ORIGINAL RESEARCH article

Front. Microbiol., 21 July 2021

Sec. Infectious Agents and Disease

Volume 12 - 2021 | https://doi.org/10.3389/fmicb.2021.715220

Clostridioides difficile Single Cell Swimming Strategy: A Novel Motility Pattern Regulated by Viscoelastic Properties of the Environment

  • Institute for Medical Microbiology and Virology, University Medical Center Göttingen, Göttingen, Germany

Abstract

Flagellar motility is important for the pathogenesis of many intestinal pathogens, allowing bacteria to move to their preferred ecological niche. Clostridioides difficile is currently the major cause for bacterial health care-associated intestinal infections in the western world. Most clinical strains produce peritrichous flagella and are motile in soft-agar. However, little knowledge exists on the C. difficile swimming behaviour and its regulation at the level of individual cells. We report here on the swimming strategy of C. difficile at the single cell level and its dependency on environmental parameters. A comprehensive analysis of motility parameters from several thousand bacteria was achieved with the aid of a recently developed bacterial tracking programme. C. difficile motility was found to be strongly dependent on the matrix elasticity of the medium. Long run phases of all four motile C. difficile clades were only observed in the presence of high molecular weight molecules such as polyvinylpyrrolidone (PVP) and mucin, which suggests an adaptation of the motility apparatus to the mucin-rich intestinal environment. Increasing mucin or PVP concentrations lead to longer and straighter runs with increased travelled distance per run and fewer turnarounds that result in a higher net displacement of the bacteria. The observed C. difficile swimming pattern under these conditions is characterised by bidirectional, alternating back and forth run phases, interrupted by a short stop without an apparent reorientation or tumbling phase. This motility type was not described before for peritrichous bacteria and is more similar to some previously described polar monotrichous bacteria.

Introduction

Clostridioides (formerly Clostridium) difficile is a spore forming, obligate anaerobic pathogen that causes CDI (C. difficile infection), which predominantly manifests as hospital-associated diarrhoea and pseudomembranous colitis (; ). Clinical symptoms are associated with toxin expression, in particular toxin A and B, which undergo a complex regulation pattern (; ; ). An intact gut microbiome is believed to protect from C. difficile infection. Dysbiosis however, for example after antibiotic treatment, favours C. difficile spore germination and subsequent colonisation (; ; ).

Flagellar motility and chemotaxis are important for successful colonisation and virulence of many gastrointestinal pathogens, for example Campylobacter jejuni, Salmonella enterica Serovar Typhimurium, Helicobacter pylori, and Vibrio chlolerae (; ; ; ). Most C. difficile strains produce peritrichous flagella, which can mediate swimming motility in soft-agar based assays (; ; ). The contribution of flagellar motility for the pathogenesis in mice was studied with the aid of C. difficile flagellar mutants, which were found to be reduced in their colonisation efficiency (; ). It is also known that the C. difficile genome contains a region which encodes for a complete set of chemotaxis genes (). Chemotaxis allows bacteria to swim up or down a chemical gradient and thus to find optimal growth conditions (). Regulation of motility by the chemotaxis system is well investigated for a variety of bacterial species. Chemical gradients are sensed by chemoreceptors, which transfer the signal via the adapter protein CheW and the histidine kinase CheA to the transducer CheY, which becomes phosphorylated (). The phosphorylated CheY finally interacts with the flagellar motor and leads to a modulation of motility characteristics.

However, in contrast to other gut pathogens, little knowledge exists on the motility of C. difficile at the single cell level. A careful analysis of the C. difficile swimming pattern and its dependency on environmental parameters would contribute to a better understanding of C. difficile motility for pathogenesis and particularly for the attachment and dissemination phases, in which flagella have previously been found to play a role (; ).

Gut pathogens like C. difficile, as well as commensalists, have to deal with the mucosal layer in the lower intestine. The colon mucus, made up mostly of Muc2, consists of a stratified inner, attached layer which is expressed by the surface goblet cells and serves as a protective barrier for the intestines (; ; ). The pore sizes of this layer allow for small molecule diffusion, but acts as a barrier to structures in the micrometre range, partially due to steric hindrance (). The mucus layer is partial degraded by commensalists and possesses then an increased pore size, allowing a better colonisation by bacteria (; ; ). As the layer is still cross-linked, it presumably still poses as a steric hindrance, influencing the rheological profile of the environment.

We report in this study on the swimming behaviour of C. difficile and its strong dependency on the viscoelastic properties of the medium across all motile clades. A comprehensive quantitative analysis of C. difficile motility parameters on the single cell level was obtained with the aid of the bacterial tracking programme YSMR (). We hypothesise that this dependency on viscoelastic properties is an adaptation to the properties of the mucin rich lower intestines, which form the habitat of C. difficile. A large-scale quantitative analysis of the swimming behaviour leads to the conclusion that the single cell motility displayed by C. difficile forms a novel pattern for peritrichous bacteria, which we describe here in detail.

Materials and Methods

Used C. difficile Strains

C. difficile 630 Δerm [Ribotype (RT) 012, DSM 28645, CP016318.1 ()], C. difficile R20291 (RT 027, DSM 27147, CP029423.1), DSM 100002 (RT 084), DSM 102978 (RT not determined, CP020380.1), DSM 28670 [RT SLO 237, CP012312.1 ()], and DSM 100005 (RT SLO 235).

Media and Strain Cultivation

C. difficile strains were grown at 37°C in BHIS (37 g/l brain heart infusion broth supplemented with 5 g/l yeast extract and 0.3 g/l cysteine) shaking at 180 rpm for liquid cultures, alternatively with 15 g/l agar for plates, or on Columbia agar with 5% sheep blood (COS, bioMérieux, Nürtingen, Germany). Cultivation was always performed under anaerobic conditions using a COY anaerobic gas chamber (COY Laboratory Products, Grass Lake, United States). The chamber was gas-flushed with 85% N2, 10% H2, and 5% CO2.

Motility With Polyvinylpyrrolidone (PVP) or Mucin

