elife01916.xml

<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.1d1 20130915//EN"  "JATS-archivearticle1.dtd"><article article-type="research-article" dtd-version="1.1d1" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"><front><journal-meta><journal-id journal-id-type="nlm-ta">elife</journal-id><journal-id journal-id-type="hwp">eLife</journal-id><journal-id journal-id-type="publisher-id">eLife</journal-id><journal-title-group><journal-title>eLife</journal-title></journal-title-group><issn publication-format="electronic">2050-084X</issn><publisher><publisher-name>eLife Sciences Publications, Ltd</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="publisher-id">01916</article-id><article-id pub-id-type="doi">10.7554/eLife.01916</article-id><article-categories><subj-group subj-group-type="display-channel"><subject>Research article</subject></subj-group><subj-group subj-group-type="heading"><subject>Biophysics and structural biology</subject></subj-group><subj-group subj-group-type="heading"><subject>Microbiology and infectious disease</subject></subj-group></article-categories><title-group><article-title><italic>Escherichia coli</italic> swimming is robust against variations in flagellar number</article-title></title-group><contrib-group><contrib contrib-type="author" id="author-9083"><name><surname>Mears</surname><given-names>Patrick J</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="aff" rid="aff2"/><xref ref-type="fn" rid="con1"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-9236"><name><surname>Koirala</surname><given-names>Santosh</given-names></name><xref ref-type="aff" rid="aff3"/><xref ref-type="fn" rid="con2"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-9237"><name><surname>Rao</surname><given-names>Chris V</given-names></name><xref ref-type="aff" rid="aff3"/><xref ref-type="other" rid="par-7"/><xref ref-type="fn" rid="con3"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" corresp="yes" id="author-9239"><name><surname>Golding</surname><given-names>Ido</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="aff" rid="aff2"/><xref ref-type="aff" rid="aff4"/><xref ref-type="corresp" rid="cor1">*</xref><xref ref-type="other" rid="par-1"/><xref ref-type="other" rid="par-4"/><xref ref-type="other" rid="par-5"/><xref ref-type="other" rid="par-6"/><xref ref-type="fn" rid="con4"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" corresp="yes" id="author-2674"><name><surname>Chemla</surname><given-names>Yann R</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="aff" rid="aff2"/><xref ref-type="corresp" rid="cor2">*</xref><xref ref-type="other" rid="par-1"/><xref ref-type="other" rid="par-2"/><xref ref-type="other" rid="par-3"/><xref ref-type="fn" rid="con5"/><xref ref-type="fn" rid="conf1"/></contrib><aff id="aff1"><institution content-type="dept">Department of Physics</institution>, <institution>University of Illinois at Urbana-Champaign</institution>, <addr-line><named-content content-type="city">Urbana</named-content></addr-line>, <country>United States</country></aff><aff id="aff2"><institution content-type="dept">Center for the Physics of Living Cells</institution>, <institution>University of Illinois at Urbana-Champaign</institution>, <addr-line><named-content content-type="city">Urbana</named-content></addr-line>, <country>United States</country></aff><aff id="aff3"><institution content-type="dept">Department of Chemical and Biomolecular Engineering</institution>, <institution>University of Illinois at Urbana-Champaign</institution>, <addr-line><named-content content-type="city">Urbana</named-content></addr-line>, <country>United States</country></aff><aff id="aff4"><institution content-type="dept">Verna and Mars McLean Department of Biochemistry and Molecular Biology</institution>, <institution>Baylor College of Medicine</institution>, <addr-line><named-content content-type="city">Houston</named-content></addr-line>, <country>United States</country></aff></contrib-group><contrib-group content-type="section"><contrib contrib-type="editor"><name><surname>Laub</surname><given-names>Michael</given-names></name><role>Reviewing editor</role><aff><institution>Massachusetts Institute of Technology</institution>, <country>United States</country></aff></contrib></contrib-group><author-notes><corresp id="cor1"><label>*</label>For correspondence: <email>igolding@illinois.edu</email> (IG);</corresp><corresp id="cor2"><label>*</label>For correspondence: <email>ychemla@illinois.edu</email> (YRC)</corresp></author-notes><pub-date date-type="pub" publication-format="electronic"><day>11</day><month>02</month><year>2014</year></pub-date><pub-date pub-type="collection"><year>2014</year></pub-date><volume>3</volume><elocation-id>e01916</elocation-id><history><date date-type="received"><day>19</day><month>11</month><year>2013</year></date><date date-type="accepted"><day>29</day><month>12</month><year>2013</year></date></history><permissions><copyright-statement>© 2013, Mears et al</copyright-statement><copyright-year>2013</copyright-year><copyright-holder>Mears et al</copyright-holder><license xlink:href="http://creativecommons.org/licenses/by/3.0/"><license-p>This article is distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/3.0/">Creative Commons Attribution License</ext-link>, which permits unrestricted use and redistribution provided that the original author and source are credited.</license-p></license></permissions><self-uri content-type="pdf" xlink:href="elife01916.pdf"/><abstract><object-id pub-id-type="doi">10.7554/eLife.01916.001</object-id><p>Bacterial chemotaxis is a paradigm for how environmental signals modulate cellular behavior. Although the network underlying this process has been studied extensively, we do not yet have an end-to-end understanding of chemotaxis. Specifically, how the rotational states of a cell’s flagella cooperatively determine whether the cell ‘runs’ or ‘tumbles’ remains poorly characterized. Here, we measure the swimming behavior of individual <italic>E. coli</italic> cells while simultaneously detecting the rotational states of each flagellum. We find that a simple mathematical expression relates the cell’s run/tumble bias to the number and average rotational state of its flagella. However, due to inter-flagellar correlations, an ‘effective number’ of flagella—smaller than the actual number—enters into this relation. Data from a chemotaxis mutant and stochastic modeling suggest that fluctuations of the regulator CheY-P are the source of flagellar correlations. A consequence of inter-flagellar correlations is that run/tumble behavior is only weakly dependent on number of flagella.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.001">http://dx.doi.org/10.7554/eLife.01916.001</ext-link></p></abstract><abstract abstract-type="executive-summary"><object-id pub-id-type="doi">10.7554/eLife.01916.002</object-id><title>eLife digest</title><p><italic>Escherichia coli</italic> is a rod-shaped bacterium commonly found in the lower intestines of humans and other warm-blooded animals. While most strains of <italic>E</italic>. <italic>coli</italic> are harmless, including most of those found in the human gut, some can cause diseases such as food poisoning. Due to its close association with humans and the fact that it is easy to grow and work with in the laboratory, <italic>E</italic>. <italic>coli</italic> has been intensively studied for over 60 years.</p><p>Many bacteria are capable of ‘swimming’ by using one or more flagella. These rotating whip-like structures are each driven by a reversible motor, and they act a bit like a propeller on a boat. While some bacteria have only a single flagellum, others, such as <italic>E</italic>. <italic>coli</italic>, have multiple flagella distributed over the cell surface. Rotating all their flagella in a counterclockwise direction allows the bacterium to swim—and it has been proposed that the clockwise movement of at least one flagellum will cause the bacterium cell to stop swimming and start tumbling.</p><p><italic>E. coli</italic> is able to control the time it spends swimming or tumbling to move towards a nutrient, such as glucose, or away from certain harmful chemicals. However, the details of how the number of flagella and the direction of rotation (clockwise or counterclockwise) influence the motion of the bacterium are not fully understood.</p><p>Now, Mears et al. have used ‘optical tweezers’ to immobilize individual <italic>E. coli</italic> cells under a microscope, and then track both their swimming behavior and the movements of their flagella. This revealed that the individual flagella on the same cell tend to move in a coordinated way. Therefore, whilst tumbling could be caused by a single flagellum stopping swimming behavior, it often involved a concerted effort by many of the cell’s flagella.</p><p>After observing that <italic>E. coli</italic> cells with more flagella spent less time tumbling than would be predicted if a single flagella always ‘vetoed’ swimming, Mears et al. propose a new mathematical relationship between the number of flagella on the cell, the direction of rotation, and the resulting probability that the cell will tumble. This work shows that swimming behavior in bacteria is less affected by variations in the number of flagella than previously thought—and this phenomenon may provide evolutionary advantages to <italic>E. coli</italic>. The next step is to explore the mechanism by which bacteria coordinate their flagella.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.002">http://dx.doi.org/10.7554/eLife.01916.002</ext-link></p></abstract><kwd-group kwd-group-type="author-keywords"><title>Author keywords</title><kwd>bacterial chemotaxis</kwd><kwd>optical tweezers</kwd><kwd>single-cell studies</kwd><kwd>flagella</kwd></kwd-group><kwd-group kwd-group-type="research-organism"><title>Research organism</title><kwd><italic>E. coli</italic></kwd></kwd-group><funding-group><award-group id="par-1"><funding-source><institution-wrap><institution>National Science Foundation</institution></institution-wrap></funding-source><award-id>PHY-082265</award-id><principal-award-recipient><name><surname>Golding</surname><given-names>Ido</given-names></name><name><surname>Chemla</surname><given-names>Yann R</given-names></name></principal-award-recipient></award-group><award-group id="par-2"><funding-source><institution-wrap><institution>Burroughs Wellcome Fund</institution></institution-wrap></funding-source><principal-award-recipient><name><surname>Chemla</surname><given-names>Yann R</given-names></name></principal-award-recipient></award-group><award-group id="par-3"><funding-source><institution-wrap><institution>Alfred P. Sloan Foundation</institution></institution-wrap></funding-source><principal-award-recipient><name><surname>Chemla</surname><given-names>Yann R</given-names></name></principal-award-recipient></award-group><award-group id="par-4"><funding-source><institution-wrap><institution>National Institutes of Health</institution></institution-wrap></funding-source><award-id>R01 GM082837</award-id><principal-award-recipient><name><surname>Golding</surname><given-names>Ido</given-names></name></principal-award-recipient></award-group><award-group id="par-5"><funding-source><institution-wrap><institution>Welch Foundation</institution></institution-wrap></funding-source><award-id>Q-1759</award-id><principal-award-recipient><name><surname>Golding</surname><given-names>Ido</given-names></name></principal-award-recipient></award-group><award-group id="par-6"><funding-source><institution-wrap><institution>National Science Foundation</institution></institution-wrap></funding-source><award-id>PHY-1147498</award-id><principal-award-recipient><name><surname>Golding</surname><given-names>Ido</given-names></name></principal-award-recipient></award-group><award-group id="par-7"><funding-source><institution-wrap><institution>National Institutes of Health</institution></institution-wrap></funding-source><award-id>R01 GM054365</award-id><principal-award-recipient><name><surname>Rao</surname><given-names>Chris V</given-names></name></principal-award-recipient></award-group><funding-statement>The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.</funding-statement></funding-group><custom-meta-group><custom-meta><meta-name>elife-xml-version</meta-name><meta-value>2</meta-value></custom-meta><custom-meta specific-use="meta-only"><meta-name>Author impact statement</meta-name><meta-value>The swimming behavior of <italic>E. coli</italic> is surprisingly robust against the natural variations in the number of flagella per cell.</meta-value></custom-meta></custom-meta-group></article-meta></front><body><sec id="s1" sec-type="intro"><title>Introduction</title><p>Many species of bacteria swim by rotating helical filaments called flagella (<xref ref-type="bibr" rid="bib3">Berg, 2004</xref>). A typical <italic>Escherichia coli</italic> cell is propelled by a bundle composed of multiple flagella. Each flagellum is controlled by a rotary motor that can switch between clockwise (CW) and counter-clockwise (CCW) rotation. When flagella on a cell rotate CCW, the cell swims along an approximately straight path called a ‘run’. When some of the flagella rotate CW, the bundle is disrupted causing an abrupt change in direction called a ‘tumble’ (<xref ref-type="bibr" rid="bib25">Macnab and Ornston, 1977</xref>). <italic>E. coli</italic> modulates the probability of being in one of these two swimming states in response to its environment, allowing it to navigate chemical, temperature, and light gradients (<xref ref-type="bibr" rid="bib4">Berg and Brown, 1972</xref>; <xref ref-type="bibr" rid="bib3">Berg, 2004</xref>). At any point in time, the probability that a flagellar motor rotates CW is determined by the concentration of phosphorylated signaling protein CheY (CheY-P). Coupling CheY phosphorylation to chemicals from the environment allows the cell to bias its random walk and migrate towards more favorable conditions. This biased random walk is called chemotaxis, and serves as a model for understanding how living organisms process information (<xref ref-type="bibr" rid="bib4">Berg and Brown, 1972</xref>; <xref ref-type="bibr" rid="bib42">Wadhams and Armitage, 2004</xref>; <xref ref-type="bibr" rid="bib32">Shimizu et al., 2010</xref>).</p><p>Tremendous progress has been made towards elucidating the mechanism of bacterial chemotaxis. The relationship between the chemotaxis signaling network and the CCW/CW rotational bias of the individual flagellar motor is now well mapped ([<xref ref-type="bibr" rid="bib6">Block et al., 1982</xref>; <xref ref-type="bibr" rid="bib10">Cluzel et al., 2000</xref>; <xref ref-type="bibr" rid="bib34">Sourjik and Berg, 2002</xref>; <xref ref-type="bibr" rid="bib44">Yuan et al., 2012</xref>]; for a review see <xref ref-type="bibr" rid="bib3">Berg, 2004</xref>), and has also been described using detailed mathematical models (<xref ref-type="bibr" rid="bib15">Emonet et al., 2005</xref>; <xref ref-type="bibr" rid="bib20">Jiang et al., 2010</xref>; <xref ref-type="bibr" rid="bib32">Shimizu et al., 2010</xref>). Despite this wealth of knowledge, how the CCW/CW states of individual motors collectively determine the run/tumble swimming behavior of the whole, multi-flagellated cell remains poorly understood. The number of flagella on an individual swimming cell can vary greatly, from one to more than ten (<xref ref-type="bibr" rid="bib11">Cohen-Ben-Lulu et al., 2008</xref>) (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1</xref>), yet very few studies are available to indicate how flagellar number affects swimming behavior. The only direct measurements of flagellar dynamics in swimming cells have been limited to short durations (∼1 s) (<xref ref-type="bibr" rid="bib38">Turner et al., 2000</xref>; <xref ref-type="bibr" rid="bib13">Darnton et al., 2007</xref>). The absence of long-term observations has precluded the development of a detailed mapping between flagellar state and cell swimming behavior. As a result, most theoretical models of bacterial chemotaxis are limited to treating an individual motor, or simply assume that all cells have a single flagellum (<xref ref-type="bibr" rid="bib8">Bray et al., 2007</xref>; <xref ref-type="bibr" rid="bib21">Kalinin et al., 2009</xref>; <xref ref-type="bibr" rid="bib26">Matthaus et al., 2009</xref>; <xref ref-type="bibr" rid="bib20">Jiang et al., 2010</xref>; <xref ref-type="bibr" rid="bib16">Flores et al., 2012</xref>). Quantifying the mapping from single-flagellum state to whole-cell swimming behavior thus remains a missing link to developing an end-to-end picture of bacterial chemotaxis.</p><p>A number of theoretical models have been put forward in an attempt to describe this mapping. One such model invokes a ‘voting’ mechanism, in which cells tumble only if a majority of flagella rotate CW (<xref ref-type="bibr" rid="bib19">Ishihara et al., 1983</xref>; <xref ref-type="bibr" rid="bib35">Spiro et al., 1997</xref>; <xref ref-type="bibr" rid="bib2">Andrews et al., 2006</xref>; <xref ref-type="bibr" rid="bib41">Vladimirov et al., 2008</xref>; <xref ref-type="bibr" rid="bib20">Jiang et al., 2010</xref>). However, by observing fluorescently labeled flagella during individual tumbles, Turner et al. established that CW rotation of a single flagellum is sufficient to ‘veto’ a run (<xref ref-type="bibr" rid="bib38">Turner et al., 2000</xref>) (<xref ref-type="fig" rid="fig1">Figure 1A</xref>). Refined versions of this ‘veto model’ were recently developed (<xref ref-type="bibr" rid="bib40">Vladimirov et al., 2010</xref>; <xref ref-type="bibr" rid="bib33">Sneddon et al., 2012</xref>), based on careful, slow-motion observations of tumbles (<xref ref-type="bibr" rid="bib13">Darnton et al., 2007</xref>). However, the extent to which these details are relevant for modeling swimming behavior is unknown, because no measurements have directly correlated long-term swimming behavior with flagellar activity in the same cell.<fig-group><fig id="fig1" position="float"><object-id pub-id-type="doi">10.7554/eLife.01916.003</object-id><label>Figure 1.</label><caption><title>Wild-type <italic>E. coli</italic> cells deviate from the ‘veto’ model.</title><p>(<bold>A</bold>) The mapping relating the run/tumble state of the cell to the CCW/CW state of its flagella according to the veto model. Schematic time trace from a cell with 3 flagella, showing CW (purple) and CCW (white) intervals for each flagellar motor and the resulting tumbles (black) and runs (white) of the cell. The veto model corresponds to an AND gate, by which cell runs only occur when all flagella rotate CCW (where CCW = 1, CW = 0, run = 1, and tumble = 0). (<bold>B</bold>) Schematic of a cell held by two optical traps (red cones) in the fluorescence excitation volume (green) within the sample chamber. (<bold>C</bold>) Representative data trace from a trapped cell with three flagella. Still images of fluorescently labeled flagella at 0.5-s intervals (top). The approximate location of the unlabeled cell body is indicated by a dashed yellow line. Flagella rotating CW (purple) and CCW (white) are numbered in frames in which they appear distinct. Corresponding cell-body rotation signal for the same cell (red line, bottom) as detected from deflections of the trapping laser. Tumbles (shaded area) were determined from the erratic cell-body rotation signal. (<bold>D</bold>) Long-term time trace of CCW/CW flagellar rotation state and run/tumble cell swimming state. CW intervals (purple, top) for each flagellum were determined from the fluorescence images. Tumbles (black, bottom) were determined from the cell-body rotation signal. (<bold>E</bold>) Mean deviation <italic>η</italic> from the veto model vs number of flagella per cell. Wild-type cells (solid black circles) with multiple flagella deviate significantly from the model (p=0.0003, <italic>N</italic> = 69 cells). CheY* cells (open gray circles; <italic>N</italic> = 46 cells) match the model (p=0.77). Error bars denote SEM. See ‘Materials and methods’ for more details.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.003">http://dx.doi.org/10.7554/eLife.01916.003</ext-link></p></caption><graphic xlink:href="elife01916f001"/></fig><fig id="fig1s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01916.004</object-id><label>Figure 1—figure supplement 1.</label><caption><title>Distribution of flagellar number.</title><p>Top, fluorescently-labeled flagella on trapped <italic>E. coli</italic> cells. Images show cells with 1, 2, 3, 4, and 6 flagella (left to right). The image contrast was adjusted to make flagella more visible (see ‘Materials and methods’ for details on the labeling, imaging and image processing protocols). Middle, transmission electron microscopy (TEM) images of cells possessing 1, 2, 3, 4 and 7 flagella (left to right). The images were enhanced to make flagella more visible (see ‘Materials and methods’ for details of the imaging and image processing). Bottom, distributions of the number of flagella per cell. Cells with zero flagella were excluded because non-swimming cells were not trapped. The distributions from trapped cells (blue, <italic>N</italic> = 86 cells) and TEM (red, <italic>N</italic> = 56 cells) have very similar means. Error bars denote SEM.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.004">http://dx.doi.org/10.7554/eLife.01916.004</ext-link></p></caption><graphic xlink:href="elife01916fs001"/></fig><fig id="fig1s2" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01916.005</object-id><label>Figure 1—figure supplement 2.</label><caption><title>Instrument layout.</title><p>Layout of the combined optical trap/epi-fluorescence microscope, showing the 1064-nm trapping laser beam (red), 532-nm fluorescence excitation laser beam (green), fluorescence emission path (yellow), objectives (01 and 02), quadrant photodiode (QPD), and fluorescence imaging, charge-coupled device camera (EMCCD). To the right is a schematic of a cell held by two optical traps (red cones) in the fluorescence excitation volume (green) within the sample chamber.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.005">http://dx.doi.org/10.7554/eLife.01916.005</ext-link></p></caption><graphic xlink:href="elife01916fs002"/></fig><fig id="fig1s3" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01916.006</object-id><label>Figure 1—figure supplement 3.</label><caption><title>Laser temporal interlacing scheme.</title><p>Schematic of the temporal interlacing of the trapping laser (red), excitation laser (green) and camera exposure (black). The 1064-nm trapping laser alternates with the 532-nm fluorescence excitation laser, such that they are never on at the same time. The camera exposure (30-μs exposure) is synchronized with the fluorescence excitation laser for stroboscopic imaging. Lower panel shows the interlacing over a longer time scale. The fluorescence excitation pulse and camera exposure occur once every 10 ms to create movies at 100 fps. See ‘Materials and methods’ for additional details.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.006">http://dx.doi.org/10.7554/eLife.01916.006</ext-link></p></caption><graphic xlink:href="elife01916fs003"/></fig><fig id="fig1s4" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01916.007</object-id><label>Figure 1—figure supplement 4.</label><caption><title>Sample data from representative cells.</title><p>Typical data from trapped cells with different numbers of flagella. A few still images (out of hundreds for each sample) of the fluorescently labeled flagella at different time points are shown. In the panels below, colors indicate the swimming state as runs/tumbles (white/black) and flagella waveforms as ‘normal’/‘semi-coiled’/‘curly-1’/‘semi-coiled’ or ‘curly-1’ (white/red/blue/purple). Light colors indicate periods when the flagellum was transitioning between two different waveforms. Light blue indicates a transition between normal and semi-coiled, light red indicates a transition between normal and curly-1. Samples A and B are wild-type cells, samples C–E are CheY* cells.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.007">http://dx.doi.org/10.7554/eLife.01916.007</ext-link></p></caption><graphic xlink:href="elife01916fs004"/></fig><fig id="fig1s5" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01916.008</object-id><label>Figure 1—figure supplement 5.</label><caption><title>Cells with flagella in the <italic>curly-1</italic> state rarely run.</title><p>Top, the fraction of time that cells spent running vs the number of flagella that were in the <italic>curly-1</italic> waveform (while all other flagella were in the <italic>normal</italic> waveform). When all flagella were in the <italic>normal</italic> waveform, cells ran 91% of the time. When the waveform of a single flagellum was <italic>curly-1</italic>, and the rest were <italic>normal</italic>, cells ran 18% of the time. Bottom, mean deviation from veto model, <italic>η</italic>, vs number of flagella per cell. Same as <xref ref-type="fig" rid="fig1">Figure 1E</xref>, with an additional theoretical model in which cells run 18% of the time when the waveform of a single flagellum is <italic>curly-1</italic> (magenta). Error bars denote SEM.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.008">http://dx.doi.org/10.7554/eLife.01916.008</ext-link></p></caption><graphic xlink:href="elife01916fs005"/></fig></fig-group></p><p>In this study, we present simultaneous, prolonged observations of individual flagella and whole-cell swimming. Using an optical trap to hold a swimming cell (<xref ref-type="bibr" rid="bib27">Min et al., 2009</xref>) while simultaneously imaging its fluorescently labeled flagella, we relate directly the number and state of each flagellum on a cell to its swimming behavior. Our measurements establish a simple mapping between CCW/CW flagellar rotation and run/tumble bacterial swimming state for cells with arbitrary numbers of flagella. Surprisingly, we find that <italic>E. coli</italic> cells wild-type for chemotaxis do not strictly comply with the veto model, because the states of individual flagella in the same cell are strongly coupled. Flagella do not switch independently, and as a result, tumbles typically involve many of the flagella on the cell. The behavior of multi-flagellated cells can still be mapped to the simple veto model by renormalizing the flagellar number to a lower effective number of independent flagella. As for the cause of inter-flagella coupling, our data strongly favor a mechanism involving fluctuations of the signaling network rather than direct, physical interactions between flagella. This model is supported by the observation that mutant cells in which flagellar switching is decoupled from the chemotaxis network do obey the simple veto model, as well as by stochastic simulations of the swimming behavior.</p></sec><sec id="s2" sec-type="results"><title>Results</title><sec id="s2-1"><title>The tumble bias of wild-type cells deviates from the predictions of the veto model</title><p>According to the veto model, a cell tumbles whenever any one flagellum rotates CW (<xref ref-type="fig" rid="fig1">Figure 1A</xref>). Thus, the probability that a cell runs equals the probability that all of its flagella remain CCW. As a consequence, cells with more flagella are expected to tumble more, since there is a higher chance that at least one flagellum will deviate from the consensus and ‘veto’ the run. These predictions can be stated mathematically, under the assumption that the rotational direction of each flagellum is independent of the other flagella. In a cell with <italic>N</italic><sub><italic>flag</italic></sub> flagella, the average tumble bias <italic>TB</italic>—the fraction of time a cell spends tumbling—will be given by<disp-formula id="equ1"><label>(1)</label><mml:math id="m1"><mml:mrow><mml:mi>T</mml:mi><mml:mi>B</mml:mi><mml:mo> </mml:mo><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mo> </mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mo> </mml:mo><mml:mi>C</mml:mi><mml:mi>B</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mrow><mml:mi>f</mml:mi><mml:mi>l</mml:mi><mml:mi>a</mml:mi><mml:mi>g</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msup></mml:mrow></mml:math></disp-formula>where <italic>CB</italic> is the average clockwise bias—the fraction of time the cell’s flagella rotate CW (‘Materials and methods’).</p><p>To test this prediction, we quantified the swimming behavior of individual <italic>E. coli</italic> cells wild-type for chemotaxis (strain HCB1660 [<xref ref-type="bibr" rid="bib39">Turner et al., 2010</xref>]; see ‘Materials and methods’, <xref ref-type="table" rid="tbl1">Table 1</xref>) using an instrument combining optical tweezers and epi-fluorescence imaging (<xref ref-type="fig" rid="fig1">Figure 1B</xref> and <xref ref-type="fig" rid="fig1s2">Figure 1—figure supplement 2</xref>). The instrument allowed us to measure simultaneously run/tumble behavior and flagellar dynamics in the same cell. The optical trap was used to hold each end of a single cell in place and the light scattered by the cell was utilized to monitor its swimming behavior, as described previously (<xref ref-type="bibr" rid="bib27">Min et al., 2009</xref>, <xref ref-type="bibr" rid="bib28">2012</xref>). As shown in <xref ref-type="fig" rid="fig1">Figure 1C</xref> (<xref ref-type="other" rid="video1">Video 1</xref>), cell runs were identified from oscillatory time signals due to cell body rotation at a frequency of ∼10 Hz. Cell tumbles were identified as periods of erratic motion during which the flagellar bundle was disrupted (<xref ref-type="bibr" rid="bib27">Min et al., 2009</xref>). Flagella were fluorescently labeled using the method of <xref ref-type="bibr" rid="bib39">Turner et al. (2010)</xref>. High speed, epi-fluorescent, stroboscopic imaging (<xref ref-type="bibr" rid="bib38">Turner et al., 2000</xref>; <xref ref-type="fig" rid="fig1s3">Figure 1—figure supplement 3</xref>) was used to resolve individual flagella (<xref ref-type="fig" rid="fig1">Figure 1C</xref>). Since the trapped cell remained in the field of view for a prolonged period, flagella were observed through multiple runs and tumbles (typically ∼5 events), limited by the time until flagella became too dim to discern due to photobleaching (∼8 to 40 s). The rotational direction of each flagellum was determined by observing its shape during 100-ms time windows. As shown by <xref ref-type="bibr" rid="bib13">Darnton et al. (2007)</xref>, flagella may take on different helical waveforms depending on their rotational state. These waveforms, termed ‘<italic>normal</italic>’, ‘<italic>semi-coiled</italic>’, ‘<italic>curly-1</italic>’, and ‘<italic>curly-2</italic>’, can be visually identified based on their pitch and wavelength (<xref ref-type="fig" rid="fig1">Figure 1C</xref>, <xref ref-type="fig" rid="fig1s4">Figure 1—figure supplement 4</xref>). CCW rotating flagella were identified based on the <italic>normal</italic> conformation, which they have been shown to adopt exclusively (<xref ref-type="bibr" rid="bib13">Darnton et al., 2007</xref>), while CW rotating flagella were identified by their <italic>curly-1</italic> or <italic>semi-coiled</italic> shape. From the identification of CCW and CW intervals, the cell’s mean CW bias was determined by averaging the fraction of time that all the flagella on the cell spent CW (‘Materials and methods’).<table-wrap id="tbl1" position="float"><object-id pub-id-type="doi">10.7554/eLife.01916.009</object-id><label>Table 1.</label><caption><p>Strains and plasmids used in this work</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.009">http://dx.doi.org/10.7554/eLife.01916.009</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><th>Strain</th><th>Genotype</th><th>Comments</th><th>Source</th></tr></thead><tbody><tr><td>HCB1660</td><td><italic>fliC::</italic>Tn5 (Kan<sup>R</sup>)</td><td>‘wild type’ Contains plasmid pBAD33-fliC<sup>S219C</sup></td><td>(<xref ref-type="bibr" rid="bib39">Turner et al., 2010</xref>) Gift of H Berg</td></tr><tr><td>PM87</td><td><italic>cheBYZ::</italic>FRT<italic>, fliC::</italic>Tn5 (Kan<sup>R</sup>)</td><td>‘CheY*’ Contains plasmids pMS164 and pPM5</td><td>This study</td></tr><tr><td>RP437</td><td/><td>Wild-type for chemotaxis</td><td>(<xref ref-type="bibr" rid="bib45">Parkinson and Houts, 1982</xref>)</td></tr><tr><td>SK109</td><td><italic>cheBYZ::</italic>Cm</td><td/><td>This study</td></tr><tr><td>SK110</td><td><italic>cheBYZ::</italic>FRT</td><td/><td>This study</td></tr><tr><td>SK112</td><td><italic>cheBYZ::</italic>FRT<italic>, fliC::</italic>Tn5 (Kan<sup>R</sup>)</td><td/><td>This study</td></tr><tr><td>Plasmids</td><td/><td/><td/></tr><tr><td> pBAD33 fliC<sup>S219C</sup></td><td><italic>fliC</italic><sup><italic>S219C</italic></sup> under P<sub>araBAD</sub> promoter, Cm<sup>R</sup>, p15a origin</td><td>Expresses mutant version of FliC for fluorescent labeling</td><td>(<xref ref-type="bibr" rid="bib39">Turner et al., 2010</xref>) Gift of H Berg</td></tr><tr><td> pPM5</td><td><italic>fliC</italic><sup><italic>S219C</italic></sup> under P<sub>araBAD</sub> promoter, Amp<sup>R</sup>, colE1 origin</td><td>Expresses mutant version of FliC for fluorescent labeling</td><td>This study</td></tr><tr><td> pMS164</td><td><italic>cheY</italic><sup><italic>D13K</italic></sup> under P<sub>lacOP</sub> promoter, Cm<sup>R</sup>, pSC101 origin</td><td>Expresses constitutively active version of CheY</td><td>(<xref ref-type="bibr" rid="bib1">Alon et al., 1998</xref>) Gift of P Cluzel</td></tr><tr><td> pDK46</td><td/><td>Helper plasmid</td><td>(<xref ref-type="bibr" rid="bib14">Datsenko and Wanner, 2000</xref>)</td></tr><tr><td> pKD3</td><td/><td>Template for Cm<sup>R</sup> cassette</td><td>(<xref ref-type="bibr" rid="bib14">Datsenko and Wanner, 2000</xref>)</td></tr><tr><td> pCP20</td><td/><td>Helper plasmid</td><td>(<xref ref-type="bibr" rid="bib9">Cherepanov and Wackernagel, 1995</xref>)</td></tr></tbody></table></table-wrap><media content-type="glencoe play-in-place height-250 width-310" id="video1" mime-subtype="mov" mimetype="video" xlink:href="elife01916v001.mov"><object-id pub-id-type="doi">10.7554/eLife.01916.010</object-id><label>Video 1.</label><caption><title>Video of trapped wild-type cell with three labeled flagella as it runs and tumbles.</title><p>Slow motion video of the wild-type cell in <xref ref-type="fig" rid="fig1">Figure 1C</xref> with three long, fluorescently labeled flagella. The approximate location of the unlabeled cell body is indicated by the dotted line. The trap signal used to determine runs and tumbles (bottom, scrolling blue curve) measures the position of the cell body in the trap as it rotates. At the beginning of the video (time stamp = 2.4 s), all three flagella are in a bundle and the cell is running. One by one, the flagella switch to CW rotation (2.7–3.0 s), which disrupts the bundle and causes the cell to tumble. Flagella can be observed in all three waveforms, <italic>normal, semi-coiled,</italic> and <italic>curly-1</italic>. Near the end of the video, the flagella all return to CCW rotation and coalesce into a bundle, causing the cell to resume running (3.4 s). In addition to the three long flagella, a short flagellar stub is visible. The stub does not affect the swimming behavior and is not analyzed. Scale bar in bottom left corner is 1 µm. Frames were recorded at 400 frames per second, video shows every other frame at 20 frames per second.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.010">http://dx.doi.org/10.7554/eLife.01916.010</ext-link></p></caption></media></p><p>Our assay allowed us to determine all the parameters in <xref ref-type="disp-formula" rid="equ1">Equation (1)</xref> directly. For each cell, we measured the tumble bias (using the optical trap), flagellar number <italic>N</italic><sub><italic>flag</italic></sub>, and CW bias (using fluorescence imaging). We used these values to compare our experimental data to the prediction of the veto model. Reorganizing <xref ref-type="disp-formula" rid="equ1">Equation (1)</xref>, we define the parameter <italic>η</italic>:<disp-formula id="equ2"><label>(2)</label><mml:math id="m2"><mml:mrow><mml:mo> </mml:mo><mml:mi>η</mml:mi><mml:mo> </mml:mo><mml:mo>≡</mml:mo><mml:mfrac><mml:mrow><mml:mtext>log</mml:mtext><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mo> </mml:mo><mml:mi>T</mml:mi><mml:mi>B</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mtext>log</mml:mtext><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>C</mml:mi><mml:mi>B</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mo> </mml:mo><mml:mo>−</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mrow><mml:mi>f</mml:mi><mml:mi>l</mml:mi><mml:mi>a</mml:mi><mml:mi>g</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></disp-formula>which quantifies the deviation of the data from the veto model. Comparing <xref ref-type="disp-formula" rid="equ2">Equation (2)</xref> to <xref ref-type="disp-formula" rid="equ1">Equation (1)</xref>, <italic>η</italic> may also be interpreted as the difference between two terms: the number of flagella estimated from the veto model based on the cell’s swimming behavior, and the number of flagella on the cell as determined by counting directly. <xref ref-type="fig" rid="fig1">Figure 1E</xref> (black circles) shows <italic>η</italic> against the flagellar number <italic>N</italic><sub><italic>flag</italic></sub>. <italic>η</italic> was calculated for each individual cell and then averaged over all cells with a given number of flagella. Unexpectedly, we found that wild-type cells with multiple flagella systematically deviated from the predicted behavior. Specifically, <italic>η</italic> was consistently negative for cells with <italic>N</italic><sub><italic>flag</italic></sub> &gt; 1 (35/48 cells), indicating that cells with multiple flagella tumbled less than expected from the model. In the context of the veto model, the cells behaved as if they had a smaller number of flagella than what they actually had.</p><p>We first considered the possibility that a more detailed version of the veto model might explain the observed behavior and reconcile this discrepancy. A recent study by <xref ref-type="bibr" rid="bib33">Sneddon et al. (2012)</xref> used the observations of <xref ref-type="bibr" rid="bib13">Darnton et al. (2007)</xref> to refine the veto model. Specifically, the Sneddon model states that a cell with a minimum of <italic>X</italic> CCW flagella will run rather than tumble, provided the remaining (<italic>N</italic><sub><italic>flag</italic></sub>–<italic>X</italic>) CW flagella are in the <italic>curly-1</italic> conformation only (<italic>X</italic> is a parameter in the model, with possible values in the range [1, <italic>N</italic><sub><italic>flag</italic></sub> − 1]). Thus, the simple veto model considered above corresponds to <italic>X</italic> = <italic>N</italic><sub><italic>flag</italic></sub>, and the least perturbative refinement to the model corresponds to <italic>X</italic> = <italic>N</italic><sub><italic>flag</italic></sub> − 1, in which a cell with a single <italic>curly-1</italic> flagellum still runs. However, in our measurements we observed that cells with a single CW flagellum in the <italic>curly-1</italic> state still tumbled 82% of the time (44 s of cumulative time in which one flagellum was in the <italic>curly-1</italic> state; <xref ref-type="fig" rid="fig1s5">Figure 1—figure supplement 5</xref>). Modifying the Sneddon model to allow runs 18% of the time was not sufficient to reproduce the trend observed in <xref ref-type="fig" rid="fig1">Figure 1D</xref>, <xref ref-type="fig" rid="fig1s5">Figure 1—figure supplement 5</xref>.</p></sec><sec id="s2-2"><title>Wild-type tumbles typically involve multiple flagella</title><p>To investigate the discrepancy between our data and the veto model, we next examined individual tumble events in greater detail. In agreement with the original observation of <xref ref-type="bibr" rid="bib38">Turner et al. (2000)</xref>, we found that CW rotation of a single flagellum was indeed sufficient to cause a tumble in multi-flagellated cells (<xref ref-type="fig" rid="fig1s4">Figure 1—figure supplement 4</xref>, samples A, D and E). However, we also observed that more than half of tumbles in multi-flagellated cells (56%, 117/210 events) actually involved multiple CW flagella. <xref ref-type="fig" rid="fig2">Figure 2A</xref> shows a representative trace from a wild-type cell with three flagella. There are times during each tumble in the trace where all three flagella are in a CW state. As shown in <xref ref-type="fig" rid="fig2">Figure 2C</xref>, the number of CW flagella ‘participating’ in a tumble (black circles) was significantly larger than would be expected if flagella were independently switching (gray dashed line, obtained from simulations of a cell with independent flagella; see ‘Materials and methods’). Our results thus suggest that while a single CW flagellum is sufficient to induce a tumble (in agreement with a simple veto model), flagella are also coupled and may thus switch in groups, in a correlated fashion. Further evidence for inter-flagella coupling was obtained by calculating the cross-correlation between pairs of flagella on a given cell. We found a significant correlation between the rotational directions of pairs of flagella on the same cell (<xref ref-type="fig" rid="fig2">Figure 2D</xref>, black data points). This correlation persisted for ∼1 s, the average duration of a tumble. Our findings are consistent with previous observations by Terasawa et al. on surface-immobilized cells (<xref ref-type="bibr" rid="bib36">Terasawa et al., 2011</xref>). There, correlations between individual motors on the same cell were detected by monitoring beads attached to flagellar stubs, as opposed to complete flagella on swimming cells in our present work.<fig-group><fig id="fig2" position="float"><object-id pub-id-type="doi">10.7554/eLife.01916.011</object-id><label>Figure 2.</label><caption><title>Tumbles in wild-type cells involve multiple CW flagella.</title><p>(<bold>A</bold>) A typical time trace from a wild-type cell with 3 flagella. Colors indicate runs/tumbles (white/black, top) and CCW/CW (white/purple, middle). The blue line (bottom) shows the corresponding number of CW flagella at each time point. (<bold>B</bold>) Same as (<bold>A</bold>) for a typical CheY* cell with 3 flagella. (<bold>C</bold>) Mean of the maximum number of CW flagella during a tumble vs number of flagella per cell. Consistently more flagella are CW during tumbles in the wild-type (black circles; <italic>N</italic> = 61 cells) compared to the CheY* strain (open gray circles; <italic>N</italic> = 24 cells). Simulations incorporating fluctuations in CheY-P (black line) and without fluctuations (gray dashed lines) reproduce the observed trends (simulations detailed in the text and <xref ref-type="fig" rid="fig4">Figure 4</xref>). (<bold>D</bold>) Cross-correlation between flagella pairs, averaged over all pairs and all cells. Wild-type (black circles) match simulations with fluctuations in CheY-P (black line). CheY* strain (open gray circles) matches simulations without fluctuations in CheY-P (gray dashed line), which exhibit almost no correlation. Error bars denote SEM. See ‘Materials and methods’ for more details.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.011">http://dx.doi.org/10.7554/eLife.01916.011</ext-link></p></caption><graphic xlink:href="elife01916f002"/></fig><fig id="fig2s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01916.012</object-id><label>Figure 2—figure supplement 1.</label><caption><title>Flagellar transition rates vs number of flagella per cell.</title><p>Flagella motor transition rates, separated into two groups corresponding to cells with a single flagellum (red) and cells with multiple flagella (blue). All data from wild-type cells (<italic>N</italic> = 52 cells). Error bars denote SEM.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.012">http://dx.doi.org/10.7554/eLife.01916.012</ext-link></p></caption><graphic xlink:href="elife01916fs006"/></fig><fig id="fig2s2" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01916.013</object-id><label>Figure 2—figure supplement 2.</label><caption><title>Flagellar transition rates vs number of flagella per cell.</title><p>Transitions between the three observed flagella waveforms (<italic>normal</italic>, <italic>semi-coiled</italic> and <italic>curly-1</italic>), separated into two groups corresponding to cells with a single flagellum (red) and cells with multiple flagella (blue). All data from wild-type cells (<italic>N</italic> = 52 cells). Error bars denote SEM.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.013">http://dx.doi.org/10.7554/eLife.01916.013</ext-link></p></caption><graphic xlink:href="elife01916fs007"/></fig></fig-group></p><p>The source of inter-flagellar correlation remains under debate. Terasawa et al. observed that mutant cells, in which the concentration of CheY-P was decoupled from the chemotaxis network, displayed no correlations (<xref ref-type="bibr" rid="bib36">Terasawa et al., 2011</xref>). This led us to likewise investigate the behavior of a strain, PM87, expressing a constitutively-active CheY (CheY<sup>D13K</sup> [<xref ref-type="bibr" rid="bib1">Alon et al., 1998</xref>], denoted CheY*, see <xref ref-type="table" rid="tbl1">Table 1</xref> and <xref ref-type="table" rid="tbl2">2</xref>). The protein was exogenously expressed, with the expression level chosen such that the population-averaged tumble bias matched that of wild-type cells (‘Materials and methods’). A representative trace from a CheY* cell with three flagella is shown in <xref ref-type="fig" rid="fig2">Figure 2B</xref> (see also <xref ref-type="fig" rid="fig1s4">Figure 1—figure supplement 4</xref>, and <xref ref-type="other" rid="video2">Video 2</xref>). Upon inspection, flagellar switching appears far less correlated than in wild-type cells (<xref ref-type="fig" rid="fig2">Figure 2D</xref>, compare black and gray data points). Comparing 54 wild-type and 24 CheY* cells with the same mean CW biases (0.11 ± 0.07 vs 0.11 ± 0.07, mean ± SD) we found that, on average, fewer CW rotating flagella participated in tumbles in the CheY* strain (<xref ref-type="fig" rid="fig2">Figure 2C</xref>, open circles). Moreover, the number of participating flagella in the CheY* strain closely matched the expectation for cells with independently switching flagella (<xref ref-type="fig" rid="fig2">Figure 2C</xref>, dashed line). (This number deviates from unity and trends upwards with number of flagella simply because of the finite probability that two tumbles overlap by chance.) Our results indicate that when the signal for flagellar motors to switch their rotational state is decoupled from the chemotaxis network, the motors switch independently. Based on our interpretation of the wild-type data, we thus expect CheY* cells to adhere to the simple veto model. As shown in the plot of <italic>η</italic> in <xref ref-type="fig" rid="fig1">Figure 1D</xref> (open circles), CheY* cells indeed match the veto model closely (<italic>η</italic> = −0.08 ± 0.15, mean ± SEM). The existence of correlations between flagella states in wild-type cells may thus explain why cells with multiple flagella deviate from the veto model.<table-wrap id="tbl2" position="float"><object-id pub-id-type="doi">10.7554/eLife.01916.014</object-id><label>Table 2.</label><caption><p>Primers used in this work</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.014">http://dx.doi.org/10.7554/eLife.01916.014</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><th>Primer</th><th>Sequence</th><th>Comments</th></tr></thead><tbody><tr><td>SK140F</td><td>TGCGTGGTCAGACGGTGTATGCGCTAAGTAAGGATTAACG GTGTAGGCTGGAGCTGCTTC</td><td><italic>cheBYZ</italic> deletion forward</td></tr><tr><td>SK140R</td><td>GCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAA CATATGAATATCCTCCTTAG</td><td><italic>cheBYZ</italic> deletion reverse</td></tr><tr><td>SK141F</td><td>CCTTAAACCCGACGGATTGC</td><td><italic>cheBYZ</italic> deletion check forward</td></tr><tr><td>SK141R</td><td>TTGCTGCCACACATCAAGC</td><td><italic>cheBYZ</italic> deletion check reverse</td></tr><tr><td>SK163F</td><td>AGGGTTATTGTCTCATGAGC</td><td>pZE11 sequencing forward</td></tr><tr><td>SK163R</td><td>GTTTTATTTGATGCCTCTAG</td><td>pZE11 sequencing reverse</td></tr><tr><td>PM7F</td><td>GGG GACGTC ATCGATGCATAATGTGCCTG</td><td>amplify P<sub>araBAD</sub> <italic>fliC</italic><sup><italic>S219C</italic></sup> forward</td></tr><tr><td>PM7R</td><td>GGG GTCGAC TTAACCCTGCAGC</td><td>amplify P<sub>araBAD</sub> <italic>fliC</italic><sup><italic>S219C</italic></sup> reverse</td></tr></tbody></table></table-wrap><media content-type="glencoe play-in-place height-250 width-310" id="video2" mime-subtype="mov" mimetype="video" xlink:href="elife01916v002.mov"><object-id pub-id-type="doi">10.7554/eLife.01916.015</object-id><label>Video 2.</label><caption><title>Video of trapped CheY* cell with two labeled flagella as it runs and tumbles.</title><p>Slow motion video, similar to <xref ref-type="other" rid="video1">Video 1</xref>, of a trapped CheY* cell with two fluorescently labeled flagella. Still images from this cell are shown in <xref ref-type="fig" rid="fig1s4">Figure 1—figure supplement 4</xref>, sample D. The approximate location of the unlabeled cell body is indicated by the dotted line. At the beginning of the video (time stamp = 7.4 s), both flagella are in a bundle, rotating CCW in the <italic>normal</italic> waveform, and the cell is running. At 7.6 s, one flagellum switches to CW rotation and transitions to the <italic>semi-coiled</italic> waveform, which disrupts the bundle and causes the cell to tumble. At the 8.5 s, the <italic>semi-coiled</italic> flagellum returns to CCW rotation and both flagella coalesce into a bundle, causing the cell to resume running. Scale bar in bottom left corner is 1 µm. Frames were recorded at 100 frames per second, video plays at 20 frames per second.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.015">http://dx.doi.org/10.7554/eLife.01916.015</ext-link></p></caption></media></p></sec><sec id="s2-3"><title>Wild-type behavior can be described by a veto model with a lower effective number of flagella</title><p>Our results so far suggest that, while wild-type cells obey the fundamental premise of the veto model—that is, a single CW flagellum is sufficient to induce a tumble—the presence of inter-flagella correlations leads to the failure of <xref ref-type="disp-formula" rid="equ1">Equation (1)</xref> in relating the observed CW bias and tumble bias. To describe the relation between single flagellar state and whole-cell behavior successfully, this expression must then be modified to account for flagellar correlations. To this end, we examined the relation between CW bias and tumble bias in all individual cells having a given flagellar number (<xref ref-type="fig" rid="fig3">Figure 3A,B</xref>). <xref ref-type="disp-formula" rid="equ1">Equation (1)</xref> defines a single curve, along which the CW bias and tumble bias of all cells with <italic>N</italic><sub><italic>flag</italic></sub> flagella should lie (dashed line in <xref ref-type="fig" rid="fig3">Figure 3A,B</xref>). As expected, the CheY* cells follow these predicted curves closely for all <italic>N</italic><sub><italic>flag</italic></sub> values (<xref ref-type="fig" rid="fig3">Figure 3B</xref>) (R<sup>2</sup> = 0.89). In contrast, wild-type cells with multiple flagella consistently fall below the predicted curves (<xref ref-type="fig" rid="fig3">Figure 3A</xref>) (35/48 cells).<fig-group><fig id="fig3" position="float"><object-id pub-id-type="doi">10.7554/eLife.01916.016</object-id><label>Figure 3.</label><caption><title>Wild-type behavior matches the veto model for cells with a lower effective number of flagella.</title><p>(<bold>A</bold>) Tumble bias vs CW bias for individual wild-type cells (<italic>N</italic> = 69), plotted separately for different numbers of flagella per cell (<italic>N</italic><sub><italic>flag</italic></sub> = 1, purple; 2, blue; 3, green; 4, red; 5, cyan). The prediction from the veto model in <xref ref-type="disp-formula" rid="equ1">Equation (1)</xref> (dashed lines) does not match the data for cells with multiple flagella (R<sup>2</sup> = 0.88, 0.60, 0.41, 0.39 for <italic>N</italic><sub><italic>flag</italic></sub> = 2, 3, 4, 5). The data were fit (solid lines) to <xref ref-type="disp-formula" rid="equ1">Equation (1)</xref>, while allowing the number of flagella to be used as a fitting parameter, <italic>N</italic><sub><italic>eff</italic></sub>. Error bars denote SD. (<bold>B</bold>) Same as (<bold>A</bold>) for CheY* (open circles, same color code as [<bold>A</bold>] <italic>N</italic> = 46 cells). The veto model prediction (dashed lines) matches the data well (R<sup>2</sup> = 0.91, 0.97, 0.93, 0.67, 0.98 for <italic>N</italic><sub><italic>flag</italic></sub> = 1, 2, 3, 4, 5). Fits (solid lines) yield <italic>N</italic><sub><italic>eff</italic></sub> values almost identical to <italic>N</italic><sub><italic>flag</italic></sub>. (<bold>C</bold>) Fitted <italic>N</italic><sub><italic>eff</italic></sub> values vs number of flagella per cell for wild-type (black circles) and CheY* (open gray circles) cells. Simulations (described in the text) reproduce the observed trends. (<bold>D</bold>) Data points from individual wild-type (solid circles) and CheY* (open circles) cells all collapse onto a single line when using <italic>N</italic><sub><italic>eff</italic></sub> from fits to wild-type data in (<bold>A</bold>) and the actual flagellar number <italic>N</italic><sub><italic>flag</italic></sub> for CheY* cells in (<bold>B</bold>). Error bars denote SEM. See ‘Materials and methods’ for more details.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.016">http://dx.doi.org/10.7554/eLife.01916.016</ext-link></p></caption><graphic xlink:href="elife01916f003"/></fig><fig id="fig3s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01916.017</object-id><label>Figure 3—figure supplement 1.</label><caption><title>Fit to <italic>N</italic><sub><italic>eff</italic></sub> vs <italic>Nf</italic><sub><italic>lag</italic></sub>.</title><p>Data points show fitted <italic>N</italic><sub><italic>eff</italic></sub> values vs number of flagella per cell for wild-type cells as in <xref ref-type="fig" rid="fig3">Figure 3C</xref>. The data was fit empirically to a power law (red line) given by <italic>N</italic><sub><italic>eff</italic></sub> = 1.27 × <italic>N</italic><sub><italic>flag</italic></sub><sup>0.5</sup>.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.017">http://dx.doi.org/10.7554/eLife.01916.017</ext-link></p></caption><graphic xlink:href="elife01916fs008"/></fig></fig-group></p><p>Based on our observation that wild-type cells exhibited <italic>η</italic> values consistent with cells with a lower number of flagella than the actual value, we hypothesized that wild-type behavior may be described within the framework of the veto model by allowing the parameter <italic>N</italic><sub><italic>flag</italic></sub> in <xref ref-type="disp-formula" rid="equ1">Equation (1)</xref> to deviate from the actual flagellar number. As shown in <xref ref-type="fig" rid="fig3">Figure 3A</xref> (solid lines), using the flagellar number as a fitting parameter (now denoted <italic>N</italic><sub><italic>eff</italic></sub>) indeed allows for a good match for the wild-type data (R<sup>2</sup> = 0.85). In this revised equation, <italic>N</italic><sub><italic>eff</italic></sub> can be thought of as the ‘effective’ number of independent flagella on a cell, which captures the fact that flagella in wild-type cells switch in a correlated manner. Consistent with this picture, the effective number of flagella <italic>N</italic><sub><italic>eff</italic></sub> is consistently smaller than the actual flagellar number <italic>N</italic><sub><italic>flag</italic></sub> (<xref ref-type="fig" rid="fig3">Figure 3C</xref>, black circles; <italic>N</italic><sub><italic>eff</italic></sub> can be approximated by <italic>N</italic><sub><italic>eff</italic></sub> = 1.27 × <italic>N</italic><sub><italic>flag</italic></sub><sup>0.5</sup>, see <xref ref-type="fig" rid="fig3s1">Figure 3—figure supplement 1</xref>). As a control, estimating <italic>N</italic><sub><italic>eff</italic></sub> for CheY* cells produces values very close to the original flagellar number <italic>N</italic><sub><italic>flag</italic></sub> (<xref ref-type="fig" rid="fig3">Figure 3B</xref>, solid line and <xref ref-type="fig" rid="fig3">Figure 3C</xref>, open circles). The introduction of the parameter <italic>N</italic><sub><italic>eff</italic></sub> allows us to formulate a generalized veto model, which describes the mapping between the CW bias and tumble bias for both wild-type and CheY* cells. The generalized model now defines a universal curve,<disp-formula id="equ3"><label>(3)</label><mml:math id="m3"><mml:mrow><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>T</mml:mi><mml:mi>B</mml:mi><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>C</mml:mi><mml:mi>B</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi mathvariant="italic">eff</mml:mi></mml:msub></mml:mrow></mml:msup></mml:mrow></mml:math></disp-formula>along which all individual cells of both genotypes should lie (using <italic>N</italic><sub><italic>eff</italic></sub> = <italic>N</italic><sub><italic>flag</italic></sub> for CheY* strain and the best fit value of <italic>N</italic><sub><italic>eff</italic></sub> for wild-type). As seen in <xref ref-type="fig" rid="fig3">Figure 3D</xref>, this expression successfully collapses all single-cell data from both strains and all flagellar numbers.</p></sec><sec id="s2-4"><title>A theoretical model incorporating CheY-P fluctuations reproduces wild-type swimming behavior</title><p>Our results show that <italic>E. coli</italic> cells adhere to the veto model, but that inter-flagellar correlations lead to a renormalization of the effect of flagellar number. What is the source of these flagellar correlations? The absence of correlations in the CheY* strain and its adherence to a simple veto model provide an important clue to understanding the mechanism of inter-flagellar coupling. In wild-type cells, CheY-P levels are subject to phosphorylation and de-phosphorylation reactions by chemotaxis network components and are believed to fluctuate in time (<xref ref-type="bibr" rid="bib23">Korobkova et al., 2006</xref>; <xref ref-type="bibr" rid="bib33">Sneddon et al., 2012</xref>). In contrast, in the mutant strain, CheY* levels are decoupled from the network and are thus expected to be constant over the timescales of interest (<xref ref-type="bibr" rid="bib23">Korobkova et al., 2006</xref>; <xref ref-type="bibr" rid="bib27">Min et al., 2009</xref>). Terasawa et al. proposed that fluctuations in CheY-P levels may thus provide a mechanism by which the CW biases of multiple flagella on a cell can be coupled (<xref ref-type="bibr" rid="bib36">Terasawa et al., 2011</xref>). To test whether such a mechanism could account for the different features of bacterial swimming observed in our study, we performed simulations of whole-cell swimming driven by the chemotaxis network. In particular, we investigated how fluctuations in CheY-P concentration could produce differences in inter-flagellar correlations between wild-type and CheY* cells and consequent differences in their respective mappings of CW bias to tumble bias.</p><p>We performed stochastic simulations of the chemotaxis network and resulting flagellar motor activity and then applied the veto rule that CW rotation of a single flagellum leads to cell tumbling (<xref ref-type="fig" rid="fig4">Figure 4A,B</xref>; ‘Materials and methods’). For wild-type cells, we incorporated fluctuations in CheY-P concentration using the approach of <xref ref-type="bibr" rid="bib33">Sneddon et al. (2012)</xref>. The model was constrained using transition rates between flagella waveforms that were directly extracted from our experiments (<xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1</xref> and <xref ref-type="table" rid="tbl3 tbl4">Tables 3 and 4</xref>). We found that, consistently, simulations with CheY-P fluctuations were able to duplicate the observed wild-type behavior, while simulations without fluctuations matched CheY* data and the predictions of the ‘naïve’ veto model. Specifically, using only two parameters—the amplitude and the characteristic timescale for CheY-P fluctuations (<xref ref-type="table" rid="tbl3">Table 3</xref>)—as fitting parameters, our simulations simultaneously reproduced (i) the relation between flagellar number and flagellar participation in tumbles (<xref ref-type="fig" rid="fig2">Figure 2C</xref>); (ii) temporal correlations between flagella (<xref ref-type="fig" rid="fig2">Figure 2D</xref>); (iii) the effective number of flagella in multi-flagellated wild-type cells (<xref ref-type="fig" rid="fig3">Figure 3C</xref>); and (iv) the degree of deviation from the veto model (<xref ref-type="fig" rid="fig4">Figure 4C</xref>).<fig-group><fig id="fig4" position="float"><object-id pub-id-type="doi">10.7554/eLife.01916.018</object-id><label>Figure 4.</label><caption><title>A theoretical model incorporating CheY-P fluctuations reproduces wild-type data.</title><p>(<bold>A</bold>) Simulated time trace for a wild-type cell. Representative simulated time trace of CheY-P concentration (gray line, top), CW bias (red), run/tumble state (white/black), CCW/CW flagellar rotational direction (white/purple) and number of CW flagella (blue line, bottom) for a cell with 3 flagella. CheY-P simulation parameters are described in the text and in <xref ref-type="table" rid="tbl3">Table 3</xref>. (<bold>B</bold>) Same as (<bold>A</bold>) from a simulated CheY* cell, in which CheY-P concentration (gray line) does not fluctuate. (<bold>C</bold>) Deviation from veto model. The theoretical model that includes CheY-P concentration fluctuations (black line) reproduces the wild-type data (black circles). A simple veto model with constant CheY-P concentration (gray dashed line) reproduces the CheY* data (open gray circles). Error bars denote SEM. See ‘Materials and methods’ for more details.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.018">http://dx.doi.org/10.7554/eLife.01916.018</ext-link></p></caption><graphic xlink:href="elife01916f004"/></fig><fig id="fig4s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01916.019</object-id><label>Figure 4—figure supplement 1.</label><caption><title>Flagellar waveform transition rates.</title><p>Diagram of the transitions rates between different flagellar waveforms: normal (CCW), semi-coiled and curly-1 (both CW). Data from wild-type cells (N = 52 cells, 203 tumbles). Values are mean ± SEM. Arrow thickness is proportional to the transition rate.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.019">http://dx.doi.org/10.7554/eLife.01916.019</ext-link></p></caption><graphic xlink:href="elife01916fs009"/></fig><fig id="fig4s2" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01916.020</object-id><label>Figure 4—figure supplement 2.</label><caption><title>Flagellar waveform transition sequences.</title><p>Histogram of the sequence of flagellar waveforms when flagella motor rotation switches from CW to CCW (N = <italic>normal</italic>, S = <italic>semi-coiled</italic>, C = <italic>curly-1</italic>). Note that prior to CW rotation all flagella had the <italic>normal</italic> waveform. The majority of CCW-to-CW switches caused flagella to transition to the <italic>semi-coiled</italic> waveform (77%, 164/213 of all flagellar switching events). Of those flagella, most remained <italic>semi-coiled</italic> until the motor returned to CCW (71%, 117/164), while some instead transitioned to the <italic>curly-1</italic> waveform (29%, 47/164). Some fraction of flagella (22%, 47/213) also skipped the <italic>semi-coiled</italic> waveform altogether and transitioned immediately from CCW to <italic>curly-1</italic>. Those flagella never switch to the <italic>semi-coiled</italic> state.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.020">http://dx.doi.org/10.7554/eLife.01916.020</ext-link></p></caption><graphic xlink:href="elife01916fs010"/></fig><fig id="fig4s3" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01916.021</object-id><label>Figure 4—figure supplement 3.</label><caption><title>Flagellar waveform transition sequences.</title><p>Decision tree showing the probability of each sequence of flagellar waveforms when flagella switch from CCW to CW rotation (N = <italic>normal</italic>, S = <italic>semi-coiled</italic>, C = <italic>curly-1</italic>). Arrow thickness is proportional to the probability of a specific transition.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.021">http://dx.doi.org/10.7554/eLife.01916.021</ext-link></p></caption><graphic xlink:href="elife01916fs011"/></fig><fig id="fig4s4" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01916.022</object-id><label>Figure 4—figure supplement 4.</label><caption><title>Tumble bias vs number of flagella per cell.</title><p>Mean tumble bias of wild-type cells vs number of flagella per cell. Data were filtered to exclude outliers with a CW bias outside the range of 0–0.25 (<italic>N</italic> = 61 cells remaining). Gray line shows the mean tumble bias from all cells. Error bars denote SEM.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.022">http://dx.doi.org/10.7554/eLife.01916.022</ext-link></p></caption><graphic xlink:href="elife01916fs012"/></fig></fig-group><table-wrap id="tbl3" position="float"><object-id pub-id-type="doi">10.7554/eLife.01916.023</object-id><label>Table 3.</label><caption><p>Model parameters</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.023">http://dx.doi.org/10.7554/eLife.01916.023</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><th>Parameter</th><th>Description</th><th>Value</th><th>Source</th></tr></thead><tbody><tr><td><italic>k</italic><sub><italic>ccw</italic>-&gt;<italic>cw</italic></sub></td><td>Motor switching rate from CCW→CW</td><td>0.26 s<sup>−1</sup></td><td>Our data</td></tr><tr><td><italic>k</italic><sub><italic>cw</italic>-&gt;<italic>ccw</italic></sub></td><td>Motor switching rate from CW→CCW</td><td>1.7 s<sup>−1</sup></td><td>Our data</td></tr><tr><td><italic>CB</italic></td><td>Average clock-wise bias of wild-type motors</td><td>0.13</td><td>Our data</td></tr><tr><td>&lt;<italic>Y</italic>&gt;</td><td>Mean concentration of CheY-P</td><td>2.59 µM</td><td>Fit to our data</td></tr><tr><td><italic>ω</italic></td><td>Characteristic motor switching time</td><td>0.5 s</td><td>Our data</td></tr><tr><td><italic>λ</italic><sup>−1</sup></td><td>Transition rate from <italic>semi-coiled</italic> to <italic>curly-1</italic> state</td><td>0.68 s<sup>−1</sup></td><td>Our data</td></tr><tr><td><italic>x</italic></td><td>From the model in Sneddon et al.; number of flagella that must be <italic>normal</italic> for a run to occur (while other flagella are <italic>curly-1</italic>) (<xref ref-type="bibr" rid="bib33">Sneddon et al., 2012</xref>)</td><td><italic>N</italic><sub><italic>flag</italic></sub> (variable in <xref ref-type="fig" rid="fig1s5">Figure 1—figure supplement 5</xref>)</td><td>Our data</td></tr><tr><td><italic>σ</italic><sup><italic>2</italic></sup></td><td>Variance in [CheY-P]</td><td>1.0 µM<sup>2</sup></td><td>Fit to our data</td></tr><tr><td><italic>τ</italic></td><td>Characteristic time-scale of fluctuations in [CheY-P]</td><td>0.2 s</td><td>Fit to our data</td></tr><tr><td><italic>K</italic><sub><italic>d</italic></sub></td><td>Midpoint of CW bias vs CheY-P response curve</td><td>3.1 µM</td><td>(<xref ref-type="bibr" rid="bib10">Cluzel et al., 2000</xref>)</td></tr><tr><td><italic>H</italic></td><td>Hill coefficient for CW bias vs CheY-P response curve</td><td>10.3</td><td>(<xref ref-type="bibr" rid="bib10">Cluzel et al., 2000</xref>)</td></tr><tr><td><italic>dt</italic></td><td>Simulation time steps</td><td>0.001 s</td><td>–</td></tr></tbody></table><table-wrap-foot><fn><p>Note: <inline-formula><mml:math id="inf1"><mml:mrow><mml:mi>ω</mml:mi><mml:mo>≡</mml:mo><mml:mfrac><mml:mrow><mml:mi>C</mml:mi><mml:mi>B</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mrow><mml:mtext>ccw</mml:mtext><mml:mo>→</mml:mo><mml:mtext>cw</mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mrow></mml:math></inline-formula></p></fn></table-wrap-foot></table-wrap><table-wrap id="tbl4" position="float"><object-id pub-id-type="doi">10.7554/eLife.01916.024</object-id><label>Table 4.</label><caption><p>Flagella waveform transition rates</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01916.024">http://dx.doi.org/10.7554/eLife.01916.024</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><th rowspan="2">Initial</th><th colspan="3">Final</th></tr><tr><th><italic>Normal</italic></th><th><italic>Semi-coiled</italic></th><th><italic>Curly-1</italic></th></tr></thead><tbody><tr><td><italic>Normal</italic></td><td/><td>0.28 ± 0.03 s<sup>−1</sup></td><td>0.08 ± 0.01 s<sup>−1</sup></td></tr><tr><td><italic>Semi-coiled</italic></td><td>1.6 ± 0.2 s<sup>−1</sup></td><td/><td>0.54 ± 0.08 s<sup>−1</sup></td></tr><tr><td><italic>Curly-1</italic></td><td>1.8 ± 0.2 s<sup>−1</sup></td><td>0.04 ± 0.02 s<sup>−1</sup></td><td/></tr></tbody></table><table-wrap-foot><fn><p>Transition rates between different flagellar waveforms: <italic>normal</italic> (CCW), <italic>semi-coiled</italic> and <italic>curly-1</italic> (both CW). Data from wild-type cells (<italic>N</italic> = 52 cells, 203 tumbles). Values are mean ± SEM.</p></fn></table-wrap-foot></table-wrap></p></sec></sec><sec id="s3" sec-type="discussion"><title>Discussion</title><p>The experimental approach described here allows the simultaneous, long-term observation of flagellar activity and swimming behavior in a single cell. By imaging many individual tumbles (<italic>N</italic> = 203 in wild-type cells), we are able to describe in great detail the underlying structure of a tumbling event. For instance, we can follow the sequence of flagellar waveforms that occurs when motors switch from CCW to CW rotation and back to CCW. <xref ref-type="fig" rid="fig4s2">Figure 4—figure supplement 2</xref> shows the distribution of possible sequences of flagellar waveforms during tumbles. In particular, the sequence of states from <italic>normal</italic> to <italic>semi-coiled</italic> to <italic>curly-1</italic> that we observed was described by <xref ref-type="bibr" rid="bib13">Darnton et al. (2007)</xref> as a ‘canonical tumble’. Although we cannot rule out that runs and tumbles in the optical trap are different in some ways than those in free swimming cells (<xref ref-type="bibr" rid="bib27">Min et al., 2009</xref>), our results are in qualitative agreement with these previous observations.</p><p>Our measurements also reveal the relationship between the cell’s run/tumble state and the CCW/CW rotational state of its flagella. In a multi-flagellated wild-type cell, a single CW flagellum (either <italic>semi-coiled</italic> or <italic>curly-1</italic>) is sufficient to induce a tumble, in agreement with the simple veto model (<xref ref-type="bibr" rid="bib38">Turner et al., 2000</xref>; <xref ref-type="bibr" rid="bib13">Darnton et al., 2007</xref>). However, the number of CW flagella during a tumble typically exceeds that expected from a cell with independently switching flagella (<xref ref-type="fig" rid="fig2">Figure 2C</xref>). The high fraction of CW flagella during tumbles in wild-type cells is in qualitative agreement with previous measurements by <xref ref-type="bibr" rid="bib38">Turner et al. (2000)</xref>, who observed that a majority of tumbles involved multiple flagella leaving the bundle. Our measurements using the CheY* strain provide an important piece of evidence linking inter-flagellar coupling to the chemotaxis network. We propose that fluctuations in the concentration of CheY-P are at the heart of wild-type <italic>E. coli</italic> behavior. In our theoretical analysis, the existence of temporal fluctuations was sufficient to explain all of our data. Stochastic simulations with and without CheY-P fluctuations (representing wild-type and CheY* cells, respectively) reproduced all of the observed differences between our two strains. <xref ref-type="fig" rid="fig4">Figure 4A</xref> summarizes how CheY-P fluctuations could lead to correlated flagellar switching. A well-known feature of the chemotaxis network is the sigmoidal relation between CW bias and CheY-P concentration (<xref ref-type="bibr" rid="bib10">Cluzel et al., 2000</xref>; <xref ref-type="bibr" rid="bib43">Yuan and Berg, 2013</xref>). A consequence of this non-linearity is that the probability of CW rotation is highly sensitive, and can respond dramatically to fluctuations in CheY-P levels, provided their amplitude is sufficiently large. As shown in <xref ref-type="fig" rid="fig4">Figure 4A</xref>, when CheY-P concentration is high, the cell experiences a near-100% probability of CW rotation and multiple motors switch CW at approximately the same time. In contrast, when CheY-P concentration is low, the probability of any motor rotating CW is essentially zero. This mechanism can explain the elevated number of CW flagella involved in tumbles (<xref ref-type="fig" rid="fig2">Figure 2C</xref>) and correlation between flagella states (<xref ref-type="fig" rid="fig2">Figure 2D</xref>). By contrast, in simulations where CheY-P level was held constant, flagellar switching was not as correlated and the majority of tumbles involved only a single CW flagellum (<xref ref-type="fig" rid="fig4">Figure 4B</xref>).</p><p>Despite the success of this model in reproducing our data, we must acknowledge that there is no direct experimental evidence to-date for the CheY-P fluctuations depicted in <xref ref-type="fig" rid="fig4">Figure 4A</xref>. Fluctuations in CheY-P have been inferred from experimental observations of CW bias in tethered-bead assays (<xref ref-type="bibr" rid="bib22">Korobkova et al., 2004</xref>). However, the fluctuations described in that study are different in their time scale and amplitude from what we found required to produce the observed correlations in flagellar rotational direction (‘Materials and methods’). Future investigation will be essential to resolving this issue and will likely have to involve direct measurements of CheY-P temporal dynamics in individual cells. Such measurements are challenging, but the development of intra-cellular fluorescence sensors for kinase activity in the network (<xref ref-type="bibr" rid="bib34">Sourjik and Berg, 2002</xref>) provides a promising approach.</p><p>As an alternative mechanism for the inter-flagellar correlations observed by <xref ref-type="bibr" rid="bib36">Terasawa et al. (2011)</xref> (and in the present work), <xref ref-type="bibr" rid="bib18">Hu and Tu (2013)</xref> proposed that hydrodynamic interactions between nearby flagella could also generate correlations in their rotational direction. One consequence of their model is that the flagellar switching rates in a cell with a single flagellum will be different than those in multi-flagellated cells. However, this prediction is not borne out by our data. In our experiments, the number of flagella per cell did not have a significant effect on the switching rates between CCW and CW states, nor on the switching rates between different flagellar waveforms (<xref ref-type="fig" rid="fig2s1 fig2s2">Figure 2—figure supplements 1 and 2</xref>). While we cannot rule out the presence of hydrodynamic interactions between flagella, these must satisfy the strict requirement that switching rates remain independent of flagellar number. In light of these constraints, we believe a mechanism in which chemotaxis network fluctuations engender inter-flagella correlations to be more plausible. For cells with <italic>N</italic><sub><italic>flag</italic></sub> &gt; 4, we note that both strains appear to deviate from the generalized veto model (<xref ref-type="fig" rid="fig3">Figure 3C</xref>). It is possible that hydrodynamic effects must be taken in consideration in cells with many flagella. Hydrodynamic coupling in the Hu model leads to a lower <italic>N</italic><sub><italic>eff</italic></sub>, in the direction of the deviation. Alternatively, a mechanism as described by <xref ref-type="bibr" rid="bib33">Sneddon et al. (2012)</xref>, in which cells with many flagella can run while some of its flagella rotate CW in the <italic>curly-1</italic> state, could lead to a similar deviation (<xref ref-type="fig" rid="fig1s5">Figure 1—figure supplement 5</xref>). Finally, we must consider the possibility that the apparent deviation is due to systematic experimental error, since it is more difficult to determine visually the state of each flagellum on cells with many flagella.</p><p>While a large number of studies have elucidated mathematical relationships between many of the components of the chemotaxis network (<xref ref-type="bibr" rid="bib6">Block et al., 1982</xref>; <xref ref-type="bibr" rid="bib10">Cluzel et al., 2000</xref>; <xref ref-type="bibr" rid="bib34">Sourjik and Berg, 2002</xref>; <xref ref-type="bibr" rid="bib32">Shimizu et al., 2010</xref>; <xref ref-type="bibr" rid="bib44">Yuan et al., 2012</xref>), there have been few experimental studies devoted to the relationship between individual flagella and whole-cell swimming (<xref ref-type="bibr" rid="bib38">Turner et al., 2000</xref>; <xref ref-type="bibr" rid="bib13">Darnton et al., 2007</xref>). As a result, existing models of bacterial chemotaxis have made drastically different assumptions in order to describe the swimming behavior of the whole cell (<xref ref-type="bibr" rid="bib8">Bray et al., 2007</xref>; <xref ref-type="bibr" rid="bib20">Jiang et al., 2010</xref>; <xref ref-type="bibr" rid="bib40">Vladimirov et al., 2010</xref>; <xref ref-type="bibr" rid="bib33">Sneddon et al., 2012</xref>). To the best of our knowledge, the current study provides the first experimentally-derived mathematical relation between flagellar and whole-cell-swimming states.</p><p>We propose that the details of this mapping are crucial for fully understanding bacterial chemotaxis. Results from recent theoretical models suggest that the details of flagellar mechanics can have significant effects on chemotactic drift. Turner et al. observed that, on average, the angular change in swimming direction upon tumbling increases as a function of the number of flagella that leave the bundle (<xref ref-type="bibr" rid="bib38">Turner et al., 2000</xref>). Vladimirov et al. showed that when this effect is incorporated into a theoretical model of bacterial chemotaxis, the chemotactic drift is nearly doubled (<xref ref-type="bibr" rid="bib40">Vladimirov et al., 2010</xref>).</p><p>Our observations that multi-flagellated wild-type cells tumble significantly less than expected also implies that the cell’s swimming behavior (and presumably its chemotactic response) is robust against variations in the number of flagella (<xref ref-type="fig" rid="fig4s4">Figure 4—figure supplement 4</xref>). We hypothesize that this phenomenon may confer evolutionary advantages, in light of the large fluctuations in flagellar numbers within a cell population (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1</xref>). If cells with many flagella did not behave like cells with fewer flagella, then they would spend the majority of their time tumbling, a behavior that would inhibit chemotaxis. <italic>E. coli</italic> thus appears to have developed a mechanism to achieve similar tumble biases with a wide range of flagellar number.</p></sec><sec id="s4" sec-type="materials|methods"><title>Materials and methods</title><sec id="s4-1"><title>Microbiology</title><sec id="s4-1-1"><title>Cell preparation</title><p>Experiments were performed using two <italic>E. coli</italic> strains (<xref ref-type="table" rid="tbl1">Table 1</xref>). The strain referred to as ‘wild-type’ is HCB1660 ([<xref ref-type="bibr" rid="bib39">Turner et al., 2010</xref>], gift of Howard Berg), <italic>ΔfliC</italic> expressing FliC<sup>S219C</sup> from a plasmid under the control of arabinose. The mutant protein FliC<sup>S219C</sup> was constructed to be specifically labeled with a maleimide functionalized fluorescent dye (<xref ref-type="bibr" rid="bib39">Turner et al., 2010</xref>). The strain referred to as ‘CheY*’ is PM87 (constructed for this study, see below), <italic>ΔfliC ΔcheBYZ</italic> expressing FliC<sup>S219C</sup> and CheY<sup>D13K</sup> from separate plasmids under the control of arabinose and IPTG, respectively. The mutant protein CheY<sup>D13K</sup> is constitutively active (<xref ref-type="bibr" rid="bib1">Alon et al., 1998</xref>), such that the CW bias was determined by the concentration of CheY<sup>D13K</sup>, decoupled from the chemotaxis network.</p><p>For each experiment, cells were picked from a single colony on an agar plate and grown overnight in 1 ml tryptone broth (1% [wt/vol] Bacto tryptone and 0.8% [wt/vol] NaCl) (<xref ref-type="bibr" rid="bib29">Saini et al., 2008</xref>; <xref ref-type="bibr" rid="bib27">Min et al., 2009</xref>) shaking at 265 RPM at 30°C with appropriate antibiotics. The overnight culture was diluted 100-fold into 12-ml tryptone broth and grown, shaking at 265 RPM at 30°C for 4.5 hr (to OD<sub>600</sub> ∼0.5) with appropriate inducers (500 µM arabinose and 50 µM IPTG). To visualize flagella, we used a fluorescence labeling protocol developed by <xref ref-type="bibr" rid="bib39">Turner et al. (2010)</xref>. The over-day culture was harvested, washed twice by slow centrifugation (1300×<italic>g</italic>, 10 min) and gently resuspended in 1 ml motility buffer (MB) (<xref ref-type="bibr" rid="bib13">Darnton et al., 2007</xref>) (10 mM KPO<sub>4</sub> (pH 7.0), 70 mM NaCl and 0.1 mM EDTA) and then in 0.5 ml MB. Flagella were specifically labeled using Alexa Fluor 532 C<sub>5</sub> Maleimide (A-10255; Life Technologies, Carlsbad, California). 1 mg of dry dye was dissolved in 300 µl H<sub>2</sub>0 by vortexing for 1 min. Aliquots containing 50 µl of dissolved dye were stored at −20°C. Cells in 500 µl MB were gently mixed with 5 µl of the dissolved dye and then incubated with slow rotation (∼10 RPM) at room temperature in the dark for 90 min. The labeled culture was washed and gently resuspended in 1 ml MB. Finally, cells were diluted 20-fold into 1 ml trap motility buffer (TMB) (70 mM NaCl, 100 mM Tris-Cl, 2% [wt/vol] glucose, and an oxygen-scavenging system [80 μg ml<sup>−1</sup> glucose oxidase and 13 μg ml<sup>−1</sup> catalase]) (<xref ref-type="bibr" rid="bib27">Min et al., 2009</xref>) and injected into the flow cell for trapping. At all stages, resuspension by pipetting was avoided to prevent shearing of the flagella (<xref ref-type="bibr" rid="bib38">Turner et al., 2000</xref>).</p></sec><sec id="s4-1-2"><title>Construction of strains</title><p>Bacterial strains and plasmids used in this study are listed in <xref ref-type="table" rid="tbl1">Table 1</xref>. Oligonucleotides used for generating mutations and creating plasmids are listed in <xref ref-type="table" rid="tbl2">Table 2</xref>. The generalized transducing phage P1<italic>vir</italic> was used in all transductional crosses (<xref ref-type="bibr" rid="bib37">Thomason et al., 2007</xref>). Chromosomal mutations were introduced using a standard λ Red recombination method as described by <xref ref-type="bibr" rid="bib14">Datsenko and Wanner (2000)</xref>. All primers were purchased from Integrated DNA Technologies (Coralville, Iowa), all sequencing was performed by ACGT (Wheeling, Illinois).</p><p>Strain PM87 (referred to as the CheY* strain in the text) was created from RP437 (wild-type for chemotaxis) in the following manner. First, <italic>cheBYZ</italic> was replaced by a chloramphenicol resistance cassette with flanking FRT sites using primers SK140F and SK140R along with pKD3 as a template (<xref ref-type="bibr" rid="bib14">Datsenko and Wanner, 2000</xref>). This deletion was then moved into a clean RP437 using P1 transduction (<xref ref-type="bibr" rid="bib37">Thomason et al., 2007</xref>), to create SK109 (<italic>cheBYZ::</italic>Cm). The chloramphenicol resistance cassette was removed from SK109 using Flp recombinase expressed from pCP20 (<xref ref-type="bibr" rid="bib9">Cherepanov and Wackernagel, 1995</xref>) to obtain strain SK110 (<italic>cheBYZ</italic>::FRT). Next, <italic>fliC</italic>::Tn5 from strain HCB1660 was moved into strain SK110 using P1 transduction to obtain strain SK112. Finally, strain PM87 was created by transforming strain SK112 with plasmids pPM5 and pMS164 to express FliC<sup>S219C</sup> and CheY<sup>D13K</sup>, respectively.</p><p>Standard molecular cloning techniques were used to construct plasmids (<xref ref-type="bibr" rid="bib31">Sambrook and Russell, 2001</xref>). Primers PM7F and PM7R were used to PCR amplify <italic>fliC</italic><sup><italic>S219C</italic></sup> and the P<sub>araBAD</sub> promoter from plasmid pBAD33-fliC<sup>S219C</sup> (gift of H Berg, [<xref ref-type="bibr" rid="bib39">Turner et al., 2010</xref>]) and to add restriction sites for AatII and SalI. The PCR product was then ligated into pZE11 (<xref ref-type="bibr" rid="bib24">Lutz and Bujard, 1997</xref>) after digesting both with restriction enzymes AatII and SalI, to obtain plasmid pPM5.</p></sec></sec><sec id="s4-2"><title>Instrument design</title><sec id="s4-2-1"><title>Optical traps</title><p>Experiments were performed using a dual optical trap instrument incorporating a custom flow cell and stroboscopic, epi-fluorescent imaging. The instrument design is shown in <xref ref-type="fig" rid="fig1s2">Figure 1—figure supplement 2</xref>. The optical trap was constructed following a previously described design (<xref ref-type="bibr" rid="bib12">Comstock et al., 2011</xref>). A 5-W, 1064-nm diode-pumped solid-state laser (BL-106C, Spectra-Physics, Santa Clara, California) was used to produce two optical traps via timesharing, by intermittently deflecting the laser with an acousto-optic deflector, AOD (DTD-274HD6, IntraAction, Bellwood, Illinois). The separation between the two traps was controlled by modulating the two deflection angles of the beams emanating from the acousto-optic deflector via the RF signal frequency driving the AOD. The IR beams were tightly focused to generate two optical traps by a 60×, water-immersion (1.2 NA) microscope objective (Nikon, Tokyo, Japan). An identical objective lens collected transmitted light for position detection and bright-field imaging, as described in <xref ref-type="bibr" rid="bib27">Min et al. (2009)</xref>. The flow chamber was positioned between the two objective lenses and was moveable relative to the two traps in all directions by a motorized three-axis translational stage (ESP301; Newport, Irvine, California). Cell motion was detected directly by the optical traps themselves, using back-focal plane interferometry, in which trap light scattered by an object relays the object’s position relative to the trap in all three directions (<xref ref-type="bibr" rid="bib17">Gittes and Schmidt, 1998</xref>). All devices and timing were controlled using custom LabVIEW (National Instruments, Austin, Texas) code (<xref ref-type="bibr" rid="bib12">Comstock et al., 2011</xref>).</p></sec><sec id="s4-2-2"><title>Fluorescence imaging</title><p>Epi-fluorescence imaging of trapped cells was achieved by excitation with a 532-nm laser (TECGL-30, World Star Tech, Toronto, Canada) and collecting backwards emitted photons with an EMCCD camera (iXon3 860 EMCCD, Andor, Belfast, Ireland). The beam size at the sample plane was approximately 20 microns in diameter. To obtain clear images of the flagella as they rotated at angular frequencies ∼100 Hz, it was necessary to take short exposures (<xref ref-type="fig" rid="fig1s3">Figure 1—figure supplement 3</xref>). To this end, we used stroboscopic illumination similarly to previous experiments (<xref ref-type="bibr" rid="bib13">Darnton et al., 2007</xref>). This was achieved by intermittently deflecting the 532-nm laser beam with a second acousto-optic modulator (AOM-802AF1, IntraAction), so as to generate short-duration excitation pulses (20 µs). The EMCCD recorded 128 × 128-pixel images (∼10 × 10 µm) synchronized with each excitation pulse (see timing in <xref ref-type="fig" rid="fig1s3">Figure 1—figure supplement 3</xref>, images <xref ref-type="fig" rid="fig1s4">Figure 1—figure supplement 4</xref>). Additionally, we pulsed the IR trapping laser out of phase with the fluorescence excitation at a rate of 16 kHz, a technique which has been shown to significantly reduce photo-bleaching, with minimal consequence to the trapping (<xref ref-type="bibr" rid="bib7">Brau et al., 2006</xref>; <xref ref-type="bibr" rid="bib12">Comstock et al., 2011</xref>) (<xref ref-type="fig" rid="fig1s3">Figure 1—figure supplement 3</xref>). The results were high speed movies showing sharp images of all the flagella on a cell as it runs and tumbles for many seconds (typically &gt;10 s), along with synchronous signal traces from the optical traps. Fluorescence movies were saved using Solis software (Andor), and subsequently analyzed manually using Matlab (Mathworks, Natick, Massachusetts). When played back in slow-motion, flagella could be reliably counted by eye during tumbles.</p></sec><sec id="s4-2-3"><title>Fluidics</title><p>Experiments were performed in a custom-built microfluidic chamber. See <xref ref-type="bibr" rid="bib28">Min et al. (2012)</xref> for a detailed description of the microfluidic chamber. Glass coverslips (12-545-M, 24 × 60−1, ThermoFisher, Waltham, Massachusetts) were sonicated in acetone for 5 min and rinsed with deionized water. The flow channel pattern was cut out from Nescofilm (Karlan, Phoenix, Arizona) using a laser engraver (Versa Laser, Scottsdale, Arizona) and placed between two coverslips, one of which had custom-drilled holes (0.05-inch diameter) for inlets and outlets. The Nescofilm flow channel pattern was bonded to coverslips by melting on a hot plate for 4 min. The completed flow chamber was inserted into a custom metal frame where inlet and outlet tubing (ABW00001; Tygon, Saint-Gobain, Paris, France) were screwed on for a tight seal. The two channels of the flow chamber were continuously injected with appropriate buffers using a syringe pump (PHD2000; Harvard Apparatus, Holliston, Massachusetts) at a linear speed of 30 μm/s, which is approximately equal to the swimming speed of a healthy cell. The upper channel was injected with TMB (see above), while the bottom channel contained cells in TMB. Cells were trapped in the bottom channel and then moved into the upper channel for observation by displacing the motorized flow chamber relative to the traps.</p></sec></sec><sec id="s4-3"><title>Electron microscopy</title><p>TEM images were recorded using the JEOL 2100 cryo-Transmission Electron Microscope (TEM) at the Frederick Seitz Materials Research Laboratory Central Facilities at the University of Illinois at Urbana–Champaign, following the protocol of <xref ref-type="bibr" rid="bib30">Saini et al. (2010)</xref>. Briefly, cells were grown as described above, and used without fluorescent labeling. Cells were fixed with glutaraldehyde and then placed on 200 Mesh Carbon Coated Copper grids (Cat. # 182; Canemco, Lakefield, Canada), which were used as supports for sample loading and imaging. Images were taken at 200 kV with a camera exposure lasting 1 s. Finally, images were contrast adjusted and an image dilation was performed in Matlab to make flagella more visible. The distribution of flagella per cell is shown in <xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1</xref>.</p></sec><sec id="s4-4"><title>Data analysis</title><sec id="s4-4-1"><title>Wavelet analysis for run-tumble detection</title><p>Determination of runs and tumbles from the optical trap signal was done as described previously (<xref ref-type="bibr" rid="bib27">Min et al., 2009</xref>). See example trace in <xref ref-type="fig" rid="fig1">Figure 1C</xref>.</p></sec><sec id="s4-4-2"><title>Image analysis</title><p>Images were contrast adjusted to make flagella easier to see using Matlab. Movies were then manually analyzed by eye to count flagella and to identify flagella waveforms. The state of each flagellum was identified during each 100-ms time window (10 movie frames at 100 f.p.s.). See example cells in <xref ref-type="fig" rid="fig1">Figure 1C</xref> and <xref ref-type="fig" rid="fig1s4">Figure 1—figure supplement 4</xref>.</p></sec></sec><sec id="s4-5"><title>Theoretical modeling</title><sec id="s4-5-1"><title>The veto model</title><p>The veto model assumes that all flagella must be rotating counter-clockwise (CCW) for the cell to run. Any flagellum can ‘veto’ the run by rotating clockwise (CW), which causes the cell to tumble (<xref ref-type="fig" rid="fig1">Figure 1A</xref>). To describe the veto model mathematically, we write the expression for the run bias (the fraction of time that a cell spends running) as a function of the CCW bias (the fraction of time that each flagellum spends rotating CCW). Assuming that the CCW bias is fixed in time, one obtains:<disp-formula id="equ4"><label>(S1)</label><mml:math id="m4"><mml:mrow><mml:mi mathvariant="italic">Run</mml:mi><mml:mo> </mml:mo><mml:mi mathvariant="italic">bias</mml:mi><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi mathvariant="italic">Time</mml:mi><mml:mo> </mml:mo><mml:mi mathvariant="italic">ALL</mml:mi><mml:mo> </mml:mo><mml:mi mathvariant="italic">motors</mml:mi><mml:mo> </mml:mo><mml:mi mathvariant="italic">CCW</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="italic">Total</mml:mi><mml:mo> </mml:mo><mml:mi mathvariant="italic">time</mml:mi></mml:mrow></mml:mfrac></mml:mrow></mml:math></disp-formula><disp-formula id="equ5"><label>(S2)</label><mml:math id="m5"><mml:mrow><mml:mi mathvariant="italic">Run</mml:mi><mml:mo> </mml:mo><mml:mi mathvariant="italic">bias</mml:mi><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mi>C</mml:mi><mml:mi>C</mml:mi><mml:mi>W</mml:mi><mml:mo> </mml:mo><mml:mi mathvariant="italic">bias</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mrow><mml:mi mathvariant="italic">flag</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msup></mml:mrow></mml:math></disp-formula>where <italic>N</italic><sub><italic>flag</italic></sub> is the number of flagella. We can write the tumble bias (TB) as a function of the run bias, and the CW bias (CB) as a function of CCW bias: <disp-formula id="equ6"><label>(S3)</label><mml:math id="m6"><mml:mrow><mml:mi>T</mml:mi><mml:mi>B</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi mathvariant="italic">Run</mml:mi><mml:mo> </mml:mo><mml:mi mathvariant="italic">Bias</mml:mi></mml:mrow></mml:math></disp-formula><disp-formula id="equ7"><label>(S4)</label><mml:math id="m7"><mml:mrow><mml:mi>C</mml:mi><mml:mi>B</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>C</mml:mi><mml:mi>C</mml:mi><mml:mi>W</mml:mi><mml:mo> </mml:mo><mml:mi mathvariant="italic">Bias</mml:mi></mml:mrow></mml:math></disp-formula></p><p>We then use these relations to solve for the tumble bias as a function of CW bias and the number of flagella:<disp-formula id="equ8"><label>(S5)</label><mml:math id="m8"><mml:mrow><mml:mi>T</mml:mi><mml:mi>B</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>C</mml:mi><mml:mi>B</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi mathvariant="italic">flag</mml:mi></mml:msub></mml:mrow></mml:msup></mml:mrow></mml:math></disp-formula></p></sec><sec id="s4-5-2"><title>Simulating the effect of CheY-P fluctuations</title><p>Stochastic simulations were performed to model the effects of fluctuations in CheY-P concentration over time, using a method similar to <xref ref-type="bibr" rid="bib33">Sneddon, et al. (2012)</xref>. Simulations are illustrated in <xref ref-type="fig" rid="fig4">Figure 4A</xref>. CheY-P fluctuations were generated using <xref ref-type="disp-formula" rid="equ9">Equation S6</xref> below (see <xref ref-type="table" rid="tbl3">Table 3</xref> for parameter values), where <italic>Y</italic> is the instantaneous concentration of CheY-P, <italic>τ</italic> is the characteristic timescale of fluctuations in [CheY-P], <italic>σ</italic><sup>2</sup> is the variance in [CheY-P], and <italic>ξ</italic> are normally distributed random numbers with unit variance. The time resolution of all simulations was one data point per millisecond.<disp-formula id="equ9"><label>(S6)</label><mml:math id="m9"><mml:mrow><mml:mi>Y</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:mi>d</mml:mi><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>Y</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>−</mml:mo><mml:mfrac><mml:mrow><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>Y</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>−</mml:mo><mml:mo>(</mml:mo><mml:mi>Y</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mrow><mml:mi>τ</mml:mi></mml:mfrac><mml:mi>d</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:mi>σ</mml:mi><mml:msqrt><mml:mrow><mml:mfrac><mml:mrow><mml:mn>2</mml:mn><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow><mml:mi>τ</mml:mi></mml:mfrac></mml:mrow></mml:msqrt><mml:mi>ξ</mml:mi></mml:mrow></mml:math></disp-formula></p><p>The CW bias was calculated from [CheY-P] using the Hill function relationship determined by <xref ref-type="bibr" rid="bib10">Cluzel et al. (2000)</xref>. After generating the CheY-P time trace, the mean CW bias of that trace was checked, to ensure that the bias was within the same range as the data mean (0.13–0.145). The rotational state of each flagellar motor was then determined stochastically, using <xref ref-type="disp-formula" rid="equ1">Equation 1</xref> in Sneddon et al., the CW bias, and characteristic motor switching rate (<italic>ω</italic>) to set the CW/CCW transition rates ([<xref ref-type="bibr" rid="bib33">Sneddon et al., 2012</xref>] see <xref ref-type="table" rid="tbl3">Table 3</xref> for parameter values).</p><p>Runs and tumbles were determined using the veto model: whenever any flagellum was rotating CW the cell tumbled, when all flagella were rotating CCW the cell ran. Finally, the simulated data were analyzed in the exact same manner as the experimental data, to extract the mean number of flagella participating in tumbles (<xref ref-type="fig" rid="fig2">Figure 2C</xref>), flagella cross-correlation (<xref ref-type="fig" rid="fig2">Figure 2D</xref>), the effective number of flagella (<italic>N</italic><sub><italic>eff</italic></sub>, <xref ref-type="fig" rid="fig3">Figure 3C</xref>) and the deviation from veto model <italic>η</italic>, (<xref ref-type="fig" rid="fig4">Figure 4C</xref>).</p><p>The values of <italic>σ</italic> and <italic>τ</italic>, which denote the amplitude and timescale of simulated CheY-P fluctuations, were determined by scanning through a range of values, and minimizing the global, reduced χ<sup>2</sup> (<xref ref-type="bibr" rid="bib5">Bevington and Robinson, 2003</xref>) from comparisons of simulations to data. The minimum reduced χ<sup>2</sup> value indicates that all of the experimental data is best reproduced by simulations in which <italic>σ</italic> = 1.0 µM<sup>2</sup> and <italic>τ</italic> = 0.2 s<sup>−1</sup>. The global reduced χ<sup>2</sup> was calculated by summing the individual reduced χ<sup>2</sup> values from fits to data in <xref ref-type="fig" rid="fig2 fig4">Figures 2C,D and 4C</xref>.</p></sec><sec id="s4-5-3"><title>Incorporating the effect of <italic>curly-1</italic> flagella during runs</title><p>To test the effect of runs involving <italic>curly-1</italic> flagella, we augmented our simulations following <xref ref-type="bibr" rid="bib33">Sneddon et al. (2012)</xref>. For each simulated flagellum, the binary CCW/CW trace was converted into a 3-state (<italic>normal</italic>, <italic>semi-coiled</italic>, <italic>curly-1</italic>) trace, using <italic>λ</italic><sup>−1</sup> = 0.68 s<sup>−1</sup> as the transition rate from <italic>semi-coiled</italic> to <italic>curly-1</italic> (<xref ref-type="table" rid="tbl3">Table 3</xref>). The veto model was then applied, with the caveat that cells ran 18% of the time when there was a single <italic>curly-1</italic> flagellum, provided the other flagella were all CCW (<italic>normal</italic>) (<xref ref-type="fig" rid="fig1s5">Figure 1—figure supplement 5</xref>).</p></sec></sec></sec></body><back><ack id="ack"><title>Acknowledgements</title><p>We are grateful to Howard Berg, Linda Turner, Vedhavalli Nathan and Philippe Cluzel for advice and for providing reagents. We thank current and former members of the Chemla, Golding, and Rao Laboratories for providing help with experiments.</p></ack><sec sec-type="additional-information"><title>Additional information</title><fn-group content-type="competing-interest"><title>Competing interests</title><fn fn-type="conflict" id="conf1"><p>The author declares that no competing interests exist.</p></fn></fn-group><fn-group content-type="author-contribution"><title>Author contributions</title><fn fn-type="con" id="con1"><p>PJM, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents</p></fn><fn fn-type="con" id="con2"><p>SK, Conception and design, Drafting or revising the article, Contributed unpublished essential data or reagents</p></fn><fn fn-type="con" id="con3"><p>CVR, Conception and design, Drafting or revising the article, Contributed unpublished essential data or reagents</p></fn><fn fn-type="con" id="con4"><p>IG, Conception and design, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con5"><p>YRC, Conception and design, Analysis and interpretation of data, Drafting or revising the article</p></fn></fn-group></sec><ref-list><title>References</title><ref id="bib1"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Alon</surname><given-names>U</given-names></name><name><surname>Camarena</surname><given-names>L</given-names></name><name><surname>Surette</surname><given-names>MG</given-names></name><name><surname>Aguera y Arcas</surname><given-names>B</given-names></name><name><surname>Liu</surname><given-names>Y</given-names></name><name><surname>Leibler</surname><given-names>S</given-names></name><name><surname>Stock</surname><given-names>JB</given-names></name></person-group><year>1998</year><article-title>Response regulator output in bacterial chemotaxis</article-title><source>The EMBO Journal</source><volume>17</volume><fpage>4238</fpage><lpage>4248</lpage><pub-id pub-id-type="doi">10.1093/emboj/17.15.4238</pub-id></element-citation></ref><ref id="bib2"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Andrews</surname><given-names>BW</given-names></name><name><surname>Yi</surname><given-names>TM</given-names></name><name><surname>Iglesias</surname><given-names>PA</given-names></name></person-group><year>2006</year><article-title>Optimal noise filtering in the chemotactic response of <italic>Escherichia coli</italic></article-title><source>PLOS Computational Biology</source><volume>2</volume><fpage>e154</fpage><pub-id pub-id-type="doi">10.1371/journal.pcbi.0020154</pub-id></element-citation></ref><ref id="bib3"><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Berg</surname><given-names>HC</given-names></name></person-group><year>2004</year><source>E. coli in motion</source><publisher-loc>New York</publisher-loc><publisher-name>Springer</publisher-name></element-citation></ref><ref id="bib4"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Berg</surname><given-names>HC</given-names></name><name><surname>Brown</surname><given-names>DA</given-names></name></person-group><year>1972</year><article-title>Chemotaxis in <italic>Escherichia coli</italic> analysed by three-dimensional tracking</article-title><source>Nature</source><volume>239</volume><fpage>500</fpage><lpage>504</lpage><pub-id pub-id-type="doi">10.1038/239500a0</pub-id></element-citation></ref><ref id="bib5"><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Bevington</surname><given-names>PR</given-names></name><name><surname>Robinson</surname><given-names>DK</given-names></name></person-group><year>2003</year><source>Data reduction and error analysis for the physical sciences</source><edition>3rd edition</edition><publisher-loc>Boston</publisher-loc><publisher-name>McGraw-Hill</publisher-name></element-citation></ref><ref id="bib6"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Block</surname><given-names>SM</given-names></name><name><surname>Segall</surname><given-names>JE</given-names></name><name><surname>Berg</surname><given-names>HC</given-names></name></person-group><year>1982</year><article-title>Impulse responses in bacterial chemotaxis</article-title><source>Cell</source><volume>31</volume><fpage>215</fpage><lpage>226</lpage><pub-id pub-id-type="doi">10.1016/0092-8674(82)90421-4</pub-id></element-citation></ref><ref id="bib7"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Brau</surname><given-names>RR</given-names></name><name><surname>Tarsa</surname><given-names>PB</given-names></name><name><surname>Ferrer</surname><given-names>JM</given-names></name><name><surname>Lee</surname><given-names>P</given-names></name><name><surname>Lang</surname><given-names>MJ</given-names></name></person-group><year>2006</year><article-title>Interlaced optical force-fluorescence measurements for single molecule biophysics</article-title><source>Biophysical Journal</source><volume>91</volume><fpage>1069</fpage><lpage>1077</lpage><pub-id pub-id-type="doi">10.1529/biophysj.106.082602</pub-id></element-citation></ref><ref id="bib8"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bray</surname><given-names>D</given-names></name><name><surname>Levin</surname><given-names>MD</given-names></name><name><surname>Lipkow</surname><given-names>K</given-names></name></person-group><year>2007</year><article-title>The chemotactic behavior of computer-based surrogate bacteria</article-title><source>Current Biology</source><volume>17</volume><fpage>12</fpage><lpage>19</lpage><pub-id pub-id-type="doi">10.1016/j.cub.2006.11.027</pub-id></element-citation></ref><ref id="bib9"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cherepanov</surname><given-names>PP</given-names></name><name><surname>Wackernagel</surname><given-names>W</given-names></name></person-group><year>1995</year><article-title>Gene disruption in <italic>Escherichia coli</italic>: TcR and KmR cassettes with the option of Flp-catalyzed excision of the antibiotic-resistance determinant</article-title><source>Gene</source><volume>158</volume><fpage>9</fpage><lpage>14</lpage><pub-id pub-id-type="doi">10.1016/0378-1119(95)00193-A</pub-id></element-citation></ref><ref id="bib10"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cluzel</surname><given-names>P</given-names></name><name><surname>Surette</surname><given-names>M</given-names></name><name><surname>Leibler</surname><given-names>S</given-names></name></person-group><year>2000</year><article-title>An ultrasensitive bacterial motor revealed by monitoring signaling proteins in single cells</article-title><source>Science</source><volume>287</volume><fpage>1652</fpage><lpage>1655</lpage><pub-id pub-id-type="doi">10.1126/science.287.5458.1652</pub-id></element-citation></ref><ref id="bib11"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cohen-Ben-Lulu</surname><given-names>GN</given-names></name><name><surname>Francis</surname><given-names>NR</given-names></name><name><surname>Shimoni</surname><given-names>E</given-names></name><name><surname>Noy</surname><given-names>D</given-names></name><name><surname>Davidov</surname><given-names>Y</given-names></name><name><surname>Prasad</surname><given-names>K</given-names></name><name><surname>Sagi</surname><given-names>Y</given-names></name><name><surname>Cecchini</surname><given-names>G</given-names></name><name><surname>Johnstone</surname><given-names>RM</given-names></name><name><surname>Eisenbach</surname><given-names>M</given-names></name></person-group><year>2008</year><article-title>The bacterial flagellar switch complex is getting more complex</article-title><source>The EMBO Journal</source><volume>27</volume><fpage>1134</fpage><lpage>1144</lpage><pub-id pub-id-type="doi">10.1038/emboj.2008.48</pub-id></element-citation></ref><ref id="bib12"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Comstock</surname><given-names>MJ</given-names></name><name><surname>Ha</surname><given-names>T</given-names></name><name><surname>Chemla</surname><given-names>YR</given-names></name></person-group><year>2011</year><article-title>Ultrahigh-resolution optical trap with single-fluorophore sensitivity</article-title><source>Nature Methods</source><volume>8</volume><fpage>335</fpage><lpage>340</lpage><pub-id pub-id-type="doi">10.1038/nmeth.1574</pub-id></element-citation></ref><ref id="bib13"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Darnton</surname><given-names>NC</given-names></name><name><surname>Turner</surname><given-names>L</given-names></name><name><surname>Rojevsky</surname><given-names>S</given-names></name><name><surname>Berg</surname><given-names>HC</given-names></name></person-group><year>2007</year><article-title>On torque and tumbling in swimming <italic>Escherichia coli</italic></article-title><source>Journal of Bacteriology</source><volume>189</volume><fpage>1756</fpage><lpage>1764</lpage><pub-id pub-id-type="doi">10.1128/Jb.01501-06</pub-id></element-citation></ref><ref id="bib14"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Datsenko</surname><given-names>KA</given-names></name><name><surname>Wanner</surname><given-names>BL</given-names></name></person-group><year>2000</year><article-title>One-step inactivation of chromosomal genes in <italic>Escherichia coli</italic> K-12 using PCR products</article-title><source>Proceedings of the National Academy of Sciences of the United States of America</source><volume>97</volume><fpage>6640</fpage><lpage>6645</lpage><pub-id pub-id-type="doi">10.1073/pnas.120163297</pub-id></element-citation></ref><ref id="bib15"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Emonet</surname><given-names>T</given-names></name><name><surname>Macal</surname><given-names>CM</given-names></name><name><surname>North</surname><given-names>MJ</given-names></name><name><surname>Wickersham</surname><given-names>CE</given-names></name><name><surname>Cluzel</surname><given-names>P</given-names></name></person-group><year>2005</year><article-title>AgentCell: a digital single-cell assay for bacterial chemotaxis</article-title><source>Bioinformatics</source><volume>21</volume><fpage>2714</fpage><lpage>2721</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/bti391</pub-id></element-citation></ref><ref id="bib16"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Flores</surname><given-names>M</given-names></name><name><surname>Shimizu</surname><given-names>TS</given-names></name><name><surname>ten Wolde</surname><given-names>PR</given-names></name><name><surname>Tostevin</surname><given-names>F</given-names></name></person-group><year>2012</year><article-title>Signaling noise enhances chemotactic drift of <italic>E. coli</italic></article-title><source>Physical Review Letters</source><volume>109</volume><fpage>148101</fpage><pub-id pub-id-type="doi">10.1103/PhysRevLett.109.148101</pub-id></element-citation></ref><ref id="bib17"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gittes</surname><given-names>F</given-names></name><name><surname>Schmidt</surname><given-names>CF</given-names></name></person-group><year>1998</year><article-title>Interference model for back-focal-plane displacement detection in optical tweezers</article-title><source>Optics Letters</source><volume>23</volume><fpage>7</fpage><lpage>9</lpage><pub-id pub-id-type="doi">10.1364/OL.23.000007</pub-id></element-citation></ref><ref id="bib18"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hu</surname><given-names>B</given-names></name><name><surname>Tu</surname><given-names>YH</given-names></name></person-group><year>2013</year><article-title>Coordinated switching of bacterial flagellar motors: evidence for direct motor-motor coupling?</article-title><source>Physical Review Letters</source><volume>110</volume><fpage>158703</fpage><pub-id pub-id-type="doi">10.1103/Physrevlett.110.158703</pub-id></element-citation></ref><ref id="bib19"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ishihara</surname><given-names>A</given-names></name><name><surname>Segall</surname><given-names>JE</given-names></name><name><surname>Block</surname><given-names>SM</given-names></name><name><surname>Berg</surname><given-names>HC</given-names></name></person-group><year>1983</year><article-title>Coordination of flagella on filamentous cells of <italic>Escherichia coli</italic></article-title><source>Journal of Bacteriology</source><volume>155</volume><fpage>228</fpage><lpage>237</lpage></element-citation></ref><ref id="bib20"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jiang</surname><given-names>LL</given-names></name><name><surname>Qi</surname><given-names>OY</given-names></name><name><surname>Tu</surname><given-names>YH</given-names></name></person-group><year>2010</year><article-title>Quantitative odeling of <italic>Escherichia coli</italic> chemotactic motion in environments varying in space and time</article-title><source>PLOS Computational Biology</source><volume>6</volume><pub-id pub-id-type="doi">10.1371/journal.pcbi.1000735</pub-id></element-citation></ref><ref id="bib21"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kalinin</surname><given-names>YV</given-names></name><name><surname>Jiang</surname><given-names>L</given-names></name><name><surname>Tu</surname><given-names>Y</given-names></name><name><surname>Wu</surname><given-names>M</given-names></name></person-group><year>2009</year><article-title>Logarithmic sensing in <italic>Escherichia coli</italic> bacterial chemotaxis</article-title><source>Biophysical Journal</source><volume>96</volume><fpage>2439</fpage><lpage>2448</lpage><pub-id pub-id-type="doi">10.1016/j.bpj.2008.10.027</pub-id></element-citation></ref><ref id="bib22"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Korobkova</surname><given-names>E</given-names></name><name><surname>Emonet</surname><given-names>T</given-names></name><name><surname>Vilar</surname><given-names>JMG</given-names></name><name><surname>Shimizu</surname><given-names>TS</given-names></name><name><surname>Cluzel</surname><given-names>P</given-names></name></person-group><year>2004</year><article-title>From molecular noise to behavioural variability in a single bacterium</article-title><source>Nature</source><volume>428</volume><fpage>574</fpage><lpage>578</lpage><pub-id pub-id-type="doi">10.1038/Nature02404</pub-id></element-citation></ref><ref id="bib23"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Korobkova</surname><given-names>EA</given-names></name><name><surname>Emonet</surname><given-names>T</given-names></name><name><surname>Park</surname><given-names>H</given-names></name><name><surname>Cluzel</surname><given-names>P</given-names></name></person-group><year>2006</year><article-title>Hidden stochastic nature of a single bacterial motor</article-title><source>Physical Review Letters</source><volume>96</volume><fpage>058105</fpage><pub-id pub-id-type="doi">10.1103/PhysRevLett.96.058105</pub-id></element-citation></ref><ref id="bib24"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lutz</surname><given-names>R</given-names></name><name><surname>Bujard</surname><given-names>H</given-names></name></person-group><year>1997</year><article-title>Independent and tight regulation of transcriptional units in <italic>Escherichia coli</italic> via the LacR/O, the TetR/O and AraC/I1-I2 regulatory elements</article-title><source>Nucleic Acids Research</source><volume>25</volume><fpage>1203</fpage><lpage>1210</lpage><pub-id pub-id-type="doi">10.1093/nar/25.6.1203</pub-id></element-citation></ref><ref id="bib25"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Macnab</surname><given-names>RM</given-names></name><name><surname>Ornston</surname><given-names>MK</given-names></name></person-group><year>1977</year><article-title>Normal-to-curly flagellar transitions and their role in bacterial tumbling. Stabilization of an alternative quaternary structure by mechanical force</article-title><source>Journal of Molecular Biology</source><volume>112</volume><fpage>1</fpage><lpage>30</lpage><pub-id pub-id-type="doi">10.1016/S0022-2836(77)80153-8</pub-id></element-citation></ref><ref id="bib26"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Matthaus</surname><given-names>F</given-names></name><name><surname>Jagodic</surname><given-names>M</given-names></name><name><surname>Dobnikar</surname><given-names>J</given-names></name></person-group><year>2009</year><article-title><italic>E. coli</italic> superdiffusion and chemotaxis-search strategy, precision, and motility</article-title><source>Biophysical Journal</source><volume>97</volume><fpage>946</fpage><lpage>957</lpage><pub-id pub-id-type="doi">10.1016/j.bpj.2009.04.065</pub-id></element-citation></ref><ref id="bib27"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Min</surname><given-names>TL</given-names></name><name><surname>Mears</surname><given-names>PJ</given-names></name><name><surname>Chubiz</surname><given-names>LM</given-names></name><name><surname>Rao</surname><given-names>CV</given-names></name><name><surname>Golding</surname><given-names>I</given-names></name><name><surname>Chemla</surname><given-names>YR</given-names></name></person-group><year>2009</year><article-title>High-resolution, long-term characterization of bacterial motility using optical tweezers</article-title><source>Nature Methods</source><volume>6</volume><fpage>831</fpage><lpage>835</lpage><pub-id pub-id-type="doi">10.1038/nmeth.1380</pub-id></element-citation></ref><ref id="bib28"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Min</surname><given-names>TL</given-names></name><name><surname>Mears</surname><given-names>PJ</given-names></name><name><surname>Golding</surname><given-names>I</given-names></name><name><surname>Chemla</surname><given-names>YR</given-names></name></person-group><year>2012</year><article-title>Chemotactic adaptation kinetics of individual <italic>Escherichia coli</italic> cells</article-title><source>Proceedings of the National Academy of Sciences of the United States of America</source><volume>109</volume><fpage>9869</fpage><lpage>9874</lpage><pub-id pub-id-type="doi">10.1073/pnas.1120218109</pub-id></element-citation></ref><ref id="bib45"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Parkinson</surname><given-names>JS</given-names></name><name><surname>Houts</surname><given-names>SE</given-names></name></person-group><year>1982</year><article-title>Isolation and behavior of Escherichia coli deletion mutants lacking chemotaxis functions</article-title><source>Journal of Bacteriology</source><volume>151</volume><fpage>106</fpage><lpage>113</lpage></element-citation></ref><ref id="bib29"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Saini</surname><given-names>S</given-names></name><name><surname>Brown</surname><given-names>JD</given-names></name><name><surname>Aldridge</surname><given-names>PD</given-names></name><name><surname>Rao</surname><given-names>CV</given-names></name></person-group><year>2008</year><article-title>FliZ is a posttranslational activator of FlhD4C2-dependent flagellar gene expression</article-title><source>Journal of Bacteriology</source><volume>190</volume><fpage>4979</fpage><lpage>4988</lpage><pub-id pub-id-type="doi">10.1128/JB.01996-07</pub-id></element-citation></ref><ref id="bib30"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Saini</surname><given-names>S</given-names></name><name><surname>Koirala</surname><given-names>S</given-names></name><name><surname>Floess</surname><given-names>E</given-names></name><name><surname>Mears</surname><given-names>PJ</given-names></name><name><surname>Chemla</surname><given-names>YR</given-names></name><name><surname>Golding</surname><given-names>I</given-names></name><name><surname>Aldridge</surname><given-names>C</given-names></name><name><surname>Aldridge</surname><given-names>PD</given-names></name><name><surname>Rao</surname><given-names>CV</given-names></name></person-group><year>2010</year><article-title>FliZ induces a Kinetic switch in flagellar gene expression</article-title><source>Journal of Bacteriology</source><volume>192</volume><fpage>6477</fpage><lpage>6481</lpage><pub-id pub-id-type="doi">10.1128/Jb.00751-10</pub-id></element-citation></ref><ref id="bib31"><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Sambrook</surname><given-names>J</given-names></name><name><surname>Russell</surname><given-names>DW</given-names></name></person-group><year>2001</year><source>Molecular cloning: a laboratory manual</source><edition>3rd edition</edition><publisher-loc>Cold Spring Harbor, NY</publisher-loc><publisher-name>Cold Spring Harbor Laboratory Press</publisher-name></element-citation></ref><ref id="bib32"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shimizu</surname><given-names>TS</given-names></name><name><surname>Tu</surname><given-names>Y</given-names></name><name><surname>Berg</surname><given-names>HC</given-names></name></person-group><year>2010</year><article-title>A modular gradient-sensing network for chemotaxis in <italic>Escherichia coli</italic> revealed by responses to time-varying stimuli</article-title><source>Molecular Systems Biology</source><volume>6</volume><fpage>382</fpage><pub-id pub-id-type="doi">10.1038/msb.2010.37</pub-id></element-citation></ref><ref id="bib33"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sneddon</surname><given-names>MW</given-names></name><name><surname>Pontius</surname><given-names>W</given-names></name><name><surname>Emonet</surname><given-names>T</given-names></name></person-group><year>2012</year><article-title>Stochastic coordination of multiple actuators reduces latency and improves chemotactic response in bacteria</article-title><source>Proceedings of the National Academy of Sciences of the United States of America</source><volume>109</volume><fpage>805</fpage><lpage>810</lpage><pub-id pub-id-type="doi">10.1073/pnas.1113706109</pub-id></element-citation></ref><ref id="bib34"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sourjik</surname><given-names>V</given-names></name><name><surname>Berg</surname><given-names>HC</given-names></name></person-group><year>2002</year><article-title>Receptor sensitivity in bacterial chemotaxis</article-title><source>Proceedings of the National Academy of Sciences of the United States of America</source><volume>99</volume><fpage>123</fpage><lpage>127</lpage><pub-id pub-id-type="doi">10.1073/pnas.011589998</pub-id></element-citation></ref><ref id="bib35"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Spiro</surname><given-names>PA</given-names></name><name><surname>Parkinson</surname><given-names>JS</given-names></name><name><surname>Othmer</surname><given-names>HG</given-names></name></person-group><year>1997</year><article-title>A model of excitation and adaptation in bacterial chemotaxis</article-title><source>Proceedings of the National Academy of Sciences of the United States of America</source><volume>94</volume><fpage>7263</fpage><lpage>7268</lpage><pub-id pub-id-type="doi">10.1073/pnas.94.14.7263</pub-id></element-citation></ref><ref id="bib36"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Terasawa</surname><given-names>S</given-names></name><name><surname>Fukuoka</surname><given-names>H</given-names></name><name><surname>Inoue</surname><given-names>Y</given-names></name><name><surname>Sagawa</surname><given-names>T</given-names></name><name><surname>Takahashi</surname><given-names>H</given-names></name><name><surname>Ishijima</surname><given-names>A</given-names></name></person-group><year>2011</year><article-title>Coordinated reversal of flagellar motors on a single <italic>Escherichia coli</italic> cell</article-title><source>Biophysical Journal</source><volume>100</volume><fpage>2193</fpage><lpage>2200</lpage><pub-id pub-id-type="doi">10.1016/j.bpj.2011.03.030</pub-id></element-citation></ref><ref id="bib37"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Thomason</surname><given-names>L</given-names></name><name><surname>Costantino</surname><given-names>N</given-names></name><name><surname>Court</surname><given-names>D</given-names></name></person-group><year>2007</year><article-title><italic>E. coli</italic> genome manipulation by P1 transduction</article-title><source>Current Protocols in Molecular Biology</source><comment>Chapter 1:Unit 1.17</comment><pub-id pub-id-type="doi">10.1002/0471142727.mb0117s79</pub-id></element-citation></ref><ref id="bib38"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Turner</surname><given-names>L</given-names></name><name><surname>Ryu</surname><given-names>WS</given-names></name><name><surname>Berg</surname><given-names>HC</given-names></name></person-group><year>2000</year><article-title>Real-time imaging of fluorescent flagellar filaments</article-title><source>Journal of Bacteriology</source><volume>182</volume><fpage>2793</fpage><lpage>2801</lpage><pub-id pub-id-type="doi">10.1128/JB.182.10.2793-2801.2000</pub-id></element-citation></ref><ref id="bib39"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Turner</surname><given-names>L</given-names></name><name><surname>Zhang</surname><given-names>R</given-names></name><name><surname>Darnton</surname><given-names>NC</given-names></name><name><surname>Berg</surname><given-names>HC</given-names></name></person-group><year>2010</year><article-title>Visualization of flagella during bacterial Swarming</article-title><source>Journal of Bacteriology</source><volume>192</volume><fpage>3259</fpage><lpage>3267</lpage><pub-id pub-id-type="doi">10.1128/JB.00083-10</pub-id></element-citation></ref><ref id="bib40"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vladimirov</surname><given-names>N</given-names></name><name><surname>Lebiedz</surname><given-names>D</given-names></name><name><surname>Sourjik</surname><given-names>V</given-names></name></person-group><year>2010</year><article-title>Predicted auxiliary navigation mechanism of peritrichously flagellated chemotactic bacteria</article-title><source>PLOS Computational Biology</source><volume>6</volume><pub-id pub-id-type="doi">10.1371/journal.pcbi.1000717</pub-id></element-citation></ref><ref id="bib41"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vladimirov</surname><given-names>N</given-names></name><name><surname>Lovdok</surname><given-names>L</given-names></name><name><surname>Lebiedz</surname><given-names>D</given-names></name><name><surname>Sourjik</surname><given-names>V</given-names></name></person-group><year>2008</year><article-title>Dependence of bacterial chemotaxis on gradient shape and adaptation rate</article-title><source>PLoS Computational Biology</source><volume>4</volume><fpage>e1000242</fpage><pub-id pub-id-type="doi">10.1371/journal.pcbi.1000242</pub-id></element-citation></ref><ref id="bib42"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wadhams</surname><given-names>GH</given-names></name><name><surname>Armitage</surname><given-names>JP</given-names></name></person-group><year>2004</year><article-title>Making sense of it all: bacterial chemotaxis</article-title><source>Nature Reviews Molecular Cell Biology</source><volume>5</volume><fpage>1024</fpage><lpage>1037</lpage><pub-id pub-id-type="doi">10.1038/nrm1524</pub-id></element-citation></ref><ref id="bib43"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yuan</surname><given-names>J</given-names></name><name><surname>Berg</surname><given-names>HC</given-names></name></person-group><year>2013</year><article-title>Ultrasensitivity of an adaptive bacterial motor</article-title><source>Journal of Molecular Biology</source><volume>425</volume><fpage>1760</fpage><lpage>1764</lpage><pub-id pub-id-type="doi">10.1016/j.jmb.2013.02.016</pub-id></element-citation></ref><ref id="bib44"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yuan</surname><given-names>J</given-names></name><name><surname>Branch</surname><given-names>RW</given-names></name><name><surname>Hosu</surname><given-names>BG</given-names></name><name><surname>Berg</surname><given-names>HC</given-names></name></person-group><year>2012</year><article-title>Adaptation at the output of the chemotaxis signalling pathway</article-title><source>Nature</source><volume>484</volume><fpage>233</fpage><lpage>236</lpage><pub-id pub-id-type="doi">10.1038/nature10964</pub-id></element-citation></ref></ref-list></back><sub-article article-type="article-commentary" id="SA1"><front-stub><article-id pub-id-type="doi">10.7554/eLife.01916.025</article-id><title-group><article-title>Decision letter</article-title></title-group><contrib-group content-type="section"><contrib contrib-type="editor"><name><surname>Laub</surname><given-names>Michael</given-names></name><role>Reviewing editor</role><aff><institution>Massachusetts Institute of Technology</institution>, <country>United States</country></aff></contrib></contrib-group></front-stub><body><boxed-text><p>eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see <ext-link ext-link-type="uri" xlink:href="http://elife.elifesciences.org/review-process">review process</ext-link>). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.</p></boxed-text><p>Thank you for sending your work entitled “A simple mapping between cell swimming behavior and single-motor state in multi-flagellated <italic>E. coli</italic>” for consideration at <italic>eLife</italic>. Your article has been favorably evaluated by a Senior editor and 3 reviewers, one of whom is a member of our Board of Reviewing Editors.</p><p>The following individuals responsible for the peer review of your submission have agreed to reveal their identity: Michael Laub (Reviewing editor); Howard Berg (peer reviewer).</p><p>The Reviewing editor and the other reviewers discussed their comments before we reached this decision, and the Reviewing editor has assembled the following comments to help you prepare a revised submission.</p><p>In general, all three reviewers enjoyed reading this manuscript and the consensus was that it should ultimately be published in <italic>eLife</italic>; we thought it was a well-executed and important study with broad appeal. We had only a few minor issues to address in a revised manuscript:</p><p>1) Change the title to something that better captures the big picture, take-home message of the paper. The current title seems tailored for chemotaxis aficionados. Perhaps something along the lines of what was submitted as the “Impact Statement” which was “The swimming behavior of <italic>E. coli</italic> is surprisingly robust against natural variations in the number of flagella per cell”. Maybe just remove the word ‘surprisingly’ and it would be a title that captures the major conclusion of the paper and might appeal to a broader audience.</p><p>2) It was noted that the generalized veto model holds only up to 4 flagella per cell. It would be prudent to note this and discuss the disagreement when there are 5 or more flagella.</p><p>3) One reviewer thought it would be useful to provide a simple empirical formula for calculating <italic>N</italic><sub><italic>eff</italic></sub> from <italic>N</italic><sub><italic>flag</italic></sub> so that modellers in particular could implement the conversion as needed.</p></body></sub-article><sub-article article-type="reply" id="SA2"><front-stub><article-id pub-id-type="doi">10.7554/eLife.01916.026</article-id><title-group><article-title>Author response</article-title></title-group></front-stub><body><p><italic>1) Change the title to something that better captures the big picture, take-home message of the paper. The current title seems tailored for chemotaxis aficionados. Perhaps something along the lines of what was submitted as the “Impact Statement” which was “The swimming behavior of E. coli is surprisingly robust against natural variations in the number of flagella per cell”. Maybe just remove the word ‘surprisingly’ and it would be a title that captures the major conclusion of the paper and might appeal to a broader audience</italic>.</p><p>We thank the reviewers for this suggestion. We have changed our title to capture the big-picture message of the paper better. The new title is “<italic>Escherichia coli</italic> swimming is robust against variations in flagellar number.”</p><p><italic>2) It was noted that the generalized veto model holds only up to 4 flagella per cell. It would be prudent to note this and discuss the disagreement when there are 5 or more flagella</italic>.</p><p>This is a good point. In <xref ref-type="fig" rid="fig3">Figure 3C</xref>, both strains deviate from the generalized veto models for cells with more than 4 flagella. We have added a discussion of the possible explanations for this deviation in the text.</p><p>It is possible that hydrodynamic effects must be taken in consideration in cells with many flagella. Hydrodynamic coupling in the Hu model leads to a lower <italic>N</italic><sub><italic>eff</italic></sub>, in the direction of the deviation. Alternatively, a mechanism as described by Sneddon et al. in which cells with many flagella can run while some of its flagella rotate CW in the <italic>curly-1</italic> state, could lead to a similar deviation (see <xref ref-type="fig" rid="fig1s5">Figure 1–figure supplement 5</xref>). Finally, we must consider the possibility that the apparent deviation is due to systematic experimental error, since it is more difficult to determine visually the state of each flagellum on cells with many flagella.</p><p><italic>3) One reviewer thought it would be useful to provide a simple empirical formula for calculating N</italic><sub><italic>eff</italic></sub> <italic>from N</italic><sub><italic>flag</italic></sub> <italic>so that modellers in particular could implement the conversion as needed</italic>.</p><p>We performed a fit to the wild-type data in <xref ref-type="fig" rid="fig3">Figure 3C</xref> to provide a simple empirical formula for calculating <italic>N</italic><sub><italic>eff</italic></sub> from <italic>N</italic><sub><italic>flag</italic></sub>. The power law <italic>N</italic><sub><italic>eff</italic></sub> = 1.27 x <italic>N</italic><sub><italic>flag</italic></sub><sup>0.5</sup> provides a reasonable match to the data and is provided in the Results section as part of the description of <xref ref-type="fig" rid="fig3">Figure 3C</xref>. Additionally, we have added <xref ref-type="fig" rid="fig3s1">Figure 3–figure supplement 1</xref> showing this empirical curve with the wild-type data.</p></body></sub-article></article>