Defined Single-Gene and Multi-Gene Deletion Mutant Collections in Salmonella enterica sv Typhimurium

Abstract

We constructed two collections of targeted single gene deletion (SGD) mutants and two collections of targeted multi-gene deletion (MGD) mutants in Salmonella enterica sv Typhimurium 14028s. The SGD mutant collections contain (1), 3517 mutants in which a single gene is replaced by a cassette containing a kanamycin resistance (Kan^R) gene oriented in the sense direction (SGD-K), and (2), 3376 mutants with a chloramphenicol resistance gene (Cam^R) oriented in the antisense direction (SGD-C). A combined total of 3773 individual genes were deleted across these SGD collections. The MGD collections contain mutants bearing deletions of contiguous regions of three or more genes and include (3), 198 mutants spanning 2543 genes replaced by a Kan^R cassette (MGD-K), and (4), 251 mutants spanning 2799 genes replaced by a Cam^R cassette (MGD-C). Overall, 3476 genes were deleted in at least one MGD collection. The collections with different antibiotic markers permit construction of all viable combinations of mutants in the same background. Together, the libraries allow hierarchical screening of MGDs for different phenotypic followed by screening of SGDs within the target MGD regions. The mutants of these collections are stored at BEI Resources (www.beiresources.org) and publicly available.

About the Authors:

Steffen Porwollik

Affiliation: Department of Microbiology and Molecular Genetics, University of California Irvine, Irvine, California, United States of America

Carlos A. Santiviago

Affiliation: Departamento de Bioquímica y Biología Molecular, Facultad de Ciencias Químicas y Farmacéuticas, Universidad de Chile, Santiago, Chile

Pui Cheng

Affiliation: Department of Microbiology and Molecular Genetics, University of California Irvine, Irvine, California, United States of America

Fred Long

Affiliation: Department of Microbiology and Molecular Genetics, University of California Irvine, Irvine, California, United States of America

Prerak Desai

Affiliation: Department of Microbiology and Molecular Genetics, University of California Irvine, Irvine, California, United States of America

Jennifer Fredlund

Current address: Institut Pasteur, Paris, France

Affiliation: Department of Microbiology and Molecular Genetics, University of California Irvine, Irvine, California, United States of America

Shabarinath Srikumar

Current address: University of Liverpool, Liverpool, United Kingdom

Affiliation: Department of Microbiology and Molecular Genetics, University of California Irvine, Irvine, California, United States of America

Cecilia A. Silva

Affiliations Department of Microbiology and Molecular Genetics, University of California Irvine, Irvine, California, United States of America, Departamento de Bioquímica y Biología Molecular, Facultad de Ciencias Químicas y Farmacéuticas, Universidad de Chile, Santiago, Chile

Weiping Chu

Affiliation: Department of Microbiology and Molecular Genetics, University of California Irvine, Irvine, California, United States of America

Xin Chen

Affiliation: Department of Microbiology and Molecular Genetics, University of California Irvine, Irvine, California, United States of America

Rocío Canals

Current address: University of Liverpool, Liverpool, United Kingdom

Affiliation: Department of Microbiology and Molecular Genetics, University of California Irvine, Irvine, California, United States of America

M. Megan Reynolds

Affiliation: Department of Microbial Pathogenesis and Immunology, Texas A&M University, College Station, Texas, United States of America

Lydia Bogomolnaya

Affiliation: Department of Microbial Pathogenesis and Immunology, Texas A&M University, College Station, Texas, United States of America

Christine Shields

Affiliation: Department of Microbial Pathogenesis and Immunology, Texas A&M University, College Station, Texas, United States of America

Ping Cui

Affiliation: Department of Microbial Pathogenesis and Immunology, Texas A&M University, College Station, Texas, United States of America

Jinbai Guo

Affiliation: Department of Microbial Pathogenesis and Immunology, Texas A&M University, College Station, Texas, United States of America

Yi Zheng

Affiliation: Department of Microbial Pathogenesis and Immunology, Texas A&M University, College Station, Texas, United States of America

Tiana Endicott-Yazdani

Affiliation: Department of Microbial Pathogenesis and Immunology, Texas A&M University, College Station, Texas, United States of America

Hee-Jeong Yang

Affiliation: Department of Microbial Pathogenesis and Immunology, Texas A&M University, College Station, Texas, United States of America

Aimee Maple

Affiliation: Department of Microbial Pathogenesis and Immunology, Texas A&M University, College Station, Texas, United States of America

Yury Ragoza

Affiliation: Department of Microbial Pathogenesis and Immunology, Texas A&M University, College Station, Texas, United States of America

Carlos J. Blondel

Current address: Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts, United States of America

Affiliation: Departamento de Bioquímica y Biología Molecular, Facultad de Ciencias Químicas y Farmacéuticas, Universidad de Chile, Santiago, Chile

Camila Valenzuela

Affiliation: Departamento de Bioquímica y Biología Molecular, Facultad de Ciencias Químicas y Farmacéuticas, Universidad de Chile, Santiago, Chile

Helene Andrews-Polymenis

* E-mail: [email protected] (HAP); [email protected] (MM)

Affiliation: Department of Microbial Pathogenesis and Immunology, Texas A&M University, College Station, Texas, United States of America

Michael McClelland

* E-mail: [email protected] (HAP); [email protected] (MM)

Affiliation: Department of Microbiology and Molecular Genetics, University of California Irvine, Irvine, California, United States of America

Introduction