For experiments with high molecular weight polymers K 90 polyvinylpyrrolidone (MW 360,000 g/mol, Carl Roth, Karlsruhe, Germany, order nr. CP15.1) and Type I-S mucin (Bovine Submaxillary glands, Merck, Darmstadt, Germany, order nr. M3895-100MG) was used. From a PVP or Mucin stock solution (see Supplementary Text 2), the desired end-concentration of the additive as stated in the experiment was prepared in Dulbecco’s phosphate buffered saline (PBS, Merck, Darmstadt, Germany, order nr. D5652) to a volume of 90 μl. The solution was then incubated anaerobically at 37°C for at least 1 h. From a mid-exponential culture (OD600 0.4–0.6) 1 ml was centrifuged at 1,500 × g for 15 min using slow acceleration and deceleration ramps during centrifugation. The supernatant was removed, and the pellet resuspended to an OD600 of 6 in anaerobised BHIS. From the bacterial suspension, 10 μl were added to the PVP/Mucin solution for a final volume of 100 μl and OD600 of 0.6.

Video Microscopy

Inside the anaerobic chamber, 4 μl of the culture were placed on an objective slide, covered with a cover slip and immediately sealed with nail polish. Except when stated otherwise, slides were recorded 15 min after sealing on a microscope at room temperature with a 10x phase contrast objective (Nikon Eclipse TE2000-S, Nikon PlanFluor 10x). Videos were recorded for 3 min at 30 fps using an Aptina CMOS Sensor 18MP 1/2.3” Colour.

For 64x phase contrast objective movies, a Leica DMR microscope with a PL APO 506082 64x objective was used. Movies were taken with a Nikon D7100 camera at 30 fps.

Single Cell Tracking With YSMR

For quantitative data of the C. difficile swimming parameters, we recorded video files and always analysed 10 s of motility parameters per detected bacteria with the aid of the recently developed bacterial tracking programme YSMR (). We used a low magnification objective (10x) for video recording in order to monitor the swimming behaviour of a large number of bacteria (100–5,000) per video. All videos were analysed using YSMR v0.1.1 using the settings as set in the “Supplementary 1_tracking.txt” file.

Turnaround Calculation

Turnarounds were calculated using YSMR v0.1.1. Briefly, turnarounds were determined by calculating the difference in heading direction for two consecutive points. Heading directions were calculated by taking the arctan2 between the original point and one 0.3 s in the future. Turnarounds were calculated to be at positions with a relative angle change of at least 30° and an average minimum speed of at least 0.03 μm/s. The position with the greatest local angle difference was chosen as the turnaround.

Displacement Calculation

The maximum Euclidean distance between any two points reached during the analysed period of the track was calculated.

Calculations, Graphs and Video Editing

All performed calculations and graphs were created using Python 3.9.1, Seaborn 0.11.1 (), Matplotlib 3.4.1 (), SciPy 1.6.3 (), Numpy 1.20.0 (), and Pandas 1.2.1 (; ). Videos were edited using openCV-contrib-python 4.5.1.48 ().

Results

C. difficile Swimming Motility in BHIS Medium

When C. difficile 630 Δerm swimming behaviour was analysed in 100% BHIS medium or in a 10% BHIS 90% PBS mixture, the bacteria could be classified into two groups. A fraction of up to 75% was completely non-motile, whereas the motile fraction displayed an unusual motility phenotype, which is characterised by alternating, short, back and forth run phases of 0.5–3 bacterial lengths, corresponding to ∼3–15 μm (Supplementary Movie 1). The frequency of directional changes is very high with ∼18 turnarounds/10 s and we thus refer to this behaviour as “jitter motility.” None of the bacteria displayed prolonged run phases and there is no substantial net displacement of the bacteria during movement. The motility of further five C. difficile strains from different clades was analysed under the same conditions. The clade 5 strain RT078, which is non-flagellated and non-motile in soft-agar assays and lacks a flagellar operon, was as expected completely immotile. All other tested strains, namely 630 Δerm (clade 1), DSM 100002 (clade 1), R20291 (clade 2), DSM 102978 (clade 3), DSM 28670 (clade 4), and DSM 100005 (clade 4) displayed the jitter-motility phenotype in BHIS medium (Supplementary Movie 1). When C. difficile 630 Δerm motility was analysed in 100% PBS instead of medium containing BHIS, the jitter-motility fraction was absent and bacteria were completely non-motile (data not shown).

The Motility Phenotype Changes in the Presence of the High Molecular Weight Polymers PVP and Mucin

It is well investigated that the motility of bacteria depends on the physical properties of the medium, particularly on the viscosity and on the matrix elasticity (; ). C. difficile 630 Δerm cells were exposed to various glycerol concentrations from 0.05 to 25% in BHIS medium in order to increase the viscosity. None of these conditions induced increased displacement and the bacteria remained in their jitter-motility stage.

As C. difficile displays motility in agarose (; ), we wondered if altering the matrix elasticity of the medium has an influence on the motility phenotype. High molecular weight polymers such as polyvinylpyrrolidone (PVP), which increase the matrix elasticity, were described to enhance motility in various bacteria in a concentration dependent manner (; ). When C. difficile was resuspended in a solution containing 3.6% PVP, bacteria displayed a dramatic change in their motility phenotype with long, alternating back and forth runs (Supplementary Movies 2, 3). In this setup, motility could be observed for at least 2 h. When motility was analysed in 3.6% PVP dissolved in pure PBS without any BHIS, cells were completely immotile, suggesting that a nutrient source is required for the observed motility phenotype. We also tested the combination of 3.6% PVP and 10% glycerol in 10% BHIS. The addition of glycerol had no influence on the motility behaviour, indicating that glycerol did not inhibit the PVP-induced motility phenotype with long, alternating back and forth runs.

To exclude the possibility that the observed motility patterns are influenced by the relatively small distance between the microscope slide and the cover slide or the wall effect (), we monitored the swimming behaviour also on cavity slides. The same motility pattern as with the original setup was observed for experiments with PVP (Supplementary Movie 9).

In order to investigate how common the observed motility phenotype is in Clostridiodes, we analyzed representative members from all four motile clades (). The PVP-induced change of the motility phenotype occurred in all motile strains that were tested, including 630 Δerm (clade 1, Supplementary Movies 47), DSM 100002 (clade 1), R20291 (clade 2), DSM 102978 (clade 3), DSM 28670 (clade 4), and DSM 100005 (clade 4) (Supplementary Movie 8).

