Growthcurver: an R package for obtaining interpretable metrics from microbial growth curves

Rapid diversification of a free-living protist is driven by adaptation to climate and habitat

Microbial eukaryotes (protists) have major functional roles in aquatic ecosystems, including the biogeochemical cycling of elements as well as occupying various roles in the food web. Despite their importance for ecosystem function, the factors that drive diversification in protists are not known. Here, we aimed to identify the factors that drive differentiation and, subsequently, speciation in a free-living protist, Synura petersenii (Chrysophyceae). We sampled five different geographic areas and utilized population genomics and quantitative trait analyses. Habitat and climate were the major drivers of diversification on the local geographical scale, while geography played a role over longer distances. In addition to conductivity and temperature, precipitation was one of the most important environmental drivers of differentiation. Our results imply that flushing episodes (floods) drive microalgal adaptation to different niches, highlighting the potential for rapid diversification in protists.

In Silico Design, In Vitro Construction, and In Vivo Application of Synthetic Small Regulatory RNAs in Bacteria

Small regulatory RNAs (sRNAs) are short non-coding RNAs in bacteria capable of post-transcriptional regulation. sRNAs have recently gained attention as tools in basic and applied sciences, for example, to fine-tune genetic circuits or biotechnological processes. Even though sRNAs often have a rather simple and modular structure, the design of functional synthetic sRNAs is not necessarily trivial. This protocol outlines how to use computational predictions and synthetic biology approaches to design, construct, and validate synthetic sRNA functionality for their application in bacteria. The computational tool, SEEDling, matches the optimal seed region with the user-selected sRNA scaffold for repression of target mRNAs. The synthetic sRNAs are assembled using Golden Gate cloning and their functionality is subsequently validated. The protocol uses the acrA mRNA as an exemplary proof-of-concept target in Escherichia coli. Since AcrA is part of a multidrug efflux pump, acrA repression can be revealed by assessing oxacillin susceptibility in a phenotypic screen. However, in case target repression does not result in a screenable phenotype, an alternative validation of synthetic sRNA functionality based on a fluorescence reporter is described.

Photo-addressable microwell devices for rapid functional screening and isolation of pathogen inhibitors from bacterial strain libraries

Discovery of new strains of bacteria that inhibit pathogen growth can facilitate improvements in biocontrol and probiotic strategies. Traditional, plate-based co-culture approaches that probe microbial interactions can impede this discovery as these methods are inherently low-throughput, labor-intensive, and qualitative. We report a second-generation, photo-addressable microwell device, developed to iteratively screen interactions between candidate biocontrol agents existing in bacterial strain libraries and pathogens under increasing pathogen pressure. Microwells (0.6 pl volume) provide unique co-culture sites between library strains and pathogens at controlled cellular ratios. During sequential screening iterations, library strains are challenged against increasing numbers of pathogens to quantitatively identify microwells containing strains inhibiting the highest numbers of pathogens. Ring-patterned 365 nm light is then used to ablate a photodegradable hydrogel membrane and sequentially release inhibitory strains from the device for recovery. Pathogen inhibition with each recovered strain is validated, followed by whole genome sequencing. To demonstrate the rapid nature of this approach, the device was used to screen a 293-membered biovar 1 agrobacterial strain library for strains inhibitory to the plant pathogen Agrobacterium tumefaciens sp. 15955. One iterative screen revealed nine new inhibitory strains. For comparison, plate-based methods did not uncover any inhibitory strains from the library (n = 30 plates). The novel pathogen-challenge screening mode developed here enables rapid selection and recovery of strains that effectively suppress pathogen growth from bacterial strain libraries, expanding this microwell technology platform toward rapid, cost-effective, and scalable screening for probiotics, biocontrol agents, and inhibitory molecules that can protect against known or emerging pathogens.

Activities of antioxidant enzymes and Hsp levels in response to elevated temperature in land snail species with varied latitudinal distribution

Land snails occupy a variety of habitats, with differing temperature and humidity regimes and exhibit a wide span of adaptations, to withstand abiotic condition changes. The present work's aim was to examine the correlation of habitat's thermal adversity in different Mediterranean type habitats with the land snail's antioxidant and heat shock responses. For this purpose, snails of different species from populations along a north-south axis from the islands and mainland of Greece were exposed to elevated temperature and antioxidant enzyme activities, and Hsp70 and Hsp90 levels were determined in their tissues. The ATP, ADP, and AMP levels and the adenylate energy charge (AEC) were also determined. The comparison of protein levels and enzymatic activities across time intervals revealed significant differences for all factors examined. While the gradation pattern over time for a given factor was similar in all populations the absolute values over time differed. Catalase activity and the Hsp90 protein levels had the higher contribution in separating the different species and populations, followed by the activity of glutathione reductase and Hsp70 protein levels which contributed to a lesser degree. In general, populations from the southern part of their distribution in Greece tend to display a faster increase than northern populations in induction levels of all factors examined. Our data seem to be in line with the concept of preparation for oxidative stress (POS) while the changes in the AEC indicate an early preparation to cover the energy demand for the induction and synthesis of antioxidant enzymes and Hsps.

A conserved transcription factor controls gluconeogenesis via distinct targets in hypersaline-adapted archaea with diverse metabolic capabilities

Timely regulation of carbon metabolic pathways is essential for cellular processes and to prevent futile cycling of intracellular metabolites. In Halobacterium salinarum, a hypersaline adapted archaeon, a sugar-sensing TrmB family protein controls gluconeogenesis and other biosynthetic pathways. Notably, Hbt. salinarum does not utilize carbohydrates for energy, uncommon among Haloarchaea. We characterized a TrmB-family transcriptional regulator in a saccharolytic generalist, Haloarcula hispanica, to investigate whether the targets and function of TrmB, or its regulon, is conserved in related species with distinct metabolic capabilities. In Har. hispanica, TrmB binds to 15 sites in the genome and induces the expression of genes primarily involved in gluconeogenesis and tryptophan biosynthesis. An important regulatory control point in Hbt. salinarum, activation of ppsA and repression of pykA, is absent in Har. hispanica. Contrary to its role in Hbt. salinarum and saccharolytic hyperthermophiles, TrmB does not act as a global regulator: it does not directly repress the expression of glycolytic enzymes, peripheral pathways such as cofactor biosynthesis, or catabolism of other carbon sources in Har. hispanica. Cumulatively, these findings suggest rewiring of the TrmB regulon alongside metabolic network evolution in Haloarchaea.

Mistranslating the genetic code with leucine in yeast and mammalian cells

Translation fidelity relies on accurate aminoacylation of transfer RNAs (tRNAs) by aminoacyl-tRNA synthetases (AARSs). AARSs specific for alanine (Ala), leucine (Leu), serine, and pyrrolysine do not recognize the anticodon bases. Single nucleotide anticodon variants in their cognate tRNAs can lead to mistranslation. Human genomes include both rare and more common mistranslating tRNA variants. We investigated three rare human tRNALeu variants that mis-incorporate Leu at phenylalanine or tryptophan codons. Expression of each tRNALeu anticodon variant in neuroblastoma cells caused defects in fluorescent protein production without significantly increased cytotoxicity under normal conditions or in the context of proteasome inhibition. Using tRNA sequencing and mass spectrometry we confirmed that each tRNALeu variant was expressed and generated mistranslation with Leu. To probe the flexibility of the entire genetic code towards Leu mis-incorporation, we created 64 yeast strains to express all possible tRNALeu anticodon variants in a doxycycline-inducible system. While some variants showed mild or no growth defects, many anticodon variants, enriched with G/C at positions 35 and 36, including those replacing Leu for proline, arginine, alanine, or glycine, caused dramatic reductions in growth. Differential phenotypic defects were observed for tRNALeu mutants with synonymous anticodons and for different tRNALeu isoacceptors with the same anticodon. A comparison to tRNAAla anticodon variants demonstrates that Ala mis-incorporation is more tolerable than Leu at nearly every codon. The data show that the nature of the amino acid substitution, the tRNA gene, and the anticodon are each important factors that influence the ability of cells to tolerate mistranslating tRNAs.

Acinetobacter baumannii Survival under Infection-Associated Stresses Depends on the Expression of Resistance-Nodulation-Division and Major Facilitator Superfamily Efflux Pumps

