Bacterial Population Dynamics Research Articles

CanKiwi: A Mechanistic Competition Model of Kiwifruit Bacterial Canker Disease Dynamics

This paper proposes a mathematical model based on a mechanistic approach and previous research findings for the bacterial canker disease development in kiwifruit vines. This disease is a leading cause of severe damage to kiwifruit vines, particularly in humid regions, and contributes to significant economic challenges for growers in many countries. The proposed model contains three parts. The first one is the model of the kiwifruit vine describing its light interception, its carbon acquisition, and the partitioning dynamics. The carbon resource represents the chemical energy required for maintaining the necessary respiration of the living organs and their growth processes. The second part of the model is the dynamics of the pathogenic bacterial population living within the vine’s tissues and competing with them for the carbon resource required for their proliferation. The third part of the model is the carbon dynamics described by a mass conservation formula which computes the remaining amount of carbon available for competition. The model was validated by comparing simulations with experimental results obtained from growth chambers. The results show that the proposed model can simulate reasonably well the functional part of the vine in both the healthy case and the disease case without plant defense mechanisms in which the bacteria are always dominant under favorable environmental conditions. They also show that the environmental effects on the vine’s growth and the infection progress are taken into account and align with the previous studies. The model can be used to simulate the infection process, predict its outcomes, test disease management techniques, and support experimental analyses.

Pneumococcal pneumonia is driven by increased bacterial turnover due to bacteriocin-mediated intra-strain competition

Using chromosomal barcoding, we observed that >97% of the Streptococcus pneumoniae (Spn) population turns over in the lung within 2 days post-inoculation in a murine model. This marked collapse of diversity and bacterial turnover was associated with acute inflammation (severe pneumococcal pneumonia), high bacterial numbers in the lungs, bacteremia, and mortality. Intra-strain competition mediated by the blp locus, which expresses bacteriocins in a quorum-sensing-dependent manner, was required for each of these effects. Bacterial turnover from the activity of Blp-bacteriocins increased the release of the pneumococcal toxin, pneumolysin (Ply), which was sufficient to account for the lung pathology. The ability of Ply to evade complement, rather than its pore-forming activity, prevented opsonophagocytic clearance of Spn enabling its multiplication in the lung, facilitating the inflammatory response and subsequent invasion into the bloodstream. Thus, our study demonstrates how an appreciation for bacterial population dynamics during infection provides new insight into pathogenesis.

Correlations in microbial abundance data reveal host-bacteria and bacteria-bacteria interactions jointly shaping the C. elegans microbiome

Compositional structure of host-associated microbiomes is potentially affected by interactions among the microbes and between the microbes and the host. To quantify the relative importance of these contributions to the microbiome composition and variation, here we analyze absolute abundance (count) data for a minimal eight-species native microbiome in the intestine. We find that a simple neutral model only considering migration, birth, death, and competition for space among the bacteria can capture the means and variances of bacterial abundance, but not the experimental bacteria-bacteria covariances. We find that either bacteria-bacteria interactions or correlations among bacterial population dynamics parameters induced by the host can qualitatively recapitulate the observed correlations among bacterial taxa. However, both models are underdetermined and not unique, and the combination of both mechanisms likely contributes. This highlights that single time-point measurements on their own are not sufficient to distinguish between two distinct biological mechanisms. Further, we observe that different interactions are required to explain (co)variance data in microbiota associated with different host genotypes, suggesting different community dynamics associated with these host types. Finally, we find that many of these signals are obscured when data are converted to proportions from counts, consistent with a growing literature on the limitations of compositional data for inference of population dynamics. We end with a discussion of the limitations of Lotka-Volterra type assumptions for microbial community data analysis revealed by our results. Published by the American Physical Society 2024

Top-down predators shape soil bacterial community composition while bottom-up nutrients drive bacterial abundance