We applied single cell tracking using the recently developed bacterial tracking programme YSMR () to analyse and quantify the motility characteristics of C. difficile (see Supplementary Videos 57). Per track we quantified the average speed, maximal displacement to body length ratio, arc-chord ratio, and calculated the turnarounds. The displacement/bacterial body length ratio indicates how much ground a bacterium covered in relation to its own length. The arc-chord ratio, calculated by dividing the length between start and end point by the total distance travelled, denotes the tortuosity of the bacterial movement, with a straight travel path resulting in an arc-chord ratio of 1.

To provide a graphical overview of the PVP-induced motility change, the track of each bacterium was displayed in a rose plot for an experiment with PVP concentrations in the range between 0.4 and 6.4% (Figure 1). More than 1,000 tracks were analysed for each PVP concentration. Each track is displayed with its starting position set to the origin of the coordinate system (0, 0). Additionally, each tracks maximum distance is indicated by its colour. For each track three selected characteristics are displayed in a bivariate scatterplot (Figure 2). The arc-chord ratio and average speed of each track are each plotted against the displacement/bacterial length ratio, with individual kernel density estimates for each condition. The experiment was repeated twice on different days and for each biological replicate independent rose and scatter plots were generated (Supplementary Figures 14).

FIGURE 1

FIGURE 2

At 0.4 and 0.8% PVP the average speed of the bacteria increases (median values 3.1–4.7 μm/s and 5.84–8.88 μm/s, respectively, Table 1) with a corresponding, albeit small increase in displacement/body length ratio (Figures 1A,B, 2). The motility phenotype at both concentrations is likely still futile, as the bacteria do not achieve notable displacement, with median displacement to length ratios below 1.5 for 0.4% and 2.5 for 0.8% PVP (Figures 1A,B and Table 1). This also is noticeable in the arc-chord ratios, which are on median between 0.08 and 0.11 (Figure 2A and Table 1). None of the bacteria revealed run phases of greater than 10 s, as can be seen in samples with higher PVP concentrations.

TABLE 1

ConditionMedian displacement/Bacterial length ratioMedian distance (μm)Median speed (μm/s)Median arc-chord ratio
0.4% K90 PVP (n = 1,271, 1,300, 1,188)1.40, 1.38, 0.8046.69, 46.96, 30.994.67, 4.70, 3.100.11, 0.10, 0.09
0.8% K90 PVP (n = 1,369, 2,486, 2,251)2.35, 2.41, 1.4874.84, 88.77, 58.367.48, 8.88, 5.840.09, 0.09, 0.08
1.6% K90 PVP (n = 4,326, 2,448, 2,034)3.77, 1.26, 2.2492.67, 38.20, 55.709.27, 3.82, 5.570.12, 0.10, 0.13
3.2% K90 PVP (n = 3,986, 2,626, 2,208)1.15, 16.57, 12.8521.24, 109.04, 93.012.12, 10.90, 9.300.22, 0.90, 0.92
3.6% K90 PVP (n = 2,580, 2,416, 2,447)21.67, 6.29, 5.60126.28, 58.68, 56.8112.63, 5.87, 5.680.90, 0.47, 0.44
6.4% K90 PVP (n = 2,221, 1,454, 328)16.80, 0.16, 0.0790.50, 7.70, 5.159.05, 0.77, 0.510.97, 0.06, 0.03

Median values for select properties for different PVP concentrations.

Values correspond to Figure 1 and Supplementary Figures 1, 3 (first, second, and third given value, respectively). For a full set, including average, 25 and 75 percentiles for each property, as well as distance, displacement, and bacterial length see Supplementary Table 1.

At 1.6% PVP the first visible shift in the motility phenotype occurs. Several individual bacteria achieve clear displacement (Figure 1C), indicating that the increase in PVP starts to form a medium in which C. difficile can be motile. Speed at 1.6% PVP is markedly increased in comparison to lower PVP concentrations (Figure 2B). On the population level though, bacteria still have a median displacement/bacterial length ratio between 1.3 and 3.8, and median arc-chord ratios between 0.1 and 0.13 (Table 1). The fraction of bacteria with continuous runs of at least 10 s was between 13 and 58%.

At PVP concentrations of 3.2% and higher, the majority of C. difficile cells display clear motility (Figures 1C–F, 2). Instead of the futile short movements, long, sustained run phases with varying curvature and intermittent turnarounds are the majority. They lead to a motility phenotype that can achieve displacement from the starting position. The motility phenotype is achieving an optimum in terms of displacement/bacterial length ratio, speed, and arc-chord ratio at around 3.2–3.6% (w/v) PVP (Figures 1D,E, 2 and Supplementary Figures 14). The fraction of bacteria with continuous run lengths of at least 10 s was between 67 and 86%.

At a concentration of 6.4% PVP the proportion of motile bacteria decreases in the majority of experiments (Figures 1F, 2A and Supplementary Figures 14). Bacteria that are still motile display very straight runs with little curvature and a median arc-chord ratio of up to 0.97 (Table 1). The fraction of bacteria showing continuous runs of at least 10 s was between 95 and 100%. The speed of motile bacteria is also markedly lower than with 3.6% PVP (Figure 2 and Supplementary Figure 2).

C. difficile inhabits areas with high molecular weight polymers in the larger intestine such as mucins, which cover the epithelial cell layer with a gel like structure (; ). We used bovine submaxillary gland mucin in concentrations ranging from 2.5 to 40 mg/ml in doubling steps (Supplementary Movie 3). The track of each bacterium was displayed in a rose plot from an experiment with mucin concentrations between 2.5 and 40 mg/ml (Supplementary Figures 79). At concentrations of 2.5–10 mg/ml, typical short back and forth motion was observed as in medium without PVP or at low PVP concentrations (Supplementary Figures 79A–C). However, at mucin concentrations of 20 and 40 mg/ml bacteria displays long, alternating back and forth run phases with few turning points and high net displacement (Supplementary Figures 79D,E).

C. difficile Displays a Distinct Motility Phenotype

