A Multiple-Relaxation-Time Lattice-Boltzmann Model for Bacterial Chemotaxis: Effects of Initial Concentration,Diffusion, and Hydrodynamic Dispersion on Traveling Bacterial Bands

Bacterial chemotaxis can enhance the bioremediation of contaminants in aqueous and subsurface environments if the contaminant is a chemoattractant that the bacteria degrade. The process can be promoted by traveling bands of chemotactic bacteria that form due to metabolism-generated gradients in chemoattractant concentration. We developed a multiple-relaxation-time (MRT) lattice-Boltzmann method (LBM) to model chemotaxis, because LBMs are well suited to model reactive transport in the complex geometries that are typical for subsurface porous media. This MRT-LBM can attain a better numerical stability than its corresponding single-relaxation-time LBM. We performed simulations to investigate the effects of substrate diffusion, initial bacterial concentration, and hydrodynamic dispersion on the formation, shape, and propagation of bacterial bands. Band formation requires a sufficiently high initial number of bacteria and a small substrate diffusion coefficient. Uniform flow does not affect the bands while shear flow does. Bacterial bands can move both upstream and downstream when the flow velocity is small. However, the bands disappear once the velocity becomes too large due to hydrodynamic dispersion. Generally bands can only be observed if the dimensionless ratio between the chemotactic sensitivity coefficient and the effective diffusion coefficient of the bacteria exceeds a critical value, that is, when the biased movement due to chemotaxis overcomes the diffusion-like movement due to the random motility and hydrodynamic dispersion.     

A material-balance approach for modeling bacterial chemotaxis to a consumable substrate in the capillary assay

A mathematical model was developed to simulate the movement of chemotactic bacteria into and within a capillary tube containing a metabolizable chemoattractant. The model was based on a material balance that accounts for the transport of bacteria and chemoattractant as well as consumption of the substrate throughout the capillary assay system. By solving the model with a numerical method, it was possible to predict the concentration of substrate and bacteria at points within the capillary and throughout the chamber. The model was tested for its ability to simulate the results of capillary assay experiments performed with Pseudomonas putida G7 and one of its chemoattractants, naphthalene, under conditions wherein naphthalene consumption was expected to affect the flux of bacteria into the capillary. While variations in the chemotactic responses were evident among different experiments, the model could simulate the accumulation of cells in the capillary using previously determined parameters, including the chemotactic sensitivity and random motility coefficients, chi(0) and mu. In particular, model predictions were consistent with the experimental observation that the chemotactic response in the capillary is diminished under conditions wherein consumption would be expected to be significant.     

Mathematical Modeling of Chemotactic Bacterial Transport through a Two-Dimensional Heterogeneous Porous Medium

The impact of bacterial chemotaxis on in situ ground-water bioremediation remains an unanswered question. Although bacteria respond to chemical gradients in aqueous environments and under no-flow conditions, it is unclear whether they can also respond in porous media with advective flow to improve overall contaminant degradation. The effect of chemotaxis is most profound in regions with sharp chemical gradients, most notably around residual nonaqueous phase liquid (NAPL) ganglia and surrounding clay lenses or aquitards with trapped contamination. The purpose of this study is to simulate bacterial transport through a two-dimensional subsurface environment, containing one region of low permeability with trapped contaminant surrounded above and below by two regions of higher permeability. Using mathematical predictions of the effect of pore size on measured bacterial transport parameters, the authors observe a 50% decrease in both motility and chemotaxis in the finer-grained, low-permeability porous medium. The authors simulate how chemotaxis affects bacterial migration to the contaminated region under various flow and initial conditions. Results indicate that bacteria traveling through a high-permeability region with advective flow can successfully migrate toward and accumulate around a contaminant diffusing from a lower permeability region.     

Mathematical Modeling of Chemotactic Bacterial Transport through a Two-Dimensional Heterogeneous Porous Medium

The impact of bacterial chemotaxis on in situ ground-water bioremediation remains an unanswered question. Although bacteria respond to chemical gradients in aqueous environments and under no-flow conditions, it is unclear whether they can also respond in porous media with advective flow to improve overall contaminant degradation. The effect of chemotaxis is most profound in regions with sharp chemical gradients, most notably around residual nonaqueous phase liquid (NAPL) ganglia and surrounding clay lenses or aquitards with trapped contamination. The purpose of this study is to simulate bacterial transport through a two-dimensional subsurface environment, containing one region of low permeability with trapped contaminant surrounded above and below by two regions of higher permeability. Using mathematical predictions of the effect of pore size on measured bacterial transport parameters, the authors observe a 50% decrease in both motility and chemotaxis in the finer-grained, low-permeability porous medium. The authors simulate how chemotaxis affects bacterial migration to the contaminated region under various flow and initial conditions. Results indicate that bacteria traveling through a high-permeability region with advective flow can successfully migrate toward and accumulate around a contaminant diffusing from a lower permeability region.     

