IAI Accepted Manuscript Posted Online 19 January 2016 Infect. Immun. doi:10.1128/IAI.01423-15 Copyright © 2016 Grant et al. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.
1
Genes Required for the fitness of Salmonella enterica serovar Typhimurium
2
During Infection of Immunodeficient gp91-/-phox Mice
3 4
Running title: Salmonella TraDIS in immunodeficient mice
5 6
Andrew J. Grant1#, Olusegun Oshota1, Roy R. Chaudhuri1a, Matthew Mayho2, Sarah
7
E. Peters1, Simon Clare2, Duncan J. Maskell1, Pietro Mastroeni1
8 9 10
1 Department of Veterinary Medicine, University of Cambridge, Madingley Road, Cambridge, UK.
11 12
2 Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton,
13
Cambridge, UK.
14 15
a Current address: Department of Molecular Biology and Biotechnology, University
16
of Sheffield, Firth Court, Western Bank, Sheffield, UK.
17 18
#
To whom correspondence should be addressed:
[email protected]
19
1
20
ABSTRACT
21
Salmonella enterica causes systemic diseases (typhoid and paratyphoid
22
fever), non-typhoidal septicaemia (NTS) and gastroenteritis in humans and
23
other animals world-wide.
24
infectious disease problem in Sub-Saharan Africa is NTS in children and
25
immunocompromised adults. A current goal is to identify Salmonella mutants
26
that are not pathogenic in the absence of key components of the immune
27
system such as might be found in immunocompromised hosts.
28
attenuated strains have the potential to be used as live vaccines. We have
29
used Transposon Directed Insertion-site Sequencing (TraDIS) to screen
30
mutants of Salmonella enterica serovar Typhimurium for their ability to infect
31
and grow in the tissues of wild type and immunodeficient mice. This was to
32
identify bacterial genes that might be deleted for the development of live
33
attenuated vaccines that would be safer to use in situations and/or
34
geographical areas where immunodeficiencies are prevalent. The relative
35
fitness of each of 9,356 transposon mutants, representing mutations in 3,139
36
different genes, was determined in gp91-/-phox mice. Mutations in certain
37
genes led to reduced fitness in both wild type and mutant mice. To validate
38
tehse results these genes were mutated by allelic replacement and resultant
39
mutants were re-tested for fitness in the mice. A defined deletion mutant of
40
cysE was attenuated in C57BL/6 wild type mice and immunodeficient gp91-/-
41
phox mice, and was effective as a live vaccine in wild type mice.
An important but under-recognized emerging
Such
2
42
INTRODUCTION
43
Salmonella enterica causes systemic diseases (typhoid and paratyphoid
44
fever), non-typhoidal septicaemia and gastroenteritis in humans and other animals
45
world-wide. Salmonella enterica serovar Typhi (S. Typhi) and S. enterica serovar
46
Paratyphi (S. Paratyphi) are restricted to humans and cause typhoid and paratyphoid
47
fever, respectively. These are important febrile illnesses often seen in crowded and
48
impoverished populations with inadequate sanitation that are exposed to unsafe
49
water and food, and in travellers visiting countries where these diseases are
50
endemic (1). In 2010, typhoid and paratyphoid fevers were together estimated to
51
account for 190,200 deaths (2). Mortality rates in untreated typhoid fever can be 10-
52
15%. In some cases patients recover but remain carriers of the bacteria for many
53
years. Invasive non-typhoidal salmonellae (iNTS) are estimated to cause over 2.1
54
million illnesses and 416,000 deaths in Africa annually despite antimicrobial
55
treatment. Vaccine development for iNTS disease is a high priority (3). iNTS disease
56
is linked with conditions which impair the immune system such as lack of antibodies
57
in young children and multiple co-morbidities such as advanced HIV infection,
58
malaria, haemolysis, sickle cell disease (SCD), malnutrition and cytokine pathway
59
deficiencies (4, 5).
60
Current measures for controlling Salmonella infections are far from ideal. The
61
emergence of new multi-drug-resistant strains has reduced the usefulness of most
62
antimicrobials (6). Prevention by implementation of hygiene measures alone is
63
proving insufficient. Currently licensed vaccines are far from optimal and no effective
64
vaccine or delivery system effective against NTS is available.
3
65
Live attenuated S. enterica vaccines are more effective in many different
66
infection systems than non-living ones due to their ability to generate protective T-
67
helper 1 type immunity in addition to antibody responses (7). New live attenuated
68
vaccines against Salmonella are currently being developed and tested (8). However,
69
concerns remain about the use of live vaccines in general when these are targeted
70
to endemic areas with a high incidence of conditions that impair the immune system,
71
and there is reluctance to use live vaccines in animals bred for food production due
72
to fears of side-effects subsequent to the unlikely, but possible, spread to humans.
73
Reversion to virulence would pose a serious threat to individuals with an impaired
74
immune system who may inadvertently receive the live vaccine. Immuno-
75
suppression can be transient, latent or undiagnosed, possibilities that are likely to be
76
more frequent in children or in adults from countries where sophisticated diagnostic
77
facilities are not always available.
78
We and others have found that some live attenuated Salmonella vaccines,
79
such as aromatic-dependent (aro) mutants, htrA mutants and aroA htrA double
80
mutants can cause lethal infections in immunodeficient mice such as those lacking in
81
the production of reactive oxgen species (ROS), IL-12 or IFNγ and in mice lacking T-
82
cells (9, 10). Bacterial mutants of these genes could therefore be dangerous if used
83
in patients with similar immune defects. We have also shown that other vaccines,
84
specifically Salmonella Pathogenicity Island-2 (SPI-2) mutants of S. enterica, similar
85
to the ones currently being tested in humans, regain virulence in ROS-deficient
86
animals (11) and could therefore prove dangerous in individuals with reduced
87
efficiency of the oxygen dependent antimicrobial mechanisms of phagocytes (e.g.
88
chronic granulomatous disease (CGD) or malaria patients) (12, 13). These
4
89
immunodeficiencies may very well be latent in the individual at the time of
90
vaccination.
91
A current goal is to identify mutations that would allow the generation of live
92
attenuated salmonellae, that would not significantly regain virulence in the absence
93
of key components of the immune system. The feasibility of this approach is
94
supported by the evidence that some Salmonella mutants can retain this attenuation
95
in situations of serious immunodeficiency in animal models of systemic infection. For
96
example, SPI-2 mutants are attenuated in gene-targeted mice lacking IFNγ, despite
97
these strains being virulent in ROS-deficient animals (14).
98
Recent advances in high-throughput techniques to assess simultaneously the
99
genotypes and relative fitness of individuals in complex pools of bacterial transposon
100
(Tn) mutants allow us to utilise a genome-wide screen for genes necessary for
101
bacterial fitness in the face of particular immunological challenges.
102
Directed Insertion-site Sequencing (TraDIS) is one such technique. TraDIS exploits
103
very high throughput Illumina sequencing to obtain sequence reads directly from the
104
region of DNA flanking each Tn insert in pools of mutants (15). This allows the
105
location of the Tns to be precisely determined in a massively parallel manner. The
106
number of reads obtained for each Tn allows the relative abundance of each mutant
107
to be assessed, and a quantitative measure of fitness can be obtained by comparing
108
data obtained before and after a selection step.
Transposon
109
In the current study we have used a stringent mouse model of
110
immunodeficiency to perform a genome-wide screen with the potential to identify
111
genes that can disable reversion to virulence in immunodeficient hosts and that
5
112
could therefore be included in new generations of live vaccines against invasive
113
salmonelloses.
6
114 115
MATERIALS AND METHODS Bacterial strains, media and growth conditions.
All wild type strains,
116
defined mutants and plasmids are summarised in Table 1. Preparation of
117
electrocompetent Escherichia coli and S. Typhimurium cells and transformations
118
were performed as previously described (21). Bacteria were grown on Luria-Bertani
119
(LB) medium. Media were supplemented with the appropriate antibiotic for selection
120
(ampicillin 100μg/ml, kanamycin 50μg/ml, chloramphenicol 10μg/ml, tetracycline
121
12.5μg/ml).
122 123
Recombinant DNA techniques: Standard methods were used for molecular
124
cloning (22).
Chromosomal and plasmid DNA purification and routine DNA
125
modifications including restriction endonuclease digestion of DNA, modifications of
126
DNA and ligations were carried out using commercial kits and supplies according to
127
the manufacturers’ instructions (QIAGEN; Thermo Fisher Scientific; New England
128
Biolabs). DNA concentration and purity were measured using a Nanodrop ND-1000
129
spectrophotometer.
130
(http://frodo.wi.mit.edu/) and purchased from Sigma (Sigma-Genosys, UK). PCRs
131
were performed in 25 μl reaction volumes in 0.2 ml Eppendorf tubes in a Perkin
132
Elmer Gene Amp 2400 thermal cycler. Reactions contained 200 M dNTPs, 2 mM
133
Mg2+, 0.01 volumes of Proof Start DNA polymerase (QIAGEN; 2.5 U μl-1), 0.1
134
volumes polymerase buffer (10x), 1 μM forward and reverse primers and template
135
DNA (~50 ng plasmid DNA or ~100 ng chromosomal DNA). Thermal cycler
136
conditions were 94ºC for 10 min, then 35 cycles of 94ºC for 1 min, 55ºC for 1 min
137
and 72ºC for 1 min, followed by a final extension at 72ºC for 10 min.
PCR primers (Table S1) were designed using Primer3
7
138 139
Generation of defined mutants of S. Typhimurium. Mutants were
140
generated using a modification of the ET-cloning procedure (23, 24) as previously
141
described (20). PCR was used to amplify the chloramphenicol resistance cassette
142
from pACYC184 (18), the kanamycin resistance cassette from pACYC177 (18) or
143
the tetracycline resistance cassette from pBR322 (19) with 5’ and 3’ 60 bp homology
144
arms complementary to the flanking regions of the gene to be deleted.
145
Approximately 1 μg of linear PCR product was used for integration into the
146
chromosome using a modification of the Lambda Red method (25), as previously
147
detailed (26).
148
chloramphenicol, kanamycin or tetracycline as appropriate. Screening for loss of the
149
pBADλred helper plasmid was as previously described (27), using MAST ID
150
Intralactam circles (MAST Diagnostics, Bootle Merseyside, UK) to screen for the
151
absence of beta-lactamase in bacterial colonies. (Further details of construction of
152
mutants are provided in Supplemental Material).
Transformants were selected by plating onto media containing
153 154
Mouse infections. All animals were handled in strict accordance with good
155
animal practice as defined by the relevant international (Directive of the European
156
Parliament and of the Council on the protection of animals used for scientific
157
purposes, Brussels 543/5) and local (Department of Veterinary Medicine, University
158
of Cambridge) animal welfare guidelines.
159
ethical review committee of the University of Cambridge and was licensed by the UK
160
Government Home Office under the Animals (Scientific Procedures) Act 1986. Sex-
161
and aged-matched 9-12 week old C57BL/6 wild type mice (Harlan Olac Ltd), and
All animal work was approved by the
8
162
gp91-/-phox mice (bred at the Wellcome Trust Sanger Institute) were infected by
163
intravenous (i.v.) injection of bacterial suspensions in a volume of 0.2 ml. Static
164
grown cultures of wild type S. Typhimurium SL1344 and defined mutants were grown
165
overnight at 37ºC from single colonies in 10 ml LB broth, and then diluted in
166
phosphate buffered saline (PBS) to the appropriate concentration for inoculation.
167
Bacterial cultures for the pools of Tn5 and Mu Tn mutants (TraDIS pools) were
168
grown as described previously (14). Inocula were enumerated by plating dilutions
169
onto LB agar plates. Mice were killed by cervical dislocation and the livers and
170
spleens were aseptically removed and homogenized in sterile water using a
171
Colworth Stomacher 80. The resulting homogenate was diluted in a 10-fold series in
172
PBS and LB agar pour plates were used to enumerate viable bacteria.
173 174
TraDIS analysis of S. Typhimurium mutant pools. The Illumina sequence
175
reads from the input (bacteria inoculated into mice) and output (bacteria recovered
176
from infected livers and spleens) pools (European Nucleotide Archive ERA201074)
177
were generated using single-end short-read sequencing with Tn specific primers and
178
by massively parallel sequencing of Tn-flanking regions in the input and output
179
mutant pools. Genomic DNA was prepared and 2 µg fragmented by Covaris to an
180
average size of ~300 bp. Following end-repair and ‘A’ tailing a modified Illumina
181
adapter, synthesized and annealed by Integrated DNA Technologies (IDT) using
182
oligonucleotides SplA5_top and SplA5_bottom (Table S1), was ligated to the
183
fragments for 40 min at 20oC. Ligated fragments were cleaned using Ampure XP
184
beads (Beckman Coulter) with a bead to sample ratio of 0.8 to 1. PCR enrichment of
185
fragments containing the Tn was performed using a primer homologous to each end
186
of the Tn (Mu_AG_5’PCR, Mu_AG_3’PCR, Tn5_AG_5’PCR or Tn5_AG_3’PCR; 9
187
Table S1) in conjunction with an adapter-specific primer containing an index tag
188
(SplAP5.x; Table S1).
189
(Beckman Coulter), quantified by qPCR and then pooled. The resultant products
190
were sequenced on a HiSeq2500 (Illumina) using a specially modified recipe to
191
overcome difficulties generated by the monotemplate Tn sequence. Briefly, a Tn-
192
specific sequencing primer (Mu_AG_5’seq, Mu_AG_3’seq, Tn5_AG_5’seq or
193
Tn5_AG_3’seq; Table S1) anneals to the Tn 10 bases away from the genomic DNA
194
junction. Sequencing takes place with no imaging for the first 10 cycles followed by
195
imaging for the next 40 cycles. This gives a 40bp genomic DNA read. The template
196
is denatured and the same sequencing primer is re-annealed and 10 cycles of
197
sequencing takes place to give a 10 bp Tn read.
PCR products were cleaned using Ampure XP beads
198 199
Sequence Analysis. Illumina reads generated from the genomic DNA of the
200
input and output pools containing Tn5 and Mu Tn insertions were stripped of the 5’
201
Tn tags using Cutadapt version 1.8.1 (28), and the remainder of each sequence read
202
was mapped to the S. Typhimurium SL1344 chromosome and plasmid sequences
203
(Genbank accession numbers FQ312003, HE654724, HE654725 and HE654726)
204
using BWA mem version 0.7.12-r1039 (29). The raw input and output read counts
205
were loaded into R version 3.0.2. A step was performed to eliminate nonspecific
206
background reads by filtering out the insertions with a total read count across all
207
samples of 94% of the Tn mutants
259
screened were mapped). The reason for the lower number of mapped mutants
260
compared with the number of mutants in the original input pools could be due to
261
mutants not surviving the initial grow-up in the input pool, and/or due to the presence
262
of duplicate mutants in the library. A total of 3,139 different genes were disrupted by
263
the 9,356 Tn insertions.
264
The sites of insertion of all Tns in the S. Typhimurium TraDIS mutant libraries,
265
the identities of the genes that are disrupted, and the fitness scores of the mutants
266
are listed in Table S2. The fitness scores of the mutants obtained during infection of
267
C57BL/6 mice showed a high proportion of seemingly attenuated mutants (Fig. S1a).
268
Additionally, the output read counts from the two replicate C57BL/6 mice did not
269
show a good correspondence (Fig. S1b), and the fitness scores did not correlate well
270
with those obtained in BALB/c mice (Fig. S1c) from our previous study (32). From
271
these observations we concluded that the C57BL/6 mice had exhibited stochastic
272
loss of mutants for reasons unrelated to their genotype (i.e. random dropout). By
273
contrast, only a small proportion of mutants were attenuated in the gp91-/-phox mice,
274
with good correspondence between the replicates and with the data obtained from
275
BALB/c mice (Fig. S2). The random dropout was not identified in the BALB/c mice,
276
probably because of the increased bacterial load used as the inoculum in this model 13
277
(~1x106 CFU in the BALB/c mice vs ~5x104 CFU in the C57BL/6 mice), or in the
278
gp91-/-phox mice, due to the absence of bacterial killing in this model. Thus, the
279
attenuated mutants in the C57BL/6 dataset are likely to include false negatives due
280
to the random dropout of mutants in this model (i.e. mutants that have been
281
assigned as attenuated, but that could be colonisation proficient). This prevented an
282
in-depth comparison of the different genes that contributed to fitness in the wild type
283
C57BL/6 and immunocompromised mice. However, the data from the screen in the
284
gp91-/-phox mice were robust and provided us with solid foundations to select
285
representative genes from the ‘attenuated mutant’ lists and testing defined allelic
286
replacement mutants of these genes in wild type and immunocompromised mice.
287 288
Analysis of defined allelic replacement S. Typhimurium mutants in
289
C57BL/6 WT and gp91-/-phox mice. An ideal live attenuated vaccine would need to
290
colonise but be attenuated in an immune-competent host and to remain attenuated
291
in conditions where the immune system may be impaired. Candidate genes for
292
generation of such vaccines will be amongst those that lead to attenuated mutant
293
phenotypes in spleens and livers of both C57BL/6 WT and gp91-/-phox mice. We
294
selected 14 genes that fit this pattern: cydC, cydD, cysE, dksA, ftsK, miaA, nuoK,
295
ptsI, recD, secG, seqA, sucA, thdF and yqiC and deleted them by allelic replacement
296
in S. Typhimurium. We also generated a double mutant knocking out both cydC and
297
cydD, a mutant containing a deletion of the entire suc operon, a mutant containing a
298
deletion of the entire nuo operon, a mutant containing a deletion of the tol operon
299
(since many different tol mutants were attenuated, and tol genes have previously
300
been shown to have a role in virulence/fitness in mice (33)), and a mutant containing
301
a deletion of the tol/pal operon. All of these were tested individually for attenuation in 14
302
C57BL/6 WT mice compared to the SL1344 wild type, administered by i.v. injection,
303
and at input doses as detailed in Table 2. We also included in the experiment
304
controls mutants lacking aroC, aroC/aroD/htrA, and aroC/ssaV, which have
305
established roles in virulence and have been or are being used as vaccine
306
candidates (34-36). Each gene, or genes, was inactivated separately by λRed
307
recombinase-mediated integration of linear PCR amplicons by homologous
308
recombination. Mice were killed at 7 days p.i. or earlier if clinical signs became
309
apparent (Table 2). Defined mutants in cydC, cydD, cydC/cydD, recD, seqA and the
310
suc operon were not attenuated in the C57BL/6 mice, and mice had to be killed at 3
311
days p.i., along with the mice infected with the wild type SL1344. Mice infected with
312
defined mutants in ftsK, nuoK and the nuo operon mutants were killed on day 4 p.i.
313
as the mice displayed clinical signs. Mice infected with a dksA mutant were killed on
314
day 5 p.i. Mice infected with sucA and thdF mutants were displaying clinical signs at
315
6 days p.i. and were killed, although interestingly viable bacterial counts in the livers
316
and spleens were not high enough to be consistent with the observed signs of
317
infection, indicating that the clinical signs might be due to focal infection at other
318
body sites. Mice infected with cysE, miaA, ptsI, secG, tol operon, tol/pal operon and
319
yqiC mutants as well as the controls infected with aroC, aroC/aroD/htrA and
320
aroC/ssaV, did not display any visible clinical signs on day 7 p.i.
321
Given their attenuation in WT mice, mutants in cysE, miaA, ptsI, secG and
322
yqiC, along with the tol operon and tol/pal operon were used to infect gp91-/-phox
323
mice. We also included aroC, aroC/aroD, and ssaV mutants as controls. Mice were
324
killed at 7 days p.i., or earlier if clinical signs became apparent (Table 3). Mice
325
infected with defined ptsI, ssaV, tol operon and tol/pal operon mutants were killed on
326
day 3 p.i. when they displayed clinical signs. Mice infected with aroC, miaA and 15
327
secG mutants were killed on day 4 p.i. Mice infected with an yqiC mutant were killed
328
on day 5 p.i. when they displayed clinical signs. Mice infected with aroC/aroD or
329
cysE mutants were killed on day 7 p.i. which was the pre-defined end point of the
330
experiment. The cysE and the aroC/aroD mutants had higher bacterial counts in the
331
organs of the gp91-/-phox mice than in the C57BL/6 WT mice. In order to determine
332
how long the gp91-/-phox mice could survive infection with a cysE mutant, we
333
performed a small pilot study, infecting gp91-/-phox mice with Log10 5.72 CFU
334
SL1344 cysE and then allowing the infection to proceed until the mice displayed
335
clinical symptoms. The mice had to be killed on day 8 p.i. when the bacterial loads in
336
the livers were Log10 8.33 + 0.11 CFU (error given as standard deviation) and the
337
bacterial loads in the spleens were Log10 7.81 + 0.05 CFU (error given as standard
338
deviation). We next proceeded to compare cysE against an aroC mutant in gp91-/-
339
phox mice (Fig. 2). At three different input doses, the time till the appearance of
340
clinical signs was consistently delayed for the cysE mutant compared to the aroC
341
mutant.
342
As a prerequisite to a vaccination study, we determined the time it took for
343
SL1344 cysE to be cleared from the systemic organs (liver and spleen) of WT
344
C57BL/6 mice (Fig. 3). We found that by day 84 p.i. mice infected i.v. had cleared all
345
SL1344 cysE from the systemic organs, while in the orally infected group it took until
346
day 91 p.i. for the organs to be clear of the SL1344 cysE.
347
We next vaccinated WT C57BL/6 mice, either i.v. or orally, with SL1344 cysE
348
and then challenged the mice, 91 days post immunisation, orally with virulent
349
SL1344 and monitored the bacterial counts in the livers and spleens at various times
350
post infection (Fig. 4). The majority of the vaccinated mice were protected from
351
challenge with the virulent SL1344. However, during the challenge experiment it was 16
352
necessary to kill six mice (five from the orally immunised group, and one from the i.v.
353
immunised group) because they were displaying signs of infection. Thus, 97% of the
354
SL1344 cysE i.v. vaccinated C57BL/6 mice, and 86% of the SL1344 cysE orally
355
vaccinated C57BL/6 mice, were protected against challenge with virulent SL1344.
17
356
CONCLUSION
357
Here we have presented a report of fitness of random Tn mutants of S.
358
Typhimurium in vivo and described attenuated mutants based on the TraDIS
359
analysis of pools of Tn mutants screened in C57BL/6 WT and gp91-/-phox mice. Our
360
study provides proof-of-principle data on the feasibility of using TraDIS analysis in
361
combination with infections of gene-targeted animals to allow identification of
362
bacterial mutants as live vaccine candidates that would be safe in situations where
363
the immune system is impaired.
364
We show that through TraDIS analysis it is possible to identify novel
365
mutations in S. Typhimurium genes that still allow colonisation in WT mice but
366
importantly also reduce bacterial growth rate in vivo in severly immunodeficient gp91-
367
/-
phox mice. It is worth noting that the variance of defined mutants from their TraDIS
368
fitness scores (Tables 2 and 3) is likely to be due to differences in competition
369
dynamics for a given mutant relative to co-screened mutant bacteria, as we
370
previously demonstrated (32), and the differences in time-points at which the
371
bacteria were enumerated p.i. between the TraDIS screen and the defined mutant
372
experiments. Thus, some of the attenuated mutants identified in the screen are
373
virulent when tested individually, while some others retain the attenuation, and
374
therefore the screen is useful, but requires further evaluation of mutations
375
individually.
376
Through the TraDIS screens and validation of defined mutants, we identified
377
that SL1344 cysE was attenuated for growth, and cleared, in WT C57BL/6 mice and
378
had reduced growth in gp91-/-phox mice, although these mice succumbed to the
379
infection 8 days p.i. SL1344 cysE offered protection when given i.v. or orally as a
380
live vaccine to C57BL/6 WT mice, however, a few of the vaccinated mice succumbed 18
381
to the challenge with the virulent SL1344, indicating that cysE does not offer full
382
protection to all vaccinated animals.
383
Our work did not identify a single-gene S. Typhimurium mutant that would
384
colonise WT animals at sufficient levels to be predicted to be immunogenic and that
385
would also be completely unable to cause disease in the severly immunodeficient
386
animals used in this study. The attractive possibility that different combinations of
387
mutations could yield this ‘ideal’ vaccine candidate remains to be explored. However,
388
even mutations that reduce growth in immunocompromised hosts would be
389
advantageous and could be used to construct live strains in a way that would
390
crucially enable more time for medical intervention should an immunocompromised
391
vaccinee develop signs of infection. Furthermore, the animal model that we used
392
was of extreme immunodeficiency, while it is likely that in the field one would
393
encounter situations where macrophage function is only partially impaired.
394
Therefore it remains to be established whether some of the mutants identified from
395
this study would be safer in models of partial immunodeficiency.
396
Extending the approach reported here to other gene-targeted mouse strains
397
has the potential in the future to generate a more comprehensive picture of the
398
functional interactions between the bacterial and host genomes and to identify
399
bacterial mutants as live vaccine candidates that would be safe across a wide range
400
of immunodeficiencies.
19
401
FUNDING INFORMATION
402
This work was supported by a Medical Research Council (MRC) grant
403
G1100102. Tn libraries were previously constructed with funding from the
404
Biotechnology and Biological Sciences Research Council grant APG19115. O.O.
405
was supported by a Newton Trust grant awarded to A.J.G. M.M. was supported by
406
the Wellcome Trust grant number WT098051. The authors have no conflicting
407
financial interests. The funders had no role in study design, data collection and
408
interpretation, or the decision to submit the work for publication.
20
409 410 411
REFERENCES 1. Crump, J.A. and Mintz, E.D. 2010. Global trends in typhoid and paratyphoid Fever. Clin. Infect. Dis. 50: 241-246.
412 413
2. Loaano, R. Hanhavi, M. Foreman, K. Lim, S. Shibuya, K. Aboyans, V.
414
Abraham, J. Adair, T. Aggarwal, R. Ahn, SY. Alvarado, M. Anderson, H.R.
415
Anderson, L.M. Andrews, K.G. Atkinson, C. Baddour, L.M. Barker-Collo,
416
S. Bartels, D.H. Bell, M.L. Benjamin, E.J. Bennett, D. Bhalla, K. Bikbov, B.
417
Bin Abdulhak, A. Birbeck, G. Blyth, F. Bolliger, I. Boufous, S. Bucello, C.
418
Burch, M. Burney, P. Carapetis, J. Chen, H. Chou, D. Chugh, S.S.
419
Coffeng, L.E. Colam, S.D. Colquhon, S. Colson, K.E. Condon, J. Connor
420
M.D. Cooper, L.T. Corriere, M. Cortinovis, M. de Vaccaro, K.C. Couser,
421
W. Cowie, B.C. Criqui, M.H. Cross, M. Dabhadkar, K.C. Dahodwala, N. De
422
Leo, D. Degenhardt, L. Delossantos, A. Denenberg, J. Des Jarlai, D.C.
423
Dharmaratne, S.D. Dorsey, E.R. Driscoll, T. Duber, H. Ebel, B. Erwin, P.J.
424
Espindola, P. Ezzati, M. Feigin, V. Flaxman, A.D. Forouzanfar, M.H.
425
Fowkes, F.G. Franklin, R. Fransen, M. Freeman, M.K. Gabriel, S.E.
426
Gakidou, E. Gaspari, F. Gillum, R.F. Gonzalez-Medina, D. Halasa, Y.A.
427
Haring, D. Harrison, J.E. Havmoeller, R. Hay, R.J. Hoen, B. Hotez, P.J.
428
Hoy, D. Jacobsen, K.H. James, S.L. Jasrasaria, R. Jayaraman, S. Johns,
429
N. Karthikeyan, G. Kassebaum, N. Keren, A. Khoo, J.P. Knowlton, L.M.
430
Kobusingye, O. Koranteng, A. Krishnamurthi, R. Lipnick, M. Lipshultz,
431
S.E. Ohno, S.L. Mabweijano, J. MacIntyre, M.F. Mallinger, L. March, L.
432
Marks, G.B. Marks, R. Matsumori, A. Matzopoulos, R. Mayosi, B.M.
433
McAnulty, J.H. McDermott, M.M. McGrath, J. Mensah, G.A. Merriman, 21
434
T.R. Michaud, C. Miller, M. Miller, T.R. Mock, C. Mocumbi, A.O. Mokdad,
435
A.A. Moran, A. Mulholland, K. Nair, M.N. Naldi, L. Narayan, K.M. Nasseri,
436
K. Norman, P. O’Donnell, M. Omer, S.B. Ortblad, K. Osborne, R. Ozgediz,
437
D. Pahari, B. Pandian, J.D. Rivero, A.P. Padilla, R.P. Perez-Ruiz, F.
438
Percio, N. Phillips, D. Pierce, K. Pope, C.A. 3rd. Porrini, E. Pourmalek, F.
439
Raju, M. Ranganathan, D. Rehm, J.T. Rein, D.B. Remuzzi, G. Rivara, F.P.
440
Roberts, T. De Leon, F.R. Rosenfeld, L.C. Rushton, L. Sacco, R.L.
441
Salomon, J.A. Sampson, U. Sanman, E. Schwebel, D.C. Segui-Gomez, M.
442
Shepard, D.S. Singh, D. Singleton, J. Sliwa, K. Smith, E. Steer, A. Taylor,
443
J.A. Thomas, B. Tleyjeh, I.M. Towbin, J.A. Truelsen, T. Undurraga, E.A.
444
venketasubramanian, N. Vijayakumar, L. Vos, T. Wagner, G.R. Wang, M.
445
Wang, W. Watt, K. Weinstock, M.A. Weintraub, R. Wilkinson, J.D. Woolf,
446
A.D. Wulf, S. Yeh, P.H. Yip, P. Zabetian, A. Zheng, Z.J. Lopez, A.D.
447
Murray, C.J. AlMazroa, M.A. Memish, Z.A. 2012. Global and regional
448
mortality from 235 causes of death for 20 age groups in 1990 and 2010: a
449
systematic analysis for the Global Burden of Disease Study 2010. Lancet.
450
380: 2095-2128.
451 452
3. Ao, T.T. Feasey, N.A. Gordon, M.A. Keddy, K.H. Angulo, F.J. and Crump,
453
J.A. 2015. Global burden of invasive nontyphoidal Salmonella disease
454
2010(1). Emerg. Infect. Dis. 21: 6.
455 456
4. Gordon, M.A. Banda, H.T. Gondwe, M. Gordon, S.B. Boeree, M.J. Walsh,
457
A.L. Corkill, J.E. Hart, C.A. Gilks, C.F. and Molneux, M.E. 2002. Non-
22
458
typhoidal salmonella bacteraemia among HIV-infected Malawian adults: high
459
mortality and frequent recrudescence. AIDS. 16: 1633-1641.
460 461 462
5. Gordon, M.A. 2008. Salmonella infections in immunocompromised adults. J. Infect. 56: 413-422.
463 464
6. Gordon, M.A. Graham, S.M. Walsh, A.L. Wilson, L. Phiri, A. Molyneux, E.
465
Zijlstra, E.E. Heyderman, R.S. Hart, C.A. and Molyneux, M.E. 2008.
466
Epidemics of invasive Salmonella enterica serovar enteritidis and S. enterica
467
Serovar typhimurium infection associated with multidrug resistance among
468
adults and children in Malawi. Clin. Infect. Dis. 46: 963-969.
469 470
7. Harrison,
J.A.
Villarreal-Ramos,
B.
Mastroeni,
P.
Demarco
de
471
Hormaeche, R. Hormaeche, C.E. 1997. Correlates of protection induced by
472
live Aro- Salmonella typhimurium vaccines in the murine typhoid model.
473
Immunology. 90: 618-625.
474 475
8. MacLennan, C.A. Martin, L.B. and Micoli, F. 2014. Vaccines against
476
invasive Salmonella disease: current status and future directions. Hum.
477
Vaccin. Immunother. 10: 1478-1493.
478 479
9. Mastroeni, P. Harrison, J.A. Robinson, J.H. Clare, S. Khan, S. Maskell,
480
D.J. Dougan, G. and Hormaeche, C.E. 1998. Interleukin-12 is required for
481
control
of
the
growth
of
attenuated
aromatic-compound-dependent 23
482
salmonellae in BALB/c mice: role of gamma interferon and macrophage
483
activation. Infect. Immun. 66: 4767-4776.
484 485
10. Sinha, K. Mastroeni, P. Harrison, J. de Hormaeche, R.D. and Hormaeche,
486
C.E. 1997. Salmonella typhimurium aroA, htrA, and aroD htrA mutants cause
487
progressive infections in athymic (nu/nu) BALB/c mice. Infect. Immun. 65:
488
1566-1569.
489 490
11. Vazquez-Torres, A. Xu, Y. Jones-Carson, J. Holden, D.W. Lucia, S.M.
491
Dinauer, M.C. Mastroeni, P. and Fang, F.C. 2000. Salmonella pathogenicity
492
island 2-dependent evasion of the phagocyte NADPH oxidase. Science. 287:
493
1655-1658.
494 495 496
12. Gilchrist, J.J. MacLennan, C.A. Hill, A.V. 2015. Genetic susceptibility to invasive Salmonella disease. Nat. Rev. Immunol. 15: 452-463.
497 498
13. Mastroeni, P. Ugrinovic, S. Chandra, A. MacLennan, C. Doffinger, R.
499
Kumararatne, D. 2003. Resistance and susceptibility to Salmonella
500
infections: lessons from mice and patients with immunodeficiencies. Rev.
501
Med. Microbiol. 14: 53-62.
502 503
14. Khan, S.A. Stratford, R. Wu, T. Mckelvie, N. Bellaby, T. Hindle, Z. Sinha,
504
K.A. Eltze, S. Mastroeni, P. Pickard, D. Dougan, G. Chatfield, S.N. and
505
Brennan, F.R. 2003. Salmonella typhi and S. typhimurium derivatives 24
506
harbouring deletions in aromatic biosynthesis and Salmonella Pathogenicity
507
Island-2 (SPI-2) genes as vaccines and vectors. Vaccine. 21: 538-548.
508 509
15. Langridge, G.C. Phan, M.D. Turner, D.J. Perkins, T.T. Parts, L. Haase, J.
510
Charles, I. Maskell, D.J. Peters, S.E. Dougan, G. Wain, J. Parkhill, J. and
511
Turner, A.K. 2009. Simultaneous assay of every Salmonella Typhi gene
512
using one million transposon mutants. Genome. Res. 19: 2308-2316.
513 514
16. Hoiseth, S.K. and Stocker, B.A. 1981. Aromatic-dependent Salmonella
515
typhimurium are non-virulent and effective as live vaccines. Nature. 21: 238-
516
239.
517 518
17. Chaudhuri, R.R. Peters, S.E. Pleasance, S.J. Northern, H. Willers, C.
519
Paterson, G.K. Cone, D.B. Allen, A.G. Owen, P.J. Shalom, G. Stekel, D.J.
520
Charles, I.G. and Maskell, D.J. 2009. Comprehensive identification of
521
Salmonella enterica serovar typhiumrium genes required for infection of
522
BALB/c mice. PLoS. Pathog. 5(7): e1000529.
523 524
18. Chang, C.Y. and Cohen, S.N. 1978. Construction and characterization of
525
amplifiable multicopy DNA cloning vehicles derived from the P15A cryptic
526
miniplasmid. J. Bacteriol. 134: 1141-1156.
527
25
528
19. Bolivar, F. Rodriguez, R.L. Greene, P.J. Betlach, M.C. Heyneker, H.L. and
529
Boyer, H.W. 1977. Construction and characterization of new cloning vehicles.
530
II. A multipurpose cloning system. Gene. 2: 95-113.
531 532
20. Mo, E. Peters, S.E. Willers, C. Maskell, D.J. and Charles, I.G. 2006. Single,
533
double and triple mutants of Salmonella enterica serovar Typhimurium degP
534
(htrA), degQ (hhoA) and degS (hhoB) have diverse phenotypes on exposure
535
to elevated temperature and their growth in vivo is attenuated to different
536
extents. Microb. Pathog. 41: 174-182.
537 538
21. Dower, W.J. Miller, J.F. and Ragsdale, C.W. 1988. High efficiency
539
transformation of E. coli by high voltage electroporation. Nucl. Acids. Res. 16:
540
6127-6145.
541 542
22. Sambrook, J. and Russell, D. 2001. Molecular cloning: a laboratory manual,
543
3rd ed. Cold Spring Harbor Laboratory Press, Cold Spring Harbour, New
544
York.
545 546
23. Muyrers, J.P.P. Zhang, Y. Testa, G. and Stewart, A.F. 1999. Rapid
547
modification of bacterial artificial chromosomes by ET-recombination. Nucl.
548
Acids. Res. 27: 1555-1557.
549
26
550
24. Yu, D. Ellis, H.M. Lee, E.C. Jenkins, N.A. Copeland, N.G. and Court, D.L.
551
2000. An efficient recombination system for chromosome engineering in
552
Escherichia coli. Proc. Natl. Acad. Sci. USA. 97: 5978-5983.
553 554
25. Datsenko, K.A. and Wanner, B.L. 2000. One-step inactivation of
555
chromosomal genes in Escherichica coli K-12 using PCR products. Proc. Natl.
556
Acad. Sci. USA. 97: 6640-6645.
557 558
26. Grant, A.J. Restif, O. McKinley, T.J. Sheppard, M. Maskell, D.J. and
559
Mastroeni, P. 2008. Modelling within-host spatiotemporal dynamics of
560
invasive bacterial disease. PLoS. Biol. 6: e74.
561 562
27. Hautefort, I. Proença, M.J. and Hinton, J.C.D. 2003. Single-copy green
563
fluorescent protein gene fusions allow accurate measurement of Salmonella
564
gene expression in vitro and during infection of mammalian cells. Appl.
565
Environ. Microbiol. 69: 7480-7491.
566 567 568
28. Martin, M. 2011. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 17: 10.
569 570 571
29. Li, H. 2013. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv preprint arXiv:1303.3997.
572
27
573
30. Love, M.I. Huber, W. and Anders, S. 2014. Moderated estimation of fold
574
change and dispersion for RNA-seq data with DESeq2. Genome. Biol. 15:
575
550
576 577
31. Benjamini, Y. and Hochberg, Y. 1995. Controlling the false discovery rate: a
578
practical and powerful approach to multiple testing. J. Royal. Stat. Soc.
579
Series. B. 57 1: 289–300.
580 581
32. Chaudhuri, R.R. Morgan, E. Peters, S.E. Pleasance, S.J. Hudson, D.L.
582
Davies, H.M. Wang, J. van Diemen, P.M. Buckley, A.M. Bowen, A.J.
583
Pullinger, G.D. Turner, D.J. Langridge, G.C. Turner, A.K. Parkhill, J.
584
Charles, I.G. Maskell, D.J. and Stevens, M.P. 2013. Comprehensive
585
assignment of roles for Salmonella typhimurium genes in intestinal
586
colonization of food-producing animals. PLoS. Genet. 9(4): e1003456.
587 588
33. Paterson, G.K. Northern, H. Cone, D.B. Willers, C. Peters, S.E. and
589
Maskell, D.J. 2009. Deletion of tolA in Salmonella Typhimurium generates an
590
attenuated strain with vaccine potential. Microbiology. 155: 220-228.
591 592
34. Tacket, C.O. Sztein, M.B. Losonsky, G.A. Wasserman, S.S. Nataro, J.P.
593
Edelman, R. Pickard, D. Dougan, G. Chatfield, S.N. and Levine, M.M.
594
1997. Safety of live oral Salmonella typhi vaccine strains with deletions in htrA
595
and aroC aroD and immune response in humans. Infect. Immun. 65: 452-456.
596
28
597
35. Dougan, G. Chatfield, S. Pickard, D. Bester, J. O’Callaghan, D. and
598
Maskell, D. 1988. Construction and characterisation of vaccine strains of
599
Salmonella harbouring mutations in two different aro genes. J. Infec. Dis. 158:
600
1329-1335.
601 602
36. Hindle, Z. Chatfield, S.N. Phillimore, J. Bentley, M. Johnson, J.
603
Cosgrove, C.A. Ghaem-Maghami, M, Sexton, A. Khan, M. Brennan, F.R.
604
Everest, P. Wu, T. Pickard, D. Holden, D.W. Dougan, G. Griffin, G.E.
605
House, D. Santangelo, J.D. Khan, S.A. Shea, J.E. Feldman, and R.G.
606
Lewis, D.J. 2002. Characterization of Salmonella enterica derivatives
607
harbouring defined aroC and Salmonella pathogenicity island 2 type III
608
secretion system (ssaV) mutations by immunization of healthy volunteers.
609
Infect. Immun. 70: 3457-3467.
29
610
FIGURES AND TABLE LEGENDS
611 612
Figure 1: Relative fitness scores (log2-fold change) plotted against normalised mean
613
read counts. Fitness scores represent the log2-fold change in the sequence read
614
counts of output pools compared with the input pools. Output pools were obtained
615
from the livers of gp91-/-phox, C57BL/6 and BALBc mice (32). Significantly
616
attenuated mutants, represented as red closed circles, are a population of mutants
617
showing negative fitness scores and an adjusted P-value (Benjamini-Hochberg) of
618
P≤0.05. Mutants with fitness scores which were not significantly different from zero
619
are shown as black circles. The wide range of fitness scores exhibited in the
620
C57B/L6 mice is indicative of large scale stochastic loss ("random dropout") of
621
mutants during those infections.
622 623
Figure 2: Time to presentation of clinical signs in gp91-/-phox mice infected with
624
various doses of SL1344 cysE or SL1344 aroC. gp91-/-phox mice were infected i.v.
625
with Log10 5.85 CFU, Log10 4.85 CFU or Log10 3.85 CFU of either SL1344 aroC (grey
626
symbols) or SL1344 cysE (black symbols) and the time to clinical signs was
627
monitored.
628 629
Figure 3: WT C57BL/6 mice were infected i.v. (circles) or orally (triangles) with Log10
630
5.80 CFU or Log10 8.68 CFU, respectively, of SL1344 cysE, and bacterial counts in
631
the systemic organs (livers and spleens) were monitored at various times post
632
infection until the bacteria had cleared from the organs. 30
633 634
Figure 4: Vaccination study, immunising C57BL/6 WT mice with SL1344 cysE and
635
subsequent challenge with WT SL1344. A) WT C57BL/6 mice were vaccinated with
636
SL1344 cysE administered at either Log10 5.7 CFU i.v. or Log10 8.5 orally. At 91 days
637
post-immunisation the vaccinated mice and an unvaccinated group of mice (naive
638
control) were challenged orally with Log10 8.8 CFU of virulent SL1344. Bacterial
639
counts in the systemic organs (livers and spleens) were monitored at various times
640
post-challenge (N.B. one mouse from the orally vaccinated group had to be killed on
641
day 9 post challenge, since it was displaying signs of infection). B) WT C57BL/6
642
mice were vaccinated with SL1344 cysE administered at either Log10 5.8 CFU i.v. or
643
Log10 8.7 orally. At 91 days post-immunisation the vaccinated mice were challenged
644
orally with Log10 8.8 CFU or virulent SL1344. Bacterial counts in the systemic organs
645
(livers and spleens) were monitored at various times post challenge (N.B. one
646
mouse from both the orally vaccinated and the i.v. vaccinated group had to be killed
647
on day 9 post challenge, one mouse from the orally vaccinated group had to be killed
648
on day 18, day 19 and day 34 post challenge, since they were displaying signs of
649
infection).
650 651
Table 1: Bacterial strains and plasmids used in this study.
652 653
Table 2: CFU in the organs of C57BL/6 WT mice infected with defined mutants of S.
654
Typhimurium, the upper and lower TraDIS fitness scores (in C57BL/6 mice) and the
655
number of mutants in each gene identified by TraDIS (in C57BL/6 mice).
31
656 657
Table 3: CFU in the organs of gp91-/-phox mice infected with defined mutants of S.
658
Typhimurium, the upper and lower TraDIS fitness scores (in gp91-/-phox mice) and
659
the number of mutants in each gene identified by TraDIS (in gp91-/-phox mice).
32
660
-/-
gp91 phox
Mu
Tn5
C57BL/6
Mu
Tn5
BALBc
Mu
Tn5
661 662
Figure 1: Relative fitness scores (log2-fold change) plotted against normalised mean
663
read counts. Fitness scores represent the log2-fold change in the sequence read
664
counts of output pools compared with the input pools. Output pools were obtained
665
from the livers of gp91-/-phox, C57BL/6 and BALBc mice (32). Significantly 33
666
attenuated mutants, represented as red closed circles, are a population of mutants
667
showing negative fitness scores and an adjusted P-value (Benjamini-Hochberg) of
668
P≤0.05. Mutants with fitness scores which were not significantly different from zero
669
are shown as black circles. The wide range of fitness scores exhibited in the
670
C57B/L6 mice is indicative of large scale stochastic loss ("random dropout") of
671
mutants during those infections.
34
672
673 674
Figure 2: Time to presentation of clinical signs in gp91-/-phox mice infected with
675
various doses of SL1344 cysE or SL1344 aroC. gp91-/-phox mice were infected i.v.
676
with Log10 5.85 CFU, Log10 4.85 CFU or Log10 3.85 CFU of either SL1344 aroC (grey
677
symbols) or SL1344 cysE (black symbols) and the time to clinical signs was
678
monitored.
35
679
680 681
Figure 3: WT C57BL/6 mice were infected i.v. (circles) or orally (triangles) with Log10
682
5.80 CFU or Log10 8.68 CFU, respectively, of SL1344 cysE, and bacterial counts in
683
the systemic organs (livers and spleens) were monitored at various times post
684
infection until the bacteria had cleared from the organs.
36
685
A
686 687
B
688 689
Figure 4: Vaccination study, immunising C57BL/6 WT mice with SL1344 cysE and
690
subsequent challenge with WT SL1344. A) WT C57BL/6 mice were vaccinated with
691
SL1344 cysE administered at either Log10 5.7 CFU i.v. (circles) or Log10 8.5 orally 37
692
(triangles). At 91 days post-immunisation the vaccinated mice and an unvaccinated
693
group of mice (naive control; squares) were challenged orally with Log10 8.8 CFU of
694
virulent SL1344. Bacterial counts in the systemic organs (livers and spleens) were
695
monitored at various times post-challenge (N.B. one mouse from the orally
696
vaccinated group had to be killed on day 9 post challenge, since it was displaying
697
signs of infection). B) WT C57BL/6 mice were vaccinated with SL1344 cysE
698
administered at either Log10 5.8 CFU i.v. (circles) or Log10 8.7 orally (triangles). At 91
699
days post-immunisation the vaccinated mice were challenged orally with Log10 8.8
700
CFU or virulent SL1344. Bacterial counts in the systemic organs (livers and spleens)
701
were monitored at various times post challenge (N.B. one mouse from both the orally
702
vaccinated and the i.v. vaccinated group had to be killed on day 9 post challenge,
703
one mouse from the orally vaccinated group had to be killed on day 18, day 19 and
704
day 34 post challenge, since they were displaying signs of infection).
38
705 Strain or Plasmid S. Typhimurium SL1344
Relevant genotype or description
Source/Reference
Virulent in mice with an i.v. LD50 of