Multidrug efflux transporters are major contributors to the antibiotic resistance of Acinetobacter baumannii in clinical settings. Previous studies showed that these transporters are tightly integrated into the physiology of A. baumannii and have diverse functions. However, for many of the efflux pumps, such functions remain poorly defined. In this study, we characterized two putative drug efflux pumps, AmfAB and AmfCD (Acinetobacter Major Facilitator), that are homologous to EmrAB-like transporters from Escherichia coli and other Gram-negative bacteria. These pumps comprise the Major Facilitator Superfamily (MFS) transporters AmfB and AmfD and the periplasmic membrane fusion proteins AmfA and AmfC, respectively. We inactivated and overproduced these pumps in the wild-type ATCC 17978 strain and its derivative strains lacking the major efflux pumps from the Resistance-Nodulation-Division (RND) superfamily and characterized antibiotic susceptibilities and growth of the strains under stresses typical during human infections. We found that neither AmfAB nor AmfCD contribute to the antibiotic non-susceptibility phenotypes of A. baumannii. The two pumps, however, are critical for the adaptation and growth of the bacterium under acidic stress, whereas AmfCD also contributes to growth under conditions of low iron, high temperature, and in the presence of bile salts. These functions are dependent on the presence of the RND pumps, the inactivation of which further diminishes A. baumannii survival and growth. Our results suggest that MFS transporters contribute to stress survival by affecting the permeability properties of the A. baumannii cell envelope.

Induction of Aspergillus fumigatus zinc cluster transcription factor OdrA/Mdu2 provides combined cellular responses for oxidative stress protection and multiple antifungal drug resistance

IMPORTANCE

An overexpression screen of 228 zinc cluster transcription factor encoding genes of A. fumigatus revealed 11 genes conferring increased tolerance to antifungal drugs. Out of these, four oxidative stress and drug tolerance transcription factor encoding odr genes increased tolerance to oxidative stress and antifungal drugs when overexpressed. This supports a correlation between oxidative stress response and antifungal drug tolerance in A. fumigatus. OdrA/Mdu2 is required for the cross-tolerance between azoles, polyenes, and oxidative stress and activates genes for detoxification. Under oxidative stress conditions or when overexpressed, OdrA/Mdu2 accumulates in the nucleus and activates detoxifying genes by direct binding at their promoters, as we describe with the mdr1 gene encoding an itraconazole specific efflux pump. Finally, this work gives new insights about drug and stress resistance in the opportunistic pathogenic fungus A. fumigatus.

The genetics of aerotolerant growth in an alphaproteobacterium with a naturally reduced genome

IMPORTANCE

The inherent complexity of biological systems is a major barrier to our understanding of cellular physiology. Bacteria with markedly fewer genes than their close relatives, or reduced genome bacteria, are promising biological models with less complexity. Reduced genome bacteria can also have superior properties for industrial use, provided the reduction does not overly restrict strain robustness. Naturally reduced genome bacteria, such as the alphaproteobacterium Zymomonas mobilis, have fewer genes but remain environmentally robust. In this study, we show that Z. mobilis is a simplified genetic model for Alphaproteobacteria, a class with important impacts on the environment, human health, and industry. We also identify genes that are only required in the absence of atmospheric oxygen, uncovering players that maintain and utilize the cellular energy state. Our findings have broad implications for the genetics of Alphaproteobacteria and industrial use of Z. mobilis to create biofuels and bioproducts.

The rise of the sourdough: Genome-scale metabolic modeling-based approach to design sourdough starter communities with tailored-made properties

Sourdough starters harbor microbial consortia that benefit the final product's aroma and volume. The complex nature of these spontaneously developed communities raises challenges in predicting the fermentation phenotypes. Herein, we demonstrated for the first time in this field the potential of genome-scale metabolic modeling (GEMs) in the study of sourdough microbial communities. Broad in-silico modeling of microbial growth was applied on communities composed of yeast (Saccharomyces cerevisiae) and different Lactic Acid Bacteria (LAB) species, which mainly predominate in sourdough starters. Simulations of model-represented communities associated specific bacterial compositions with sourdough phenotypes. Based on ranking the phenotypic performances of different combinations, Pediococcus spp. - Lb. sakei group members were predicted to have an optimal effect considering the increase in S. cerevisiae growth abilities and overall CO2 secretion rates. Flux Balance Analysis (FBA) revealed mutual relationships between the Pediococcus spp. - Lb. sakei group members and S. cerevisiae through bidirectional nutrient dependencies, and further underlined that these bacteria compete with the yeast over nutrients to a lesser extent than the rest LAB species. Volatile compounds (VOCs) production was further modeled, identifying species-specific and community-related VOCs production profiles. The in-silico models' predictions were validated by experimentally building synthetic sourdough communities and assessing the fermentation phenotypes. The Pediococcus spp. - Lb. sakei group was indeed associated with increased yeast cell counts and fermentation rates, demonstrating a 25 % increase in the average leavening rates during the first 10 fermentation hours compared to communities with a lower representation of these group members. Overall, these results provide a possible novel strategy towards the de-novo design of sourdough starter communities with tailored-made characterizations, including a shortened leavening period.

Jumbo phages are active against extensively drug-resistant eyedrop-associated Pseudomonas aeruginosa infections

Antibiotic-resistant bacteria present an emerging challenge to human health. Their prevalence has been increasing across the globe due in part to the liberal use of antibiotics that has pressured them to develop resistance. Those bacteria that acquire mobile genetic elements are especially concerning because those plasmids may be shared readily with other microbes that can then also become antibiotic resistant. Serious infections have recently been related to the contamination of preservative-free eyedrops with extensively drug-resistant (XDR) isolates of Pseudomonas aeruginosa, already resulting in three deaths. These drug-resistant isolates cannot be managed with most conventional antibiotics. We sought to identify alternatives to conventional antibiotics for the lysis of these XDR isolates and identified multiple bacteriophages (viruses that attack bacteria) that killed them efficiently. We found both jumbo phages (>200 kb in genome size) and non-jumbo phages that were active against these isolates, the former killing more efficiently. Jumbo phages effectively killed the three separate XDR P. aeruginosa isolates both on solid and liquid medium. Given the ongoing nature of the XDR P. aeruginosa eyedrop outbreak, the identification of phages active against them provides physicians with several novel potential alternatives for treatment.

Alleviation of DSS-induced colitis via bovine colostrum-derived extracellular vesicles with microRNA let-7a-5p is mediated by regulating Akkermansia and -hydroxybutyrate in gut environments

IMPORTANCE

Even though studying on the possible involvement of extracellular vesicles (EVs) in host-microbe interactions, how these relationships mediate host physiology has not clarified yet. Our current findings provide insights into the encouraging benefits of dietary source-derived EVs and microRNAs (miRNAs) on organic acid production and ultimately stimulating gut microbiome for human health, suggesting that supplementation of dietary colostrum EVs and miRNAs is a novel preventive strategy for the treatment of inflammatory bowel disease.

Bifidobacterium infantis utilizes N-acetylglucosamine-containing human milk oligosaccharides as a nitrogen source

Bifidobacterium longum subsp. infantis (B. infantis) utilizes oligosaccharides secreted in human milk as a carbohydrate source. These human milk oligosaccharides (HMOs) integrate the nitrogenous residue N-acetylglucosamine (NAG), although HMO nitrogen utilization has not been described to date. Herein, we characterize the B. infantis nitrogen utilization phenotype on two NAG-containing HMO species, LNT and LNnT. This was characterized through in vitro growth kinetics, incorporation of isotopically labeled NAG nitrogen into the proteome, as well as modulation of intracellular 2-oxoglutarate levels while utilizing HMO nitrogen. Further support is provided by comparative transcriptomics and proteomics that identified global regulatory networks deployed during HMO nitrogen utilization. The aggregate data demonstrate that B. infantis strains utilize HMO nitrogen with the potential to significantly impact fundamental and clinical studies, as well as enable applications.

Parsley: a web app for parsing data from plate readers

SUMMARY