Representative tracks from a video (Supplementary Movie 4) of C. difficile 630 Δerm, as identified by the YSMR software, are shown in Figure 3A. The swimming motility phenotype in PVP is characterised by alternating back and forth run phases, which are indistinguishable in length and duration. There is a short stop between two run phases, but there is neither an obvious tumbling nor a sudden reorientation during the stop (Figures 3A,B). Nevertheless, the reverse run is not exactly mirroring the previous forward run, but their curvatures can be different. This leads to a significant net displacement of the bacterium during movement.

FIGURE 3

A manual inspection of bacterial tracks revealed no obvious difference between a forward and a backward run phase. In addition, we performed a thorough quantitative comparison of two motility parameter, speed and arc-chored ratio, before and after a stop phase. We calculated the difference in average speed and in average arc-chored ratio before and after the stop for > 300 of these events and plotted the obtained difference values (Figure 4). If back and forth runs possess non-identical motility parameter, a histogram of difference values would show two peaks. In contrast, the histogram shows a typical Gaussian distribution, which is indicative for identical speed values and arc-chord ratios before and after a stop. This confirms that the two run phases between a stop are indistinguishable in their motility parameters.

FIGURE 4

In summary, this characterisation of the C. difficile swimming phenotype reveals a novel motility type for peritrichous bacteria.

Discussion

The Motility of C. difficile Is Dependent on Matrix Elasticity

We used single cell video microscopy in conjunction with a software-based quantification of motility parameters to characterise the motility profile of C. difficile under various conditions. The applied YSMR software allowed the unbiased and automated analysis of a large number of tracks per frame. It is well investigated that the motility of bacteria depends on the physical properties of the medium in which they move (; ; ; ). When analysed in regular, water-based media without increased viscosity or matrix elasticity, C. difficile displayed a seemingly useless motility phenotype. Moving rapidly back and forth while achieving no notable displacement, the motility seemed to be highly energy intensive, yet ineffectual. It is unclear whether this motility phenotype has any relevance for C. difficile physiology and it is possible that it does not occur under natural environmental conditions.

Important physical characteristics shaping the environment include the viscosity, as well as the matrix elasticity of the medium (; ). Viscosity is generally accepted to be a major criterion for bacterial movement, as it plays an overwhelming role in the physics governing motility at the Reynolds number at which bacteria exist (). Increasing the medium viscosity by adding glycerol was therefore a logical first step to study motility. However, glycerol, which increases medium viscosity without affecting the matrix elasticity, did not strongly affect the swimming behaviour of C. difficile and effective motility could not be observed. In contrast, increase of matrix elasticity by PVP, a long-chained polymer at 360 MDa, or by mucins, both of which can form gel like structures, has a dramatic, concentration dependent effect on the swimming motility of C. difficile.

At PVP concentrations of 1.6%, a marked change in motility behaviour becomes visible with long, sustained run phases with varying curvature and intermittent turnarounds starting to appear that result in a net displacement of the bacterium. An optimum in terms of displacement/bacterial length ratio, speed, and arc-chord ratio is achieved at 3.2–3.6% (w/v) PVP. At 6.4% a decrease or complete cessation of movement can be observed, in which the strongly reduced number of motile cells perform very straight runs.

Our findings lead to the conclusion that the motility apparatus in C. difficile, across all clades, is adapted to an environment of large, high molarity molecules, which increase the matrix elasticity. Such molecules can be found in the natural environment of C. difficile, for example in the mucin-rich lower intestines. Association and interaction of C. difficile with the mucosal layer has been researched previously. Its ability to adhere to, and survive in the mucosal layer, as well as its demonstrated presence in the outer mucosal layer during infection, further indicates that the C. difficile motility apparatus is specialised for this environment (; ; ). It should be noted that Mucin types differ in their protein and glycosylation patterns. The bovine type I-S mucin used here is different from the human gastrointestinal mucin.

C. difficile Displays a New Swimming Strategy

Depending on the species and the number and location of flagella across the cell body, bacteria use different swimming strategies. At least four different swimming types are described in detail in the literature, but none of them matches the observed C. difficile swimming strategy. The swimming pattern of the peritrichous model organism E. coli and its regulation are particularly well investigated. E. coli switches between a “run” mode and a “tumble” mode (). During the run mode, the counter-clockwise rotating flagella form a bundle that pushes the cell forward, leading to a straight swimming phase. In the “tumble” mode the clockwise rotation of flagella leads to a breakup of the flagellar bundle, followed by a random change of the cell orientation. Bacteria with a single flagellum, such as Vibrio alginolyticus and Rhodobacter sphaeroides, possess different swimming patterns. V. alginolyticus displays a three-step pattern, in which a forward swimming phase is followed by a reverse swimming phase and a reorientation (“flick”) phase (; ). The flagellum of R. sphaeroides rotates in only one direction and the motility pattern follows a “stop and coil” pattern (). Random reorientation occurs when the flagellar rotation stops. Pseudomonas aeruginosa performs a run-reverse-pause pattern, where reversals are marked with ∼180° angles. After pauses, the swimming direction does generally not change. Forward and backward travel directions do not differ and are achieved by clockwise or counter-clockwise turning of the single flagellum (). Recently, a similar swimming strategy was identified in the singly flagellated proteobacterium Caulobacter crescentus, suggesting that the forward and backward motility pattern is more widely spread than previously expected (). In fact, the “run and tumble” strategy of E. coli has not been identified for any polar flagellated bacterium (). Pseudomonas putida produces a tuft of flagella on one cell pole. It can switch between a pushing mode, a pulling mode and a wrapping mode, with reorientation taking place during switches in modes ().

The PVP-induced C. difficile motility pattern appears to be markedly different from those previously described. In particular, an obvious reorientation phase caused by a “tumble” or a “flick” appears to be absent. Instead, C. difficile stops in place at the end of a run phase and performs a ∼180° turnaround by changing the movement direction from forward to backward. Turnarounds therefore lead to a direct reversal in direction, but do not play an immediate role in the random reorientation process. While at a first glance this is similar to the polar monotrichous Pseudomonas aeruginosa, C. difficile has peritrichous flagella, which leads to the question how the reversal is achieved (; ). At the beginning of a backward run C. difficile travel back on a very similar path as the previous forward run. However, the curvature can change during the course of the backward run, leading to a significant net displacement of the bacterium during movement.