Salmonella enterica are found in many hosts including humans and their agricultural and companion animals. This bacterial species consists of a few known typhoidal serovars, such as S. Typhi or S. Paratyphi, that induce life-threatening enteric fever, and a large number of non-typhoidal serovars (NTS). In humans, NTS cause approximately 94 million cases of gastroenteritis and 155,000 deaths each year [1] and are the leading bacterial pathogens responsible for foodborne gastroenteritis [2], [3].

Salmonella enterica serovar Typhimurium ATCC14028 is a widely studied strain that was originally isolated by the Center of Disease Control in 1960 from pools of hearts and livers from 4-week-old chickens. Although many strains of serovar Typhimurium, including this strain, cause a systemic disease in BALB/c mice [4]–[6], this strain causes self-limiting gastroenteritis in virtually all hosts with an intact immune system. The strain persists in the intestine of immunologically intact mice including 129Sv and CBA/J [7] and in chicks [8]–[11]. The genome of this strain has been sequenced and annotated [12].

A widely used approach to identify a gene function is to generate a mutant with a deleted individual gene, followed by screening of the mutant for a phenotype. Generation of these mutants one at a time is time-consuming, so a resource to quickly obtain such mutants is of considerable interest to the research community. With this in mind, collections of defined deletion mutants have been established for several organisms in the past, including Escherichia coli K-12 [13], [14], Acinetobacter baylyi ADP1 [15] and several strains of the yeast Saccharomyces cerevisiae [16], [17]. The importance of establishing collections of such defined mutants is also recognized for metazoans such as Arabidopsis thaliana [18], Caenorhabditis elegans [19] and zebrafish [20]. These resources are invaluable for screening individual mutants for phenotypes [21], and for studies of epistasis and synthetic lethality [22]. As Salmonella Typhimurium is a widely used model organism, building a similar resource for Salmonella makes sense.

We report on the generation of the first comprehensive collections of both single-gene deletion for Salmonella enterica sv Typhimurium. We also present a first-of-its-kind set of systematic multi-gene deletion mutants. These collections are constructed with two markers to facilitate the subsequent combination of multiple mutations into the same genetic background. We further discuss how these resources are useful in infection models, where a high complexity of a population of different mutants is difficult to maintain throughout the entire process [23], and in situations where a host response to Salmonella must be screened, one mutant at a time [24].

Results

Generation of specific single-gene deletions in strain 14028

The GenBank annotation of Salmonella enterica sv Typhimurium 14028 contains 5372 protein coding sequences (CDSs) (GenBank GI:378448274). We excluded deletion targets for the following classes: rRNAs; tRNAs; ∼150 CDSs with syntenic orthologs known to be essential in three Escherichia coli K12 isolates [13], [14], [25], [26]; ∼100 structural elements of active lysogenic phage; most CDSs that occurred in more than one copy in the genome; and ∼750 CDSs under 100 amino acids in size that had no orthologous annotated counterpart in S. Typhimurium LT2 [27]. The genes not attempted for deletion are annotated with “nd” in Table S1.

We targeted the remaining genes including sRNAs for single-gene deletion (overall 4203 genes from strain 14028s). Primers were designed such that the resulting mutants would keep the first and last 30 bases of each targeted gene intact, to minimize interference with neighboring genes, while replacing the rest of the gene with an antibiotic resistance cassette. Primer sequences and positions can be found in Table S1 for all experiments that succeeded in mutant generation. For 26 CDSs, the mutant was generated in more than one experiment and listed twice.

The procedure for sequence replacement with an antibiotic resistance cassette followed the Lambda-Red technique [28]. The few differences from the original plasmid are described in detail elsewhere [23]. Briefly, we added a T7 RNA polymerase promoter, positioned such that it would generate a unique transcript from the Salmonella genome directly downstream of each mutant, thereby creating a mutant-specific RNA “signature”.

Each gene was deleted as both a Cam^R and a Kan^R construct in the same electroporation, and plated on separate selective plates (Figure 1A). The Kan^R and Cam^R resistance genes were oriented in the sense and anti-sense direction of the gene, respectively, to maximize any differences in polarity due to the active antibiotic resistance promoter. Deletion attempts in 4203 genes obtained at least one kanamycin-resistant (Kan^R) or chloramphenicol-resistant (Cam^R) mutant in 90% of cases. In most cases, two mutants of each type (two Cam^R and two Kan^R) were saved separately in individual wells in 96-well plates. Table 1 summarizes the overall results of all mutant construction attempts, and Table S1 presents the details of each mutant. Note that we also targeted a few genes annotated in S. Typhimurium LT2, but not annotated in 14028s (22 CDSs and 22 sRNA candidates), and included mutagenesis attempts for a few regions in the 14028s genome where alternative annotations had suggested a CDS to be present (31 mutants). We did not include the mutants that resulted from those efforts in the summary in Table 1, but they are displayed in detail in Table S1.

[Figure omitted. See PDF.]

Figure 1. Workflow of collection construction and verification.

A. General procedure and logistics of mutant generation, storage and characterization. B. Identification of obtained mutants and mapping of mutants to specific wells by a “Hypercube” Combinatorial Pooling approach, essentially as described in [30].

https://doi.org/10.1371/journal.pone.0099820.g001

[Figure omitted. See PDF.]

Table 1. Summary of collection statistics.

https://doi.org/10.1371/journal.pone.0099820.t001

Because both the Kan^R and the Cam^R PCR products to target each gene were transformed simultaneously, it was possible for both resistance cassettes to integrate into a single genome during the recombination event. Only nine such double mutants (5 in the SGD-K collection and 4 in the SGD-C collection) were identified by screening on LB agar plates supplemented with both Kan and Cam among over 8000 deletions. These double resistant mutants are annotated in Table S1.