Although the top-down and bottom-up concept in microbial food-webs has been a primary interest in ecology, less is still known about it in soil ecosystems. Protists are the primary top-down predators of bacterial communities, altering their compositions, while the bottom-up resources are the primary factors limiting bacterial growth. Here, we hypothesized that the top-down predators modulate soil bacterial community composition, while the bottom-up nutrients control the bacterial growth and population. To precisely control nutrient levels, we used an inert soil substitute consisting of a combination of calcined clay and sand. Nutrients equivalent to the reference paddy field soil were added to microcosms as a control treatment. To investigate the effects of C, N, and P, six additional bottom-up treatments in the absence and double amounts of the nutrients were prepared. Four top-down treatments (no protist addition, Acanthamoeba castellanii, Vermamoeba vermiformis, and Heteromita globosa) were set up for each bottom-up treatment. A total of 252 microcosms under 28 treatments were incubated. Bacterial communities were analyzed using high-throughput sequencing and real-time PCR in the 1st, 3rd, and 5th weeks. The results revealed that the top-down predators significantly altered the bacterial community composition, and the bacterial population was predominantly controlled by the bottom-up nutrients. Analysis of absolute abundance data demonstrated that both top-down and bottom-up factors shaped the bacterial community structure (community composition and population). Random forest analysis classified the amplicon sequence variants associated with the treatments, showing that mostly similar families were affected by both bottom-up and top-down factors. In conclusion, the results of this study fully supported our hypothesis that top-down predators alter community composition, while bottom-up factors influence bacterial population dynamics.

Succession of Bacteria and Archaea Within the Soil Micro-Food Web Shifts Soil Respiration Dynamics.

Bacterivorous nematodes are important grazers in the soil micro-food web. Their trophic regulation shapes the composition and ecosystem services of the soil microbiome, but the underlying population dynamics of bacteria and archaea are poorly understood. We followed soil respiration and 221 dominant bacterial and archaeal 16S rRNA gene amplicon sequencing variants (ASVs) in response to top-down control by a common bacterivorous soil nematode, Acrobeloides buetschlii, bottom-up control by maize litter amendment and their combination over 32 days. Maize litter amendment significantly increased soil respiration, while A. buetschlii addition caused an earlier peak in soil respiration. Underlying bacterial and archaeal population dynamics separated into five major response types, differentiating in their temporal abundance maxima and minima. In-depth analysis of these population dynamics identified a broad imprint of A. buetschlii grazing on dominant bacterial (Acidobacteriota, Bacteroidota, Gemmatimonadota, Pseudomonadota) and archaeal (Nitrososphaerota) ASVs. Combined bottom-up control by maize litter and top-down control by A. buetschlii grazing caused a succession of soil microbiota, driven by population changes first in the Bacteroidota, then in the Pseudomonadota and finally in the Acidobacteriota and Nitrososphaerota. Our results are an essential step forward in understanding trophic modulation of soil microbiota and its feedback on soil respiration.

Diverse morphology and motility induced emergent order in bacterial collectives.

Microbial communities exhibit complex behaviors driven by species interactions and individual characteristics. In this study, we delve into the dynamics of a mixed bacterial population comprising two distinct species with different morphology and motility aspects. Employing agent-based modeling and computer simulations, we analyze the impacts of size ratios and packing fractions on dispersal patterns, aggregate formation, clustering, and spatial ordering. Notably, we find that motility and anisotropy of elongated bacteria significantly influence the distribution and spatial organization of nonmotile spherical species. Passive spherical cells display a superdiffusive behavior, particularly at larger size ratios in the ballistic regime. As the size ratio increases, clustering of passive cells is observed, accompanied by enhanced alignment and closer packing of active cells in the presence of higher passive cell area fractions. In addition, we identify the pivotal role of passive cell area fraction in influencing the response of active cells toward nematicity, with its dependence on size ratio. These findings shed light on the significance of morphology and motility in shaping the collective behavior of microbial communities, providing valuable insights into complex microbial behaviors with implications for ecology, biotechnology, and bioengineering.

Using bacterial population dynamics to count phages and their lysogens

Traditional assays for counting bacteriophages and their lysogens are labor-intensive and perturbative to the host cells. Here, we present a high-throughput infection method in a microplate reader, where the growth dynamics of the infected culture is measured using the optical density (OD). We find that the OD at which the culture lyses scales linearly with the logarithm of the initial phage concentration, providing a way of measuring phage numbers over nine orders of magnitude and down to single-phage sensitivity. Interpreting the measured dynamics using a mathematical model allows us to infer the phage growth rate, which is a function of the phage-cell encounter rate, latent period, and burst size. Adding antibiotic selection provides the ability to measure the rate of host lysogenization. Using this method, we found that when E. coli growth slows down, the lytic growth rate of lambda phages decreases, and the propensity for lysogeny increases, demonstrating how host physiology influences the viral developmental program.

Quantifying trade-offs between therapeutic efficacy and resistance dissemination for enrofloxacin dose regimens in cattle