Generally, flagella have an energetically favourable turning direction (). In some species a difference in torque profile based on rotational direction, likely due to stator-rotor interactions, could be shown, which leads to different motility parameter for clockwise/counter-clockwise rotation (; ). The energetically favourable turning direction can translate into a preferred direction of travel. In C. difficile, forward and backward runs are indistinguishable, suggesting that C. difficile has no preferred direction of travel. C. difficile is able to traverse bi-directional, which appears to be a novel mode of motility for peritrichous bacteria (). These results suggest that matrix elasticity affects the C. difficile flagellar motor rotation state with a strong reduction of switching frequency in an environment of high matrix elasticity. Our findings lead to the conclusion that the motility apparatus in C. difficile is adapted to an environment of large, high molarity molecules with the ability for bi-directional motility.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Author contributions

JS, WB, and IO: conceptualisation, methodology, and validation. JS: software, formal analysis, data curation, and visualisation. JS and IO: investigation. UG, AZ, and WB: resources, supervision, and project administration. JS and WB: writing—original draft preparation. JS, UG, AZ, and WB: writing—review and editing. UG: funding acquisition. All authors contributed to the article and approved the submitted version.

Funding

This work was funded by the Federal State of Lower Saxony, Niedersächsisches Vorab (VWZN2889/3215/3266).

Acknowledgments

We would like to thank Dr. Joachim Schwanbeck for kindly loaning us the camera with which all images and videos were taken.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.715220/full#supplementary-material

Supplementary Figure 1

Single cell motility dependency on PVP of C. difficile 630 Δerm. The motility of C. difficile cells was tracked in the presence of six increasing K90 PVP concentrations (w/v) from 0.4 to 6.4% (A–F). All observation durations were limited to 10 s for better comparability. The number of displayed tracks is stated in the title of each experiment. The starting position of each track is at the origin of the plot at position 0 μm/0 μm. The scale is identical for all six plots. Track colours denote the total travelled distance, with a colour bar on the right-hand side as a legend. A 5 μm scale bar is depicted in the lower left of subfigure (A). The median cell lengths of C. difficile was determined to be between 4.81 and 5.28 μm (Supplementary Table 1).

Supplementary Figure 2

Direct comparison of motility characteristics and statistical distribution within experiments for different PVP concentrations. The motility of C. difficile was tracked in the presence of six increasing K90 polyvinylpyrrolidone (PVP) concentrations. The PVP concentration starts at 0.4% and increase at a doubling rate up to 6.4%, with an additional step of 3.6% PVP. All observations were limited to 10 s for better comparability. Average speed (A) and arc-chord ratio (B, denoting tortuosity of path) are plotted against the displacement/bacterial length ratio (showing the maximal distance in body lengths each bacterium has covered) for each tracked bacterium. The kernel density estimate is plotted independently for each experiment and variable, denoting the distribution for each characteristic.

Supplementary Figure 3

Single cell motility dependency on PVP of C. difficile 630 Δerm. The motility of C. difficile cells was tracked in the presence of six increasing K90 PVP concentrations (w/v) from 0.4–6.4% (A–F). All observation durations were limited to 10 s for better comparability. The number of displayed tracks is stated in the title of each experiment. The starting position of each track is at the origin of the plot at position 0 μm/0 μm. The scale is identical for all six plots. Track colours denote the total travelled distance, with a colour bar on the right-hand side as a legend. A 5 μm scale bar is depicted in the lower left of subfigure (A). The median cell lengths of C. difficile was determined to be between 4.99 and 5.74 μm (Supplementary Table 1).

Supplementary Figure 4

Direct comparison of motility characteristics and statistical distribution within experiments for different PVP concentrations. The motility of C. difficile was tracked in the presence of six increasing K90 polyvinylpyrrolidone (PVP) concentrations. The PVP concentration starts at 0.4% and increase at a doubling rate up to 6.4%, with an additional step of 3.6% PVP. All observations were limited to 10 s for better comparability. Average speed (A) and arc-chord ratio (B, denoting tortuosity of path) are plotted against the displacement/bacterial length ratio (showing the maximal distance in body lengths each bacterium has covered) for each tracked bacterium. The kernel density estimate is plotted independently for each experiment and variable, denoting the distribution for each characteristic.

Supplementary Figure 5

Comparison of speed and arc-chord ratio before and after a turnaround. Events with more than one calculated turnaround within a 4 s window were excluded. Turnarounds are set to the point with local greatest change in angle. Tracks are from Figure 1C. (A) speed profile of 30 example tracks before and after a turnaround (black line, T = 0). (B) Five additional detailed speed graphs before and after turnaround (black line, T = 0). The colouration indicates the relative angle change between positions. The speed maximum at or up to 0.5 s after the turnaround ranges between 7.6 and 19.3 μm/s. (C,D) Histograms and corresponding box plots on the speed delta (C) and arc-chord ratio delta (D) before and after turnaround. The deltas are calculated between the time intervals [-2 s, -1 s] and [1 s, 2 s] (grey regions in A) around the turnaround.

Supplementary Figure 6

Comparison of speed and arc-chord ratio before and after a turnaround. Events with more than one calculated turnaround within a 4 s window were excluded. Turnarounds are set to the point with local greatest change in angle. Tracks are from Figure 1C. (A) speed profile of 30 example tracks before and after a turnaround (black line, T = 0). (B) Five additional detailed speed graphs before and after turnaround (black line, T = 0). The colouration indicates the relative angle change between positions. The speed maximum at or up to 0.5 s after the turnaround ranges between 2.5 and 20.0 μm/s. (C,D) Histograms and corresponding box plots on the speed delta (C) and arc-chord ratio delta (D) before and after turnaround. The deltas are calculated between the time intervals [-2 s, -1 s] and [1 s, 2 s] (grey regions in A) around the turnaround.

Supplementary Figure 7

Single cell motility dependency on mucin concentration of C. difficile 630 Δerm. The motility of C. difficile cells was tracked in the presence of 5 increasing mucin concentrations. The mucin concentration starts at 2.5 mg/ml and increases in doubling steps up to 40 mg/ml (A–E). The number of observations is stated in the title of each experiment. All observation durations were limited to 10 s for better comparability. The starting position of single C. difficile tracks are centred on the origin. The scale is identical for all six plots. Track colours denote the total travelled distance, with a colour bar on the right-hand side as a legend. A 5 μm scale bar is depicted in the lower left of subfigure (A).