As demand for the automation of biological assays has increased over recent years, the range of measurement types implemented by multiwell plate readers has broadened and the list of published software packages that caters to their analysis has grown. However, most plate readers export data in esoteric formats with little or no metadata, while most analytical software packages are built to work with tidy data accompanied by associated metadata. 'Parser' functions are therefore required to prepare raw data for analysis. Such functions are instrument- and data type-specific, and to date, no generic tool exists that can parse data from multiple data types or multiple plate readers, despite the potential for such a tool to speed up access to analysed data and remove an important barrier for less confident coders. We have developed the interactive web application, Parsley, to bridge this gap. Unlike conventional programmatic parser functions, Parsley makes few assumptions about exported data, instead employing user inputs to identify and extract data from data files. In doing so, it is designed to enable any user to parse plate reader data and can handle a wide variety of instruments (10+) and data types (53+). Parsley is freely available via a web interface, enabling access to its unique plate reader data parsing functionality, without the need to install software or write code.

AVAILABILITY AND IMPLEMENTATION

The Parsley web application can be accessed at: https://gbstan.shinyapps.io/parsleyapp/. The source code is available at: https://github.com/ec363/parsleyapp and is archived on Zenodo: https://zenodo.org/records/10011752.

Strain-level diversity of symbiont communities between individuals and populations of a bioluminescent fish

The bioluminescent symbiosis involving the urchin cardinalfish, Siphamia tubifer, and Photobacterium mandapamensis, a luminous member of the Vibrionaceae, is highly specific compared to other bioluminescent fish-bacteria associations. Despite this high degree of specificity, patterns of genetic diversity have been observed for the symbionts from hosts sampled over relatively small spatial scales. We characterized and compared sub-species, strain-level symbiont diversity within and between S. tubifer hosts sampled from the Philippines and Japan using PCR fingerprinting. We then carried out whole genome sequencing of the unique symbiont genotypes identified to characterize the genetic diversity of the symbiont community and the symbiont pangenome. We determined that an individual light organ contains six symbiont genotypes on average, but varied between 1-13. Additionally, we found that there were few genotypes shared between hosts from the same location. A phylogenetic analysis of the unique symbiont strains indicated location-specific clades, suggesting some genetic differentiation in the symbionts between host populations. We also identified symbiont genes that were variable between strains, including luxF, a member of the lux operon, which is responsible for light production. We quantified the light emission and growth rate of two strains missing luxF along with the other strains isolated from the same light organs and determined that strains lacking luxF were dimmer but grew faster than most of the other strains, suggesting a potential metabolic trade-off. This study highlights the importance of strain-level diversity in microbial associations and provides new insight into the underlying genetic architecture of intraspecific symbiont communities within a host.

Copper reduces the virulence of bacterial communities at environmentally relevant concentrations

Increasing environmental concentrations of metals as a result of anthropogenic pollution are significantly changing many microbial communities. While there is evidence metal pollution can result in increased antibiotic resistance, the effects of metal pollution on the virulence of bacterial communities remains largely undetermined. Here, we experimentally test whether metal stress alters the virulence of bacterial communities. We do this by incubating three wastewater influent communities under different environmentally relevant copper concentrations for three days. We then quantify the virulence of the community phenotypically using the Galleria mellonella infection model, and test if differences are due to changes in the rate of biomass accumulation (productivity), copper resistance, or community composition (quantified using 16S amplicon sequencing). The virulence of the communities was found to be reduced by the highest copper concentration, but not to be affected by the lower concentration. As well as reduced virulence, communities exposed to the highest copper concentration were less diverse and had lower productivity. This work highlights that metal pollution may decrease virulence in bacterial communities, but at a cost to diversity and productivity.

Thermal adaptation of soil microbial growth traits in response to chronic warming

IMPORTANCE

Soils are the largest terrestrial carbon sink and the foundation of our food, fiber, and fuel systems. Healthy soils are carbon sinks, storing more carbon than they release. This reduces the amount of carbon dioxide released into the atmosphere and buffers against climate change. Soil microbes drive biogeochemical cycling and contribute to soil health through organic matter breakdown, plant growth promotion, and nutrient distribution. In this study, we determined how soil microbial growth traits respond to long-term soil warming. We found that bacterial isolates from warmed plots showed evidence of adaptation of optimum growth temperature. This suggests that increased microbial biomass and growth in a warming world could result in greater carbon storage. As temperatures increase, greater microbial activity may help reduce the soil carbon feedback loop. Our results provide insight on how atmospheric carbon cycling and soil health may respond in a warming world.

Unlocking the secrets of peptide transport in wine yeast: insights into oligopeptide transporter functions and nitrogen source preferences

IMPORTANCE

Limited nitrogen supply can prevent the completion of alcoholic fermentation. Supplementation through peptides as an alternative, natural source of nitrogen for yeast offers an interesting solution for this issue. In this work, the S. cerevisiae peptide transporters of the Opt and Fot families were studied. We demonstrated that Fot and Opt2 have a broader peptide length preference than previously reported, enabling yeasts to acquire sufficient nitrogen from peptides without requiring additional ammonia or amino acids to complete fermentation. On the contrary, Opt1 was unable to consume any peptide in the given conditions, whereas it has been described elsewhere as the main peptide transporter for peptides longer than three amino acid residues in experiments in laboratory conditions. This controversy signifies the need in applied sciences for approaching experimental conditions to those prevalent in the industry for its more accurate characterization. Altogether, this work provides further evidence of the importance of peptides as a nitrogen source for yeast and their consequent positive impact on fermentation kinetics.

Acquired Triazole Resistance Alters Pathogenicity-Associated Features in Candida auris in an Isolate-Dependent Manner

Fluconazole resistance is commonly encountered in Candida auris, and the yeast frequently displays resistance to other standard drugs, which severely limits the number of effective therapeutic agents against this emerging pathogen. In this study, we aimed to investigate the effect of acquired azole resistance on the viability, stress response, and virulence of this species. Fluconazole-, posaconazole-, and voriconazole- resistant strains were generated from two susceptible C. auris clinical isolates (0381, 0387) and compared under various conditions. Several evolved strains became pan-azole-resistant, as well as echinocandin-cross-resistant. While being pan-azole-resistant, the 0381-derived posaconazole-evolved strain colonized brain tissue more efficiently than any other strain, suggesting that fitness cost is not necessarily a consequence of resistance development in C. auris. All 0387-derived evolved strains carried a loss of function mutation (R160S) in BCY1, an inhibitor of the PKA pathway. Sequencing data also revealed that posaconazole treatment can result in ERG3 mutation in C. auris. Despite using the same mechanisms to generate the evolved strains, both genotype and phenotype analysis highlighted that the development of resistance was unique for each strain. Our data suggest that C. auris triazole resistance development is a highly complex process, initiated by several pleiotropic factors.

Design, construction, and functional characterization of a tRNA neochromosome in yeast

Here, we report the design, construction, and characterization of a tRNA neochromosome, a designer chromosome that functions as an additional, de novo counterpart to the native complement of Saccharomyces cerevisiae. Intending to address one of the central design principles of the Sc2.0 project, the ∼190-kb tRNA neochromosome houses all 275 relocated nuclear tRNA genes. To maximize stability, the design incorporates orthogonal genetic elements from non-S. cerevisiae yeast species. Furthermore, the presence of 283 rox recombination sites enables an orthogonal tRNA SCRaMbLE system. Following construction in yeast, we obtained evidence of a potent selective force, manifesting as a spontaneous doubling in cell ploidy. Furthermore, tRNA sequencing, transcriptomics, proteomics, nucleosome mapping, replication profiling, FISH, and Hi-C were undertaken to investigate questions of tRNA neochromosome behavior and function. Its construction demonstrates the remarkable tractability of the yeast model and opens up opportunities to directly test hypotheses surrounding these essential non-coding RNAs.

Integrated genomic and functional analyses of human skin-associated Staphylococcus reveal extensive inter- and intra-species diversity

Human skin is stably colonized by a distinct microbiota that functions together with epidermal cells to maintain a protective physical barrier. Staphylococcus, a prominent genus of the skin microbiota, participates in colonization resistance, tissue repair, and host immune regulation in strain-specific manners. To unlock the potential of engineering skin microbial communities, we aim to characterize the diversity of this genus within the context of the skin environment. We reanalyzed an extant 16S rRNA amplicon dataset obtained from distinct body sites of healthy volunteers, providing a detailed biogeographic depiction of staphylococcal species that colonize our skin. S. epidermidis, S. capitis, and S. hominis were the most abundant staphylococcal species present in all volunteers and were detected at all body sites. Pan-genome analysis of isolates from these three species revealed that the genus-core was dominated by central metabolism genes. Species-restricted-core genes encoded known host colonization functions. The majority (~68%) of genes were detected only in a fraction of isolate genomes, underscoring the immense strain-specific gene diversity. Conspecific genomes grouped into phylogenetic clades, exhibiting body site preference. Each clade was enriched for distinct gene sets that are potentially involved in site tropism. Finally, we conducted gene expression studies of select isolates showing variable growth phenotypes in skin-like medium. In vitro expression revealed extensive intra- and inter-species gene expression variation, substantially expanding the functional diversification within each species. Our study provides an important resource for future ecological and translational studies to examine the role of shared and strain-specific staphylococcal genes within the skin environment.