Generation of specific multi-gene deletions in strain 14028

To establish which regions of the genome would survive deletion when growing in rich media, we analyzed the locations of 230,000 independent EZ-Tn5 <T7/KAN-2> (Epicentre Biotechnologies) insertion sites in the genome of strain 14028 [29]. Analysis of these locations resulted in identification of regions that would not tolerate disruptions (“essential regions”). We attempted deletion of 412 intervening “non-essential” regions of which 198 MGD-K and 251 MGD-C mutants were successfully constructed. Primer pairs that resulted in a successful multi-gene deletion in at least one variety (Kan or Cam) are listed in Table S1. Primers flanking the MGD regions to be deleted were offset from the outer ends of the 5′ and 3′ bordering genes by 60 nucleotides, allowing clear distinction from the primers used in the SGD collections where an offset by 30 nucleotides was used.

Identification of mutants and their locations in the library of 96-well plates

We used a combinatorial pooling approach [30] followed by Illumina sequencing to unambiguously verify the existence of mutants and simultaneously identify their plate and well locations in the collection (outlined in Figure 1B). Every Kan^R mutant (including both SGD-K and MGD-K mutants) was pipetted into a unique combination of seven out of twenty-four pools. The same procedure was performed for every Cam^R mutant (an additional twenty-four pools). Following a specific library preparation protocol (see Materials and Methods), sequence analysis then identified the presence of each mutant in each of the pools, thereby revealing the zip-code which pointed to a specific well in the collection.

During analysis of expected versus obtained zip-codes for each mutant, perfect concordance was found in the vast majority of cases. However, the robot occasionally failed to deposit a droplet into a pool or deposition of a droplet occurred during movement of the robotic arm. If over 75% of the digits of the obtained zip-codes concurred with the addresses in the expected zip code for that mutant, we assigned the mutant to that expected well. This was possible because we maximized the distance between the patterns used (error encoding), so that imperfect matches to the expected patterns still allowed mapping to a most probable well. Mutants that were present in the collection but remained un-mappable were also recorded.

This analysis identified a few systematic handling errors such as rotated plates that had occurred during collection establishment, leading to some mutants being identified in a well other than the one originally expected. These mis-assignments were corrected and the correct assignments are shown in Table S1.

For deletion of individual genes we attempted generation of mutants in 4203 genes of strain 14028s. Sequencing confirmed mutant generation for 3517 SGD-K and 3376 SGD-C clones, covering a total of 3773 genes in at least one variety. Analysis of the obtained zip-codes for each mutant successfully mapped 3331 SGD-K and 3175 SGD-C clones (covering a total of 3700 genes).

For multi-gene deletions, we targeted 412 regions spanning, in total, 4085 annotated features of the S. enterica sv Typhimurium 14028s genome. Sequencing of the junction of the mutants confirmed mutant generation for 198 MGD-K and 251 MGD-C clones, covering a total of 3476 genes. Zip-codes mapped 162 MGD-K and 234 MGD-C clones (covering a total of 3340 genes). Failed mutants may have been unfit, and smaller deletions in these regions can subsequently be attempted. Overall, the well locations of 93.8% of Kan^R mutants and 94.1% of Cam^R mutants were assigned with high probability to a single well (Table 1, Table S1).

Our sequencing approach was ideal for detection of systematic plate handling errors during assembly and storage of the collection and for detection of mutants in unexpected locations. The method was not designed to detect low levels of cross-contamination events where multiple mutants reside in the same well. A few such instances were detected through manual curation of the sequencing data. We tested 304 Kan^R mutants and 231 Cam^R mutants by PCR using primers to the flanking genomic regions of the expected mutant [23]. A cross-contaminant as well as a duplicated gene (see below) would have generated a second visible PCR product, particularly if the targeted mutant gene was less than 2 kb in size. We did not detect cross-contamination in any of these mutants. The few suspected and known contaminants are annotated in Table S1. Users should be aware that it is possible for additional wells to contain both an identified expected mutant and an unidentified contaminating mutant, most likely from an adjacent well.

Identification of insertions with functional T7 promoters

We investigated whether the T7 promoter introduced during mutagenesis was functional, using a previously described procedure [23] (see also Materials and Methods). DNA that included the T7 promoter in the mutant and flanking Salmonella genome sequence was PCR amplified. The resulting DNA was used to create fluorescently labeled cRNA with T7 RNA polymerase. The cRNA was hybridized to an array consisting of about 385,000 unique oligonucleotides, representing the entire genome of strain ATCC14028 on both strands in overlapping tiles at 12 bp intervals [31]. In three biological replicates and a total of five hybridizations, we repeatedly identified labeled transcript in the correct orientation and the correct location for the 3′ end of at least 3436 Kan^R mutants from our SGD-K collection where an annotated 14028s gene was deleted, and from at least 3280 such expected insertions from our SGD-C collection. We interrogated the presence of functional T7 promoters in our MGD-K collection by five microarray hybridization experiments. A transcript was detected in all 198 Kan^R mutants in the majority of these hybridizations. The status of detection of the expected transcripts for the mutants in all collections is listed in Table S1. The MGD-C collection was not checked for T7 promoter activity.

Evasion of mutagenesis by gene duplication

During establishment of the E. coli KEIO collection of targeted deletions in individual genes [13], [14], 25 of 7728 successfully established mutants evaded a gene knockout through duplication of the targeted gene. Most of these genes are now known to be essential or near-essential and we did not attempt to delete orthologs of this class of genes in Salmonella.