Supplementary Figure 8

Single cell motility dependency on mucin concentration of C. difficile 630 Δerm. The motility of C. difficile cells was tracked in the presence of 5 increasing mucin concentrations. The mucin concentration starts at 2.5 mg/ml and increases in doubling steps up to 40 mg/ml (A–E). The number of observations is stated in the title of each experiment. All observation durations were limited to 10 s for better comparability. The starting position of single C. difficile tracks are centred on the origin. The scale is identical for all six plots. Track colours denote the total travelled distance, with a colour bar on the right-hand side as a legend. A 5 μm scale bar is depicted in the lower left of subfigure (A).

Supplementary Figure 9

Single cell motility dependency on mucin concentration of C. difficile 630 Δerm. The motility of C. difficile cells was tracked in the presence of 5 increasing mucin concentrations. The mucin concentration starts at 2.5 mg/ml and increases in doubling steps up to 40 mg/ml (A–E). The number of observations is stated in the title of each experiment. All observation durations were limited to 10 s for better comparability. The starting position of single C. difficile tracks are centred on the origin. The scale is identical for all six plots. Track colours denote the total travelled distance, with a colour bar on the right-hand side as a legend. A 5 μm scale bar is depicted in the lower left of subfigure (A).

Supplementary Table 1

Various statistics on individual experiments. Values are given as: average (25 % quantile, median, 75 % quantile). N: number of observations.

Supplementary Movie 1

Motility of various C. difficile strains in BHIS without PVP. Video was taken with a 10x objective and converted to inverted grey-scale while filming. Shown are 630 Δerm (clade 1), DSM 100002 (clade 1), R20291 (clade 2), DSM 102978 (clade 3), DSM 28670 (clade 4), and DSM 100005 (clade 4) in 100% BHIS.

Supplementary Movie 2

Motility in 10% BHIS/PBS with 3.6% PVP. The video of C. difficile 630 Δerm was taken with a 64x objective.

Supplementary Movie 3

Motility in 10% BHIS/PBS with 2.5–40 mg/ml Mucin. Shown are excerpts of movies used in Supplementary Figure 7. The video of C. difficile 630 Δerm was taken with a 10x objective and converted to inverted grey-scale while filming.

Supplementary Movie 4

Motility in 10% BHIS/PBS with 1.8% PVP and annotated tracks. The video of C. difficile 630 Δerm was taken with a 10x objective and converted to inverted grey-scale while filming. Examples of tracks are depicted in Figure 3.

Supplementary Movie 5

Parts of Movies used in Figures 1, 2. The videos of C. difficile 630 Δerm were taken with a 10x objective and converted to inverted grey-scale while filming. Shown are excerpts of movies of C. difficile 630 Δerm in increasing PVP concentrations.

Supplementary Movie 6

Parts of Movies used in Supplementary Figures 1, 2. The videos of C. difficile 630 Δerm were taken with a 10x objective and converted to inverted grey-scale while filming. Shown are excerpts of movies of C. difficile 630 Δerm in increasing PVP concentrations.

Supplementary Movie 7

Parts of Movies used in Supplementary Figures 3, 4. The videos of C. difficile 630 Δerm were taken with a 10x objective and converted to inverted grey-scale while filming. Shown are excerpts of movies of C. difficile 630 Δerm in increasing PVP concentrations.

Supplementary Movie 8

Motility of various C. difficile strains in BHIS with PVP. The videos were taken with a 10x objective and converted to inverted grey-scale while filming. Shown are DSM 100002 (clade 1), R20291 (clade 2), DSM 102978 (clade 3), DSM 28670 (clade 4), and DSM 100005 (clade 4) with 2% (w/v) PVP.

Supplementary Text 1

Settings-file used for YSMR. “Supplementary_text 1_tracking.txt”.

Supplementary Text 2

Mucin/PVP handling and storage.