The use of antimicrobial drugs in food-producing animals contributes to the selection pressure on pathogenic and commensal bacteria to become resistant. This study aims to evaluate the existence of trade-offs between treatment effectiveness, cost, and the dynamics of resistance in gut commensal bacteria. We developed a within-host ordinary differential equation model to track the dynamics of antimicrobial drug concentrations and bacterial populations in the site of infection (lung) and the gut. The model was parameterized to represent enrofloxacin treatment for bovine respiratory disease (BRD) caused by Pastereulla multocida in cattle. Three approved enrofloxacin dosing regimens were compared for their effects on resistance on P. multocida and commensal E. coli: 12.5 mg/kg and 7.5 mg/kg as a single dose, and 5 mg/kg as three doses. Additionally, we explored non-FDA-approved regimes. Our results indicated that both 12.5 mg/kg and 7.5 mg/kg as a single dose scenario increased the most the treatment costs and prevalence of P. multocida resistance in the lungs, while 5 mg/kg as three doses increased resistance in commensal E. coli bacteria in the gut the most out of the approved scenarios. A proposed non-FDA-approved scenario (7.5 mg/kg, two doses 24 h apart) showed low economic costs, minimal P. multocida, and moderate effects on resistant E. coli. Overall, the scenarios that decrease P. multocida, including resistant P. multocida did not coincide with those that decrease resistant E. coli the most, suggesting a trade-off between both outcomes. The sensitivity analysis suggests that bacterial populations were the most sensitive to drug conversion factors into plasma (β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\beta}$$\\end{document}), elimination of the drug from the colon (ϑ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\vartheta$$\\end{document}), fifty percent sensitive bacteria (P. multocida) killing effect (Ls50\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\ ext{L}}_{\ ext{s50}}$$\\end{document}), fifty percent of bacteria (E. coli) above ECOFF killing effect (Cr50\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\ ext{C}}_{\ ext{r50}}$$\\end{document}), and net drug transfer rate in the lung (γ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\gamma$$\\end{document}) parameters.

Antibiotic resistance gene dynamics in the commensal infant gut microbiome over the first year of life

Colonization of the infant gut is an important developmental process characterized by high carriage of antimicrobial resistance genes (ARGs) and high abundances of pathobionts. The horizontal transfer of ARGs to pathogenic bacteria represents a major public health concern. However, there is still a paucity of longitudinal studies surveilling ARGs in healthy infant guts at high temporal resolution. Furthermore, we do not yet have a clear view of how temporal variation in ARG carriage relates to the dynamics of specific bacterial populations, as well as community virulence potential. Here, we performed deep shotgun metagenomic sequencing of monthly fecal samples from a cohort of 12 infants, covering the first year of life to interrogate the infant gut microbiome for ARG content. We further relate ARG dynamics to the dynamics of taxa, virulence potential, as well as the potential for ARG mobilization. We identify a core resistome dominated by efflux systems typically associated with Enterobacteriaceae. Overall ARG carriage declined over the first year of life and showed strong contemporaneous correlation with the population dynamics of Proteobacteria. Furthermore, the majority of ARGs could be further mapped to metagenome-assembled genomes (MAGs) classified to this phylum. We were able to assign a large number of ARGs to E. coli by correlating the temporal dynamics of individual genes with species dynamics, and we show that the temporal dynamics of ARGs and virulence factors are highly correlated, suggesting close taxonomic associations between these two gene classes. Finally, we identify ARGs linked with various categories of mobile genetic elements, demonstrating preferential linkage among mobility categories and resistance to different drug classes. While individual variation in ARG carriage is substantial during infancy there is a clear reduction over the first year of life. With few exceptions, ARG abundances closely track the dynamics of pathobionts and community virulence potential. These findings emphasize the potential for development of resistant pathogens in the developing infant gut, and the importance of effective surveillance in order to detect such events.

Investigation of the growth kinetics of Streptococcus mutans bacteria on Zr-C coating surfaces deposited on 316L stainless steel substrates