The secondary metabolite hydrogen cyanide protects Pseudomonas aeruginosa against sodium hypochlorite-induced oxidative stress

The high pathogenicity of Pseudomonas aeruginosa is attributed to the production of many virulence factors and its resistance to several antimicrobials. Among them, sodium hypochlorite (NaOCl) is a widely used disinfectant due to its strong antimicrobial effect. However, bacteria develop many mechanisms to survive the damage caused by this agent. Therefore, this study aimed to identify novel mechanisms employed by P. aeruginosa to resist oxidative stress induced by the strong oxidizing agent NaOCl. We analyzed the growth of the P. aeruginosa mutants ΔkatA, ΔkatE, ΔahpC, ΔahpF, ΔmsrA at 1 μg/mL NaOCl, and showed that these known H2O2 resistance mechanisms are also important for the survival of P. aeruginosa under NaOCl stress. We then conducted a screening of the P. aeruginosa PA14 transposon insertion mutant library and identified 48 mutants with increased susceptibility toward NaOCl. Among them were 10 mutants with a disrupted nrdJa, bvlR, hcnA, orn, sucC, cysZ, nuoJ, PA4166, opmQ, or thiC gene, which also exhibited a significant growth defect in the presence of NaOCl. We focussed our follow-up experiments (i.e., growth analyzes and kill-kinetics) on mutants with defect in the synthesis of the secondary metabolite hydrogen cyanide (HCN). We showed that HCN produced by P. aeruginosa contributes to its resistance toward NaOCl as it acts as a scavenger molecule, quenching the toxic effects of NaOCl.

Vesicle-enriched secretomes alter bacterial competitive abilities and are drivers of evolution in microbial communities

Microbial membrane vesicles can carry compounds that inhibit bacterial growth, but how they impact the fitness of the vesicle-producing bacterial species and influence community dynamics remain unexplored questions. To address these questions, we examined the effect of vesicle-enriched secretomes (VESs) in different single-species and multi-species systems. Effects of VESs on single-species growth dynamics were determined for nine bacterial species belonging to four genera (Escherichia, Salmonella, Pseudomonas and Bacillus) in nutrient-rich and poor growth media. Results showed both species-specific and nutrient-dependent effects of the VESs on bacterial growth. The strongest antagonistic effects were observed for VES isolated from the natural isolates of E. coli, while those isolated from P. aeruginosa PA14 affected the highest number of species. We further demonstrated that these VESs altered the competitive abilities of the species involved in two-species (S. Typhimurium LT2 and S. arizonae) and three-species systems (E. coli, S. Typhimurium LT2 and B. subtilis). Finally, using experimental evolution we showed that different bacterial species could rapidly acquire mutations that abrogated the antagonistic effects of VESs. This study demonstrates how VESs can contribute in shaping microbial communities, both by increasing the competitive ability of a given bacterial species and as a driver of genetic adaptation.

Characterization of potential probiotic starter cultures of lactic acid bacteria isolated from Ethiopian fermented cereal beverages, Naaqe, and Cheka

AIMS

To test the in vitro probiotic potential and starter culture capacity of lactic acid bacteria (LAB) isolated from Naaqe and Cheka, cereal-based Ethiopian traditional fermented beverages.

METHODS AND RESULTS

A total of 44 strains were isolated from spontaneously fermented Ethiopian cereal-based beverages, Naaqe and Cheka with 24 putatively identified as LAB and 14 identified up to the species level. The species Limosilactobacillus fermentum (6/12; 50%) and Weissella confusa (5/12, 41.67%) were the predominant species identified from Naaqe, while the two Cheka isolates were L. fermentum and Pediococcus pentosaceus. Six LAB strains inhibited eight of the nine gastrointestinal indicator key pathogens in Ethiopia, including Escherichia coli, Salmonella enterica subsp. enterica var. Typhimurium, Staphylococcus aureus, Shigella flexneri, and Listeria monocytogenes. Three of the LAB isolates exhibited strain-specific immunostimulation in human monocytes. Based on these probiotic properties and growth, six strains were selected for in situ evaluation in a mock fermentation of Naaqe and Cheka. During primary fermentations, L. fermentum 73B, P. pentosaceus 74D, L. fermentum 44B, W. confusa 44D, L. fermentum 82C, and Weissella cibaria 83E and their combinations demonstrated higher pH-lowering properties and colony-forming unit counts compared to the control spontaneous fermentation. The same pattern was also observed in the secondary mock fermentation by the Naaqe LAB isolates.

CONCLUSIONS

In this study, we selected six LAB strains with antipathogenic, immunostimulatory, and starter culture potentials that can be used as autochthonous probiotic starters for Naaqe and Cheka fermentations once their health benefit is ascertained in a clinical trial as a next step.

The toxicological spectrum of the Bacillus cereus toxin cereulide points towards niche-specific specialisation

Most microbes share their environmental niches with very different forms of life thereby engaging in specialised relationships to enable their persistence. The bacterium Bacillus cereus occurs ubiquitously in the environment with certain strain backgrounds causing foodborne and opportunistic infections in humans. The emetic lineage of B. cereus is capable of producing the toxin cereulide, which evokes emetic illnesses. Although food products favouring the accumulation of cereulide are known, the ecological role of cereulide and the environmental niche of emetic B. cereus remain elusive. To better understand the ecology of cereulide-producing B. cereus, we systematically assayed the toxicological spectrum of cereulide on a variety of organisms belonging to different kingdoms. As cereulide is a potassium ionophore, we further tested the effect of environmental potassium levels on the action of cereulide. We found that adverse effects of cereulide exposure are species-specific, which can be exacerbated with increased environmental potassium. Additionally, we demonstrate that cereulide is produced within an insect cadaver indicating its potential ecological function for a saprophytic lifestyle. Collectively, distinct cereulide susceptibilities of other organisms may reflect its role in enabling competitive niche specialization of emetic B. cereus.

Metabolic similarity and the predictability of microbial community assembly

When microbial communities form, their composition is shaped by selective pressures imposed by the environment. Can we predict which communities will assemble under different environmental conditions? Here, we hypothesize that quantitative similarities in metabolic traits across metabolically similar environments lead to predictable similarities in community composition. To that end, we measured the growth rate and by-product profile of a library of proteobacterial strains in a large number of single nutrient environments. We found that growth rates and secretion profiles were positively correlated across environments when the supplied substrate was metabolically similar. By analyzing hundreds of in-vitro communities experimentally assembled in an array of different synthetic environments, we then show that metabolically similar substrates select for taxonomically similar communities. These findings lead us to propose and then validate a comparative approach for quantitatively predicting the effects of novel substrates on the composition of complex microbial consortia.

Fitness Costs of Antibiotic Resistance Impede the Evolution of Resistance to Other Antibiotics

Antibiotic resistance is a major threat to global health, claiming the lives of millions every year. With a nearly dry antibiotic development pipeline, novel strategies are urgently needed to combat resistant pathogens. One emerging strategy is the use of sequential antibiotic therapy, postulated to reduce the rate at which antibiotic resistance evolves. Here, we use the soft agar gradient evolution (SAGE) system to carry out high-throughput in vitro bacterial evolution against antibiotic pressure. We find that evolution of resistance to the antibiotic chloramphenicol (CHL) severely affects bacterial fitness, slowing the rate at which resistance to the antibiotics nitrofurantoin and streptomycin emerges. In vitro acquisition of compensatory mutations in the CHL-resistant cells markedly improves fitness and nitrofurantoin adaptation rates but fails to restore rates to wild-type levels against streptomycin. Genome sequencing reveals distinct evolutionary paths to resistance in fitness-impaired populations, suggesting resistance trade-offs in favor of mitigation of fitness costs. We show that the speed of bacterial fronts in SAGE plates is a reliable indicator of adaptation rates and evolutionary trajectories to resistance. Identification of antibiotics whose mutational resistance mechanisms confer stable impairments may help clinicians prescribe sequential antibiotic therapies that are less prone to resistance evolution.

