Host-selected mutations converging on a global regulator drive an adaptive leap towards symbiosis in bacteria

13 figures, 6 tables and 6 additional files


Host selection mechanisms that shape adaptive evolution by V. fischeri.

(A) Dorsal view of juvenile host E. scolopes (left) with box indicating the relative position of the ventrally situated symbiotic light organ. On the right, a schematic illustrating the stages at which host-imposed selection occurs during squid–V. fischeri symbiosis: host recruitment (mucus entrapment, aggregation at light organ pores), initiation of symbiosis (host defenses, including hemocyte engulfment and oxidative stress), and colonization and maintenance (nutrient provisioning, sanctioning of non-luminous cheaters, continued hemocyte patrolling, and daily purging). (B) Symbiont population growth modeled for a single passage on the basis of growth dynamics of V. fischeri ES114. Light-organ populations are initiated with as few as ~10 cells (Wollenberg and Ruby, 2009; Altura et al., 2013) or as much as 1% of the inoculum, but are reduced by 95% following venting of the light organ at dawn (every 24 hr) (Boettcher et al., 1996). Shaded areas represent night periods whereas light areas represent daylight, which induces the venting behavior. (C) Experimental evolution of V. fischeri under host selection as described in Schuster et al. (2010). Each ancestral V. fischeri population was prepared by recovering cells from five colonies, growing them to mid-log phase, and sub-culturing them into 100 mL filtered seawater at a concentration sufficient to colonize squid (≤20,000 CFU/mL). On day 1, ten un-colonized (non-luminous) juvenile squid were communally inoculated by overnight incubation, during which bacteria were subjected to the first host-selective bottleneck. Following venting of ~95% of the light organ population, the squid were separated into isolated lineages in individual wells of a 24-well polystyrene plate containing filtered sea water with intervening rows of squid from an un-inoculated control cohort, the aposymbiotc control (‘apo control’). Note that only two of the ten passage squid populations are shown. On days 2, 3, and 4, after venting, squid were rinsed and transferred into 2 mL fresh filtered seawater. Luminescence was measured at various intervals for each squid to monitor colonization and the absence of contamination in aposymbiotic control squid. On the fourth day, the squid and half of the ventate were frozen at −80°C to preserve bacteria, and the remaining 1 mL ventate was combined with 1 mL of fresh filtered seawater, and used to inoculate a new uncolonized 24-hr-old juvenile squid. The process continued for 15 squid only for those lineages in which squid were detectably luminous at 48 hr post inoculation.
Figure 2 with 1 supplement
Experimental evolution of Vibrio fischeri produced multiple alleles in the sensor kinase BinK.

(A). Phylogenetic relationship, symbiotic capacity, and mutations accrued during squid experimental evolution of ecologically diverse Vibrio fischeri strains. Strain relationships were inferred under maximum likelihood using whole genomes with RealPhy (Bertels et al., 2014) and with node supports calculated from 1,000 bootstraps. Graphic symbols for ecological niches represent the source of isolation. Intrinsic squid symbiotic capacities of the five experimentally evolved strains, as determined by the minimum inoculum concentration required for successful colonization of 90% of squid with a 3 hr (ES114, EM17, and WH1) or over-night (H905 and MJ11) inoculum, are represented by color spectrum. Consensus genomes for each of the parallel V. fischeri populations evolved through E. scolopes are shown on the right, with variants indicated by circles. Mutation details are shown in Table 2. The mutations that were selected in host-passaged populations improved symbiotic capacity rather than general vigor. (B) BinK mutations arising in squid-evolved populations of MJ11 occurred in the HAMP and HATPaseC domains. A homo-dimer structural model for BinK using TMPRed and hybrid histidine kinase domain modelling (Anantharaman and Aravind, 2000; Stewart and Chen, 2010) predicts that the accessory sensory Cache1 domain localizes to the periplasm whereas the remaining four functional domains (accessory HAMP, and conserved HisKA, HATPaseC, and REC phosphorelay domains) are cytoplasmic (shown as gray band). A position-specific scoring matrix (PSSM) analysis for each of the squid-evolved BinK positions indicates whether a given amino acid is more (positive) or less (negative) likely to be functionally neutral. Scores for the substitutions incurred at these sites are shown in bold. Please refer to Figure 2—figure supplement 1 for a phylogenetic assessment of BinK orthology across Aliivibrio and V. fischeri strains.
Figure 2—figure supplement 1
BinK orthology, conserved domains and squid-adapted binK alleles.

(A) Unrooted maximum-likelihood (ML) phylogeny of all of the hybrid histidine kinases identified in V. fischeri genomes. Gene families were phylogenetically annotated using Escherichia coli references where possible (not shown), otherwise using the ES114 locus tag.
Figure 3 with 2 supplements
Evolved binK alleles enhanced host colonization and conferred a fitness tradeoff in non-host environments.

(A) Symbiotic colonization efficiency of MJ11 and derivatives in squid. Percentage of squid colonized by culture-evolved (c1–c5) and squid-evolved (binK1- binK4, bolded isolates in Table 2) derivatives of MJ11. Three hours after a cohort of 10–20 squid were inoculated with 3000 CFU/mL of each MJ11 strain, the squid were separated into individual vials, and colonization percentages determined by detectable luminescence at 24 hr. Bars: 95% CI. (B) Growth rates of MJ11 and evolved strains during competition in laboratory culture. Average growth rates (realized Malthusian parameters) of ΔbinK, squid-evolved binK and culture-evolved flagellar mutants (fliA and fliP variants, see Table 2) following in vitro culture competition in minimal media with ancestral binK+ MJ11, estimated using CFU yields of each competitor recovered at regular intervals. Bars: 95% CI. The diagonal line indicates 1:1 growth. Please refer to Figure 3—figure supplement 1 for data on the competitive abilities of binK1 and binK3 during colonization. Please refer to Figure 3—figure supplement 2 for symbiotic yields (CFU) of ES114 and MJ11 strains after 24 and 48 hr.
Figure 3—figure supplement 1
Relative competitive ability of binK1 and binK3 variants to colonize squid.

