Rolling controls sperm navigation in response to the dynamic rheological properties of the environment
Abstract
Mammalian sperm rolling around their longitudinal axes is a longobserved component of motility, but its function in the fertilization process, and more specifically in sperm migration within the female reproductive tract, remains elusive. While investigating bovine sperm motion under simple shear flow and in a quiescent microfluidic reservoir and developing theoretical and computational models, we found that rolling regulates sperm navigation in response to the rheological properties of the sperm environment. In other words, rolling enables a sperm to swim progressively even if the flagellum beats asymmetrically. Therefore, a rolling sperm swims stably along the nearby walls (walldependent navigation) and efficiently upstream under an external fluid flow (rheotaxis). By contrast, an increase in ambient viscosity and viscoelasticity suppresses rolling, consequently, nonrolling sperm are less susceptible to nearby walls and external fluid flow and swim in twodimensional diffusive circular paths (surface exploration). This surface exploration mode of swimming is caused by the intrinsic asymmetry in flagellar beating such that the curvature of a sperm’s circular path is proportional to the level of asymmetry. We found that the suppression of rolling is reversible and occurs in sperm with lower asymmetry in their beating pattern at higher ambient viscosity and viscoelasticity. Consequently, the rolling component of motility may function as a regulatory tool allowing sperm to navigate according to the rheological properties of the functional region within the female reproductive tract.
Introduction
In mammals, sperm must migrate through the female reproductive tract to fertilize an egg (Suarez and Pacey, 2006; Suarez, 2016). During this migration, sperm require navigational mechanisms to swim in the correct direction (Eisenbach and Giojalas, 2006; Kaupp et al., 2008). These navigational mechanisms rely on external and dynamic biochemical and biophysical cues that are present in the female reproductive tract (Kaupp et al., 2008; Bahat et al., 2003; Tung et al., 2015a). Although the role of biochemical cues in mammalian sperm navigation remains poorly understood (Kaupp et al., 2008; Suarez, 2008), in vitro and in vivo studies have provided evidence for two navigational mechanisms that rely on external biophysical cues, namely rheotaxis (Miki and Clapham, 2013; Kantsler et al., 2014; Bukatin et al., 2015; Tung et al., 2015b) and walldependent navigation (Guidobaldi et al., 2014; Nosrati et al., 2016; Denissenko et al., 2012; Zaferani et al., 2019; Wang and Larina, 2018). Rheotaxis, as an upstream swimming in response to an external fluid flow, has been observed and quantitatively studied for bovine, human, and mouse sperm (Miki and Clapham, 2013; Kantsler et al., 2014). Walldependent navigation, as sperm response to the nearby physical boundaries such as walls of the female reproductive tract, has also been observed and characterized for bovine and human sperm (Tung et al., 2015a; Denissenko et al., 2012).
Although walldependent navigation combined with rheotaxis may characterize regulatory mechanisms of sperm navigation within the complex geometry of the female reproductive tract and under dynamic fluid flow, it remains unclear how sperm rolling contributes to these navigational mechanisms (Gadadhar et al., 2021; Miller et al., 2018; Drake, 1974; Babcock et al., 2014; Schiffer et al., 2020). Furthermore, the dynamic biophysical factors of the female reproductive tract are not limited to its complex geometry and the varying flow of the mucus within it. The rheological properties of the mucus also change according to the functional region within the female reproductive tract (Suarez and Pacey, 2006; Carlstedt and Sheehan, 1984; Tung et al., 2017). It also remains poorly understood how the rheology of the environment influences sperm swimming behavior, and in particular, rheotaxis and walldependent navigation.
To address these questions, we investigated bovine sperm motion in a microfluidic device to identify how the rolling component of its motility contributes to navigation under simple shear flow and within a quiescent reservoir. To avoid errors arising from studying sperm motility under external fluid flow, such as experimental inaccuracies caused by decoupling the effect of flow on sperm motion from active swimming, we employed a multistep approach. In this approach, we first isolated sperm within a quiescent reservoir using a rheotaxisbased method to ensure that sperm within the reservoir could swim upstream before entering the reservoir (Zaferani et al., 2018). We then characterized the components of sperm motility, including flagellar beating and rolling, in viscous and viscoelastic media, in the absence of external flow. Finally, we studied walldependent navigation within the reservoir and, by tracking the sperm prior to their entry into the reservoir, we evaluated rheotaxis to identify the function of rolling in sperm navigation under fluid flow.
We found that rolling enables sperm to swim progressively even when the sperm flagellar beating pattern is intrinsically asymmetric, which subsequently promotes rheotaxis and walldependent navigation. Sperm that lack rolling swim along twodimensional (2D) diffusive circular paths and are less susceptible to being influenced by the nearby walls and external flow. We observed that the suppression of rolling occurs by increasing ambient viscosity or viscoelasticity, such that an increase in ambient rheological properties transitions progressive motion into 2D diffusive circular surface exploration. Such diffusive circular surface exploration is caused by the intrinsic asymmetry of flagellation and the curvature of the circular path is proportional to the level of asymmetry. We noticed that suppression of rolling was reversible, as a decrease in ambient viscosity or viscoelasticity resulted in reactivation of rolling. Furthermore, we found out that the level of flagellar asymmetry in a sperm population forms a continuum, and suppression of rolling in sperm with lower flagellar asymmetry occurs at higher viscosity or viscoelasticity.
Sperm swimming behavior transitions between progressive and diffusive circular motions after each incidence of suppression or reactivation of rolling. Since the suppression or reactivation of rolling relies on changes in the viscosity or viscoelasticity of the media, sperm swimming behavior manifests differently in response to the rheological properties of the environment. Because the characteristics of these swimming behaviors (circular versus progressive) are different and the viscosity and viscoelasticity of the mucus in the female reproductive tract varies according to functional regions (Suarez and Pacey, 2006; Tung et al., 2017; Ishimoto and Gaffney, 2018; Ishimoto and Gaffney, 2016; Gaffney et al., 2011), our results suggest that rolling potentially enables sperm to regulate its navigation in response to the dynamic rheological properties of the mucus in the tract.
Results
To select sperm based on their rheotactic behavior and to isolate them inside the quiescent reservoir, we used a microfluidic corral system. Within the reservoir filled with standard Tyrode’s albumin lactate pyruvate medium (TALP), a subpopulation of sperm (<5%) exhibited an inplane 2D asymmetric flagellar beating pattern, where the midpiece of each flagellum was consistently bent more significantly to one side (Figure 1A and Video 1). This asymmetric beating pattern results in circular swimming at an angular velocity of $\mathrm{\Omega}$, whereas other sperm within the population swim progressively (Figure 1B). Unlike circular motion, progressive motion is not just produced by 2D inplane symmetric flagellation, rather frequent but irregular rollings contribute to sperm motility (Figure 1C and Video 2). The rolling component was detected when a change in the light intensity of the sperm heads was visualized under a phase contrast microscope (Figure 1C). Rollings occur rapidly and discontinuously (Drake, 1974; Babcock et al., 2014), such that the elapsed time between two consecutive rolling events $\left({T}_{SR}\right)$ is not constant. While tracking the sperm within the reservoir and under the flow prior to their entry into the reservoir, we found that progressive motility included the rolling component, both under the flow and within the quiescent reservoir, whereas sperm exhibiting a circular motion did not exhibit rolling under either condition (Figure 1D). Although rolling occurs independently of external fluid flow, it does depend on the rheological properties of the medium, as the percentage of rolling sperm decreases with an increase in the viscosity and viscoelasticity of the medium (Figure 1E, Video 3). Our results indicate that the suppression of rolling is more sensitive to viscoelasticity than to viscosity. We used varying concentrations of polyvinylpyrrolidone (PVP) and polyacrylamide (PAM) to prepare viscous and viscoelastic solutions (Tung et al., 2017), the rheological characteristics of which are presented in Section I of Appendix 1.
We observed that, in all the solutions we used, rolling sperm swam progressively, whereas nonrolling sperm exhibited 2D asymmetric flagellar beatings and swam in circles, such that the curvature of the circular path was found to be proportional to the asymmetry level of flagellar beating. Furthermore, an increase in viscosity decreased the frequency of rolling $\left({\omega}_{SR}=\raisebox{1ex}{$1$}\!\left/ \!\raisebox{1ex}{${T}_{SR}$}\right.\right)$ in rolling sperm (Figure 1F) as well as their progressive velocities, such that rolling sperm continue to swim progressively but slowly, with rolling occurring at lower frequency in a viscous solution. An increase in viscoelasticity did not, however, change the frequency of rolling and the average path velocity of rolling sperm (Figure 1F). We also observed that suppression of rolling in the viscous solution (4% PVP) occurred with a significant decrease in the propulsive velocity of nonrolling sperm, whereas rolling suppression in the viscoelastic solution (1% PAM) did not result in a decrease in the propulsive velocity of nonrolling sperm (Figure 1G). Considering that viscosity of 1% PAM solution is an order of magnitude higher than that of 4% PVP solution (Appendix I, Section I), these results seem to contradict our claim that an increase of viscosity leads to suppression of rolling as well as decrease in the propulsive velocity. However, we stress that the storage modulus $\left({G}^{\text{'}}\right)$ of 1% PAM solution which represents the elastic properties of the fluid was two orders of magnitude higher than that of 4% PVP solution (SI, Section I). The higher elasticity of 1% PAM indicates that this increase contributed significantly to the suppression of rolling, while it did not result in decrease of the sperm propulsive velocity.
We also noticed that the range of curvature $\left(\kappa \right)$ of the circular path in the nonrolling sperm population depended on ambient viscoelasticity (Figure 1H). As shown in Figure 1H, the suppression of rolling in sperm with a higher degree of asymmetry in their beating pattern occurred at lower ambient viscoelasticity, whereas higher ambient viscoelasticity was needed to suppress the rolling of sperm that exhibited lower asymmetry in their beating patterns. Similar behavior was observed when we increased the viscosity of the solution.
In TALP + 1% and 2% PVP solutions, some sperm exhibiting circular motion (<10%) also exhibited infrequent rollings, such that $\mathrm{\Omega}{T}_{SR}>2\pi $. Such infrequent rollings changed the direction of the circular motion (Figure 1I and Video 4) and resulted in abrupt relocations of the circular path’s center without significantly changing the curvature of the path.
Although our results indicate that rolling is a key contributor to progressive motility, more evidence is needed to validate the hypothesis. Therefore, we decreased the concentration of bovine serum albumin to 0.5% in TALP and tethered sperm heads to the glass surface upon their entry into the reservoir (Saggiorato et al., 2017). Tethering the sperm heads to the glass surface suppressed the rolling component, and the 2D flagellar beating pattern was observed separately from rolling. We observed that, without rolling, flagellation was not necessarily symmetric and >80% of the sperm began rotating around the sperm head upon tethering, although their motion was progressive prior to the surface binding (Figure 1J, Videos 5 and 6).
Another piece of evidence that validates our hypothesis was observed by tracking single sperm migrating between TALP and 1% PAM or 4% PVP solutions (Figure 1K). The representative sperm trajectories shown in Figure 1K indicate that suppression of rolling upon exiting TALP and entering into the viscous or viscoelastic solutions (blue trajectories) resulted in an instantaneous transition in sperm swimming behavior from progressive to circular. The transition in sperm swimming behavior was found to be reversible, as reactivation of rolling upon exiting the viscous or viscoelastic solutions and entering into TALP (red trajectories) caused sperm to start swimming progressively. This reversible transition was observed in more than 95% of sperm cells (total count = 109) migrating between TALP and 1% PAM or 4% PVP solutions. This finding agrees with the results shown in Figure 1G, as sperm velocity of average path for the blue and red trajectories was measured to be (70 ± 5) µms^{−1} in TALP, (45 ± 5) µms^{−1} in 4% PVP, and (70 ± 5) µms^{−1} in 1% PAM solutions.
We measured the angular velocities of all tethered sperm and compared them with that of untethered sperm, which were swimming freely in circles (Figure 1L). Consistent with Figure 1H, tethered sperm exhibited lower angular velocities than untethered sperm. Thus, rolling sperm exhibit lower asymmetry than their nonrolling counterparts and, accordingly, higher ambient viscosity and viscoelasticity is needed to suppress the rolling motion of sperm with lower asymmetry in flagellation and vice versa.
Rolling and progressive motion
To quantify the relationship between rolling and progressive motion, we first characterized asymmetry in the flagellar beating pattern. In agreement with Friedrich et al., 2010, measuring bending in the midpiece of the flagellum in time (Figure 2A) and the corresponding normalized fast Fourier transform $\left({P}^{*}\left(t\right)\right)$ revealed the presence of a zeroth harmonic within the frequency domain of flagellation (Figure 2B). Note that the constant offset in the normalized power spectrum is white noise in our measurement system. Furthermore, the experimental noise coming from our measurement system (e.g., image processing) is also included in the peak width around the first harmonic frequency. That is, the peak width centered at $\omega $ includes the intrinsic noise originated from flagellar sources, as well as the noise associated with our measurement system.
Writing the flagellar beating pattern in Fourier series ansatz (Equation 1) and further modeling the flagellar beating pattern (SI, Section II), we found that $\sigma =\frac{{a}_{1}{a}_{0}}{{a}_{1}}$ approximately presents the asymmetry level in the beating pattern. That is, $\sigma =1$ presents symmetric beating whereas lower values correspond to higher asymmetry in the beating pattern:
In Equation 1, the $x$ axis was set parallel to the flagellum at its straightline form, $k=\frac{2\pi}{\lambda}$ (with $\lambda \approx L$) is the wave number, $\omega $ is the main frequency, and ${a}_{n}$ is the amplitude of the nth harmonic. Note that the amplitude and frequencies of beating include deltacorrelated Gaussian noise that can be described as ${a}_{n}={\stackrel{~}{a}}_{n}\left(1+{\eta}_{n}\left(t\right)\right)$ and ${\omega}_{n}={\stackrel{~}{\omega}}_{n}\left(1+\xi \left(t\right)\right)$, respectively, with $\u27e8{\eta}_{n}\left(t\right)\u27e9=\u27e8\xi \left(t\right)\u27e9=0$, $\u27e8{\eta}_{n}\left(t\right){\eta}_{{n}^{\text{'}}}\left(t\text{'}\right)\u27e9={D}_{n}{\delta}_{n{n}^{\text{'}}}\delta (tt\text{'})$, and $\u27e8\xi \left(t\right)\xi \left(t\text{'}\right)\u27e9={D}_{\omega}\delta (tt\text{'})$. The tilde sign above a quantity represents its time average value. Following Saggiorato et al., 2017; Friedrich et al., 2010; Elgeti et al., 2015, we applied small amplitude (i.e., $\forall n\to {a}_{n}\ll L$) and length preservation constraints and used resistive force theory (SI, Sections II–IV), to find that the zeroth harmonic in the beating pattern yields a toque as follows:
Note that ${\xi}_{T}$ and ${\xi}_{N}$ are drag coefficients in the tangential and normal directions, respectively, and $L$ is the sperm length. Although the propulsive force produced by the flagellum correlates with the characteristics of the first and higher harmonics, the amplitude of the zeroth harmonic is involved solely in the torque produced (Equation 2). Applying the zero net torque and force constraints, we found that propulsive $\left({\stackrel{~}{V}}_{P}\right)$ and angular $\left(\stackrel{~}{\mathrm{\Omega}}\right)$ velocities were related through ${a}_{0}$:
To verify Equation 3, we measured the angular velocity of the sperm and plotted a normalized path curvature $\left(L\mathrm{\Omega}{V}_{P}^{1}\right)$ with respect to the normalized amplitude of the zeroth harmonic $\left({a}_{0}{L}^{1}\right)$ in TALP, TALP + 4% PVP (viscous), and TALP + 1% PAM (viscoelastic) solutions (Figure 2C and D). The linear correlation between $L\mathrm{\Omega}{V}_{P}^{1}$ and ${a}_{0}{L}^{1}$ is consistent with our mathematical arguments (Equations 1–3), confirming that ${a}_{0}$ modulates the curvature of the nonrolling sperm path. Furthermore, the linear relationship between $L\mathrm{\Omega}{V}_{P}^{1}$ and ${a}_{0}{L}^{1}$ predicted by the resistive force theory is preserved for sperm motion in a viscoelastic solution whereas the slope differs from that observed for sperm motion in standard and viscous solutions (Zhang and Goldman, 2014; Teran et al., 2010; Li and Ardekani, 2015; Riley and Lauga, 2017; Spagnolie et al., 2013).
The resistive force theory, as a mean field approach, cannot explain the observed inconsistency in the circular path. The inconsistency in the circular path, however, can be characterized by quantifying the random fluctuations at the center of the circular path (Ma et al., 2014). We noticed that fluctuations at the center, and thus circular motion, are diffusive in character as the mean square displacement (MSD) of the center is proportional to the elapsed time: $\mathrm{M}\mathrm{S}\mathrm{D}~t$ (Figure 2E). This diffusive motion can be quantified by the diffusion coefficient of the center determined by the intercept of MSD, or alternatively, the center’s speed distribution (Figure 2F).
We measured the head light intensity (HLI) of the sperm over time to characterize rollings during progressive motility (Figure 2G). HLI is a pulsetype quantity, with the pulse duration shorter than the time that elapses between two consecutive pulses, ${T}_{SR}$. Therefore, we defined the edgesensitive function $\mathrm{\Pi}\left(t\right)$ with $\mathrm{\Pi}\left(0\right)=1$, such that $\mathrm{\Pi}\left(t\right)$ is multiplied by −1 at each positive edge of HLI (Figure 2G). Note that $\mathrm{\Pi}\left(t\right)$ captures the rolling component as a rapid switch in the direction of asymmetry (Figure 1I); rolling can therefore be incorporated into Equation 1 with $\mathrm{\Pi}\left(t\right)$:
Solving the equations of motion using Equation 4, we found that, depending on $\stackrel{~}{\mathrm{\Omega}{T}_{SR}}$, sperm swim along varying pathways with average progressive velocity of:
(Figure 2—figure supplement 4A and B). At the frequent rolling limit $\left(\stackrel{~}{\mathrm{\Omega}{T}_{SR}}\to 0\right)$, $\stackrel{}{V}$ approaches ${\stackrel{~}{V}}_{P}$, which means that frequent rolling asymptotically yields progressive motion with the average path velocity equal to the propulsive velocity even though $\stackrel{~}{\mathrm{\Omega}}\ne 0$. Our experimental measurements are consistent with this finding because $\stackrel{~}{\mathrm{\Omega}{T}_{SR}}$ is less than $2\pi $ for rolling sperm, whereas for circular motion with infrequent rolling it is greater than $2\pi $ (Figure 2—figure supplement 4C). Note that for nonrolling sperm, $\stackrel{~}{\mathrm{\Omega}{T}_{SR}}\to \mathrm{\infty}$.
Even though $\stackrel{~}{\mathrm{\Omega}{T}_{SR}}$ determines the average path, how do deviations of ${T}_{SR}$ from the mean during motion influence the trajectory? Our measurements of ${T}_{SR}$ for one rolling sperm for ~40 s indicate that the mean and standard deviation of this random variable are $\stackrel{~}{T}}_{SR}\approx 0.11\text{}\mathrm{s},{\sigma}_{SR}\approx 0.02\text{}\mathrm{s$, respectively (Figure 2H). Combining arbitrary values for $\stackrel{~}{\mathrm{\Omega}}$ with ${T}_{SR}$ distribution, we found that the direction of the average path obeys a similar distribution with mean of 0 and standard deviation of ${\sqrt{{n}_{SR}}\sigma}_{SR}\stackrel{~}{\mathrm{\Omega}}$ (Figure 2I), in which ${n}_{SR}$ is the number of rolling occurrences.
Considering rolling as a switch in the direction of asymmetry, we investigated how rolling influenced sperm rheotaxis. The spermrheotactic behavior and angular velocity of upstream orientation is modeled using an Adlertype Equation 10:
in which $\gamma $ is the shear rate, $A$ is a constant, and $\theta $ is the relative orientation of the sperm with respect to the external fluid stream. For sperm with intrinsic angular velocity $\left(\mathrm{\Omega}\right)$, the net angular velocity under fluid flow is ${\Omega}_{RH}+\mathrm{\Omega}$, which orients the sperm with respect to the flow during upstream motion $\left({\theta}_{UP}\right)$:
Including the rolling component in Equation 7 using $\Pi \left(t\right)$, the average orientation of the sperm with respect to flow during upstream motion is:
Notably, $\stackrel{~}{\mathrm{\Pi}\left(t\right)}\ll 1$. Equation 8 predicts that rollings result in significant decay in $\stackrel{~}{\theta}}_{UP$. Consequently, Equations 9 and 10 respectively present sperm net velocity in the upstream direction (${\stackrel{~}{V}}_{UP}$) in the presence or absence of frequent rollings:
In Equations 9 and 10, ${V}_{P}$ and ${V}_{N}$ are sperm propulsive and perpendicular velocities (corresponding to the angular velocity), respectively, whereas ${V}_{F}$ is the external flow velocity (SI, Section V). To experimentally confirm the prediction obtained from Equations 9 and 10, we backtracked and analyzed sperm motion under flow prior to their entry into the quiescent zone. We observed that $\stackrel{~}{\theta}}_{UP$ is greater for nonrolling sperm $\left({\stackrel{~}{\theta}}_{UP}=40\xb0\pm 10\xb0\right)$ than for those exhibiting the rolling motion $\left({\stackrel{~}{\theta}}_{UP}=10\xb0\pm 5\xb0\right)$. Furthermore, the average upstream velocity of rolling sperm was much higher $\left({\stackrel{~}{V}}_{UP}=60\pm 10\text{}\mu \mathrm{m}/\mathrm{s}\right)$ than that of nonrolling sperm $\left({\stackrel{~}{V}}_{UP}=30\pm 10\text{}\mu \mathrm{m}/\mathrm{s}\right)$ (SI, Section V). These results indicate that, even though the rolling component is not required for sperm rheotaxis, it facilitates rheotactic behavior by minimizing the angle between sperm orientation and the external fluid flow, thus maximizing the upstream component of their motion.
Sperm surface exploration
Nonrolling sperm swim along diffusive circular paths and explore the surface (Ma et al., 2014). To characterize the diffusivity in circular motion (Figure 2E and F), we overlayed consecutive images of one sperm with 0.08 s intervals over different times (Figure 3A). The thickness of the combined circular paths $\left({\delta}_{14}\right)$ increases over time in all directions, ${\delta}_{14}\propto \sqrt{t}$ (Figure 3B) as predicted by the MSD of the circular path’s center (Figure 2E and F). For any sperm that diffused more rapidly, cell tracking followed by circle fitting was used to characterize the motion of circular path’s center (Figure 3C).
To identify the flagellar source of the center’s diffusion, we solved equations of motion where the amplitude and phase of all harmonics included white Gaussian noises, as described in Equation 1 (the results are shown in Figure 3D). Our simulations suggest that noise in the amplitudes and the phases of first and higher harmonics yielded a noisy ${V}_{P}$ without resulting in fluctuations of the center (SI, Section VI). By contrast, the noise in the amplitude of the zeroth harmonic yielded a similar noise in $\mathrm{\Omega}$ but not in ${V}_{P}$. Therefore, $\mathrm{\Omega}$ and $\kappa $ can be expressed as $\mathrm{\Omega}\left(t\right)=\stackrel{~}{\mathrm{\Omega}}(1+{\eta}_{0}(t\left)\right)$ and $\kappa \left(t\right)=\stackrel{~}{\kappa}(1+{\eta}_{0}(t\left)\right)$, respectively, such that lower signaltonoise ratios $\left(\mathrm{S}\mathrm{N}\mathrm{R}=\frac{{\stackrel{~}{\kappa}}^{2}}{\u27e8{\kappa}^{2}\u27e9{\stackrel{~}{\kappa}}^{2}}\right)$ for identical mean curvatures yielded faster diffusion, whereas higher SNRs yielded more consistent pathways (Figure 3D). We measured the normalized mean step size $\u27e8{\delta}^{2}{\kappa}^{2}\u27e9$ and the corresponding SNR and found that our experimental data were consistent with the simulationbased results, except for a constant as shown in Figure 3E and F. That is, as suggested by simulations, the normalized mean step size was found to be inversely correlated with the SNRs: $\u27e8{\delta}^{2}{\kappa}^{2}\u27e9\propto {\mathrm{S}\mathrm{N}\mathrm{R}}^{1}$.
As previously published (Ma et al., 2014; Goldstein et al., 2009), these fluctuations in the circular path’s center correspond to nonthermal noise in the sperm flagellar beating, rather than thermally driven fluctuations. Approximating sperm as a rod with length of ∼80 µm and diameter of ∼5 µm, the rotational diffusion caused by thermal fluctuations is in the other of 10^{−6}s^{−1}, which yields a $SNR~{10}^{7}$ for $\stackrel{~}{\kappa}~(5\times {10}^{3}){\mu m}^{1}$. This order of $SNR$ estimated for thermal fluctuations are far greater than our measured values and the experimental resolution; therefore, the contribution of thermal noise to the sperm diffusive circular motion is negligible.
Diffusive circular motion might be abruptly interrupted by infrequent rollings (Figure 1I, Figure 3G and H). In this case, the center of the circular path not only diffuses in time but also relocates ballistically upon rolling in a random direction with ${V}_{r}$ (Figure 3H,I), such that $\frac{D}{{V}_{r}}$∼10^{−6} −1 µm (SI, Section VII). This twophase motion is an intermittent search, with greater efficiency in surface exploration than with normal diffusion (Bénichou et al., 2011). This increase in the efficiency of exploration may be interpreted as a decrease in the probability of revisiting previously swept spots. More precisely, the average area swept by the sperm through diffusion is proportional to the square root of diffusion time $\u27e8A\left(t\right)\u27e9\propto \sqrt{t}$. Assuming that ${T}_{i}$ is the time frame between the (i−1)th and ith relocations, the areas swept with and without relocation are proportional to $\sum _{i}\sqrt{{T}_{i}}$ and $\sqrt{\sum _{i}{T}_{i}}$, respectively. Because $\sqrt{\sum _{i}{T}_{i}}\mathrm{}\le \sum _{i}\sqrt{{T}_{i}}$, infrequent rolling increases the area swept by the sperm (see SI, Section VII).
Rolling and walldependent navigation
In addition to active swimming and external flow, the surrounding walls contribute to sperm motion (Elgeti et al., 2010). For a sperm that is positioned away from the walls (distance from the wall > the sperm length), the wall’s effect on the swimmer is known to be a drift velocity that can be either attractive or repulsive depending on the swimmer’s angle with respect to the wall (Elgeti et al., 2010). The model proposed for studying this drift velocity is based on positing the swimmer as a force dipole that includes the propulsive force provided by the flagellum as well as the corresponding drag force (Elgeti et al., 2010). However, this model can be used when the sperm exhibits fully symmetric flagellation. We therefore developed a new swimmer model that includes the torque caused by asymmetric beating. The proposed model is depicted in Figure 4—figure supplement 1A, in which $f$ is the propulsive force, ${f}^{\text{'}\text{'}}$ is the perpendicular force corresponding to the torque, and ${f}^{\text{'}}$ is the drag force required for the torquefree condition. We then carried out finite element method simulations in a cylindrical domain, like our microfluidic quiescent reservoir, and solved the Stokes and mass conservation equations for our model to find the velocity field imposed by the sperm active swimming.
The velocity field imposed by flagellar beating for $\stackrel{~}{\mathrm{\Delta}\theta}=0{15}^{\circ}$ is shown in Figure 4—figure supplement 1B. Integrating net flow in the y direction imposed on the sperm body that is caused by noslip walls, we calculated the drift velocity toward the wall (Figure 4—figure supplement 1C). The simulation results indicated that the motion of the microswimmer is influenced to a lesser degree by nearby boundaries as the circular component emerges in the motility (SI, Section VIII). This decay in the drift velocity caused by the circular component of motion was also predicted by the analytical solution derived from a Stokeslet (Ardekani and Stocker, 2010; Li and Ardekani, 2014) description (SI, Section IX). The attraction of the sperm toward the wall in the presence of circular behavior can be described by the following equation:
where ${U}_{w}$ and ${U}_{w}^{p}$ are the drift velocities with and without the circular motion. Calculating the average farfield drift velocity imposed on the sperm during one round of circulation, we found that the average magnitude of the drift velocity in one round (<1 µm/s) is much smaller than random fluctuations ((>540 µm/s) of the circular path’s center (Figure 4—figure supplement 1D and SI, Section X); therefore, the sperm’s circular motion is more strictly controlled by the level of asymmetry in its flagellum and diffusivity than by distant walls.
Frequent rollings, however, negate the contribution of circular components of the motion, producing a drift velocity as if the flagellation is symmetric and the sperm behaves like a dipole swimmer,
We reiterate that $\stackrel{~}{\mathrm{\Pi}\left(\mathrm{t}\right)}\ll 1$. Equation 12 states that frequent rollings make the sperm more readily disposed to physical boundaries at the far field.
Sperm nearfield interactions with its nearby walls fit in one of the four categories shown in Figure 4A. A rolling sperm reorients upon wall contact and swims along the wall (${S}_{1}$, Video 7), whereas with a nonrolling sperm, depending on $\mathrm{\Omega}$ and the location of the circular path’s center at the contact point, three other behaviors were observed. At a high magnitude of $\mathrm{\Omega}$, the sperm do not contact the wall, maintaining their circular motion (${S}_{2}$, Video 8). At lower magnitudes of $\mathrm{\Omega}$, the sperm contact the wall and, depending on the location of the path’s center relative to the contact point, either swim near the wall temporarily and detach (${S}_{3}$, Video 9) or swim more slowly along the wall (compared with ${S}_{1}$) with a tilted orientation with respect to the wall (${S}_{4}$, Video 10). We compared these four types of sperm–wall interactions quantitatively using the time of sperm detention on the wall (${T}_{D}$, Figure 4B), sperm velocity on the wall divided by its velocity before wall contact (${V}^{*}$, Figure 4C), and the normalized length of detention ${\rho}^{*}$ (Figure 4D,E). Based on the measurements shown in Figure 4B–E, we found that the ${S}_{1}$ and ${S}_{3}$ categories yield similar ${V}^{*}$ values (close to 1), whereas for ${S}_{4}$, ${V}^{*}$ is smaller than 1 and for $\mathrm{\Omega}>1.5\text{}{\mathrm{s}}^{1}$ the sperm do not migrate along the wall. Furthermore, ${S}_{1}$ and ${S}_{4}$ yield similar ${\rho}^{*}$ values, close to 1, whereas for ${S}_{3}$, ${\rho}^{*}$ is much smaller than 1, indicating temporary detention on the wall. A simple yet insightful approach to understanding these four categories involves surface contact force analysis (Marion, 2013), where the sperm contact with the wall can be modeled by a positive force that is perpendicular to the surface $N$.
Suppose that a sperm swims along a wall at an angle $\beta $ (Figure 4F). Under a zero net force constraint, the normal surface force becomes $N={F}_{T}\mathrm{sin}\left(\beta \right){F}_{N}\mathrm{cos}\left(\beta \right)$, where ${F}_{T}$ and ${F}_{N}$ are the propulsive and perpendicular forces produced by the sperm, respectively. The threshold angle (${\beta}_{th}$) that corresponds to the $N=0$ situation is equal to ${\mathrm{tan}}^{1}\gamma $, where $\gamma =\frac{{F}_{N}}{{F}_{T}}$. For $\beta <{\beta}_{th}$, $N$ becomes negative and no contact occurs (${S}_{2}$). Because an increase in $\gamma $ leads to higher ${\beta}_{th}$, sperm with greater $\gamma $ values are less likely to contact and follow the wall. For $\beta $ values greater than ${\beta}_{th}$, where sperm–wall contact occurs, greater ${F}_{N}$ yields a smaller $N$ and leads to easier detachment from the wall (${S}_{3}$). When the direction of the sperm perpendicular force at the contact point is against the wall $\left(N={F}_{T}\mathrm{sin}\left(\beta \right)+{F}_{N}\mathrm{cos}\left({\beta}_{s}\right)\right)$, a greater surface force is exerted, yielding a stronger attachment to the wall (Figure 4G). In this configuration, however, the paralleltothewall components of the perpendicular and propulsive forces thwart each other and cause slower sperm motion along the wall. Therefore, asymmetric flagellation perturbs sperm motion along the wall by decreasing either its detention time (${S}_{3}$) or velocity (${S}_{4}$) on the wall. However, the rolling component functions as a switch between ${S}_{3}$ and ${S}_{4}$, and not only guarantees longer detention but also increases sperm velocity along the wall (Figure 4H,I).
Based on previous theoretical and computational studies on hydrodynamic interactions of a sperm (or a bacteria) with a wall (Drescher et al., 2011; Schaar et al., 2015; Elgeti and Gompper, 2013; Rode et al., 2019; Elgeti and Gompper, 2016), we developed a simple hydrodynamic model of sperm–wall interaction at the lubrication limit. Solving Stokes and mass conservation equations (SI, Section XI), we plotted the phase curve $\left(\dot{\beta}\mathrm{v}\mathrm{s}\beta \right)$ of sperm dynamics after contacting the wall (Figure 4J). Note that $\dot{\beta}$ is the effect of the wall superimposed with intrinsic $\mathrm{\Omega}$. For $\mathrm{\Omega}=0$, stability ($\dot{\beta}=0$) occurs at ${\beta}_{S}=0$, implying that the final orientation of a symmetrically beating sperm with respect to the wall is 0. For $\mathrm{\Omega}>0$, corresponding to ${S}_{4}$, ${\beta}_{S}>0$ is the final angle between the sperm and the wall. For $\mathrm{\Omega}<0$, corresponding to ${S}_{3}$, ${\beta}_{S}<0$ and $\frac{d\beta}{dt}>0$ at $\beta =0$, suggesting that swimming parallel to the wall is unstable, and the sperm detaches from the wall with ${\beta}_{S}$. Figure 4J shows that the instability that occurs at ${\beta}_{U}=\frac{\pi}{2}$ (for $\mathrm{\Omega}=0$) changes with $\mathrm{\Omega}$ as well.
To include the contribution of rolling in sperm dynamics after wall contact, we modeled rollings as transitions between two curves in the phase space with positive and negative values of $\mathrm{\Omega}$ (Figure 4K). We then can rewrite $\dot{\beta}$ and include $\mathrm{\Pi}\left(t\right)$ such that:
Note that $g\left(\beta \right)$ is the curve in the phase space that corresponds to $\mathrm{\Omega}=0$. Insofar as stability occurs at $\dot{\beta}=0$,
in which ${{\beta}_{S}}^{*}$ and ${\beta}_{S}$ are the stable points with and without taking rolling into account. Because $\stackrel{~}{\mathrm{\Pi}\left(\mathrm{t}\right)}\ll 1$, the average of ${\beta}_{S}^{*}\left(t\right)$ is:
Equation 15 suggests that, at the frequent rolling limit where $\stackrel{~}{\mathrm{\Pi}\left(\mathrm{t}\right)}$ approaches 0, $\stackrel{~}{{\beta}_{S}^{\ast}}$ approaches 0 as well. Consequently, frequent rollings mitigate the destructive role of $\mathrm{\Omega}$ in sperm motion along the wall, thereby yielding faster and longerlasting swimming along the wall (SI, Section XII). Because at the frequent rolling limit $\stackrel{~}{{\beta}_{S}^{\ast}}$ is close to 0, we posit that $g\left(\stackrel{~}{{\beta}_{S}^{\ast}}\right)\approx 0$, and thus $\dot{\beta}$ near the stable point can be written as:
Based on Equation 16, $\beta \left(t\right)$ is a triangular function near the stable point in the form of:
We measured $\dot{\beta}$ experimentally to verify the results obtained from the lubrication theory (Figure 4J). The net rotation ($\dot{\beta}$) with respect to the angle between the sperm and wall ($\beta $) that corresponds to ${S}_{1}$ and ${S}_{4}$ is consistent with the results obtained from our model (Figure 4J,L,M). Furthermore, the frequent rollings observed in the ${S}_{1}$ category decreased $\stackrel{~}{{\beta}_{S}}$ significantly to ∼10° (Figure 4N) in comparison with $\stackrel{~}{{\beta}_{S}}$ for ${S}_{4}$, which was greater and linearly correlated with $\frac{\mathrm{\Omega}}{V}$ (Figure 4N, inset plot). This measurement supports our prediction that rolling enables the sperm to swim stably along the wall by decreasing $\stackrel{~}{{\beta}_{s}}$. Measuring $\beta \left(t\right)$ during stable swimming along the wall for 25 sperm that belong to the ${S}_{1}$ category, we found that $\beta \left(t\right)$ is indeed a triangular function near the stable point, as predicted by Equation 17 (SI, Section XII).
A unified picture of sperm motion within the quiescent reservoir can be obtained by developing a state diagram and identifying the transitions between the ${S}_{14}$ states, possibly through the diffusivity of the circular motion and rolling. Suppose that, within the quiescent reservoir (radius of $R$), the nonrolling sperm swims in circle (with a radius of $R}^{\prime$), such that the distance between the centers of the two circles is $d\left(t\right)$, which evolves through diffusion (Figure 5A). Defining $s\left(t\right)=\frac{1}{2}\left({R}^{2}+{R}^{\prime 2}d{\left(t\right)}^{2}\right)$, the circular path does not intersect the reservoir for $s\left(t\right)>RR\text{'}$ (category ${S}_{2}$), whereas for $0<s\left(t\right)\le R{R}^{\prime}$ the angle between the two circles at the intersection point (sperm incidence angle) is greater than ${\beta}_{u}$ (category ${S}_{3}$); for $s<0$, the angle at the intersection point is less than ${\beta}_{u}$ (category ${S}_{4}$). Assuming that at $t=0$ and the sperm is in the ${S}_{2}$ state, the circular path starts to diffuse over time, which leads to transitions between the ${S}_{2}$ and ${S}_{3}$ states. The average time for the first occurrence of a transition is $\u27e8T\u27e9=\frac{{\left(R{R}^{\prime}\right)}^{2}}{2D}$ (Figure 5B). One might expect the diffusion process to continue until $s\left(t\right)$ becomes negative and the transition to ${S}_{4}$ occurs (Figure 5B). However, the angle at which the sperm detaches from the wall after the first contact $\left({\beta}_{s}\right)$ depends on $\mathrm{\Omega}$ rather than the location and angle at which contact occurs (Figure 5C). Therefore, after first contact, the sperm returns to the wall at an incidence angle that is identical to the angle at which it detaches from the wall ${\beta}_{s})$. Therefore, the transition from ${S}_{3}$ to ${S}_{4}$ occurs if ${\beta}_{s}>{\beta}_{u}$. Given that $\beta}_{s}\approx \frac{\mathrm{\Omega}}{\left{g}^{\prime}\left(0\right)\right$, ${\beta}_{u}\approx \frac{\pi}{2}+\frac{\mathrm{\Omega}}{{g}^{\text{'}}\left(\frac{\pi}{2}\right)}$, and $\left{g}^{\prime}(0)\right\ll {g}^{\prime}\left(\frac{\pi}{2}\right)$, the condition of such a transition reduces to:
The prime indicates derivative with respect to $\beta $. The corresponding curvature that satisfies Equation 18 is κ ∼ 0.1 µm^{−1}, which is much greater than the experimentally observed values. The state diagram of nonrolling sperm–wall interactions and the possible transitions are summarized in Figure 5D.
For rolling sperm, the state diagram shown in Figure 5D alters significantly. First, infrequent rollings result in abrupt changes in $d\left(t\right)$ and thus $s\left(t\right)$, causing reversible transitions between ${S}_{2}$, ${S}_{3}$, and ${S}_{4}$ (Figure 5F). Second, at the frequent rollings limit, one may write $s\left(t\right)\approx \mathrm{\Pi}\left(t\right)\lefts\right$, so that $\stackrel{~}{s\left(t\right)}$ approaches 0 (Figure 5E). Consequently, in the presence of frequent rollings, all states transition to ${S}_{1}$, where the sperm swims stably along the wall (Figure 5F).
Discussion
In summary, rolling is a component of mammalian sperm motility that is sensitive to ambient viscosity and viscoelasticity. In a solution with low viscosity and viscoelasticity, most sperm exhibit a rolling and progressive motion that is susceptible to external fluid flow (rheotaxis) and rigid physical boundaries (walldependent navigation). As ambient viscosity or viscoelasticity increases, the rolling component becomes suppressed and subsequently the sperm swim in diffusive circular paths (surface exploration), a type of motion that is less susceptible to being influenced by external fluid flow or nearby walls. Suppression of rolling was found to be reversible, as sperm migrate into medium with low viscosity or viscoelasticity, rolling reactivates, and thus sperm swimming transitions from surface exploration to progressive motion. Furthermore, we demonstrated that the suppression of rolling in sperm with lower asymmetry in their flagellar beating pattern occurs at higher viscosity or viscoelasticity. Therefore, the suppression of rolling, and thus the onset of surface exploration, depends on both ambient rheological properties and the level of asymmetry in sperm flagellation.
Our results evidenced sperm flagellar beating is intrinsically asymmetric, but frequent rolling counteracts this asymmetric flagellation by alternating the direction of asymmetry, and results in a progressive motion. This progressive motion resulted from frequent rollings found to be key to the sperm rheotaxis and walldependent navigations. But why is susceptibility to external fluid flow and nearby walls under dynamic conditions not controlled solely by the level of asymmetry in flagellation so that maximum susceptibility appears with the fully symmetric beating pattern?
We argue that while fully symmetric beating pattern yields an efficient walldependent navigation, it does not result in an efficient rheotaxis, because the tilted orientation of the sperm (caused by rolling) with respect to the surface is needed for efficient rheotaxis, as demonstrated previously (Kantsler et al., 2014). More importantly, it is possible that, unlike the asymmetric beating pattern that is modulated by chemical factors (Suarez, 2008; Quill et al., 2003; RamírezGómez et al., 2020), rolling depends on the ambient viscosity and viscoelasticity of the medium alone (Schiffer et al., 2020). Because the viscosity and viscoelasticity of the fluid within the female reproductive tract varies across functional regions, the tract possibly regulates sperm navigation by independently controlling: (1) the rolling component through regulating the rheology of the environment and (2) the asymmetry level in the flagellar beating pattern through secreting chemical factors. To show how regulating of asymmetrical beating caused by chemical factors works together with rolling, which is impacted by ambient viscosity and viscoelasticity, more studies regulating both aspects are needed.
Our results also demonstrate that the transition found between fast progressive and slow diffusive circular motions is reversible and occurs through suppression or reactivation of rolling. This finding, in particular, suggests that sperm motion during migration within the female reproductive tract is possibly bimodal. The fast progressive mode is an appropriate swimming behavior for sperm to migrate long distances between different functional regions within the female reproductive tract. This mode of motion may be regulated by the tract through rheotaxis and walldependent navigation. Whereas the slow diffusive motion is an appropriate swimming behavior for exploring the functional regions within the female reproductive tract to possibly receive physiological signals from these regions that are essential for the fertilization process. These physiological signals may include specific ligands secreted by the tract (Chang and Suarez, 2010), or pH of the functional region (Marquez and Suarez, 2007).
Our findings also suggest that elastic properties of the swimming media are key to the suppression of rolling and thus sperm motion. Characterizing the rheological properties of our viscous (4% PVP) and viscoelastic (1% PAM) solutions, we noticed that viscosity of 1% PAM is two orders of magnitude greater than that of 4% PVP solution. Rolling suppression in 1% PAM occurred without a loss in sperm propulsive velocity, whereas rolling suppression in the 4% PVP occurred with a significant loss in the propulsive velocity. Since the storage modulus of 1% PAM solution was two orders of magnitude greater than that of 4% PVP, we, therefore, conclude that the absence of a loss in the sperm propulsive velocity was due to the elastic properties of the solution. Accordingly, the elasticity of the swimming media is a key contributing factor to the sperm navigation within the female reproductive tract.
Our findings revealed the potential role of sperm rolling in mammalian fertilization, which results in a holistic and fundamental understanding of the fertilization process. Furthermore, these results are useful for more practical purposes such as designing technologies to improve fertility in cattle industry, as well as diagnosing and treating human male infertility. Furthermore, our results can be used to design new types of synthetic microswimmers that are responsive to dynamic physical environments, and thus more sensitive while exploring confined spaces.
Materials and methods
Sperm sample and culture media preparation
Request a detailed protocolCommercially available cryopreserved bovine semen samples taken from two mature black and white Holstein bulls (5.5 and 6 years of age) were kindly donated by Genex Cooperative (Ithaca, NY) in milk and egg yolk–based extender in plastic straws. The ejaculate concentration was 2.9 and 3 billion cells/mL, respectively, and had a prefreeze motility of 65%. The semen was thawed at 37°C in a water bath and diluted in a 1:4 ratio with TALP. After dilution, the viscosity of the samples was ~5 mPas. The initial sperm concentration in the thawed semen samples was ~200 million/mL, which was diluted to ~40 million/mL with TALP. The motility of the semen sample after dilution decreased to 20–30%. We used 10 separate semen samples in both the milk and egg yolk–based extender. At least three replicates were performed for each experiment to validate the accuracy of data and obtain a valid value for error bars.
TALP was prepared as follows: NaCl (110 mM), KCl (2.68 mM), NaH_{2}PO_{4} (0.36 mM), NaHCO_{3} (25 mM), MgCl_{2} (0.49 mM), CaCl_{2} (2.4 mM), HEPES buffer (25 mM), glucose (5.56 mM), pyruvic acid (1.0 mM), penicillin G (0.006% or 3 mg/500 mL), and bovine serum albumin (20 mg/mL). To tether the sperm head to the glass surface, we reduced the concentration of bovine serum albumin to 5 mg/mL. To increase the viscosity and viscoelasticity of TALP, we added 1–4% of PVP (weight percent) and 0.25–1% PAM (weight percent).
Rheological measurements
Request a detailed protocolDynamic rheological measurements were performed with a rheometer (MCR 501, Anton Paar, Stuttgart, Germany) with a 50 mm parallel plate at a gap of 0.5 mm. The amplitude sweep was conducted from 0.01% to 100% strain with a 1 Hz angular frequency to identify the linear viscoelastic region. The frequency sweep was performed from 1 to 100 s^{−1} (Tung et al., 2017) with a constant 1% strain (within the linear viscoelastic region). The viscosity was measured using the steady shear mode with the shear rate from 0.01 to 100 s^{−1}. Samples were characterized at 37°C.
Microfabrication and semen injection
Request a detailed protocolThe microfluidic device was made of polydimethylsiloxane using a standard soft lithography protocol. The diameter of the circular quiescent zone was 500 µm and the height of the chamber was 25 µm. Diluted semen was injected into the microfluidic device using gravity and the flow generated in the channel was controlled by changing the height of the semen container. Because sperm rheotaxis occurs under a very low shear rate (0.6 s^{−1}), using gravity instead of conventional syringe pumps is more efficient for obtaining and controlling low flow rates.
Rheotaxisbased sperm isolation and phasecontrast microscopy
Request a detailed protocolTo isolate motile bovine sperm inside the quiescent reservoir, we used a microfluidic corral system that isolated motile swimmers based on their ability to move upstream. As we injected the sample at an injection rate of 1.2 µLh^{−1}, sperm with motilities higher than 53.2 µms^{−1} could swim upstream and enter the quiescent zone, which was filled with TALP, allowing us to study sperm movement with minimal fluid mechanical noise. Sperm movement was observed with a Nikon Eclipse TE300 inverted phasecontrast microscope (20× and 40× magnifications) and recorded with an Andor Zyla 5.5 sCMOS camera (25 and 50 frames/s).
Cell tracking and zeroth harmonic measurement
Request a detailed protocolSperm trajectories and other motilityrelated characteristics were analyzed using ImageJ and MATLAB. To identify the harmonics of midpiece bending, we first removed noise and background using Gaussian filter and image subtraction. We then binarized the images taken from the sperm at 25 frames/s and measured the deviation of the midpiece (i.e., the segment located at 10 ± 1 µm from the head) from the centerline using the optical flow Flareback method. Taking the fast Fourier transform of the bending, we identified the zeroth, first, and second harmonics of the bending signal. A simple method for measuring the amplitude of the zeroth harmonic involves measuring the maximum bending toward the left $\left({y}_{L}\right)$ and right $\left({y}_{R}\right)$ sides of the swimmer. Therefore, the magnitude of the zeroth harmonic can be calculated using the following equation:
The main advantage of our method is its simplicity, as tracking the whole flagellum was not required for measurement of zeroth harmonic.
Numerical simulation
Beating pattern
Request a detailed protocolTo model the beating pattern of a sperm, we posited that the flagellation in one beat could be described by a sine wave at a temporal interval of $\left[\pi {\varphi}_{0},{\varphi}_{0}\right]$, such that ${\varphi}_{0}\in \left[\pi ,2\pi \right]$. ${\varphi}_{0}$ determines the level of asymmetry in the beating pattern. Thereafter, by evenly extending the function, we obtained the beating patterns, where ${\varphi}_{0}$ determined the asymmetry in the beating. For example, ${\varphi}_{0}=2\pi $ corresponds to $\sigma =1,$ and thus symmetric beating, whereas ${\varphi}_{0}<2\pi $ results in $\sigma <1$, and thus asymmetric beating. We then applied a fast Fourier transform on the beating patterns to determine their temporal frequencies. These steps were performed using MATLAB (version R2017a).
Finite element method simulations
Request a detailed protocolTo obtain the velocity field imposed by the swimmer model shown in Figure 4—figure supplement 1 and determine the farfield hydrodynamic interactions, we first imported the cylindrical structure of the quiescent zone to the COMSOL MULTIPHYSICS (version 5.2) platform. Two orthogonal Gaussian pulse functions (defined in the $x$ and $y$ directions) were used to define each point force in the swimmer model. The mathematical form of the pulse is a 2D Gaussian distribution, as follows:
We used ${x}_{0}$,${y}_{0}$ to move and ${\sigma}_{x}$, ${\sigma}_{y}$ to focus the point forces arbitrarily. This strategy was chosen to lower the computational cost and avoid issues related to using small volumetric forces and their associated meshing problems in the finite element method. Finally, assuming that sperm swim in a quasi2D plane that is located 5 μm below and parallel to the top surface, we solved the Stokes (Equation 21) and mass conservation equations for varying $\mathrm{\Delta}\theta $ values:
In Equation 21, p is the pressure, μ is the dynamic viscosity of the TALP medium ($3.2\text{}\mathrm{m}\mathrm{P}\mathrm{a}\mathrm{s}$), v is the fluid velocity, r is the position, ${r}_{i}$ is the position of the point force ${f}_{i}$, and $\delta $ is the Dirac delta function. The results obtained from this section are demonstrated in Figure 4—figure supplement 1B. Then, by integrating the velocity field imposed by the sperm, we obtained the drift velocity toward the wall, as shown in Figure 4—figure supplement 1C.
To find the torque imposed on the sperm in nearfield conditions through lubrication approximation, we used a finite element method to solve Stokes and mass conservation equations for the configuration shown in Appendix 1. Given that the contribution of pressure in the stress tensor dominates that of viscous stress (SI, Section XI) ($pI\gg \mu (\nabla v+{\nabla v}^{T})$), we extracted the pressure exerted on the sperm (Figure 4—figure supplement 5) at varying incident angles ($90\xb0<\beta <90\xb0$) for a constant progressive velocity (${V}_{p}=80\mathrm{\mu}\mathrm{m}/\mathrm{s}$). The torque exerted on the sperm from the wall and the corresponding angular velocity was then calculated:
In Equation 22, $p$ is pressure, ${x}_{cm}$ is the coordinate of the sperm center of mass, $L$ is the sperm length, and $\alpha $ is a fitting parameter. Note that $\dot{\beta}$ is linearly correlated to ${V}_{p}$.
Appendix 1
I Rheology of the sperm media
To increase the viscosity and viscoelasticity of the standard TALP medium, we added 1–4% of PVP and 0.25–1% of PAM. The rheological properties of the prepared solutions are shown in Figure 1—figure supplement 1. Note that ${G}^{\text{'}}$ of PVPbased solutions with percentages lower than 4% was very low and not detected by the rheometer.
II Modeling the sperm beating patterns and Fourier analysis
To model beating patterns that resemble those of sperm flagella, we studied the pattern in one cycle of flagellar beating using a traveling sine wave with a temporal phase of $\varphi \left(t\right)=\omega t$ (with $\omega =40$ Hz) in the range of $\pi {\varphi}_{0}\le \varphi \left(t\right)\le {\varphi}_{0}$, such that ${\varphi}_{0}\in [\pi ,2\pi ]$. Notably, ${\varphi}_{0}=2\pi $ corresponds to a fully symmetric beating pattern and thus, $\sigma =1$, whereas lower values of ${\varphi}_{0}$ correspond to lower $\sigma $ values and higher asymmetry in the beating pattern.
Positing that the beating pattern will identically repeat in time, we constructed the even extension of the partial sine wave to form the flagellar beating function over time. Figure 2—figure supplement 1 shows the modeled beating patterns of untethered sperm, which resemble the flagellar beating observed inside the quiescent zone, especially that of the mid and principal piece, which we are interested in. To model the beating pattern, we considered the length preservation constraint as well because the sine wave form does not preserve the length in one beat.
Length preservation constraint
Knowing that the sperm length (${L}_{0}$) is constant over time, we must find the position of its tip on the x axis in time ($L\left(t\right)$), to satisfy the length preservation constraint:
Finding $L\left(t\right)$ enables us to define the beating pattern in $x\in \left[0,L\left(t\right)\right]$ and thereby solve the problem of length preservation.
Assuming that $\forall i,{\left(k{a}_{i}\right)}^{2}\ll 1$:
Plugging Equation S2 into Equation S1 yields Equation S3:
The following calculations simplify Equation S3 and break it down into ${I}_{1}$ and ${I}_{2}$:
${I}_{1}$ has no $x$dependent variable and can be calculated as follows:
To calculate ${I}_{2}$, we first calculate the first term of the equation:
Knowing that
we can use the following approximation:
Using the same method to calculate the second and third terms of ${I}_{2}$ shows that:
Given that $\delta \left(t\right)\ll L$, and $L\left(t\right)$, ${I}_{2}$ can be neglected compared to ${I}_{1}$.
Therefore, we may write $L\left(t\right)$ as follows:
After modeling the flagellar beating pattern and considering length preservation and small amplitude constraints, we analyzed the beating pattern and the resulting sperm motions by applying the Fourier transform to identify temporal frequencies (Figure 2—figure supplement 2). Interestingly, with increasing temporal asymmetry in the beating pattern and thus $\sigma $, the frequency of the main (first) harmonic increases, whereas its amplitude decreases. Moreover, the zeroth and second harmonics simultaneously appear in the frequency domain; however, the amplitudes of the second and higher harmonics are much smaller than that of the zeroth and first harmonics, which enables us to approximate the level of asymmetry as follows:
III Resistive force theory
In this section, we used resistive force theory to derive equations describing the forces produced by each segment of the flagellum and follow the presentation used by Friedrich et al. and Saggiorato et al. The velocity of each segment in the y direction ($V$) can be decomposed into its tangential and normal components ${V}_{T}$ and ${V}_{N}$, respectively, using $\alpha $, which is the tangential angle (Equations S14–S17):
According to the resistive force theory, the forces produced by each element in the tangential and normal direction are linearly related to the velocity in those directions:
where ${\xi}_{T}$ and ${\xi}_{N}$ are drag coefficients in the tangential and normal direction, respectively. Because the amplitudes of all harmonics are small in comparison to sperm length, we may posit the following assumptions:
Using the approximations in Equations S20 and S21, we can write the tangential and normal velocities and forces using Equations S22–S25:
The force produced by each segment in the x and y directions can be described as follows:
Plugging Equations S24 and S25 into Equations S26 and S27:
yields Equations S30 and S31, which describe the forces produced by the flagellum in the x and y directions.
IV Force and torque produced by the flagellum
Force in the x direction
The forces produced by a segment of the flagellum moving with $y(x,t)$ in the tangential and normal directions can be described by Equations S31 and S32, where ${\xi}_{T}$ and ${\xi}_{N}$ are the corresponding drag coefficients, respectively. Plugging Equation S32 into Equation S30 yields Equation S34.
For $i\ne j$, the time average of ${f}_{x}$ would be zero, thus we only retain terms with $i=j$. This yields Equations S35 and S36.
By calculating the time average of Equation S36, the average force produced by each segment of the flagellum in the x direction can be described by Equation S37.
Integrating $\stackrel{~}{{f}_{x}}$ over the flagellum, the total force and velocity produced in the x direction are:
Force in the y direction and the corresponding torque
Plugging Equation S32 into Equation S31 yields Equation S40.
One may rewrite Equation S40 in the form of Equation S41.
Because we are interested in the time average values of ${f}_{y}$, the following terms of Equation S41 are nonzero:
Therefore, Equation S41 is reduced to Equations S42 and S43:
Taking the average of Equation S43, the average force produced by each segment in the y direction is:
The time average of the force produced by each segment of the flagellum in the y direction is not constant and is a function of x, meaning that the effect of the zeroth harmonic can be seen in the force produced in the y direction. Integrating the forces produced by each segment in the y direction over the flagellum, the total force produced in the y direction becomes 0:
However, the magnitude of the force produced in the front half of the sperm is nonzero and equal to the force produced in the rear half of the sperm:
We will use two forces for the $\gamma $ ratio, which is required for the Stokeslet description of the microswimmer model (Section VIII of Appendix 1). Although the total force produced in the y direction is 0, the torque produced by the flagellum is not:
To find the angular velocity of the sperm ($\stackrel{~}{\Omega}$), we need to calculate the torque produced by drag as well:
Considering the zero net torque constraint, we calculated the angular velocity of the sperm:
Because the curvature of the sperm trajectory is $\stackrel{~}{\Omega}{{\stackrel{~}{V}}_{p}}^{1}$ (where $\stackrel{~}{\Omega}$ is the angular velocity of the sperm and ${\stackrel{~}{V}}_{p}$ is the sperm velocity), the curvature of sperm path is:
We may write the $\gamma $ ratio as follows:
V Rolling facilitates sperm rheotaxis
Sperm rheotaxis and the subsequent angular velocity imposed by external shear flow, which aligns the sperm in the upstream direction, can be described by an Adlertype equation:
where $\theta $ is the angle between sperm orientation and flow stream $\left({\mathrm{cos}}^{1}\left(\frac{{V}_{P}\bullet {V}_{F}}{\Vert {V}_{P}\Vert \Vert {V}_{F}\Vert}\right)\right)$, $\gamma $ is the shear rate, and A is a constant. Therefore, sperm orientation with respect to the flow evolves with ${\Omega}_{Rh}$ superimposed with ${\Omega}_{Asy}$, which is caused by the asymmetry in the beating pattern (Figure 2—figure supplement 5A):
Considering that the angle between sperm orientation and flow while swimming upstream (${\theta}_{UP}$) is the angle at which $\dot{\theta}=0$:
Now, we include the influence of rolling by replacing $\Omega $ with $\Pi \left(t\right)\Omega $:
Taking the average value of ${\theta}_{UP}\left(t\right)$ and assuming that $\stackrel{~}{\mathrm{\Pi}\left(t\right)}\ll 1$:
Equation S57 suggests that rolling results in a significant decrease in ${\stackrel{~}{\theta}}_{UP}$. To see how this decrease in ${\stackrel{~}{\theta}}_{UP}$ facilitates sperm rheotaxis, we calculated the net upstream component of sperm motion with and without the rolling component.
According to the velocity components shown in Figure 2—figure supplement 5B, the net upstream velocity is:
where ${\theta}_{UP}$ includes $\Pi \left(t\right):$
Taking the average value of ${V}_{\mathrm{U}\mathrm{P}}\left(t\right)$ and assuming that $\stackrel{~}{\Pi \left(t\right)}\approx 0$, upstream velocity can be written as follows:
At high frequencies of rolling, $\stackrel{~}{\theta}}_{UP$ approaches 0 (as discussed previously) and therefore,
Thus, at the frequent rolling limit, the tangential velocity produced by flagellation is inclined toward upstream of the flow and thus rheotactic behavior is efficient, whereas without the rolling component ${V}_{\mathrm{U}\mathrm{P}}$ is:
which is smaller than ${V}_{P}{V}_{F}$.
We experimentally measured ${\stackrel{~}{\theta}}_{UP}$ for rolling and nonrolling sperm. For rolling sperm, $\stackrel{~}{\theta}}_{UP}=10\pm {5}^{\circ$, whereas for nonrolling sperm, $\stackrel{~}{\theta}}_{UP}=40\pm {10}^{\circ$, which is consistent with our theoretical analysis. By tracking the sperm under flow (Figure 2—figure supplement 6), we noticed that the swimming velocity was much higher $\left({\stackrel{~}{V}}_{UP}=60\pm 10\mathrm{\mu}\mathrm{m}/\mathrm{s}\right)$ when the sperm exhibited rolling compared to nonrolling $\left({\stackrel{~}{V}}_{UP}=30\pm 10\mathrm{\mu}\mathrm{m}/\mathrm{s}\right)$.
VI Noise in the amplitude and phase of the first and higher harmonics
To identify the source of diffusive motion at the center of the circular path, we employed results obtained from the resistive force theory (Equations S38 and S50), and included white Gaussian noises in the amplitudes and phases of all harmonics:
where $\u27e8{\eta}_{n}\left(t\right)\u27e9=\u27e8\xi \left(t\right)\u27e9=0$, $\u27e8{\eta}_{n}\left(t\right){\eta}_{n\text{'}}\left(t\text{'}\right)\u27e9={D}_{n}{\delta}_{n{n}^{\text{'}}}\delta \left(t{t}^{\text{'}}\right)$, and $\u27e8\xi \left(t\right)\xi \left(t\text{'}\right)\u27e9={D}_{\omega}\delta \left(t{t}^{\text{'}}\right)$. Inserting the noise terms, we solved the equations of motion for a sperm with arbitrary and constant progressive and angular velocities and noticed that except for the noise in the amplitude of the zeroth harmonics, all other sources of noise did not change the pathway and hence do not result in a diffusive motion of the center (Figure 3—figure supplement 1). By contrast, the noise in the amplitude of the zeroth harmonic produces a diffusive component during circular motion. Moreover, the SNR of the zeroth harmonic determines the diffusion coefficient of the center.
VII Infrequent rolling and intermittent search
Our experimental measurement of the thickness of the layer shown in Figure 3A, that is, $\delta $, indicates that it correlates with $\sqrt{t}$, as in a normal diffusion process (Figure 3—figure supplement 2). Furthermore, we observed that the sperm with circular diffusive motion at certain ambient viscosities might show infrequent rolling that relocates the center. Measuring the velocity of the center, we demonstrated that relocation happens much faster than diffusion, and subsequently, one can posit the relocations as abrupt interruptions to the diffusion process (Figure 3—figure supplement 3). Because sperm circular motion includes a diffusion process with $\delta \propto \sqrt{t}$, and infrequent rolling results in abrupt relocations of the center, we now propose a model to investigate how infrequent rolling increases the area swept by circular motion.
The area swept by the sperm through a sole diffusion process is equivalent to the area of a ring with the thickness δ and curvature $\kappa $, which is approximately $2\pi \delta {\kappa}^{1}$ for $\delta \ll {\kappa}^{1}$. Assuming that $\kappa $ is constant in time, the area of the ring increases in time with the square root of time ($\u27e8A\left(t\right)\u27e9\propto \sqrt{t}$) (Figure 3—figure supplement 4). Now, suppose that the time between (i−1)th and ith rolling is ${T}_{i}$, during which the sperm swims in diffusive circular trajectories. The area that the sperm sweeps after $N$ rolling occurrences is proportional to $\sum _{i=1}^{N}\sqrt{{T}_{i}}$, which is larger than the area swept without rolling, $\sqrt{\sum _{i=1}^{N}{T}_{i}}$. This twophase motion in the sperm is reminiscent of intermittent search algorithms that have been demonstrated to be more efficient than a simple diffusion process.
VIII The flow field produced by the propulsive and circular components of motility
The flow field produced by the swimmer model in Figure 4—figure supplement 1 is the superposition of the flow field produced by the propulsive $\left(\overrightarrow{f}\right)$ and circular $\left(\overrightarrow{{f}^{\u2033}}\text{}\mathrm{a}\mathrm{n}\mathrm{d}\text{}\overrightarrow{{f}^{\u2033}}\right)$ components of motility and their corresponding drags $\left(\overrightarrow{f},\overrightarrow{{f}^{\prime}},\text{}\mathrm{a}\mathrm{n}\mathrm{d}\text{}\overrightarrow{{f}^{\prime}}\right)$. To identify the mechanism of reduction in the farfield drift velocity caused by circular components, we separately simulated the propulsive term with its corresponding drag (Figure 4—figure supplement 2A) and the circular term with its corresponding drag (Figure 4—figure supplement 2B). The arrows show the normalized vector field, and the magnitude of the flow is represented in color. The flow field produced by the circular terms suppresses the flow field produced by the propulsive term, and thus, reduces the farfield drift velocity.
IX The analytic solution for the swimmer model
The swimmer model for calculating the farfield drift velocity is shown in Figure 4—figure supplement 3. The velocity field imposed by a source point is known as the Stokeslet (Equation S64), which is the most fundamental solution for the Stokes equation. Based on the linearity in the Stokes equation, the contribution of the actively swimming sperm on the fluid flow can be described by superimposing the flow fields produced by each point force. For the sake of simplicity, we write the imposed velocity field in three terms, including:
The velocity field imposed by $\overrightarrow{f}$ and $\overrightarrow{f}$ (Equation S65),
$\overrightarrow{f\text{'}}$ and $\overrightarrow{f\text{'}}$ (Equation S66), and
$\overrightarrow{f\text{'}\text{'}}$ and $\overrightarrow{f\text{'}\text{'}}$ (Equation S67), in which the magnitudes of $\overrightarrow{f\text{'}}$ and $\overrightarrow{f\text{'}\text{'}}$ are equal to $\gamma \overrightarrow{f}$.
Note that $\overrightarrow{d}$ is a vector with a magnitude equal to the sperm length in $\widehat{d}$ direction and $\overrightarrow{d\text{'}}=\left[\begin{array}{cc}cos\u2206\theta & sin\u2206\theta \\ sin\u2206\theta & cos\u2206\theta \end{array}\right]\overrightarrow{d}$.
Using Equations S79, S80, and S81, the velocity field imposed by the swimmer model is:
with:
$\stackrel{=}{B}=\left[\begin{array}{cc}\mathrm{cos}\mathrm{\Delta}\theta & \mathrm{sin}\mathrm{\Delta}\theta \\ \mathrm{sin}\mathrm{\Delta}\theta & \mathrm{cos}\mathrm{\Delta}\theta \end{array}\right]$ and $\stackrel{=}{I}=\left[\begin{array}{cc}1& 0\\ 0& 1\end{array}\right].$
To calculate the farfield drift velocity caused by a wall with noslip boundary condition, we replaced the wall using the mirror image of the swimmer in the boundary. Accordingly, the velocity field imposed on the sperm that causes the farfield drift toward the wall (${U}_{w}$) is ${u}_{T}$ with $r=2h$, where $h$ is the distance between the sperm and the wall. For $\phi =\pi /2$, $A=1+\frac{3}{2}\mathrm{s}\mathrm{i}\mathrm{n}\left(2\mathrm{\Delta}\theta \right)$ and $\gamma \ll 1$:
This equation can also be written as:
where ${U}_{w}^{p}$ is the drift velocity when the circular terms are 0.
Interestingly, at ${\stackrel{~}{\u2206\theta}}_{n}=\frac{1}{2}{\mathrm{sin}}^{1}\frac{2}{3}\approx 21\xb0$, the swimmer experiences no attraction toward the walls and becomes neutral. Given the relation between $\stackrel{~}{\mathrm{\Delta}\theta}$, $\stackrel{~}{\kappa}$, and ${\stackrel{~}{V}}_{p}$, the corresponding neutral curvature is:
X The effect of the wall on the center of the circular path
Now, using Equation S70, we develop a model to calculate the average drift velocity imposed by the wall on the circular path’s center in one round of circulation (Figure 4—figure supplement 4):
with
where N is the number of beats required to swim a complete circle.
Equations S73 and S74 are derived by plugging ${\phi}_{k}$ and $\u2206\theta $ into Equation S72 and ${h}_{k}$:
Using Equation S73 and Equation S74, the average farfield drift velocity imposed on the sperm by the wall can be written as follows:
The calculations based on Equation S75 are shown in Figure 4—figure supplement 1D, which indicates that the average drift velocity imposed by the wall on the circular path’s center is much lower than the one produced by noise in the amplitude of the zeroth harmonic. Therefore, the circular path is solely evolving with diffusion of the center caused by the noise in the amplitude of the zeroth harmonic, rather than by the farfield influence of the wall.
XI Nearfield interactions and lubrication approximation
To develop a hydrodynamic model and describe sperm nearfield interactions with the wall at the lubrication limit, we solved the Stokes equation for a single sperm swimming progressively near a wall with the noslip boundary condition (Figure 4—figure supplement 5). The stress tensor is:
At distances adequately close to the wall, the contribution of the pressure dominates the viscous term, such that $\sigma \approx p\stackrel{\u033f}{I}$. Accordingly, the torque exerted on the sperm by the wall can be written as follows:
Neglecting sperm mass, the net torque applied on the sperm is 0, meaning that the drag torque cancels out the torque exerted by the wall. This constraint gives us the following equation:
where $\beta $ is the angle of the sperm swimming direction with respect to the wall. To find $\dot{\beta}$, we performed finite element simulations to find the pressure distribution between the sperm and wall, and the results are shown in Figure 4—figure supplement 5. Note that the sperm tangential velocity in our simulations was arbitrary and constant, the pressure was found to be linearly correlated to tangential velocity.
XII The effect of rolling on sperm–wall interaction
In our finite element simulations, we did not include the role of rolling. However, a simple way to include rolling is to include $\mathrm{\Pi}\left(t\right)$ in the rotation caused by asymmetric flagellar beating. Therefore, we can write $\dot{\beta}$ as follows:
To find ${\beta}_{s}$ in the presence of rolling, we set $\dot{\beta}=0$, which yields:
Taking the average value of ${\beta}_{S}\left(t\right)$, and assuming that $\stackrel{~}{\mathrm{\Pi}\left(t\right)}\ll 1$, $\stackrel{~}{{\beta}_{S}^{\ast}}\ll {\beta}_{S}$, suggests that a rolling component in motility significantly decreases the average angle of the sperm with respect to the wall. To find $\beta \left(t\right)$ at $\beta \to \stackrel{~}{{\beta}_{S}^{\ast}}$ limit, we plug $\stackrel{~}{{\beta}_{S}}$ into Equation S79:
where $g\left(\stackrel{~}{{\beta}_{S}^{\ast}}\right)\approx 0$. Relying upon Equation S81, $\beta \left(t\right)$ is a triangular function with an offset of $\stackrel{~}{{\beta}_{S}^{\ast}}$:
Our experimental measurement of $\beta \left(t\right)$ for the rolling sperm swimming along the wall agrees with Equation S82. Our results demonstrate that the angle between the rolling sperm and wall is a triangular function with an offset of $\stackrel{~}{{\beta}_{S}^{\ast}}\sim {10}^{\circ}$ (Figure 4—figure supplement 6).
Data availability
All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. All data related to this paper are deposited in https://doi.org/10.5061/dryad.ngf1vhhtd.

Dryad Digital RepositoryData from: Rolling controls sperm navigation in response to the dynamic rheological properties of the environment.https://doi.org/10.5061/dryad.ngf1vhhtd
References

Stratlets: low Reynolds number pointforce solutions in a stratified fluidPhysical Review Letters 105:084502.https://doi.org/10.1103/PhysRevLett.105.084502

Intermittent search strategiesReviews of Modern Physics 83:81–129.https://doi.org/10.1103/RevModPhys.83.81

Rethinking the relationship between hyperactivation and chemotaxis in mammalian spermBiology of Reproduction 83:507–513.https://doi.org/10.1095/biolreprod.109.083113

Observations on bull sperm rotationBiology of Reproduction 10:78–84.https://doi.org/10.1095/biolreprod10.1.78

Sperm guidance in mammals  an unpaved road to the eggNature Reviews. Molecular Cell Biology 7:276–285.https://doi.org/10.1038/nrm1893

Hydrodynamics of sperm cells near surfacesBiophysical Journal 99:1018–1026.https://doi.org/10.1016/j.bpj.2010.05.015

Physics of microswimmers—single particle motion and collective behavior: a reviewReports on Progress in Physics 78:56601.https://doi.org/10.1088/00344885/78/5/056601

Microswimmers near surfacesThe European Physical Journal Special Topics 225:2333–2352.https://doi.org/10.1140/epjst/e2016600706

Highprecision tracking of sperm swimming fine structure provides strong test of resistive force theoryThe Journal of Experimental Biology 213:1226–1234.https://doi.org/10.1242/jeb.039800

Mammalian sperm motility: observation and theoryAnnual Review of Fluid Mechanics 43:501–528.https://doi.org/10.1146/annurevfluid121108145442

Noise and synchronization in pairs of beating eukaryotic flagellaPhysical Review Letters 103:168103.https://doi.org/10.1103/PhysRevLett.103.168103

Geometrical guidance and trapping transition of human sperm cellsPhysical Eeview. E, Statistical, Nonlinear, and Soft Matter Physics 89:32720.https://doi.org/10.1103/PhysRevE.89.032720

Mechanical tuning of mammalian sperm behaviour by hyperactivation, rheology and substrate adhesion: a numerical explorationJournal of the Royal Society Interface 13:20160633.https://doi.org/10.1098/rsif.2016.0633

Hydrodynamic clustering of human sperm in viscoelastic fluidsScientific Reports 8:15600.https://doi.org/10.1038/s41598018335848

Mechanisms of sperm chemotaxisAnnual Review of Physiology 70:93–117.https://doi.org/10.1146/annurev.physiol.70.113006.100654

Hydrodynamic interaction of microswimmers near a wallPhysical Review. E, Statistical, Nonlinear, and Soft Matter Physics 90:13010.https://doi.org/10.1103/PhysRevE.90.013010

Undulatory swimming in nonNewtonian fluidsJournal of Fluid Mechanics 784:R4.https://doi.org/10.1017/jfm.2015.595

Active phase and amplitude fluctuations of flagellar beatingPhysical Review Letters 113:48101.https://doi.org/10.1103/PhysRevLett.113.048101

Bovine sperm hyperactivation is promoted by alkalinestimulated Ca2+ influxBiology of Reproduction 76:660–665.https://doi.org/10.1095/biolreprod.106.055038

Rheotaxis guides mammalian spermCurrent Biology : CB 23:443–452.https://doi.org/10.1016/j.cub.2013.02.007

Predominance of sperm motion in cornersScientific Reports 6:26669.https://doi.org/10.1038/srep26669

Sperm motility in modulated microchannelsNew Journal of Physics 21:013016.https://doi.org/10.1088/13672630/aaf544

Human sperm steer with second harmonics of the flagellar beatNature Communications 8:1415.https://doi.org/10.1038/s4146701701462y

Control of hyperactivation in spermHuman Reproduction Update 14:647–657.https://doi.org/10.1093/humupd/dmn029

Mammalian sperm interactions with the female reproductive tractCell and Tissue Research 363:185–194.https://doi.org/10.1007/s0044101522442

Sperm transport in the female reproductive tractHuman Reproduction Update 12:23–37.https://doi.org/10.1093/humupd/dmi047

Viscoelastic fluid response can increase the speed and efficiency of a free swimmerPhysical Review Letters 104:038101.https://doi.org/10.1103/PhysRevLett.104.038101

Emergence of upstream swimming via a hydrodynamic transitionPhysical Review Letters 114:108102.https://doi.org/10.1103/PhysRevLett.114.108102

The effectiveness of resistive force theory in granular locomotionPhysics of Fluids 26:101308.https://doi.org/10.1063/1.4898629
Decision letter

Raymond E GoldsteinReviewing Editor; University of Cambridge, United Kingdom

Anna AkhmanovaSenior Editor; Utrecht University, Netherlands
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
Acceptance summary:
Mammalian sperm cells achieve locomotion in a liquid environment by beating their flagellum in wavelike pattern. It is wellknown that the resulting forward motion is often accompanied by a rolling motion around the cell's longitudinal axis. However, the function of this rolling motion in the context of navigation is not well understood. The authors combine experiments and simulations to show how rolling aids sperm navigation along surfaces and external fluid flows.
Decision letter after peer review:
Thank you for submitting your article "Rolling controls sperm navigation in response to the dynamic rheological properties of the environment" for consideration by eLife. Your article has been reviewed by 2 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Anna Akhmanova as the Senior Editor. The reviewers have opted to remain anonymous.
The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.
Essential revisions:
1. The conclusion that rolling is "in response to the dynamic rheological properties of the environment", while compelling, could be better supported. To show that sperm swimming responds to the rheology of the medium, the authors would first need to track the same sperm in media of different rheologies. Second, it is unclear if the sperm is responding to the rheology, or the swimming merely manifests differently in response to the rheology. Third, the data presented for the rheological characterization (Figure S1) do not support the clean delineation between the PVP media being "viscous" and the PAM media being "viscoelastic": the PVP media are 12 orders of magnitude less viscous than the PAM, and in both cases, the viscosity appears to sharply rise at low shear rates (although this is likely noise as it also manifests for the TALP). The PVP media appear to show weak shear thinning, similar to the PAM. Furthermore, only the shear viscosity is presented for both solutions; to adequately characterize viscoelasticity, the authors need to show data for e.g. the first normal stress difference vs shear rate, or the elastic and viscous moduli vs oscillation frequency, for example.
2. The writing in general is dense and highly technical; it would be useful to the reader if the authors could provide broader introduction and Discussion sections that clearly discuss the context of the work, the findings, and the implications of the findings.
3. The authors say in lines 9899 that an increase in viscoelasticity slightly increased the propulsive velocity. I do not see this in the data.
4. In Figure 1F, what do the tildes indicate? And what are the units of the quantities plotted?
5. Lines 162164: The authors discuss the diffusion coefficient, but I do not see a quantitative description in the actual data. Can they plot e.g. a MSD vs time?
6. Figure 2AB: What are the units? What is shown in blue? What is the influence of experimental noise (e.g. in the imaging) on these results?
7. Figure 2G: It is unclear what the points vs the lines are.
8. Figure 2H: The authors claim this is a normal distribution, but it does not appear to be so.
9. Experimental uncertainties are improperly reported throughout. E.g. 38.7 +/ 12.2 should be reported as 40 +/ 10, given the uncertainty in the measurement.
10. What is the influence of thermal diffusion in the results throughout?
11. Line 46: The phrase "To avoid errors arising from studying sperm motility under external fluid flow, …" is unclear – which specific potential "errors" do the authors mean?
12. Line 119: It is not obvious how "Thus, rolling causes progressive motion." follows from the preceding sentences.
13. Line 134: period before equation 1 seems unnecessary?
14. Line 135: How exactly is x defined?
15. Line 142: "(24, 29)(24, 29)(24, 29)(24, 29)"?
16. Line 343 344: the notation is unclear. Does the tilde indicate a time or ensemble average ? Assuming the former, why is the time average of \Pi still timedependent ? Shouldn't it read "average of \β_s(t)" in line 344?
17. The authors should carefully check and correct the reference list. The listing of doi 's seems inconsistent and several references are incomplete. For example:
a. Line 524: doi necessary?
b. Line 525: reference incomplete.
c. Line 559: reference incomplete.
d. Line 587: reference incomplete.
e. Line 601: page number seems incorrect.
f. Line 605: page number missing?
g. Line 611: page number seems incorrect.
h. Line 611: page number seems incorrect?
18. Bold font on rhs. of equation S44 should be removed, similar in S51, S52.
19. Please check constant use of multiplicative dots in SI (sometimes lower dots and sometimes centered dots are used).
20. Major point: The definition of eta in equation S65S68 is either incomplete or incorrect.
21. What exactly do the authors mean by Gaussian noise in S6667?
22. S67 does not look standard Gaussian white noise, which should be δ correlated. Are the eta(t)'s time correlated or simply random numbers of variance D? In the latter case, the stochastic process eta(t) would be discontinuous?
23. Below equation S77 the authors mention "white Gaussian noise" which seems inconsistent with S67.
24. Also, if eta(t) has a Gaussian distribution, then we do not see how S73 can be satisfied, since the tails of a Gaussian distribution extend to infinity – in particular, when using D=1, as stated below equation S77.
25. Which type of stochastic calculus do the authors use to arrive at S74? It looks like Stratonovich but this is not explicitly stated.
26. Based on the above, we are not convinced that equations S65S77 are correct.
This part of the SI should be carefully revised and corrected as needed, or omitted.Reviewer #2 (Recommendations for the authors):
The authors provide a functional explanation for the observation that mammalian sperm roll around their longitudinal axis: it enables a sperm to continue to move progressively. The authors further sought to provide evidence that rolling responds to the rheological properties of the environment.
The use of imaging to track single sperm and directly observe rolling events is a major strength.
The rheological characterization of the media requires further analysis.
The authors nicely show that rolling enables a sperm to continue to move progressively. However, the conclusion that rolling is "in response to the dynamic rheological properties of the environment", while compelling, is not well supported. The authors tested this hypothesis by studying sperm swimming in a "viscous medium" of TALP + PVP and a "viscoelastic medium" of TALP + PAM, although the rheological characterization of these media does not support this clean delineation.
The observation that sperm rolling has a key functional role broadens understanding of sperm swimming. If the authors could more rigorously support the claim that rolling is modulated in response to the dynamic rheological properties of the environment, then this would broaden understanding of sperm swimming in complex environments during e.g. fertilization.Reviewer #3 (Recommendations for the authors):
The wavelike beating patterns of sperm flagella can vary between cells, with some cells exhibiting more symmetric beat patterns than others. To investigate the dependence of rolling motion on beat asymmetry and rheological properties, the author studied bovine sperm motion in a microfluidic device under simple shear flow and within a quiescent reservoir.
The authors' results suggest that rolling becomes suppressed at relatively lower viscosities/viscoelasticity for sperm exhibiting less symmetric beat patterns, with suppression of rolling more sensitive to viscoelasticity than to viscosity. They found that rolling sperm swim on mostly straight trajectories, whereas nonrolling sperm display 2D planar beat forms resulting in circular trajectories.
Building on previous theoretical work, the authors formulate a minimal 2Dprojected model that accounts for rolling through a timedependent prefactor \Pi(t) in equation 4. They find that this model recapitulates their main experimental observations. In particular, their joint theoretical and experimental analysis suggests that, although the rolling component is not required for sperm rheotaxis, it facilitates rheotactic behavior by minimizing the angle between sperm orientation and externally applied fluid flow, thus maximizing the upstream component of sperm motion.
Overall, I find this comprehensive characterization of sperm rolling interesting and it appears to close an existing gap in the literature.
https://doi.org/10.7554/eLife.68693.sa1Author response
Essential revisions:
1. The conclusion that rolling is "in response to the dynamic rheological properties of the environment", while compelling, could be better supported. To show that sperm swimming responds to the rheology of the medium, the authors would first need to track the same sperm in media of different rheologies.
We thank the reviewers for this helpful comment. To address it, we have conducted new experiments and demonstrated reversibility in the rolling suppression, as we observed reactivation of rolling upon leaving the viscoelastic (1% PAM) or viscous (4% PVP) solutions and entering less viscous and viscoelastic (TALP) solution. This suppression and reactivation of sperm rolling while migrating between TALP and 4% PVP or 1% PAM solutions resulted in transition of sperm swimming behavior between progressive and circular motion. To fully illustrate this point, we have added a new subfigure (Figure 1K) to Figure 1 of this new version of the manuscript, then explained these new findings in the results (page 7, lines 138149) and further discussed the implications of this observation in the Discussion section (pages 2223 lines 458478). Due to the high significance of this reversibility in rolling suppression, we have revised the overall manuscript to reflect this concept throughout.
Second, it is unclear if the sperm is responding to the rheology, or the swimming merely manifests differently in response to the rheology.
Thank you for this helpful comment. We do not have any evidence for the existence of an active mechanism through which sperm “responds” to the rheology of the environment. We replaced “respond” to “in response” throughout the manuscript to clarify. Furthermore, we clearly stated on page 4, lines 7274 that “sperm swimming behavior manifests differently in response to the rheological properties of the environment”. Unlike small molecules such as caffeine and 4AP as hyperactivation agonists, there is no evidence of a biochemical response to the PVP or PAM, which are large molecules with molecular weight measured in MDa range, so it is likely that the change in sperm swimming behavior is in response to the rheological properties of the solution.
Third, the data presented for the rheological characterization (Figure S1) do not support the clean delineation between the PVP media being "viscous" and the PAM media being "viscoelastic": the PVP media are 12 orders of magnitude less viscous than the PAM, and in both cases, the viscosity appears to sharply rise at low shear rates (although this is likely noise as it also manifests for the TALP). The PVP media appear to show weak shear thinning, similar to the PAM. Furthermore, only the shear viscosity is presented for both solutions; to adequately characterize viscoelasticity, the authors need to show data for e.g. the first normal stress difference vs shear rate, or the elastic and viscous moduli vs oscillation frequency, for example.
Thank you for this very useful comment. We have now measured the viscosity, storage modulus and loss modulus for our polymerbased solutions using a new device which was equipped with large parallel plates (50 mm in diameter). Using such large parallel plates helped us to measure viscosities, as well as G’ and G’’ more precisely. The amplitude sweep was conducted from 0.01% to 100% strain with a 1 Hz angular frequency to identify the linear viscoelastic region. The frequency sweep was performed from 1 to 100 s1 with a constant 1% strain (within the linear viscoelastic region). The viscosity was measured using the steady shear mode with the shear rate from 0.01 to 100 s1. Samples were characterized at 37°C. The results of rheological measurements are now available in Supplementary Information, section I (Figure 1—figure supplement 1). These measurements helped us to make more accurate statements in the main text as well.
Based on these new rheological measurements, we noticed that the viscosity of 1% PAM is two orders of magnitude greater than that of 4% PVP. Therefore, both solutions are more viscous than TALP, however, our new rheological measurements also indicate that the storage modulus (G’) of 1% PAM is also two orders of magnitude greater than that of 4% PVP. Therefore, we refer to PVP as the viscous solution and PAM as the viscoelastic solution.
These new measurements helped us to explain the results shown in Figure 1F and G. As demonstrated in Figure 1 G, the propulsive velocity of nonrolling sperm in viscous solution is very low, indicating that increase of viscosity decreases sperm velocity. Further, in PVP, we also observed a similar decay in the frequency of rolling as we increased the concentration (Figure 1F). However, in PAM, we did not observe a decrease in frequency of rolling as we increased the concentration of PAM up to 1%. We also observed two things in 1% PAM: first, that sperm propulsive velocity was measured to be higher than that of the 4% PVP solution, and second, that more than 90% of sperm cells did not show rolling while in 4% PVP only 40% of sperm cells did not roll. Given that the viscosity of 1% PAM is two orders of magnitude higher than that of 4% PVP, these results, based solely on viscosity, seemed contradictory. Therefore, we needed additional parameters beyond viscosity to explain our observations.
To solve this apparent anomaly, we used the additional rheological measurements asked by the reviewers and found that the storage modulus of 1% PAM was two orders of magnitude greater than that of 4% PVP. Therefore, the elasticity of the swimming media plays a major role in the suppression of rolling without loss in the sperm propulsive velocity. This finding was key to our discussion regarding the sensitivity of rolling suppression to viscoelasticity in the previous version, but in this new version this claim is now supported by the rheological measurements.
Including section I of the Supplementary Information, we have also changed Figure 1G according to your comment and our new experiments. We also reflected these findings on pages 56 (lines 105119), pages 23 (lines 469478), and page 24 (lines 503509).
2. The writing in general is dense and highly technical; it would be useful to the reader if the authors could provide broader introduction and Discussion sections that clearly discuss the context of the work, the findings, and the implications of the findings.
We thank the reviewers for this comment. According to your comment, we revised the introduction and Discussion sections of the manuscript. Please see pages 24 for the revised introduction. Please also see pages 2123 for the revised discussion.
3. The authors say in lines 9899 that an increase in viscoelasticity slightly increased the propulsive velocity. I do not see this in the data.
Thank you for calling our attention to this issue. After careful examination of our results, we noticed that our previous claim was not crucial to our discussion, and given the error in our experimental measurements, the slight increase in sperm propulsive velocity was not significant. However, we added Figure 1G to the updated version of the manuscript. This data is important as it demonstrates that suppression of rolling in the viscous solution comes with a significant loss in the propulsive velocity, while suppression of rolling in the viscoelastic solution is not only more efficient but also does not accompany a loss in the sperm propulsive velocity.
4. In Figure 1F, what do the tildes indicate? And what are the units of the quantities plotted?
Thank you for calling our attention to this issue. The tilde signs are now omitted, and the units are added to this figure.
5. Lines 162164: The authors discuss the diffusion coefficient, but I do not see a quantitative description in the actual data. Can they plot e.g. a MSD vs time?
Thank you for raising this issue. In response, we have plotted the MSD (mean square displacement) of two sperm swimming in diffusive circles and the results are shown in Figure 2E. We noticed that for both sperm, $MSD\sim t$, indicating that the circular motion is indeed diffusive in character.
6. Figure 2AB: What are the units? What is shown in blue? What is the influence of experimental noise (e.g. in the imaging) on these results?
Thank you for calling our attention to this issue. The units are now added to Figure 2AB, the blue line in the previous figure was a fitting curve, but it is omitted in this version and the signal itself is shown in blue in Figure 2A.
We discussed the effect of experimental noise in the main text page 8 (lines 161165). The constant offset in the normalized power spectrum is white noise in our measurement system. Furthermore, the experimental noise coming from our measurement system (image processing) is also included in the peak width around the first harmonic frequency. That is, the peak width centered at ω includes the intrinsic noise originated from flagellar sources, as well as the noise associated with our measurement system.
7. Figure 2G: It is unclear what the points vs the lines are.
Thank you for this comment. To increase clarity, we omitted the points and changed the style of Figure 2G. In HLI graph, the lines are showing the times at which rollings occurred and $\Pi (t)$ is the defined function.
8. Figure 2H: The authors claim this is a normal distribution, but it does not appear to be so.
Thank you for this comment. We agree with the reviewers, as according to the Kolmogorov–Smirnov test (limiting form, Stephens modification, Marsaglia method, and Lilliefors modification), the distribution shown in Figure 2H is not normal. We omitted this claim in the updated version. Note that simulations performed to obtain sperm trajectories (Figure 2I) was based on the actual distribution of ${T}_{\mathrm{\text{SR}}}$, so we did not change Figure 2I in the updated version.
9. Experimental uncertainties are improperly reported throughout. E.g. 38.7 +/ 12.2 should be reported as 40 +/ 10, given the uncertainty in the measurement.
Thank you for calling our attention to this issue. We carefully checked the experimental errors throughout and modified the reported values.
10. What is the influence of thermal diffusion in the results throughout?
Thank you for this question. In response we have added a new paragraph to the manuscript on page 13 (lines 271277). Using rough estimations, we argue that thermal fluctuations are negligible in our study. While our estimation is particular to our study, its main idea was borrowed from two previously published studies (ref. 36 and 37).
11. Line 46: The phrase "To avoid errors arising from studying sperm motility under external fluid flow, …" is unclear – which specific potential "errors" do the authors mean?
Thank you. To specify, we revised the sentence according to your comment. The revised sentence on page 3 (lines 5052) is:
“To avoid errors arising from studying sperm motility under external fluid flow, such as experimental inaccuracies caused by decoupling the effect of flow on sperm motion from active swimming, we employed a multistep approach.”
12. Line 119: It is not obvious how "Thus, rolling causes progressive motion." follows from the preceding sentences.
Thank you for your comment. We omitted this sentence. Instead, we discussed the reversibility of rolling suppression in the next paragraph (page 7, lines 138149) and introduced it as another evidence to support our claim that the rolling component is a key contributor to progressive motility.
13. Line 134: period before equation 1 seems unnecessary?
Thank you! We have replaced the period with a colon.
14. Line 135: How exactly is x defined?
Thank you for your comment. We modified the text on page 8 (line 171) to clarify that the $x$ axis was set parallel to the flagellum at its straightline form.
15. Line 142: "(24, 29)(24, 29)(24, 29)(24, 29)"?
Thank you for calling our attention to this issue. Unfortunately, our reference manager did not work properly. We resolved the issue and deleted this part from the text.
16. Line 343 344: the notation is unclear. Does the tilde indicate a time or ensemble average ? Assuming the former, why is the time average of \Pi still timedependent ? Shouldn't it read "average of \β_s(t)" in line 344?
Thank you for this comment. We have modified the text (pages 1819, lines 377387) and made the mathematical exposition as clear and easy to follow as possible.
The tilde sign indicates time average, so the tilde sign should be over $\Pi (t)$, and not just $\Pi $.
From the text:
“We then can rewrite $\stackrel{\u0307}{\beta}$ and include $\Pi \left(t\right)$ such that:
Note that $g\left(\beta \right)$ is the curve in the phase space that corresponds to $\Omega =0$. Insofar as stability occurs at $\stackrel{\u0307}{\beta}=0$,
in which ${{\beta}_{S}}^{*}$ and ${\beta}_{S}$ are the stable points with and without taking rolling into account. Because $\stackrel{\u0303}{\Pi (t)}\ll 1$, the average of ${\beta}_{S}^{*}(t)$ is:
Equation 15 suggests that, at the frequent rolling limit where $\stackrel{\u0303}{\Pi (t)}$ approaches zero, $\stackrel{\u0303}{{\beta}_{S}^{*}}$ approaches zero as well. Consequently, frequent rollings mitigate the destructive role of $\Omega $ in sperm motion along the wall, thereby yielding faster and longerlasting swimming along the wall (SI, section XII). Because at the frequent rolling limit $\stackrel{\u0303}{{\beta}_{S}^{*}}$ is close to zero, we posit that $g\left(\stackrel{\u0303}{{\beta}_{S}^{*}}\right)\approx 0$, and thus $\stackrel{\u0307}{\beta}$ near the stable point can be written as
$\stackrel{\u0307}{\beta}=g\left(\stackrel{\u0303}{{\beta}_{S}^{*}}\right)+\Pi \left(t\right)\Omega \approx \Pi \left(t\right)\Omega .\phantom{\rule{0.222em}{0ex}}\#\left(16\right)\phantom{\rule{0.222em}{0ex}}$“
17. The authors should carefully check and correct the reference list. The listing of doi 's seems inconsistent and several references are incomplete. For example:
a. Line 524: doi necessary?
b. Line 525: reference incomplete.
c. Line 559: reference incomplete.
d. Line 587: reference incomplete.
e. Line 601: page number seems incorrect.
f. Line 605: page number missing?
g. Line 611: page number seems incorrect.
h. Line 611: page number seems incorrect?
Thank you for calling our attention to this issue. Due to the flexibilities offered by eLife, we are using the PNAS citation style in this version. We resolved the issue and carefully checked and corrected the reference list. Furthermore, we noticed that the following paper was retracted during the revision process, so we omitted this reference from our reference list.
“Human sperm uses asymmetric and anisotropic flagellar controls to regulate swimming symmetry and cell steering”, (DOI: 10.1126/sciadv.aba5168)
18. Bold font on rhs. of equation S44 should be removed, similar in S51, S52
According to your comment, we removed the bold fonts in equations and mathematical expression throughout the main text and Supplementary Information.
19. Please check constant use of multiplicative dots in SI (sometimes lower dots and sometimes centered dots are used).
According to your comment, we carefully revised the Supplementary Information and removed all multiplicative dots, except for the dot product of two vectors. Accordingly, all dots are simple periods.
20. Major point: The definition of eta in equation S65S68 is either incomplete or incorrect.
21. What exactly do the authors mean by Gaussian noise in S6667?
22. S67 does not look standard Gaussian white noise, which should be δ correlated. Are the eta(t)'s time correlated or simply random numbers of variance D? In the latter case, the stochastic process eta(t) would be discontinuous?
23. Below equation S77 the authors mention "white Gaussian noise" which seems inconsistent with S67.
24. Also, if eta(t) has a Gaussian distribution, then we do not see how S73 can be satisfied, since the tails of a Gaussian distribution extend to infinity – in particular, when using D=1, as stated below equation S77.
25. Which type of stochastic calculus do the authors use to arrive at S74? It looks like Stratonovich but this is not explicitly stated.
26. Based on the above, we are not convinced that equations S65S77 are correct.
This part of the SI should be carefully revised and corrected as needed, or omitted.
Thank you for your comments. As suggested, we have chosen to omit this section because it was tangential to our main findings in this work. In the future we plan to revisit this part of this project and to make sure that we can accurately and convincingly address the reviewer’s concerns.
https://doi.org/10.7554/eLife.68693.sa2Article and author information
Author details
Funding
Cornell University
 Alireza Abbaspourrad
National Science Foundation (I20252)
 Alireza Abbaspourrad
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors would like to thank SS Suarez for fruitful discussions about the biological aspects of the work. We also thank DL Koch, CK Tung, J Fan, and M Esmaily for helpful discussions about hydrodynamic interactions, the dipole swimmer model, and the surface force analysis used here. The authors also thank our dearest KJ Donaghy for helping us edit this manuscript. This work was performed in part at the Cornell NanoScale Facility, a member of the National Nanotechnology Coordinated Infrastructure (NNCI), which is supported by the National Science Foundation (Grant NNCI2025233).
Senior Editor
 Anna Akhmanova, Utrecht University, Netherlands
Reviewing Editor
 Raymond E Goldstein, University of Cambridge, United Kingdom
Publication history
 Received: March 23, 2021
 Accepted: August 3, 2021
 Accepted Manuscript published: August 4, 2021 (version 1)
 Version of Record published: August 25, 2021 (version 2)
Copyright
© 2021, Zaferani 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,309
 Page views

 207
 Downloads

 5
 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

 Developmental Biology
 Stem Cells and Regenerative Medicine
Basal cells are multipotent stem cells of a variety of organs, including the respiratory tract, where they are major components of the airway epithelium. However, it remains unclear how diverse basal cells are, and how distinct subpopulations respond to airway challenges. Using single cell RNAsequencing and functional approaches, we report a significant and previously underappreciated degree of heterogeneity in the basal cell pool, leading to identification of six subpopulations in the adult murine trachea. Among these, we found two major subpopulations collectively comprising the most uncommitted of all the pool, but with distinct gene expression signatures. Notably, these occupy distinct ventral and dorsal tracheal niches and differ in their ability to selfrenew and initiate a program of differentiation in response to environmental perturbations in primary cultures and in mouse injury models in vivo. We found that such heterogeneity is acquired prenatally, when the basal cell pool and local niches are still being established, and depends on the integrity of these niches, as supported by the altered basal cell phenotype of tracheal cartilagedeficient mouse mutants. Lastly, we show that features that distinguish these progenitor subpopulations in murine airways are conserved in humans. Together, the data provide novel insights into the origin and impact of basal cell heterogeneity on the establishment of regionally distinct responses of the airway epithelium during injuryrepair and in disease conditions.

 Developmental Biology
 Stem Cells and Regenerative Medicine
Planarians have become an established model system to study regeneration and stem cells, but the regulatory elements in the genome remain almost entirely undescribed. Here, by integrating epigenetic and expression data we use multiple sources of evidence to predict enhancer elements active in the adult stem cell populations that drive regeneration. We have used ChIPseq data to identify genomic regions with histone modifications consistent with enhancer activity, and ATACseq data to identify accessible chromatin. Overlapping these signals allowed for the identification of a set of highconfidence candidate enhancers predicted to be active in planarian adult stem cells. These enhancers are enriched for predicted transcription factor (TF) binding sites for TFs and TF families expressed in planarian adult stem cells. Footprinting analyses provided further evidence that these potential TF binding sites are likely to be occupied in adult stem cells. We integrated these analyses to build testable hypotheses for the regulatory function of TFs in stem cells, both with respect to how pluripotency might be regulated, and to how lineage differentiation programs are controlled. We found that our predicted GRNs were independently supported by existing TF RNAi/RNAseq datasets, providing further evidence that our work predicts active enhancers that regulate adult stem cells and regenerative mechanisms.