Transposon-insertion sequencing screens unveil requirements for EHEC growth and intestinal colonization.

Enterohemorrhagic Escherichia coli O157:H7 (EHEC) is an important food-borne pathogen that colonizes the colon. Transposon-insertion sequencing (TIS) was used to identify genes required for EHEC and E. coli K-12 growth in vitro and for EHEC growth in vivo in the infant rabbit colon. Surprisingly, many conserved loci contribute to EHEC's but not to K-12's growth in vitro. There was a restrictive bottleneck for EHEC colonization of the rabbit colon, which complicated identification of EHEC genes facilitating growth in vivo. Both a refined version of an existing analytic framework as well as PCA-based analysis were used to compensate for the effects of the infection bottleneck. These analyses confirmed that the EHEC LEE-encoded type III secretion apparatus is required for growth in vivo and revealed that only a few effectors are critical for in vivo fitness. Over 200 mutants not previously associated with EHEC survival/growth in vivo also appeared attenuated in vivo, and a subset of these putative in vivo fitness factors were validated. Some were found to contribute to efficient type-three secretion while others, including tatABC, oxyR, envC, acrAB, and cvpA, promote EHEC resistance to host-derived stresses. cvpA is also required for intestinal growth of several other enteric pathogens, and proved to be required for EHEC, Vibrio cholerae and Vibrio parahaemolyticus resistance to the bile salt deoxycholate, highlighting the important role of this previously uncharacterized protein in pathogen survival. Collectively, our findings provide a comprehensive framework for understanding EHEC growth in the intestine.

Microbiota Metabolism Promotes Synthesis of the Human Ah Receptor Agonist 2,8-Dihydroxyquinoline.

The aryl hydrocarbon receptor (AHR) is a major regulator of immune function within the gastrointestinal tract. Resident microbiota are capable of influencing AHR-dependent signaling pathways via production of an array of bioactive molecules that act as AHR agonists, such as indole or indole-3-aldehyde. Bacteria produce a number of quinoline derivatives, of which some function as quorum-sensing molecules. Thus, we screened relevant hydroxyquinoline derivatives for AHR activity using AHR responsive reporter cell lines. 2,8-Dihydroxyquinoline (2,8-DHQ) was identified as a species-specific AHR agonist that exhibits full AHR agonist activity in human cell lines, but only induces modest AHR activity in mouse cells. Additional dihydroxylated quinolines tested failed to activate the human AHR. Nanomolar concentrations of 2,8-DHQ significantly induced CYP1A1 expression and, upon cotreatment with cytokines, synergistically induced IL6 expression. Ligand binding competition studies subsequently confirmed 2,8-DHQ to be a human AHR ligand. Several dihydroxyquinolines were detected in human fecal samples, with concentrations of 2,8-DHQ ranging between 0 and 3.4 pmol/mg feces. Additionally, in mice the microbiota was necessary for the presence of DHQ in cecal contents. These results suggest that microbiota-derived 2,8-DHQ would contribute to AHR activation in the human gut, and thus participate in the protective and homeostatic effects observed with gastrointestinal AHR activation.

Genetic analysis of Vibrio parahaemolyticus intestinal colonization.

Vibrio parahaemolyticus is the most common cause of seafood-borne gastroenteritis worldwide and a blight on global aquaculture. This organism requires a horizontally acquired type III secretion system (T3SS2) to infect the small intestine, but knowledge of additional factors that underlie V. parahaemolyticus pathogenicity is limited. We used transposon-insertion sequencing to screen for genes that contribute to viability of V. parahaemolyticus in vitro and in the mammalian intestine. Our analysis enumerated and controlled for the host infection bottleneck, enabling robust assessment of genetic contributions to in vivo fitness. We identified genes that contribute to V. parahaemolyticus colonization of the intestine independent of known virulence mechanisms in addition to uncharacterized components of T3SS2. Our study revealed that toxR, an ancestral locus in Vibrio species, is required for V. parahaemolyticus fitness in vivo and for induction of T3SS2 gene expression. The regulatory mechanism by which V. parahaemolyticus ToxR activates expression of T3SS2 resembles Vibrio cholerae ToxR regulation of distinct virulence elements acquired via lateral gene transfer. Thus, disparate horizontally acquired virulence systems have been placed under the control of this ancestral transcription factor across independently evolved human pathogens.