As mentioned above, we checked the first 304 mutants of our SGD-K collection and the first 231 mutants of our SGD-C collection by PCR using primers flanking the gene and found no evidence of duplication [23]. Next, we looked for mutants that appeared to be successfully deleted in our collection but might be near-essential, i.e. those that displayed a low frequency of transposon insertions in our previous work [29]. We investigated 30 Kan^R mutants and 33 Cam^R mutants for possible gene duplications including mutants in five genes with very few transposon insertions, indicating that these genes might be difficult to delete. We tested these candidates by qPCR amplification using primers inside the target genes (see Materials and Methods). In five mutants (two Kan^R and three Cam^R) a duplication of the gene was present along with the mutant. Two of these genes were among those selected because they had few transposon insertions and were, thus, perhaps hard to delete. The results of this assay are included in the comments column of Table S1. It is likely that a few other mutants in the collection contain duplicates. In some cases it is possible that these duplications fully restore fitness.

Collection usage alerts

Cross-contamination of some wells, as reported for the E. coli mutant collection [14], was also observed in our collection, and is likely to continue to occur when using the collection. It is also possible that a few successful deletion mutants may have second site mutations that alter their phenotypes, especially if these second site mutations improve growth. As is always good practice, we recommend mutant verification and transfer into a clean background for mechanistic follow-up studies.

Discussion

We present collections of defined single-gene and multi-gene deletion mutants in Salmonella enterica sv Typhimurium ATCC14028. Two collections are single-gene deletion (SGD) mutants in genes not required for growth in rich liquid media (LB), marked with either kanamycin resistance (Kan^R) or chloramphenicol resistance (Cam^R). The two additional collections of mutants are multi-gene deletion (MGD) mutants in chromosomal regions encompassing three or more genes that are not needed for growth of Salmonella enterica sv Typhimurium in LB, also marked with Kan^R or Cam^R. The generation of libraries with two different markers allows all viable pairwise combinations of deletions to be readily assembled in the same strain by generalized transduction.

The selectable marker within each insertion can be removed by FLP recombination leaving an in-frame scar, to diminish or eliminate polar effects on adjacent genes caused by the selective marker and its promoter [23], [28]. FLP recombination also allows three or more deletions to be assembled in the same isolate by sequential removal of antibiotic resistances followed by addition of a new marked mutation.

There are other advantages of these collections. Researchers that are engaged in functional genomics sometimes want to screen for as many phenotype-associated genes as possible in a single experiment. For this purpose, various techniques have been devised that employ pools of strains carrying random transposon mutant integrations, and monitor the change in representation of each mutant in the population [29], [30], [32]–[42]. To guarantee transposon insertions in almost all genes on the genome, including short genes, the number of random mutants that needs to be generated and screened is very high. Screening pools of large numbers of mutants is trivial in most in vitro conditions. However, in complex environments, such as infections in animals, the number of wild type bacterial cells surviving various stresses may be low at one or more bottlenecks, even if a huge initial dose is used.

To overcome random loss of mutants in complex environments, alternative strategies that use decreased numbers of mutants must be employed. First, one can simply reduce the number of random transposon insertion mutants in the pool used for infection [42]–[45]. This strategy lowers the number of genes screened or requires screening of separate sub-pools in different animals. Alternatively, one can use a subset of individual mapped transposon insertions, an approach that has been used in Staphylococcus aureus [46], Vibrio cholerae [47], Pseudomonas aeruginosa PAO1 [48], [49], and Burkholderia thailandensis [50]. Finally, mutants can be built with a defined deletion in each gene. Because only a single mutant is generated for each gene, the total number of different mutants required in the pool for maximum coverage of the genome is minimized.

Using defined mutants has a number of advantages: (a) Defined deletions usually involve deletion of the whole gene, instead of an insertion at an arbitrary location. (b) The mutants can be made in a way that the selectable marker and promoter can subsequently be removed in frame to downstream minimize any polar effects. (c) Different selectable markers can be introduced at the same location. Multiple markers permits construction of mutations in multiple genes in the same strain to measure, for example, redundancy or epistatic effects [14]. (d) Groups of adjacent, continuous, non-essential genes can be deleted, en bloc, so that a library of a few hundred mutants can encompass thousands of genes.

Large deletions are of utility in certain situations, such as screening organs after oral infection, where severe bottlenecks occur that cause random loss of mutants from complex pools. When a bottleneck allows only 20,000 bacteria to survive a step during infection, our MGD libraries of a few hundred bacteria can be used successfully. After such an MGD study identifies candidate genetic regions with functions in the investigated process, a smaller pool of only those SGD mutants representing the candidate regions can then be used to pinpoint the exact gene(s) with a function. We have used parts of the MGD and SGD libraries to validate this approach [51], [52].

The MGD collections are also useful in conditions where a host response elicited by specific genes of the bacterium is of interest. A pool of mutants cannot be used in such a scenario, because most mutants in the mixture will retain a functioning trigger and will elicit the response. Instead, such studies must be performed one mutant at a time. Part of our MGD library has been used in a hierarchical screen to identify regions in the Salmonella genome involved in down-modulation of expression of the human T-Cell Receptor, inhibition of blastogenesis, and proliferation of T cells, with subsequent mapping to single genes using SGDs [24].

Collection availability

All four collections have been deposited with BEI Resources (www.beiresources.org), with the following catalog numbers: SGD-K: NRC29399–29409, NRS42824–42858, NRS42859 (partial); SGD-C: NRC29410–29420, NRS42859 (partial), NRS42865–42899; MGD-K: NRS42861–42864; MGD-C: NRS42900–42902. Note that BEI will provide a maximum of 10 mutants at a time. Researchers interested in studies of mutant pools of our collections should contact MM directly.