In vivo competitions suggest no competitive advantage in squid colonization between evolved V. fischeri MJ11 variants carrying either HAMP or HATPaseC domain mutations. Relative competitive indices for binK1 and binK3 MJ11 variants (carrying HATPaseC and HAMP domain mutations, respectively) used to co-inoculate squid across a range of inoculum densities. Points above or below zero represent squid light organs that are dominated by bink3 or bink1, respectively.
Figure 3—figure supplement 2
Growth of strain ES114 and strain MJ11 and its binK variants in squid light organs 24 or 48 hr after inoculation.

Yields of symbionts determined by plating serial dilutions of squid homogenate as described previously (Whistler and Ruby, 2003). Note: the Y-axis is log-scaled. Bars: 95% CI.
Figure 4 with 1 supplement
Empirical and modeled estimates of selective advantage in evolving V. fischeri symbiont populations.

(A) Conceptual overview of symbiont population dynamics during growth in inoculum and following host colonization (black line), including daily host-imposed bottlenecks. (B) Comparison of the selection coefficients conferred by binK1 in strain MJ11EP2-4-1 (harboring no other mutations) relative to binK+ from co-inoculated squid light organs after 24 or 48 hr. The selective advantage (i.e., relative competitiveness) of the evolved allele increased significantly during this period from 1.1 to 1.8 (Fisher-Pitman permutation test, **p=0.0011). Each circle represents the selective advantage of each strain measured from the strain ratios recovered in an individual hatchling. Please refer to Figure 4—figure supplement 1 for the effect of starting binK1 frequencies and inoculum densities on estimates of selective advantage. (C) Modeled survival probabilities for new beneficial alleles arising in a growing symbiont population facing host-imposed bottlenecks. The gray shaded curves estimate the survival probability of new mutants following the subsequent population bottleneck, which depends on both the generation of growth in the inoculum or host in which they arise (x-axis) and the selective advantage (s) conferred by mutation (gray shading). Notably, beneficial variants that arise early in inoculum culture are likely to survive extinction at the subsequent bottleneck, and this probability of survival rapidly decreases even when conferring a large selective coefficient. On the basis of this model, for example, a mutation conferring a large selective advantage (s ~2) would have less than a 10% chance of surviving the subsequent colonization bottleneck if it arose during the tenth generation of inoculum growth (red line).
Figure 4—figure supplement 1
Estimates of the selective advantage of the binK1 allele during squid colonization across a range of starting frequencies and inoculum densities.

Comparison of selection coefficients conferred by binK1 in strain MJ11EP2-4-1 (‘Evo’) (harboring no other mutations) relative to binK+ (‘Anc’) from co-inoculated squid light organs. Each point represents the selective advantage of each strain measured from the strain ratios recovered in an individual hatchling. The estimated selective advantage conferred by the evolved binK1 allele was not influenced by starting frequency (A) (R2 = 0.025, pfrequency = 0.62), but it was marginally influenced by density (B) (R2 = 0.025, pdensity = 0.03), based on a multiple regression analysis.
Figure 5 with 3 supplements
Host-adapted binK1 improved initiation phenotypes through enhanced biofilm.

(A) V. fischeri MJ11 aggregate formation near light-organ ducts. Host tissue stained with CellTracker Orange. Symbionts carry GFP plasmids (pKV111) (Nyholm et al., 2000). Micrographs show representative V. fischeri aggregates following the dissection of 30 newly hatched animals incubated with each strain. Aggregates were visualized between 2 and 3hr after of inoculation using a Zeiss LSM 510 Meta laser-scanning confocal microscope. Please refer to Figure 5—figure supplement 1 for additional views of aggregate formation. (B) Biofilm production (crystal violet staining relative to MJ11) by wild-type MJ11(binK+), squid-adaptive binK1 and ∆binK variants in the presence of either empty vector (EV, pVSV105) (white fill), Syp biofilm repressor sypE (pCLD48) (hatched fill), or cellulose repressor binA (pRF2A3) (gray fill). n = 12–16 biological replicates. See Figure 5—figure supplement 2 for evidence of increased cellulose in binK variants, and Figure 5—figure supplement 3 for biofilm repressor schematic. Followed by influence of a sypK deletion on biofilm production of MJ11 and binK variants. n = 10 biological replicates. (C) Binomial mean of survival following exposure to hydrogen peroxide of wild-type MJ11(binK+), squid-adaptive binK1 and ∆binK variants in the presence of either empty vector (EV, pVSV105) (white fill), sypE (pCLD48) (hatched fill), or binA (pRF2A3) (gray fill). n = 20–50 biological replicates. Followed by influence of a sypK deletion (diagonal line overlay) on population survival of MJ11 and binK variants (color fill). n = 15–106 biological replicates. Error bars 95% CI. Significant p values (p<0.05) are indicated above each comparison. *p<2.2e-16. Although the effects of overexpression of binA and deletion of sypK on oxidative resistance in the ∆binK variant followed the same trends as these genes in binK1, the reductions were only marginally significant (p=0.051 and 0.15, respectively). Please refer to Figure 5—figure supplement 2 for transcriptomic evidence of reduced expression of two cellulose loci in the ∆binK mutant. A schematic of the impact of the BinA and SypE repressors on biofilm substrates is available as Figure 5—figure supplement 3.
Figure 5—figure supplement 1
In vivo aggregation behavioral changes conferred by evolved binK1 variant.

