# Genome-wide CRISPR-dCas9 screens in E. coli identify essential genes and phage host factors

François Rousset, Lun Cui, Elise Siouve, Christophe Bécavin, Florence

Depardieu, David Bikard

► To cite this version:

François Rousset, Lun Cui, Elise Siouve, Christophe Bécavin, Florence Depardieu, et al.. Genome-wide CRISPR-dCas9 screens in E. coli identify essential genes and phage host factors. PLoS Genetics, 2018, 14 (11), pp.e1007749. 10.1371/journal.pgen.1007749. pasteur-01975438

HAL Id: pasteur-01975438

https://pasteur.hal.science/pasteur-01975438

Submitted on 9 Jan 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Genome-wide CRISPR-dCas9 screens in *E. coli* identify essential genes and phage host factors

François Rousset<sup>1,2,●</sup>, Lun Cui<sup>1,●</sup>, Elise Siouve<sup>1,■</sup>, Christophe Becavin<sup>3</sup>,
Florence Depardieu<sup>1,□</sup>, David Bikard<sup>1,*</sup>

1 Synthetic Biology Group, Microbiology Department, Institut Pasteur, Paris, France, 2 Sorbonne Université, Collège Doctoral, Paris, France, 3 Hub Bioinformatique et Biostatistique, Institut Pasteur - C3BI, USR 3756 IP CNRS, Paris, France

● These authors contributed equally to this work.

■ Current address: Department of Chemical Engineering and Biotechnology, University of Cambridge, Cambridge, United Kingdom

* david.bikard@pasteur.fr

<figure><img src="image_2.png" /></figure>

## Abstract

High-throughput genetic screens are powerful methods to identify genes linked to a given phenotype. The catalytic null mutant of the Cas9 RNA-guided nuclease (dCas9) can be conveniently used to silence genes of interest in a method also known as CRISPRi. Here, we report a genome-wide CRISPR-dCas9 screen using a starting pool of ~ 92,000 sgRNAs which target random positions in the chromosome of *E. coli*. To benchmark our method, we first investigate its utility to predict gene essentiality in the genome of *E. coli* during growth in rich medium. We could identify 79% of the genes previously reported as essential and demonstrate the non-essentiality of some genes annotated as essential. In addition, we took advantage of the intermediate repression levels obtained when targeting the template strand of genes to show that cells are very sensitive to the expression level of a limited set of essential genes. Our data can be visualized on CRISPRbrowser, a custom web interface available at crispr.pasteur.fr. We then apply the screen to discover *E. coli* genes required by phages λ, T4 and 186 to kill their host, highlighting the involvement of diverse host pathways in the infection process of the three tested phages. We also identify colanic acid capsule synthesis as a shared resistance mechanism to all three phages. Finally, using a plasmid packaging system and a transduction assay, we identify genes required for the formation of functional λ capsids, thus covering the entire phage cycle. This study demonstrates the usefulness and convenience of pooled genome-wide CRISPR-dCas9 screens in bacteria and paves the way for their broader use as a powerful tool in bacterial genomics.

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Rousset F, Cui L, Siouve E, Becavin C, Depardieu F, Bikard D (2018) Genome-wide CRISPR-dCas9 screens in *E. coli* identify essential genes and phage host factors. PLoS Genet 14(11): e1007749. https://doi.org/10.1371/journal.pgen.1007749

**Editor:** Melanie Blokesch, Swiss Federal Institute of Technology Lausanne (EPFL), SWITZERLAND

**Received:** June 1, 2018

**Accepted:** October 9, 2018

**Published:** November 7, 2018

**Copyright:** © 2018 Rousset et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** The sequencing data is available from the European Nucleotide Archive with the accession number PRJEB28256. The processed screen results are provided as supplementary tables. The code and data tables used to perform the analyses are available at https://gitlab.pasteur.fr/dbikard/dCas9_genome_wide-screen. The data from the essentiality screen can also be visualized on CRISPRbrowser, a web-based genome browser available at crisprgo.pasteur.fr. Other relevant data supporting the

## Author summary

Over the past few years, CRISPR-Cas technologies have emerged as powerful tools to edit genomes and modulate gene expression. They have been applied to perform high-throughput genetic screens with the purpose to understand the function of genes in a
findings of the study are available in this article and
its Supporting Information files.

**Funding:** This work was supported by the European Research Council (ERC) under the Europe Union's Horizon 2020 research and innovation program (grant agreement No [677823]), by the French Government's Investissement d'Avenir program and by Laboratoire d'Excellence 'Integrative Biology of Emerging Infectious Diseases' (ANR-10-LABX-62-IBEID), and by the Pasteur-Weizmann consortium; F.R is supported by a doctoral scholarship from Ecole Normale Supérieure. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

systematic manner, but the application of these screens to bacteria have so far remained limited. Here, we present the use of a library of ~92,000 guide RNAs directing the dCas9 protein to silence one by one all the genes in the chromosome of *E. coli*. To benchmark our method, we first investigate the performance of the technique to identify essential genes, highlighting several non-essential genes also found to be essential by other methods. We then apply our method to detect bacterial genes required by three different bacteriophages to kill *E. coli* and for the production of functional progeny by phage λ. Our screens highlight previously known and new genetic interactions between phages and their host's pathways and emphasize the importance of bacterial capsule in the resistance to multiple phages. Altogether, our results demonstrate the usefulness of genome-wide CRISPR-dCas9 screens in bacteria to uncover genes involved in various phenotypes.

Introduction

The technological applications of RNA-guided nucleases derived from the Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR) prokaryotic immune system [1-3] represents a true paradigm shift in our ability to manipulate cells at the genetic level [4]. In particular, the Cas9 nuclease from type II systems can be guided by a chimeric single-guide RNA (sgRNA) which directs it to specifically cleave target sequences [5]. The ease with which these tools can be reprogrammed enables the development of powerful high-throughput screens. Recent studies have shown how an engineered CRISPR system packaged in lentiviral vectors can be used to deliver libraries of guide RNAs to human cells and create libraries of gene knockouts [6,7]. Such libraries can target most genes in the genome and were used, among other things, to identify essential genes in the human genome [8-12] and genetic requirements for different human viruses [13-15].

This approach is not directly applicable to bacteria where Cas9 cleavage leads to cell death rather than gene knockout [16-20]. The catalytic dead variant of Cas9, known as dCas9, carries mutations inactivating its two catalytic domains. It binds specific targets without introducing a double-strand break and can conveniently be used to silence genes in bacteria [21,22]. The strand orientation of dCas9 binding does not seem to matter when blocking transcription initiation by targeting a promoter sequence; however binding of the guide RNA to the non-template (coding) strand is required to block the running RNA polymerase.

Arrayed libraries of a few hundreds of guide RNAs have already proven useful to decipher the function of essential genes in *Bacillus subtilis* and *Streptococcus pneumoniae* [23,24]. While arrayed libraries enable to access many phenotypes and easily link them to genotypes, they are cumbersome to work with and maintain. Pooled screens present the major advantage of enabling the study of much larger libraries at a low cost, using well-established high-throughput sequencing methods [8-15,25].

We recently constructed a pooled library of ~ 92,000 sgRNAs targeting random positions in the genome of *E. coli* MG1655 with the sole requirement of a proper NGG protospacer adjacent motif (PAM) [26]. This library enabled us to investigate the properties of CRISPR-dCas9 screens in *E. coli* in an unbiased way and identify several important design rules to avoid off-target effects and toxicity of dCas9. In this study, we now analyze the results of the screen taking these design rules into account to investigate what it can teach us about gene essentiality in *E. coli* during growth in rich medium. We show that this method can be used to confidently predict gene essentiality in most cases despite the polar effect produced by dCas9. Among other findings, we further reveal the importance of some repeated elements, identify essential
genes that cannot tolerate small reductions in their expression levels, and challenge the essenti-
ality of a few genes.

As a proof of concept of the usefulness of this approach, we then applied the screen to iden-
tify *E. coli* genes required for bacteriophage infection. This knowledge can help identify genes
that bacteria can mutate to become phage-resistant, an interesting insight for the design of
improved phage therapies. More generally, the study of phage-host interactions over the last
century has led to many basic discoveries in genetics and molecular biology, as well as to the
development of powerful molecular biology tools [27,28]. In addition to classical genetic
approaches, genome-wide screens have already been performed to identify host dependencies
of *E. coli* phages T7 and λ using the Keio collection, an in-frame single-gene knockout strain
collection [29-31]. While powerful, this method shows several limitations, including the time-
consuming development of such strain collection, the laborious screening process and the
focus on nonessential genes only. Here, the application of our CRISPR-dCas9 screen to the
study of phages λ, T4 and 186 revealed various host factors including phage receptors and lipo-
polysaccharide (LPS) requirements. It also highlighted capsule synthesis as a shared resistance
mechanism to the three phages. We finally took advantage of the ability of phage λ to package
plasmids carrying a cos site to perform a pooled transduction assay of the library, enabling to
distinguish host genes used by phage λ to lyse its host from host genes required for the produc-
tion of functional λ particles.

**Results**

Screen design

In our previous work, we performed CRISPR-dCas9 screens using a library 92,000 guide RNAs grown over 17 generations in rich medium. We constructed to this end strain LC-E75, a derivative of *E. coli* MG1655 carrying dCas9 on its chromosome under the control of an aTc-inducible promoter [26]. In this strain, the expression level of dCas9 was fine-tuned to limit a surprising phenomenon that we termed the “bad-seed” effect: dCas9 appears to be toxic to *E. coli* when guided by sgRNAs sharing some specific 5 nt PAM-proximal sequences. This effect is particularly pronounced at high dCas9 concentrations and could be alleviated by tuning dCas9 levels while maintaining strong on-target repression. Tuning dCas9 concentration also enabled to limit off-target repression of genes, which we found could occur with as little as 9 nt of identity between the PAM-proximal region of a guide and an off-target. We now analyze the data obtained from this screen to study gene essentiality in *E. coli*. During this experiment, guides that reduce the cell fitness for instance by blocking the expression of essential genes are depleted from the library (Fig 1A). The DNA from the sgRNA library was extracted and sequenced before and after dCas9 induction. We used the number of reads as a measure of the abundance of each guide in the library and computed the log2-transformed fold change (log2FC) from DESeq2 R package [32] as a measure of relative sgRNA fitness.

Effect of guides targeting multiple positions

While most guides in our screen target unique positions in the genome, our library also
includes guides targeting multiple positions. These guides can be conveniently used to investi-
gate the role of genes or genomics regions present in several copies in the genome, which can-
not be easily achieved by other methods. Interestingly, 18.0% (348/1,932) of sgRNAs with
multiple targets around the genome are strongly depleted (log2FC < -2) against 7.4% (6,117/
82,816) for guides having a single target (p < 10⁻⁴¹, Fischer’s exact test). Most of these depleted
guides target genes involved in translation: ribosomal RNAs, tRNAs, and the elongation factor
Tu encoded by highly similar *tufA* and *tufB* genes (S1 Fig). More surprisingly, 37 of these
<figure><img src="image_4.png" /><figcaption>**Fig 1. Predicting *E. coli* essential genes from the CRISPR-dCas9 screen data.** (A) A sgRNA library was transformed into cells expressing dCas9 under the control of an aTc-inducible promoter. The distribution of sgRNAs was retrieved by deep sequencing before dCas9 induction and after 17 generations of induction. (B) Genome-wide visualization of the screen results. The median log2FC of sgRNAs targeting the coding strand of protein-coding genes is displayed in grey. The outer red track represents genes annotated as essential in the EcoGene database while blue track shows genes selected as candidate essential in our screen (median log2FC < -2) (C) Genes were ranked according to the median log2FC of the sgRNAs targeting the coding strand of protein-coding genes. Genes annotated as essential in the EcoGene database are shown in red. Dashed line represents the chosen essentiality threshold (median log2FC = -2). Error bars represent median absolute deviation. (D) The median log2FC of sgRNAs targeting the coding strand of protein-coding genes was used to predict gene essentiality. Receiver Operating Characteristic (ROC) curve of the prediction model</figcaption></figure>
is plotted (AUC = 0.943). The black dot represents the chosen threshold. The index score from the TraDIS dataset was also used to predict the EcoGene essential gene set (AUC = 0.957). (E) 379 candidate essential genes were chosen with the threshold. Polar effect indicates genes located upstream essential or near-essential genes. Slow-growth phenotype indicates non-essential genes known to have a growth defect when mutated or deleted. (F) 52 genes annotated as essential were not classified as essential in our screen. Genes found essential only in the Keio collection but not by other sources are termed as "likely not essential". Conflicting results indicates genes with discrepancies between various datasets. (G) Venn diagram comparing essential genes selected in our screen, in the TraDIS study [35] and in the EcoGene database [34]. (H) Within an operon, different scenarios are displayed as well as the corresponding conclusion. If a fitness defect is induced by guides targeting the first gene but not the second gene, then the first gene can be confidently classified as essential and the second as non-essential (top panel). If a fitness defect is induced when targeting both genes, then the second gene can confidently be classified as essential while the first remains uncertain (middle panel). Finally, if a fitness defect is induced by guides targeting the second gene but not the first gene, then an internal promoter likely drives the expression of the second gene (bottom panel).