Chemotactic Sensing towards Ambient and Secreted Attractant Drives Collective Behaviour of E. coli

We simulate the dynamics of a suspension of bacterial swimmers, which chemotactically sense gradients in either ambient or self-secreted attractants (e.g. nutrient or aspartate respectively), or in both. Unlike previous mean field models based on a set of continuum partial differential equations, our model resolves single swimmers and therefore incorporates stochasticity and effects due to fluctuations in the bacterial density field. The algorithm we use is simple enough that we can follow the evolution of colonies of up to over a million bacteria for timescales relevant to pattern formation for E. coli growing in semisolid medium such as agar, or in confined geometries. Our results confirm previous mean field results that the patterns observed experimentally can be reproduced with a model incorporating chemoattractant secretion, chemotaxis (towards gradients in the chemoattractant field), and bacterial reproduction. They also suggest that further experiments with bacterial strains chemotactically moving up both nutrient and secreted attractant field may yield yet more dynamical patterns.     

A microscale model of bacterial and biofilm dynamics in porous media

A microscale model for the transport and coupled reaction of microbes and chemicals in an idealized two-dimensional porous media has been developed. This model includes the flow, transport, and bioreaction of nutrients, electron acceptors, and microbial cells in a saturated granular porous media. The fluid and chemicals are represented as a continuum, but the bacterial cells and solid granular particles are represented discretely. Bacterial cells can attach to the particle surfaces or be advected in the bulk fluid. The bacterial cells can also be motile and move preferentially via a run and tumble mechanism toward a chemoattractant. The bacteria consume oxygen and nutrients and alter the profiles of these chemicals. Attachment of bacterial cells to the soil matrix and growth of bacteria can change the local permeability. The coupling of mass transport and bioreaction can produce spatial gradients of nutrients and electron acceptor concentrations. We describe a numerical method for the microscale model, show results of a convergence study, and present example simulations of the model system.     

Cell balance equation for chemotactic bacteria with a biphasic tumbling frequency

Alts three-dimensional cell balance equation characterizing the chemotactic bacteria was analyzed under the presence of one-dimensional spatial chemoattractant gradients. Our work differs from that of others who have developed rather general models for chemotaxis in the use of a non-smooth anisotropic tumbling frequency function that responds biphasically to the combined temporal and spatial chemoattractant gradients. General three-dimensional expressions for the bacterial transport parameters were derived for chemotactic bacteria, followed by a perturbation analysis under the planar geometry. The bacterial random motility and chemotaxis were summarized by a motility tensor and a chemotactic velocity vector, respectively. The consequence of invoking the diffusion-approximation assumption and using intrinsic one-dimensional models with modified cellular swimming speeds was investigated by numerical simulations. Characterizing the bacterial random orientation after tumbles by a turn angle probability distribution function, we found that only the first-order angular moment of this turn angle probability distribution is important in influencing the bacterial long-term transport. Mathematics Subject Classification (2000):60G05, 60J60, 82A70     

Finite time blow-up in some models of chemotaxis

We consider a class of models of chemotactic bacterial populations, introduced by Keller-Segel. For those models, we investigate the possibility of chemotactic collapse, in other words, the possibility that in finite time the population of predators aggregates to form a delta-function. To study this phenomenon, we construct self-similar solutions, which may or may not blow-up (in finite time), depending on the relative strength of three mechanisms in competition: (i) the chemotactic attraction of bacteria towards regions of high concentration in substrate (ii) the rate of consumption of the substrate by the bacteria and (iii) (possibly) the diffusion of bacteria. The solutions we construct are radially symmetric, and therefore have no relation with the classical traveling wave solutions. Our scaling can be justified by a dimensional analysis. We give some evidence of numerical stability.     

Kinetic models for chemotaxis: Hydrodynamic limits and spatio-temporal mechanisms