The conducted research aims to describe the kinetics of Streptococcus mutans bacteria growth on the surface of Zr-C coatings with varying carbon content deposited on 316L medical-grade stainless steel substrates.The coatings were deposited using the MS PVD (Magnetron Sputtering Physical Vapour Deposition) technique at a constant temperature of 400C. Varying carbon contents in the coatings were achieved by changing the flow rate of acetylene during the process. The dynamics of bacteria growth cultivated at 37C were examined on the surface of the coatings for incubation times ranging from 0 to 72 hours. Generalised logistic functions were utilised for the mathematical description of the obtained results.The produced coatings exhibited varying carbon content, determining the structural composition ranging from a mixture of metallic Zr and ZrC through stoichiometric zirconium carbide to nanocrystalline grains of ZrC in an amorphous carbon matrix. The obtained results unequivocally demonstrate that in the case of ZrC coatings, the carbon content influences both the bacterial growth rate during the proliferation phase and the final population of bacteria.In the conducted research on bacterial population dynamics, only a single strain of Streptococcus mutans was considered, and the mathematical description was limited to reaching a pseudo-steady state by the population.The proposed model can successfully be adapted to describe the dynamics of other individual strains of bacteria dwelling on metallic surfaces as well as coatings.The proposed model can serve as a tool for studying, among other things, inflexion points of logistic curves, which are considered as the boundary between the dominance of anabolic and catabolic processes in the bacterial population.

Phacelia and Buckwheat Cover Crops’ Effects on Soil Quality in Organic Vegetable Production in a High Tunnel System

Cover crops (CCs) are regarded as beneficial to agricultural practice as an option for soil quality improvement in field production systems. The main goal of this study was to assess the impact of spring phacelia (Phacelia tanacetifolia Benth.) and buckwheat (Fagopyrum Mill.) in a crop rotation (CC–leek–parsley, 2020–2021) on the physicochemical and biological properties of the soil in an organic high tunnel system. Soil analyses involved measurements of bulk density, water capacity, soil aggregation, soil organic carbon (SOC), available soil nutrients, as well as microbial abundance and diversity. Phacelia generated more aboveground biomass (58.2 t fresh matter ha−1) than buckwheat (33.0 t ha−1), and their biomass contained 161 kg N ha−1 and 67 kg N ha−1, respectively. A large quantity of elements, such as N, Ca, P, S, B, and Cu, were found in phacelia biomass. More Mg and Na were found in buckwheat plants. The results showed that CC biomass significantly improved some of the soil physical and chemical properties, such as soil organic carbon stock and wet aggregate stability, and decreased soil bulk density. Cover crop treatments changed the dynamics of soil bacterial and fungus populations in a high tunnel system. Phacelia increased the quantity of ammonifiers and nitrifiers in the soil substantially. Further research with a long-term focus is needed to assess the impact of cover crops on soil properties, soil quality, and subsequent crop yields in high tunnel crop rotation and management systems.

A benchmarked comparison of software packages for time-lapse image processing of monolayer bacterial population dynamics.

Time-lapse microscopy provides a detailed window into the world of bacterial behavior. However, the vast amount of data produced by these techniques is difficult to analyze manually. We have analyzed four software tools designed to process such data and compared their performance, using populations of commonly studied bacterial species as our test subjects. Our findings offer a roadmap to scientists, helping them choose the right tool for their research. This comparison bridges a gap between microbiology and computational analysis, streamlining research efforts.

Bradyrhizobium and the soybean rhizosphere: Species level bacterial population dynamics in established soybean fields, rhizosphere and nodules

Bradyrhizobium and the soybean rhizosphere: Species level bacterial population dynamics in established soybean fields, rhizosphere and nodules

A phage tail-like bacteriocin suppresses competitors in metapopulations of pathogenic bacteria.

Bacteria can repurpose their own bacteriophage viruses (phage) to kill competing bacteria. Phage-derived elements are frequently strain specific in their killing activity, although there is limited evidence that this specificity drives bacterial population dynamics. Here, we identified intact phage and their derived elements in a metapopulation of wild plant-associated Pseudomonas genomes. We discovered that the most abundant viral cluster encodes a phage remnant resembling a phage tail called a tailocin, which bacteria have co-opted to kill bacterial competitors. Each pathogenic Pseudomonas strain carries one of a few distinct tailocin variants that target the variable polysaccharides in the outer membrane of co-occurring pathogenic Pseudomonas strains. Analysis of herbarium samples from the past 170 years revealed that the same tailocin and bacterial receptor variants have persisted in Pseudomonas populations. These results suggest that tailocin genetic diversity can be mined to develop targeted "tailocin cocktails" for microbial control.

Genome-resolved metagenomics reveals the effect of nutrient availability on bacterial genomic properties across 44 European freshwater lakes.