The PitA protein contributes to colistin susceptibility in Pseudomonas aeruginosa

Pseudomonas aeruginosa is an opportunistic pathogen that causes a wide range of problematic infections in individuals with predisposing conditions. Infections can be treated with colistin but some isolates are resistant to this antibiotic. To better understand the genetic basis of resistance, we experimentally evolved 19 independent resistant mutants from the susceptible laboratory strain PAO1. Whole genome sequencing identified mutations in multiple genes including phoQ and pmrB that have previously been associated with resistance, pitA that encodes a phosphate transporter, and carB and eno that encode enzymes of metabolism. Individual mutations were engineered into the genome of strain PAO1. Mutations in pitA, pmrB and phoQ increased the minimum inhibitory concentration (MIC) for colistin 8-fold, making the bacteria resistant. Engineered pitA/phoQ and pitA/pmrB double mutants had higher MICs than single mutants, demonstrating additive effects on colistin susceptibility. Single carB and eno mutations did not increase the MIC suggesting that their effect is dependent on the presence of other mutations. Many of the resistant mutants had increased susceptibility to β-lactams and lower growth rates than the parental strain demonstrating that colistin resistance can impose a fitness cost. Two hundred and fourteen P. aeruginosa isolates from a range of sources were tested and 18 (7.8%) were colistin resistant. Sequence variants in genes identified by experimental evolution were present in the 18 resistant isolates and may contribute to resistance. Overall our results identify pitA mutations as novel contributors to colistin resistance and demonstrate that resistance can reduce fitness of the bacteria.

Cultivating efficiency: High-throughput growth analysis of anaerobic bacteria in compact microplate readers

Anaerobic microbes play crucial roles in environmental processes, industry, and human health. Traditional methods for monitoring the growth of anaerobes, including plate counts or subsampling broth cultures for optical density measurements, are time and resource intensive. The advent of microplate readers revolutionized bacterial growth studies by enabling high-throughput and real-time monitoring of microbial growth kinetics but their use in anaerobic microbiology has remained limited. Here, we present a workflow for using small-footprint microplate readers and the Growthcurver R package to analyze the kinetic growth metrics of anaerobic bacteria. We benchmarked the small-footprint Cerillo Stratus microplate reader against a BioTek Synergy HTX microplate reader in aerobic conditions using Escherichia coli DSM 28618 cultures. The growth rates and carrying capacities obtained from the two readers were statistically indistinguishable. However, the area under the logistic curve was significantly higher in cultures monitored by the Stratus reader. We used the Stratus to quantify the growth responses of anaerobically grown E. coli and Clostridium bolteae DSM 29485 to different doses of the toxin sodium arsenite. The growth of E. coli and C. bolteae was sensitive to arsenite doses of 1.3 μM and 0.4 μM, respectively. Complete inhibition of growth was achieved at 38 μM arsenite for C. bolteae, and 338 μM in E. coli. These results show that the Stratus performs similarly to a leading brand of microplate reader and can be reliably used in anaerobic conditions. We discuss the advantages of the small format microplate readers and our experiences with the Stratus.

Characterization of Trichoderma species from forest ecosystems by high-throughput phenotypic microarray

The use of beneficial organisms for the biocontrol of soil-borne pathogens in forestry is still poor explored. In this work, the nutritional demands of 10 previously selected isolates of Trichoderma for the biocontrol of forest soil-borne pathogens have been tested by Phenotype Microarray technology, to investigate about their C-source utilization and exploring the possibility to obtain a microbial consortia (SynCom), an innovative strategy for the biocontrol of plant disease. All Trichoderma isolates tested in this study showed a high spore germination percentage within 3 d and evidenced nutritional preference regardless of the species they belong to, and unrelated to their soil of origin. Results of growth curve analysis and MANOVA test revealed that all isolates assimilate a broad range of substrates, generally preferring complex compounds such as monosaccharides related compounds, nitrogen compounds, carboxylic acids and esters. No evidence of competition for nutritional resources have been observed among isolates of this study. As a result, a combination of different isolates could be proposed to obtain a SynCom useful for the practice of phytopathogen biocontrol in forestry. The addition of i-erythritol, adenosine and turanose to a growth substrate could be suggested as stimulating compounds for the growth of the selected Trichoderma isolates.

Evolutionary innovation through transcription factor rewiring in microbes is shaped by levels of transcription factor activity, expression, and existing connectivity

The survival of a population during environmental shifts depends on whether the rate of phenotypic adaptation keeps up with the rate of changing conditions. A common way to achieve this is via change to gene regulatory network (GRN) connections-known as rewiring-that facilitate novel interactions and innovation of transcription factors. To understand the success of rapidly adapting organisms, we therefore need to determine the rules that create and constrain opportunities for GRN rewiring. Here, using an experimental microbial model system with the soil bacterium Pseudomonas fluorescens, we reveal a hierarchy among transcription factors that are rewired to rescue lost function, with alternative rewiring pathways only unmasked after the preferred pathway is eliminated. We identify 3 key properties-high activation, high expression, and preexisting low-level affinity for novel target genes-that facilitate transcription factor innovation. Ease of acquiring these properties is constrained by preexisting GRN architecture, which was overcome in our experimental system by both targeted and global network alterations. This work reveals the key properties that determine transcription factor evolvability, and as such, the evolution of GRNs.

Mutational hotspots lead to robust but suboptimal adaptive outcomes in certain environments

The observed mutational spectrum of adaptive outcomes can be constrained by many factors. For example, mutational biases can narrow the observed spectrum by increasing the rate of mutation at isolated sites in the genome. In contrast, complex environments can shift the observed spectrum by defining fitness consequences of mutational routes. We investigate the impact of different nutrient environments on the evolution of motility in Pseudomonas fluorescens Pf0-2x (an engineered non-motile derivative of Pf0-1) in the presence and absence of a strong mutational hotspot. Previous work has shown that this mutational hotspot can be built and broken via six silent mutations, which provide rapid access to a mutation that rescues swimming motility and confers the strongest swimming phenotype in specific environments. Here, we evolved a hotspot and non-hotspot variant strain of Pf0-2x for motility under nutrient-rich (LB) and nutrient-limiting (M9) environmental conditions. We observed the hotspot strain consistently evolved faster across all environmental conditions and its mutational spectrum was robust to environmental differences. However, the non-hotspot strain had a distinct mutational spectrum that changed depending on the nutrient environment. Interestingly, while alternative adaptive mutations in nutrient-rich environments were equal to, or less effective than, the hotspot mutation, the majority of these mutations in nutrient-limited conditions produced superior swimmers. Our competition experiments mirrored these findings, underscoring the role of environment in defining both the mutational spectrum and the associated phenotype strength. This indicates that while mutational hotspots working in concert with natural selection can speed up access to robust adaptive mutations (which can provide a competitive advantage in evolving populations), they can limit exploration of the mutational landscape, restricting access to potentially stronger phenotypes in specific environments.

Multispecies interactions shape the transition to multicellularity

The origin of multicellularity transformed the adaptive landscape on Earth, opening diverse avenues for further innovation. The transition to multicellular life is understood as the evolution of cooperative groups which form a new level of individuality. Despite the potential for community-level interactions, most studies have not addressed the competitive context of this transition, such as competition between species. Here, we explore how interspecific competition shapes the emergence of multicellularity in an experimental system with two yeast species, Saccharomyces cerevisiae and Kluyveromyces lactis, where multicellularity evolves in response to selection for faster settling ability. We find that the multispecies context slows the rate of the transition to multicellularity, and the transition to multicellularity significantly impacts community composition. Multicellular K. lactis emerges first and sweeps through populations in monocultures faster than in cocultures with S. cerevisiae. Following the transition, the between-species competitive dynamics shift, likely in part to intraspecific cooperation in K. lactis. Hence, we document an eco-evolutionary feedback across the transition to multicellularity, underscoring how ecological context is critical for understanding the causes and consequences of innovation. By including two species, we demonstrate that cooperation and competition across several biological scales shapes the origin and persistence of multicellularity.

Unveiling metabolic pathways involved in the extreme desiccation tolerance of an Atacama cyanobacterium