(A–B) Aggregation of ancestral (A) and evolved (B) MJ11 on host mucosal epithelium prior to colonization. Host tissue stained with CellTracker Orange. Symbionts carry GFP plasmids (pKV111) (Nyholm et al., 2000). Micrographs show representative V. fischeri aggregates following the dissection of 30 newly hatched animals incubated with each strain. Aggregates were visualized between 2 and 3 hr after inoculation using a Zeiss LSM 510 Meta laser scanning confocal microscope. Scale bars: 24 μm.
Figure 5—figure supplement 2
Transcriptional shifts associated with binK variants.

Transcriptomic differences between wild-type MJ11 (binK+), squid-adapted MJ11 binK1, and MJ11 ∆binK for the coding loci in the MJ11 genome as determined by RNA-Seq. Variants were sampled during early log growth (OD600 ~0.25) in rich media (SWTO) prior to detectable biofilm production from four biological replicates for each strain. Green indicates increased expression; red indicates reduced expression relative to mean expression per locus (i.e., read counts z-scaled relative to mean logCPM). The heat map only displays loci for which mean expression in a binK variant differed from that in MJ11 at a FDR significance threshold of 0.05 (Table 3). The colored labels refer to compounds whose metabolism, transport, or synthesis are affected by the expression of these genes. Genes involved in cellulose synthesis are indicated with arrows (VFMJ11_A1000- cellulose synthase operon C protein, and VFMJ11_A1007- cellulose synthase operon protein YhjU) and log fold change (logFC) relative to wild-type is indicated for binK1 and then ∆binK.
Figure 5—figure supplement 3
Schematic of regulation by the biofilm repressors SypE and BinA.

SypE represses Syp biofilm production post transcriptionally (Morris and Visick, 2013). BinA represses cellulose, but not Syp, biofilm formation by increasing phosphodiesterase activity (Bassis and Visick, 2010). Black-capped lines indicate negative regulation. Gray arrows indicate transcription/translation. Note that because binA is expressed from a syp locus promoter, activation of the syp locus leading to Syp production also leads to repression of cellulose.
Figure 6 with 1 supplement
Biofilm production by squid-adaptive binK1 variants mediates hemocyte evasion.

(A) Relative efficiency of squid hemocyte binding of GFP-labelled V. fischeri strains including: squid-native symbiont ES114, binK+ MJ11, ∆binK MJ11 (RF1A4), binK1 MJ11, and shellfish pathogen V. harveyi B392. (B) Relative efficiency of squid hemocyte binding of squid-native symbiont ES114 and squid-adapted bink1 MJ11 carrying the empty vector (pVSV104), sypE (pRF2A1) or binA (pRF2A4). N = 30–52 hemocytes quantified per strain. Error bars: 95% CI. Significant p-values (p<0.05) are indicated above each comparison. Please refer to Figure 6—figure supplement 1 for micrographs of Vibrio–hemocyte interactions.
Figure 6—figure supplement 1
In vitro response of squid hemocytes to wild, squid-evolved and mutant Vibrio.

The micrographs show examples of hemocyte-bound non-symbiotic (A: Vibrio harveyi), squid-symbiotic (B: V. fischeri ES114), squid-naive (C: V. fischeri MJ11 binK+) and squid-adapted (D: MJ11 binK1) cells. The mean number of GFP-labelled Vibrio cells bound by hemocytes was quantified relative to total bacterial count in a 60 µm radius using confocal microscopy at 63X magnification, following one hour of bacterial exposure. Squid hemocytes in red (CellTracker Orange), Vibrio in green (GFP). Scale bars: 12 μm.
Contribution of Syp and cellulose to improved squid colonization by binK variants.

(A) Colonization efficiency (% colonized squid at 24 hr) by wild-type MJ11 (binK+), squid-adaptive binK1 and ∆binK variants in the presence of empty vector (EV, pVSV105) (white fill), the Syp repressor sypE (pCLD48) (hatched fill), or the cellulose repressor binA (pRF2A3) (gray fill). n = 15–20 biological replicates. (B) Influence of a sypK deletion on colonization efficiency of MJ11 and binK variants. n = 31–52 biological replicates. Error bars: 95% CI. Significant p-values (p<0.05) are indicated above each comparison. *p<2.2e-16.
Host-adapted binK1 attenuates quorum-sensing regulation of luminescence.

(A) Supernatant concentrations (nM/OD600) of N-(3-oxohexanoyl) homoserine lactone (C6-HSL), as quantified against synthetic standards (Schaefer et al., 2000; Pearson et al., 1994; Duerkop et al., 2007) and corresponding luminescence (Lum/OD600) of 10 independent cultures each for wild-type MJ11, binK1 and ΔbinK derivatives during quorum-sensing induction of luminescence determined from cultures grown to early log (Average OD600 1.1, range 0.9–1.4,). (B) Average cell density as measured by absorbance (OD600), colony-forming units (CFU)/mL/OD600, N-(3-oxohexanoyl) homoserine lactone (C6) nM concentration, N-octanoyl homoserine lactone (C8) nM concentration, and luminescence (Lum)/1 mL culture for ten biological replicates of each variant relative to wild-type MJ11. Error bars: 95% CI. Significant p-values (p<0.05) are indicated above each comparison. *p<2.2e-16.
Effect of binK on squid colonization and biofilm production.

(A) Improvement in colonization by multi-copy in trans expression of the evolved binK1 allele and decreased colonization by expression of the ancestral binK+ allele. Colonization assessed by percentage of squid that are luminous after 24 hr. Error bars: 95% CI. N = 15–25. (B) Increased biofilm production resulting from in trans expression of the binK1 allele, and decreased biofilm production resulting from expression of the ancestral binK+. Comparisons of biofilm production in control-plasmids (pVSV105= EV) with that in multi-copy plasmids carrying binK suggest an inhibitory role for BinK in biofilm production, presumably alleviated by the dominance of the binK1 allele. Biofilm production was quantified by absorbance of crystal violet at A550. Background color depicts strain background in which multicopy plasmid effects were measured, mirroring those used throughout where blue is wild-type MJ11, green is the evolved binK1 variant and salmon is the ∆binK derivative. Error bars: 95% CI; non-overlap indicates significance. N = 7–8. Significant p-values (p<0.05) are indicated above each comparison. *p<0.05, **p<0.005, ***p<0.005.
Model of BinK regulation of traits adaptive during squid symbiosis.