Chemoproteomic profiling of host and pathogen enzymes active in cholera.

Activity-based protein profiling (ABPP) is a chemoproteomic tool for detecting active enzymes in complex biological systems. We used ABPP to identify secreted bacterial and host serine hydrolases that are active in animals infected with the cholera pathogen Vibrio cholerae. Four V. cholerae proteases were consistently active in infected rabbits, and one, VC0157 (renamed IvaP), was also active in human choleric stool. Inactivation of IvaP influenced the activity of other secreted V. cholerae and rabbit enzymes in vivo, and genetic disruption of all four proteases increased the abundance of intelectin, an intestinal lectin, and its binding to V. cholerae in infected rabbits. Intelectin also bound to other enteric bacterial pathogens, suggesting that it may constitute a previously unrecognized mechanism of bacterial surveillance in the intestine that is inhibited by pathogen-secreted proteases. Our work demonstrates the power of activity-based proteomics to reveal host-pathogen enzymatic dialog in an animal model of infection.

Ligand-mediated cytoplasmic retention of the Ah receptor inhibits macrophage-mediated acute inflammatory responses.

The Ah receptor (AHR) has been shown to exhibit both inflammatory and anti-inflammatory activity in a context-specific manner. In vivo macrophage-driven acute inflammation models were utilized here to test whether the selective Ah receptor modulator 1-allyl-7-trifluoromethyl-1H-indazol-3-yl]-4-methoxyphenol (SGA360) would reduce inflammation. Exposure to SGA360 was capable of significantly inhibiting lipopolysaccharide (LPS)-mediated endotoxic shock in a mouse model, both in terms of lethality and attenuating inflammatory signaling in tissues. Topical exposure to SGA360 was also able to mitigate joint edema in a monosodium urate (MSU) crystal gout mouse model. Inhibition was dependent on the expression of the high-affinity allelic AHR variant in both acute inflammation models. Upon peritoneal MSU crystal exposure SGA360 pretreatment inhibited neutrophil and macrophage migration into the peritoneum. RNA-seq analysis revealed that SGA360 attenuated the expression of numerous inflammatory genes and genes known to be directly regulated by AHR in thioglycolate-elicited primary peritoneal macrophages treated with LPS. In addition, expression of the high-affinity allelic AHR variant in cultured macrophages was necessary for SGA360-mediated repression of inflammatory gene expression. Mechanistic studies revealed that SGA360 failed to induce nuclear translocation of the AHR and actually enhanced cytoplasmic localization. LPS treatment of macrophages enhanced the occupancy of the AHR and p65 to the Ptgs2 promoter, whereas SGA360 attenuated occupancy. AHR ligand activity was detected in peritoneal exudates isolated from MSU-treated mice, thus suggesting that the anti-inflammatory activity of SGA360 is mediated at least in part through AHR antagonism of endogenous agonist activity. These results underscore an important role of the AHR in participating in acute inflammatory signaling and warrants further investigations into possible clinical applications.

Divergent Ah Receptor Ligand Selectivity during Hominin Evolution.

We have identified a fixed nonsynonymous sequence difference between humans (Val381; derived variant) and Neandertals (Ala381; ancestral variant) in the ligand-binding domain of the aryl hydrocarbon receptor (AHR) gene. In an exome sequence analysis of four Neandertal and Denisovan individuals compared with nine modern humans, there are only 90 total nucleotide sites genome-wide for which archaic hominins are fixed for the ancestral nonsynonymous variant and the modern humans are fixed for the derived variant. Of those sites, only 27, including Val381 in the AHR, also have no reported variability in the human dbSNP database, further suggesting that this highly conserved functional variant is a rare event. Functional analysis of the amino acid variant Ala381 within the AHR carried by Neandertals and nonhuman primates indicate enhanced polycyclic aromatic hydrocarbon (PAH) binding, DNA binding capacity, and AHR mediated transcriptional activity compared with the human AHR. Also relative to human AHR, the Neandertal AHR exhibited 150-1000 times greater sensitivity to induction of Cyp1a1 and Cyp1b1 expression by PAHs (e.g., benzo(a)pyrene). The resulting CYP1A1/CYP1B1 enzymes are responsible for PAH first pass metabolism, which can result in the generation of toxic intermediates and perhaps AHR-associated toxicities. In contrast, the human AHR retains the ancestral sensitivity observed in primates to nontoxic endogenous AHR ligands (e.g., indole, indoxyl sulfate). Our findings reveal that a functionally significant change in the AHR occurred uniquely in humans, relative to other primates, that would attenuate the response to many environmental pollutants, including chemicals present in smoke from fire use during cooking.

