Noise and crosstalk in two quorum-sensing inputs of Vibrio fischeri Pablo Delfino Prez, Joel T. Weiss, Stephen J. Hagen* Department of Physics, University of Florida, Gainesville FL 32611-8440, USA *Corresponding Author Email addresses: SJH email@example.com PDP firstname.lastname@example.org JTW email@example.com
Noise and crosstalk in Vibrio fischeri Abstract Background One of the puzzles in bacterial quorum sensing is understanding how an organism integr ates the information gained from multiple input signals. The marine bacterium Vibrio fischeri regulates its bioluminescence through a quorum sensing mechanism that receives input from t hree pheromone signals, including two acyl homoserine lactone (HSL) signals. While t he role of the 3-oxo-C6 homoserine lactone (3OC6HSL) signal in activating the lux genes has been extensively studied and modeled, the role of the C8 homoserine lactone (C8HSL) is less obvious, as it can either activate luminescence or block its activation. It remains unclear how c rosstalk between C8HSL and 3OC6HSL affects the information that the bacterium obtains through quorum sensing. Results We have used microfluidic methods to measure the response of individual V.fischeri cells to combinations of C8HSL and 3OC6HSL. By measuring the fluorescence of individual V.fischeri cells containing a chromosomal gfp -reporter for the lux genes, we study how combinations of exogenous HSLs affect both the population average and the cell-to-cell variability of lux activation levels. At the level of a population average, the crosstalk between the C8H SL and 3OC6HSL inputs is well-described by a competitive inhibition model. At the level of i ndividual cells, the heterogeneity in the lux response depends only on the average degree of activation, so that the noise in the output is not reduced by the presence of the second HSL signal. Overal l we find that the mutual information between the signal inputs and the lux output is less than one bit. A nonlinear correlation between fluorescence and bioluminescence outputs from lux leads to different noise properties for these reporters. Conclusions The lux genes in V.fischeri do not appear to distinguish between the two HSL inputs, and even with two signal inputs the regulation of lux is extremely noisy. Hence the role of crosstalk from
Noise and crosstalk in Vibrio fischeri the C8HSL input may not be to improve sensing precision, but rather to suppress the sensiti vity of the switch for as long as possible during colony growth.
Noise and crosstalk in Vibrio fischeri Background Quorum sensing is a mechanism of bacterial gene regulation that is based on the rel ease and detection of diffusible chemical signals. It is classically desc ribed as a population-sensing scheme: the bacteria release a pheromone (autoinducer) into their environment, a nd the accumulation of this autoinducer is an indicator of a high population density, triggerin g changes in phenotype. However it has become increasingly apparent that bacterial quor um sensing (QS) behaviors are often more complex than simple population-counting [1-4]. Many QS regul atory networks employ multiple receptors that receive signals from different auto inducers, forming interacting detectors that may act sequentially or in parallel to regul ate downstream genes. The design principles of these multi-input systems remain mysterious: One of the interesting puzzles in the study of QS is to understand what benefit or information an organism can gain from combining multiple autoinducer inputs, and how information from different inputs is processe d to generate a useful output [5, 6]. Here we investigate this question for two autoinducer inputs in Vibrio fischeri a g proteobacterium that uses QS to regulate bioluminescence as well as other behavior s that are important to colonization of its symbiotic host animal. These two autoinducers exhibi t a competitive or antagonistic interaction in regulating the lux operon that controls bioluminescence. We use microfluidic and single-cell methods to observe how combinati ons of the two autoinducer signals affect the bulk or average output of the QS network, as well a s the cell-to-cell variability in the activation of lux. We ask whether different combinations of signal inputs that produce the same average response across a population also produce the same response from individual V.fischeri, and therefore whether the lux system gains additional information from the presence of an additional signal. Bioluminescence in Vibrio fischeri is generated by the lux operon luxICDABEG which encodes the bacterial luciferase as well as enzymes for production of the lucifera se substrate . It is regulated by three QS channels  (Figure 1). Most well known is the LuxI/R me chanism. LuxI is the synthase of the autoinducer N -3-oxohexanoylLhomoserine lactone (3OC6HSL), which
Noise and crosstalk in Vibrio fischeri interacts with its cognate receptor LuxR to form a transcriptional act ivator for the lux operon. Threshhold concentrations (nM) of 3OC6HSL induce V.fischeri bioluminescence. The second QS system in V.fischeri is the AinS/R system (Figure 1). It uses the synthase AinS to produce the signal N -octanoylL -homoserine lactone (C8HSL). C8HSL interacts with its cognate receptor kinase AinR to initiate a phosphorelay signaling cascade (involving LuxU, LuxO, and a small RNA) that interrupts negative regulation of LitR, an activat or of LuxR. In addition to regulating luminescence via LuxR, AinS/R also regulates a numbe r of other behaviors, such as motility and acetate utilization , that are important to s uccessful colonization of the symbiotic host [10, 11]. The third QS system in V.fischeri differs from the first and second because it does not employ an acyl homoserine lactone (HSL) autoinducer. Instead the autoinducer is a furanosyl borate diester (AI2) that is synthesized by LuxS and detected by LuxP and LuQ. The signal fe eds into the same phosphorelay channel that detects the C8HSL autoinducer of the AinS/R syst em. AI2 influences luminescence (and presumably also noise in luminescence) via its dow nstream effect on LuxR expression. However the AI2 input makes a relatively small contribution to luminescence regulation and colonization , especially in comparison to the HS L autoinducers C8HSL and 3OC6HSL. Mutants deficient in production of the HSL autoinducers produce only very low, basal luminescence if any [13, 14]. Therefore, although AI2 may be import ant in interspecies communication , we have not included it in this study. Interestingly, the second autoinducer, C8HSL, also acts on luminescence through an additional route, where it bypasses the phosphorelay and interacts directly with LuxR to a ctivate lux expression (Figure 1). Therefore, while QS control of V.fischeri luminescence primarily occurs through the C8HSL and 3OC6HSL autoinducers, the effect of the two HSLs is not simply additive. The two routes of C8HSL action lead to a complex crosstalk between the AinS /R and LuxI/R systems. Generally, in a V.fischeri culture lacking 3OC6HSL, the addition of C8HSL induces bioluminescence. By contrast, in the presence of 3OC6HSL, the addition of C8HSL suppresses bioluminescence. Furthermore C8HSL appears to influence luminesce nce largely through direct interaction with the lux operon, rather than through the phosphorelay channel [13,
Noise and crosstalk in Vibrio fischeri 15]. These findings suggested a competitive inhibition model [11, 13] in which C8HSL modulates the bioluminescence by competing for the 3OC6HSL receptor LuxR (Fi gure 1): both C8HSL and 3OC6HSL are capable of binding to the receptor LuxR and activating transc ription of the lux genes , but the C8HSL LuxR complex is a less effective activator than the 3OC6HSL LuxR complex. The sensitivity of LuxR (as an activator of the lux genes) to C8HSL vs. 3OC6HSL is readily tunable through single-residue mutations , indicating that crosstalk could be m inimized if it impaired optimal regulation of bioluminescence. In fact, interaction between Ain S/R and LuxI/R not only exists but is strain-dependent, as the luminescence of V.fischeri mutants lacking the C8HSL synthase ( ainS mutants) behaves differently for strains derived from different symbiotic host animals. While the ainS mutation suppressed the luminescence of a strain extracted from the squid host Euprymna scolopes , the ainS mutation accelerated the induction of luminescence in a strain gathered from the fish host Monocentris japonicus . Furthermore, genomic analysis of V.fischeri strains derived from squid versus fish hosts showed that while most genes are highly conserved, luminescence activation and the regulatory targets of the LuxI/LuxR system exhibit significant divergence between strains [18, 19]. The strength of the crosstalk between C8HSL and 3OC6HSL, and the tuning of this inter action in strains that occupy different symbiotic environments, suggests that it is not i ncidental but rather that it provides an adaptive benefit. We should ask how a system of two HSLs, wo rking in opposition to each other, improves the regulation of bioluminescence. For example, Kuo et al noted that synthesis of C8HSL will delay the induction of luminescence early in gr owth, conserving the energy resources of the organism . Lupp and Ruby , noting that C8HS L regulates colonization factors in addition to luminescence, suggested that AinS /R and LuxI/R act sequentially so that the maximum induction of luminescence occurs after host col onization is initiated. It is still puzzling however that AinS/R should exhibit such strong cr osstalk with LuxI/R in regulating bioluminescence, since the same average delay in lumi nescence could presumably be achieved through a higher activation threshold for 3OC6HSL.
Noise and crosstalk in Vibrio fischeri Here we have investigated how combinations of C8HSL and 3OC6HSL signals affect the response of the lux operon at the individual cell level. Our recent study of the bioluminescent emission from individual V.fischeri found that the response of individual cells to defined concentrations of exogenous 3OC6HSL (alone) was extremely heterogeneous in the overall magnitude of luminescent emission and in the time scale for response to the HSL si gnal . Although the LuxI/R system exerts good control of the average bioluminescence of a V.fischeri population, it provides only weak control of individual cell behavior. Therefore we ask whet her the presence of two signals, C8HSL and 3OC6HSL, provides an additional dimension of control at the individual cell level, i.e. whether combinations of HSLs elicit a less noisy and more precise response in individual cells and therefore whether the second HSL signal i mproves the sensing precision of the QS system. Although we measured individual cell bioluminescence directly in our previous study, he re we use a fluorescent reporting strain (JB10) of V.fischeri containing a chromosomal gfp reporter of lux operon activation. We first show that the activation of JB10 GFP fluorescence in the presence of C8HSL and 3OC6HSL is described quantitatively by the competitive i nhibition model of Figure 1 The empirical parameters from the model define the average response and provide the basis for microfluidic single cell studies, in which we apply combina tions of C8HSL and 3OC6HSL autoinducers and measure the cell-to-cell variation in the lux response. From the observed distributions we can determine whether combinations of HSL inputs improve t he precision of individual cell response. Combining the single cell observations with the parametrized model we can also estimate the throughput of information from the H SL signal inputs to the overall lux output. Results Modeling the lux response Here our goal is to construct a mathematical representation of the interac tion between HSL signals by fitting bulk (well-plate) data to the competitive-binding model shown i n Figure 1 (see Methods ). The HSL-induced bioluminescence of V.fischeri strain JB10 observed in the wellplate assay is very similar to that observed in wild type strains [11, 13, 15]: C8HSL w eakly
Noise and crosstalk in Vibrio fischeri activates the bioluminescence in the absence of 3OC6HSL, while it represses bioluminescence in the presence of 3OC6HSL. 3OC6HSL consistently activates bioluminescence. Fi gure 2 shows the GFP fluorescence of JB10 as a function of HSL inputs. Since the gfp reporter is inserted into the lux operon of JB10, we expect the GFP fluorescence to correlate closely with bioluminescence over the full range of C8HSL and 3OC6HSL inputs. In principle we then ha ve the choice of fitting the model to either the luminescence L or the fluorescence F, as fluorescence and luminescence reporters are both regarded as reliable measures of ge ne expression . However Figure 3 shows that the correlation between the fluorescence and luminescence of JB10 is strong but it is not linear: at t = 10 hrs, when the bioluminescence response spans a dynamic range of ~3-4 decades, the fluorescence spans only ~ 20-fold. In fact, the empirica l relationship between fluorescence F ([3OC6HSL], [C8HSL]) and luminescence L ([3OC6HSL], [C8HSL]) is nearer to a power law L 1/2 a ( F F 0 ). Eqn. (1) Here a is a proportionality constant and F 0 describes a baseline fluorescence that is present even at the lowest activation levels where the luminescence is undetectable ( L 0). As the baseline F 0 grows with the population density (Figure 3) we interpret it as either a baseli ne expression of lux occuring independent of HSL activation  or autofluorescence of the cells. The power law in Eqn. (1) suggests that the luminescence intensity i s affected by the association equilibrium of LuxA and LuxB, which form the bacterial luciferase heterodimer : In a simple dimer association model, the concentration of enzymatically active luciferase ( L ) should scale as the product of the LuxA and LuxB concentrations. We also expect t hat luxA luxB and gfp should all be expressed at similar levels as they are all under the control of the same promoter. Therefore L should correlate with the square of the GFP concentration ( F-F 0 ), leading to Eqn. (1) From this perspective fluorescence F (relative to its background F 0 ) is preferable to L as a reporter of lux activation and should give a better fit to a physical model. Therefo re we based our analysis on the fluorescence F data. However the model also fits L 1/2 with virtually the same parameters as it fits F (Table 1).
Noise and crosstalk in Vibrio fischeri Fitting the competitive inhibition model to GFP fluorescence data at a range of optical densities early in growth (OD = 0.05-0.15 cm -1 ) gives a very satisfactory fit with minimal spread in the parameter values (Figure 2 and Table 1). Averaging the fit parameters obtaine d for F data over the OD range produces the surface F ([3OC6HSL],[C8HSL]) shown in Figure 2 Using this model surface we can predict the (average) lux response under any HSL condition and identify HSL combinations of interest for single-cell studies. Microfluidic studies of individual cells By loading V.fischeri cells into a three-channel microfluidic device (Figure 4) on a fluorescence microscope we can observe simultaneously three groups of cells subject to diffe rent combinations of HSL inputs and characterize the heterogeneity of their GFP response. The contour map of Figure 4 shows the HSL combinations used in four such experiments. Two experiments explore the heterogeneity along contours of near-constant lux activation while two experiments use signal combinations that cross contour lines. Experiment #1 applied three different C8HSL concentrations along with 100 nM 3OC6HS L to the three device channels at time t = 0. Initially all cells exhibit a weak fluorescence with a narrow distribution. In the absence of C8HSL, the response evolves over ~3-5 hrs to give a broad and distinctly non-Gaussian distribution that extends over an order of magnitude in GF P fluorescence, with a minority of cells becoming far brighter than the aver age. By contrast, when C8HSL was present at 500 nM or 1000 nM, the average cell fluorescence increased s lightly over time, but the distribution still remained narrower than in the absence of C8HSL. Sinc e the high C8HSL conditions ( i.e. C8HSL = 1000 nM and 3OC6HSL = 100 nM) induced virtually no lux response in our bulk experiments, we interpret the weak response as the baseline autofluorescence F 0 of Eqn. (1) Figure 5 shows that C8HSL does not simply reduce the average fluorescence, but rather reshapes the distribution by suppressing the development of the highly heterogeneous (noisy) activated response. Experiments #2 and #3 examined whether different combinations of C8HSL and 3OC6HSL that induce the same average response also elicit the same degree of heterogen eity. The selected
Noise and crosstalk in Vibrio fischeri n HSL conditions for #2 and #3 (Figure 6) follow two contours in Figure 4 corresponding to roughly 25% and 60% of full activation, respectively. Although the variance s 2 in GFP expression increases relative to the mean m at higher activation levels, we find that different combinations of HSLs that produce similar overall average fluorescence F also produce similar distributions. That is, for a given degree of activation, the distribution of individual c ell responses does not appear sensitive to the particular combination of HSL signals that induced that response. We also find that the lux response develops on the same time scale in all three channels, regardless of the relative proportions of 3OC6HSL and C8HSL. This is cons istent with prior findings that C8HSL signaling through the AinR/S route is not essential for luminescence . C8HSL acts on the same time scale as 3OC6HSL as it primarily regula tes luminescence through a direct association with LuxR. Mutual information between inputs and output The cell fluorescence histograms can be expressed as probability distribut ions P ( F | ([C8HSL], [3OC6HSL])) for the cell fluorescence F given the HSL inputs. The finding that these distributions are a function of the overall lux activation F ([C8 HSL], [3OC6 HSL]) rather than depending on the C8HSL and 3OC6HSL levels separately (Figures 6 and 7), strongly su ggests that the additional (C8HSL) input does not induce a response from the lux genes of the individual cell that 3OC6HSL alone could not extract. In this sense the lux system does not gain additional information by employing two HSL autoinducers. (This would not be true of other regul atory targets of the phosphorelay controlled by AinS/R and LuxS/P/Q, which is not regulate d by 3OC6HSL.) We can quantify the information that is gained by calculating the mut ual information, which measures the regulatory precision, or the number of practically distinguishable input/output states, of this regulatory system. One may think of mutua l information as measuring the amount by which the uncertainty in the (GFP) output is reduc ed by knowledge of the HSL inputs (and vice versa) [23-25]. Mutual information has been illuminating in recent studies of information throughput in other single-cell che mical sensing systems such as the multiple-autoinducer QS scheme of V.harveyi  and the chemotaxis of Dictyostelium discoideum .
Noise and crosstalk in Vibrio fischeri nn Calculating the mutual information between the cells environment, as defined by the HSL inputs, and the lux output requires a mathematical model for P ( F ). For higher activation levels we found that P ( F ) is too broad to be satisfactorily represented as a Gaussian, although it is reasonably well-described by a gamma distribution : P ( F ) dF = dF F (n-1) exp( F / b ) / G (n) bn Eqn. (2) The gamma distribution depends on two independent parameters, n and b, that depend on the variance s 2 and mean m of the distribution. The dimensionless parameter n = m 2 / s 2 completely determines the shape of P ( F ). The Fano factor b = s 2 / m does not affect the shape of the distribution, but it scales the horizontal axis and normalization according to the uni ts of measurement of F ( e.g. protein copy number, GFP fluorescence counts, etc.). Eqn. (2) is an appealing model for P ( F ) because it arises naturally in an intrinsic noise model, where mRNAs from one gene are synthesized in a Poisson process, with each mRNA leading to a burs t of protein expression [27, 28]. Therefore we modeled the observed heterogeneity P ( F ) as a gamma distribution that has the same relationship between the ratio s / m and the average activation level as we observed in our experiments. For this calculation we also define two new coordinates, X and Y to represent the progress of the two LuxR-HSL binding equilibria: X = [3OC6] n /([3OC6] n + k 2 n ) Y = [C8] m /([C8] m + k 1 m ) (See Methods ). X and Y are both confined to the interval 0 1. When plotted in terms of these coordinates, the response surface F (X,Y) generated by the competitive model has a curved shape (Figure 7) that highlights the asymmetry in the lux response to Y ( i.e. to [C8 HSL]) versus X ([3OC6 HSL]): the sensitivity to Y depends on the value of X The models for F ( X,Y ) and P ( F ) together provide an accurate mathematical representation of our bulk and single-cell data, from which we can calculate the mutual information I ( F ; ( X,Y ))
Noise and crosstalk in Vibrio fischeri n between the input combination ( X,Y ) and the output F. The calculation, described in Methods, leads to a very modest I ( F ; ( X,Y )) 0.53 bits. By contrast, a simple noiseless ON / OFF switch would transmit precisely one bit. Therefore, while the LuxI/LuxR system pr ovides precise control of the average lux response of a population over a range of HSL input concentrations, that response activates so gradually with respect to the input levels, and with such a heterogeneous output, that the individual cell cannot be said to exhibit a clean switch ing between the OFF and ON states of its lux response. Discussion The functional advantages of multi-input quorum sensing (QS) systems are not gene rally understood, even though many such architectures are known [5, 6]. Because our previous study of V.fischeri bioluminescence showed an extremely noisy response to 3OC6HSL alone, we investigated whether the presence of the C8HSL signal in addition to 3OC6HSL aff ects the noise performance of the LuxI/R system. To accomplish this we first constructed a data-based model for competitive inhibition of 3OC6HSL by C8HSL and used this model to draw a contour map of lux activation by the two signals. We then used microfluidic devices to control the che mical environment while we measured the heterogeneity of lux response among individual cells. The microfluidic flow chamber allows for extended observation of individual cells as w ell as precise definition of the exogenous HSL levels, eliminating the possibility of QS circ uit autoactivation. V.fischeri luminescence is a model system in quorum regulation, and accordingly it has been the subject of mathematical studies of both deterministic and stochastic QS behavior [29-36]. Most studies have emphasized dynamics and steady states of the lux system and 3OC6HSL alone, without considering the role of AinS/R. One exception is the work of Kuttler and Hense  who presented a detailed dynamical model for the combined ain and lux signaling pathway as outlined by Lupp et al  and others. The resulting system of ordinary differential equations displays a variety of possible stationary states and dynamics, although thes e outcomes depend on some poorly-known microscopic parameters that characterize transport, producti on and degradation of the HSLs, synthase production, and the kinetics of HSL-LuxR complex
Noise and crosstalk in Vibrio fischeri n formation. Interestingly however, those authors found that the experimentally-obs erved differences in the effect of the ainS mutation in squid-derived  versus fish-derived  strains of V.fischeri could arise from relatively minor differences in parameters describing the C8HSL and 3OC6HSL competition for LuxR. Specifically, the relative affinity of the C8HSLLuxR vs. 3OC6HSL-LuxR complex for the lux box, and the relative strength of lux activation by C8HSL-LuxR vs 3OC6HSL-LuxR, can determine whether an ainS mutant will be dark (as in ) or show accelerated luminescence response (as in ). These parameter s influence dynamical effects such as the role of ainS in the autoactivation of the LuxI system. Because autoactivation is not possible in the microfluidic chamber, and because C8HS L primarily affects bioluminescence through its direct interaction with L uxR, we set aside many of these dynamical complexities in modeling our population-averaged data on the ef fect of different HSL combinations. Instead we used the four-parameter competitive inhibition model (Figure 1) described in Methods The model provides a satisfactory fit to experimental data collected on the fluorescent reporting strain (JB10), although the same model will not describe the bioluminescence as accurately (see Eqn. (1)) From the fit parameters we calculated the fluorescence response surface F ((C8HSL), (3OC6HSL)), which we used as the basis for microfluidic studies of the heterogeneity in lux activation in the presence of multiple autoinducers. We studied the heterogeneity in the lux response along contours (of constant fluorescence F ) or along slices (of varying F ) in the [C8HSL],[3OC6HSL] plane. In general the heterogeneity varies with the average degree of lux expression, with the coefficient of variation ( s / m ) trending downward as activation increases (Figure 7). The variance in individua l cell responses appears much less sensitive to the particular combination of C8HSL and 3OC6H SL inputs than to the overall degree of activation. The relative proportions of 3OC6HSL and C8HSL do not appear to influence the time scale for development of the lux response. These findings are fully consistent with the competitive inhibition model, if the noise in L uxI/R is controlled by a low copy number of the HSL receptor LuxR . One of the puzzles in our previous study of bioluminescence was that the noise in the 3OC6H SL response was quite large i.e. CV = s / m 1, especially as a GFP-reporter study found significantly lower noise levels, CV ~ 0.15-0.4 in V.harveyi QS . Here it is interesting the
Noise and crosstalk in Vibrio fischeri n empirical Eqn. (1) predicts different noise levels for luminescence vs. fluorescence reporters of lux If F = F-F 0 is the level of fluorescence activation above threshhold, then Eqn. (1) predicts d L / L 2 d F / F so that the coefficient of variation (CV = s / m ) should be roughly twice as large for the luminescence as for the fluorescence. Therefore the heterogeneity seen in the single cell fluorescence, which is characterized by CV 0.4-0.6, is fully consistent with our single cell bioluminescence data. Nevertheless the regulation of lux still appears noisier in V.fischeri than in V.harveyi [37, 38]. The antagonistic interaction between HSL signals in V.fischeri is an intriguing contrast to the additive signaling found in V. harveyi bioluminescence, a model system for QS regulation that lacks the LuxI/LuxR mechanism. In V.harveyi three distinct autoinducers are detected by three membrane-bound histidine kinases that feed into the same phosphotransferase LuxU. LuxU controls a phosphorelay cascade that regulates the single output LuxR VH (unrelated to V.fischeri LuxR). The fact that the circuit merges inputs from both an HSL and a furanosyl borat e diester (AI2) autoinducer suggests that it senses both intra-species and interspecies QS signals, possibly functioning as a coincidence detector for the input signals . This could increas e the systems resistance to crosstalk from other bacterial QS signals or prevent it from responding in certain habitats. Moreover the autoinducer response is additive and symmetric in the sense t hat all three receptors contribute positively and in parallel to the output, with two having equal kina se activities, so that the output responds in the same way to activation of each recepto r [5, 25, 37]. Those authors suggested that equal sensitivity to each autoinducer benefits the orga nism by providing a graded, sequential activation of bioluminescence during growth. Lupp et al proposed a similar, sequential interpretation  for the role of C8HSL and 3OC6HSL in V.fischeri. They suggested that C8HSL acts first to stimulate luminescence at intermediate cell densities (as in cultures), activating luxR expressi on through the AinS/R route and also interacting directly with LuxR. At higher cell densities (as late r during colonization) luxR remains activated by C8HSL but 3OC6HSL accumulates to sufficient concentra tions to interact with LuxR and activate lux. It would indeed be remarkable if both V.fischeri and V.harveyi used multiple autoinducers to achieve sequential activation of lux yet only V.fischeri did so by using antagonistic signal inputs.
Noise and crosstalk in Vibrio fischeri n Alternatively, Kuo et al suggested that the suppressing role of C8HSL served a different function, conserving the energy resources of the organism by delaying the inducti on of luminescence early in V.fischeri growth . One puzzle however is that the same delayed outcome could presumably be achieved by setting a higher threshold for induction by 3OC6HSL, making the second signal unnecessary. We cannot fully interpret bioluminescence regulation in V.fischeri without considering its symbiotic context, as the full QS network that regulates both bioluminescence and host colonization receives input from many environmental factors [22, 40]. However we ca n still ask which properties of the V.fischeri LuxI/LuxR system could make an antagonistic interaction between 3OC6HSL and C8HSL advantageous. We used our experimental and modeling result s to quantify the signal-transmission property of the two-HSL system. We calcula ted the mutual information between lux output and the signal inputs [23-25] by modeling the populationaveraged lux activation F ( X,Y ) with the competitive inhibition model (where X and Y are scaled variables corresponding to the relative saturation of LuxR by 3OC6HSL and C8HSL respectively), and then modeling the noise in F by a gamma distribution that captures the coefficient of variation observed in our single-cell experiments (Figure 7 a nd Methods ) In the absence of any other information about signal inputs i.e. using the simplest assumption that all input combinations ( X,Y ) are equally likely a priori the calculation leads to a surprisingly low estimate for the mutual information, I ( Z ,( X,Y )) 0.53 bits. Even with its two signal inputs, the output F ( X,Y ) of the LuxI/LuxR system transmits less information about its inputs than would a simple ON/OFF switch. By contrast, Mehta et al. estimated ~1.2-1.7 bits of mutual information between the output and two inputs (AI1 and AI2) of the phosphorelay system in V.harveyi QS. The noisy performance (CV ~ 0.5) and gradual switching of LuxI/R significa ntly degrades its sensing capability, in comparison to the V.harveyi circuit. Therefore we find no indication that the second (C8HSL) autoinducer enhances the prec ision of signal response in the V.fischeri LuxI/R system. However the poor information throughput of this system does suggest a different perspective on the idea  that C8HSL c onserves energy resources by delaying induction of lux at low population densities. At the population-average
Noise and crosstalk in Vibrio fischeri n level, the same delay in activation could be achieved by raising the activation t hreshold for the 3OC6HSL signal. However at the single cell level, the presence of any 3OC6HSL induces a highly heterogeneous response, with some cells luminescing much more brightly than a verage. Thus, even if the threshold is set very high, a few cells will waste energy by e mitting light during early growth. One benefit of producing a small concentration of C8HSL is that it collapses the bioluminescence distribution, suppressing the most active emitters and conserving m etabolic energy. Simultaneous synthesis of C8HSL and 3OC6HSL may therefore reduce lumine scent output by virtually all cells, at least until 3OC6HSL attains high concentrat ion. In this sense the crosstalk from the C8HSL signal does not improve the environment-sensing precision of LuxI/R at steady state, but it may tend to compensate for the noisy performance of the LuxI/R switch by suppressing the switch for as long as possible during growth and colonization. It woul d be intriguing to see if dynamical models that accurately capture the noise in the c ircuit and the temporal accumulation of HSL can characterize this behavior quantitatively Conclusions Although multiple-input quorum-sensing systems are widespread in the microbial w orld, the mechanisms by which they combine and process information from parallel sig nal inputs are in general poorly understood. One of the intriguing properties of the V.fischeri QS network is that it employs two autoinducer signals that can act competitively or antagonis tically in regulating the lux genes. In order to understand the possible advantages of this competitive interacti on we have studied the response of individual V.fischeri to combinations of HSL signals The populationaveraged, steady state activation of lux by the two HSL signals is readily described by a quantitative, competitive inhibition model. Our measurements of lux activation in individual cells show a noisy response, with the LuxI/R circuit conveying less than one bi t of mutual information between its HSL signal inputs and its lux output. Further the data provide no indication that either the dynamics of the lux response or the heterogeneity in that response are sensitive to different combinations of signals that generate the same populati on-averaged output. In this sense the second HSL signal input appears to provide little if any additiona l information to the lux system. These findings may instead suggest a dynamical role in which the produc tion of C8HSL signal provides an energetic advantage by suppressing sensitivity of t he luminescence switch during the growth of a population.
Noise and crosstalk in Vibrio fischeri n Methods Fluorescence and luminescence response of bulk cult ure V.fischeri mutant JB10 is a derivative of the ES114 strain in which a chromosomal gfp reporter is inserted into the lux operon by allele exchange, producing luxI-gfp-luxCDABEG . We prepared JB10 from a glycerol stock and grew the cells to exponential phase in define d artificial seawater medium  to which was added 0.3% casamino acids. Cells were then dilute d and regrown to OD ~ 0.1-0.3 in fresh medium, washed three times, and then rediluted 100 into a 96-well assay plate containing fresh medium. The individual wells were preloaded w ith an 118 array of concentrations of the two HSL autoinducers N -3-oxohexanoylL -homoserine lactone (3OC6HSL, Sigma #K3007) and NoctanoylL -homoserine lactone (C8HSL, Cayman Chemical Co. #10011199). The well plate was then incubated in a Biotek Synergy 2 plate reader at 25 C, giving a growth rate 1.1 0.1 hr -1 Optical density was measured at 600 nm, and GFP fluorescence was measured using a 485/20 nm excitation filter and a 528/20 nm emissi on filter. The optical density, luminescence and GFP fluorescence values for each well we re recorded at regular intervals during exponential growth (Figure 2). Data collected ea rly in growth ( t < 12 hrs) showed a sensitive dependence on the exogenous levels of both HSLs, indicating that endogenous HSL did not accumulate significantly during this interval. Competitive inhibition model for bulk response In order to generate a mathematical representation of the lux response, as a function of the 3OC6HSL and C8HSL signals, we fit the JB10 well-plate data (fluorescence vs HSL concentrations) to the competitive inhibition model of Figure 1 [11, 13]. In this model lux is regulated primarily through competition between C8HSL and 3OC6HSL to form LuxR complexes that act as transcriptional activators for the lux genes. The action of C8HSL on LuxR synthesis through AinR and the phosphorelay is not considered. We assume that 3OC6HSL and C8HSL diffuse freely across the cell envelope and form multimeric complex es with LuxR. We allow an arbitrary degree of multimerization but we do not consider heterocomplexes ( i.e. involving both C8HSL and 3OC6HSL). Although it is simple to include the weak activation of lux by C8HSL-LuxR, which is evident in the bioluminescence data at low 3OC6HSL
Noise and crosstalk in Vibrio fischeri n concentrations, this activation is scarcely visible in the GFP fluorescence d ata that is the target of our modeling. Therefore we omitted this mechanism from our model and considered C8HSL only in its role as a competitor for LuxR. That is, we assume that the GFP fluor escence is proportional to the concentration of the 3OC6HSL-LuxR complex. More biochemical accurac y could be included by introducing extra parameters, but the simpler model appears suffi cient to describe the JB10 data. The model allows C8HSL and 3OC6HSL to form multimeric complexes (of degree m and n respectively) with LuxR, where the Hill coefficients n and m are not assumed to be integers. These equilibria are characterized by two dissociation constants, K 1 and K 2 : K 1 and K 2 are defined so as to have units of concentration, regardless of the values of m and n If [LuxR 0 ] is the average total concentration of LuxR, including complexes, then [LuxR 0 ] = [LuxR] + n [(3OC6HSL-LuxR) n ] +
+ m [(C8HSL-LuxR) m ] Eqn. ( 3) As we do not measure the actual LuxR copy number (although see ), it is convenient t o redefine the dissociation constants in terms of [LuxR 0 ] and a dimensionless concentration r : r = [LuxR]/[LuxR 0 ] [ ] [ ] () m m m m LuxR HSL C HSL C LuxR K = 8 8 12 1 [ ] [ ] () n n n n LuxR HSL OC HSL OC LuxR K = 6 3 6 3 12 2 ( ) mLuxR HSL C LuxR m HSL C m + 8 8 ( ) nLuxR HSL OC LuxR n HSL OC n + 6 3 6 3
Noise and crosstalk in Vibrio fischeri n Eqn. ( 4) Here k 1 and k 2 have dimensions of (autoinducer) concentration. Then Eqn. (3) becomes Eqn. (5) Starting from the HSL concentrations and an initial guess for the paramete rs ( k 1 k 2 m, n ) we solve Eqn. (5) to find r. Then Eqn. (4) gives the concentrations (relative to LuxR 0 ) of the two multimer species. We compare the model to the well-plate data by assuming t hat the GFP fluorescence F is a linear, non-saturating function of the two multimer concentrations : Eqn. (6) Here F 0 a 1 and a 2 are positive constants (see Results ). As explained above, a 1 is evident in luminescence but is scarcely detectable in the fluorescence; setting a 1 = 0 does not impair the fit. Then the shape of the 2D surface F (3OC6HSL, C8HSL) is determined solely by the four parameters k 1 k 2 n, and m, while the parameters F 0 and a 2 provide an instrument-dependent offset and amplitude that scale the 2d model F surface onto the measured values. We estimate the four model parameters through a nonlinear least squares fit of the fluorescence r esponse surfaces F (3OC6HSL, C8HSL) measured at optical densities 0.05 0.15 cm -1 to Eqn. (6), with the scale parameters a 2 and F 0 determined by linear regression. This provides a parametrization of the average response F as a function of the two HSL inputs (Table 1). [ ] () 0 1 8 8 LuxR LuxR HSL C HSL C r k m m m m = [ ] () 0 2 6 3 6 3 LuxR LuxR HSL OC HSL OC r k n n n n = [ ] n n n m m mk HSL OC nr k HSL C mr r2 1] 6 3[ 8 1 + + = ( ) [ ]  ()  0 2 0 1 0 6 3 8 LuxR LuxR HSL OC a LuxR LuxR HSL C a F F n m + + =
Noise and crosstalk in Vibrio fischeri The data do not require that m and n are different. As Table 1 indicates, the fit yield s similar values for the two Hill coefficients ( m = 1.1 0.4 and n = 1.35 0.05), and in fact we obtain a very similar fit if we assume that the same coeffic ient applies for both autoinducers ( m = n = 1.2 0.2). Table 1 also shows (as expected from Eqn. ( 1)) that we obtain similar parameters when we fit Eqn. (6) to the square root of the measured luminescence L1/2 rather than to the GFP fluorescence F. Microfluidic studies of individual cells To measure the effect of exogenous HSL signals on lux expression in individual JB10 cells we loaded cells into microfluidic perfusion chambers t hat supplied a flow of medium containing exogenous 3OC6HSL and C8HSL. Each microfluidic dev ice consisted of three parallel and unconnected channels (Figure 4), with each channel having width 400 m (parallel to the observation window but perpendicular to the fluid f low), depth 10-15 m (perpendicular to the observation window), and length (parallel to observ ation window and to fluid flow) 10 mm. The devices were fabricated from PDMS silicone elastomer (Sylgard 184, Dow-Corning Corporation) by a standard soft-lithographic method in which a PDMS replica is cast from a reactive ion-etched silicon master . The devic e channels were sealed by a glass coverslip bonded to the PDMS. In order to promote cell adhesi on to the interior of the glass window, we coated the interior of the device by filling it wit h a solution of polyL -lysine (1 mg/ml, MW 300 000) and incubating it for 24 hours at 5 C, prior to cell injection. This provided stable a dhesion of the V.fischeri to the glass window. JB10 cells for microfluidic studies were prepared i n exponential phase as for the 96-well assay above: We grew cells to exponential phase in define d artificial seawater medium with casamino acids , then washed (3) and rediluted the cell s, and then regrew them to OD (600 nm) = 0.015 0.03 cm-1 in fresh medium. Once the cells and the microflui dic device were prepared, we flushed the polyL -lysine solution by pumping the JB10 culture into a ll three parallel channels at 1-2 ml/hr with a syringe pump. We then placed the device (with glass window facing downward) on the stage of a Nikon TE2000U microscop e and reduced the flow rate to ~0.02 ml/hr. At this slow flowrate the cells gradually s ettle and adhere to the glass window. Once a sufficient number of cells had adhered to the windo w (requiring 15-30 minutes), we supplied
Noise and crosstalk in Vibrio fischeri n autoinducer by connecting the device inputs to syri nge pumps that delivered defined medium containing exogenous 3OC6 HSL and/or C8HSL. Each o f the three channels was supplied with a different combination of HSLs, flowing at a rate ~0.02 ml/hr during fluorescence measurements. The 0.02 ml/hr flow rate of medium corresponds to a n average flow velocity of ~1 mm/s within each channel. Both the device design and experimen tal testing ensured that this flow was sufficiently uniform and rapid to wash away endogen ous (natively produced) autoinducer that might otherwise affect activation of the lux genes. First of all, control experiments in our fl ow system showed that in the absence of any exogenou s autoinducer (HSL) gfp expression from the lux reporter strain was at its baseline level (and lumi nescence was unobservable). Moreover, the physical parameters of the flow system make it highly implausible that spatial heterogeneity in the flow could develop or allow experimentally r elevant concentrations of HSL to accumulate near any of the cells under observation: First, th e dimensions of the device and the flow rate of growth medium lead to fluid flow at a very low Reyn olds number ( Re ~ 0.03). At this Re the flow velocity profile is highly uniform across the width and length of the flow chamber, up to within ~10 m m of the chamber edges . Second, the only sign ificant heterogeneity in this flow velocity profile occurs along the depth of the channel ( i.e perpendicular to the window), which is 10-15 m m. However HSL requires only ~0.1 s to diffuse thi s distance. This is so much faster than other relevant time scales in the exper iment that a meaningful HSL gradient cannot be established in this direction. Third, the chamber v olume and the 1 mm/s flow rate together indicate that the entire volume of the cell chamber region (10 mm length) is completely flushed every ~10 seconds. However the only cells occupyin g the chamber (and producing HSL) are those forming a sparse single layer (cells are typi cally spaced > 20 m m apart) on the chamber window. Literature estimates of HSL production rat es in V.fischeri indicate that such a sparse layer of individual cells, within a chamber that is flushed at this rate, would not be able to generate an endogenous HSL concentration above ~100 pM . This concentration is at least two orders of magnitude smaller than the exogenous HSL concentrations that we are providing.
Noise and crosstalk in Vibrio fischeri Finally, if the cells did generate enough HSL to af fect local concentrations, we would expect that cells downstream would in general express more GFP than cells upstream. More generally we would expect the correlation Cij between the GFP fluorescence Fi, Fj of a pair of cells i, j Cij = ( Fim )( Fj m )/ s2 Eqn. (7) to depend on their spatial separations xij, yij, or rij. (Here m is the mean cell fluorescence and s2 is the variance in F .) We analyzed our data for such spatial correlati ons and found none. For example, Figure 4 shows no relationship between Cij and rij: the gfp expression of two neighboring cells is no more similar than that of t wo distant cells. In short the data and the system design argue strongly against any autoactiva tion of (or local crosstalk between) the individual cells under observation. Characterizing heterogeneity in lux activation The three-channel device allowed us to collect the fluorescence histogram of cells under three different HSL signal combinations, as it evolved ov er 4-5 hours. Once HSLs were introduced to the device at t= 0 we collected phase contrast and fluorescence image pairs for each channel (HSL combination) at intervals of 20 minutes, using a 20/0.50NA phase objective and a GFP filter cube. Images were recorded by a Coolsnap HQ 2 camera (Photometrics) at 30 C and corrected in software for dark current and flat-fie ld. For each experimental condition we evaluated the fl uorescent emission from (typically) ~200 individual cells by first determining the physical locations (pixel coordinates) of single cells in a phase contrast image. We then used a homemade Matla b code to evaluate the fluorescence per cell pixel in the associated fluorescence image by summing the fluorescence emission (relative to background) of the contiguous bright pixels associa ted with the cells pixel coordinates. Normalizing the histogram of individual cell fluore scence values gives a distribution P ( F | ([3OC6HSL],[C8HSL])), representing the probabilit y of cell fluorescence F given the HSL input concentrations. Calculating the mutual information By combining mathematical parametrizations of both the lux response, F ([3OC6HSL],[C8HSL]) and the probability distribution P ( F | ([3OC6HSL],[C8HSL])) we calculate the mutual
Noise and crosstalk in Vibrio fischeri information  between the HSL signal inputs and lux output. The signal concentrations are inconvenient parameters for this calculation becaus e only an infinite concentration of autoinducer can saturate the response. In their an alysis of V.harveyi QS, Mehta et al.  defined new coordinates that describe the state of saturation of the autoinducer receptors. In similar fashion we replace [3OC6HSL] and [C8HSL] wi th coordinates X and Y that describe the state of the association equilibria for the HSL-Lux R complexes: Eqn. (8) X =1 or Y = 1 corresponds to complete saturation of the 3OC6H SL-LuxR or C8HSL-LuxR binding equilibrium respectively. With these coord inates the response surface F ( X,Y ) for the competitive inhibition model has a simple shape (Fi gure 7) that is independent of the parameters k1, n k2, m. The mutual information between a combination of inp uts ( X,Y ) and the output F is then calculated as: Eqn. (9) Here P ( F ) is the probability of finding output F, in the absence of any knowledge of the input ( X,Y ). P ( F |( X,Y )) is the probability of F, given the combination ( X,Y ). P ( F, ( X,Y )) is the probability of observing the particular combination F, ( X,Y ): These probability distributions are normalized as f ollows: [ ]  n n n k HSL C HSL OC X 2 6 30 6 3 + = [ ]  m m m k HSL C HSL C Y 1 8 8 + =() ()() () ( ) ( ) ()() () ()() () () () = = F P Y X F P Y X P Y X F P dFdXdY Y X P F P Y X F P Y X F P dFdXdY Y X F I | log , | , log , ; 2 2() ()() () () () = = = = = 1 | 1 , | )), (,( ) ( Y X F P dF F P dF Y X P Y X F P dXdY Y X F P dXdY F P ( ) ( ) ( ) ( ) ( ) Y X P Y X F P Y X F P , , =
Noise and crosstalk in Vibrio fischeri To evaluate Eqn. (9) we model P ( F ,( X,Y )) as the gamma distribution that has the same mean and variance as observed in the bulk and single-cell me asurements respectively. The calculation also requires an estimate of P ( X,Y ), the prior probability of a particular combinatio n ( X,Y ). P ( X,Y ) is not so easily predicted. However, given that X and Y are both bounded by 0 and 1 we made the straightforward assumption that P ( X,Y ) = constant. The mutual information Eqn. (9) is then found to be I 0.53 bits. However this result is not sensitive t o our assumptions about the prior probability: various P ( X,Y ) functions that were strongly bimodal in both X and Y, and either symmetric or asymmetric in X vs. Y , all gave similar values of I 0.5 bits. List of Abbreviations C8HSL C8 homoserine lactone, NoctanoylL -homoserine lactone; GFP green fluorescent protein; HSL homoserine lactone; JB10 Vibrio fischeri strain JB10; QS quorum sensing; 3OC6HSL 3-oxo-C6 homoserine lactone, N -3-oxohexanoyl homoserine lactone; Authors Contributions PDP Acquired and analyzed data; JTW Acquired da ta; SJH Designed research, analyzed data, and drafted the manuscript. All authors read and approved the final manuscript. Acknowledgment Prof. Eric Stabb and Prof. Rahul Kulkarni provided valuable advice and many suggestions during this work. Prof. Eric Stabb also provided t he strain JB10 used in this study. The silicon master for the microfluidic devices was designed an d fabricated by Minjun Son. The authors gratefully acknowledge funding support from the Nat ional Science Foundation under award DMR-0851707 and the National Institutes of Health N IDCR under award 1R21DE018826. Publication of this article was funded in part by t he University of Florida Open-Access () = = 1 1 )), (,( Y X P dXdY Y X F P dFdXdY
Noise and crosstalk in Vibrio fischeri Publishing Fund. The funding agencies played no ro le in the research design and activity or in the preparation and submission of this manuscript.
Noise and crosstalk in Vibrio fischeri rrrr nrr nnrnnn r !"#r$%"r%& nnnrnnr r '(r)$r*r) % n +rr%',(##r$(rr#-"r+" #."',r# nnn nn n nrnnn /"#(rr)r012 *(',3#"45 nnn#$%r rrnr '"%""$-r. &'( n /1!&%"r%& "rrnnn)%r# 'r)6rr# *& n 62"-r)'% nn#rrnn nn)% !r2 7#r8"2r)r032#r.%12"r r n r1-r4' +r rrnn %r)nn !! nn n 59,& ,r nrnnr"r-rrn %"#r '( nn 3#(r35"r(246 .n//nn0r#1#. r/)# %"#r 2 nn n&(88$(246 n )nn #rn r rn rn3rnr %"#r '( nn nn&(88$:"096rrr14;(246 4# n nn r rn nr nnnn$nrn rrnnn# # nn n n&(88$(246 n ,$/rn.n/rr5 )nr n #r %"#r '6 n n n,('$""-"3*("8;5 +rnnnn$n 7 rn #nrnnrn/rn1# %"#rn (' n n59,&<#r*<" /1"( 246 n 8n9r rnrn0 nn3rnrn5n#-,# rn %"#r n
Noise and crosstalk in Vibrio fischeri n,('%(1-/5*("8;5 +7. ,-n,rn.n ,nnn#+rn/ "r n %"#rn (6 n3-"rr'&+"=r9"%&4r-"'6rrr1 4; nnn n 9n rn,$0nrn)#rnrnr %"#rn (' n n n$$+'<+&r"r##r 5n1#,$0n r n5 rrr #n rn !! n n"r!rr133#"4559,& (246 .nr nn rrr#rn /"#(r *!' n : n%r&!rr13$#*"r3 r8#r'/*(',3#"45 nn 0r:5n#$0$nn0n #5"nn1 #/rnrn:n5nn '884)nn (( ;rr=;*+"1r3 -n0nr/nn/nr n#,nn n5r n ;&3>rn nnrn nr1+"?(r#$8r*;r$6rr ." :$nrrnrnr5rrn nrnnn nnrr %.32#r.%12n &2r/&*(',%r&3#"45 "#+rn1#:/*05rn nrn 0r4#rn"nnrn:$ n#$n %"#rn 2 nn nn $)r-."' nnr r+9r/!r2 7#rrr@ 9"96$""$6%"r9! nrn)rn3rnnrnnr rn ;/"#'"3:3' n n r-#";62"3&1%"r%&!1rr /3 nrnnrnnrnrnn rrnn r("32#r.%12 <(r*$-r!'r6."'&r)r+ "88r!&.!< :$nrrnnnr nrnn%r #r$ ;rr1#-r/"#"'"r.23rrn ( n $"&
Noise and crosstalk in Vibrio fischeri ".r3/;".r6,Crrr13<"1r #. ,nnnn#rn r n .nrnr # nr rn % 26 nn nn/;>'<"1r<"1r#. r3,Crrr133#rr1; ;n# .-,r nr )rn rrnr<< %n &2 n n$D$*;r#r6*'r3&""#r $(.'(#*B"&3"2r633.8& .nr 7n/nn >7$3'("7#r1"#)r%12 ( n E-(3$-r&/'-"", +rnrn###r n#n 3rnn $rrrn ;-2r)&r## 2! nnn (rr,(##r$+rr%'#-"r+ "#."' nrn nnrnnnn "#-% !& ,(##r$+rr%' nr )nnrn n -r % n n !".!$(A&r)-r9'3#r)r' 0rnn5nrnn )n%r% r32#% 'r1-3;r4$-5$"#3r <)"#'3?("#' nnr5 rnnnnnrnrr nrn <43&r## 2 n nn &1(,$!"1Br-#";>1/;%"r %&!1rr/3 rn n#nrn /nn/nr)#/n0n ;&3% ( 'r#=r1r$;-(1, -n nnnrn n n (& 9,$!1rr/3%"r%& 1#r5 nnrnn ) rnnn$n 4%> n %r&,.:%"#909!6("(;>)r r2'&2r/&59,&3#"45 "nnn n n #$n5r. 6! n(246/r",+ / .rn nn )##+rn,n =# rrr .+/ "rn"rr/ %%(n r13!"1B(,$&1r-#";!1rr /3%"r%>1/; +rn# 7#+r/nn0rr "rr %8-2n 2' n
Noise and crosstalk in Vibrio fischeri 3"3,!-#rr6 +5rrn ># $rn?r 4r#8-r n $#"#r(5/ rrn /r0B9381r 5r"1@n
Noise and crosstalk in Vibrio fischeri Figure Legends Figure 1. Schematic of QS regulation of V.fischeri bioluminescence and competitive model. ( A ) QS regulation of bioluminescence in V.fischeri uses three autoinducer channels . The autoinducer 3OC6HSL is synthesized by LuxI and bind s to LuxR to form a transcriptional activator for the bioluminescence genes luxCDABEG. Two more autoinducers (C8HSL and AI2) drive the phosphorelay that regulates producti on of LuxR as well as other colonization behaviors. ( B ) A simplified model considers only competitive int eraction between 3OC6HSL and C8HSL, as proposed by Kuo et al  and Lupp et al . The receptor LuxR binds the autoinducers C8HSL and 3OC6HSL to form multimeric c omplexes (of degree m and n respectively) which activate lux transcription. In our fit we omit lux activation by the C8HSL complex of LuxR. Figure 2. Data and fit for lux activation of JB10 strain by two autoinducers The figure shows the population-average (bulk measu rement) response of combinations of exogenous C8HSL and 3OC6HSL autoinducers, as measur ed by GFP fluorescence of strain JB10. ( A ) Data at OD = 0.15 cm -1 ( B ) fit to competitive inhibition model, ( C ) residual from fit, and ( D ) Data and fit sectioned at constant [3OC6HSL]. Figure 3. Correlation between luminescence and GFP fluorescence of JB10 strain ( A ) Correlation between bioluminescence L and GFP fluorescence F of the JB10 strain, as observed in a bulk (well-plate) experiment. Data w ere collected at OD = 0.10 in the presence of exogenous C8HSL and 3OC6HSL spanning a range of con centrations from zero to ~ 1000 nM. ( B ) Same data plotted as L 1/2 vs F showing the approximate power law of Eqn. (1). ( C ) The HSL-independent intercept F 0 in Eqn. (1) grows in proportion to optical density during early growth. Figure 4. Microfluidic measurements of GFP fluoresc ence ( A ) Schematic of microfluidic device for measuring lux expression in individual V.fischeri Growth medium containing exogenous HSL flows throug h three parallel rectangular channels that are cast into the lower surface of a PDMS bloc k. Live cells in channels adhere to the glass
Noise and crosstalk in Vibrio fischeri n window (coverslip) that seals the channels from ben eath, and are observed in an inverted microscope. The shallow ratio of channel height to width ( h/w ~ 0.02-0.04) ensures a uniform flow velocity profile across the width and length o f each channel. ( B ) Contour map of lux activation F versus HSL input. The white circles show HSL combi nations applied to the cells during the microfluidic experiments described in th e text. The contour labels show the activation fraction above the base level, i.e. ( F-F 0 )/max( F-F 0 ), as derived from the bulk measurements and competitive inhibition model. ( C ) Histogram comparing the correlation C ij (Eqn. (7)) in gfp expression of a pair of cells ( i,j ) to the physical separation r ij = ( x ij 2 + y ij 2 ) between those cells. The color of each bin indicates the number of cell pairs ( i,j ) whose physical separation and brightness correlation fall within that bin. Pairs of near-neighbor cells are not more correlated in their lux activation than pairs of distant cells. Figure 5. Temporal response of the lux activation histogram ( A )-( E ) Progression of the GFP expression histogram for i ndividual cells in the microfluidic chamber, as measured for the three signal combinati ons indicated as Experiment #1 ( i.e. 100 nM 3OC6HSL and variable C8HSL) in Figure 4 B ; ( F )-( J ) Percentiles in gfp fluorescence of individual cells under different HSL conditions. E ach curve shows the cumulative area under the corresponding histogram in panels ( A )-( E ) (and with the same color scheme). In the absence of C8HSL the distribution broadens (and increases in m ean value) over 3-4 hrs, while the presence of C8HSL suppresses this response. Figure 6. Cell-to-cell heterogeneity at fixed avera ge activation Cell fluorescence brightness distributions measured for experiments #2 and #3 in Figure 4 B ( A )-( C ) Distributions collected for two-HSL combinations that generated ~25% of full lux activation (Experiment #2) and ( D )-( F ) distributions collected for combinations that gen erated ~60% of full lux activation (Experiment #3). The dashed curves show maximum likelihood fits to a gamma distribution. Each histogram is derived from roughly 200 individual cells. Figure 7. Response surface F ( X,Y ) and the coefficient of variation ( A ) The (population-average) response surface F ( X,Y ) generated by the competitive model with the change of variable, Eqn. (8). ( B ) The measured coefficient of variation CV = s / m in lux
Noise and crosstalk in Vibrio fischeri activation of individual cells declines at higher a ctivation levels. Relative activation is defined as ( F-min ( F ))/( max ( F )-min ( F )). The dashed line shows a linear fit, used to pa rameterize CV in the calculation of the mutual information I ( F; ( X,Y )) Table 1. Parameter values for the competitive inhi bition model Parameter values for the competitive model were obt ained by fitting Eqn. (6) to the GFP expression data (Figure 2), for a bulk culture of Vibrio fischeri strain JB10 at density OD = 0.05 0.15 cm -1 The first row contains parameters obtained by fi tting the GFP fluorescence data ( F ) with the four-parameter model described in the text The second row contains parameters obtained by fitting the fluorescence data under the constraint that m = n. The third row contains fit parameters obtained when the four-parameter mod el is fit to the luminescence data (actually L 1/2 ) rather than to the fluorescence F. (See Eqn. (1).) Fit k1 (nM) m k2 (nM) n 1 39 8 1.1 0.4 163 15 1.35 0.05 2 46 2 1.2 0.2 179 6 (= m ) 3 38 5 0.8 0.2 118 11 1.6 0.06
! # Figure 2
012 0 Fluorescence (10 5 cts) 012 000.050.10.15 OD 01234 F0(104cts) Fluorescence (10 5 cts) 200 400Luminescence 1/2 0.5 1 1.5 2Luminescence (105cts) OD = 0.10 OD = 0.10 F 0 Figure 3
0 0.5 1 0 0.5 0 0.5 0 0.5 0501000 0.5 lux activation (GFP)0 nM 500 nM 1000 nM C8HSL t = 0.7 hrst = 2.1 hrst = 3.1 hrs t = 4.1 hrst = 4.5 hrs # $lux activation (GFP)100 25102050 t = 0.7 hrs t = 2.1 hrs t = 3.1 hrs t = 4.1 hrs 00000 5050505050 100100100100100 t = 4.5 hrs% & ( )Histograms in lux activation Percentiles in lux activation Figure 5
GFP fluorescence 0 50 100 150 0 0.05 0.1 0 nM C8HSL200 nM 3OC6HSL 0 50 100 150 0 0.05 0.1 20 nM C8HSL500 nM 3OC6HSL 0 50 100 150 0 0.05 0.1 100 nM C8HSL800 nM 3OC6HSL 0 20 40 60 0 0.05 0.1 300 nM C8HSL500 nM 3OC6HSL 0 20 40 60 0 0.05 0.1 1000 nM C8HSL1000 nM 3OC6HSL 0 20 40 60 0 0.05 0.1 0 nM C8HSL50 nM 3OC6HSLGFP fluorescence Figure 6
0 0.5 1 0 0.5 1 0 1 2 3 10 5 Y X Fluorescence (model) 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 Relative ActivationCoefficient of variation Figure 7