Arrows originating from BinK point to characteristics that are activated or enhanced, and blocked lines point to those that are repressed or blocked by BinK. Hashed lines point to polysaccharides that contribute to biofilm.
Appendix 1—figure 1
Siderophore production in MJ11 and binK variants.

(A) Squid-native ES114, (B) squid-naïve MJ11 binK+ and (c) squid-evolved binK1 plated on CAS agar.
Appendix 2—figure 1
Metabolic shifts associated with binK variants.

Significantly differing metabolic responses to BIOLOG compounds for wild-type MJ11 (binK+), squid-adapted MJ11 binK1, MJ11 ∆binK and squid-native ES114. Responses to all tested compounds are reported in the Figure Supplement.
Appendix 2—figure 2
Metabolic profiles using BIOLOG phenotyping assays.

Plots enclosed by boxes indicate substrates that are significantly differentially metabolized across strains (listed in Table 2). X-axis represents time (0–48 hr); Y-axis represents metabolic activity as detected by BIOLOG redox (tetrazolium) dye absorbance (OD490).


Table 1

Strains and plasmids used in this study.
Strain nameDescription*Reference/source
Vibrio fischeri strains
ES114Isolated from Euprymna scolopes(Boettcher and Ruby, 1990)
MJ11Isolated from Monocentris japonica light-organ(Haygood et al., 1984)
EM17Isolated from Euprymna morseii light-organ(Ruby and Lee, 1998)
H905Isolated from Hawaiian plankton(Lee and Ruby, 1992)
WH1Isolated from Massachusetts plankton(Lee, 1994)
RF1A4MJ11 ∆binK::ermB; EmRThis study
RF1A5MJ11 ∆sypK::aphA1; KmRThis study
RF1A6MJ11 ∆binK::ermB ∆sypK::aphA1; EmR KmRThis study
RF1A7MJ11 binK1 ∆sypK::aphA1; KmRThis study
MJ11EP2-3-2MJ11 binK4This study
MJ11EP2-3-3MJ11 binK4This study
MJ11EP2-3-4MJ11 binK4This study
MJ11EP2-3-5MJ11 binK4This study
MJ11EP2-3-6MJ11 binK4This study
MJ11EP2-3-7MJ11 binK4This study
MJ11EP2-3-8MJ11 binK4This study
MJ11EP15-3-1MJ11 binK4This study
MJ11EP15-3-3MJ11 binK4This study
MJ11EP15-3-4MJ11 binK4This study
MJ11EP15-3-7MJ11 binK4This study
MJ11EP15-3-8MJ11 binK4This study
MJ11EP2-4-1MJ11 binK1This study
MJ11EP2-4-3MJ11 binK1This study
MJ11EP2-4-4MJ11 binK1This study
MJ11EP2-4-5MJ11 binK1This study
MJ11EP2-4-6MJ11 binK1This study
MJ11EP15-4-1MJ11 binK1 tadC1G593T(Schuster et al., 2010)
MJ11EP15-4-6MJ11 binK1This study
MJ11EP15-4-7MJ11 binK1This study
MJ11EP15-4-8MJ11 binK1This study
MJ11EP2-5-2MJ11 binK3This study
MJ11EP2-5-3MJ11 binK3This study
MJ11EP2-5-4MJ11 binK3This study
MJ11EP2-5-5MJ11 binK3This study
MJ11EP2-5-6MJ11 binK3This study
MJ11EP15-5-2MJ11 binK4This study
MJ11EP15-5-3MJ11 binK3This study
MJ11EP15-5-4MJ11 binK3This study
MJ11EP15-5-5MJ11 binK3This study
MJ11EP2-6-1MJ11 binK2This study
MJ11EP15-6-1MJ11 binK2(Schuster et al., 2010)
MJ11EP15-6-2MJ11 binK2This study
MJ11EP15-6-3MJ11 binK2This study
MJ11EP15-6-4MJ11 binK2This study
MJ11EP15-6-5MJ11 binK2This study
MJ11CE4-1MJ11 fliAG80DThis study
MJ11CE5-1MJ11 fliP∆476This study
Strain nameDescription*Reference/source
Escherichia coli strains
DH5αF recA1 endA1 hsdR17 supE44 thi-1 gyrA96
relA1Δ (argF-lacZYA) U169φ 80lacZΔM15λ
DH5αλpirsupE44 ΔlacU169lacZΔM15) recA1 endA1 hsdR17
thi-1 gyrA96 relA1; λpir phage lysogen
(Kolter and Helinski, 1978)
CC118λpirΔ(arg-leu) araD ΔlacX74 galE galK phoA20 thi-1 rpsE
rpoB argE(Am) recA1, lysogenized with λ pir dam dcm
(Martín-Mora et al., 2016)
NEB 10-betaΔ(ara-leu)7697 araD139 fhuA ΔlacX74 galK16 galE15 e14- Φ80dlacZΔM15 recA1
relA1 endA1 nupG rpsL (SmR) rph spoT1 Δ(mrr-hsdRMS-mcrBC)
New England Biolabs, Ipswich, MA
TOP10F- mcrA Δ(mrr-hsdRMS-mcrBC) Φ80lacZΔM15 ΔlacX74 recA1
araD139 Δ(ara-leu)7697 galU galK rpsL (SmR) endA1 nupG
Invitrogen, Carlsbad, CA
pCR2.1-TOPOCommercial cloning vector; ApR KmRInvitrogen, Carlsbad, CA
pVSV105Mobilizable vector; ChR(Dunn et al., 2006)
pRAD2E1pVSV105 carrying wild-type binK; ChRThis study
pRF2A2pVSV105 carrying binK1; ChRThis study
pCLD48pVSV105 carrying ES114 sypE; ChR(Hussa et al., 2008)
pRF2A3pVSV105 carrying MJ11 binA; ChRThis study
pVSV104Mobilizable vector; KmR(Stabb and Ruby, 2002)
pRF2A1pVSV104 carrying sypE; KmRThis study
pRF2A4pVSV104 carrying binA; KmRThis study
pKV111Mobilizable vector containing gfp; ChR(Nyholm et al., 2000)
pRF2B7pCR2.1-TOPO containing MJ11 ∆sypK::aph1 SOE fragment; KmRThis study
pVSV103Mobilizable vector containing lacZ; KmR(Dunn et al., 2006)
pCAW7B1pVSV103 containing lacZ∆147–1080 bp; KmRThis study
  1. *ApR, ampicillin resistance; ChR, chloramphenicol resistance; EmR, erythromycin resistance; KmR, kanamycin resistance;SmR streptomycin resistance.

  2. Experimentally evolved strains are designated ‘MJ11EP#-#-#’, where the first and second numbers after the ‘P’ designates the squid passage and population from which the strain was isolated, and the third number designates isolate number; strains derived from evolution in culture are designated ‘MJ11CE’.