References

  • 1

    Anjuwon-FosterB. R.TamayoR. (2017). A genetic switch controls the production of flagella and toxins in Clostridium difficile.PLoS Genet.13:e1006701. 10.1371/journal.pgen.1006701

  • 2

    AntaniJ. D.SumaliA. X.LeleT. P.LeleP. P. (2021). Asymmetric random walks reveal that the chemotaxis network modulates flagellar rotational bias in Helicobacter pylori.eLife10:e63936. 10.7554/eLife.63936

  • 3

    ArmitageJ. P.MacnabR. M. (1987). Unidirectional, intermittent rotation of the flagellum of Rhodobacter sphaeroides.J. Bacteriol.169514518. 10.1128/JB.169.2.514-518.1987

  • 4

    BabanS. T.KuehneS. A.Barketi-KlaiA.CartmanS. T.KellyM. L.HardieK. R.et al (2013). The Role of Flagella in Clostridium difficile Pathogenesis: Comparison between a Non-Epidemic and an Epidemic Strain.PLoS One8:e73026. 10.1371/journal.pone.0073026

  • 5

    BartlettT. M.BrattonB. P.DuvshaniA.MiguelA.ShengY.MartinN. R.et al (2017). A periplasmic polymer curves Vibrio cholerae and promotes pathogenesis.Cell168172.e185.e. 10.1016/j.cell.2016.12.019

  • 6

    BatahJ.KobeissyH.Bui PhamP. T.Denève-LarrazetC.KuehneS.CollignonA.et al (2017). Clostridium difficile flagella induce a pro-inflammatory response in intestinal epithelium of mice in cooperation with toxins.Sci. Rep.7:3256. 10.1038/s41598-017-03621-z

  • 7

    BoinM. A.AustinM. J.Hã¤seC. C. (2004). Chemotaxis in Vibrio cholerae.FEMS Microbiol. Lett.23918. 10.1016/j.femsle.2004.08.039

  • 8

    BradskiG. (2000). The OpenCV Library.Dr. Dobbs J. Softw. Tools120122125.

  • 9

    BuffieC. G.PamerE. G. (2013). Microbiota-mediated colonization resistance against intestinal pathogens.Nat. Rev. Immunol.13790801. 10.1038/nri3535

  • 10

    CaiQ.LiZ.OuyangQ.LuoC.GordonV. D. (2016). Singly Flagellated Pseudomonas aeruginosa Chemotaxes Efficiently by Unbiased Motor Regulation.mBio7e00013e16. 10.1128/mBio.00013-16

  • 11

    ChandrasekaranR.LacyD. B. (2017). The role of toxins in Clostridium difficile infection.FEMS Microbiol. Rev.41723750. 10.1093/femsre/fux048

  • 12

    CoursonD. S.PokhrelA.ScottC.MadrillM.RineholdA. J.TamayoR.et al (2019). Single cell analysis of nutrient regulation of Clostridioides (Clostridium) difficile motility.Anaerobe59205211. 10.1016/j.anaerobe.2019.102080

  • 13

    DannheimH.RiedelT.Neumann-SchaalM.BunkB.SchoberI.SpröerC.et al (2017). Manual curation and reannotation of the genomes of Clostridium difficile 630Δerm and C. difficile 630.J. Med. Microbiol.66286293. 10.1099/jmm.0.000427

  • 14

    GaoY.NeubauerM.YangA.JohnsonN.MorseM.LiG.et al (2014). Altered motility of Caulobacter crescentus in viscous and viscoelastic media.BMC Microbiol.14:322. 10.1186/s12866-014-0322-3

  • 15

    GómezS.ChavesF.OrellanaM. A. (2017). Clinical, epidemiological and microbiological characteristics of relapse and re-infection in Clostridium difficile infection.Anaerobe48147151. 10.1016/j.anaerobe.2017.08.012

  • 16

    GreenbergE. P.Canale-ParolaE. (1977). Motility of flagellated bacteria in viscous environments.J. Bacteriol.132356358. 10.1128/JB.132.1.356-358.1977

  • 17

    GriffithsD.FawleyW.KachrimanidouM.BowdenR.CrookD. W.FungR.et al (2010). Multilocus Sequence Typing of Clostridium difficile.J. Clin. Microbiol.48770778. 10.1128/JCM.01796-09

  • 18

    GrognotM.TauteK. M. (2021a). A multiscale 3D chemotaxis assay reveals bacterial navigation mechanisms.Commun. Biol.4:669. 10.1038/s42003-021-02190-2

  • 19

    GrognotM.TauteK. M. (2021b). More than propellers: how flagella shape bacterial motility behaviors.Curr. Opin. Microbiol.617381. 10.1016/j.mib.2021.02.005

  • 20

    HarrisC. R.MillmanK. J.van der WaltS. J.GommersR.VirtanenP.CournapeauD.et al (2020). Array programming with NumPy.Nature585357362. 10.1038/s41586-020-2649-2

  • 21

    HintscheM.WaljorV.GroßmannR.KühnM. J.ThormannK. M.PeruaniF.et al (2017). A polar bundle of flagella can drive bacterial swimming by pushing, pulling, or coiling around the cell body.Sci. Rep.7:16771. 10.1038/s41598-017-16428-9

  • 22

    HunterJ. D. (2007). Matplotlib: A 2D Graphics Environment.Comput. Sci. Eng.99095. 10.1109/MCSE.2007.55

  • 23

    JohanssonM. E. V.LarssonJ. M. H.HanssonG. C. (2011). The two mucus layers of colon are organized by the MUC2 mucin, whereas the outer layer is a legislator of host-microbial interactions.Proc. Natl. Acad. Sci.10846594665. 10.1073/pnas.1006451107

  • 24

    JustI.SelzerJ.WilmM.Eichel-StreiberC.von, MannM.et al (1995). Glucosylation of Rho proteins by Clostridium difficile toxin B.Nature375500503. 10.1038/375500a0

  • 25

    KorolikV. (2019). The role of chemotaxis during Campylobacter jejuni colonisation and pathogenesis.Curr. Opin. Microbiol.473237. 10.1016/j.mib.2018.11.001

  • 26

    KuehneS. A.CartmanS. T.HeapJ. T.KellyM. L.CockayneA.MintonN. P. (2010). The role of toxin A and toxin B in Clostridium difficile infection.Nature467711713. 10.1038/nature09397

  • 27

    LaiS. K.WangY.-Y.WirtzD.HanesJ. (2009). Micro- and macrorheology of mucus.Adv. Drug Deliv. Rev.6186100. 10.1016/j.addr.2008.09.012

  • 28

    LawsonP. A.CitronD. M.TyrrellK. L.FinegoldS. M. (2016). Reclassification of Clostridium difficile as Clostridioides difficile (Hall and O’Toole 1935) Prévot 1938.Anaerobe409599. 10.1016/j.anaerobe.2016.06.008

  • 29

    LefflerD. A.LamontJ. T. (2015). Clostridium difficile Infection.N. Engl. J. Med.37215391548. 10.1056/NEJMra1403772

  • 30

    LertsethtakarnP.OttemannK. M.HendrixsonD. R. (2011). Motility and Chemotaxis in Campylobacter and Helicobacter.Annu. Rev. Microbiol.65389410. 10.1146/annurev-micro-090110-102908

  • 31

    LessaF. C.MuY.BambergW. M.BeldavsZ. G.DumyatiG. K.DunnJ. R.et al (2015). Burden of Clostridium difficile Infection in the United States.N. Engl. J. Med.372825834. 10.1056/NEJMoa1408913

  • 32

    MagariyamaY.IchibaM.NakataK.BabaK.OhtaniT.KudoS.et al (2005). Difference in Bacterial Motion between Forward and Backward Swimming Caused by the Wall Effect.Biophys. J.8836483658. 10.1529/biophysj.104.054049

  • 33

    MartinezV. A.Schwarz-LinekJ.ReuferM.WilsonL. G.MorozovA. N.PoonW. C. K. (2014). Flagellated bacterial motility in polymer solutions.Proc. Natl. Acad. Sci.1111777117776. 10.1073/pnas.1415460111

  • 34

    MatillaM. A.KrellT. (2018). The effect of bacterial chemotaxis on host infection and pathogenicity.FEMS Microbiol. Rev.42:fux052. 10.1093/femsre/fux052

  • 35

    McGuckinM. A.LindénS. K.SuttonP.FlorinT. H. (2011). Mucin dynamics and enteric pathogens.Nat. Rev. Microbiol.9265278. 10.1038/nrmicro2538

  • 36

    McKinneyW. (2010). “Data Structures for Statistical Computing in Python,” in Proceedings of the Python in Science Conference, (Austin: Python in Science Conference), 10.25080/Majora-92bf1922-00a

  • 37

    MinaminoT.KinoshitaM.NambaK. (2019). Directional Switching Mechanism of the Bacterial Flagellar Motor.Comput. Struct. Biotechnol. J.1710751081. 10.1016/j.csbj.2019.07.020

  • 38

    PorterS. L.WadhamsG. H.ArmitageJ. P. (2011). Signal processing in complex chemotaxis pathways.Nat. Rev. Microbiol.9153165. 10.1038/nrmicro2505

  • 39

    PurcellE. M. (1977). Life at low Reynolds number.Am. J. Phys.45311. 10.1119/1.10903

  • 40

    RebackJ.McKinneyW.Jbrockmendel, BosscheJ. V. D.AugspurgerT.CloudP.et al (2021). pandas-dev/pandas: Pandas 1.2.1.San Francisco: GitHub, 10.5281/ZENODO.4452601

  • 41

    RiedelT.WetzelD.HofmannJ. D.PlorinS. P. E. O.DannheimH.BergesM.et al (2017). High metabolic versatility of different toxigenic and non-toxigenic Clostridioides difficile isolates.Int. J. Med. Microbiol.307311320. 10.1016/j.ijmm.2017.05.007

  • 42

    SarkarM. K.PaulK.BlairD. (2010). Chemotaxis signaling protein CheY binds to the rotor protein FliN to control the direction of flagellar rotation in Escherichia coli.Proc. Natl. Acad. Sci.10793709375. 10.1073/pnas.1000935107

  • 43

    SchneiderW. R.DoetschR. N. (1974). Effect of Viscosity on Bacterial Motility.J. Bacteriol.117696701. 10.1128/JB.117.2.696-701.1974

  • 44

    SchwanbeckJ.OehmigI.DretzkeJ.ZautnerA. E.GroßU.BohneW. (2020). YSMR: a video tracking and analysis program for bacterial motility.BMC Bioinformatics21:166. 10.1186/s12859-020-3495-9

  • 45

    SemenyukE. G.PoroykoV. A.JohnstonP. F.JonesS. E.KnightK. L.GerdingD. N.et al (2015). Analysis of Bacterial Communities during Clostridium difficile Infection in the Mouse.Infect. Immun.8343834391. 10.1128/IAI.00145-15

  • 46

    SoavelomandrosoA. P.GaudinF.HoysS.NicolasV.VedantamG.JanoirC.et al (2017). Biofilm Structures in a Mono-Associated Mouse Model of Clostridium difficile Infection.Front. Microbiol.8:2086. 10.3389/fmicb.2017.02086

  • 47

    StecherB.HapfelmeierS.MüllerC.KremerM.StallmachT.HardtW.-D. (2004). Flagella and Chemotaxis Are Required for Efficient Induction of Salmonella enterica Serovar Typhimurium Colitis in Streptomycin-Pretreated Mice.Infect. Immun.7241384150. 10.1128/IAI.72.7.4138-4150.2004

  • 48

    StockerR. (2011). Reverse and flick: Hybrid locomotion in bacteria.Proc. Natl. Acad. Sci.10826352636. 10.1073/pnas.1019199108

  • 49

    TasteyreA.BarcM.-C.CollignonA.BoureauH.KarjalainenT. (2001). Role of FliC and FliD Flagellar Proteins of Clostridium difficile in Adherence and Gut Colonization.Infect. Immun.6979377940. 10.1128/IAI.69.12.7937-7940.2001

  • 50

    TaylorJ. A.SichelS. R.SalamaN. R. (2019). Bent Bacteria: A Comparison of Cell Shape Mechanisms in Proteobacteria.Annu. Rev. Microbiol.73457480. 10.1146/annurev-micro-020518-115919

  • 51

    TheriotC. M.YoungV. B. (2015). Interactions Between the Gastrointestinal Microbiome and Clostridium difficile.Annu. Rev. Microbiol.69445461. 10.1146/annurev-micro-091014-104115

  • 52

    TwineS. M.ReidC. W.AubryA.McMullinD. R.FultonK. M.AustinJ.et al (2009). Motility and Flagellar Glycosylation in Clostridium difficile.J. Bacteriol.19170507062. 10.1128/JB.00861-09

  • 53

    VirtanenP.GommersR.OliphantT. E.HaberlandM.ReddyT.et al (2020). SciPy 1.0: fundamental algorithms for scientific computing in Python.Nat. Methods17261272. 10.1038/s41592-019-0686-2

  • 54

    WaskomM. (2021). seaborn: statistical data visualization.J. Open Source Softw.6:3021. 10.21105/joss.03021

  • 55

    XieL.AltindalT.ChattopadhyayS.WuX.-L. (2011). Bacterial flagellum as a propeller and as a rudder for efficient chemotaxis.Proc. Natl. Acad. Sci.10822462251. 10.1073/pnas.1011953108

  • 56

    YuanJ.FahrnerK. A.TurnerL.BergH. C. (2010). Asymmetry in the clockwise and counterclockwise rotation of the bacterial flagellar motor.Proc. Natl. Acad. Sci.1071284612849. 10.1073/pnas.1007333107

Summary

Keywords

Clostridioides difficile, motility, viscoelastic medium, video microscopy, motility tracking, bacterial swimming strategy

Citation

Schwanbeck J, Oehmig I, Groß U, Zautner AE and Bohne W (2021) Clostridioides difficile Single Cell Swimming Strategy: A Novel Motility Pattern Regulated by Viscoelastic Properties of the Environment. Front. Microbiol. 12:715220. doi: 10.3389/fmicb.2021.715220

Received

26 May 2021

Accepted

29 June 2021

Published

21 July 2021

Volume

12 - 2021

Edited by

Caterina Guzmán-Verri, National University of Costa Rica, Costa Rica

Reviewed by

Jonathan David Partridge, University of Texas at Austin, United States; Mindy Engevik, Medical University of South Carolina, United States; David Courson, Old Dominion University, United States

Updates

Copyright

*Correspondence: Wolfgang Bohne,

This article was submitted to Infectious Diseases, a section of the journal Frontiers in Microbiology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics