Abstract:An important application of CRISPR interference (CRISPRi) technology is for identifying chemical-genetic interactions (CGIs). Discovery of genes that interact with exposure to antibiotics can yield insights to drug targets and mechanisms of action or resistance. The objective is to identify CRISPRi mutants whose relative abundance is suppressed (or enriched) in the presence of a drug when the target protein is depleted, reflecting synergistic behavior. Different sgRNAs for a given target can induce a wide range of protein depletion and differential effects on growth rate. The effect of sgRNA strength can be partially predicted based on sequence features. However, the actual growth phenotype depends on the sensitivity of cells to depletion of the target protein. For essential genes, sgRNA efficiency can be empirically measured by quantifying effects on growth rate. We observe that the most efficient sgRNAs are not always optimal for detecting synergies with drugs. sgRNA efficiency interacts in a non-linear way with drug sensitivity, producing an effect where the concentration-dependence is maximized for sgRNAs of intermediate strength (and less so for sgRNAs that induce too much or too little target depletion). To capture this interaction, we propose a novel statistical method called CRISPRi-DR (for Dose-Response model) that incorporates both sgRNA efficiencies and drug concentrations in a modified dose-response equation. We use CRISPRi-DR to re-analyze data from a recent CGI experiment in Mycobacterium tuberculosis to identify genes that interact with antibiotics. This approach can be generalized to non-CGI datasets, which we show via an CRISPRi dataset for E . coli growth on different carbon sources. The performance is competitive with the best of several related analytical methods. However, for noisier datasets, some of these methods generate far more significant interactions, likely including many false positives, whereas CRISPRi-DR maintains higher precision, which we observed in both empirical and simulated data. CRISPRi technology is revolutionizing research in various areas of the life sciences, including microbiology, affording the ability to partially deplete the expression of target proteins in a specific and controlled way. Among the applications of CRISPRi, it can be used to construct large (even genome-wide) libraries of knock-down mutants for profiling antibacterial inhibitors and identifying chemical-genetic interactions (CGIs), which can yield insights on drug targets and mechanisms of action and resistance. The data generated by these experiments (i.e., sgRNA counts from high throughput sequencing) is voluminous and subject to various sources of noise. The goal of statistical analysis of such data is to identify significant CGIs, which are genes whose depletion sensitizes cells to an inhibitor. In this paper, we show how to incorporate both sgRNA efficiency and drug concentration simultaneously in a model (CRISPRi-DR) based on an extension of the classic dose-response (Hill) equation in enzymology. This model has advantages over other analytical methods for CRISPRi, which we show using empirical and simulated data.

Exponential family measurement error models for single-cell CRISPR screens

Linear mixed effects models for non‐Gaussian continuous repeated measurement data

A model for accurate quantification of CRISPR effects in pooled FACS screens

A statistical simulation model to guide the choices of analytical methods in arrayed CRISPR screen experiments

Pooled Genome-Scale CRISPR Screens in Single Cells

Modeling variable guide efficiency in pooled CRISPR screens with ContrastiveVI+

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

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

CRISPR-MIP replaces PCR and reveals GC and oversampling bias in pooled CRISPR screens

A dose-response model for statistical analysis of chemical genetic interactions in CRISPRi screens

A new Bayesian factor analysis method improves detection of genes and biological processes affected by perturbations in single-cell CRISPR screening

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

Robust differential expression testing for single-cell CRISPR screens at low multiplicity of infection

Inference of heterogeneous effects in single-cell genetic perturbation screens

Linking Genotypes with Multiple Phenotypes in Single-Cell CRISPR Screens

Improving Accuracy in Cell-Perturbation Experiments by Leveraging Auxiliary Information

CEDA: integrating gene expression data with CRISPR-pooled screen data identifies essential genes with higher expression

Integrative Analysis of Pooled CRISPR Genetic Screens Using MAGeCKFlute.

Reliable and efficient parameter estimation using approximate continuum limit descriptions of stochastic models