Table 2

Summary of mutations detected following experimental evolution of V. fischeri using Illumina genome resequencing and targeted Sanger sequencing. For culture-evolved populations of V. fischeri MJ11, five isolates from each evolved population were combined to generate five metagenomes. For squid-evolved populations of MJ11, EM17, WH1 and H905, individual isolates were sequenced from lineages that ultimately survived 15 host passages. Isolates saved from early evolutionary time-points (host passage 2) are shown along with isolate genomes from the endpoint (host passage 15). Mean read depth and genome coverage for isolates analyzed with WGS are also provided.
AncestorEvolved Passage (EP)PopulationIsolateDetected mutationsIllumina sequencing statistics
binK (VFMJ11_A0397)tadC1
All other mutations
detectected by WGS
gene (locus);
mutation (reads)
Reads% Mapped to ancestorCoverage
allele/mutationMethod (reads)§ChIChII
MJ11211binK3/S311LWGS (35)375335299.5135.2118
MJ11213binK3/S311LWGS (32)371708899.6134.2113.5
MJ111514binK3/S311LWGS (17)171614499.546.842.5
MJ111531binK4/N292TWGS (42)303114998.9104.393.5
MJ111533binK4/N292TWGS (63)377771499.4114.6105.2
MJ111534binK4/N292TWGS (42)342021299.5106.497.1
MJ111537binK4/N292TWGS (41)330489199.590.382.5
MJ111538binK4/N292TWGS (63)294874399.685.581.2
MJ11241binK1/R537CWGS (62)2511256998478
MJ111541binK1/R537CWGS (131)G198V (85)412614999.4117.8106.1
MJ111546binK1/R537CWGS (61)G198V (55)226682199.260.852.5
MJ111547binK1/R537CWGS (89)G198V (93)307443799.69283.6
MJ111548binK1/R537CWGS (47)G198V (96)290297799.58477.5
MJ11252binK3/S311LWGS (26)377104899.6132.4123.7
MJ11253binK3/S311LWGS (46)259551899.684.283.7
MJ11254binK3/S311LWGS (20)178571399.560.657.2
MJ11255binK3/S311LWGS (62)364134699.6117.4113.1
MJ11256binK3/S311LWGS (81)412875199.6141.1134.8
MJ111552binK4/N292TWGS (89)443082399.1152.3138.4
MJ111553binK3/S311LWGS (10)324858099.38881.1
MJ111554binK3/S311LWGS (59)360938299.5106.897.1
MJ111555binK3/S311LWGS (28)291557099.587.482.6
MJ11261binK2/K482NWGS (104)474856999.1164.6147
MJ111561binK2/K482NWGS (75)276491099.483.275.5
MJ111562binK2/K482NWGS (63)324096899.28872.6
MJ111563binK2/K482NWGS (93)381436799.5108.1101.7
MJ111564binK2/K482NWGS (108)371463899.5121.485.7
MJ111565binK2/K482NWGS (90)300636299.485.572
MJ1115Culture4mgfliA (VF_1834); G80D (63)289416098.376.667.4
MJ1115Culture5mgfliP (VF_1842); ∆1 @ 476/870nt (110)557143997.9148.5132.1
WH11551NADH oxidase (VF_A0027); A402T (62)173210699.560.864.7
WH11552NADH oxidase (VF_A0027); A402T (61)173709599.461.964.9
WH11553NADH oxidase (VF_A0027); A402T (80)21948479660.863.4
EM171575icmF (VF_0992);S171N, (72)241428899.571.671.5
EM171593gdh2 (VF_1284); E732D (72)348418899.595.793.2
H9051511(∆37168 bp/25 genes)WGS (230)IlvY (VF_2529); M25I (233)764550894.2250.4222.1
H9051512(∆37168 bp/25 genes)WGS (167)IlvY (VF_2529); M25I (112)353111496.8117.5104.4
H9051513(∆37168 bp/25 genes)WGS (175)IlvY (VF_2529); M25I (97)359668997122.3109.1
H9051522∆16 bp@ 498/2595WGS (75)purR (VF_1572); N71T (60)281938797.691.479.6
H9051524∆16 bp@ 498/2595WGS (94)purR (VF_1572); N71T (52)299297896.9103.391.4
H9051525∆16 bp@ 498/2595WGS (90)purR (VF_1572); N71T (95)384483096.3123.6109
H905231tadF2 (VF_A0228); G21D (68)339361190.799.592.2
H9051531tadF2 (VF_A0228); G21D (140)797477391.5147.9143.9
H9051532T195IWGS (65)tadF2 (VF_A0228); G21D (28)198987595.565.458.2
H9051533tadF2 (VF_A0228); G21D (77)325389996.7103.894.4
H9051534tadF2 (VF_A0228); G21D (58)324274997.1103.394.7
H9051535tadF2 (VF_A0228); G21D (25)219077195.967.559
H9051541E43*WGS (102)665138592125.1130
H9051543E43*WGS (111)403237396.4135.9120.4
H9051544E43*WGS (187)612216895.8203.4179.4
H9051545E43*WGS (90)317781796.7100.890.6
H9051551∆1 bp @ 2325/2595ntWGS (113)716687090.4134.5130.9
H9051552∆1 bp @ 2325/2595ntWGS (94)370394696.7118.6108.3
H9051553∆1 bp @ 2325/2595ntWGS (66)282810297.498.690.4
H9051554∆1 bp @ 2325/2595ntWGS (109)472157597158.9143.8
H905261T195IWGS (105)tadF2 (VF_A0228); G21D (28)27436939483.373.6
H9051563T195IWGS (142)tadF2 (VF_A0228); G21D (49)559477197.5191.7175.3
H9051564T195IWGS (105)tadF2 (VF_A0228); G21D (37)336120696115.9101.4
  1. Individual characterized strain collection names assigned to isolates were derived from their ancestral lineage (e.g. MJ11), evolved passage (e.g. EP2), the population (e.g. 1), and isolate number (e.g. 1), which in the preceding example would give rise to strain collection name of MJ11EP2-1-1. Isolates in bold served as allelic binK representatives for further assays. mg: metagenome sequencing by pooling five isolates from a population.

  2. The presence of mutations was determined from Illumina short read (100PE) whole genome sequencing (WGS), by allele-specific PCR (PCR), and/or by locus PCR-amplification, followed by Sanger sequencing (SS). ‘–' indicates that no mutations were identified at this locus by breseq (Deatherage and Barrick, 2014) in this isolate using WGS. ‘n.d.’ indicates that the presence of mutations at this locus was not determined.

  3. §The number of reads supporting the mutation call from WGS data is provided. Mutations were called for sites with minimum coverage of 20 mappable reads. Mutations identified by Sanger sequencing (SS) of PCR-generated amplicons were confirmed from alignments of both forward and reverse reads. Coding genes reference V. fischeri ES114 locus tags.