We study kinetic models for chemotaxis, incorporating the ability of cells to assess temporal changes of the chemoattractant concentration as well as its spatial variations. For prescribed smooth chemoattractant density, the macroscopic limit is carried out rigorously. It leads to a drift equation with a chemotactic sensitivity depending on the time derivative of the chemoattractant density. As an application it is shown by numerical experiments that the new model can resolve the chemotactic wave paradox. For this purpose, the macroscopic equation is coupled to a simple activation-inhibition model for the chemoattractant which produces the chemoattractant waves typical for the slime mold Dictyostelium discoideum.     

Microbes to-go: slugs as source for Caenorhabditis elegans microbiota acquisition

Research on the Caenorhabditis elegans microbiota only recently started, with little known about how C. elegans acquires its microbiota. Slugs live in the same habitat as C. elegans and are known vectors for the worm. Hence, we wondered how the passage through a slug affects the C. elegans gut microbiota and whether worms can acquire bacteria from the slug. Using fluorescently labelled microbiota and 16S rRNA gene amplicon sequencing, we evaluated microbiota persistence and acquisition in C. elegans after slug passage. We compared C. elegans gut microbiomes isolated from wild-caught slugs to the microbiomes of worms after experimental slug passage to compare similarities and differences in microbiome composition. We found that microbiota persists in C. elegans while passing the slug gut and that worms simultaneously acquire additional bacteria species from the slug. Although the amplicon sequencing variant (ASV) richness of worms from the experiment did not exceed the richness of worms that naturally occur in slugs, we found a high number of shared ASVs indicating the importance of commonly associated microbiota. We demonstrate that C. elegans can take advantage of its passage through the slug by acquiring new potential microbiota without losing its native microbiota.     

Comparative toxicity of nanoparticulate CuO and ZnO to soil bacterial communities

The increasing industrial application of metal oxide Engineered Nano-Particles (ENPs) is likely to increase their environmental release to soils. While the potential of metal oxide ENPs as environmental toxicants has been shown, lack of suitable control treatments have compromised the power of many previous assessments. We evaluated the ecotoxicity of ENP (nano) forms of Zn and Cu oxides in two different soils by measuring their ability to inhibit bacterial growth. We could show a direct acute toxicity of nano-CuO acting on soil bacteria while the macroparticulate (bulk) form of CuO was not toxic. In comparison, CuSO(4) was more toxic than either oxide form. Unlike Cu, all forms of Zn were toxic to soil bacteria, and the bulk-ZnO was more toxic than the nano-ZnO. The ZnSO(4) addition was not consistently more toxic than the oxide forms. Consistently, we found a tight link between the dissolved concentration of metal in solution and the inhibition of bacterial growth. The inconsistent toxicological response between soils could be explained by different resulting concentrations of metals in soil solution. Our findings suggested that the principal mechanism of toxicity was dissolution of metal oxides and sulphates into a metal ion form known to be highly toxic to bacteria, and not a direct effect of nano-sized particles acting on bacteria. We propose that integrated efforts toward directly assessing bioavailable metal concentrations are more valuable than spending resources to reassess ecotoxicology of ENPs separately from general metal toxicity.     

生姜根系不同生态位细菌群落多样性特征、组成及结构差异