Understanding intricate microbial interactions in the environment is crucial. This is especially true for the relationships between nutrients and bacteria, as phosphorus, nitrogen and organic carbon availability are known to influence bacterial population dynamics. It has been suggested that low nutrient conditions prompt the evolutionary process of genome streamlining. This process helps conserve scarce nutrients and allows for proliferation. Genome streamlining is associated with genomic properties such as %GC content, genes encoding sigma factors, percent coding regions, gene redundancy, and functional shifts in processes like cell motility and ATP binding cassette transporters, among others. The current study aims to unveil the impact of nutrition on the genome size, %GC content, and functional properties of pelagic freshwater bacteria. We do this at finer taxonomic resolutions for many metagenomically characterized communities. Our study confirms the interplay of trophic level and genomic properties. It also highlights that different nutrient types, particularly phosphorus and nitrogen, impact these properties differently. We observed a covariation of functional traits with genome size. Larger genomes exhibit enriched pathways for motility, environmental interaction, and regulatory genes. ABC transporter genes reflect the availability of nutrients in the environment, with small genomes presumably relying more on metabolites from other organisms. We also discuss the distinct strategies different phyla adopt to adapt to oligotrophic environments. The findings contribute to our understanding of genomic adaptations within complex microbial communities.

Bacterial imaging in tumour diagnosis.

Some bacteria, such as Escherichia coli (E. coli) and Salmonella typhimurium (S. typhimurium), have an inherent ability to locate solid tumours, making them a versatile platform that can be combined with other tools to improve the tumour diagnosis and treatment. In anti-cancer therapy, bacteria function by carrying drugs directly or expressing exogenous therapeutic genes. The application of bacterial imaging in tumour diagnosis, a novel and promising research area, can indeed provide dynamic and real-time monitoring in both pre-treatment assessment and post-treatment detection. Different imaging techniques, including optical technology, acoustic imaging, magnetic resonance imaging (MRI) and nuclear medicine imaging, allow us to observe and track tumour-associated bacteria. Optical imaging, including bioluminescence and fluorescence, provides high-sensitivity and high-resolution imaging. Acoustic imaging is a real-time and non-invasive imaging technique with good penetration depth and spatial resolution. MRI provides high spatial resolution and radiation-free imaging. Nuclear medicine imaging, including positron emission tomography (PET) and single photon emission computed tomography (SPECT) can provide information on the distribution and dynamics of bacterial population. Moreover, strategies of synthetic biology modification and nanomaterial engineering modification can improve the viability and localization ability of bacteria while maintaining their autonomy and vitality, thus aiding the visualization of gut bacteria. However, there are some challenges, such as the relatively low bacterial abundance and heterogeneously distribution within the tumour, the high dimensionality of spatial datasets and the limitations of imaging labeling tools. In summary, with the continuous development of imaging technology and nanotechnology, it is expected to further make in-depth study on tumour-associated bacteria and develop new bacterial imaging methods for tumour diagnosis.

Fire-Induced Changes in Soil Properties and Bacterial Communities in Rotational Shifting Cultivation Fields in Northern Thailand.

Fire is a common practice in rotational shifting cultivation (RSC), but little is known about the dynamics of bacterial populations and the impact of fire disturbance in northern Thailand. To fill the research gap, this study aims to investigate the dynamics of soil bacterial communities and examine how the soil's physicochemical properties influence the bacterial communities in RSC fields over a period of one year following a fire. Surface soil samples (0-2 cm depth) were collected from sites with 6 (RSC-6Y) and 12 (RSC-12Y) years of fallow in Chiang Mai Province, northern Thailand at six different time points: before burning, 5 min after burning (summer), 3 months after burning (rainy season), 6 months after burning (rainy season), 9 months after burning (winter), and 12 months after burning (summer). The results revealed a reduction in the soil bacterial communities' diversity and an increase in soil nutrient levels immediately after the fire. The fire significantly influenced the abundance of Firmicutes, Proteobacteria, Acidobacteria, and Planctomycetes, but not that of Actinobacteria. At the genus level, Bacillus, Conexibacter, and Chthoniobacter showed increased abundance following the fire. During the rainy season, a recovery in the abundance of the soil bacterial communities was observed, although soil nutrient availability declined. Soil physicochemical properties such as pH, organic matter, organic carbon, electrical conductivity, cation exchange capacity, nitrate-nitrogen, available phosphorus, exchangeable potassium, total nitrogen, bulk density, sand, and silt contents significantly influenced the composition of bacterial communities. Alpha diversity indices indicated a decrease in diversity immediately after burning, followed by an increase from the early rainy season until the summer season, indicating that seasonal variation affected the composition of the soil bacterial communities. After one year of burning, an increase in bacterial richness was observed, while the diversity of the bacterial communities reverted to pre-burning levels.