https://doi.org/10.1371/journal.pgen.1007749.g001

guides target Repetitive Extragenic Palindromic (REP) elements. Some REP elements were reported to play a regulatory role at the translational level by modulating mRNA stability [33], but their exact role remains elusive. For many of these guides, none of the target positions are in the vicinity of essential genes, suggesting that the fitness defect caused by dCas9 binding at REP sequences is not simply due to an effect on the expression of essential genes. Guides targeting several repeat regions at the same time had a significantly higher chance of inducing a strong fitness defect (log2FC < -2) than guides targeting only one repeat region (p < 2.10⁻⁵, Fischer's exact test). Further work will be required to understand the mechanism responsible for the fitness defect produced by REP-targeting guides, which will likely reveal interesting biology.

## Identification of essential genes

Starting from the initial library of ~ 92,000 sgRNAs targeting random positions along the chromosome of *E. coli*, we filtered sgRNAs with either a bad seed effect, an off-target, multiple targets in the genome or a low number of reads (see Methods), yielding a library of ~ 59,000 guides used to perform the analyses below (S1 Table). To investigate the usefulness of such CRISPR-dCas9 screens to identify essential genes, we first ranked genes according to the median log2FC of sgRNAs targeting the coding strand (S2 Table). Mapping our data to the genome of *E. coli* highlighted a good concordance with previously known essential genomic regions (Fig 1B). As expected, genes with the most depleted sgRNAs included a large majority of genes annotated as essential in the EcoGene database [34] (Fig 1C). A gene network analysis of the 100 top-scoring genes in our screen highlighted the main essential functions, namely ribosome assembly, peptidoglycan synthesis, DNA replication and transcription, fatty acid metabolism and tRNA metabolism (S2 Fig). We then assessed the ability to predict gene essentiality from the median log2FC value, using the EcoGene database as a gold standard. The receiver operating characteristic (ROC) curve showed an area under the curve (AUC) of 0.943 (Fig 1D). Our CRISPR-dCas9 screen achieved a similar prediction performance as a recent screen based on Transposon-Directed Insertion-site Sequencing (TraDIS) using the largest ever built transposon library (901,383 insertion sites, AUC = 0.957) [35]. This suggests that CRISPR-dCas9 screens can reach the efficacy of transposon-derived methods but with much smaller libraries (~ 15 fold in this study). In particular, during the review process, a new study using genome-wide CRISPRi screening in *E. coli* reported similar conclusions by achieving a comparable prediction performance with a rationally designed library of ~ 60,000 sgRNAs, thus highlighting the reproducibility of CRISPRi screens [25].

Candidate essential genes were selected when the median log2FC of sgRNAs targeting their coding strand was lower than -2. This arbitrary threshold enabled to select a consistent number of candidate genes while limiting the selection of false positives (see Fig 1D). We selected 379 genes including 235 genes annotated as essential (62%), representing 79% of the annotated essential genes (Fig 1E). Among the 144 remaining genes, 49 are known to show a slow-growth
phenotype when mutated or deleted [36] and can thus be termed as near-essential. The effect of 57/144 genes can be explained by a polar effect leading to the repression of a downstream essential or near-essential gene in the same operon. All in all, polar effects account for a false-positive discovery rate of essential genes of 15%. In addition, 2/144 genes, *hipB* and *ratA* are part of toxin/antitoxin systems. Blocking the expression of *hipB* will also silence the downstream *hipA* toxin. This is expected to be toxic as degradation of HipB by the Lon protease will free HipA [37]. The antitoxin of *ratA* was not described and its regulation is poorly understood but we can make the hypothesis that guides targeting this gene also block the expression of the antitoxin. Among the 36 remaining genes, 5 genes can be explained by off-target effects. When filtering the data to eliminate guides with off-targets, we made a compromise between the number of guides to keep in the analysis and the stringency of the filter. As a result, a few guides with off-target effects remain in the library. Since the median log2FC was used to assess gene essentiality, a gene targeted by only one or two guides might be classified as essential if one of the guides has an off-target. Two such examples are given in S3 Fig. Finally, 31 genes remain ambiguous either because they are targeted by a very low number of sgRNAs in our library or because these genes may be essential or near-essential in our experimental conditions. Note that 7 of these genes (*aceF*, *lpd*, *dcd*, *hemE* and *ihfA*) were found to be essential in the recent TraDIS screen [35]. Overall comparison of our data with the EcoGene database and with the TraDIS screen revealed 227 essential genes common to all datasets while the genes unique to our dataset are mostly explained by polar effects (44/99) and slow-growth phenotypes (25/99) (Fig 1G).

Note that the analysis above was performed by computing a median log2FC value for each gene without taking into account their genomic context. Silencing a gene in an operon with dCas9 leads to the repression of all the downstream genes, and in many cases, several genes with a fitness defect are adjacent within an operon. In this scenario, only the last gene with a fitness defect can confidently be classified as essential or near-essential due to polar effects. Upstream genes might or might not cause a fitness defect when silenced alone, and can be classified as potentially essential (Fig 1H). When doing the analysis in this manner rather than gene by gene, only 155 genes can be confidently classified as essential while 224 genes remain uncertain.

Another interesting scenario when looking at operons is when silencing a gene upstream of an essential or near-essential gene induces no fitness defect despite the theoretical polar effect. A likely explanation is that an unidentified promoter internal to the operon drives the expression of the essential gene independently of the upstream gene (Fig 1H). We identified 7 such cases (S4 Fig), and in all instances recent data supports the existence of an internal promoter [38].

Conversely, some genes annotated as essential do not show a strong fitness defect in our screen (median log2FC > -2) (Fig 1F, S3 Table). Out of 52 essential genes that we fail to identify, 4 genes (*chpS*, *mazE*, *yafN* and *yefM*) are part of toxin-antitoxin systems which are encoded in operons with the antitoxin gene preceding the toxin gene. By definition, an antitoxin is essential to the cell as silencing it leads the cognate toxin to kill the cell. In the four cases above, dCas9 blocks the expression of both the antitoxin and the toxin gene, which we expected to be toxic if the antitoxin is more labile. Surprisingly, unlike *hipB*, no effect on fitness was measured here. Another 13/52 genes annotated as essential in the Ecogene database were also supposed to be non-essential by other studies [35,39,40]. We successfully replaced three of these genes, *alsK*, *bcsB* and *entD* by a kanamycin resistance cassette, confirming that they are indeed not essential (S5 Fig). In 10/53 other cases, gene essentiality has also been challenged by conflicting results [31,35,39–41]. For the remaining 25 genes, the absence of fitness defect could be explained by an inefficient repression by dCas9, a strong robustness of the cell to low
levels of the protein or could indicate genes that are actually not essential in our experimental conditions. In addition, weak repression could arise from negative regulatory feedback loops. We recently showed that silencing a gene regulated by a negative feedback leads to an increased transcription initiation rate, ultimately leading to weak silencing [42]. Here, we specifically investigated two essential genes, *lexA* and *rho*, which are not detected as essential in our screen and whose expression is known to be controlled by a negative feedback loop: *lexA* binds its own promoter while the expression of *rho* is attenuated by the presence of Rho-dependent terminators when Rho levels are high. As expected, targeting these genes only led to a weak repression level as measured by RT-qPCR (S6 Fig). When cloning a sfGFP reporter under the control of the *lexA* promoter or the *rho* leader sequence, we could confirm that blocking *lexA* or *rho* led to an increased sfGFP expression (S6 Fig). We can make the prediction that guide RNAs binding in the promoter region of *lexA* should not be affected by the negative feedback loop since they inhibit the initiation of transcription. Indeed, a guide targeting the promoter of *lexA* does show a strong fitness defect (no guide targeting the promoter of *rho* are present in the screen).

## Unexpected effects on the template strand of target genes

We further investigated the effect of sgRNAs targeting the template strand of genes. Guides in this orientation should only have a moderate effect on gene expression. As expected, the vast majority of them do not produce any fitness defect. We compared the median log2FC of sgRNAs targeting the coding strand and the template strand (Fig 2A). Interestingly, we observed that a set of essential genes that we term “sensitive” produced a strong fitness defect (median log2FC < -2) when targeted on the template strand (including *accA*, *alaS*, *ftsA*, *ftsQ*,

<figure><img src="image_5.png" /><figcaption>**Fig 2.** Identification of essential genes that do not tolerate small changes in their expression level. (A) Comparison between median log2FoldChange of sgRNAs targeting the coding strand or the template strand highlights unexpected effects. (B) Effect of guide RNAs targeting the *glyQS* operon. (C) Effect of guide RNAs targeting *gyrA*. (B,C) Relative expression of *glyQS* and *gyrA* when targeted by guides T-glyQ1, T-glyQ2, T-gyrA1 and T-gyrA2 as measured by RT-qPCR. Black bars show the mean of 3 independent biological replicates and 2 technical replicates.</figcaption></figure>

https://doi.org/10.1371/journal.pgen.1007749.g002
glyQ, glyS, lolB, lptC, rplX, rpmB, rpsL and yejM) whereas targeting the template strand of other essential genes did not induce any fitness defect. To verify that this effect does not occur from a surprisingly high repression level on the template strand of these genes, we measured the repression efficiency of 2 guides targeting the template strand of the "sensitive" gene *glyQ* (Fig 2B). As a control, we also measured the repression efficiency of 2 guides targeting the template strand of the essential gene *gyrA* which does not show this sensitivity phenotype (Fig 2C). The results show similar intermediate repression levels in both cases (33-57% of the initial expression level, see Fig 2B and 2C) suggesting that the cell indeed seems very sensitive to the expression level of some essential genes whereas a weak repression of most essential genes can be tolerated.

The other genes that display a strong fitness defect when targeted on the template strand can be explained by polar effects on nearby essential or near-essential genes and shed light on atypical gene organizations (S7 Fig). Gene *yejL* is located upstream of the sensitive essential gene *yejM*. Gene *ypaB* contains the promoter of the essential gene *nrdA*. Gene *sroG* is actually a riboswitch controlling the expression of the essential gene *ribB* involved in riboflavin synthesis [43]. One sgRNA binding to the template strand in the very beginning of the gene was strongly depleted, suggesting that it blocks transcription initiation of the *ribB* mRNA. *yceQ* is a small open reading frame of just 321bp in the promoter region of the essential RNAse E gene (*rne*) and in the opposite orientation. Guide RNAs directing dCas9 to bind the template strand of *yceQ* are thus in the good orientation to effectively block the expression of *rne*. Conversely, guides targeting the coding strand of *yceQ* show no fitness defect, suggesting that *yceQ* is actually not essential, which is also supported by recent TraDIS data [35]. A similar case is that of *psrO* which encodes a small RNA located in the promoter region of *pnp*.

Overall, these results show the performance of CRISPRi screens to assess gene essentiality in E. coli but also to highlight diverse genomic organizations. Our results can be visualized on CRISPRbrowser, a web-based genome browser available at crispr.pasteur.fr.

## Identification of genes providing phage resistance when silenced

To demonstrate the broad usefulness of this strategy to the study of different phenotypes, we applied the method to unveil bacterial genes required for successful phage infection, also known as host factors. Considering previous results, we now focused on sgRNAs targeting genes on the coding strand, yielding a library of ~ 17,220 sgRNAs after removing guides with an insufficient number of reads. As a proof of concept, we used temperate phage λ whose host requirements are well documented. We also used temperate phage 186cIts (a thermosensitive strain of phage 186) and virulent phage T4 in order to compare their host requirements. As strain LC-E75 carries the *dcas9* expression cassette integrated at the *attB* site of phage 186, we constructed a new strain, FR-E01, with the same cassette integrated at the *attB* site of phage HK022 in order to avoid any interference with phage 186. Both strains expressed dCas9 at the same level (S8 Fig).

A culture of strain FR-E01 carrying the library was grown with aTc to stationary phase allowing for dCas9 to be expressed and for the target gene products to be diluted and/or degraded. Cells were then diluted 100-fold and grown to exponential phase, still with aTc, followed by infection with phage λ, T4 or 186cIts at a multiplicity of infection (MOI) of 1 (Fig 3A). During the experiment, phages will lyse bacteria unless the sgRNA they carry makes them resistant to infection. The pool of sgRNAs was sequenced before and after 2 h of infection to allow phages to lyse a maximum of sensitive cells while limiting the rise of resistant mutants in the population. Log2FC values were computed for each sgRNA as previously (S4 Table and
A

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

B

Phage λ

Phage T4

Phage 186

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

C

Phage λ

Phage T4

Phage 186

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

D

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

**Fig 3.** The CRISPRi screen reveals phage receptors and LPS requirements of phages λ, T4 and 186. (A) Overview of the experimental procedure. dCas9 expression was induced for ~ 8h before infection by λ, T4 or 186 at a MOI of 1. The distribution of sgRNAs was retrieved before and after infection by sequencing. (B) The log2FoldChange was computed for each sgRNA. For phages λ and T4, sgRNAs targeting the phage receptor (respectively LamB and OmpC) are highlighted in red as well as genes targeting waa genes for phage 186. Dashed line represents adjusted p-value = 10⁻³. (C) For each gene, a resistance score was computed as the median log2FoldChange of sgRNAs targeting the coding strand. Error bars represent median absolute deviation. 5 top-scoring genes are highlighted in red in each panel. (D) Infection time courses of MG1655 (top panel) or Keio collection strain ΔwaaJ without phage or after infection by λ, T4 or 186 (MOI ~ 0.01).
Fig 3B). For each gene, a resistance score was computed as the median log2FC of the guides targeting the coding strand (S5 Table).

Host genes are required at several steps of the phage infectious cycle and the genes involved differ depending on each phage. The first step of phage infection is adsorption to the cell followed by receptor binding and injection of the phage DNA. Surface proteins or specific structures of the LPS are typically used as receptors. The gpJ tail tip protein of phage λ binds to the maltose outer membrane porin LamB [44], while T4 recognizes both the OmpC protein and the core-lipid A region of the LPS that includes the heptose region [45]. In addition to surface receptors, in the case of λ, the ManXYZ mannose permease helps the phage DNA cross the inner membrane [46]. The surface receptor of phage 186 or other genes required for DNA entry are not known. As expected, the highest-scoring genes in our screen correspond to phage receptors and their positive regulators (Fig 3B and 3C). The phage λ receptor gene *lamB* and its regulators *malT*, *cyaA* and *crp* respectively ranked 3rd, 2st, 6th and 9th after infection by λ. The gene with the 1st highest resistance score was *malK* which is expressed upstream of *lamB* in an operon structure. Guides that block the expression of this gene therefore also block the expression of *lamB*, explaining its high ranking although it is not known to participate in the infection process. Many genes involved in LPS biosynthesis (*lpcA*, *gmhB*, *rfaD*, *rfaE*, *waaA*, *waaC*, *waaF*, *waaG*) were also identified in our screen. Mutations affecting LPS synthesis were previously associated with reduced adsorption of λ, presumably linked to a reduced availability of LamB on the outer membrane [47]. Guides that target the manXYZ operon where also significantly albeit weakly enriched (p < 0.0034, Mann-Whitney U test). In the case of phage T4, gene *ompC* and its regulator *ompR* ranked 1st and 2nd while genes *rfaD* and *waaF* involved in the synthesis of glycerol-D-manno-heptose and its addition to the LPS, ranked 5th and 8th. Considering that the receptor for phage 186 is not clearly established, we used this approach to identify candidate host components. Out of the 20 highest scoring genes, 16 were involved in LPS biosynthesis while no surface proteins were present, suggesting that 186 uses the LPS as a receptor (Fig 3C). To validate this observation, we performed infection time courses with a strain deleted for *waaJ*, one of the last genes in the LPS biosynthesis pathway, in the presence of each of the three phages. The absence of gene *waaJ* induced a strong resistance to 186 but not to λ and T4 (Fig 3D) suggesting that the outer core of the K-12 strain LPS is required for adhesion of phage 186. Taken together, these observations validate the ability of our method to identify phage receptors which are the most common bacterial components giving rise to phage resistance.

Comparison between phages λ, T4 and 186 revealed closer host requirements between the two temperate phages, λ and 186, than between λ and T4 or between 186 and T4 (Fig 4A). Genes whose silencing provide resistance to the different phages were selected if their resistance score was at least 20% of the maximum score for each phage (see Methods, Fig 4B). Interestingly, we could identify 3 genes with a high resistance score to all 3 phages. These include *rfaD* and *waaF* discussed above for their involvement in LPS biosynthesis. In addition, we identified gene *yrfF*, an essential gene indirectly linked to capsule synthesis. YrfF inhibits the Rcs signaling pathway by an unknown mechanism [48]. This signaling pathway senses and responds to damage to the membrane or to the peptidoglycan by activating genes involved in colanic acid capsule synthesis [49,50]. We hypothesized that silencing *yrfF* would lead to the activation of the Rcs signaling pathway, induce capsule synthesis and prevent phage adsorption to the cell surface. The synthesis of capsule was indeed previously reported to provide broad resistance to phages [51–54]. Since *yrfF* gene is essential, we reduced its expression to intermediate levels using a sgRNA bearing 4 mismatches at the 5'-end (*yrfF*-4m) instead of using a fully matched sgRNA [21,42]. The resulting strain only showed a slight growth defect and became resistant to all three phages (Fig 4C). Deletion of gene *rcsB*, a central actor of the Rcs
<figure><img src="image_16.png" /><figcaption>**Fig 4.** Colanic acid capsule synthesis is a common resistance mechanism to phages λ, T4 and 186. (A) Pearson correlation coefficient of resistance scores between phages λ, T4 and 186-cIts. (B) For each phage, genes targeted by > 1 guide with a resistance score greater than 20% of the maximum score were selected and compared in a Venn diagram. (C) Infection time courses of FR-E01 carrying an empty vector, FR-E01 expressing *yrfF* sgRNA with 4 mismatches (*yrfF*-4m) or FR-E01Δ*rcsB* expressing *yrfF*-4m without infection or after infection by λ, T4 or 186-cI (MOI ~ 0.01). Experiments were performed in triplicates. (D) Expression of genes *wza* and *wzb* involved in colanic acid capsule synthesis in strain FR-E01 carrying an empty vector or expressing the *yrfF*-4m sgRNA, as measured by RT-qPCR. Results are shown for 3 biological replicates and 2 technical replicates. (E) Schematic of colanic acid capsule-mediated phage resistance. *YrfF* inhibits the *Rcs* phosphorelay. Inhibition of *YrfF* expression activates the *Rcs* phosphorelay which triggers the expression of genes involved in colanic acid capsule synthesis.</figcaption></figure>

https://doi.org/10.1371/journal.pgen.1007749.g004

pathway, restored sensitivity to the 3 phages, showing that the resistance phenotype that results from *yrfF* silencing is mediated by the *Rcs* pathway. We further confirmed this by showing that silencing *yrfF* induces a ~ 140-fold increase in the expression of the *wza-wzb* operon directly involved in colanic acid capsule synthesis (Fig 4D and 4E). Interestingly, deletion of *rfaD* or *waaF* (but not other genes of the LPS biosynthesis pathway) was shown to induce a mucoid phenotype through the synthesis of a capsule [55]. This suggests that they might provide resistance to all three phages through this pathway rather than through their role in LPS synthesis. Consistently with this hypothesis, silencing genes upstream of *rfaD* and *waaF* in the LPS biosynthesis pathway did not provide resistance to T4 in our screen.

After DNA entry, temperate phages make a decision to either enter the lytic or lysogenic cycle. In the lysogenic cycle, the phage integrates in the host genome and the cell survives. Host factors that increase the rate of lysogeny when repressed are thus expected to be enriched in our screen. As a matter of fact, we identify genes *hflD* and *ftsH* (*hflB*) after λ infection. These genes were historically described through the identification of mutants that affect the lysis/lysogeny decision of λ (*hfl* stands for “high-frequency of lysogeny”) [56]. HflD is known to
directly interact with the transcription factor CII of λ and facilitate its degradation by the FtsH-HflKC protease [57]. Repressing *hfI*D or *ftsH* is thus expected to increase the level of CII which will push the phage toward lysogeny and increase survival. Genes *hfIKC* did not pass our selection threshold for an unknown reason. Interestingly, *hfI*D is also found for phage 186 suggesting a similar control of the lysogeny decision for both phages.

In lytic cycle, phages replicate their DNA to high copy number and produce the phage capsids into which the DNA is packaged. Finally, the cell is lysed leading to phage release. Host genes that are required for any of the steps of the lytic cycle are expected to be enriched in our screen as long as silencing their expression prevents cell death. After λ infection, we could identify genes *dnaKJ* encoding a chaperone known to be required for DNA replication of λ [58]. These genes are also required for phage 186 which suggests a similar replication mechanism. Genes involved in two tRNA modification pathways were also identified after λ infection: uridine 2-thiolation (*tusA*, *tusBCD*, *tusE* and *mnmA*) and the essential N⁶-threonylcarbamoyladenosine modification (*tsaB*, *tsaC*, *tsaD* and *tsaE*). Uridine 2-thiolation has been previously described to be necessary for the proper translation of λ proteins gpG and gpGT involved in tail assembly through its impact on programmed ribosomal slippage [29,59], but the essential N⁶-threonylcarbamoyladenosine modification of tRNAs was not previously associated with λ infection. We also identified RNAses *rnt*, *rnpA* and *rnd* involved in tRNA maturation. Altogether, 27 out the 100 top-scoring genes after λ infection are involved in translation either directly or through ribosome or tRNA synthesis. These include genes also found after infection by phage 186 such as *deaD*, *hpf*, *prfA*, *prmC*, *rluD*, *rnpA*, *tilS* and *trmH*. This highlights the importance of an intact translation machinery for the completion of the lytic cycle for both λ and 186. Finally, genes linked to DNA repair including *recA*, *recBCD*, *recG*, and *ligA* seem to be important for 186 infection but not for λ infection. To our knowledge, no recombinase has been identified in the genome of phage 186 suggesting that it may rely on the host recombination system while phage λ carries its own recombination machinery (*gam*, *exo* and *bet*) [60].

## Identification of genes involved in later steps of λ infection

The screen performed here is thus a powerful method to identify genes required by phages to kill the cell. However, this first strategy cannot identify genes necessary for the synthesis of functional phage capsids if blocking the expression of these genes does not prevent the phage from killing the cell. One can expect that this will be the case of any host gene involved in late stages of the infectious cycle when the host cell is already doomed. In order to get better insights into the genes affecting the production of functional phages, we implemented a second step focusing on phage λ (Fig 5A). The vector carrying the sgRNA library was designed to carry a λ packaging site (*cos* site), thus forming a cosmid [61]. Upon infection, phage λ is thus able to package the plasmid carrying the guide RNA only if a functional capsid was produced (S9 Fig). After infection, the cell lysate containing a mixture of λ and cosmid particles was used to transduce strain MG1655::λ and thus recover guides in the library that do not affect the infection process. The distribution of sgRNAs recovered after transduction was compared to the initial pool to identify depleted sgRNAs corresponding to bacterial genes required for the production of functional capsids (S6 Table).

Interestingly, 45 guides in our library which target prophages of *E. coli* MG1655 also have a match with 100% nucleotide identity and a proper PAM sequence in the genome of phage λ. Among these guides, those predicted to block the expression of the structural genes of phage λ should inhibit the production of functional capsids. As expected, these guides were strongly
<figure><img src="image_17.png" /><figcaption>**Fig 5.** Transduction of the library reveals host genes necessary for the production of functional capsids. (A) Overview of the experimental system to screen for λ host factors. Strain FR-E01 expressing dCas9 under aTc-inducible promoter and carrying the sgRNA library on a cosmid was infected by λ. The library can be packaged into capsids when the sgRNA does not disrupt phage propagation. After transduction into a lysogenic strain, the distribution of sgRNAs is measured by sequencing. (B) log2FC values after transduction are strongly correlated to the effect of guides on cell growth (ANOVA, F = 27263, df = 17198, p < 10⁻¹⁶). (C) For each gene, a linear model was built to explain log2FoldChange after transduction while taking the effect on cell growth into account. Estimate and corrected p-value (FDR) are shown for each gene. Genes previously described as host factors are shown in red. (D) Comparison of results after infection and after transduction. Genes targeted by >1 guide which significantly decrease the production of functional phages (ANOVA, FDR < 0.05) are displayed in red or blue when they show a negative or positive resistance score respectively, as computed in the first screen. (E) Growth of strain FR-E01 carrying a sgRNA targeting *nusG* without infection or after infection by λ at MOI = 0.01 or MOI = 1. (F) Phage titer during infection of strain FR-E01 carrying an empty vector or expressing a sgRNA targeting *nusG* (MOI = 0.01). Error bars show standard deviation of three independent experiments.</figcaption></figure>

https://doi.org/10.1371/journal.pgen.1007749.g005

depleted after transduction (S10 Fig). Most of these guides do not block the expression of the lysis genes of λ, explaining why they were not identified in our first screen.

When comparing our data with the essentiality data from our previous screen, we observed a strong correlation between log2FC after transduction and the effect of guides on the cell fitness (Fig 5B). One might indeed expect that genes whose repression slows down cell growth will also slow down phage production, leading to fewer particles being released after 2h of infection. In order to identify genes which disproportionately affect phage production over cell growth, we performed statistical analyses taking the effect on growth into account (see Methods, S7 Table). We identified 57 genes (46 targeted by > 1 guide) which significantly
decrease the amount of transduced particles when repressed (FDR < 0.05) (Fig 5C). These genes could be separated into two groups according to their resistance score (Fig 5D): genes with a positive resistance score provided resistance to lysis when silenced, while genes with a negative resistance score led to an increased sensitivity to the phage but a reduced production of functional phage particles. The first group includes genes involved in the expression of LamB (*lamB*, *malK*, *malT*, *crp*, *cyaA*) as well as other genes identified in the first step, such as *rpoH*, *dnaJ*, *prtC*, *rnt* and *rnpA*. The second group includes genes that could not be identified in the first step of our screen based on their resistance score as they decrease or have little effect on cell survival to λ infection. These genes include well-characterized genes involved in the transcription of the λ genome: RNA polymerase subunits β *rpoB* and β' *rpoC*, RNA polymerase σ⁷⁰ factor *rpoD* and transcription termination factors *rho*, *nusA* and *nusG* [62–65]. These factors interact with λ protein N to avoid transcription termination at terminator sites present on its genome, thus allowing the expression of downstream genes in a process called transcription antitermination [63]. The identification of *rho* is somewhat surprising. Indeed, while it was found to be associated to the antitermination complex [66], it was not previously described as essential for phage λ transcription. While antitermination at Rho-dependent terminators has been extensively studied, note that *rho* is an essential protein and λ replication was never tested in a null mutant of *rho* to our knowledge.

The screen also identified *groL* encoding the GroEL chaperonin involved in the assembly of the λ head [67–69]. Its transcription is controlled by the heat shock-specific RNA polymerase σ³² factor *rpoH* found in the first group [70]. Genes with little effect on survival after λ infection include DNA polymerase III subunit α encoded by *dnaE*, as well as many genes encoding aminoacyl-tRNA synthetases, namely *argS*, *pheS*, *pheT*, *leuS*, *hisS*, *proS*, *trpS*, *valS*, *cysS*, and *metG*, suggesting that these essential elements are a limiting resource for the production of phage particles.

This analysis highlights the fact that silencing some genes inhibits the formation of phage particles while leaving the cell susceptible to lysis. To confirm this observation, infection time courses and phage titer measurements were performed in strain FR-E01 expressing a sgRNA targeting *nusG* which has a negative resistance score (Fig 5E). Infection at a MOI of 1 led to a complete lysis of the cell population showing that silencing *nusG* does not provide resistance to infection. However, infection at MOI = 0.01 did not lead to lysis of the population, presumably because the first infected cells did not release enough active phages to lyse the rest of the population. Measuring phage concentration during infection indeed showed a ~ 10-fold reduction in burst size when *nusG* is silenced (Fig 5F).

Discussion

This work establishes a proof of concept of the performance of genome-wide CRISPRi screens in bacteria for the identification of essential genes and phage host factors. During the review process, another study reported the use of a library of ~60,000 guides (~15 guides per gene) to identify essential genes in *E. coli* [25]. In contrast, our library was designed to target ~92,000 randomly chosen positions with a proper PAM along the chromosome of *E. coli* in order to investigate the properties of CRISPRi in an unbiased way. A previous screen performed with this library enabled to identify an intriguing toxicity of dCas9 mediated by 5 nucleotides in the seed sequence of the guide RNA, as well as the existence of off-target effects with as little as 9nt of identity between the guide and the off-target position [26]. After filtering our guides to avoid these effects, ~23,000 guides targeting the coding strand of genes were included in our analyses. Despite the much lower effective size of our library compared to the recently published study, both libraries achieved similar results. This highlights the importance of the
design rules established in our previous work. A current limitation of CRISPRi is the heterogeneity in sgRNA repression efficiency. Efforts have been made in Eukaryotic systems to predict sgRNA activity [71–75] but such models remain to be developed in bacteria. The ability to confidently design efficient guides will in the future enable to decrease the necessary library size even more, thus increasing the number of experiments that can be performed simultaneously while decreasing DNA synthesis and sequencing costs.

A notable feature to consider during CRISPR-dCas9 experiments is that when series of genes in operons are found to have an effect, it is theoretically only possible to conclude that the last gene in the series is essential. A gene by gene analysis of our data was nonetheless able to reliably predict gene essentiality with a false positive discovery rate due to polar effects of only 15%. This surprisingly accurate performance is due to the fact that essential genes tend to cluster in the same operons. Yet, this is not always the case and polar effects represent a limitation of what can be done with CRISPRi screens: one can easily identify operons of interest but a detailed understanding of the involvement of each gene requires further investigation with other techniques.

Overall, CRISPRi screens present several key advantages: (i) The expression of dCas9 is inducible, which enables to maintain guides targeting essential genes in the library, allowing their study; (ii) Duplicated regions can easily be studied since sgRNAs can be designed to target identical sequences present in several copies along the genome; (iii) Intermediate repression levels can be obtained either by targeting the template strand of genes or by using sgRNAs with a variable number of mismatches [21,42]; (iv) libraries can be rationally designed to target specific locations or subsets of genes, as opposed to the random insertion of transposons. (v) Preparing libraries for sequencing only requires a simple PCR reaction to amplify the guide RNAs.

In this study, we used these key advantages to investigate essential genes and other interesting features of the *E. coli* genome. Overall, 79% (235/298) of the genes previously annotated as essential were correctly identified. A large part of the genes that we failed to identify might not actually be essential. Indeed, our screen showed some discrepancies with annotated databases, particularly with the Keio collection, which were also observed in the other recent CRISPRi screen [25]. We experimentally confirmed that some genes were wrongly annotated as essential, corroborating other recent work [35]. Nonetheless, our screen failed to identify a few well-known essential genes. This might be due to a strong robustness of the cell to low levels of some essential proteins or to a weak dCas9 repression caused by negative feedback loops as we could indeed demonstrate in the case of genes *lexA* and *rho*. Note that the design of our library being random, some genes are only targeted by a small number of guides in our screen, reducing the likelihood to correctly identify them. This will easily be corrected in future screens using rationally designed libraries.

In addition to the investigation of gene essentiality, our results highlighted the importance of REP elements and the usefulness of CRISPRi to investigate such repeated loci. We could also identify genes that cannot tolerate small reductions of their expression level. Proteins encoded by these genes could be promising antibiotic targets as a partial inhibition of their activity is likely lethal. Our screen also highlighted some atypical genomic organizations such as internal promoters within operons or nonessential genes antisense to neighboring essential genes.

It is interesting to discuss the advantages and drawbacks of CRISPR-dCas9 screens in comparison with other genome-wide approaches. Arrayed collections of gene deletions have proven extremely useful for the study of the few bacterial species where they have been constructed. The *E. coli* Keio collection [31,41] is used as a reference for the classification of *E. coli* genes as essential, and has been employed in many different types of screens, including
for the identification of phage host factors [29,30]. However, the development of such strain
collections and the screening process are laborious and focused on nonessential genes only.
The Keio collection is also known to contain a few errors arising from the failure to obtain
some gene knockouts for other reasons than gene essentiality, or from gene duplications
occurring during the construction of the collection [41]. Transposon-based strategies such as
Tn-seq and TraDIS are powerful methods to identify essential and conditionally-essential
genes in a wide range of bacteria [35,76,77] by looking for genes in which no transposition
events are recovered. With these approaches, false negatives may arise when the insertion of
the transposon does not disrupt the functionality of the protein, or when only a small region
of a gene is essential. False positives might also arise from biases in the number of transposon
insertions along the genome. Several data analysis pipelines have been developed to take
these effects into account [78-80], but recent results still highlighted the need to combine
statistical analysis with manual curation in transposon-based screens [35]. Finally, both
methods can cause polar effects on the expression of adjacent genes which can be hard to
predict as opposed to the predictable effects of dCas9 binding on the expression of down-
stream genes.

The application of our CRISPRi screen to identify *E. coli* genes required for phage infection enabled us to identify well known host factors while providing novel insights into phage-host interactions. Among the highlighted host factors, phage receptors and other genes involved in phage adhesion had the strongest effect. Interestingly, some host factors were common between λ and 186 suggesting some similarities between the infectious cycles of these two temperate phages. Phage T4 relies on very few host components in comparison to λ and 186. This is consistent with the fact that T4 has a much larger genome (169 kbp) than the two other phages (48kbp for λ and 31kbp for 186) and encodes most of the functions it needs to complete its cycle [81]. *E. coli* was reported to become resistant to phage T4 only by acquiring mutations that block its adsorption [82,83], whereas *E. coli* can also become resistant to λ by acquiring mutations in intracellular components [64]. A TraDIS screen was recently performed in *E. coli* O157 to investigate genes involved in susceptibility and resistance to phage T4 [82]. This screen identified a few host factors putatively involved in DNA entry that were not found in the CRISPR-dCas9 screen, suggesting the existence of experiment or strain-specific effects. Note that T4 is known to degrade the host DNA [84]. As a consequence, pooled screens might not be able to identify host factors that act at late stages of the infection as the guide RNA or transposons inserted in the genome will be degraded. Note also that our screen failed to identify a few previously described host factors for λ infection [29], likely for the same reasons that we failed to identify a few essential genes.

A notable result of our screen is the identification of colanic acid capsule synthesis as a
shared resistance mechanism to phages λ, 186 and T4. The capsule of *E. coli* has been shown to
protect against external aggressions such as antibiotics or desiccation and helps to evade the
host immune system [51]. While previous studies have reported that bacterial capsule can also
provide resistance to several phages [52–54], this is rarely seen as the main function of the cap-
sule. Our findings argue for a broader role of colanic acid capsule synthesis as a defense mech-
anism against a wide range of phages, acting as a physical barrier by masking phage receptors
on the cell surface. Accordingly, some phages have evolved the capacity to degrade capsule
through the action of enzymes such as endosialidases which are frequently carried by the tail
fibers of caudovirales [52,85–87].

While genome-wide CRISPR screens have already demonstrated their usefulness in eukary-
otic systems [8–15], this study should now set the stage for their broader adoption as a power-
ful tool in bacterial genomics.
# Methods

## Bacterial strains and media

*E. coli* strains were grown in Luria-Bertani (LB) broth. LB-Agar 1.5% was used as a solid medium. Antibiotics were used in standard concentrations (50 µg / mL kanamycin (Sigma), 100 µg / mL carbenicillin (Euromedex)). *E. coli* DH5α (New England Biolabs) was used as a cloning strain and *E. coli* K12 MG1655, a gift from Didier Mazel, was used for screening experiments.

## Phage strains and stocks

The MG1655::λ lysogeny was a gift from Luciano Marraffini, the MG1655::186cI-ts lysogen was a gift from Keith Shearwin and phage T4 was a gift from Laurent Debarbieux. Overnight cultures of MG1655::λ and C600::186cIts were harvested and the supernatant was filtered (0.22 µm). All the liquid stocks were further propagated in MG1655 grown in LB supplemented with maltose 0.2% (Sigma) and CaCl₂ 5 mM (Sigma) at a multiplicity of infection (MOI) of 1. Phage titers were measured by spotting 2 µL-drops of 10-fold serial dilutions in λ dilution buffer (TrisHCl 20 mM, NaCl 0.1M and MgSO₄ 10 mM) on a bacterial lawn of *E. coli* MG1655 in LB-agar (0.5%) containing 5 mM CaCl₂.

## *E. coli* strain construction

Strain LC-E75 used for the essentiality screen was described in our previous study [26]. This strain expresses an optimized *dcas9* cassette under the control of an aTc-inducible pTet promoter integrated at the phage 186 *attB* site. In order to study phage 186 with our screen, a new strain FR-E01 was constructed with the same cassette integrated at the HK022 *attB* site to avoid any interference. A high-copy pOSIP backbone [88] containing the HK022 integrase, HK022 *attB* site and a Kanamycin-resistance cassette was digested by EcoRI and PstI (New England Biolabs). A fragment containing *dcas9* under the control of pTet promoter was amplified by Phusion PCR (ThermoScientific) with primers LC124/LC125 to introduce homology regions with the backbone. The 2 fragments were assembled by the Gibson method [89]. The resulting vector was electroporated into strain MG1655. The backbone containing the HK022 integrase and a Kan<sup>R</sup> selection marker was then removed using the pE-FLP (Amp<sup>R</sup>) plasmid able to recombine FRT sites flanking the backbone [88]. pE-FLP was then cured through serial restreaks on LB plates, yielding strain FR-E01. Primers used for cloning are listed in S8 Table. Strain LC-E75 and FR-E01 are available on AddGene with the accession numbers 115925 and 118727 respectively.

## CRISPRi library design and assembly

A library of ~ 92,000 sgRNAs was designed previously [26]. These sgRNAs target 20-nt regions adjacent to NGG sites in *E. coli* K-12 MG1655 (NC_000913.2) and were chosen randomly among the total pool of possible sgRNAs in this strain. Briefly, the library was generated through on-chip oligo synthesis (CustomArray), amplified with primers LC296 and LC297 using Phusion DNA polymerase (ThermoScientific) and assembled into the psgRNAcos backbone using the Gibson method [89]. The resulting plasmid library DNA was transferred to strains LC-E75 and FR-E01 by electroporation yielding > 10⁸ colonies, thus ensuring a ~1000X coverage of the library. Primers used for cloning are listed in S8 Table.

## High-throughput screens

The data for the screen performed with strain LC-E75 grown in rich medium was obtained from our previous study [26]. This screen was performed over 17 generations in triplicates
from independent aliquots of the library generated from 3 independent transformations into
strain LC-E75.

The phage screen was performed in triplicates as follows: FR-E01 was grown at 37°C from 1 mL aliquots stored at -80 °C into 500 mL LB. At OD₆₀₀ = 0.2, dCas9 expression was induced by addition of 1 µM aTc (Acros Organics) to trigger the silencing of the target genes. The culture was grown to stationary phase (OD₆₀₀ = 2) and diluted 100-fold in LB containing 1 µM aTc, 0.2% Maltose and 5 mM CaCl₂. At OD₆₀₀ = 0.4, 20 mL of the culture was sampled and the library was extracted by miniprep (Nucleospin Plasmid, Macherey-Nage) to obtain the sgRNA distribution before infection. The culture was then infected with 1 mL of λ, T4 or 186cI-ts stocks (10⁷ pfu / µL) to reach a MOI of 1 ensuring a high infection rate while limiting double infections. After 2 h at 37°C, the cultures were harvested (7 min–4,000 g). Pellets were washed twice in PBS and the library was extracted by miniprep to obtain the sgRNA distribution after infection. For phage λ, the supernatant containing a mixture of λ and packaged library was collected and filtered (0.22 µm). Transduction was performed with 25 mL of lysate and 75 mL of stationary phase culture of MG1655 carrying λ prophage (providing resistance to super-infection) grown in LB supplemented with 0.2% maltose. After 30 min at 37°C, cells were harvested and resuspended in 1 mL LB and transduced cells were selected on 4 12x12 cm Petri dishes containing kanamycin. After 4h at 37°C, nascent colonies were washed and harvested in 5 mL LB-Kan. Each sample was split into 6 tubes from which library plasmids were extracted by miniprep (Nucleospin Plasmid, Macherey-Nage) before pooling back the 6 extractions into one sample for sequencing.

## Illumina sample preparation and sequencing

Library sequencing was performed as previously described [26]. Briefly, a customized Illumina sequencing method was designed to avoid problems arising from low library diversity when sequencing PCR products. Two Phusion PCR reactions (ThermoScientific) were used to generate the sequencing library with primers listed in S8 Table. The 1st PCR adds the 1st index. The 2nd PCR adds the 2nd index and flow cell attachment sequences. Sequencing is then performed using primer LC609 as a custom read 1 primer. Custom index primers were also used: LC499 reads index 1 and LC610 reads index 2. Sequencing was performed on a NextSeq 500 benchtop sequencer (Illumina). The first 2 cycles which read bases common to all clusters were set as dark cycles, followed by 20 cycles to read the guide. Using this strategy, we obtained on average 17 million reads per sample for the growth-based screen in LC-E75 and 4.6 million reads per sample for the phage screens in FR-E01 respectively.

## Data analysis

Raw sequencing files are available on the European Nucleotide Archive (https://www.ebi.ac.uk/ena) with the accession number PRJEB28256. Indexes were used to de-multiplex the data. Guides were filtered for potential off-target effect by discarding guides when the 9 PAM-proximal bases have a perfect match next to an NGG PAM in the promoter region of a gene (loosely defined as a window of 100bp before to 20bp after the start codon of all genes in the genome), or when the 11 PAM-proximal bases have a perfect match next to an NGG PAM allowing binding to the coding strand of a gene. Guides were also filtered to remove those carrying one of the 10 strongest bad seeds described in our previous study (AGGAA, TAGGA, ACCCA, TTGGA, TATAG, GAGGC, AAAGG, GGGAT, TAGAC and GTCCT) [26]. Guides with fewer than 20 reads in total were discarded.

Statistical analysis was performed from count data using the DESeq2 package [32] in R.
This package allows comparison of expression data by modeling read counts using a negative
binomial generalized linear model. The read counts of a non-targeting control guide RNA present in the library (5'-TGAGACCAGTCTAGGTCTCG-3') were used as the normalization factor. A paired analysis was performed to compare each sample to its initial condition. The log2FoldChange (log2FC) value represents the enrichment or depletion of each sgRNA. For the phage screen, guides targeting the template strand of genes or outside of genes were excluded from the analysis as well as well guides with insufficient number of reads (BaseMean < 10), yielding a library of ~ 17,200 sgRNAs. The lists of all sgRNAs with computed log2FC values after the growth-based screen, after phage screens and after transduction assay are provided as S1, S4 and S6 Tables respectively. For each gene, the median log2FC value of the sgRNAs targeting the coding strand was used for ranking (S2 and S5 Tables). For the λ transduction experiment, a linear model was built for each gene to predict log2FC using the log2FC values of the essentiality screen and a Boolean factor considering the gene, or not, as a host factor. A given gene was regarded as a hit when the Boolean factor significantly improved the model after correction for multiple testing (ANOVA, FDR < 0.05), showing that the log2FC values obtained for this gene cannot be explained only by the fitness defect induced by the guides. The regression coefficient of this parameter termed as "estimate" was used for gene ranking. Genes whose silencing decreases the production of functional capsids have a negative estimate after transduction. A list of estimates and FDR for each gene is provided as S7 Table. Transcriptional units (operons) were obtained from RegulonDB [90]. The log2FC values after the essentiality screen can be visualized on our web-based tool CRISPRbrowser (crispr.pasteur.fr). The code used for data analysis is available at gitlab.pasteur.fr/dbikard/dCas9_genome_wide_screen.

## Infection dynamics

sgRNAs were cloned into a psgRNAcos backbone through Golden Gate assembly using BsaI [91] and were electroporated into strain FR-E01. For essential gene *yrfF*, 4 mismatches were inserted at the 5'-end of the sgRNA to decrease expression of the target gene to intermediate levels. Constructions were validated by Sanger sequencing. The list of individual sgRNAs is provided in S9 Table. Strains were grown overnight and diluted 100-fold in LB medium containing 0.2% maltose, 1 µM aTc, 5 mM CaCl₂ and kanamycin. At OD₆₀₀ = 0.4, 90 µL of cultures were mixed in 96-well plates with 100 µL of fresh medium and 10 µL of a 8.10⁴ pfu/µL stock of the appropriate phage (MOI ~ 0.01). Infection dynamics were monitored in three replicates on an Infinite M200Pro (Tecan) at 37°C with shaking for 8h. Strains Δ*waaJ* and Δ*rcsB* were obtained from the Keio collection. The *rcsB* deletion was transferred to strain FR-E01 by P1 transduction.

## Gene deletions

Plasmid pKOBEG-A [92] carrying the λ-red system components under the control of an arabinose-inducible promoter was transformed into strain MG1655. A kanamycin resistance cassette was amplified by PCR with Phusion polymerase (ThermoScientific) from plasmid pKD4 with primers introducing homology regions with sequences flanking genes *alsK*, *bcsB* or *entD*. Electrocompetent MG1655-pKOBEG-A cells were prepared with arabinose (Sigma) before transformation of 1 µg of the DNA fragments. Recovery and overnight incubation were performed at 30°C with arabinose and kanamycin. The next day, colonies were restreaked and incubated at 37°C. Primers used for pKD4 amplification and colony screening are provided in S8 Table.
RT-qPCR

Overnight cultures were diluted 1:100 in 3 mL LB containing 1 µM aTc. Cells were further grown for 2h before RNA extraction using Direct-zol RNA MiniPrep (Zymo Research) followed by DNAsse treatment using TURBO DNA-free Kit (Thermo Fisher Scientific). All RNAs were reverse transcribed into cDNA using the Transcription First Strand cDNA Synthesis Kit (Roche) using 500 ng RNA. qPCR was performed with the FastStart Essential DNA Green master mix (Roche) in a LightCycler 96 (Roche) following the manufacturer's instructions. qPCR was performed in two technical replicates and three biological replicates. Relative gene expression was computed on LightCycler 96 software (Roche) using the ΔΔCq method [93] after normalization by 5S rRNA (rrsA). qPCR primers are listed in S10 Table.

**Supporting information**

**S1 Fig.** Target regions of guides that produce a strong fitness defect while targeting multiple positions. (A) 348/1932 sgRNAs that simultaneously target several positions induce a fitness defect (log2FC < -2). These guides mostly target rRNAs, tRNAs, repeat regions or elongation factor genes (tufA-tufB). (B) Among sgRNAs targeting REP elements, sgRNAs simultaneously targeting several regions have more chance of inducing a fitness defect than sgRNAs targeting a single region (Fisher's exact test, odds ratio = 0.33, p < 2.10-5). (C) Example of a sgRNA simultaneously targeting 3 repeat regions. sgRNAs targeting the +1 or -1 strand are dotted in red or blue respectively. Repeat regions are represented as grey boxes.

(PNG)

**S2 Fig.** Gene interaction network of the top 100 genes having the lowest median log2Fold-Change. The STRING database was used to compute a gene interaction network [94]. Genes were colored by function. Line thickness indicates confidence of the interaction.

(PNG)

**S3 Fig.** sgRNAs with off-targets in essential genes induce a strong fitness defect. (A) A sgRNA targeting *mdfA* has a 10-nt perfect match to *rpsL*. (B) A sgRNA targeting *flgC* has a 10-nt perfect match to *lptG*. (A,B) Matched base pairs are shown in blue. On the plots, sgRNAs targeting the +1 strand or the -1 strand are shown in red or blue respectively. A black arrow indicates the sgRNA with an off-target activity. A dashed line indicates the off-target position.

(PNG)

**S4 Fig.** Predicted internal promoters within operons. We identified 7 operons in which an expected polar effect is not observed, i.e guides targeting a gene upstream of an essential or near-essential gene are not depleted. This suggests that the downstream gene can be expressed from an internal promoter. Promoters predicted from a recent transcription start site dataset are shown as dashed arrow [38].

(PNG)

**S5 Fig.** Replacement of genes *alsK*, *bcsB* and *entD* by a Kanamycin resistance cassette through λ-red recombination shows that these genes are not essential. The kanamycin resistance cassette from plasmid pKD4 was amplified with primers designed to introduce 50 bp-long homologies with regions flanking genes *alsK*, *bcsB* and *entD*. Primer couples P1 + P2, P1' + P2' and P1" + P2" were designed to flank the genetic region of *alsK*, *bcsB* and *entD* respectively, while primer couples P3 + P4, P3' + P4' and P3" + P4" were designed to amplify within *alsK*, *bcsB* and *entD* respectively to demonstrate that gene duplication did not occur
during the experiment.

(PNG)

**S6 Fig.** *Inefficient dCas9-mediated repression of genes regulated by negative feedback loops.* *lexA* and *rho* are two well-known essential genes whose product inhibits their own expression and which are classified as nonessential in our screen. (A) Relative *lexA* or *rho* expression was measured in presence of a *lexA*- or *rho*-targeted sgRNA respectively, with or without dCas9 repression (± aTc), showing a low repression activity (66.4% and 85.2% respectively). RT-qPCR results are shown for 3 biological replicates and 2 technical replicates. (B) To measure the activity of the *lexA* and *rho* promoters while targeting the respective gene with dCas9, we built plasmids pFR42 which expresses sfGFP from the *lexA* promoter and pFR43 which expresses sfGFP from the *rho* promoter. These plasmids also express the corresponding sgRNA constitutively. (C) pFR42 and pFR43 were transformed into strain LC-E75 used in the screen. An overnight culture was diluted 100-fold with or without aTc and OD600 and GFP fluorescence were measured over time. Raw GFP fluorescence after 12h was normalized by OD600 and the normalized fluorescence of the control (LC-E75 with psgRNAcos) was subtracted. Bar plot shows mean ± standard deviation (n = 3).

(PNG)

**S7 Fig.** *Unexpected fitness defects shed light on atypical gene organizations.* (A) *yejL* is located upstream of the sensitive essential gene *yejM*. (B) *ypaB* contains the promoter of the essential gene *nrdA*. (C) *sroG* is a riboswitch controlling the expression of the essential gene *ribB*. (D) *yceQ* is annotated as essential but is actually located upstream *rne* in the opposite direction and contains the promoter driving the expression of *rne*. (E) *psrO* encodes a small RNA and is located in the promoter region of near-essential gene *pnp*. sgRNAs targeting the +1 or -1 strand are shown as red or blue dots respectively.

(PNG)

**S8 Fig.** *Strains LC-E75 and FR-E01 produce similar concentrations of dCas9.* Strains LC-E75 and FR-E01 were grown overnight and diluted 1:100 with aTc for 2 h before harvesting. Samples were run in NuPAGE Novex Bis-Tris gels in reducing condition before transfer to PVDF membranes. Rabbit monoclonal CRISPR-Cas9 antibody and rabbit polyclonal RecA antibody were used.

(PNG)

**S9 Fig.** *Phage λ can package plasmid psgRNAcos upon infection.* (A) MG1655 carrying psgRNAcos was infected with λ at MOI = 1. The lysate containing a mix of phage and packaged cosmid was extracted after 2 h and the relative concentrations of phage and cosmid were measured by plaque assay and by transduction into strain MG1655::λ respectively. This strain carries the λ lysogen and is thus resistant to superinfection by the λ particles present in the lysate. (B,C) Bar plot and dot plot show mean ± standard deviation (n = 9).

(PNG)

**S10 Fig.** *sgRNAs in the library with a perfect match in λ phage genome.* sgRNAs targeting the positive or negative strand are dotted in red or blue respectively. Genes colored in yellow, orange and red respectively correspond to early right, early left and late operon.

(PNG)

**S1 Table.** *List of sgRNAs with log2FC values after the growth-based screen in strain LC-E75.* Each sgRNA is indexed with its position, orientation (ori), target strand (coding) and computed fold change (log2FC, padj and gamma). Information on the targeted gene is also
provided: name (gene), essentiality (essential), position (gene_left, gene_right) and orientation (gene_ori).
(CSV)

**S2 Table. List of genes with median log2FC values after growth-based screen.** The following information is provided. Gene name (gene), essentiality (essential), orientation (gene_ori), position (gene_left, gene_right), median log2FC for each strand (median_coding and median_template) and median absolute deviation (mad_coding and mad_template), number of guides targeting the coding or template strand (coding and template), as well as information on operon structure (operon). This table contains all genes targeted in the screen including protein-coding and RNAs.
(CSV)

**S3 Table. List of essential genes not found essential in our screen.**
(XLSX)

**S4 Table. List of sgRNAs with log2FC values after phage screens.** Each sgRNA is indexed with its position, orientation (ori), strand (coding), information on the targeted gene (gene, essential, gene_left, gene_right, gene_ori) and computed fold change for each phage (log2FCent λ, log2FC_T4 and log2FC_186).
(CSV)

**S5 Table. List of genes with resistance scores after phage screens.** Each gene is indexed with its information (gene, essential), median log2FC (Resistance_Scoreent λ, Resistance_Score_T4 and Resistance_Score_186) and median absolute deviation of guides targeting the coding strand for each phage, number of guides targeting the coding strand used for calculations (sgRNAs), as well as information on operon structures (operon).
(CSV)

**S6 Table. List of sgRNAs with log2FC values after transduction screen.** Each sgRNA is indexed with its position, orientation (ori), strand (coding), information on the targeted gene (gene, essential, gene_left, gene_right, gene_ori) and computed fold change after transduction (log2FC).
(CSV)

**S7 Table. List of genes with estimates and FDR after transduction screen.** Each gene is indexed with its information (gene, essential), operon information (operon) and number of guides targeting the coding strand used for calculations (sgRNAs). For each gene, a linear model was built to take the effect on cell fitness into account (see Methods). The regression coefficient of the Boolean parameter termed as "estimate" is used for gene ranking. Genes whose silencing decreases the production of functional capsids have a negative estimate after transduction. P-value and FDR associated to this parameter are also displayed for each gene.
(CSV)

**S8 Table. List of primers used for cloning and sequencing.**
(DOCX)

**S9 Table. List of individual sgRNAs used in this study.**
(DOCX)

**S10 Table. List of primers used for qPCR.**
(DOCX)
**Acknowledgments**

We thank Alicia Calvo-Villamañán and Olaya Rendueles-Garcia for helpful advice.

**Author Contributions**

**Conceptualization:** François Rousset, Lun Cui, Elise Siouve, David Bikard.

**Data curation:** François Rousset, Lun Cui.

**Formal analysis:** François Rousset, Lun Cui, David Bikard.

**Funding acquisition:** David Bikard.

**Investigation:** François Rousset, Lun Cui, Elise Siouve, Florence Depardieu.

**Methodology:** François Rousset, Lun Cui.

**Project administration:** David Bikard.

**Software:** Christophe Becavin.

**Supervision:** David Bikard.

**Visualization:** François Rousset, Christophe Becavin.

**Writing – original draft:** François Rousset, David Bikard.

**References**

1. Brouns SJJ, Jore MM, Lundgren M, Westra ER, Slijkhuis RJH, Snijders APL, et al. Small CRISPR RNAs Guide Antiviral Defense in Prokaryotes. Science (80-). 2008; 321: 960–964. https://doi.org/10.1126/science.1159689 PMID: 18703739

2. Barrangou R, Fremaux C, Deveau H, Richards M, Boyaval P, Moineau S, et al. CRISPR Provides Acquired Resistance Against Viruses in Prokaryotes. Science (80-). 2007; 315: 1709–1712. https://doi.org/10.1126/science.1138140 PMID: 17379808

3. Garneau JE, Dupuis M-È, Villion M, Romero DA, Barrangou R, Boyaval P, et al. The CRISPR/Cas bacterial immune system cleaves bacteriophage and plasmid DNA. Nature. Nature Research; 2010; 468: 67–71. https://doi.org/10.1038/nature09523 PMID: 21048762

4. Barrangou R, Doudna JA. Applications of CRISPR technologies in research and beyond. Nat Biotechnol. 2016; 34: 933–941. https://doi.org/10.1038/nbt.3659 PMID: 27606440

5. Jinek M, Chylinski K, Fonfara I, Hauer M, Doudna JA, Charpentier E, et al. A programmable dual-RNA-guided DNA endonuclease in adaptive bacterial immunity. Science. American Association for the Advancement of Science; 2012; 337: 816–21. https://doi.org/10.1126/science.1225829 PMID: 22745249

6. Shalem O, Sanjana NE, Hartenian E, Shi X, Scott DA, Mikkelson T, et al. Genome-scale CRISPR-Cas9 knockout screening in human cells. Science. NIH Public Access; 2014; 343: 84–87. https://doi.org/10.1126/science.1247005 PMID: 24336571

7. Wang T, Wei JJ, Sabatini DM, Lander ES. Genetic screens in human cells using the CRISPR-Cas9 system. Science. NIH Public Access; 2014; 343: 80–4. https://doi.org/10.1126/science.1246981 PMID: 24336569

8. Wang T, Birsoy K, Hughes NW, Krupczak KM, Post Y, Wei JJ, et al. Identification and characterization of essential genes in the human genome. Science (80-). American Association for the Advancement of Science; 2015; 350: 1096–1101. https://doi.org/10.1126/SCIENCE.AAC7041 PMID: 26472758

9. Bertomeu T, Coulombe-Huntington J, Chatr-Aryamontri A, Bourdages KG, Coyaud E, Raught B, et al. A High-Resolution Genome-Wide CRISPR/Cas9 Viability Screen Reveals Structural Features and Contextual Diversity of the Human Cell-Essential Proteome. Mol Cell Biol. American Society for Microbiology; 2018; 38: e00302–17. https://doi.org/10.1128/MCB.00302-17 PMID: 29038160

10. Evers B, Jastrzebski K, Heijmans JPM, Grenrum W, Beijersbergen RL, Bernards R. CRISPR knockout screening outperforms shRNA and CRISPRi in identifying essential genes. Nat Biotechnol. Nature Publishing Group; 2016; 34: 631–633. https://doi.org/10.1038/nbt.3536 PMID: 27111720
11. Morgens DW, Deans RM, Li A, Bassik MC. Systematic comparison of CRISPR/Cas9 and RNAi screens for essential genes. Nat Biotechnol. Nature Publishing Group; 2016; 34: 634–636. https://doi.org/10.1038/nbt.3567 PMID: 27159373

12. Hart T, Chandrashekhar M, Aregger M, Steinhart Z, Brown KR, MacLeod G, et al. High-Resolution CRISPR Screens Reveal Fitness Genes and Genotype-Specific Cancer Liabilities. Cell. 2015; 163: 1515–1526. https://doi.org/10.1016/j.cell.2015.11.015 PMID: 26627737

13. Ma H, Dang Y, Wu Y, Jia G, Anaya E, Zhang J, et al. A CRISPR-based screen identifies genes essential for west-nile-virus-induced cell death. Cell Rep. The Authors; 2015; 12: 673–683. https://doi.org/10.1016/j.celrep.2015.06.049 PMID: 26190106

14. Zhang R, Miner JJ, Gorman MJ, Rausch K, Ramage H, White JP, et al. A CRISPR screen defines a sig- nal peptide processing pathway required by flaviviruses. Nature. Nature Publishing Group; 2016; 535: 164–168. https://doi.org/10.1038/nature18625 PMID: 27383988

15. Park RJ, Wang T, Koundakjian D, Hultquist JF, Lamothe-Molina P, Monel B, et al. A genome-wide CRISPR screen identifies a restricted set of HIV host dependency factors. Nat Genet. Nature Publishing Group; 2017; 49: 193–203. https://doi.org/10.1038/ng.3741 PMID: 27992415

16. Bikard D, Euler CW, Jiang W, Nussenzweig PM, Goldberg GW, Duportet X, et al. Exploiting CRISPR-Cas nucleases to produce sequence-specific antimicrobials. Nat Biotechnol. Nature Publishing Group; 2014; 32: 1146–50. https://doi.org/10.1038/nbt.3043 PMID: 25282355

17. Gomaa AA, Klumpe HE, Luo ML, Selle K, Barrangou R, Beisel CL. Programmable removal of bacterial strains by use of genome-targeting CRISPR-Cas systems. MBio. American Society for Microbiology; 2014; 5: e00928–13. https://doi.org/10.1128/mBio.00928-13 PMID: 24473129

18. Cui L, Bikard D. Consequences of Cas9 cleavage in the chromosome of Escherichia coli. Nucleic Acids Res. 2016; 44: 4243–4251. https://doi.org/10.1093/nar/gkw223 PMID: 27060147

19. Bikard D, Hatoum-Aslan A, Mucida D, Marraffini LA. CRISPR Interference Can Prevent Natural Trans- formation and Virulence Acquisition during In Vivo Bacterial Infection. Cell Host Microbe. 2012; 12: 177–186. https://doi.org/10.1016/j.chom.2012.06.003 PMID: 22901538

20. Citorik RJ, Mimee M, Lu TK. Sequence-specific antimicrobials using efficiently delivered RNA-guided nucleases. Nat Biotechnol. Nature Research; 2014; 32: 1141–1145. https://doi.org/10.1038/nbt.3011 PMID: 25240928

21. Bikard D, Jiang W, Samai P, Hochschild A, Zhang F, Marraffini LA. Programmable repression and acti- vation of bacterial gene expression using an engineered CRISPR-Cas system. Nucleic Acids Res. 2013; 41: 7429–7437. https://doi.org/10.1093/nar/gkt520 PMID: 23761437

22. Qi LS, Larson MH, Gilbert LA, Doudna JA, Weissman JS, Arkin AP, et al. Repurposing CRISPR as an RNA-Guided Platform for Sequence-Specific Control of Gene Expression. Cell. 2013; 152: 1173–1183. https://doi.org/10.1016/j.cell.2013.02.022 PMID: 23452860

23. Peters JM, Colavin A, Shi H, Czarny TL, Larson MH, Wong S, et al. A Comprehensive, CRISPR-based Functional Analysis of Essential Genes in Bacteria. Cell. 2016; 165: 1493–1506. https://doi.org/10.1016/j.cell.2016.05.003 PMID: 27238023

24. Liu X, Gallay C, Kjos M, Domenech A, Slager J, van Kessel SP, et al. High-throughput CRISPRi pheno- typing identifies new essential genes in Streptococcus pneumoniae. Mol Syst Biol. EMBO Press; 2017; 13: 931. https://doi.org/10.15252/msb.20167449 PMID: 28490437

25. Wang T, Guan C, Guo J, Liu B, Wu Y, Xie Z, et al. Pooled CRISPR interference screening enables genome-scale functional genomics study in bacteria with superior performance. Nat Commun. Nature Publishing Group; 2018; 9: 2475. https://doi.org/10.1038/s41467-018-04899-x PMID: 29946130

26. Cui L, Vigouroux A, Rousset F, Varet H, Khanna V, Bikard D. A CRISPRi screen in E. coli reveals sequence-specific toxicity of dCas9. Nat Commun. Nature Publishing Group; 2018; 9: 1912. https://doi.org/10.1038/s41467-018-04209-5 PMID: 29765036

27. Henry M, Debarbieux L. Tools from viruses: Bacteriophage successes and beyond. Virology. 2012; 434: 151–161. https://doi.org/10.1016/j.virol.2012.09.017 PMID: 23063405

28. Keen EC. A century of phage research: bacteriophages and the shaping of modern biology. Bioessays. NIH Public Access; 2015; 37: 6–9. https://doi.org/10.1002/bies.201400152 PMID: 25521633

29. Maynard ND, Birch EW, Sanghvi JC, Chen L, Gutschow M V., Covert MW. A Forward-Genetic Screen and Dynamic Analysis of Lambda Phage Host-Dependencies Reveals an Extensive Interaction Net- work and a New Anti-Viral Strategy. Matic I, editor. PLoS Genet. Public Library of Science; 2010; 6: e1001017. https://doi.org/10.1371/journal.pgen.1001017 PMID: 20628568

30. Qimron U, Marintcheva B, Tabor S, Richardson CC. Genomewide screens for Escherichia coli genes affecting growth of T7 bacteriophage. Proc Natl Acad Sci U S A. 2006; 103: 19039–44. https://doi.org/10.1073/pnas.0609428103 PMID: 17135349
31. Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, Baba M, et al. Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol. EMBO Press; 2006; 2: 2006.0008. https://doi.org/10.1038/msb4100050 PMID: 16738554

32. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15.

33. Liang W, Rudd KE, Deutscher MP. A role for REP sequences in regulating translation. Mol Cell. Elsevier; 2015; 58: 431-9. https://doi.org/10.1016/j.molcel.2015.03.019 PMID: 25891074

34. Zhou J, Rudd KE. EcoGene 3.0. Nucleic Acids Res. Oxford University Press; 2012; 41: D613-D624. https://doi.org/10.1093/nar/gks1235 PMID: 23197660

35. Goodall ECA, Robinson A, Johnston IG, Jabbari S, Turner KA, Cunningham AF, et al. The Essential Genome of Escherichia coli K-12. MBio. American Society for Microbiology; 2018; 9: e02096-17. https://doi.org/10.1128/mBio.02096-17 PMID: 29463657

36. Keseler IM, Mackie A, Santos-Zavaleta A, Billington R, Bonavides-Martínez C, Caspi R, et al. The Eco-Cyc database: reflecting new knowledge about *Escherichia coli* K-12. Nucleic Acids Res. Oxford University Press; 2017; 45: D543-D550. https://doi.org/10.1093/nar/gkw1003 PMID: 27899573

37. Hansen S, Vulić M, Min J, Yen T-J, Schumacher MA, Brennan RG, et al. Regulation of the Escherichia coli HipBA Toxin-Antitoxin System by Proteolysis. de Crécy-Lagard V, editor. PLoS One. Public Library of Science; 2012; 7: e39185. https://doi.org/10.1371/journal.pone.0039185 PMID: 22720069

38. Thomason MK, Bischler T, Eisenbart SK, Förstner KU, Zhang A, Herbig A, et al. Global transcriptional start site mapping using differential RNA sequencing reveals novel antisense RNAs in Escherichia coli. J Bacteriol. American Society for Microbiology; 2015; 197: 18-28. https://doi.org/10.1128/JB.02096-14 PMID: 25266388

39. Gerdes SY, Scholle MD, Campbell JW, Balázsi G, Ravasz E, Daugherty MD, et al. Experimental determination and system level analysis of essential genes in Escherichia coli MG1655. J Bacteriol. American Society for Microbiology; 2003; 185: 5673-84. https://doi.org/10.1128/JB.185.19.5673-5684.2003

40. Yamazaki Y, Niki H, Kato J. Profiling of Escherichia coli Chromosome Database. Methods in molecular biology (Clifton, NJ). 2008. pp. 385-389. https://doi.org/10.1007/978-1-59745-321-9_26 PMID: 18392982

41. Yamamoto N, Nakahigashi K, Nakamichi T, Yoshino M, Takai Y, Touda Y, et al. Update on the Keio collection of Escherichia coli single-gene deletion mutants. Mol Syst Biol. EMBO Press; 2009; 5: 335. https://doi.org/10.1038/msb.2009.92 PMID: 20029369

42. Vigouroux A, Oldewurtel E, Cui L, Bikard D, van Teeffelen S. Tuning dCas9's ability to block transcription enables robust, noiseless knockdown of bacterial genes. Mol Syst Biol. EMBO Press; 2018; 14: e7899-e7899. https://doi.org/10.15252/msb.20177899 PMID: 29519933

43. Richter G, Volk R, Krieger C, Lahm HW, Röthlisberger U, Bacher A. Biosynthesis of riboflavin: cloning, sequencing, and expression of the gene coding for 3,4-dihydroxy-2-butanone 4-phosphate synthase of Escherichia coli. J Bacteriol. 1992; 174: 4050-6. PMID: 1597419

44. Randall-Hazelbauer L, Schwartz M. Isolation of the bacteriophage lambda receptor from *Escherichia coli*. J Bacteriol. 1973; 116: 1436-1446. PMID: 4201774

45. Yu F, Mizushima S. Roles of lipopolysaccharide and outer membrane protein OmpC of Escherichia coli K-12 in the receptor function for bacteriophage T4. J Bacteriol. American Society for Microbiology; 1982; 151: 718-22.

46. Elliott J, Arber W. E. coli K-12 pel mutants, which block phage λ DNA injection, coincide with ptsM, which determines a component of a sugar transport system. MGG Mol Gen Genet. Springer-Verlag; 1978; 161: 1-8. https://doi.org/10.1007/BF00266608

47. Randall LL. Quantitation of the loss of the bacteriophage lambda receptor protein from the outer membrane of lipopolysaccharide-deficient strains of Escherichia coli. J Bacteriol. American Society for Microbiology; 1975; 123: 41-6.

48. Clarke DJ. The Rcs phosphorelay: more than just a two-component pathway. Future Microbiol. 2010; 5: 1173-84. https://doi.org/10.2217/fmb.10.83 PMID: 20722597

49. Gottesman S, Stout V. Regulation of capsular polysaccharide synthesis in Escherichia coli K12. Mol Microbiol. Blackwell Publishing Ltd; 1991; 5: 1599-1606. https://doi.org/10.1111/j.1365-2958.1991.tb01906.x

50. Stout V. Regulation of capsule synthesis includes interactions of the RcsC/RcsB regulatory pair. Res Microbiol. Elsevier Masson; 1994; 145: 389-392. https://doi.org/10.1016/0923-2508(94)90086-8

51. Rendueles O, Garcia-Garcerà M, Néron B, Touchon M, Rocha EPC. Abundance and co-occurrence of extracellular capsules increase environmental breadth: Implications for the emergence of pathogens. Bessen DE, editor. PLOS Pathog. Public Library of Science; 2017; 13: e1006525. https://doi.org/10.1371/journal.ppat.1006525 PMID: 28742161
52. Kim MS, Kim YD, Hong SS, Park K, Ko KS, Myung H. Phage-encoded colanic acid-degrading enzyme permits lytic phage infection of a capsule-forming resistant mutant Escherichia coli strain. Appl Environ Microbiol. American Society for Microbiology; 2015; 81: 900-9. https://doi.org/10.1128/AEM.02606-14 PMID: 25416767

53. Scholl D, Adhya S, Merril C. Escherichia coli K1's capsule is a barrier to bacteriophage T7. Appl Environ Microbiol. American Society for Microbiology; 2005; 71: 4872-4. https://doi.org/10.1128/AEM.71.8.4872-4874.2005

54. Paynter MJB, Bungay HR. Capsular protection against virulent coliphage infection. Biotechnol Bioeng. Wiley Subscription Services, Inc., A Wiley Company; 1970; 12: 341-346. https://doi.org/10.1002/bit.260120303 PMID: 4920047

55. Joloba ML, Clemmer KM, Sledjeski DD, Rather PN. Activation of the gab operon in an RpoS-dependent manner by mutations that truncate the inner core of lipopolysaccharide in Escherichia coli. J Bacteriol. American Society for Microbiology; 2004; 186: 8542-6. https://doi.org/10.1128/JB.186.24.8542-8546.2004

56. Belfort M, Wulff DL. A Mutant of Escherichia coli That Is Lysogenized with High Frequency. Cold Spring Harb Monogr Arch. 1971; 739-742.

57. Kihara A, Akiyama Y, Ito K. Revisiting the Lysogenization Control of Bacteriophage λ. J Biol Chem. 2001; 276: 13695-13700. https://doi.org/10.1074/jbc.M011699200 PMID: 11278968

58. Yochem J, Uchida H, Sunshine M, Saito H, Georgopoulos C, Feiss M. Genetic Analysis of Two Genes, dnaJ and dnaK, Necessary for Escherichia coli and Bacteriophage Lambda DNA Replication. Mol Genet. 1978; 164: 9-14. PMID: 360041

59. Maynard ND, Macklin DN, Kirkegaard K, Covert MW. Competing pathways control host resistance to virus via tRNA modification and programmed ribosomal frameshifting. Mol Syst Biol. EMBO Press; 2012; 8: 567. https://doi.org/10.1038/msb.2011.101 PMID: 22294093

60. Takahashi N, Kobayashi I. Evidence for the double-strand break repair model of bacteriophage lambda recombination. Proc Natl Acad Sci U S A. 1990; 87: 2790-4. PMID: 2138786

61. Cronan JE. Improved plasmid-based system for fully regulated off-to-on gene expression in Escherichia coli: Application to production of toxic proteins. Plasmid. Elsevier Inc.; 2013; 69: 81-89. https://doi.org/10.1016/j.plasmid.2012.09.003 PMID: 23022297

62. Li J, Horwitz R, McCracken S, Greenblatt J. NusG, a new Escherichia coli elongation factor involved in transcriptional antitermination by the N protein of phage λ. J Biol Chem. 1992; 267: 6012-6019. PMID: 1532577

63. Gottesman ME, Adhya S, Das A. Transcription antitermination by bacteriophage lambda N gene product. J Mol Biol. Academic Press; 1980; 140: 57-75. https://doi.org/10.1016/0022-2836(80)90356-3

64. Friedman DI. Interaction between bacteriophage λ and its Escherichia coli host. Curr Opin Genet Dev. Elsevier Current Trends; 1992; 2: 727-738. https://doi.org/10.1016/S0959-437X(05)80133-9

65. Ghysen A, Pironio M. Relationship between the N function of bacteriophage λ and host RNA polymerase. J Mol Biol. Academic Press; 1972; 65: 259-272. https://doi.org/10.1016/0022-2836(72)90281-1

66. Gottesman ME, Weisberg RA. Little lambda, who made thee? Microbiol Mol Biol Rev. American Society for Microbiology; 2004; 68: 796-813. https://doi.org/10.1128/MMBR.68.4.796-813.2004

67. Sternberg N. Properties of a mutant of Escherichia coli defective in bacteriophage λ head formation (groE). I. Initial characterization. J Mol Biol. 1973; 76: 1-23. https://doi.org/10.1016/0022-2836(73)90078-8 PMID: 4578098

68. Georgopoulos CP, Hendrix RW, Casjens SR, Kaiser AD. Host participation in bacteriophage lambda head assembly. J Mol Biol. 1973; 76: 45-60. https://doi.org/10.1016/0022-2836(73)90080-6 PMID: 4578100

69. Sternberg N. Properties of a mutant of Escherichia coli defective in bacteriophage lambda head formation (groE). II. The propagation of phage lambda. J Mol Biol. 1973; 76: 25-44. https://doi.org/10.1016/0022-2836(73)90079-X PMID: 4578099

70. Chuang SE, Blattner FR. Characterization of twenty-six new heat shock genes of Escherichia coli. J Bacteriol. 1993; 175: 5242-52. PMID: 8349564

71. Kim HK, Min S, Song M, Jung S, Choi JW, Kim Y, et al. Deep learning improves prediction of CRISPR-Cpf1 guide RNA activity. Nat Biotechnol. Nature Publishing Group; 2018; 36: 239-241. https://doi.org/10.1038/nbt.4061 PMID: 29431740

72. Chuai G, Wang Q-L, Liu Q. In Silico Meets In Vivo: Towards Computational CRISPR-Based sgRNA Design. Trends Biotechnol. Elsevier Current Trends; 2017; 35: 12-21. https://doi.org/10.1016/J.TIBTECH.2016.06.008 PMID: 27418421
73. Doench JG, Fusi N, Sullender M, Hegde M, Vaimberg EW, Donovan KF, et al. Optimized sgRNA design to maximize activity and minimize off-target effects of CRISPR-Cas9. Nat Biotechnol. Nature Publishing Group; 2016; 34: 1–12. https://doi.org/10.1038/nbt.3437 PMID: 26780180

74. Chari R, Yeo NC, Chavez A, Church GM. sgRNA Scorer 2.0: A Species-Independent Model To Predict CRISPR/Cas9 Activity. ACS Synth Biol. American Chemical Society; 2017; 6: 902–904. https://doi.org/10.1021/acssynbio.6b00343 PMID: 28146356

75. Labuhn M, Adams FF, Ng M, Knoess S, Schambach A, Charpentier EM, et al. Refined sgRNA efficacy prediction improves large- and small-scale CRISPR-Cas9 applications. Nucleic Acids Res. Oxford University Press; 2018; 46: 1375–1385. https://doi.org/10.1093/nar/gkx1268 PMID: 29267886

76. Price MN, Wetmore KM, Waters RJ, Callaghan M, Ray J, Liu H, et al. Mutant phenotypes for thousands of bacterial genes of unknown function. Nature. Nature Publishing Group; 2018; 1. https://doi.org/10.1038/s41586-018-0124-0 PMID: 29769716

77. Gallagher LA, Shendure J, Manoil C. Genome-scale identification of resistance functions in Pseudomonas aeruginosa using Tn-seq. MBio. American Society for Microbiology; 2011; 2: e00315–10. https://doi.org/10.1128/mBio.00315-10 PMID: 21253457

78. Barquist L, Mayho M, Cummins C, Cain AK, Boinett CJ, Page AJ, et al. The TraDIS toolkit: sequencing and analysis for dense transposon mutant libraries. Bioinformatics. Oxford University Press; 2016; 32: 1109–1111. https://doi.org/10.1093/bioinformatics/btw022 PMID: 26794317

79. Zomer A, Burghout P, Bootsma HJ, Hermans PWM, van Hijum SAFT. ESSENTIALS: Software for Rapid Analysis of High Throughput Transposon Insertion Sequencing Data. Provart NJ, editor. PLoS One. Public Library of Science; 2012; 7: e43012. https://doi.org/10.1371/journal.pone.0043012 PMID: 22900082

80. Solaimanpour S, Sarmiento F, Mrázek J. Tn-Seq Explorer: A Tool for Analysis of High-Throughput Sequencing Data of Transposon Mutant Libraries. Antoniewski C, editor. PLoS One. Public Library of Science; 2015; 10: e0126070. https://doi.org/10.1371/journal.pone.0126070 PMID: 25938432

81. Miller ES, Kutter E, Mosig G, Arisaka F, Kunisawa T, Rüger W. Bacteriophage T4 genome. Microbiol Mol Biol Rev. American Society for Microbiology; 2003; 67: 86–156, table of contents. https://doi.org/10.1128/MMBR.67.1.86-156.2003

82. Cowley LA, Low AS, Pickard D, Boinett CJ, Dallman TJ, Day M, et al. Transposon Insertion Sequencing Elucidates Novel Gene Involvement in Susceptibility and Resistance to Phages T4 and T7 in Escherichia coli O157. MBio. American Society for Microbiology; 2018; 9: e00705–18. https://doi.org/10.1128/mBio.00705-18 PMID: 30042196

83. Lenski RE. Experimental studies of pleiotropy and epistasis in Escherichia coli. I. Variation in competitive fitness among mutants resistant to virus T4. Evolution (N Y). 1988; 42: 425–432. https://doi.org/10.1111/j.1558-5646.1988.tb04149.x PMID: 28564005

84. Warren RJ, Bose SK. Bacteriophage-induced inhibition of host functions. I. Degradation of Escherichia coli deoxyribonucleic acid after T4 infection. J Virol. American Society for Microbiology Journals; 1968; 2: 327–34.

85. Bull JJ, Vimr ER, Molineux J. A tale of tails: Sialidase is key to success in a model of phage therapy against K1-capsulated Escherichia coli. Virology. Academic Press; 2010; 398: 79–86. https://doi.org/10.1016/J.VIROL.2009.11.040 PMID: 20006993

86. Cornelissen A, Ceyssens P-J, Krylov VN, Noben J-P, Volckaert G, Lavigne R. Identification of EPS-degrading activity within the tail spikes of the novel Pseudomonas putida phage AF. Virology. Academic Press; 2012; 434: 251–256. https://doi.org/10.1016/J.VIROL.2012.09.030 PMID: 23084421

87. Scholl D, Rogers S, Adhya S, Merril CR. Bacteriophage K1-5 encodes two different tail fiber proteins, allowing it to infect and replicate on both K1 and K5 strains of Escherichia coli. J Virol. American Society for Microbiology (ASM); 2001; 75: 2509–15. https://doi.org/10.1128/JVI.75.6.2509-2515.2001

88. St-Pierre F, Cui L, Priest DG, Endy D, Dodd IB, Shearwin KE. One-step cloning and chromosomal integration of DNA. ACS Synth Biol. 2013; 2: 537–541. https://doi.org/10.1021/sb400021j PMID: 24050148

89. Gibson DG, Young L, Chuang R-Y, Venter JC, Hutchison CA, Smith HO. Enzymatic assembly of DNA molecules up to several hundred kilobases. Nat Methods. Nature Publishing Group; 2009; 6: 343–345. https://doi.org/10.1038/nmeth.1318 PMID: 19363495

90. Gama-Castro S, Salgado H, Santos-Zavaleta A, Ledezma-Tejeida D, Muñiz-Rascado L, Garcia-Sotelo JS, et al. RegulonDB version 9.0: high-level integration of gene regulation, coexpression, motif clustering and beyond. Nucleic Acids Res. Oxford University Press; 2016; 44: D133–D143. https://doi.org/10.1093/nar/gkv1156 PMID: 26527724

91. Engler C, Kandzia R, Marillonnet S. A One Pot, One Step, Precision Cloning Method with High Throughput Capability. PLoS One. 2008; 3: e3647. https://doi.org/10.1371/journal.pone.0003647 PMID: 18985154
92. Chaveroche M-K, Ghigo J-M, d' Enfert C. A rapid method for efficient gene replacement in the filamen-
tous fungus Aspergillus nidulans. Nucleic Acids Res. 2000; 28: e97. PMID: 11071951

93. Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. Nature Publishing Group; 2008; 3: 1101–1108. https://doi.org/10.1038/nprot.2008.73

94. Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, et al. The STRING database in 2017: quality-controlled protein–protein association networks, made broadly accessible. Nucleic Acids Res. 2017; 45: D362–D368. https://doi.org/10.1093/nar/gkw937 PMID: 27924014