Assessment of Ah receptor transcriptional activity mediated by halogenated dibenzo-p-dioxins and dibenzofurans (PXDD/Fs) in human and mouse cell systems.

Polybrominated and mixed bromo/chloro dibenzo-p-dioxins and dibenzofurans (PXDD/Fs) are emerging environmental contaminants of concern. Thus far, an understanding of the toxicological behavior of these chemical species and their impact upon human health is incomplete. Here we utilized human and mouse hepatocellular carcinoma cell lines to examine the ability of differentially halogenated PXDD/F congeners to induce aryl hydrocarbon receptor (AHR)-mediated transcriptional activity. Dose-response experiments in reporter cell lines identified varied potencies among differentially halogenated PXDD/F isomers by comparison of EC50 values relative to the prototypical AHR agonist, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Brominated PXDD/F species displayed reduced capacity to activate the mouse AHR, compared to TCDD. Only BrCl3 dibenzo-p-dioxin was found to have a greater relative potency than TCDD to induce human AHR transcriptional activity. Human cells required ∼10-29-fold higher ligand concentrations to induce analogous AHR activity, relative to mouse cells. Decreased sensitivity of the human AHR to brominated dibenzofuran congeners directly corresponded to the number of bromine functional groups. Mixtures of these compounds exhibited an additive effect on AHR activation. The data also support the inclusion of mixed halogenated dibenzo-p-dioxins and dibenzofurans into routine environmental screening procedures as well as more thorough toxicological characterization of PXDD/Fs.

RpoS and quorum sensing control expression and polar localization of Vibrio cholerae chemotaxis cluster III proteins in vitro and in vivo.

The diarrheal pathogen Vibrio cholerae contains three gene clusters that encode chemotaxis-related proteins, but only cluster II appears to be required for chemotaxis. Here, we present the first characterization of V. cholerae's 'cluster III' chemotaxis system. We found that cluster III proteins assemble into foci at bacterial poles, like those formed by cluster II proteins, but the two systems assemble independently and do not colocalize. Cluster III proteins are expressed in vitro during stationary phase and in conjunction with growth arrest linked to carbon starvation. This expression, as well as expression in vivo in suckling rabbits, is dependent upon RpoS. V. cholerae's CAI-1 quorum sensing (QS) system is also required for cluster III expression in stationary phase and modulates its expression in vivo, but is not required for cluster III expression in response to carbon starvation. Surprisingly, even though the CAI-1 and AI-2 QS systems are thought to feed into the same signaling pathway, the AI-2 system inhibited cluster III gene expression, revealing that the outputs of the two QS systems are not always the same. The distinctions between genetic determinants of cluster III expression in vitro and in vivo highlight the distinctive nature of the in vivo environment.

Comparative RNA-Seq based dissection of the regulatory networks and environmental stimuli underlying Vibrio parahaemolyticus gene expression during infection.