Gloeocapsopsis dulcis strain AAB1 is an extremely xerotolerant cyanobacterium isolated from the Atacama Desert (i.e., the driest and oldest desert on Earth) that holds astrobiological significance due to its ability to biosynthesize compatible solutes at ultra-low water activities. We sequenced and assembled the G. dulcis genome de novo using a combination of long- and short-read sequencing, which resulted in high-quality consensus sequences of the chromosome and two plasmids. We leveraged the G. dulcis genome to generate a genome-scale metabolic model (iGd895) to simulate growth in silico. iGd895 represents, to our knowledge, the first genome-scale metabolic reconstruction developed for an extremely xerotolerant cyanobacterium. The model's predictive capability was assessed by comparing the in silico growth rate with in vitro growth rates of G. dulcis, in addition to the synthesis of trehalose. iGd895 allowed us to explore simulations of key metabolic processes such as essential pathways for water-stress tolerance, and significant alterations to reaction flux distribution and metabolic network reorganization resulting from water limitation. Our study provides insights into the potential metabolic strategies employed by G. dulcis, emphasizing the crucial roles of compatible solutes, metabolic water, energy conservation, and the precise regulation of reaction rates in their adaptation to water stress.

Regulatory response to a hybrid ancestral nitrogenase in Azotobacter vinelandii

Biological nitrogen fixation, the microbial reduction of atmospheric nitrogen to bioavailable ammonia, represents both a major limitation on biological productivity and a highly desirable engineering target for synthetic biology. However, the engineering of nitrogen fixation requires an integrated understanding of how the gene regulatory dynamics of host diazotrophs respond across sequence-function space of its central catalytic metalloenzyme, nitrogenase. Here, we interrogate this relationship by analyzing the transcriptome of Azotobacter vinelandii engineered with a phylogenetically inferred ancestral nitrogenase protein variant. The engineered strain exhibits reduced cellular nitrogenase activity but recovers wild-type growth rates following an extended lag period. We find that expression of genes within the immediate nitrogen fixation network is resilient to the introduced nitrogenase sequence-level perturbations. Rather the sustained physiological compatibility with the ancestral nitrogenase variant is accompanied by reduced expression of genes that support trace metal and electron resource allocation to nitrogenase. Our results spotlight gene expression changes in cellular processes adjacent to nitrogen fixation as productive engineering considerations to improve compatibility between remodeled nitrogenase proteins and engineered host diazotrophs. IMPORTANCE Azotobacter vinelandii is a key model bacterium for the study of biological nitrogen fixation, an important metabolic process catalyzed by nitrogenase enzymes. Here, we demonstrate that compatibilities between engineered A. vinelandii strains and nitrogenase variants can be modulated at the regulatory level. The engineered strain studied here responds by adjusting the expression of proteins involved in cellular processes adjacent to nitrogen fixation, rather than that of nitrogenase proteins themselves. These insights can inform future strategies to transfer nitrogenase variants to non-native hosts.

High quality de novo genome assembly of the non-conventional yeast Kazachstania bulderi describes a potential low pH production host for biorefineries

Kazachstania bulderi is a non-conventional yeast species able to grow efficiently on glucose and δ-gluconolactone at low pH. These unique traits make K. bulderi an ideal candidate for use in sustainable biotechnology processes including low pH fermentations and the production of green chemicals including organic acids. To accelerate strain development with this species, detailed information of its genetics is needed. Here, by employing long read sequencing we report a high-quality phased genome assembly for three strains of K. bulderi species, including the type strain. The sequences were assembled into 12 chromosomes with a total length of 14 Mb, and the genome was fully annotated at structural and functional levels, including allelic and structural variants, ribosomal array and mating type locus. This high-quality reference genome provides a resource to advance our fundamental knowledge of biotechnologically relevant non-conventional yeasts and to support the development of genetic tools for manipulating such strains towards their use as production hosts in biotechnological processes.

A generalizable Cas9/sgRNA prediction model using machine transfer learning with small high-quality datasets

The CRISPR/Cas9 nuclease from Streptococcus pyogenes (SpCas9) can be used with single guide RNAs (sgRNAs) as a sequence-specific antimicrobial agent and as a genome-engineering tool. However, current bacterial sgRNA activity models struggle with accurate predictions and do not generalize well, possibly because the underlying datasets used to train the models do not accurately measure SpCas9/sgRNA activity and cannot distinguish on-target cleavage from toxicity. Here, we solve this problem by using a two-plasmid positive selection system to generate high-quality data that more accurately reports on SpCas9/sgRNA cleavage and that separates activity from toxicity. We develop a machine learning architecture (crisprHAL) that can be trained on existing datasets, that shows marked improvements in sgRNA activity prediction accuracy when transfer learning is used with small amounts of high-quality data, and that can generalize predictions to different bacteria. The crisprHAL model recapitulates known SpCas9/sgRNA-target DNA interactions and provides a pathway to a generalizable sgRNA bacterial activity prediction tool that will enable accurate antimicrobial and genome engineering applications.

Grainyhead-like 2 is required for morphological integrity of mouse embryonic stem cells and orderly formation of inner ear-like organoids

Mutations in the transcription factor gene grainyhead-like 2 (GRHL2) are associated with progressive non-syndromic sensorineural deafness autosomal dominant type 28 (DFNA28) in humans. Since complete loss of Grhl2 is lethal in mouse embryos, we studied its role during inner ear pathology and hearing loss in vitro. To this end, we generated different homozygous deletions to knockout Grhl2 in mouse embryonic stem cells (Grhl2-KO ESCs), including some mimicking naturally occurring truncations in the dimerisation domain related to human DFNA28. Under naïve culture conditions, Grhl2-KO cells in suspension were more heterogenous in size and larger than wild-type controls. Adherent Grhl2-KO cells were also larger, with a less uniform shape, flattened, less circular morphology, forming loose monolayer colonies with poorly defined edges. These changes correlated with lower expression of epithelial cadherin Cdh1 but no changes in tight junction markers (Ocln, Tjp2) or other Grhl isoforms (Grhl1, Grhl3). Clonogenicity from single cells, proliferation rates of cell populations and proliferation markers were reduced in Grhl2-KO ESCs. We next induced stepwise directed differentiation of Grhl2-KO ESCs along an otic pathway, giving rise to three-dimensional inner ear-like organoids (IELOs). Quantitative morphometry revealed that Grhl2-KO cells initially formed larger IELOs with a less compacted structure, more eccentric shape and increased surface area. These morphological changes persisted for up to one week. They were partially rescued by forced cell aggregation and fully restored by stably overexpressing exogenous Grhl2 in Grhl2-KO ESCs, indicating that Grhl2 alters cell-cell interactions. On day 8, aggregates were transferred into minimal maturation medium to allow self-guided organogenesis for another two weeks. During this period, Grhl2-KO cells and wild-type controls developed similarly, expressing neural, neuronal and sensory hair cell markers, while maintaining their initial differences in size and shape. In summary, Grhl2 is required for morphological maintenance of ESCs and orderly formation of IELOs, consistent with an essential role in organising epithelial integrity during inner ear development. Our findings validate quantitative morphometry as a useful, non-invasive screening method for molecular phenotyping of candidate mutations during organoid development.

Metabolic resource overlap impacts competition among phyllosphere bacteria

The phyllosphere is densely colonised by microbial communities, despite sparse and heterogeneously distributed resources. The limitation of resources is expected to drive bacterial competition resulting in exclusion or coexistence based on fitness differences and resource overlap between individual colonisers. We studied the impact of resource competition by determining the effects of different bacterial colonisers on the growth of the model epiphyte Pantoea eucalypti 299R (Pe299R). Resource overlap was predicted based on genome-scale metabolic modelling. By combining results of metabolic modelling and pairwise competitions in the Arabidopsis thaliana phyllosphere and in vitro, we found that ten resources sufficed to explain fitness of Pe299R. An effect of both resource overlap and phylogenetic relationships was found on competition outcomes in vitro as well as in the phyllosphere. However, effects of resource competition were much weaker in the phyllosphere when compared to in vitro experiments. When investigating growth dynamics and reproductive success at the single-cell resolution, resource overlap and phylogenetic relationships are only weakly correlated with epiphytic Pe299R reproductive success, indicating that the leaf's spatial heterogeneity mitigates resource competition. Although the correlation is weak, the presence of competitors led to the development of Pe299R subpopulations that experienced different life histories and cell divisions. In some in planta competitions, Pe299R benefitted from the presence of epiphytes despite high resource overlap to the competitor strain suggesting other factors having stronger effects than resource competition. This study provides fundamental insights into how bacterial communities are shaped in heterogeneous environments and a framework to predict competition outcomes.