Materials and Methods

Strain

ATCC14028s, a commonly used laboratory strain, was obtained from Dr. Don Guiney (UCSD), but was originally obtained from ATCC. We re-sequenced this strain to over 10× coverage [23] and found no differences with the published sequence [12].

Mutant generation

Primers for deletions were based on the published GenBank annotation of S. Typhimurium strain ATCC14028s [12] (GenBank GI:378448274) supplemented with in-house annotations using RAST and JCVI tools and the annotation of S. Typhimurium strain LT2, a very close relative [27]. Primers used in this project, along with their targeted coordinates in the 14028 chromosome and virulence plasmid (pSLT), mutant generation success and locations of the mutants in the collection plates are shown in Table S1. Wells without mutants either failed to produce a mutant or were blank.

The generation of targeted deletion mutants using the Lambda-Red mutagenesis is described in detail elsewhere [23], [28]. In brief, we amplified resistance cassettes from slightly modified pKD3 (Cam^R) and pKD4 (Kan^R) plasmids (named pCLF3 and pCLF4, respectively), using gene-specific primers. These primers are listed in Table S1. Unlike pKD3 and pKD4, plasmids pCLF3 and pCLF4 contain a T7 promoter to enable mutant-specific RNA signature production. All other features of pKD3 and pKD4, such as the FLP recombination sites, the ribosome binding site (RBS) and the start codon at the 3′ end of the mutagenic cassette for translational coupling, are retained (GenBank accession numbers EU629213 and EU629214 for pCLF3 and pCLF4, respectively). PCR products were used to transform competent S. Typhimurium ATCC14028 expressing the lambda Red recombinase by electroporation. Successful mutants were harvested from selective agar plates (supplemented with either 20 µg/ml Cam or 60 µg/ml Kan) after growth for 16–48 h. Two colonies were stocked for each mutant, only one of which was subsequently verified and deposited at BEI.

Array-based analysis of presence and functionality of mutants in pools

Separate pools for the different mutant classes (SGD-K, SGD-C, MGD-K), containing all generated mutants in each class, were grown overnight in selective LB broth. Genomic DNA was isolated using the GenElute Bacterial Genomic DNA kit (Sigma). This was done in three replicates for the SGD-K and MGD-K pools and in two replicates for the SGD-C pool. The region directly 3′ from each mutation in the pools was labeled and hybridized as described in detail previously [23]. Essentially, after shearing and polyA-tailing of the genomic DNA using 40 U terminal transferase in a 50 µl reaction, nested PCR amplified the subset of polyA-tailed DNA fragments that contained the T7 promoter encoded in the cassettes integrated into each mutant, along with the flanking genomic region. In the first PCR reaction, 50 ng of purified polyA-tailed DNA was used in a PCR with 0.05 U Taq (Promega, WI) and 0.2 µM of primers FRT Out 3-1 (TTCCTATACTTTCTAGAGAA), and PolyA-P (CCT₂₄VN) in 25 µl. The PCR reaction was performed under the following conditions: 94°C for 1 min; 30 cycles of 94°C for 10 s, 50°C for 10 s, and 72°C for 5 s; 72°C for 3 min. A subsequent nested PCR used 1 µl amplified product from the initial PCR in a total volume of 50 µl. Internal primer FRT Out 3-2 (TAGGAACTTCGGAATAGGAA) and primer CCT₂₄VN were used under identical cycling conditions. In vitro transcription using the AmpliScribe T7 transcription kit (Epicentre) and labeled Cy5- or Cy3-UTP (GE Healthcare) generated cRNA that was hybridized to a 384 k Nimblegen oligo tiling array covering the entire genome of S. Typhimurium 14028 in 12 bp overlaps. In each hybridization, 4 µg of labeled RNA was mixed with alignment oligo, hybridization buffer and Component A according to the Roche/Nimblegen hybridization protocol, and hybridized at 42°C for 16 hours. Arrays were washed according to the manufacturer's protocol, and scanned using a GenePix 4000B laser scanner (Molecular Devices, Sunnyvale, California) at 5 µm resolution. Each replicate gDNA for each pool was hybridized at least once. Array signals were quantified and analyzed using NimbleScan v2.4, WebarrayDB [53] and custom Perl scripts. For the SGD pools, a mutant was deemed to be present with an intact T7 promoter if the median signal intensity of the five closest strand-specific probes downstream from the expected insertion site of that respective mutant was higher than the 90^th percentile of the intensity values of all nonsense probe intensities on the array (i.e. the probes on the array that were hybridizing to the opposite strand of the expected T7 transcripts in the collection, which should reflect background signal). For the MGD pool, a plot of signal intensities for probes along the genome was visually inspected and all expected T7 transcript sites scored manually.

Identification of the location of all mutants in the collections using combinatorial pooling

Following a combinatorial pooling approach [30], every Kan^R mutant (including both SGD-K and MGD-K mutants) was initially pipetted robotically into a unique combination of seven out of twenty-four pools. The same procedure was performed for every Cam^R mutant. During this pipetting step, 24 pools of Kan^R mutants (P01–P24) and 24 pools of Cam^R mutants (P25–P48) were generated that each contained between 1325 and 1496 individual mutants, with every mutant being deposited in 7 different pools. A custom Perl script identified sufficient unique 7-out-of-24 pool combinations that each varied from all others by at least two pools. The resulting well-specific seven-pool “zip code” therefore differed from any other code by at least two numbers (error encoding). As an example, the Kan^R mutant of gene STM14_0001, which resides in plate SGD_067/068_Kan, well G11, was deposited in pools P01, P04, P06, P10, P17, P23 and P24, and is therefore assigned the corresponding “zip code” of 01-04-06-10-17-23-24. No other mutant was deposited in a pattern that contained more than 5 of these seven pools.