Vibrio parahaemolyticus is the leading worldwide cause of seafood-associated gastroenteritis, yet little is known regarding its intraintestinal gene expression or physiology. To date, in vivo analyses have focused on identification and characterization of virulence factors--e.g. a crucial Type III secretion system (T3SS2)--rather than genome-wide analyses of in vivo biology. Here, we used RNA-Seq to profile V. parahaemolyticus gene expression in infected infant rabbits, which mimic human infection. Comparative transcriptomic analysis of V. parahaemolyticus isolated from rabbit intestines and from several laboratory conditions enabled identification of mRNAs and sRNAs induced during infection and of regulatory factors that likely control them. More than 12% of annotated V. parahaemolyticus genes are differentially expressed in the intestine, including the genes of T3SS2, which are likely induced by bile-mediated activation of the transcription factor VtrB. Our analyses also suggest that V. parahaemolyticus has access to glucose or other preferred carbon sources in vivo, but that iron is inconsistently available. The V. parahaemolyticus transcriptional response to in vivo growth is far more widespread than and largely distinct from that of V. cholerae, likely due to the distinct ways in which these diarrheal pathogens interact with and modulate the environment in the small intestine.

Indole and Tryptophan Metabolism: Endogenous and Dietary Routes to Ah Receptor Activation.

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor recognized for its role in xenobiotic metabolism. The physiologic function of AHR has expanded to include roles in immune regulation, organogenesis, mucosal barrier function, and the cell cycle. These functions are likely dependent upon ligand-mediated activation of the receptor. High-affinity ligands of AHR have been classically defined as xenobiotics, such as polychlorinated biphenyls and dioxins. Identification of endogenous AHR ligands is key to understanding the physiologic functions of this enigmatic receptor. Metabolic pathways targeting the amino acid tryptophan and indole can lead to a myriad of metabolites, some of which are AHR ligands. Many of these ligands exhibit species selective preferential binding to AHR. The discovery of specific tryptophan metabolites as AHR ligands may provide insight concerning where AHR is activated in an organism, such as at the site of inflammation and within the intestinal tract.

Short-term toxicity studies of thallium (I) sulfate administered in drinking water to Sprague Dawley rats and B6C3F1/N mice.

Thallium is a heavy metal that is known to induce a broad spectrum of adverse health effects in humans including alopecia, neurotoxicity, and mortality following high dose acute poisoning events. Widespread human exposure to thallium may occur via consumption of contaminated drinking water; limited toxicity data are available to evaluate the corresponding public health risk. To address this data gap, the Division of Translational Toxicology conducted short-term toxicity studies of a monovalent thallium salt, thallium (I) sulfate. Thallium (I) sulfate was administered via dosed drinking water to time-mated Sprague Dawley (Hsd:Sprague Dawley® SD®) rats (F0 dams) and their offspring (F1) from gestation day (GD) 6 until up to postnatal day (PND) 28 at concentrations of 0, 3.13, 6.25, 12.5, 25, or 50 mg/L, and adult male and female B6C3F1/N mice for up to 2 weeks at concentrations of 0, 6.25, 12.5, 25, 50, or 100 mg/L. Rat dams in the 50 mg/L exposure group were removed during gestation, and dams and offspring in the 25 mg/L exposure group were removed on or before PND 0 due to overt toxicity. Exposure to thallium (I) sulfate at concentrations ≤ 12.5 mg/L did not impact F0 dam body weights, maintenance of pregnancy, littering parameters, or F1 survival (PND 4-28). However, in F1 pups, exposure to 12.5 mg/L thallium (I) sulfate resulted in decreased body weight gains relative to control rats and onset of whole-body alopecia. Measurement of thallium concentrations in dam plasma, amniotic fluid, fetuses (GD 18), and pup plasma (PND 4) indicated marked maternal transfer of thallium to offspring during gestation and lactation. Mice exposed to 100 mg/L thallium (I) sulfate were removed early due to overt toxicity, and mice exposed to ≥ 25 mg/L exhibited exposure concentration-related decreases in body weight. Lowest-observed-effect levels of 12.5 mg/L (rats) and 25 mg/L (mice) were determined based on the increased incidence of clinical signs of alopecia in F1 rat pups and significantly decreased body weights for both rats and mice.

High-throughput fitness screening and transcriptomics identify a role for a type IV secretion system in the pathogenesis of Crohn's disease-associated Escherichia coli.