Airway metabolic profiling during Streptococcus pneumoniae infection identifies branched chain amino acids as signatures of upper airway colonisation

Streptococcus pneumoniae is a leading cause of community-acquired pneumonia and bacteraemia and is capable of remarkable phenotypic plasticity, responding rapidly to environmental change. Pneumococcus is a nasopharyngeal commensal, but is responsible for severe, acute infections following dissemination within-host. Pneumococcus is adept at utilising host resources, but the airways are compartmentalised and those resources are not evenly distributed. Challenges and opportunities in metabolite acquisition within different airway niches may contribute to the commensal-pathogen switch when pneumococcus moves from nasopharynx into lungs. We used NMR to characterise the metabolic landscape of the mouse airways, in health and during infection. Using paired nasopharynx and lung samples from naïve animals, we identified fundamental differences in metabolite bioavailability between airway niches. Pneumococcal pneumonia was associated with rapid and dramatic shifts in the lung metabolic environment, whilst nasopharyngeal carriage led to only modest change in upper airway metabolite profiles. NMR spectra derived from the nasopharynx of mice infected with closely-related pneumococcal strains that differ in their colonisation potential could be distinguished from one another using multivariate dimensionality reduction methods. The resulting models highlighted that increased branched-chain amino acid (BCAA) bioavailability in nasopharynx is a feature of infection with the high colonisation potential strain. Subsequent analysis revealed increased expression of BCAA transport genes and increased intracellular concentrations of BCAA in that same strain. Movement from upper to lower airway environments is associated with shifting challenges in metabolic resource allocation for pneumococci. Efficient biosynthesis, liberation or acquisition of BCAA is a feature of adaptation to nasopharyngeal colonisation.

Evidence of host specificity in Lactobacillus johnsonii genomes and its influence on probiotic potential in poultry

To date, the selection of candidate strains for probiotic development in production animals has been largely based upon screens for desired phenotypic traits. However, increasing evidence indicates that the use of host-specific strains may be important, because coevolution with the animal host better prepares a bacterial strain to colonize and succeed in its respective host animal species. This concept was applied to Lactobacillus johnsonii in commercial poultry production because of its previous correlation with enhanced bird performance. Using 204 naturally isolated chicken- and turkey-source L. johnsonii, we demonstrate that there is a strong phylogenetic signal for coevolution with the animal host. These isolates differ phenotypically, even within host source, and these differences can be correlated with certain L. johnsonii phylogenetic clades. In commercial turkey poults, turkey-specific strains with strong in vitro phenotypes performed better early in life than strains lacking those phenotypes. A follow-up performance trial in broiler chickens demonstrated that chicken-specific strains result in better overall bird performance than nonchicken-specific strains. Collectively, this work provides evidence for the impact of host adaptation on a probiotic strain's potential. Furthermore, this top-down approach is useful for screening larger numbers of isolates for probiotic candidates.

Transcriptome-guided metabolic network analysis reveals rearrangements of carbon flux distribution in Neisseria gonorrhoeae during neutrophil co-culture

The ability of bacterial pathogens to metabolically adapt to the environmental conditions of their hosts is critical to both colonization and invasive disease. Infection with Neisseria gonorrhoeae (the gonococcus, Gc) is characterized by the influx of neutrophils [polymorphonuclear leukocytes (PMNs)], which fail to clear the bacteria and make antimicrobial products that can exacerbate tissue damage. The inability of the human host to clear Gc infection is particularly concerning in light of the emergence of strains that are resistant to all clinically recommended antibiotics. Bacterial metabolism represents a promising target for the development of new therapeutics against Gc. Here, we generated a curated genome-scale metabolic network reconstruction (GENRE) of Gc strain FA1090. This GENRE links genetic information to metabolic phenotypes and predicts Gc biomass synthesis and energy consumption. We validated this model with published data and in new results reported here. Contextualization of this model using the transcriptional profile of Gc exposed to PMNs revealed substantial rearrangements of Gc central metabolism and induction of Gc nutrient acquisition strategies for alternate carbon source use. These features enhanced the growth of Gc in the presence of neutrophils. From these results, we conclude that the metabolic interplay between Gc and PMNs helps define infection outcomes. The use of transcriptional profiling and metabolic modeling to reveal new mechanisms by which Gc persists in the presence of PMNs uncovers unique aspects of metabolism in this fastidious bacterium, which could be targeted to block infection and thereby reduce the burden of gonorrhea in the human population. IMPORTANCE The World Health Organization designated Gc as a high-priority pathogen for research and development of new antimicrobials. Bacterial metabolism is a promising target for new antimicrobials, as metabolic enzymes are widely conserved among bacterial strains and are critical for nutrient acquisition and survival within the human host. Here we used genome-scale metabolic modeling to characterize the core metabolic pathways of this fastidious bacterium and to uncover the pathways used by Gc during culture with primary human immune cells. These analyses revealed that Gc relies on different metabolic pathways during co-culture with human neutrophils than in rich media. Conditionally essential genes emerging from these analyses were validated experimentally. These results show that metabolic adaptation in the context of innate immunity is important to Gc pathogenesis. Identifying the metabolic pathways used by Gc during infection can highlight new therapeutic targets for drug-resistant gonorrhea.

Tools for Genetic Engineering and Gene Expression Control in Novosphingobium aromaticivorans and Rhodobacter sphaeroides

Alphaproteobacteria have a variety of cellular and metabolic features that provide important insights into biological systems and enable biotechnologies. For example, some species are capable of converting plant biomass into valuable biofuels and bioproducts have the potential to form the backbone of the sustainable bioeconomy. Among the Alphaproteobacteria, Novosphingobium aromaticivorans, Rhodobacter sphaeroides, and Zymomonas mobilis, show particular promise as organisms that can be engineered to convert extracted plant lignin or sugars into bioproducts and biofuels. Genetic manipulation of these bacteria is needed to introduce engineered pathways and modulate expression of native genes with the goal of enhancing bioproduct output. Although recent work has expanded the genetic toolkit for Z. mobilis, N. aromaticivorans and R. sphaeroides still need facile, reliable approaches to deliver genetic payloads to the genome and to control gene expression. Here, we expand the platform of genetic tools for N. aromaticivorans and R. sphaeroides to address these issues. We demonstrate that Tn7 transposition is an effective approach for introducing engineered DNA into the chromosome of N. aromaticivorans and R. sphaeroides. We screen a synthetic promoter library to identify inducible promoters with strong, regulated activity in both organisms. Combining Tn7 integration with promoters from our library, we establish CRISPR interference systems for N. aromaticivorans and R. sphaeroides that can target essential genes and modulate engineered pathways. We anticipate that these systems will greatly facilitate both genetic engineering and gene function discovery efforts in these industrially important species and other Alphaproteobacteria.

Adaptation versus plastic responses to temperature, light, and nitrate availability in cultured snow algal strains

Snow algal blooms are widespread, dominating low temperature, high light, and oligotrophic melting snowpacks. Here, we assessed the photophysiological and cellular stoichiometric responses of snow algal genera Chloromonas spp. and Microglena spp. in their vegetative life stage isolated from the Arctic and Antarctic to gradients in temperature (5 - 15°C), nitrate availability (1 - 10 µmol L-1), and light (50 and 500 µmol photons m-2 s-1). When grown under gradients in temperature, measured snow algal strains displayed Fv/Fm values increased by ∼115% and electron transport rates decreased by ∼50% at 5°C compared to 10 and 15°C, demonstrating how low temperatures can mimic high light impacts to photophysiology. When using carrying capacity as opposed to growth rate as a metric for determining the temperature optima, these snow algal strains can be defined as psychrophilic, with carrying capacities ∼90% higher at 5°C than warmer temperatures. All strains approached Redfield C:N stoichiometry when cultured under nutrient replete conditions regardless of temperature (5.7 ± 0.4 across all strains), whereas significant increases in C:N were apparent when strains were cultured under nitrate concentrations that reflected in situ conditions (17.8 ± 5.9). Intra-specific responses in photophysiology were apparent under high light with Chloromonas spp. more capable of acclimating to higher light intensities. These findings suggest that in situ conditions are not optimal for the studied snow algal strains, but they are able to dynamically adjust both their photochemistry and stoichiometry to acclimate to these conditions.