Table 3

Genomes used in phylogenetic analyses. This table lists GenBank accessions for nucleotide genomes used in strain phylogeny and source for gene models used in hybrid histidine kinase phylogeny.
StrainNCBI accession/de novoProkka/NCBI gene models
Escherichia coliNC_000913NCBI
Aliivibrio wodanisLN554846-51NCBI
A. salmonicidaNC_011311–6NCBI
A. logeiNZ_AJYJ00000000Prokka
Vibrio furnissiiNC_016602, NC_016628NCBI
Vibrio parahaemolyticusNC_004603–5NCBI
Vibrio fischeri SR5NZ_AHIH00000000Prokka
Vibrio fischeri ES114NC_006840–2NCBI
Vibrio fischeri MJ11NC_011184–6NCBI
Vibrio fischeri EM17De novoProkka
Vibrio fischeri WH1De novoProkka
Vibrio fischeri ZF211AJYI01Prokka
Vibrio fischeri WH4De novoProkka
Vibrio fischeri SA1De novoProkka
Vibrio fischeri CG101De novoProkka
Vibrio fischeri H905De novoProkka
Vibrio fischeri PP3De novoProkka
Vibrio fischeri VLS2De novoProkka
Table 4

DNA oligonucleotide primers used in this study.
Primer namePrimer DNA sequence (5’−3’)Annealing temperatureSource
 A0397 4+ RGTCCCACCAAATTGACA53°CThis study
Appendix 1—table 1

Transcript expression differences between wild-type binK+ (ancestral MJ11) and binK mutants (∆binK and binK1) as detected by RNAseq under Fisher’s Exact test in edgeR. Positive fold-changes (logFC) indicate elevated expression in the wild-type relative to the indicated mutant; negative values indicate reduced expression in the wild-type relative to mutants. Loci with similar and significant expression changes in both ∆binK and binK1 relative to the wild-type are listed in bold. Only loci showing significant differences in transcript abundance compared with ancestral MJ11 are listed (FDR < 0.05).
Wildtype vs mutant ∆binKWildtype vs evolved binK1