Adherent-invasive Escherichia coli (AIEC) are pathogenic bacteria frequently isolated from patients who have Crohn's disease (CD). Despite the phenotypic differences between AIEC and commensal E. coli, comparative genomic approaches have been unable to differentiate these two groups, making the identification of key virulence factors a challenge. Here, we conduct a high-resolution, in vivo genetic screen to map AIEC genes required for intestinal colonization of mice. In addition, we use in vivo RNA-sequencing to define the host-associated AIEC transcriptome. We identify diverse metabolic pathways required for efficient gut colonization by AIEC and show that a type IV secretion system (T4SS) is required to form biofilms on the surface of epithelial cells, thereby promoting AIEC persistence in the gut. E. coli isolated from CD patients are enriched for a T4SS, suggesting a possible connection to disease activity. Our findings establish the T4SS as a principal AIEC colonization factor and highlight the use of genome-wide screens in decoding the infection biology of CD-associated bacteria that otherwise lack a defined genetic signature.

Development of an Analytical Method for Quantitation of 2,2'-Dimorpholinodiethyl Ether (DMDEE) in Rat Plasma, Amniotic Fluid and Fetal Homogenate by UPLC-MS-MS for Determination of Gestational and Lactational Transfer in Rats.

2,2'-Dimorpholinodiethyl ether (DMDEE) is a specialty amine catalyst used in the production of flexible foams, adhesives and coatings. The potential for occupational exposure to DMDEE is high, but toxicity data are very limited. The objective of this work was to develop a method to quantitate DMDEE in biological matrices to assess gestational and lactational transfer of DMDEE in rats following exposure of dams The method used protein precipitation, followed by removal of phospholipids and analysis of supernatant by ultra-performance liquid chromatography-tandem mass spectrometry. Rat fetuses were homogenized in water prior to protein precipitation and delipidation procedures. The method was evaluated in male Sprague Dawley rat plasma over the concentration range 5 to 1000 ng/mL. The method was linear (r ≥ 0.99), accurate (mean relative error (RE) ≤ ±11.9%) and precise (relative standard deviation (RSD) ≤ 2.7%). The mean absolute recovery was 106%. The limit of detection was 0.262 ng/mL. Standards as high as ∼100,000 ng/mL could be successfully diluted into the calibration range (mean %RE = -14.9; %RSD = 0.5). The method was evaluated in Sprague Dawley rat dam plasma, post-natal day 4 pup plasma, gestational day (GD) 18 amniotic fluid and fetal homogenate (mean %RE ≤ ±11.9; %RSD ≤ 2.3). Concentrations of DMDEE in rat dam plasma, amniotic fluid and fetal homogenate stored for at least 29 days and in pup plasma for at least 18 days at -80°C were within 87.7 to 99.5% of Day 0 concentrations, demonstrating that DMDEE is stable in these matrices. The method was used to quantitate DMDEE in rat plasma, amniotic fluid and fetus samples from a dose range finding toxicology study in which dams were dosed via gavage with DMDEE from GD 6 at doses of 0 (control), 62.5 and 250 mg/kg/day. DMDEE concentration increased with the dose in all matrices examined. The concentration in GD 18 fetuses was almost 2-fold higher than GD 18 dams demonstrating gestational transfer of DMDEE. However, the concentration in post-natal day 4 pup plasma was more than an order of magnitude lower than corresponding dam plasma suggesting less potential for transfer of DMDEE from dams to pups via lactation. There was no significant difference in concentration for male and female pup plasma.

The aryl hydrocarbon receptor activates ceramide biosynthesis in mice contributing to hepatic lipogenesis.

Aryl hydrocarbon receptor (AHR) activation via 2,3,7,8-tetrachlorodibenzofuran (TCDF) induces the accumulation of hepatic lipids. Here we report that AHR activation by TCDF (24  µg/kg body weight given orally for five days) induced significant elevation of hepatic lipids including ceramides in mice, was associated with increased expression of key ceramide biosynthetic genes, and increased activity of their respective enzymes. Results from chromatin immunoprecipitation (ChIP), electrophoretic mobility shift assay (EMSA) and cell-based reporter luciferase assays indicated that AHR directly activated the serine palmitoyltransferase long chain base subunit 2 (Sptlc2, encodes serine palmitoyltransferase 2 (SPT2)) gene whose product catalyzes the initial rate-limiting step in de novo sphingolipid biosynthesis. Hepatic ceramide accumulation was further confirmed by mass spectrometry-based lipidomics. Taken together, our results revealed that AHR activation results in the up-regulation of Sptlc2, leading to ceramide accumulation, thus promoting lipogenesis, which can induce hepatic lipid accumulation.