Rapid Evolution of Multidrug Resistance in a Candida lusitaniae Infection during Micafungin Monotherapy

Candida (Clavispora) lusitaniae is a rare, emerging non-albicans Candida species that can cause life-threatening invasive infections, spread within hospital settings, and rapidly acquire antifungal drug resistance, including multidrug resistance. The frequency and spectrum of mutations causing antifungal drug resistance in C. lusitaniae are poorly understood. Analyses of serial clinical isolates of any Candida species are uncommon and often analyze a limited number of samples collected over months of antifungal therapy with multiple drug classes, limiting the ability to understand relationships between drug classes and specific mutations. Here, we performed comparative genomic and phenotypic analysis of 20 serial C. lusitaniae bloodstream isolates collected daily from an individual patient treated with micafungin monotherapy during a single 11-day hospital admission. We identified isolates with decreased micafungin susceptibility 4 days after initiation of antifungal therapy and a single isolate with increased cross-resistance to micafungin and fluconazole, despite no history of azole therapy in this patient. Only 14 unique single nucleotide polymorphisms (SNPs) were identified between all 20 samples, including three different FKS1 alleles among isolates with decreased micafungin susceptibility and an ERG3 missense mutation found only in the isolate with increased cross-resistance to both micafungin and fluconazole. This is the first clinical evidence of an ERG3 mutation in C. lusitaniae that occurred during echinocandin monotherapy and is associated with cross-resistance to multiple drug classes. Overall, the evolution of multidrug resistance in C. lusitaniae is rapid and can emerge during treatment with only first-line antifungal therapy.

How Do Prebiotics Affect Human Intestinal Bacteria?-Assessment of Bacterial Growth with Inulin and XOS In Vitro

Prebiotics are believed to exhibit high specificity in stimulating the growth or activity of a limited number of commensal microorganisms, thereby conferring health benefits to the host. However, the mechanism of action of prebiotics depends on multiple factors, including the composition of an individual's gut microbiota, and is therefore difficult to predict. It is known that different bacteria can utilize inulin and xylooligosaccharides (XOS), but an overview of which bacteria in the human gut may be affected is lacking. Detailed knowledge of how bacterial growth is affected by prebiotics is furthermore useful for the development of new synbiotics, which combine a living microorganism with a selective substrate to confer a health benefit to the host. Hence, we developed a statistical model to compare growth in vitro among typical human gut bacteria from different phylogenetic lineages. Based on continuous observation of the optical density (OD600), we compare maximal growth rates (rmax), maximal attained OD600 (ODmax), and area under the growth curve (AUC) of bacteria grown on inulin or XOS. The consideration of these three parameters suggests strain-specific preferences for inulin or XOS and reveals previously unknown preferences such as Streptococcus salivarius growth on XOS.

The conserved iol gene cluster in Pseudomonas is involved in rhizosphere competence

The Pseudomonas genus has shown great potential as a sustainable solution to support agriculture through its plant-growth-promoting and biocontrol activities. However, their efficacy as bioinoculants is limited by unpredictable colonization in natural conditions. Our study identifies the iol locus, a gene cluster in Pseudomonas involved in inositol catabolism, as a feature enriched among superior root colonizers in natural soil. Further characterization revealed that the iol locus increases competitiveness, potentially caused by an observed induction of swimming motility and the production of fluorescent siderophore in response to inositol, a plant-derived compound. Public data analyses indicate that the iol locus is broadly conserved in the Pseudomonas genus and linked to diverse host-microbe interactions. Together, our findings suggest the iol locus as a potential target for developing more effective bioinoculants for sustainable agriculture.

SIMBA Method-Simultaneous Detection of Antimicrobial and Anti-biofilm Activity of New Compounds Using Salmonella Infantis

The development of antimicrobial resistance and the formation of Salmonella biofilms are serious public health problems. For this reason, new natural compounds with antimicrobial and anti-biofilm activity are being sought, and wild fungi represent an untapped potential. Various extraction agents, including organic solvents and aqueous buffers, can be used to obtain bioactive compounds from natural sources. To evaluate their bioactivity, extensive screening studies are required to determine antimicrobial and anti-biofilm activity using methods such as broth microdilution or crystal violet assay, respectively, but none of these methods allow simultaneous evaluation of both activities against bacteria. Cold water extraction from wild fungi offers the advantage of extracting water-soluble compounds. The SIMultaneous detection of antiMicrobial and anti-Biofilm Activity (SIMBA) method combines the testing of both types of activity against bacteria with the evaluation of the 20 h growth curve of the Salmonella Infantis ŽM9 strain determined with absorbance measurements at 600 nm in a 96-well plate. SIMBA method thus shortens the time to determine the bioactivity of extracts, reduces material consumption, and eliminates the need for additional reagents. SIMBA enables rapid selection of bioactive extracts for their fractionation and shortens the time to determine new natural products with antimicrobial and anti-biofilm activity. Graphical overview.

The Ypk1 protein kinase signaling pathway is rewired and not essential for viability in Candida albicans

Protein kinases are central components of almost all signaling pathways that control cellular activities. In the model organism Saccharomyces cerevisiae, the paralogous protein kinases Ypk1 and Ypk2, which control membrane lipid homeostasis, are essential for viability, and previous studies strongly indicated that this is also the case for their single ortholog Ypk1 in the pathogenic yeast Candida albicans. Here, using FLP-mediated inducible gene deletion, we reveal that C. albicans ypk1Δ mutants are viable but slow-growing, explaining prior failures to obtain null mutants. Phenotypic analyses of the mutants showed that the functions of Ypk1 in regulating sphingolipid biosynthesis and cell membrane lipid asymmetry are conserved, but the consequences of YPK1 deletion are milder than in S. cerevisiae. Mutational studies demonstrated that the highly conserved PDK1 phosphorylation site T548 in its activation loop is essential for Ypk1 function, whereas the TORC2 phosphorylation sites S687 and T705 at the C-terminus are important for Ypk1-dependent resistance to membrane stress. Unexpectedly, Pkh1, the single C. albicans orthologue of Pkh1/Pkh2, which mediate Ypk1 phosphorylation at the PDK1 site in S. cerevisiae, was not required for normal growth of C. albicans under nonstressed conditions, and Ypk1 phosphorylation at T548 was only slightly reduced in pkh1Δ mutants. We found that another protein kinase, Pkh3, whose ortholog in S. cerevisiae cannot substitute Pkh1/2, acts redundantly with Pkh1 to activate Ypk1 in C. albicans. No phenotypic effects were observed in cells lacking Pkh3 alone, but pkh1Δ pkh3Δ double mutants had a severe growth defect and Ypk1 phosphorylation at T548 was completely abolished. These results establish that Ypk1 is not essential for viability in C. albicans and that, despite its generally conserved function, the Ypk1 signaling pathway is rewired in this pathogenic yeast and includes a novel upstream kinase to activate Ypk1 by phosphorylation at the PDK1 site.

Solving the stochastic dynamics of population growth

Population growth is a fundamental process in ecology and evolution. The population size dynamics during growth are often described by deterministic equations derived from kinetic models. Here, we simulate several population growth models and compare the size averaged over many stochastic realizations with the deterministic predictions. We show that these deterministic equations are generically bad predictors of the average stochastic population dynamics. Specifically, deterministic predictions overestimate the simulated population sizes, especially those of populations starting with a small number of individuals. Describing population growth as a stochastic birth process, we prove that the discrepancy between deterministic predictions and simulated data is due to unclosed-moment dynamics. In other words, the deterministic approach does not consider the variability of birth times, which is particularly important with small population sizes. We show that some moment-closure approximations describe the growth dynamics better than the deterministic prediction. However, they do not reduce the error satisfactorily and only apply to some population growth models. We explicitly solve the stochastic growth dynamics, and our solution applies to any population growth model. We show that our solution exactly quantifies the dynamics of a community composed of different strains and correctly predicts the fixation probability of a strain in a serial dilution experiment. Our work sets the foundations for a more faithful modeling of community and population dynamics. It will allow the development of new tools for a more accurate analysis of experimental and empirical results, including the inference of important growth parameters.