Genomic DNA of the resulting 24 Kan^R mutant pools and the 24 Cam^R mutant pools was prepared using the GenElute Bacterial Genomic DNA kit (Sigma). These gDNAs were subsequently denatured for 5 min at 95°C, randomly primed using 0.2 uM Random_Multiplexing_Primer (5′- GTGACTGGAGTTCAGACGTGTGCTCTTCCGATCTNNNNNNNNN) and extended, using 50 U exo- Klenow enzyme (New England Biolabs) and a gradual increase of temperature over 6 min to 37°C, where the reaction proceeded for 30 min. Following enzyme deactivation at 75°C for 20 min and subsequent QIAquick purification (Qiagen), the 3′ ends of products were tractably amplified by PCR using 2.5 U Taq enzyme (Invitrogen) and primers Read2_Illumina (5′ GTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT) and SGD_Index (5′ACACTCTTTCCCTACACGACGCTCTTCCGATCT(X₇)GACCTAAGGAGGATATTCA, where X₇ indicates a pool-specific known “barcode” sequence). The 5′ ends were amplified using primers Read2_Illumina and SGD_5′end_code_2 (5′ ACACTCTTTCCCTACACGACGCTCTTCCGATCTAGGAGTCTCGAAGCAGCTCCAGCCTG). In both reactions, primer Read2_Illumina was added after an initial extension cycle, following this regimen: 1 cycle of 94°C 3 min - 56°C 3 min - 72°C 5 min; addition of Read2_Illumina; 25 cycles of 94°C 15 sec - 60°C 30 sec - 72°C 1 min; final extension at 72°C 5 min. As a result of this approach, the 3′ends of the DNAs in every pool were now tagged with a specific “barcode”. Products were QIAquick purified (Qiagen) and assembled at roughly equimolar amounts.

The mixtures were then paired-end sequenced on the Genome Analyzer IIx (Illumina), and sequences were analyzed following a pipeline that involved Bowtie [54], de-duplication and custom Perl scripts to improve alignments to the genome of S. Typhimurium 14028. Reads that aligned within a 10-base window to expected genome locations were deemed to have been generated from that expected mutant. The “barcode” reads obtained for each mutant were then used to identify the exact physical plate/well location of each mutant in the collection. Data interpretation is reported in Table S1.

Test for target gene duplication

Primers were designed to amplify 140–170 bp internal fragments of the targeted genes. A qPCR assay was then performed, in duplicate, on boiled material from each mutant well, using Evagreen dye (Biotium) and the Kapa 2G Robust Hotstart ReadyMix PCR kit (Kapa Biosystems). For each mutant, a control reaction was included using an internal primer pair for a gene that had been targeted in a different mutant. The obtained C_Ts of control PCR versus internal PCR were compared, and a mean ΔΔC_Ts above 3 indicated lack of duplication, successful mutant generation, and a lack of significant contamination with another mutant. Mutants where the difference between the C_Ts was lower than 3 were purified to single colony and tested further by PCR of the mutant junction in the genome, performed with primer SGD_5_prime_out (5′-TCGAAGCAGCTCCAGCCTG), facing outward towards the 5′ end of the resistance cassette, and a primer adjacent to each targeted gene, facing towards the 5′ end of the resistance cassette. Successful PCR amplification indicated accurate integration of the resistance cassette into the target gene. Internal PCR amplification was repeated in this sample to reveal if there was target gene duplication.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

List of mutants present in our collections.

https://doi.org/10.1371/journal.pone.0099820.s001

(XLSX)

Acknowledgments

The authors thank Howard Brickner for technical support.

Author Contributions

Conceived and designed the experiments: MM HAP. Performed the experiments: SP PC FL PD JF SS CAS WC XC RC MMR LB CS PC JG YZ TEY HJY AM YR CJB CV CAS HAP. Analyzed the data: MM SP PD HAP. Wrote the paper: MM SP HAP.

Citation: Porwollik S, Santiviago CA, Cheng P, Long F, Desai P, Fredlund J, et al. (2014) Defined Single-Gene and Multi-Gene Deletion Mutant Collections in Salmonella enterica sv Typhimurium. PLoS ONE 9(7): e99820. https://doi.org/10.1371/journal.pone.0099820

References

1. Majowicz SE, Musto J, Scallan E, Angulo FJ, Kirk M, et al. (2010) The global burden of nontyphoidal Salmonella gastroenteritis. Clin Infect Dis 50: 882–889.

2. EFSA (2013) The European Union Summary Report on Trends and Sources of Zoonoses, Zoonotic Agents and Food-borne Outbreaks in 2011. EFSA Journal 11: 3129.

3. Gould LH, Walsh KA, Vieira AR, Herman K, Williams IT, et al. (2013) Surveillance for foodborne disease outbreaks - United States, 1998–2008. MMWR Surveill Summ 62: 1–34.

4. Baumler AJ, Tsolis RM, Bowe FA, Kusters JG, Hoffmann S, et al. (1996) The pef fimbrial operon of Salmonella typhimurium mediates adhesion to murine small intestine and is necessary for fluid accumulation in the infant mouse. Infect Immun 64: 61–68.

5. Bowe F, Lipps CJ, Tsolis RM, Groisman E, Heffron F, et al. (1998) At least four percent of the Salmonella typhimurium genome is required for fatal infection of mice. Infect Immun 66: 3372–3377.