Butylparaben multigenerational reproductive assessment by continuous breeding in Hsd:Sprague Dawley SD rats following dietary exposure.

Butylparaben (BP) is an antimicrobial agent utilized for decades as a preservative in numerous consumer products. The safety of parabens has recently come under scrutiny based on reports of estrogenic activity and suggested adverse effects upon the reproductive system. Due to the limited availability of studies that address the potential for BP exposure to induce reproductive toxicity, and clear evidence of human exposure, the National Toxicology Program conducted a multigenerational continuous breeding study to evaluate the impact of dietary BP-exposure at 0, 5000, 15,000, or 40,000 ppm on reproductive and developmental parameters in Hsd:Sprague Dawley SD rats. BP-exposure was not associated with adverse alterations of fertility, fecundity, pubertal attainment, or reproductive parameters in F0, F1, or F2 generations. Exposure-dependent increases in liver weights, and incidences of non-neoplastic liver lesions suggest the liver is a target organ of BP toxicity. No findings were observed that would support the purported mechanism of BP-induced endocrine disruption in perinatally-exposed rodents.

Unsupervised Learning Approach for Comparing Multiple Transposon Insertion Sequencing Studies.

Transposon insertion sequencing (TIS) is a widely used technique for conducting genome-scale forward genetic screens in bacteria. However, few methods enable comparison of TIS data across multiple replicates of a screen or across independent screens, including screens performed in different organisms. Here, we introduce a post hoc analytic framework, comparative TIS (CompTIS), which utilizes unsupervised learning to enable meta-analysis of multiple TIS data sets. CompTIS first implements screen-level principal-component analysis (PCA) and clustering to identify variation between the TIS screens. This initial screen-level analysis facilitates the selection of related screens for additional analyses, reveals the relatedness of complex environments based on growth phenotypes measured by TIS, and provides a useful quality control step. Subsequently, PCA is performed on genes to identify loci whose corresponding mutants lead to concordant/discordant phenotypes across all or in a subset of screens. We used CompTIS to analyze published intestinal colonization TIS data sets from two vibrio species. Gene-level analyses identified both pan-vibrio genes required for intestinal colonization and conserved genes that displayed species-specific requirements. CompTIS is applicable to virtually any combination of TIS screens and can be implemented without regard to either the number of screens or the methods used for upstream data analysis.IMPORTANCE Forward genetic screens are powerful tools for functional genomics. The comparison of similar forward genetic screens performed in different organisms enables the identification of genes with similar or different phenotypes across organisms. Transposon insertion sequencing is a widely used method for conducting genome-scale forward genetic screens in bacteria, yet few bioinformatic approaches have been developed to compare the results of screen replicates and different screens conducted across species or strains. Here, we used principal-component analysis (PCA) and hierarchical clustering, two unsupervised learning approaches, to analyze the relatedness of multiple in vivo screens of pathogenic vibrios. This analytic framework reveals both shared pan-vibrio requirements for intestinal colonization and strain-specific dependencies. Our findings suggest that PCA-based analytics will be a straightforward widely applicable approach for comparing diverse transposon insertion sequencing screens.

Using Tox21 High-Throughput Screening Assays for the Evaluation of Botanical and Dietary Supplements.