Fold-change in expressionAverage
Corr. PFold-change in expressionAverage
Corr. P
Locus TaglogFClogCPMFDRlogFClogCPMFDRGene description
VFMJ11_0008−1.46.790.0027−0.35.940.6192Cystine-binding protein
VFMJ11_00132.77.870.00001.38.140.3143Small heat shock protein
VFMJ11_01951.07.430.00150.47.580.6346Co-chaperonin GroES
VFMJ11_02601.24.860.00030.74.980.2770Universal stress protein family protein
VFMJ11_0297−1.45.510.0063−0.64.810.4928Sulfite reductase NADPH flavoprotein alpha-component
VFMJ11_0307−−0.54.460.5758Sulfate adenylyltransferase subunit 2
VFMJ11_0421−0.86.960.0011−0.66.680.2201mutY AG-specific adenine glycosylase
VFMJ11_05781.19.320.00050.79.370.2725ATP-dependent chaperone ClpB
VFMJ11_0653−1.16.700.0000−0.25.950.8114Aminobenzoyl-glutamate transport protein
VFMJ11_06901.53.250.00010.63.530.5043urease accessory protein UreE
VFMJ11_06911.65.180.00000.45.550.5650Urease subunit alpha UreC
VFMJ11_11331.83.600.00000.63.940.3689Peptidase T pepT_1
VFMJ11_11371.18.560.00000.68.690.2725Glutamate decarboxylase
VFMJ11_11381.57.030.00000.47.340.4196TrkA domain integral membrane protein
VFMJ11_1253−0.89.480.0100− transporter
VFMJ11_1268−1.58.450.0000−0.57.570.3861Insulin-cleaving metalloproteinase outer membrane protein
VFMJ11_1269−−0.65.650.3122Thiol oxidoreductase
VFMJ11_1270−1.06.400.0062−0.45.870.4916Imelysin superfamily protein
VFMJ11_1305−1.43.640.0082− protein
VFMJ11_13171.28.360.00000.18.770.9186Hemin receptor
VFMJ11_1370−1.27.830.0055−1.47.820.05303-hydroxydecanoyl-ACP dehydratase fabA
VFMJ11_1398−0.66.430.0634−1.06.500.0335Na-dependent nucleoside transporter family protein
VFMJ11_1464−0.87.740.0088− U32
VFMJ11_1477−0.98.480.0065−1.08.430.1214Glycine betaine transporter
VFMJ11_1534−− RNA helicase RhlE
VFMJ11_1637−1.66.350.0000−1.15.790.1727Long-chain fatty acid transport protein
VFMJ11_1945−1.110.990.0001−1.010.790.1727Long-chain fatty acid transport protein
VFMJ11_2039−0.610.700.0678−1.010.900.0335Nitrate reductase catalytic subunit NapA
VFMJ11_2103−−0.68.780.2201Queuine tRNA-ribosyltransferase tgt
VFMJ11_21111.53.600.00021.03.740.2591Protein YgiW
VFMJ11_2127−−0.78.930.1727Peptidase U32
VFMJ11_21651.24.760.00620.54.970.4470DNA-binding transcriptional activator CadC
VFMJ11_22211.59.280.00000.19.800.9638Autonomous glycyl radical cofactor GrcA
VFMJ11_22230.76.550.00790.46.620.4482Homoserine kinase thrB
VFMJ11_22311.06.470.00020.56.580.3889Glutamate synthase subunit beta gltD
VFMJ11_2259−1.09.660.0011−0.79.320.2410IronIII ABC transporter periplasmic binding protein
VFMJ11_23941.15.830.00580.46.030.6194Succinylglutamic semialdehyde dehydrogenase astD
VFMJ11_24161.09.760.0032−0.310.320.6532Argininosuccinate synthase argG
VFMJ11_24561.310.160.00000.110.580.9472Fumarate reductase flavoprotein subunit frdA
VFMJ11_24571.48.350.00000.08.881.0000Fumarate reductase iron-sulfur subunit
VFMJ11_24581.76.080.00000.46.500.7063Fumarate reductase subunit C
VFMJ11_24591.46.680.0069− reductase subunit D
VFMJ11_2504−0.95.360.0048−0.54.930.39403-dehydroquinate dehydratase aroQ
VFMJ11_2505−1.08.610.0015−0.98.470.1727Acetyl-CoA carboxylase biotin carboxyl carrier protein subunit accB
VFMJ11_2506−1.210.800.0001−1.010.540.1384Acetyl-CoA carboxylase biotin carboxylase subunit accC
VFMJ11_26930.95.740.00830.55.800.3940Branched-chain amino acid aminotransferase ilvE
VFMJ11_2696−−1.28.520.1793Cold-shock DNA-binding domain
VFMJ11_A01061.77.960.00000.18.500.9472Anaerobic dimethyl sulfoxide reductase chain a
VFMJ11_A01111.07.200.00000.17.490.8523Outer membrane protein RomA
VFMJ11_A01510.97.290.0001−0.27.720.7881Putative tripeptide transporter permease tppB
VFMJ11_A02001.56.460.00000.76.680.3479L-lysine 6-monooxygenase
VFMJ11_A02011.67.970.00000.78.180.3687Ferric aerobactin receptor
VFMJ11_A02221.36.350.00000.66.540.3940Trimethylamine-n-oxide reductase 2
VFMJ11_A0224−−1.05.370.1748FhuE receptor
VFMJ11_A02801.23.750.00640.24.080.8702Methyl-accepting chemotaxis protein
VFMJ11_A03171.47.280.00000.27.650.6304Anaerobic ribonucleoside triphosphate reductase
VFMJ11_A0388−0.95.940.00490.45.090.4808Cyclic nucleotide binding protein
VFMJ11_A0389−2.09.660.00000.28.000.7303Sodium glucose cotransporter
VFMJ11_A0390−3.17.570.00000.44.860.5426UDP-glucose 4-epimerase
VFMJ11_A0391− uridylyltransferase
VFMJ11_A0393−1.56.890.00000.05.761.0000Aldose 1-epimerase
VFMJ11_A0394−6.39.930.00000.34.230.7720Transporter AcrB-D-F
VFMJ11_A0395− resistance periplasmic protein
VFMJ11_A0396−5.26.940.00000.02.500.9472Transcriptional regulator TetR family
VFMJ11_A0398− transcriptional regulator GalR
VFMJ11_A0487−1.47.990.0000− system N-acetylglucosamine-specific iibc component
VFMJ11_A0619−0.78.610.0099− multidrug transport system ATPase and permease component
VFMJ11_A0620−1.77.950.0000− antiporter
VFMJ11_A07551.23.780.00200.83.870.2319Restriction endonuclease
VFMJ11_A0885−1.47.510.0000−0.66.790.3687Cyclopropane-fatty-acyl-phospholipid synthase
VFMJ11_A0887−1.35.630.0000−0.44.890.6304Amine oxidase
VFMJ11_A0888−−0.54.350.4916Short chain dehydrogenase
VFMJ11_A0890−−0.25.300.9034Transcriptional activator ChrR
VFMJ11_A0891−1.25.700.0013−0.24.940.8635RNA polymerase sigma factor
VFMJ11_A0909−1.55.520.0000− receptor
VFMJ11_A1000−0.76.900.0075−0.36.490.5856Cellulose synthase operon C protein
VFMJ11_A1007−0.94.950.0040−0.54.550.4061Cellulose synthase operon protein YhjU
VFMJ11_A10381.03.890.00400.73.940.2725Alkanal monooxygenase beta chain
VFMJ11_A10391.44.020.00000.74.220.4313Alkanal monooxygenase alpha chain
VFMJ11_A10401.73.160.00051.13.330.2320Acyl transferase
VFMJ11_A10411.53.800.00030.74.010.3851Acyl-CoA reductase
VFMJ11_A10481.07.250.0069−0.47.870.6194Carboxypeptidase G2
VFMJ11_A1058−2.88.750.0000−1.77.690.0887Pts system fructose-specific eiibc component
VFMJ11_A1060−2.97.770.0000−1.66.520.0335Bifunctional PTS system fructose-specific transporter subunit IIA Hpr protein
VFMJ11_A1061−2.14.340.0000−1.23.440.1793DNA-binding transcriptional regulator FruR
VFMJ11_A12561.08.210.00000.48.380.4150Iron-regulated protein FrpC
Appendix 1—table 1—source data 1

