Coordination of two opposite flagella allows highspeed swimming and active turning of individual zoospores
Abstract
Phytophthora species cause diseases in a large variety of plants and represent a serious agricultural threat, leading, every year, to multibillion dollar losses. Infection occurs when their biflagellated zoospores move across the soil at their characteristic high speed and reach the roots of a host plant. Despite the relevance of zoospore spreading in the epidemics of plant diseases, individual swimming of zoospores have not been fully investigated. It remains unknown about the characteristics of two opposite beating flagella during translation and turning, and the roles of each flagellum on zoospore swimming. Here, combining experiments and modeling, we show how these two flagella contribute to generate thrust when beating together, and identify the mastigonemesattached anterior flagellum as the main source of thrust. Furthermore, we find that turning involves a complex active process, in which the posterior flagellum temporarily stops, while the anterior flagellum keeps on beating and changes its gait from sinusoidal waves to power and recovery strokes, similar to Chlamydomonas’s breaststroke, to reorient its body to a new direction. Our study is a fundamental step toward a better understanding of the spreading of plant pathogens’ motile forms, and shows that the motility pattern of these biflagellated zoospores represents a distinct eukaryotic version of the celebrated ‘runandtumble’ motility class exhibited by peritrichous bacteria.
Editor's evaluation
The authors present a study of the swimming behaviour of the zoospores of the water mold Phytophthora (Greek "Plant Destroyer"), which is responsible for significant crop damage worldwide. The motility of the zoospores is likely a significant contributor to the successful spread of the disease, and as such its study has potential wide impact. The authors suggest using a model that the anterior "hairy" (covered with mastigonemes) flagellum is the primary contributor to motility, and show with highspeed imaging that the microorganism is able to turn on the spot by stopping its posterior flagellum, and changing the beatpattern of its anterior flagellum from "spermlike" to "Chlamydomonaslike".
https://doi.org/10.7554/eLife.71227.sa0eLife digest
Microorganisms of the Phytophthora genus are serious agricultural pests. They cause diseases in many crops, including potato, onion, tomato, tobacco, cotton, peppers, and citrus. These diseases cause billions of dollars in losses each year. Learning more about how the tiny creatures disseminate and reach host plants could help scientists develop new ways to prevent such crop damage.
The spore cells of Phytophthora, also known as zoospores, have two appendages called flagella on their bodies. A tinselshaped flagellum is near the front of the creature and a long smooth filamentlike flagellum is near the posterior. Zoospores use their flagella to swim at high speeds through liquid toward potential plant hosts. Their complex swimming patterns change in response to different physical, chemical, and electrical signals in the environment. But exactly how they use their flagella to generate these movements is not clear.
Tran et al. reveal new details about zoospore locomotion. In the experiments, Tran et al. recorded the movements of zoospores in a tiny ‘swimming pool’ of fluid on top of a glass slide and analyzed the movements using statistical and mathematical models. The results uncovered coordinated actions of the flagella when zoospores swim in a straight line and when they turn. The tinsellike front flagellum provides most of the force that propels the zoospore forward. To do this, it beats with an undulating wave pattern. It shifts the beating to a breaststroke pattern to change direction. The posterior flagellum provides a smaller forward thrust and temporarily pauses during turns.
The study provides new details about zoospore’s movements that may help scientists develop new strategies to control these pests. It also offers more information about how flagella coordinate their actions to switch speeds or change directions that may be of interest to other scientists studying organisms that use flagella to move.
Introduction
Life of swimming microorganisms in viscositydominant world has been of great interest in biophysics research. The problems on microbial locomotion of those tiny individual flagellated swimmers are still far to be fully understood. There have been a multitude of theoretical and experimental models of microswimmers that study the hydrodynamics of the individual and collective motions of those cells (Ghanbari, 2020; Lauga and Goldstein, 2012). These swimming cells can be categorized into two groups: eukaryotes (having nuclei) and prokaryotes (no nuclei). Escherichia coli is one of the most studied prokaryotic swimmers, which possesses a bundle of passive helical flagella controlled by a rotary motor attached to the cell body (Turner et al., 2000). Eukaryotic microswimmers, such as green algae Chlamydomonas (Mitchell, 2001; Polin et al., 2009) and spermatozoa (Ishijima et al., 1986; Mortimer et al., 1997), have active and flexible flagella along which molecular motors are distributed. Here, we introduce a new type of microswimmer, named Phytophthora zoospores, which has two different flagella collaborating for unique swimming and turning mechanisms (Figure 1(A)).
Phytophthora is a genus of eukaryotic and filamentous microorganisms. They are classified as oomycetes and grouped in the kingdom of the Stramenopiles with the heterokont algae (such as diatoms and brown algae) (Burki, 2014; Derelle et al., 2016). A number of Phytophthora species are plant pathogens and cause tremendous damages to agro and ecosystems (Derevnina et al., 2016; Kamoun et al., 2015). Nowadays, Phytophthora diseases are responsible for a big impact on economies with billions of dollars of damages each year and remain a threat to the food security worldwide (Drenth and Sendall, 2004; Judelson and Blanco, 2005; Strange and Scott, 2005). The diseases are pervasive as they release swimming biflagellated spores called ‘zoospores’ which initiates the spreading through water. These zoospores are able to achieve speed up to $250\mathrm{\mu}\mathrm{m}{\mathrm{s}}^{1}$(Appiah et al., 2005) through thin water films, water droplets on leaves, or through pores within moist soils. To facilitate the spreading, their cell bodies store an amount of energy (mycolaminarin, lipid) allowing them to swim continuously for several hours (Judelson and Blanco, 2005). In natural ecosystems and even more in agrosystems, putative host plants are usually close. This proximity makes the distance to find a plant relatively short and it is compatible to the timeability of zoospores to swim. When the zoospores reach plant roots, they stop swimming and release their flagella to produce a primary cell wall and become germinative cysts which are able to penetrate into the host tissue. Then, they start a hyphal growth inside the infected plant. In this study, we investigate the telluric species P. parasitica, a polyphagous pathogen attacking a wide range of hosts such as potato, onion, tomato, tobacco, ornamentals, cotton, pepper, citrus plants, and forest ecosystems (Panabieres et al., 2016).
Previous studies have shown that during the spreading and approaching the host, zoospores can have complicated swimming patterns and behaviors as they experience multiple interactions with environmental signals, both physical, electrical and chemical, in soil and hostroot surface (Bassani et al., 2020a). Near the plantroot, zoospores can perceive various stimuli from the environment, such as ion exchange between soil particles and plant roots, the chemical gradients generated by root exudates, which activate cell responses. This results in coordinated behaviors of zoospores, allowing them to preferentially navigate to the water film at the interface between soil particles and plant roots. For instance, potassium, which is uptaken by roots in the soil, reduces zoospore swimming speed, causes immediate directional changes and also results in perpetual circle trajectories (Appiah et al., 2005; Galiana et al., 2019). Bassani et al. provide transcriptomic studies showing that potassium induces zoospore aggregation, which facilitates the advantages for zoospores to attack the hostroot (Bassani et al., 2020b). Experimental evidence has demonstrated that zoosporezoospore interaction can lead to ‘pattern swimming’, a microbial bioconvection happened without the appearance of chemical or electrical signals (Ochiai et al., 2011; Savory et al., 2014). These findings urge for a better understanding of the swimming physics of the individual zoospores and how the combination of their two heterogeneous flagella results in those complex swimming behaviors.
A zoospore is usually about 10 µm in size and has a kidneylike cell body (Mitchell et al., 2002; Walker and van West, 2007). At least two unique traits distinguish zoospores from other prokaryotic and eukaryotic microswimmers currently studied using physical approaches: (i) The two flagella beat longitudinally along the anteriorposterior axis of the cell body and not laterally as in the case of the green algae Chlamydomonas; (ii) two flagella distinguished from each other as the anterior flagellum has a tinsellike structure, while the posterior flagellum has a smooth whiplash one, both beating periodically with wave propagation directions outwards the body. At first sight, these two flagella seem to be competing each other due to the opposite wave propagation directions. Contrarily, multiple mastigoneme structures on the anterior flagellum of zoospores are shown to have thrust reversal ability, which makes both flagella generate thrust in the same direction and propel the cell body forwards (Bassani et al., 2020b; Cahill et al., 1996). Although it has been known about how zoospores swim, characteristics of the swimming and the beating flagella have not been statistically reported. The effects of mastigonemes on zoospore swimming also need to be carefully investigated since the mechanical properties of mastigonemes such as size, rigidity, density, can affect the swimming differently (Namdeo et al., 2011). For instance, although mastigonemes are shown to generate thrust reversal in P. palmivora zoospores (Cahill et al., 1996), they do not contribute to enhance swimming of C. reinhardtii (Amador et al., 2020).
In other microswimmers, their flagella are often synchronized to perform a cooperative swimming when they are a few microns away from each other (Friedrich and Jülicher, 2012; Lauga and Goldstein, 2012). For examples, C. reinhardtii performs breaststroke swimming by two flagella drawing away and back to each other (Polin et al., 2009) or E. coli, B. subtilis bacteria’s flagella form a bundle and rotate together like a corkscrew to propel the cell body (Turner et al., 2000). The question that whether zoospores, as eukaryotic swimmers, possess the similar cooperative behaviors of their flagella is of good interest. It has been previously claimed that zoospore flagella are independent on each other and able to perform different tasks. Carlile, 1983 describes that the anterior flagellum is responsible for pulling the zoospore through water whereas the posterior flagellum acts as a rudder for steering the cell. However, Morris et al., 1995 observe P. palmivora zoospores stop momentarily and then selforientate their bodies to a new direction relatively to the posterior flagellum. Nevertheless, the cooperative actions of the motor and rudder have not been carefully observed nor investigated, which remains unclear about how zoospores change direction either by random walks or in response to chemical and physical environment. This motivates us to unveil the physics behind individual swimming of zoospores.
In this article, we first investigate characteristics of zoospore trajectories at a global scale, then focus on the flagella scale’s swimming mechanisms. We observe that zoospores can perform long and stable straight runs, discontinued by active turning events. We obtain statistics of the trajectories and develop a numerical model to study and extrapolate the zoospore spreading characteristics solely by random walks. Then, we detail an indepth study on the hydrodynamics of P. parasitica’s flagella and acquire a mathematical model to correlate the functions of two flagella on the motion of straight runs. Although theoretical models for microswimmers with single mastigonemesattached flagella have been formulated (Brennen, 1975; Namdeo et al., 2011), models for microswimmers with two heterokont flagella have yet been considered as in case of zoospores. Here, we use Resistive Force Theory and further develop the model of a single flagellum with mastigonemes (Brennen, 1975; Namdeo et al., 2011) to adapt it with another smooth flagellum and a cell body, using a hypothesis of no interactions between two flagella. Moreover, we discover a unique active turning mechanism of zoospores including a body rotation then steering to a new direction, which results from the instantanous gait changing ability of their anterior flagellum. Our study reveals the mechanism and characteristics of zoospore spreading, which provides better insights on understanding and control of Phytophthora diseases.
Results and discussion
Characteristics of P. parasitica’s cell body and flagella
To understand the swimming, we first look at the cell body and flagellar structures of P. parasitica. By using Scanning Electron Microscopy (EM), we are able to observe the shape of the cell body and the positions of the flagellar base (Figure 1(B)). The cell in general has an ellipsoidal shape with tapered heads and a groove along the body. The size of the body is measured to be $8.8\pm 0.4\mathrm{\mu}\mathrm{m}$ (SEM) in length, and $4.7\pm 0.1\mathrm{\mu}\mathrm{m}$ (SEM) in width. The anterior flagellum attaches to the cell body in a narrow hole at one side of the groove, possessing an average length of $15.5\pm 0.1\mathrm{\mu}\mathrm{m}$ (SEM). The posterior flagellum has the same diameter as the anterior’s ($0.3\mathrm{\mu}\mathrm{m}$) but it is longer ($20.3\pm 0.76\mathrm{\mu}\mathrm{m}$ (SEM)), and attaches directly to the surface of the groove. The roots of two flagella are apart from each other with a distance of $2.9\pm 0.1\mathrm{\mu}\mathrm{m}$ (SEM). Some mastigoneme structures were observed on the anterior flagellum, but could only be distinguished with difficulty by Scanning EM technique.
Observing the zoospores with Transmission Electron Microscopy (TEM) after negative staining, we discover multiple mastigonemes on both flagella (Figure 1 (C13)). There are two different types of mastigonemes on the anterior flagellum: (type1) straight and tubular shape, high density (∼ 13 per $\mathrm{\mu}\mathrm{m}$), $0.03\mathrm{\mu}\mathrm{m}$ diameter, $1.5\mathrm{\mu}\mathrm{m}$ long; (type2) curved and irregular shape, longer and thicker in size ($0.1\mathrm{\mu}\mathrm{m}$ diameter, $1.8\mathrm{\mu}\mathrm{m}$ long), and randomly distributed (Figure 1 (C2)). The posterior flagellum instead has a smooth whip shape with plenty of very fine hairs on the surface (Figure 1 (C3)). These hairs wrap around the flagellum to increase the contact surface, thus increasing the propulsion efficiency (Lee, 2018). We also see a few type2 mastigonemes on the posterior flagellum but they only appear near the root. The function of type2 mastigonemes (Figure 1 (C4)) is unknown, but their flexibility and random arrangement suggest that they might not contribute to generate drag. In contrary, the type1 mastigonemes (Figure 1 (C5)) are tripartite hairs that occur in most of Stramenopiles kingdom. They are known to be able to generate increased drag and reverse thrust for the anterior flagellum (Cahill et al., 1996; Dodge, 2012).
Statistics of individual swimming patterns
We investigate the characteristics of zoospore swimming by analyzing their trajectories and behaviors in water. To facilitate that, we perform microscopic assays where a low concentration of individual zoospores are released to an open thin film of water with thickness of ∼100 µm on a glass slide. The setup of the water thin film can be visualized as a ‘swimming pool’ that is not covered as we want to avoid the unwanted physical interactions of zoospores with the top when they experience aerotaxis. The zoospore swimming is captured at 60 fps (interval time between two consecutive frames, $\mathrm{\Delta}t\approx 0.0167\mathrm{s}$). The images are processed by Fiji (Schindelin et al., 2012) and Trackmate plugin (Tinevez et al., 2017) to semiautomatically track the positions of zoospores during the experiment duration (see Video 1). Figure 2(A) illustrates the trajectories of zoospores captured from the microscopic assay. These trajectories indicate that zoospore can perform long and straight runs, some can even cross the whole observed region. The straight runs are separated by multiple turning events when zoospores randomly change directions. With this swimming strategy, zoospores can be categorized as runandtumble active particles (Khatami et al., 2016). From the position data of each zoospore over time, we achieve its movement characteristics defined by two parameters: magnitude of speed $U$ and moving directions θ (Figure 2(B–C)), after applying moving average method with step length $n=12$ to improve the accuracy of moving direction and instantaneous speed of the zoospore. From $U$ values, we can separate the movement of zoospores into 2 states: running state during straight runs and stopping state at turning events. While running, $U$ and θ vary around a constant value. At turning events, $U$ drops drastically, (occasionally close to 0) then quickly recovers, θ also rapidly changes to a new value. ${U}_{\mathrm{th}}$ is defined as the threshold speed that separates the two states of running and turning. With ${U}_{\mathrm{th}}$, we determine two important parameters of zoospore swimming: (i) running time ${\tau}_{\mathrm{r}}$ as the duration when $U\ge {U}_{\mathrm{th}}$, deciding how long a zoospore is able to travel without turning; (ii) stopping time ${\tau}_{\mathrm{s}}$ as the duration, when $U\le {U}_{\mathrm{th}}$, for a zoospore to perform a turn.d
We plot the distribution of $U$ for all the trajectories of zoospores swimming in the observed region for duration of 60 s (total number of zoospores $N=58$) in Figure 2(D). The speed distribution $p(U)$ exhibits a combination of two different normal distributions: ${f}_{1}(U)=0.082{e}^{{(U{\mathrm{\mu}}_{1})}^{2}/2{\sigma}_{1}^{2}}$ with ${\mathrm{\mu}}_{1}=176.6$, ${\sigma}_{1}=22.3\mathrm{\mu}\mathrm{m}{\mathrm{s}}^{1}$, and ${f}_{2}(U)=0.043{e}^{{(U{\mathrm{\mu}}_{2})}^{2}/2{\sigma}_{2}^{2}}$ with ${\mathrm{\mu}}_{2}=110$, ${\sigma}_{2}=53\mathrm{\mu}\mathrm{m}{\mathrm{s}}^{1}$. This bimodal distribution of $U$ indicates that zoospore speed fluctuates around two speed values $U={\mathrm{\mu}}_{1}$ and $U={\mathrm{\mu}}_{2}$, corresponding to two behavioral states of running and turning. The distribution f_{1} is associated with the running state where zoospores experience stable moving speed, while f_{2} represents the speed at turning events where zoospores reduce their speed from a stable running speed to 0 then quickly recover. We achieve the fitting curve for $p(U)$, resulting from the sum of two Gaussian fits ${f}_{1}+{f}_{2}$, and choose the speed at inflection point of the fitting curve where ${\mathrm{\mu}}_{\mathrm{s}}\le U\le {\mathrm{\mu}}_{\mathrm{r}}$ as the speed threshold to separate two behavioral states, $U}_{\mathrm{t}\mathrm{h}}=111.5\phantom{\rule{thinmathspace}{0ex}}\mu \mathrm{m}\phantom{\rule{thinmathspace}{0ex}}{\mathrm{s}}^{1$ at the inflection point determines a turning point for a significant change of the speed from running to turning state. The sensitivity of our ${U}_{\mathrm{th}}$ selection can be tolerated by ±10 % of the chosen value, ranging from 100 to 122.5 $\mathrm{\mu}\mathrm{m}{\mathrm{s}}^{1}$ (See Appendix 5). We also plot the polarity distribution of zoospores in Figure 2(E) based on the moving direction θ and acquire an equally distributed in all directions. With the defined ${U}_{\mathrm{th}}$, we calculate and plot the distributions of ${\tau}_{\mathrm{r}}$ and ${\tau}_{\mathrm{s}}$ in form of survival curves $p(\tau \ge t)$ in Figure 2(F). Both survival curves show complex behaviors of zoospores during running and turning. The statistics $p({\tau}_{\mathrm{s}}\ge t)$ can be considered as sum of two exponential decays with average stopping time ${\overline{\tau}}_{\mathrm{s}}=0.37\mathrm{s}$. Also, $p({\tau}_{\mathrm{r}}\ge t)$ is in form of two exponential decays with average running time ${\overline{\tau}}_{\mathrm{r}}=1.0\mathrm{s}$. We can estimate that zoospores stop and turn with the frequency 1/${\overline{\tau}}_{\mathrm{r}}$ greater than $1\mathrm{Hz}$ ($p({\tau}_{\mathrm{r}}<1.0\phantom{\rule{thinmathspace}{0ex}}\mathrm{s})=0.82$). Based on the moving direction over time, we calculate the average turning speed of zoospores at $\overline{\dot{\theta}}\approx 0.6\pi \mathrm{rad}{\mathrm{s}}^{1}$. At each stopping time they perform a turning angle $\mathrm{\Delta}\theta ={\theta}_{\mathrm{i}}{\theta}_{\mathrm{e}}$, where ${\theta}_{\mathrm{i}}$ and ${\theta}_{\mathrm{e}}$ is the moving direction right before and after each turning event, respectively. The distribution of $\mathrm{\Delta}\theta $ is shown in Figure 2(G), demonstrating the equal preference of turning directions. It is also shown that zoospores preferentially turn with the angle around 0°, which we speculate that it results from the failed outofplane movement when zoospores swim near the water/air interface during their aerotaxis (See Video 2).
Since the motion of zoospore is characterized by the succession of straight runs and turning events, as illustrated in Figure 2(H), in order to quantify their largescale transport properties, we assemble all previous measurements in the following way. Each straight run is characterized by a speed $U}_{\mathrm{r}}>{U}_{\mathrm{t}\mathrm{h}$ and a duration ${\tau}_{\mathrm{r}}$ drawn from the distributions in Figure 2(D), and (F), respectively. After a run phase, an idle phase of duration ${\tau}_{\mathrm{s}}$, drawn from Figure 2(F) follows. The moving direction of the $r$th run phase is given by $\mathrm{cos}({\theta}_{r})\widehat{x}+\mathrm{sin}({\theta}_{r})\widehat{y}$, i.e. it is parameterized by an angle ${\theta}_{r}$. Note that moving direction of two consecutive $r$th and $(r+1)$th run phases are correlated. Moreover, ${\theta}_{r+1}={\theta}_{r}+\mathrm{\Delta}\theta $, where $\mathrm{\Delta}\theta $ is a random angle drawn from the distribution in Figure 2(G). Mathematically, the position $\mathbf{x}}_{m$ of the zoospore after $m$ run phases is given by ${\mathbf{x}}_{m}=\sum _{r=1}^{m}{U}_{r}\phantom{\rule{thinmathspace}{0ex}}{\tau}_{r}[\mathrm{cos}({\theta}_{r})\hat{x}+\mathrm{sin}({\theta}_{r})\hat{y}]$, and its meansquare displacement is $\mathrm{M}\mathrm{S}\mathrm{D}(m)=\u27e8({\mathbf{x}}_{m}\u27e8{\mathbf{x}}_{m}\u27e9{)}^{2}\u27e9$, where $\u27e8\mathrm{\cdots}\u27e9$ denotes average over realizations of the process. Our simulation results in a diffusive behavior of zoospores, with the MSD proportional to $t$ (see Video 3). The diffusion coefficient is then obtained from $D=\underset{m\to \mathrm{\infty}}{lim}\frac{\mathrm{M}\mathrm{S}\mathrm{D}(m)}{4m(\u27e8{\tau}_{r}\u27e9+\u27e8{\tau}_{s}\u27e9)}$. In the computation of $\mathrm{MSD}$ and $D$, we assume that the only random variable exhibiting correlations is ${\theta}_{r}$, while ${U}_{r}$, ${\tau}_{r}$, and ${\tau}_{s}$ are uncorrelated. This procedure allows us to obtain an estimate of $D$, with $D=3.5\times {10}^{4}{\mathrm{cm}}^{2}{\mathrm{s}}^{1}$, which is in the same order of magnitude as the diffusion coefficient of C. reinhardtii’s (Polin et al., 2009). We stress that direct measurements of $D$ based on experimental MSD data are highly unreliable given the relatively small number of trajectories and their short duration. In our case, the data are enough to obtain reliable estimates on the distributions of $\mathrm{\Delta}\theta $, ${U}_{r}$, ${\tau}_{r}$, and ${\tau}_{s}$ from which, as explained above, $D$ can be reliably estimated from simulations. Similar methods of using a theoretical random walk model to estimate the macroscopic parameters from the microscopic experiment have been previously developed for C. reinhardtii (Hill and Hader, 1997; Vladimirov et al., 2004). We emphasize that $D$ represents the estimation of diffusion coefficient of individual swimming of zoospores from random walks. This is more to show the intrinsic ability of individual zoospores to perform spatial exploration, rather than to quantify the bulk diffusivity where the collective swimming behaviors, which involve zoosporezoospore interactions, play a major role.
Role of two flagella in swimming motions of zoospores
Our statistics study on swimming trajectories of zoospores has delivered characteristics of their movement at largescale, including the straight runs and turning events. These motions are controlled by two flagella oriented in opposite directions along the cell body’s anteriorposterior axis. In this section, we look indepth to how these two flagella together generate speed and perform turning for zoospores by conducting microscopic assays at small length scale, of which the flagella are visible and in very short time scale.
Straight runs
We record movement of zoospores during their straight runs with visible flagella by conducting brightfield microscopy with 40× objective and a highspeed camera capturing at 2000 fps at exposure time 200 $\mathrm{\mu}\mathrm{s}$. In Figure 3(A), we show images of a P. parasitica zoospore swimming by two flagella beating in sinusoidal shapes with the wave propagation in opposite directions. While translating, the cell body gyrates around the moving direction simultaneously, which results in a helical swimming trajectory (see Video 4 for a long run of a zoospore swimming in water). We believe that this gyrational motion might result from the intrinsic chiral shape of the zoospore body and offaxis arrangement of their flagella (Figure 1(B)). Indeed, previous studies have shown that chirality of a microswimmer’s body induces spontaneous axial rotation resulting from the translational motion (Keaveny and Shelley, 2009; Löwen, 2016; Namdeo et al., 2014). From multiple observations, we obtain the pitch and radius of the helical trajectories at $p=130\pm 8\mathrm{\mu}\mathrm{m}$ (SEM) and $R=4.0\pm 0.2\mathrm{\mu}\mathrm{m}$ (SEM), respectively (data presented in Appendix). We then estimate the gyrational speed of the cell body $\dot{\varphi}=2\pi /\mathrm{\Delta}{t}_{p}$, where $\mathrm{\Delta}{t}_{p}$ is the duration the zoospore travels through a full turn of the helical path (Figure 3(B)). We obtain $\dot{\varphi}=(3.6\pm 0.3)\pi \mathrm{rad}{\mathrm{s}}^{1}$ (SEM) (see Appendix). The observations of helical trajectories also confirm that each flagellum of zoospores beat as a flexible oar in a 2D plane as we notice the two flagella flattened into two straight lines during the gyration. Thus, zoospores are not expected to swim in circles when interacting with noslip boundaries as seen in case of E. coli with a rotating flagellar bundle. The ‘curved straight runs’ of zoospores that we observed in Figure 2(A) might result from rotational diffusion and thermal fluctuations.
We retrieve the parameters of the beating flagella including beating frequencies $f$ and wavelengths λ, by applying kymographs on the crosssections within the two flagella and normal to the moving direction of the zoospore. We present more details of the parameter retrieval with kymographs in Appendix. Our data show that when swimming in straight runs, the anterior flagellum of zoospores usually beats at ${f}_{1}\approx 70\mathrm{Hz}$ while the posterior flagellum beats with the frequency ${f}_{2}\approx 120\mathrm{Hz}$ that is approximately 1.7fold faster than the anterior’s. We focus on a very short duration of less than 50 ms, which is equivalent to a translation of less than 10 µm. Compared with the rotation motion at time scale of 1 s, we neglect the effect of the cell body’s rotation in this short duration.
We then develop a mathematical model to study how the dynamics of beating flagella helps generating thrust for the cell to move forward. We assume that the gyration of the cell body does not affect the shapes and motions of the flagella since the beating frequencies of the two flagella are much higher than the gyrational speed. The gyration also does not contribute to the translation as we consider it as a passive motion resulting from the chirality of zoospores. Thus, the swimming zoospore can be considered as a 2D model (Figure 3(C)) in which the cell body is an ellipse defined as ${x}^{2}/{a}^{2}+{y}^{2}/{b}^{2}=1$ in its bodyfixed frame ($xOy$), with the anterior and posterior flagellum having sine waveform shapes defined as
as ${(1)}^{k}x\le c$, where ${A}_{k}$ is the amplitude, ${\omega}_{k}$ is the angular speed, ${\lambda}_{k}$ is the wavelength of the anterior ($k=1$) and the posterior flagellum ($k=2$). The two flagella are attached to the cell body at two points both lying on $x$axis and distanced to the origin $O$ a gap of $c$, and beating with the wave propagation ${\overrightarrow{v}}_{\mathrm{w}}$ having directions against each other. The anterior has length L_{1} and diameter d_{1} while those of the posterior are L_{2} and d_{2}. Additionally, there are multiple tubular type1 mastigonemes with length $h$ and diameter ${d}_{\mathrm{m}}$ attached to the surface of the anterior flagellum with density ${N}_{\mathrm{m}}$ indicating the number of mastigonemes attached on a unit length of the flagellum. It is important to determine the flexibility of these mastigonemes as it would impact the ability of the mastigonemes to produce drag. We estimate the flexibility by a dimensionless parameter, which has been carefully characterized in previous studies (Khaderi et al., 2009; Namdeo et al., 2011), ${F}_{\mathrm{m}}=12\mu KA\omega {h}^{3}/(E{d}_{\mathrm{m}}^{3})$, where μ is the fluid viscosity, $K=2\pi /\lambda $ is the wave number, $E$ is the Young modulus of the mastigonemes. With this estimation, if ${F}_{\mathrm{m}}<0.1$, mastigonemes are considered as fully rigid. In case of zoospores’ mastigonemes, we achieve ${F}_{\mathrm{m}}$ at order of 10^{4}, which is much lower than 0.1. Thus, we can assume that the mastigonemes of zoospores are nondeformable and rigidly attached to the anterior flagellum. As a result, hydrodynamic interactions between neighboring mastigonemes can also be neglected. Additionally, we ignore the effects of the type2 mastigonemes in producing drag due to their nontubular and random structures.
Zoospores swim in water with very low Reynolds number ($Re<<1$), resulting in negligible inertia, dominant viscous force and the kinetic reversibility (Purcell, 1977; Qiu et al., 2014). Microswimmers with flexible flagella generate thrust from drag force acted by fluid on the flagellum segments. In our model, we use Resistive Force Theory (RFT) to deal with the calculation of fluid’s drag force on the two flagella of the zoospore. RFT has proven to be an effective and accurate method to predict the propulsive force and velocity of microswimmers regardless of the interactions of flagellumflagellum or flagellumbody (Gray and Hancock, 1955; Lauga and Powers, 2009; Marcos et al., 2014). In case of zoospores where two flagella are in opposite directions, and the flagellumbody interaction is insignificant, RFT is a suitable solution to apply. Following this method, each flagellum is divided into an infinite numbers of very small segments with length ds_{k}, and each segment is located in the bodyfixed frame ($xOy$) by a position vector ${\overrightarrow{r}}_{k}={x}_{k}\overrightarrow{i}+{y}_{k}\overrightarrow{j}$, where $\overrightarrow{i}$, $\overrightarrow{j}$ are the unit vectors in $x$ and $y$direction, respectively; x_{k} and y_{k} satisfy the shape equation (Equation 1).
RFT states that the drag force by fluid acting on an infinitesimal segment ds of the flagellum is proportional to the relative velocity of fluid to the flagellum segment (Gray and Hancock, 1955; Hancock, 1953; Koh and Shen, 2016), as follows
where ${V}_{\mathrm{N}}$ and ${V}_{\mathrm{L}}$ are two components of relative velocity of fluid in normal and tangent direction to the flagellum segment, ${K}_{\mathrm{N}}$ and ${K}_{\mathrm{L}}$ are drag coefficients of the flagellum in normal and tangent to the flagellum segment, $\hat{n}$ and $\hat{l}$ are the unit vectors normal and tangent to the flagellum segment, respectively. The drag coefficients ${K}_{N}$ and ${K}_{L}$ are estimated by Brennen and Winet, 1977, which depends on fluid viscosity, the wavelength and diameter of a flagellum.
We then apply RFT on each flagellum of the zoospore to calculate the total drag force acting on it. For the posterior flagellum, each segment ds_{2} is a simple smooth and slender filament (see inset ds_{2} in Figure 3(C)), having drag coefficients ${K}_{\mathrm{N2}}$ and ${K}_{\mathrm{L2}}$. For the anterior flagellum, each segment ds_{1} contains additional ${N}_{\mathrm{m}}\mathrm{d}{s}_{1}$ mastigonemes. Using a strategy from previous models for flagella with mastigonemes (Brennen, 1975; Namdeo et al., 2011), we consider these mastigonemes stay perpendicular to the segment itself (see inset ds_{1} of Figure 3(C)) and also act as slender filaments experienced drag from water. Interestingly, due to the direction arrangement, the relative velocity normal to the flagellum segment results in drag force in tangent direction to the mastigonemes, and subsequently, the relative velocity tangent to the flagellum segment results in drag force in normal direction to the mastigonemes. In another perspective, we can consider the anterior flagellum receives additional drag from the mastigonemes, which is presented by two increased drag coefficients in normal and tangent direction, defined as
and
respectively. Here, ${K}_{\mathrm{Nf1}}$ and ${K}_{\mathrm{Lf1}}$ are the drag coefficients in normal and tangent direction of the flagellum filament, respectively; ${K}_{\mathrm{Nm1}}$ and ${K}_{\mathrm{Lm1}}$ are the drag coefficients in normal and tangent direction of the mastigonemes, respectively.
In low Reynolds number condition, total forces equate to zero due to approximately zero inertia. Thus, we derive translational velocity ${U}_{X}$ of the zoospore as shown in Equation 5
where ${v}_{\mathrm{wk}}={\lambda}_{k}{f}_{k}$ is the wave propagation velocity of the flagellum, $\gamma}_{k}={K}_{\mathrm{L}k}/{K}_{\mathrm{N}k$ is the drag coefficient ratio, ${\beta}_{k}={A}_{k}/{\lambda}_{k}$ is the flagellar shape coefficient, and ${\xi}_{\mathrm{e}}$ is the shape coefficient of the ellipse cell body (see Appendix for the detailed derivatives).
We first study the effects of mastigonemes on zoospore speed by ploting the value of ${U}_{X}$ with different density of the mastigonemes, and varying the beating frequency of the posterior flagellum f_{2} while the anterior flagellum beats with a usual ${f}_{1}=70\mathrm{Hz}$ in Figure 3(D). The dimensions and physical parameters of the zoospore’s cell body and two flagella are shown in Appendix. The plot shows that the appearance of mastigonemes results in a reversed thrust from the anterior flagellum. To illustrate this, when there are no mastigonemes on the front flagellum (${N}_{\mathrm{m}}=0$) and the posterior flagellum is excluded (${f}_{2}=0$, ${L}_{2}=0$), the anterior flagellum generates a thrust in negative $X$direction (${U}_{X}<0$). This agrees well with previous studies modeling a smooth reciprocal beating flagellum (Gray and Hancock, 1955; Hancock, 1953; Lauga and Powers, 2009). But since ${N}_{\mathrm{m}}>2$, the velocity of the zoospore is reversed to positive $X$direction due to the extra drag from the mastigonemes. This phenomenon was also described in hydrodynamics by Namdeo et al., 2011, and observed in experiment of Cahill et al., 1996. Interestingly, when changing directions, zoospores are observed to swim with the solely beating anterior flagellum to pull the body forwards while the posterior flagellum is immobile (Figure 4(A)) and Video 5, which confirms the thrust reversal ability of the mastigonemes. The beating frequency of the anterior flagellum in this case increases to ${f}_{1}\approx 110\mathrm{Hz}$. This finding also reassures the importance of mastigonemes on zoospore swimming, which is not similar to those of C. reinhardtii (Amador et al., 2020). The front flagellum of zoospores has fibrillar mastigonemes, similar to Chlamydomonas, but at higher density with tubular shape and larger in size, that could render into account of the different beating properties from the smooth posterior flagellum. However, high mastigoneme density (more than 20 per 1 µm flagellum length) shows mild effect on speed. From TEM images taken at the anterior flagellum, we estimate the mastigoneme density by averaging the number of mastigonemes manually counted over a flagellum length (See Appendix). We obtain ${N}_{\mathrm{m}}=13.0\pm 0.8\mathrm{\mu}{\mathrm{m}}^{1}$ (SD), which falls between the optimum range to generate speed (Figure 3(D)).
To understand how the coordination of two flagella influences zoospore speed, we vary the beating frequency of one flagellum while the other’s remains constant and obtain the resultant speed (Figure 3(E)). We find that although both flagella contribute to zoospore speed, the anterior flagellum has larger impact on speed than the posterior does. For instance, the anterior flagellum can singly generate a speed threefold higher than the posterior flagellum can do at the same beating frequency. Moreover, the additional speed contributed by the anterior flagellum remains almost the same regardless of the varation of frequency of the posterior flagellum, while the speed contribution of the posterior flagellum decreases as the anterior flagellum increases its frequency. Since the contribution to zoospore speed is different between two flagella, we ask whether the energy consumption of each flagellum might also be different or equally distributed. In our model, each flagellar segment consumes a power deriving from the dot product between the drag force of water acting on the segment and the relative velocity of water to the segment, which can be written as $P\equiv \int {\overrightarrow{v}}_{\mathrm{w}/\mathrm{f}}\cdot d\overrightarrow{F}$ (see Appendix for detailed derivatives). We achieve that, the power consumption of each flagellum depends on both flagella at low frequencies ($f$ ≤ 40 Hz), and becomes solely dependent on its own beating at high frequencies ($f$ ≥ 40 Hz) (Figure 3(F)). At high frequencies, the anterior flagellum consumes approximately fivefold more power than the posterior, given the same beating frequency. As a result, the total power consumed for both flagella becomes less dependent on the posterior flagellum as the anterior flagellum beats faster (Figure 3(G)). In Figure 3(H), we show the fraction of power consumed by the anterior flagellum over the total power consumption of the zoospore. We notice that, at the same beating frequency with the posterior flagellum, the anterior flagellum accounts for ∼80 % of the total power. Interestingly, when the frequency of the posterior flagellum is ∼1.7fold higher than that of the anterior flagellum, both flagella consume the same amount of power. Indeed, this result agrees well with our experimental data, in which we obtain the anterior flagellum normally beats at 70 Hz and the posterior flagellum at 120 Hz. Thus, we can speculate that the energy is equally distributed for both flagella. In addition, we also estimate the propelling efficiency of the zoospore $\eta =\frac{{P}_{0}}{{P}_{1}+{P}_{2}}$, with P_{0} as required power to move the cell body forward at speed ${U}_{X}$ (see Appendix for derivatives). We achieve that η is higher as the anterior flagellum increases its frequency (Figure 3(I)). The efficiency reaches its maximum value at ∼1.2 % when the beating frequency of the posterior is ∼1.7fold higher than that of the anterior flagellum. Overall, we show that the energy is shared in comparable manner between two flagella, but the anterior flagellum has more influence on zoospore speed, power consumption and propelling efficiency. On the other hand, the posterior flagellum provides a modest contribution to zoospore speed despite beating at higher frequency and consuming half of the energy. Taking the fact that the anterior flagellum can beat singly with the posterior flagellum staying on pause during turning events, we speculate that the anterior flagellum is the main motor of zoospores. Nevertheless, the function of the posterior flagellum remains cryptic and we investigate more on its role during turning events.
Turning events
We capture zoospores changing their directions by a unique active turning mechanism and an interesting coordination between two flagella. The videos of zoospore turning are recorded at 2000 fps with the same setup as the microscopic assay (Figure 4(A) and Video 5). Using Fiji and Trackmate, we track the positions of the zoospore and plot the trajectory, smoothed by moving average with step length $n=40$, during its turning event (Figure 4(B)). We then compute the zoospore speed $U$ (Figure 4(C)), and the moving directions θ from the trajectory while manually measure body orientation ψ over time (Figure 4(D)). We observe that at first, to prepare for the turn, the zoospore reduces its speed as both flagella beat with smaller amplitudes. Then, the posterior flagellum instantaneously stops beating, which marks the beginning of a turning event (Figure 4 (A12)). The anterior flagellum now takes full control of zoospore motion during turning. The zoospore then perform two distinct sequences of motions right after the posterior stops beating: (i) rotation of the cell body out of the old direction resulting from a few repetitive strokelike beating of the anterior flagellum (Figure 4 (A3)), then (ii) steering towards a new direction as the anterior flagellum switches back to normal beating with sinusoidal waveform to propel the body (Figure 4 (A45)). These two sequences of turning can be distinguished by two different patterns of trajectory during turning. While the rotation is indicated by multiple circular curves, in which each of them corresponds to a strokelike back and forth beating motions, the steering results in straight line (Figure 4(B)). Additionally, the zoospore speed $U$ during this period also consists of two different patterns of fluctuations (Figure 4(C)), while the moving direction θ changes from large fluctuations (for rotation) to stable direction (for steering) (Figure 4(C)). The turning event ends when the cell body stably moves in the new direction, which is indicated by the recovery of speed and the overlap of the moving direction and body orientation. The anterior flagellum continues to propel the cell body out of the location of the turning event (∼5 µm away), which we call ‘stabilize’ step (Figure 4 (A6)). We notice that in this step, the anterior flagellum beats at a higher frequency than usual (∼110 Hz, compared to normally at ∼70 Hz) and achieves a highspeed of ∼250 µm s^{1}. Finally, the posterior flagellum resumes its beating and the zoospore returns to the normal straight run state.
It is striking that the anterior flagellum is able to completely change its gait from sinusoidal traveling wave to strokelike beating during the rotation step of the turning event. Thus, we need to have experimental evidence from direct observation of the shape of the anterior flagellum during this gait changing period, which is sometimes hidden underneath the cell body due to the offaxis postion of the flagellar base, to extend our understanding on this turning behavior. We capture new videos of zoospore turning events, but the focus plane is slightly offset from the focus plane of the cell body (Figure 4(E) and Video 6). We observe that the anterior flagellum instantaneously perform continuous power and recovery strokes, similar to the breaststroke beating of Chlamydomonas’s flagella. As a result, the cell body can rotate quickly about a fixed point, which helps the zoospore make a sharp turn. We summarize the turning of zoospores by a schematic in Figure 4(F). Indeed, this rotation motion is similar to that of uniflagellate C. reinhardtii (Bayly et al., 2011), which strengthens the evidence for the actively switchable beating pattern of zoospores’ anterior flagellum. While the anterior flagellum plays a major role in turning events on top of straight runs, the posterior flagellum only stays immobile throughout the process. Here, we can confirm that the posterior flagellum does not act as a rudder to steer the direction. Instead, it is fully stretched during turning events and might contribute to increase drag at one end of the cell body. The function of the posterior flagellum is not completely clear regarding thrust production during translation or turning. Thus, we speculate that it might contribute to other nonphysical activities such as chemical and electrical sensing. The advantage of zoospore turning mechanism is that it allows them to actively and quickly achieve a new direction, which is not the case of other microswimmers such as the tumbling of E. coli’s (Darnton et al., 2007; Perez Ipiña et al., 2019) or Chlamydomonas’s (Bennett and Golestanian, 2015; Polin et al., 2009).
Conclusion
We have performed the first systematic study of the swimming pattern and spreading features of P. parasitica zoospores, a plant pathogen, which is considered a major agricultural threat. Combining highspeed imaging and Resistive Force Theory, we show how the two opposite flagella are coordinated in producing thrust by beating together, allowing the microorganism to achieve highspeed swimming during straight runs. Furthermore, we find that turning is a coordinated process, in which the the posterior flagellum stops beating, while the anterior flagellum actively moves causing the cell body to rotate. Finally, we explain how fastswimming periods and active turning events combine to produce a diffusion coefficient of $D=3.5\times {10}^{4}{\mathrm{cm}}^{2}{\mathrm{s}}^{1}$, a quantity that characterizes spatiotemporal spreading of this pathogen during plant epidemics.
It is worth stressing the motility pattern exhibited by the zoospores represents an Eukaryotic version of the ‘runandtumble’ motility class exhibited by bacteria peritrichous bacteria. Several eukaryotic swimmers, for example spermatozoa, do not exhibit such kind of motility pattern, but C. reinhardtii have been reported to also fall into this category (Polin et al., 2009). There are, however, important differences. Chlamydomonas possess two identical flagella located at one tip of the cell. Reorientation events occur during asynchronous beating periods of the two identical flagella, while straight runs require synchronous beating. In sharp contrast to this picture, we show that in zoospores while straight runs also involve a coordinated beating of the opposite and different flagella, turning involves temporary halting of posterior flagellum, while the anterior flagellum continues beating. This strongly suggests that P. parasitica navigates using a fundamentally different internal regulation mechanism to control swimming, than C. reinhardtii, a mechanism that is likely to be present in other Eukaryotic swimmers with two opposite and different flagella.
We believe our findings on the coordination of two flagella bring more insights on zoospore swimming dynamics. It was also not known from the literature about the role of each flagellum on the straight runs, and on turning events in particular. We show that although the energy is shared in comparable manner between both flagella, the anterior flagellum contributes more to zoospore speed. Zoospores actively change directions thanks to the sole beating of the anterior flagellum. We believe that anterior flagellum could be the main motor of zoospores that is in charge of generating speed, changing beating patterns from sinusoidal wave to power and recover stroke to quickly rotate the body, while the posterior flagellum might play a role in chemical/electrical sensing and providing an anchorlike turning point for zoospores, instead of acting like a rudder as previously hypothesized. These findings pave new ways for controlling the disease since now we can have different strategies on targeting one of the flagella.
Materials and methods
P. parasitica mycelium culture and zoospore release
Request a detailed protocolWe culture mycelium of Phytophthora parasitica (isolate 310, Phytophthora INRAE collection, SophiaAntipolis, France) (Bassani et al., 2020b; Galiana et al., 2019) routinely on malt agar at 24 °C in the dark. To produce zoospores, we prepare the mycelium which is grown for one week in V8 liquid medium at 24 °C under continuous light. The material is then drained, macerated and incubated for a further four days on water agar (2 %) to induce sporangiogenesis. Zoospores are released from sporangia by a heat shock procedure. We place a petridish of mycelium inside a refrigerator at 4 °C for 30 min, then pour 10 mL of 37 °C distilled water on top of the mycelium and continue to incubate it at room temperature (25 °C) for another 30 min. Zoospores then escape from sporangia and swim up to the water. The zoospore suspension is then collected for further experiments.
Scanning electron microscopy and transmission electron microscopy with negative staining
Request a detailed protocolFor Scanning Electron Microscopy and Transmission Electron Microscopy, cell pellets are fixed in a 2.5 % glutaraldehyde solution in 0.1 M sodium cacodylate buffer (pH 7.4) at room temperature (∼25 °C) for 1 h and then stored at 4 °C. For Scanning EM observations, after three rinsing in distilled water, protists are filtered on a 0.2 µm isopore filter. Samples on filters are subsequently dehydrated in a series of ethanol baths (70 %, 96 %, 100 % three times, 15 min each). After a final bath in hexamethyldisilazane (HMDS, 5 min), samples are left to dry overnight. Samples on filters are mounted on Scanning EM stubs with silver paint and coated with platinum (3 nm) prior to observing. The Scanning EM observations are performed with a Jeol JSM6700F scanning electron microscope at an accelerating voltage of 3 kV.
For TEM observations, samples are prepared using the negative staining method. After three rinsing in distilled water, a drop of cells suspension (∼10 µL) is left for 5 min on a TEM copper grid (400 mesh) with a carbon support film. The excess liquid is removed with a filter paper. Subsequently, staining is done by adding a drop of 0.5 % (w/v) aqueous solution of uranyl acetate on the grid for 1.5 min, followed by removal of excess solution. The TEM observations are carried out with a JEOL JEM1400 transmission electron microscope equipped with a Morada camera at 100 kV.
Microscopic assays of zoospores
Request a detailed protocolWe pipette a droplet of 10 µL water containing zoospores onto a microscopic glass slide and spread the droplet to thoroughly cover the marked area of 1 × 1 cm. We then achieve a thin water film of approximately 100 µm thickness. We do not put coverslips on the water film to prevent the unwanted interactions between the zoospores and rigid surface of coverslips. We observe the swimming of individual zoospores inside the water film under a bright field transmission microscope (Nikon Eclipse Ti2, Minato, Tokyo, Japan) at 40× objective with the highspeed camera Phantom v711 (Vision Research, NJ, USA). For the experiment to observe the swimming trajectories of the zoospores, we use 4× objective to capture a large swimming region of 5000 × 4000 µm. The captured images are processed by Fiji with Trackmate plugin.
Estimation of trajectory parameters
Request a detailed protocolPositions of an individual zoospore are captured at each time frame $\mathrm{\Delta}t$. At each ${t}_{j}=j\mathrm{\Delta}t$ ($j=1,2,3,\mathrm{\dots}$), the zoospore has a position ${z}_{j}=(x({t}_{j}),y({t}_{j}))$. First, we smooth the trajectory by moving average with step $n$. The smoothed positions ${Z}_{j,n}$ are calculated as,
Each new position ${Z}_{j,n}$ possesses a velocity vector with speed,
and moving direction (angle between velocity vector and $x$axis).
Appendix 1
Obtaining characteristics of beating flagella by kymograph
We first use Fiji with StackReg plugin to fix the position of the cell body in the swimming videos, and only the two flagella remain beating.
We retrieve the parameters of the beating flagella by applying kymographs on multiple crosssections which are along to or normal to the moving direction at the positions shown in the figure below. These kymographs help converting the complex timedependent 2D image data of the flagella from ($X,Y,t$) to separate signals ($X,t$) and ($Y,t$). As a result, the kymographs at crosssection (1) and (2) give us the average amplitudes $A$ of the waveform shapes of the two flagella, and the ones at (3) and (4) provide us with information of beating frequencies $f$ and wavelengths λ.
Appendix 2
Variation of ${U}_{\mathrm{th}}$
To verify if the way we chose the value of ${U}_{\mathrm{th}}$ will affect the result of turning angles and other swimming parameters, we vary the ${U}_{\mathrm{th}}$ value by ±10 % of the chosen value, which ranges from 100 to 122.5 µm s^{1}. We plot the running time ${\tau}_{\mathrm{r}}$, stopping time ${\tau}_{\mathrm{s}}$ and turning angles $\mathrm{\Delta}\theta $ according to the varied ${U}_{\mathrm{th}}$. Data are shown in the figure below. We see that the changes do not affect the results of ${\tau}_{\mathrm{r}}$, ${\tau}_{\mathrm{s}}$ and $\mathrm{\Delta}\theta $ estimation. Thus, we believe the sensitivity of ${U}_{\mathrm{th}}$ selection can be tolerated by ±10%.
Appendix 3
Characteristics of helical trajectories
From the microscopic assays with highspeed camera and 40× objective, we observe zoospores swimming in helical trajectories. We manually measure the characteristics of these helical trajectories and present the data in the table below.
Appendix 4
Estimation of mastigoneme density
Mastigoneme density is an important parameter, which determines the thrust reversal effect of the anterior flagellum and generates speed for zoospores. Using TEM, we capture images of anterior flagellum containing the mastigonemes. We first select a longitudinal flagellum length and manually count the total number of mastigonemes attached along that distance (See figure below). The mastigoneme density ${N}_{\mathrm{m}}$ is calculated as the ratio between the number of counted mastigonemes and the selected flagellum length.
Appendix 5
Hydrodynamics of two flagella in a swimming zoospore using Resistive Force Theory
Following the schematics of Figure 3(a) in the main text, we consider a parameter s_{k} as the distance along the flagellum from the root to the segment ds_{k}, then we define dimensionless constants ${\alpha}_{k}$ as the ratios of the wavelength to the arc length of the anterior ($k=1$) and posterior flagellum ($k=2$), respectively. Whereas
and
Thus, we can express the segment position ${\overrightarrow{r}}_{k}$ in terms of s_{k}, as follows
and
We achieve the velocity of the flagella relative to the cell body as
The beating of flagella results in the movement of the cell body relative to water, with a velocity:
We neglect the velocity of the cell body in $Y$direction with the assumption that the zoospore swimming is very directional and $U}_{X}>>{U}_{Y$.
Resistive Force Theory (RFT) states that the drag force by fluid acting on an infinitesimal segment d $s$ of the flagellum is proportional to the relative velocity of fluid to the flagellum segment, as follows
where ${V}_{\mathrm{N}}$ and ${V}_{\mathrm{L}}$ are two components of relative velocity of fluid in normal and tangent direction to the flagellum segment, ${K}_{\mathrm{N}}$ and ${K}_{\mathrm{L}}$ are the drag coefficients of the flagellum in normal and tangent to the flagellum segment, $\hat{n}$ and $\hat{l}$ are the unit vectors normal and tangent to the flagellum segment, respectively.
Equation (12) can be expressed as:
where ${\overrightarrow{v}}_{\mathrm{w}/\mathrm{f}}$ is the relative velocity of water to the flagellum, which is achieved as
and the tangent unit vector $\hat{l}$ is expressed as
During the derivative, we use an assumption that the flagella have small amplitude deflection ($A<<L$). Thus, we can approximate that ${\mathrm{cos}}^{2}\left(\frac{2\pi {s}_{k}}{{L}_{k}}{\omega}_{k}t\right)\approx 0.5$. This helps to simplify the term $\frac{\partial \overrightarrow{r}}{\partial s}=\sqrt{{\alpha}_{k}^{2}+4{\pi}^{2}\frac{{A}_{k}^{2}}{{L}_{k}^{2}}{\mathrm{cos}}^{2}\left(\frac{2\pi {s}_{k}}{{L}_{k}}{\omega}_{k}t\right)}=\sqrt{{\alpha}_{k}^{2}+2{\pi}^{2}\frac{{A}_{k}^{2}}{{L}_{k}^{2}}}$.
The drag coefficients ${K}_{N}$ and ${K}_{L}$ is estimated by Brennen and Winet, as follows:
where $\mu =8.9\times {10}^{4}\mathrm{Pa}\mathrm{s}$ is the water viscosity at $25{\phantom{\rule{thinmathspace}{0ex}}}^{\circ}\mathrm{C}$, λ is the wavelength and $d$ is the diameter of the flagellum. With ${U}_{Y}$ being neglected, we also consider that the two flagella do not generate force in $Y$direction, but only induce thrust in $X$direction. Hence,
We then apply RFT on each flagellum of the zoospore to calculate the total drag force acting on it. For the posterior flagellum, each segment ds_{2} is a simple smooth and slender filament, having drag coefficients ${K}_{\mathrm{N2}}$ and ${K}_{\mathrm{L2}}$. following equation (16) and (17). From equation (13) and (18), we derive the drag force of water acting on posterior flagellum as
where ${v}_{\mathrm{w2}}={\lambda}_{2}{f}_{2}$ is the wave propagation velocity, ${\gamma}_{2}={K}_{\mathrm{L2}}/{K}_{\mathrm{N2}}$, ${\beta}_{2}={A}_{2}/{\lambda}_{2}$. For the anterior flagellum, each segment ds_{1} contains additional $N\mathrm{d}{s}_{1}$ mastigonemes which are considered perpendicular to the segment itself. These mastigonemes also act as slender filaments experienced drag from water. Interestingly, due to the direction arrangement, the relative velocity normal to the flagellum segment results in drag force in tangent direction to the mastigonemes, and subsequently, the relative velocity tangent to the flagellum segment results in drag force in normal direction to the mastigonemes. Hence, the total drag force acting on a segment of the anterior flagellum is derived as:
where $N$ is the density of mastigonemes, $h$ is the length of each mastigoneme, ${V}_{\mathrm{N1}}$ and ${V}_{\mathrm{L1}}$ are two components of relative velocity of fluid in normal and tangent direction to the segment ds_{1}, ${K}_{\mathrm{Nf1}}$ and ${K}_{\mathrm{Lf1}}$ are the drag coefficients in normal and tangent direction of the flagellum filament, respectively; ${K}_{\mathrm{Nm1}}$ and ${K}_{\mathrm{Lm1}}$ are the drag coefficients in normal and tangent direction of the mastigonemes, respectively. ${K}_{\mathrm{Nf1}}$, ${K}_{\mathrm{Lf1}}$, ${K}_{\mathrm{Nm1}}$ and ${K}_{\mathrm{Lm1}}$ are also calculated from equation (16) and (17). In another perspective, we can consider the anterior flagellum receives additional drag from the mastigonemes, which is presented by two increased drag coefficients in normal and tangent direction defined as
and
respectively. Here, ${N}_{\mathrm{m}}$ is the density of mastigonemes; $h$ is the length of each mastigoneme; ${K}_{\mathrm{Nf1}}$ and ${K}_{\mathrm{Lf1}}$ are the drag coefficients in normal and tangent direction of the flagellum filament, respectively; ${K}_{\mathrm{Nm1}}$ and ${K}_{\mathrm{Lm1}}$ are the drag coefficients in normal and tangent direction of the mastigonemes, respectively. ${K}_{\mathrm{Nf1}}$, ${K}_{\mathrm{Lf1}}$, ${K}_{\mathrm{Nm1}}$ and ${K}_{\mathrm{Lm1}}$ are also calculated from equation (16) and (17). The total fluid drag force acting on the anterior flagellum is also derived from equation 13 and 18:
where ${v}_{\mathrm{w1}}={\lambda}_{1}{f}_{1}$ is the wave propagation velocity of the anterior flagellum, ${\gamma}_{1}={K}_{\mathrm{L1}}/{K}_{\mathrm{N1}}$, ${\beta}_{1}={A}_{1}/{\lambda}_{1}$.
At the same time, the ellipsoidal cell body moving with velocity ${U}_{X}$ also experiences a drag force from water (following Happel and Brenner)
where μ is water viscosity, ${\xi}_{e}$ is the shape coefficient of the ellipse cell body in 2D. ${\xi}_{e}$ is estimated as
with $\kappa =a/b$ being the ellipse body ratio. In low Reynolds number condition, total forces equate to zero due to approximately zero inertia. Thus,
Combine equation (19), (23) and (24), we achieve translational velocity ${U}_{X}$ of the zoospore as shown in Equation 5.
Appendix 6
Power and efficiency of two flagella during in zoospore swimming
We derive the power that each flagellum generates during zoospore swimming to investigate how the two flagella cooperate and contribute to the swimming. The power generated by the anterior flagellum P_{1} is calculated as
where ${\overrightarrow{v}}_{1}={\overrightarrow{v}}_{\mathrm{w}/\mathrm{f}(1)}$ is the relative velocity vector of water to each flagellar segment, and $\mathrm{d}{\overrightarrow{F}}_{1}$ is the drag force of water acting on each flagellar segment. From (13), we have
Then, we can derive P_{1} from (14), (15) and (28) as
Similarly, we can also derive P_{2} as
The useful power P_{0} required to propel the cell body with a speed ${U}_{X}$ can be derived as
The efficiency of the two flagella in propelling the cell body is estimated as
Appendix 7
Dimensions and physical parameters of zoospore’s flagella
With the characteristics of beating flagella obtained by kymograph, we calculate the values of drag coefficients ${K}_{\mathrm{N}}$, ${K}_{\mathrm{L}}$, flagellum shape ratio β, and drag coefficient ratio γ. The results are presented in the table below.
Appendix 8
Data publication
Our experimental data are accessible via Zenodo (https://doi.org/10.5281/zenodo.4710633).
In the data, we include
Datasets of all zoospore positions along multiple trajectories in the experiment of Figure 2,
A MATLAB file to compute all the statistical results in Figure 2(D–G),
A MATLAB file containing the simulation model presented in Figure 2(H),
Datasets of zoospore positions, speed, moving directions, body orientations during the turning, presented in Figure 4(A–D).
Data availability
All data generated and simulation files are available via Zenodo using this URL: https://doi.org/10.5281/zenodo.4710633https://doi.org/10.5281/zenodo.4710633. In the data, we include: (1) datasets of all zoospore positions along multiple trajectories in the experiment of Figure 2, (2) a MATLAB file to compute all the statistical results in Figure 2(DG), (3) a MATLAB file containing the simulation model presented in Figure 2(H), (4) datasets of zoospore positions, speed, moving directions, body orientations during the turning, presented in Figure 4(AD).

ZenodoCoordination of two opposite flagella allows highspeed swimming and active turning of individual zoospores.https://doi.org/10.5281/zenodo.4710633
References

Fibrous Flagellar Hairs of Chlamydomonas reinhardtii Do Not Enhance SwimmingBiophysical Journal 118:2914–2925.https://doi.org/10.1016/j.bpj.2020.05.003

Potassium homeostasis influences the locomotion and encystment of zoospores of plant pathogenic oomycetesFungal Genetics and Biology 42:213–223.https://doi.org/10.1016/j.fgb.2004.11.003

Phytophthora zoospores: From perception of environmental signals to inoculum formation on the hostroot surfaceComputational and Structural Biotechnology Journal 18:3766–3773.https://doi.org/10.1016/j.csbj.2020.10.045

Propulsive forces on the flagellum during locomotion of Chlamydomonas reinhardtiiBiophysical Journal 100:2716–2725.https://doi.org/10.1016/j.bpj.2011.05.001

A steering mechanism for phototaxis in ChlamydomonasJournal of the Royal Society, Interface 12:20141164.https://doi.org/10.1098/rsif.2014.1164

Locomotion of flagellates with mastigonemesJournal of Mechanochemistry & Cell Motility 3:207–217.

Fluid Mechanics of Propulsion by Cilia and FlagellaAnnual Review of Fluid Mechanics 9:339–398.https://doi.org/10.1146/annurev.fl.09.010177.002011

The eukaryotic tree of life from a global phylogenomic perspectiveCold Spring Harbor Perspectives in Biology 6:a016147.https://doi.org/10.1101/cshperspect.a016147

Motility, taxis and tropism in phytophthoraPhytophthora: Its Biology, Taxonomy, Ecology and Pathology 1:95–107.

On torque and tumbling in swimming Escherichia coliJournal of Bacteriology 189:1756–1764.https://doi.org/10.1128/JB.0150106

A Phylogenomic Framework to Study the Diversity and Evolution of Stramenopiles (=HeterokontsMolecular Biology and Evolution 33:2890–2898.https://doi.org/10.1093/molbev/msw168

Emerging oomycete threats to plants and animalsPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 371:1709.https://doi.org/10.1098/rstb.2015.0459

Economic Impact of Phytophthora Diseases in Southeast AsiaDiversity and Management of Phytophthora in Southeast Asia ACIAR Monograph 114:10–28.

Flagellar synchronization independent of hydrodynamic interactionsPhysical Review Letters 109:13.https://doi.org/10.1103/PhysRevLett.109.138102

Guidance of zoospores by potassium gradient sensing mediates aggregationJournal of the Royal Society, Interface 16:20190367.https://doi.org/10.1098/rsif.2019.0367

Bioinspired reorientation strategies for application in micro/nanorobotic controlJournal of MicroBio Robotics 16:173–197.https://doi.org/10.1007/s12213020001307

The Propulsion of SeaUrchin SpermatozoaJournal of Experimental Biology 32:802–814.https://doi.org/10.1242/jeb.32.4.802

A biased random walk model for the trajectories of swimming microorganismsJournal of Theoretical Biology 186:503–526.https://doi.org/10.1006/jtbi.1997.0421

Flagellar movement of human spermatozoaGamete Research 13:185–197.https://doi.org/10.1002/mrd.1120130302

The spores of Phytophthora: weapons of the plant destroyerNature Reviews. Microbiology 3:47–58.https://doi.org/10.1038/nrmicro1064

The Top 10 oomycete pathogens in molecular plant pathologyMolecular Plant Pathology 16:413–434.https://doi.org/10.1111/mpp.12190

Hydrodynamic mobility of chiral colloidal aggregatesPhysical Review. E, Statistical, Nonlinear, and Soft Matter Physics 79:051405.https://doi.org/10.1103/PhysRevE.79.051405

Natureinspired microfluidic propulsion using magnetic actuationPhysical Review. E, Statistical, Nonlinear, and Soft Matter Physics 79:046304.https://doi.org/10.1103/PhysRevE.79.046304

Theoretical modeling in microscale locomotionMicrofluidics and Nanofluidics 20:e1761.https://doi.org/10.1007/s104040161761y

The hydrodynamics of swimming microorganismsReports on Progress in Physics 72:e96601.https://doi.org/10.1088/00344885/72/9/096601

Chirality in microswimmer motion: From circle swimmers to active turbulenceThe European Physical Journal Special Topics 225:2319–2331.https://doi.org/10.1140/epjst/e2016600546

Analysis of a swimming sperm in a shear flowMicrofluidics and Nanofluidics 17:809–819.https://doi.org/10.1007/s1040401413715

Chlamydomonas flagellaJournal of Phycology 36:261–273.https://doi.org/10.1046/j.15298817.2000.99218.x

Characterisation of Phytophthora nicotianae zoospore and cyst membrane proteinsMycological Research 106:1211–1223.https://doi.org/10.1017/S0953756202006512

Tactic response of zoospores of the fungus Phytophthora palmivora to solutions of different pH in relation to plant infectionMicrobiology (Reading, England) 141:1231–1237.https://doi.org/10.1099/1350087214151231

Quantitative observations of flagellar motility of capacitating human spermatozoaHuman Reproduction (Oxford, England) 12:1006–1012.https://doi.org/10.1093/humrep/12.5.1006

Swimming direction reversal of flagella through ciliary motion of mastigonemesBiomicrofluidics 5:34108–3410815.https://doi.org/10.1063/1.3608240

Numerical modelling of chiralityinduced bidirectional swimming of artificial flagellaProceedings. Mathematical, Physical, and Engineering Sciences 470:20130547.https://doi.org/10.1098/rspa.2013.0547

Phytophthora nicotianae diseases worldwide: new knowledge of a longrecognised pathogenPhytopathologia Mediterranea 55:20–40.https://doi.org/10.14601/Phytopathol_Mediterr16423

Bacteria display optimal transport near surfacesNature Physics 15:610–615.https://doi.org/10.1038/s4156701904605

Chlamydomonas swims with two “gears” in a eukaryotic version of runandtumble locomotionScience (New York, N.Y.) 325:487–490.https://doi.org/10.1126/science.1172667

Swimming by reciprocal motion at low Reynolds numberNature Communications 5:5119.https://doi.org/10.1038/ncomms6119

Autoaggregation in zoospores of Phytophthora infestans: the cooperative roles of bioconvection and chemotaxisJournal of the Royal Society, Interface 11:20140017.https://doi.org/10.1098/rsif.2014.0017

Fiji: an opensource platform for biologicalimage analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019

Plant disease: A threat to global food securityAnnual Review of Phytopathology 43:83–116.https://doi.org/10.1146/annurev.phyto.43.113004.133839

TrackMate: An open and extensible platform for singleparticle trackingMethods (San Diego, Calif.) 115:80–90.https://doi.org/10.1016/j.ymeth.2016.09.016

Realtime imaging of fluorescent flagellar filamentsJournal of Bacteriology 182:2793–2801.https://doi.org/10.1128/JB.182.10.27932801.2000

Measurement of cell velocity distributions in populations of motile algaeThe Journal of Experimental Biology 207:1203–1216.https://doi.org/10.1242/jeb.00881

Zoospore development in the oomycetesFungal Biology Reviews 21:10–18.https://doi.org/10.1016/j.fbr.2007.02.001
Article and author information
Author details
Funding
French Government  Investissements d'Avenir  ANR (UCAJEDI ANR15IDEX01)
 Eric Galiana
 Xavier Noblin
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work has been supported by the French Government, managed by the National Research Agency under the Project UCA^{JEDI} "Investissements d'Avenir Investments" bearing the reference n° ANR15IDEX01. CCMA electron microscopy equipments have been funded by the Région Sud  ProvenceAlpesCôte d’Azur, the Conseil Général des Alpes Maritimes, and the GISIBiSA. We thank Ilaria Bassani, Nicolas Bruot, Xinhui Shen and NgocPhu Tran for fruitful discussions, as well as Marcos, David GonzalezRodriguez and the three reviewers for their careful reading and giving valuable feedback to improve the manuscript.
Version history
 Preprint posted: April 23, 2021 (view preprint)
 Received: June 12, 2021
 Accepted: March 25, 2022
 Accepted Manuscript published: March 28, 2022 (version 1)
 Accepted Manuscript updated: April 4, 2022 (version 2)
 Version of Record published: May 3, 2022 (version 3)
Copyright
© 2022, Tran et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics

 1,638
 Page views

 256
 Downloads

 3
 Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading

 Microbiology and Infectious Disease
A productive HIV1 infection in humans is often established by transmission and propagation of a single transmitted/founder (T/F) virus, which then evolves into a complex mixture of variants during the lifetime of infection. An effective HIV1 vaccine should elicit broad immune responses in order to block the entry of diverse T/F viruses. Currently, no such vaccine exists. An indepth study of escape variants emerging under host immune pressure during very early stages of infection might provide insights into such a HIV1 vaccine design. Here, in a rare longitudinal study involving HIV1 infected individuals just days after infection in the absence of antiretroviral therapy, we discovered a remarkable genetic shift that resulted in near complete disappearance of the original T/F virus and appearance of a variant with H173Y mutation in the variable V2 domain of the HIV1 envelope protein. This coincided with the disappearance of the first wave of strictly H173specific antibodies and emergence of a second wave of Y173specific antibodies with increased breadth. Structural analyses indicated conformational dynamism of the envelope protein which likely allowed selection of escape variants with a conformational switch in the V2 domain from an αhelix (H173) to a βstrand (Y173) and induction of broadly reactive antibody responses. This differential breadth due to a single mutational change was also recapitulated in a mouse model. Rationally designed combinatorial libraries containing 54 conformational variants of V2 domain around position 173 further demonstrated increased breadth of antibody responses elicited to diverse HIV1 envelope proteins. These results offer new insights into designing broadly effective HIV1 vaccines.

 Microbiology and Infectious Disease
For most retroviruses, including HIV, association with the plasma membrane (PM) promotes the assembly of immature particles, which occurs simultaneously with budding and maturation. In these viruses, maturation is initiated by oligomerization of polyprotein precursors. In contrast, several retroviruses, such as MasonPfizer monkey virus (MPMV), assemble in the cytoplasm into immature particles that are transported across the PM. Therefore, protease activation and specific cleavage must not occur until the preassembled particle interacts with the PM. This interaction is triggered by a bipartite signal consisting of a cluster of basic residues in the matrix (MA) domain of Gag polyprotein and a myristoyl moiety Nterminally attached to MA. Here, we provide evidence that myristoyl exposure from the MA core and its insertion into the PM occurs in MPMV. By a combination of experimental methods, we show that this results in a structural change at the Cterminus of MA allowing efficient cleavage of MA from the downstream region of Gag. This suggests that, in addition to the known effect of the myristoyl switch of HIV1 MA on the multimerization state of Gag and particle assembly, the myristoyl switch may have a regulatory role in initiating sequential cleavage of MPMV Gag in immature particles.