6. Swearingen MC, Porwollik S, Desai PT, McClelland M, Ahmer BM (2012) Virulence of 32 Salmonella strains in mice. PLoS One 7: e36043.

7. Katribe E, Bogomolnaya LM, Wingert H, Andrews-Polymenis H (2009) Subspecies IIIa and IIIb Salmonellae are defective for colonization of murine models of salmonellosis compared to Salmonella enterica subsp. I serovar typhimurium. J Bacteriol 191: 2843–2850.

8. Pezoa D, Yang HJ, Blondel CJ, Santiviago CA, Andrews-Polymenis HL, et al. (2013) The type VI secretion system encoded in SPI-6 plays a role in gastrointestinal colonization and systemic spread of Salmonella enterica serovar Typhimurium in the chicken. PLoS One 8: e63917.

9. Andrews-Polymenis HL, Rabsch W, Porwollik S, McClelland M, Rosetti C, et al. (2004) Host restriction of Salmonella enterica serotype Typhimurium pigeon isolates does not correlate with loss of discrete genes. J Bacteriol 186: 2619–2628.

10. Izat AL, Adams MH, Cabel MC, Colberg M, Reiber MA, et al. (1990) Effects of formic acid or calcium formate in feed on performance and microbiological characteristics of broilers. Poult Sci 69: 1876–1882.

11. Waldroup AL, Skinner JT, Hierholzer RE, Waldroup PW (1993) An evaluation of fructooligosaccharide in diets for broiler chickens and effects on salmonellae contamination of carcasses. Poult Sci 72: 643–650.

12. Jarvik T, Smillie C, Groisman EA, Ochman H (2010) Short-term signatures of evolutionary change in the Salmonella enterica serovar typhimurium 14028 genome. J Bacteriol 192: 560–567.

13. Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, et al. (2006) Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol 2: 2006 0008.

14. Typas A, Nichols RJ, Siegele DA, Shales M, Collins SR, et al. (2008) High-throughput, quantitative analyses of genetic interactions in E. coli. Nat Methods 5: 781–787.

15. de Berardinis V, Vallenet D, Castelli V, Besnard M, Pinet A, et al. (2008) A complete collection of single-gene deletion mutants of Acinetobacter baylyi ADP1. Mol Syst Biol 4: 174.

16. Winzeler EA, Shoemaker DD, Astromoff A, Liang H, Anderson K, et al. (1999) Functional characterization of the S. cerevisiae genome by gene deletion and parallel analysis. Science 285: 901–906.

17. Mulleder M, Capuano F, Pir P, Christen S, Sauer U, et al. (2012) A prototrophic deletion mutant collection for yeast metabolomics and systems biology. Nat Biotechnol 30: 1176–1178.

18. Bolle C, Huep G, Kleinbolting N, Haberer G, Mayer K, et al. (2013) GABI-DUPLO: a collection of double mutants to overcome genetic redundancy in Arabidopsis thaliana. Plant J 75: 157–171.

19. Huang T-T, Gonzales CB, Gu F, Hsu Y-T, Jadhav RR, et al. (2013) Epigenetic deregulation of the anaplastic lymphoma kinase gene modulates mesenchymal characteristics of oral squamous cell carcinomas. Carcinogenesis

20. Varshney GK, Huang H, Zhang S, Lu J, Gildea DE, et al. (2013) The Zebrafish Insertion Collection (ZInC): a web based, searchable collection of zebrafish mutations generated by DNA insertion. Nucleic Acids Res 41: D861–864.

21. Nichols RJ, Sen S, Choo YJ, Beltrao P, Zietek M, et al. (2011) Phenotypic landscape of a bacterial cell. Cell 144: 143–156.

22. Butland G, Babu M, Diaz-Mejia JJ, Bohdana F, Phanse S, et al. (2008) eSGA: E. coli synthetic genetic array analysis. Nat Methods 5: 789–795.

23. Santiviago CA, Reynolds MM, Porwollik S, Choi SH, Long F, et al. (2009) Analysis of pools of targeted Salmonella deletion mutants identifies novel genes affecting fitness during competitive infection in mice. PLoS Pathog 5: e1000477.

24. Kullas AL, McClelland M, Yang HJ, Tam JW, Torres A, et al. (2012) L-asparaginase II produced by Salmonella typhimurium inhibits T cell responses and mediates virulence. Cell Host Microbe 12: 791–798.

25. Kato J, Hashimoto M (2007) Construction of consecutive deletions of the Escherichia coli chromosome. Mol Syst Biol 3: 132.

26. Hashimoto M, Ichimura T, Mizoguchi H, Tanaka K, Fujimitsu K, et al. (2005) Cell size and nucleoid organization of engineered Escherichia coli cells with a reduced genome. Mol Microbiol 55: 137–149.

27. McClelland M, Sanderson KE, Spieth J, Clifton SW, Latreille P, et al. (2001) Complete genome sequence of Salmonella enterica serovar Typhimurium LT2. Nature 413: 852–856.

28. Datsenko KA, Wanner BL (2000) One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc Natl Acad Sci U S A 97: 6640–6645.

29. Canals R, Xia XQ, Fronick C, Clifton SW, Ahmer BM, et al. (2012) High-throughput comparison of gene fitness among related bacteria. BMC Genomics 13: 212.

30. Goodman AL, Wu M, Gordon JI (2011) Identifying microbial fitness determinants by insertion sequencing using genome-wide transposon mutant libraries. Nat Protoc 6: 1969–1980.