Read counts estimated by RSEM for chromosome I transcript abundance (Appendix 1, Figure 5—figure supplement 2).
Appendix 1—table 1—source data 2

Read counts estimated by RSEM for chromosome II transcript abundance (Appendix 1, Figure 5—figure supplement 2).
Appendix 2—table 1

Metabolic convergence between squid native V. fischeri ES114 and squid-evolved binK1. The net changes in metabolic activity (as indicated by change in absorption of the Biolog tetrazolium redox dye) are shown for each V. fischeri strain after 48 hr of exposure to each substrate. Only substrates which induced significant (FDR < 0.05) differences across strains are listed. Metabolic changes in each strain relative to wild-type MJ11 binK+ are colored to indicate relatively increased or decreased activity. Of the 190 substrates tested, 44 substrates yielded significant differences across strain, including 39 which indicate congruent metabolic responses by ES114 and binK1 (Exact binomial test, p=1.405e-7).
Metabolic activity
(∆A550 over 48 hr)
Metabolic activity change
relative to remS+ MJ11
G10Methyl pyruvate0.0190.6770.4620.6390.00033.7022.7031.78+
H08Pyruvic acid0.0130.1870.2760.3950.00013.0919.7628.71+
E02m-Tartaric acid0.0260.4240.3380.4510.00015.1011.8416.13+
D02D-Aspartic acid0.0400.4590.3630.7350.00010.568.1617.55+
H02p-Hydroxy phenyl acetic acid0.0270.0630.0370.6360.0001.270.3622.11+
B06D-Gluconic acid0.0480.3000.3240.6280.0005.245.7512.06+
B09L-Lactic acid0.0290.0680.0400.6470.0001.300.3720.93+
C05Tween 200.0250.0650.0010.5170.0001.60−0.9519.54+
H09L-Galactonic acid—Lactone0.0890.2750.3790.4500.0002.103.274.06+
F8Mucic acid0.0260.0720.0440.0350.0001.780.680.36+
G11D-Malic acid0.0280.0620.0360.0290.0001.210.310.05+
F06Bromo succinic acid0.0330.0610.0370.0350.0000.820.110.04+
D06α-Keto-glutaric acid0.0450.0740.0430.0420.0000.65−0.04−0.08-
F9Glycolic acid0.0390.0620.0320.0400.0000.61−0.170.04+
E062-Hydroxy benzoic acid0.0890.0700.0800.0930.000−0.21−0.100.04+
E02Caproic acid0.2410.1010.0070.2150.000−0.58−0.97−0.11+
G02Tricarballylic acid0.0290.018−0.0080.0060.000−0.37−1.27−0.80+
D08α-Methyl-D- Galactoside−0.0110.0180.0080.0000.000−2.56−1.68−0.99+
Appendix 2—table 1—source data 1

Raw data for redox activity over 48 hr in BIOLOG plate PM1 (Appendix 2).
Appendix 2—table 1—source data 2

Raw data for redox activity over 48 hr in BIOLOG plate PM2A (Appendix 2).

Additional files

Source code 1

Statistical analysis of transcriptome changes in R (Appendix 1, Figure 5—figure supplement 2).
Source code 2

Statistical analysis of metabolic differences in BIOLOG assays in R (Appendix 2).
Appendix 1—table 1—source data 1

Read counts estimated by RSEM for chromosome I transcript abundance (Appendix 1, Figure 5—figure supplement 2).
Appendix 1—table 1—source data 2

Read counts estimated by RSEM for chromosome II transcript abundance (Appendix 1, Figure 5—figure supplement 2).
Appendix 2—table 1—source data 1

Raw data for redox activity over 48 hr in BIOLOG plate PM1 (Appendix 2).
Appendix 2—table 1—source data 2

Raw data for redox activity over 48 hr in BIOLOG plate PM2A (Appendix 2).

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. M Sabrina Pankey
  2. Randi L Foxall
  3. Ian M Ster
  4. Lauren A Perry
  5. Brian M Schuster
  6. Rachel A Donner
  7. Matthew Coyle
  8. Vaughn S Cooper
  9. Cheryl A Whistler
Host-selected mutations converging on a global regulator drive an adaptive leap towards symbiosis in bacteria
eLife 6:e24414.