生姜作为常见的调味品和传统中药材,是我国重要的经济作物之一。作为取食部分的生姜块茎与根系直接相连,其产量、品质与根相关细菌群落密切相关。然而,关于生姜根系微环境中细菌群落的特点仍鲜有报道,土壤环境能否衍生出宿主特异性内生菌群落尚不清楚。以生姜根系不同生态位细菌群落为研究对象,采用高通量测序技术,对非根际、根际及根内细菌进行16S rRNA基因测序。结果表明,不同生态位细菌群落多样性存在显著差异,其中非根际及根际细菌群落多样性(Shannon index, Observed species, Faith′s PD)显著高于内生菌群落。同时,各生态位共现网络稳定性和复杂度表现为非根际>根际>根内细菌群落。而在组成上,细菌群落在不同生态位差异显著(R²=0.57,P=0.001)。其中变形菌门(Proteobacteria)是根内的优势门,该门类下假单胞菌属(Pseudomonas)、短波单胞菌属(Brevundimonas)、寡养单胞菌属(Stenotrophomonas)及泛菌属(Pantoea)在根内显著富集。在根际细菌中,拟杆菌门(Bacteroid...     

Analysis of the distribution and evolution of the ATP‐dependent DNA ligases of bacteria delineates a distinct phylogenetic group ‘Lig E’

Prior to the discovery of a minimal ATP‐dependent DNA ligase in Haemophilus influenzae, bacteria were thought to only possess a NAD‐dependent ligase, which was involved in sealing of Okazaki fragments. We now know that a diverse range of bacterial species possess up to six of these accessory bacterial ATP‐dependent DNA ligases (b‐ADLs), which vary in size and enzymatic domain associations. Here we compare the domain structure of different types of b‐ADLs and investigate their distribution among the bacterial domain to describe possible evolutionary trajectories that gave rise to the sequence and structural diversity of these enzymes. Previous biochemical and genetic analyses have delineated three main classes of these enzymes: Lig B, Lig C and Lig D, which appear to have descended from a common ancestor within the bacterial domain. In the present study, we delineate a fourth group of b‐ADLs, Lig E, which possesses a number of unique features at the primary and tertiary structural levels. The biochemical characteristics, domain structure and inferred extracellular location sets this group apart from the other b‐ADLs. The results presented here indicate that the Lig E type ligases were horizontally transferred into bacteria in a separate event from other b‐ADLs possibly from a bacteriophage.     

NDP52, a novel autophagy receptor for ubiquitin-decorated cytosolic bacteria

Autophagy functions as a cell-autonomous effector mechanism of innate immunity by separating bacteria from cytosolic resources and delivering them for lysosomal destruction. How cytosolic bacteria are targeted for autophagy is incompletely understood. We recently discovered that Salmonella enterica serotype Typhimurium and Streptococcus pyogenes are detected by NDP52 (nuclear dot protein 52kDa), after these bacteria enter the cytosol of human cells and become decorated with poly-ubiquitinated proteins. NDP52 binds the bacterial ubiquitin coat as well as ATG8/LC3 and delivers cytosolic bacteria into autophagosomes. In the absence of NDP52 ubiquitin-coated bacteria accumulate outside ATG8/LC3⁺ autophagosomes. Cells lacking NDP52 fail to restrict bacterial proliferation, as do cells depleted of TBK1, an IKK family kinase colocalizing with NDP52 at the bacterial surface. Our findings demonstrate the existence of a receptor for the selective autophagy of cytosolic bacteria, suggesting that cells are able to differentiate between anti-bacterial and other forms of autophagy.     

Pathogenicity of Moraxella osloensis, a bacterium associated with the nematode Phasmarhabditis hermaphrodita, to the slug Deroceras reticulatum.

Moraxella osloensis, a gram-negative bacterium, is associated with Phasmarhabditis hermaphrodita, a nematode parasite of slugs. This bacterium-feeding nematode has potential for the biological control of slugs, especially the grey garden slug, Deroceras reticulatum. Infective juveniles of P. hermaphrodita invade the shell cavity of the slug, develop into self-fertilizing hermaphrodites, and produce progeny, resulting in host death. However, the role of the associated bacterium in the pathogenicity of the nematode to the slug is unknown. We discovered that M. osloensis alone is pathogenic to D. reticulatum after injection into the shell cavity or hemocoel of the slug. The bacteria from 60-h cultures were more pathogenic than the bacteria from 40-h cultures, as indicated by the higher and more rapid mortality of the slugs injected with the former. Coinjection of penicillin and streptomycin with the 60-h bacterial culture reduced its pathogenicity to the slug. Further work suggested that the reduction and loss of pathogenicity of the aged infective juveniles of P. hermaphrodita to D. reticulatum result from the loss of M. osloensis from the aged nematodes. Also, axenic J1/J2 nematodes were nonpathogenic after injection into the shell cavity. Therefore, we conclude that the bacterium is the sole killing agent of D. reticulatum in the nematode-bacterium complex and that P. hermaphrodita acts only as a vector to transport the bacterium into the shell cavity of the slug. The identification of the toxic metabolites produced by M. osloensis is being pursued.     

Pathogenicity of Moraxella osloensis, a Bacterium Associated with the Nematode Phasmarhabditis hermaphrodita, to the Slug Deroceras reticulatum

Moraxella osloensis, a gram-negative bacterium, is associated with Phasmarhabditis hermaphrodita, a nematode parasite of slugs. This bacterium-feeding nematode has potential for the biological control of slugs, especially the grey garden slug, Deroceras reticulatum. Infective juveniles of P. hermaphrodita invade the shell cavity of the slug, develop into self-fertilizing hermaphrodites, and produce progeny, resulting in host death. However, the role of the associated bacterium in the pathogenicity of the nematode to the slug is unknown. We discovered that M. osloensis alone is pathogenic to D. reticulatum after injection into the shell cavity or hemocoel of the slug. The bacteria from 60-h cultures were more pathogenic than the bacteria from 40-h cultures, as indicated by the higher and more rapid mortality of the slugs injected with the former. Coinjection of penicillin and streptomycin with the 60-h bacterial culture reduced its pathogenicity to the slug. Further work suggested that the reduction and loss of pathogenicity of the aged infective juveniles of P. hermaphrodita to D. reticulatum result from the loss of M. osloensis from the aged nematodes. Also, axenic J1/J2 nematodes were nonpathogenic after injection into the shell cavity. Therefore, we conclude that the bacterium is the sole killing agent of D. reticulatum in the nematode-bacterium complex and that P. hermaphrodita acts only as a vector to transport the bacterium into the shell cavity of the slug. The identification of the toxic metabolites produced by M. osloensis is being pursued.     

Intracellular periodontal pathogen exploits recycling pathway to exit from infected cells

Although human gingival epithelium prevents intrusions by periodontal bacteria, Porphyromonas gingivalis, the most well‐known periodontal pathogen, is able to invade gingival epithelial cells and pass through the epithelial barrier into deeper tissues. We previously reported that intracellular P. gingivalis exits from gingival epithelial cells via a recycling pathway. However, the underlying molecular process remains unknown. In the present study, we found that the pathogen localized in early endosomes recruits VAMP2 and Rab4A. VAMP2 was found to be specifically localized in early endosomes, although its localization remained unclear in mammalian cells. A single transmembrane domain of VAMP2 was found to be necessary and sufficient for localizing in early endosomes containing P. gingivalis in gingival epithelial cells. VAMP2 forms a complex with EXOC2/Sec5 and EXOC3/Sec6, whereas Rab4A mediates dissociation of the EXOC complex followed by recruitment of RUFY1/Rabip4, Rab4A effector, and Rab14. Depletion of VAMP2 or Rab4A resulted in accumulation of bacteria in early endosomes and disturbed bacterial exit from infected cells. It is suggested that these novel dynamics allow P. gingivalis to exploit fast recycling pathways promoting further bacterial penetration of gingival tissues.     

Phytopathogenic bacteria utilize host glucose as a signal to stimulate virulence through LuxR homologues

Chemical signal-mediated biological communication is common within bacteria and between bacteria and their hosts. Many plant-associated bacteria respond to unknown plant compounds to regulate bacterial gene expression. However, the nature of the plant compounds that mediate such interkingdom communication and the underlying mechanisms remain poorly characterized. Xanthomonas campestris pv. campestris (Xcc) causes black rot disease on brassica vegetables. Xcc contains an orphan LuxR regulator (XccR) which senses a plant signal that was validated to be glucose by HPLC-MS. The glucose concentration increases in apoplast fluid after Xcc infection, which is caused by the enhanced activity of plant sugar transporters translocating sugar and cell-wall invertases releasing glucose from sucrose. XccR recruits glucose, but not fructose, sucrose, glucose 6-phosphate, and UDP-glucose, to activate pip expression. Deletion of the bacterial glucose transporter gene sglT impaired pathogen virulence and pip expression. Structural prediction showed that the N-terminal domain of XccR forms an alternative pocket neighbouring the AHL-binding pocket for glucose docking. Substitution of three residues affecting structural stability abolished the ability of XccR to bind to the luxXc box in the pip promoter. Several other XccR homologues from plant-associated bacteria can also form stable complexes with glucose, indicating that glucose may function as a common signal molecule for pathogen–plant interactions. The conservation of a glucose/XccR/pip-like system in plant-associated bacteria suggests that some phytopathogens have evolved the ability to utilize host compounds as virulence signals, indicating that LuxRs mediate an interkingdom signalling circuit.     

Two continuum models for the spreading of myxobacteria swarms

We analyze the phenomenon of spreading of a Myxococcus xanthus bacterial colony on plates coated with nutrient. The bacteria spread by gliding on the surface. In the first few hours, cell growth is irrelevant to colony spread. In this case, bacteria spread through peninsular protrusions from the edge of the initial colony. We analyze the diffusion through the narrowing reticulum of cells on the surface mathematically and derive formulae for the spreading rates. On the time scale of tens of hours, effective diffusion of the bacteria, combined with cell division and growth, causes a constant linear increase in the colony's radius. Mathematical analysis and numerical solution of reaction-diffusion equations describing the bacterial and nutrient dynamics demonstrate that, in this regime, the spreading rate is proportional to the square root of both the effective diffusion coefficient and the nutrient concentration. The model predictions agree with the data on spreading rate dependence on the type of gliding motility.