Introduction: Recent nationwide surveys found that natural products, including botanical dietary supplements, are used by ∼18% of adults. In many cases, there is a paucity of toxicological data available for these substances to allow for confident evaluations of product safety. The National Toxicology Program (NTP) has received numerous nominations from the public and federal agencies to study the toxicological effects of botanical dietary supplements. The NTP sought to evaluate the utility of in vitro quantitative high-throughput screening (qHTS) assays for toxicological assessment of botanical and dietary supplements. Materials and Methods: In brief, concentration-response assessments of 90 test substances, including 13 distinct botanical species, and individual purported active constituents were evaluated using a subset of the Tox21 qHTS testing panel. The screen included 20 different endpoints that covered a broad range of biologically relevant signaling pathways to detect test article effects upon endocrine activity, nuclear receptor signaling, stress response signaling, genotoxicity, and cell death signaling. Results and Discussion: Botanical dietary supplement extracts induced measurable and diverse activity. Elevated biological activity profiles were observed following treatments with individual chemical constituents relative to their associated botanical extract. The overall distribution of activity was comparable to activities exhibited by compounds present in the Tox21 10K chemical library. Conclusion: Botanical supplements did not exhibit minimal or idiosyncratic activities that would preclude the use of qHTS platforms as a feasible method to screen this class of compounds. However, there are still many considerations and further development required when attempting to use in vitro qHTS methods to characterize the safety profile of botanical/dietary supplements.

A live vaccine rapidly protects against cholera in an infant rabbit model.

Outbreaks of cholera, a rapidly fatal diarrheal disease, often spread explosively. The efficacy of reactive vaccination campaigns-deploying Vibrio cholerae vaccines during epidemics-is partially limited by the time required for vaccine recipients to develop adaptive immunity. We created HaitiV, a live attenuated cholera vaccine candidate, by deleting diarrheagenic factors from a recent clinical isolate of V. cholerae and incorporating safeguards against vaccine reversion. We demonstrate that administration of HaitiV 24 hours before lethal challenge with wild-type V. cholerae reduced intestinal colonization by the wild-type strain, slowed disease progression, and reduced mortality in an infant rabbit model of cholera. HaitiV-mediated protection required viable vaccine, and rapid protection kinetics are not consistent with development of adaptive immunity. These features suggest that HaitiV mediates probiotic-like protection from cholera, a mechanism that is not known to be elicited by traditional vaccines. Mathematical modeling indicates that an intervention that works at the speed of HaitiV-mediated protection could improve the public health impact of reactive vaccination.

Time-Resolved Transposon Insertion Sequencing Reveals Genome-Wide Fitness Dynamics during Infection.

Transposon insertion sequencing (TIS) is a powerful high-throughput genetic technique that is transforming functional genomics in prokaryotes, because it enables genome-wide mapping of the determinants of fitness. However, current approaches for analyzing TIS data assume that selective pressures are constant over time and thus do not yield information regarding changes in the genetic requirements for growth in dynamic environments (e.g., during infection). Here, we describe structured analysis of TIS data collected as a time series, termed pattern analysis of conditional essentiality (PACE). From a temporal series of TIS data, PACE derives a quantitative assessment of each mutant's fitness over the course of an experiment and identifies mutants with related fitness profiles. In so doing, PACE circumvents major limitations of existing methodologies, specifically the need for artificial effect size thresholds and enumeration of bacterial population expansion. We used PACE to analyze TIS samples of Edwardsiella piscicida (a fish pathogen) collected over a 2-week infection period from a natural host (the flatfish turbot). PACE uncovered more genes that affect E. piscicida's fitness in vivo than were detected using a cutoff at a terminal sampling point, and it identified subpopulations of mutants with distinct fitness profiles, one of which informed the design of new live vaccine candidates. Overall, PACE enables efficient mining of time series TIS data and enhances the power and sensitivity of TIS-based analyses.IMPORTANCE Transposon insertion sequencing (TIS) enables genome-wide mapping of the genetic determinants of fitness, typically based on observations at a single sampling point. Here, we move beyond analysis of endpoint TIS data to create a framework for analysis of time series TIS data, termed pattern analysis of conditional essentiality (PACE). We applied PACE to identify genes that contribute to colonization of a natural host by the fish pathogen Edwardsiella piscicida. PACE uncovered more genes that affect E. piscicida's fitness in vivo than were detected using a terminal sampling point, and its clustering of mutants with related fitness profiles informed design of new live vaccine candidates. PACE yields insights into patterns of fitness dynamics and circumvents major limitations of existing methodologies. Finally, the PACE method should be applicable to additional "omic" time series data, including screens based on clustered regularly interspaced short palindromic repeats with Cas9 (CRISPR/Cas9).