31. Navarre WW, Porwollik S, Wang Y, McClelland M, Rosen H, et al. (2006) Selective silencing of foreign DNA with low GC content by the H-NS protein in Salmonella. Science 313: 236–238.

32. Sassetti CM, Boyd DH, Rubin EJ (2001) Comprehensive identification of conditionally essential genes in mycobacteria. Proc Natl Acad Sci U S A 98: 12712–12717.

33. Langridge GC, Phan MD, Turner DJ, Perkins TT, Parts L, et al. (2009) Simultaneous assay of every Salmonella Typhi gene using one million transposon mutants. Genome Res 19: 2308–2316.

34. Skurnik D, Roux D, Aschard H, Cattoir V, Yoder-Himes D, et al. (2013) A comprehensive analysis of in vitro and in vivo genetic fitness of Pseudomonas aeruginosa using high-throughput sequencing of transposon libraries. PLoS Pathog 9: e1003582.

35. van Opijnen T, Camilli A (2013) Transposon insertion sequencing: a new tool for systems-level analysis of microorganisms. Nat Rev Microbiol 11: 435–442.

36. Mazurkiewicz P, Tang CM, Boone C, Holden DW (2006) Signature-tagged mutagenesis: barcoding mutants for genome-wide screens. Nat Rev Genet 7: 929–939.

37. Barquist L, Langridge GC, Turner DJ, Phan MD, Turner AK, et al. (2013) A comparison of dense transposon insertion libraries in the Salmonella serovars Typhi and Typhimurium. Nucleic Acids Res 41: 4549–4564.

38. Pickard D, Kingsley RA, Hale C, Turner K, Sivaraman K, et al. (2013) A genomewide mutagenesis screen identifies multiple genes contributing to Vi capsular expression in Salmonella enterica serovar Typhi. J Bacteriol 195: 1320–1326.

39. Sabbagh SC, Lepage C, McClelland M, Daigle F (2012) Selection of Salmonella enterica serovar Typhi genes involved during interaction with human macrophages by screening of a transposon mutant library. PLoS One 7: e36643.

40. Silva CA, Blondel CJ, Quezada CP, Porwollik S, Andrews-Polymenis HL, et al. (2012) Infection of mice by Salmonella enterica serovar Enteritidis involves additional genes that are absent in the genome of serovar Typhimurium. Infect Immun 80: 839–849.

41. Khatiwara A, Jiang T, Sung SS, Dawoud T, Kim JN, et al. (2012) Genome scanning for conditionally essential genes in Salmonella enterica Serotype Typhimurium. Appl Environ Microbiol 78: 3098–3107.

42. Chaudhuri RR, Morgan E, Peters SE, Pleasance SJ, Hudson DL, et al. (2013) Comprehensive assignment of roles for Salmonella typhimurium genes in intestinal colonization of food-producing animals. PLoS Genet 9: e1003456.

43. Parsons DA, Heffron F (2005) sciS, an icmF homolog in Salmonella enterica serovar Typhimurium, limits intracellular replication and decreases virulence. Infect Immun 73: 4338–4345.

44. Lawley TD, Chan K, Thompson LJ, Kim CC, Govoni GR, et al. (2006) Genome-wide screen for Salmonella genes required for long-term systemic infection of the mouse. PLoS Pathog 2: e11.

45. Chaudhuri RR, Peters SE, Pleasance SJ, Northen H, Willers C, et al. (2009) Comprehensive identification of Salmonella enterica serovar typhimurium genes required for infection of BALB/c mice. PLoS Pathog 5: e1000529.

46. Fey PD, Endres JL, Yajjala VK, Widhelm TJ, Boissy RJ, et al. (2013) A genetic resource for rapid and comprehensive phenotype screening of nonessential Staphylococcus aureus genes. MBio 4: e00537–00512.

47. Cameron DE, Urbach JM, Mekalanos JJ (2008) A defined transposon mutant library and its use in identifying motility genes in Vibrio cholerae. Proc Natl Acad Sci U S A 105: 8736–8741.

48. Held K, Ramage E, Jacobs M, Gallagher L, Manoil C (2012) Sequence-verified two-allele transposon mutant library for Pseudomonas aeruginosa PAO1. J Bacteriol 194: 6387–6389.

49. Feinbaum RL, Urbach JM, Liberati NT, Djonovic S, Adonizio A, et al. (2012) Genome-wide identification of Pseudomonas aeruginosa virulence-related genes using a Caenorhabditis elegans infection model. PLoS Pathog 8: e1002813.

50. Gallagher LA, Ramage E, Patrapuvich R, Weiss E, Brittnacher M, et al. (2013) Sequence-defined transposon mutant library of Burkholderia thailandensis. MBio 4: e00604–00613.

51. Reynolds MM, Bogomolnaya L, Guo J, Aldrich L, Bokhari D, et al. (2011) Abrogation of the twin arginine transport system in Salmonella enterica serovar Typhimurium leads to colonization defects during infection. PLoS One 6: e15800.

52. Zheng Y, Sambou T, Bogomolnaya LM, Cirillo JD, McClelland M, et al. (2013) The EAL domain containing protein STM2215 (rtn) is needed during Salmonella infection and has cyclic di-GMP phosphodiesterase activity. Mol Microbiol 89: 403–419.

53. Wang Y, McClelland M, Xia XQ (2009) Analyzing microarray data using WebArray. Cold Spring Harb Protoc 2009: pdb prot5260.

54. Langmead B, Schatz MC, Lin J, Pop M, Salzberg SL (2009) Searching for SNPs with cloud computing. Genome Biol 10: R134.