Development, feeding, and sex shape the relative quantity of the nutritional obligatory symbiont Wolbachia in bed bugs.

The common bed bug, Cimex lectularius, is a hemipteran insect that feeds only on blood, and whose bites cause public health issues. Due to globalization and resistance to insecticides, this pest has undergone a significant and global resurgence in recent decades. Blood is an unbalanced diet, lacking notably sufficient B vitamins. Like all strict hematophagous arthropods, bed bugs host a nutritional symbiont supplying B vitamins. In C. lectularius, this nutritional symbiont is the intracellular bacterium Wolbachia (wCle). It is located in specific symbiotic organs, the bacteriomes, as well as in ovaries. Experimental depletion of wCle has been shown to result in longer nymphal development and lower fecundity. These phenotypes were rescued by B vitamin supplementation. Understanding the interaction between wCle and the bed bug may help to develop new pest control methods targeting the disruption of this symbiotic interaction. The objective of this work was thus to quantify accurately the density of wCle over the life cycle of the host and to describe potential associated morphological changes in the bacteriome. We also sought to determine the impact of sex, feeding status, and aging on the bacterial population dynamics. We showed that the relative quantity of wCle continuously increases during bed bug development, while the relative size of the bacteriome remains stable. We also showed that adult females harbor more wCle than males and that wCle relative quantity decreases slightly in adults with age, except in weekly-fed males. These results are discussed in the context of bed bug ecology and will help to define critical points of the symbiotic interaction during the bed bug life cycle.

Interspecies interactions in dairy biofilms drive community structure and response against cleaning and disinfection

Interspecies interactions within a biofilm community influence population dynamics and community structure, which in turn may affect the bacterial stress response to antimicrobials. This study was conducted to assess the impact of interactions between Kocuria salsicia and a three-species biofilm community (comprising Stenotrophomonas rhizophila, Bacillus licheniformis, and Microbacterium lacticum) on biofilm mass, the abundance of individual species, and their survival under a laboratory-scale cleaning and disinfection (C&D) regime. The presence of K. salsicia enhanced the cell numbers of all three species in pairwise interactions. The outcomes derived from summing up pairwise interactions did not accurately predict the bacterial population dynamics within communities of more than two species. In four-species biofilms, we observed the dominance of S. rhizophila and B. licheniformis, alongside a concurrent reduction in the cell counts of K. salsicia and M. lacticum. This pattern suggests that the underlying interactions are not purely non-transitive; instead, a more complex interplay results in the dominance of specific species. We observed that bacterial spatial organization and matrix production in different mixed-species combinations affected survival in response to C&D. Confocal microscopy analysis of spatial organization showed that S. rhizophila localized on the biofilm formed by B. licheniformis and M. lacticum, and S. rhizophila was more susceptible to C&D. Matrix production in B. licheniformis, evidenced by alterations in biofilm mass and by scanning electron microscopy, demonstrated its protective role against C&D, not only for this species itself, but also for neighbouring species. Our findings emphasise that various social interactions within a biofilm community not only affect bacterial population dynamics but also influence the biofilm community's response to C&D stress.

Higher dose makes higher lethality? A dose–response model of pulsed electric fields inactivation from multiscale coarse-graining method

As an emerging technology in liquid inactivation, one of the main challenges of pulsed electric fields (PEFs) inactivation lies in quantitatively describing and predicting its lethality to microorganisms. However, due to its cross-scaled complexity and the consequent numerous regulatory factors, there is currently still no unified framework to understand the PEF dose–response relationship and the population dynamics theoretically. In this study, a simple yet powerful model from multiscale coarse-graining method is proposed to simulate the bacterial inactivation in suspensions during PEF processing. The complex dose–response effects at the macroscale are successfully reconstructed from simple evolution rules and several coarse-graining parameters, while considering the damage and death of a single bacterium at the microscale. Our model uncovers the seemingly chaotic and even controversial dose–response relationship of PEF in literatures and systematically explores the regulatory effect of experimental parameters in a unified framework. One of the interesting findings is that PEF with shorter pulsed width enhances lethality and reduces the minimal inhibitory time at a constant energy output per pulse, owing to the phase transitions in three bacterial population dynamics (Bistability mode, Avalanche mode, and Hybrid mode). Our study provides a new insight for numerically modeling PEF lethality in liquid inactivation and could serve as a guide for dosage management in practical applications.