High-throughput protein characterization by complementation using DNA barcoded fragment libraries.

Our ability to predict, control, or design biological function is fundamentally limited by poorly annotated gene function. This can be particularly challenging in non-model systems. Accordingly, there is motivation for new high-throughput methods for accurate functional annotation. Here, we used complementation of auxotrophs and DNA barcode sequencing (Coaux-Seq) to enable high-throughput characterization of protein function. Fragment libraries from eleven genetically diverse bacteria were tested in twenty different auxotrophic strains of Escherichia coli to identify genes that complement missing biochemical activity. We recovered 41% of expected hits, with effectiveness ranging per source genome, and observed success even with distant E. coli relatives like Bacillus subtilis and Bacteroides thetaiotaomicron. Coaux-Seq provided the first experimental validation for 53 proteins, of which 11 are less than 40% identical to an experimentally characterized protein. Among the unexpected function identified was a sulfate uptake transporter, an O-succinylhomoserine sulfhydrylase for methionine synthesis, and an aminotransferase. We also identified instances of cross-feeding wherein protein overexpression and nearby non-auxotrophic strains enabled growth. Altogether, Coaux-Seq's utility is demonstrated, with future applications in ecology, health, and engineering.

Potential applications of microbial genomics in nuclear non-proliferation.

As nuclear technology evolves in response to increased demand for diversification and decarbonization of the energy sector, new and innovative approaches are needed to effectively identify and deter the proliferation of nuclear arms, while ensuring safe development of global nuclear energy resources. Preventing the use of nuclear material and technology for unsanctioned development of nuclear weapons has been a long-standing challenge for the International Atomic Energy Agency and signatories of the Treaty on the Non-Proliferation of Nuclear Weapons. Environmental swipe sampling has proven to be an effective technique for characterizing clandestine proliferation activities within and around known locations of nuclear facilities and sites. However, limited tools and techniques exist for detecting nuclear proliferation in unknown locations beyond the boundaries of declared nuclear fuel cycle facilities, representing a critical gap in non-proliferation safeguards. Microbiomes, defined as "characteristic communities of microorganisms" found in specific habitats with distinct physical and chemical properties, can provide valuable information about the conditions and activities occurring in the surrounding environment. Microorganisms are known to inhabit radionuclide-contaminated sites, spent nuclear fuel storage pools, and cooling systems of water-cooled nuclear reactors, where they can cause radionuclide migration and corrosion of critical structures. Microbial transformation of radionuclides is a well-established process that has been documented in numerous field and laboratory studies. These studies helped to identify key bacterial taxa and microbially-mediated processes that directly and indirectly control the transformation, mobility, and fate of radionuclides in the environment. Expanding on this work, other studies have used microbial genomics integrated with machine learning models to successfully monitor and predict the occurrence of heavy metals, radionuclides, and other process wastes in the environment, indicating the potential role of nuclear activities in shaping microbial community structure and function. Results of this previous body of work suggest fundamental geochemical-microbial interactions occurring at nuclear fuel cycle facilities could give rise to microbiomes that are characteristic of nuclear activities. These microbiomes could provide valuable information for monitoring nuclear fuel cycle facilities, planning environmental sampling campaigns, and developing biosensor technology for the detection of undisclosed fuel cycle activities and proliferation concerns.

Fitness factors impacting survival of a subsurface bacterium in contaminated groundwater.

Many factors contribute to the ability of a microbial species to persist when encountering complexly contaminated environments including time of exposure, the nature and concentration of contaminants, availability of nutritional resources, and possession of a combination of appropriate molecular mechanisms needed for survival. Herein we sought to identify genes that are most important for survival of Gram-negative Enterobacteriaceae in contaminated groundwater environments containing high concentrations of nitrate and metals using the metal-tolerant Oak Ridge Reservation (ORR) isolate, Pantoea sp. MT58 (MT58). Survival fitness experiments in which a randomly barcoded transposon insertion (RB-TnSeq) library of MT58 was exposed directly to contaminated ORR groundwater samples from across a nitrate and mixed metal contamination plume were used to identify genes important for survival with increasing exposure times and concentrations of contaminants, and availability of a carbon source. Genes involved in controlling and using carbon, encoding transcriptional regulators, and related to Gram-negative outer membrane processes were among those found to be important for survival in contaminated ORR groundwater. A comparative genomics analysis of 75 Pantoea genus strains allowed us to further separate the survival determinants into core and non-core genes in the Pantoea pangenome, revealing insights into the survival of subsurface microorganisms during contaminant plume intrusion.

Interactive tools for functional annotation of bacterial genomes.

Automated annotations of protein functions are error-prone because of our lack of knowledge of protein functions. For example, it is often impossible to predict the correct substrate for an enzyme or a transporter. Furthermore, much of the knowledge that we do have about the functions of proteins is missing from the underlying databases. We discuss how to use interactive tools to quickly find different kinds of information relevant to a protein's function. Many of these tools are available via PaperBLAST (http://papers.genomics.lbl.gov). Combining these tools often allows us to infer a protein's function. Ideally, accurate annotations would allow us to predict a bacterium's capabilities from its genome sequence, but in practice, this remains challenging. We describe interactive tools that infer potential capabilities from a genome sequence or that search a genome to find proteins that might perform a specific function of interest. Database URL: http://papers.genomics.lbl.gov.

Nitrogen accountancy in space agriculture.

Food production and pharmaceutical synthesis are posited as essential biotechnologies for facilitating human exploration beyond Earth. These technologies not only offer critical green space and food agency to astronauts but also promise to minimize mass and volume requirements through scalable, modular agriculture within closed-loop systems, offering an advantage over traditional bring-along strategies. Despite these benefits, the prevalent model for evaluating such systems exhibits significant limitations. It lacks comprehensive inventory and mass balance analyses for crop cultivation and life support, and fails to consider the complexities introduced by cultivating multiple crop varieties, which is crucial for enhancing food diversity and nutritional value. Here we expand space agriculture modeling to account for nitrogen dependence across an array of crops and demonstrate our model with experimental fitting of parameters. By adding nitrogen limitations, an extended model can account for potential interruptions in feedstock supply. Furthermore, sensitivity analysis was used to distill key consequential parameters that may be the focus of future experimental efforts.

High-throughput genetics enables identification of nutrient utilization and accessory energy metabolism genes in a model methanogen.

Archaea are widespread in the environment and play fundamental roles in diverse ecosystems; however, characterization of their unique biology requires advanced tools. This is particularly challenging when characterizing gene function. Here, we generate randomly barcoded transposon libraries in the model methanogenic archaeon Methanococcus maripaludis and use high-throughput growth methods to conduct fitness assays (RB-TnSeq) across over 100 unique growth conditions. Using our approach, we identified new genes involved in nutrient utilization and response to oxidative stress. We identified novel genes for the usage of diverse nitrogen sources in M. maripaludis including a putative regulator of alanine deamination and molybdate transporters important for nitrogen fixation. Furthermore, leveraging the fitness data, we inferred that M. maripaludis can utilize additional nitrogen sources including ʟ-glutamine, ᴅ-glucuronamide, and adenosine. Under autotrophic growth conditions, we identified a gene encoding a domain of unknown function (DUF166) that is important for fitness and hypothesize that it has an accessory role in carbon dioxide assimilation. Finally, comparing fitness costs of oxygen versus sulfite stress, we identified a previously uncharacterized class of dissimilatory sulfite reductase-like proteins (Dsr-LP; group IIId) that is important during growth in the presence of sulfite. When overexpressed, Dsr-LP conferred sulfite resistance and enabled use of sulfite as the sole sulfur source. The high-throughput approach employed here allowed for generation of a large-scale data set that can be used as a resource to further understand gene function and metabolism in the archaeal domain.IMPORTANCEArchaea are widespread in the environment, yet basic aspects of their biology remain underexplored. To address this, we apply randomly barcoded transposon libraries (RB-TnSeq) to the model archaeon Methanococcus maripaludis. RB-TnSeq coupled with high-throughput growth assays across over 100 unique conditions identified roles for previously uncharacterized genes, including several encoding proteins with domains of unknown function (DUFs). We also expand on our understanding of carbon and nitrogen metabolism and characterize a group IIId dissimilatory sulfite reductase-like protein as a functional sulfite reductase. This data set encompasses a wide range of additional conditions including stress, nitrogen fixation, amino acid supplementation, and autotrophy, thus providing an extensive data set for the archaeal community to mine for characterizing additional genes of unknown function.

Molecular mechanisms and environmental adaptations of flagellar loss and biofilm growth of Rhodanobacter under environmental stress.

Biofilms aid bacterial adhesion to surfaces via direct and indirect mechanisms, and formation of biofilms is considered as an important strategy for adaptation and survival in suboptimal environmental conditions. However, the molecular underpinnings of biofilm formation in subsurface sediment/groundwater ecosystems where microorganisms often experience fluctuations in nutrient input, pH, and nitrate or metal concentrations are underexplored. We examined biofilm formation under different nutrient, pH, metal, and nitrate regimens of 16 Rhodanobacter strains isolated from subsurface groundwater wells spanning diverse levels of pH (3.5 to 5) and nitrates (13.7 to 146 mM). Eight Rhodanobacter strains demonstrated significant biofilm growth under low pH, suggesting adaptations for survival and growth at low pH. Biofilms were intensified under aluminum stress, particularly in strains possessing fewer genetic traits associated with biofilm formation, findings warranting further investigation. Through random barcode transposon-site sequencing (RB-TnSeq), proteomics, use of specific mutants, and transmission electron microscopy analysis, we discovered flagellar loss under aluminum stress, indicating a potential relationship between motility, metal tolerance, and biofilm growth. Comparative genomic analyses revealed the absence of flagella and chemotaxis genes and the presence of a putative type VI secretion system in the highly biofilm-forming strain FW021-MT20. In this study we identified genetic determinants associated with biofilm growth under metal stress in a predominant environmental genus, Rhodanobacter, and identified traits aiding survival and adaptation to contaminated subsurface environments.

Barcoded overexpression screens in gut Bacteroidales identify genes with roles in carbon utilization and stress resistance.

A mechanistic understanding of host-microbe interactions in the gut microbiome is hindered by poorly annotated bacterial genomes. While functional genomics can generate large gene-to-phenotype datasets to accelerate functional discovery, their applications to study gut anaerobes have been limited. For instance, most gain-of-function screens of gut-derived genes have been performed in Escherichia coli and assayed in a small number of conditions. To address these challenges, we develop Barcoded Overexpression BActerial shotgun library sequencing (Boba-seq). We demonstrate the power of this approach by assaying genes from diverse gut Bacteroidales overexpressed in Bacteroides thetaiotaomicron. From hundreds of experiments, we identify new functions and phenotypes for 29 genes important for carbohydrate metabolism or tolerance to antibiotics or bile salts. Highlights include the discovery of a D-glucosamine kinase, a raffinose transporter, and several routes that increase tolerance to ceftriaxone and bile salts through lipid biosynthesis. This approach can be readily applied to develop screens in other strains and additional phenotypic assays.

Mixed waste contamination selects for a mobile genetic element population enriched in multiple heavy metal resistance genes.

Mobile genetic elements (MGEs) like plasmids, viruses, and transposable elements can provide fitness benefits to their hosts for survival in the presence of environmental stressors. Heavy metal resistance genes (HMRGs) are frequently observed on MGEs, suggesting that MGEs may be an important driver of adaptive evolution in environments contaminated with heavy metals. Here, we report the meta-mobilome of the heavy metal-contaminated regions of the Oak Ridge Reservation subsurface. This meta-mobilome was compared with one derived from samples collected from unimpacted regions of the Oak Ridge Reservation subsurface. We assembled 1615 unique circularized DNA elements that we propose to be MGEs. The circular elements from the highly contaminated subsurface were enriched in HMRG clusters relative to those from the nearby unimpacted regions. Additionally, we found that these HMRGs were associated with Gamma and Betaproteobacteria hosts in the contaminated subsurface and potentially facilitate the persistence and dominance of these taxa in this region. Finally, the HMRGs were associated with conjugative elements, suggesting their potential for future lateral transfer. We demonstrate how our understanding of MGE ecology, evolution, and function can be enhanced through the genomic context provided by completed MGE assemblies.

Contribution of Microorganisms with the Clade II Nitrous Oxide Reductase to Suppression of Surface Emissions of Nitrous Oxide.

The sources and sinks of nitrous oxide, as control emissions to the atmosphere, are generally poorly constrained for most environmental systems. Initial depth-resolved analysis of nitrous oxide flux from observation wells and the proximal surface within a nitrate contaminated aquifer system revealed high subsurface production but little escape from the surface. To better understand the environmental controls of production and emission at this site, we used a combination of isotopic, geochemical, and molecular analyses to show that chemodenitrification and bacterial denitrification are major sources of nitrous oxide in this subsurface, where low DO, low pH, and high nitrate are correlated with significant nitrous oxide production. Depth-resolved metagenomes showed that consumption of nitrous oxide near the surface was correlated with an enrichment of Clade II nitrous oxide reducers, consistent with a growing appreciation of their importance in controlling release of nitrous oxide to the atmosphere. Our work also provides evidence for the reduction of nitrous oxide at a pH of 4, well below the generally accepted limit of pH 5.

A fast comparative genome browser for diverse bacteria and archaea.

Genome sequencing has revealed an incredible diversity of bacteria and archaea, but there are no fast and convenient tools for browsing across these genomes. It is cumbersome to view the prevalence of homologs for a protein of interest, or the gene neighborhoods of those homologs, across the diversity of the prokaryotes. We developed a web-based tool, fast.genomics, that uses two strategies to support fast browsing across the diversity of prokaryotes. First, the database of genomes is split up. The main database contains one representative from each of the 6,377 genera that have a high-quality genome, and additional databases for each taxonomic order contain up to 10 representatives of each species. Second, homologs of proteins of interest are identified quickly by using accelerated searches, usually in a few seconds. Once homologs are identified, fast.genomics can quickly show their prevalence across taxa, view their neighboring genes, or compare the prevalence of two different proteins. Fast.genomics is available at https://fast.genomics.lbl.gov.

Genomic and environmental controls on Castellaniella biogeography in an anthropogenically disturbed subsurface.

Castellaniella species have been isolated from a variety of mixed-waste environments including the nitrate and multiple metal-contaminated subsurface at the Oak Ridge Reservation (ORR). Previous studies examining microbial community composition and nitrate removal at ORR during biostimulation efforts reported increased abundances of members of the Castellaniella genus concurrent with increased denitrification rates. Thus, we asked how genomic and abiotic factors control the Castellaniella biogeography at the site to understand how these factors may influence nitrate transformation in an anthropogenically impacted setting. We report the isolation and characterization of several Castellaniella strains from the ORR subsurface. Five of these isolates match at 100% identity (at the 16S rRNA gene V4 region) to two Castellaniella amplicon sequence variants (ASVs), ASV1 and ASV2, that have persisted in the ORR subsurface for at least 2 decades. However, ASV2 has consistently higher relative abundance in samples taken from the site and was also the dominant blooming denitrifier population during a prior biostimulation effort. We found that the ASV2 representative strain has greater resistance to mixed metal stress than the ASV1 representative strains. We attribute this resistance, in part, to the large number of unique heavy metal resistance genes identified on a genomic island in the ASV2 representative genome. Additionally, we suggest that the relatively lower fitness of ASV1 may be connected to the loss of the nitrous oxide reductase (nos) operon (and associated nitrous oxide reductase activity) due to the insertion at this genomic locus of a mobile genetic element carrying copper resistance genes. This study demonstrates the value of integrating genomic, environmental, and phenotypic data to characterize the biogeography of key microorganisms in contaminated sites.

Environmental stress mediates groundwater microbial community assembly.

Community assembly describes how different ecological processes shape microbial community composition and structure. How environmental factors impact community assembly remains elusive. Here we sampled microbial communities and >200 biogeochemical variables in groundwater at the Oak Ridge Field Research Center, a former nuclear waste disposal site, and developed a theoretical framework to conceptualize the relationships between community assembly processes and environmental stresses. We found that stochastic assembly processes were critical (>60% on average) in shaping community structure, but their relative importance decreased as stress increased. Dispersal limitation and 'drift' related to random birth and death had negative correlations with stresses, whereas the selection processes leading to dissimilar communities increased with stresses, primarily related to pH, cobalt and molybdenum. Assembly mechanisms also varied greatly among different phylogenetic groups. Our findings highlight the importance of microbial dispersal limitation and environmental heterogeneity in ecosystem restoration and management.

Eliminating genes for a two-component system increases PHB productivity in Cupriavidus basilensis 4G11 under PHB suppressing, nonstress conditions.

Species of bacteria from the genus Cupriavidus are known, in part, for their ability to produce high amounts of poly-hydroxybutyrate (PHB) making them attractive candidates for bioplastic production. The native synthesis of PHB occurs during periods of metabolic stress, and the process regulating the initiation of PHB accumulation in these organisms is not fully understood. Screening an RB-TnSeq transposon library of Cupriavidus basilensis 4G11 allowed us to identify two genes of an apparent, uncharacterized two-component system, which when omitted from the genome enable increased PHB productivity in balanced, nonstress growth conditions. We observe average increases in PHB productivity of 56% and 41% relative to the wildtype parent strain upon deleting each gene individually from the genome. The increased PHB phenotype disappears, however, in nitrogen-free unbalanced growth conditions suggesting the phenotype is specific to fast-growing, replete, nonstress growth. Bioproduction modeling suggests this phenotype could be due to a decreased reliance on metabolic stress induced by nitrogen limitation to initiate PHB production in the mutant strains. Due to uncertainty in the two-component system's input signal and regulon, the mechanism by which these genes impart this phenotype remains unclear. Such strains may allow for the use of single-stage, continuous bioreactor systems, which are far simpler than many PHB bioproduction schemes used previously, given a similar product yield to batch systems in such a configuration. Bioproductivity modeling suggests that omitting this regulation in the cells may increase PHB productivity up to 24% relative to the wildtype organism when using single-stage continuous systems. This work expands our understanding of the regulation of PHB accumulation in Cupriavidus, in particular the initiation of this process upon transition into unbalanced growth regimes.

A bacterial sensor taxonomy across earth ecosystems for machine learning applications.

Microbial communities have evolved to colonize all ecosystems of the planet, from the deep sea to the human gut. Microbes survive by sensing, responding, and adapting to immediate environmental cues. This process is driven by signal transduction proteins such as histidine kinases, which use their sensing domains to bind or otherwise detect environmental cues and "transduce" signals to adjust internal processes. We hypothesized that an ecosystem's unique stimuli leave a sensor "fingerprint," able to identify and shed insight on ecosystem conditions. To test this, we collected 20,712 publicly available metagenomes from Host-associated, Environmental, and Engineered ecosystems across the globe. We extracted and clustered the collection's nearly 18M unique sensory domains into 113,712 similar groupings with MMseqs2. We built gradient-boosted decision tree machine learning models and found we could classify the ecosystem type (accuracy: 87%) and predict the levels of different physical parameters (R2 score: 83%) using the sensor cluster abundance as features. Feature importance enables identification of the most predictive sensors to differentiate between ecosystems which can lead to mechanistic interpretations if the sensor domains are well annotated. To demonstrate this, a machine learning model was trained to predict patient's disease state and used to identify domains related to oxygen sensing present in a healthy gut but missing in patients with abnormal conditions. Moreover, since 98.7% of identified sensor domains are uncharacterized, importance ranking can be used to prioritize sensors to determine what ecosystem function they may be sensing. Furthermore, these new predictive sensors can function as targets for novel sensor engineering with applications in biotechnology, ecosystem maintenance, and medicine.IMPORTANCEMicrobes infect, colonize, and proliferate due to their ability to sense and respond quickly to their surroundings. In this research, we extract the sensory proteins from a diverse range of environmental, engineered, and host-associated metagenomes. We trained machine learning classifiers using sensors as features such that it is possible to predict the ecosystem for a metagenome from its sensor profile. We use the optimized model's feature importance to identify the most impactful and predictive sensors in different environments. We next use the sensor profile from human gut metagenomes to classify their disease states and explore which sensors can explain differences between diseases. The sensors most predictive of environmental labels here, most of which correspond to uncharacterized proteins, are a useful starting point for the discovery of important environment signals and the development of possible diagnostic interventions.

Systematic and scalable genome-wide essentiality mapping to identify nonessential genes in phages.

Phages are one of the key ecological drivers of microbial community dynamics, function, and evolution. Despite their importance in bacterial ecology and evolutionary processes, phage genes are poorly characterized, hampering their usage in a variety of biotechnological applications. Methods to characterize such genes, even those critical to the phage life cycle, are labor intensive and are generally phage specific. Here, we develop a systematic gene essentiality mapping method scalable to new phage-host combinations that facilitate the identification of nonessential genes. As a proof of concept, we use an arrayed genome-wide CRISPR interference (CRISPRi) assay to map gene essentiality landscape in the canonical coliphages λ and P1. Results from a single panel of CRISPRi probes largely recapitulate the essential gene roster determined from decades of genetic analysis for lambda and provide new insights into essential and nonessential loci in P1. We present evidence of how CRISPRi polarity can lead to false positive gene essentiality assignments and recommend caution towards interpreting CRISPRi data on gene essentiality when applied to less studied phages. Finally, we show that we can engineer phages by inserting DNA barcodes into newly identified inessential regions, which will empower processes of identification, quantification, and tracking of phages in diverse applications.

Evaluating E. coli genome-scale metabolic model accuracy with high-throughput mutant fitness data.

The Escherichia coli genome-scale metabolic model (GEM) is an exemplar systems biology model for the simulation of cellular metabolism. Experimental validation of model predictions is essential to pinpoint uncertainty and ensure continued development of accurate models. Here, we quantified the accuracy of four subsequent E. coli GEMs using published mutant fitness data across thousands of genes and 25 different carbon sources. This evaluation demonstrated the utility of the area under a precision-recall curve relative to alternative accuracy metrics. An analysis of errors in the latest (iML1515) model identified several vitamins/cofactors that are likely available to mutants despite being absent from the experimental growth medium and highlighted isoenzyme gene-protein-reaction mapping as a key source of inaccurate predictions. A machine learning approach further identified metabolic fluxes through hydrogen ion exchange and specific central metabolism branch points as important determinants of model accuracy. This work outlines improved practices for the assessment of GEM accuracy with high-throughput mutant fitness data and highlights promising areas for future model refinement in E. coli and beyond.

Mixed nitrate and metal contamination influences operational speciation of toxic and essential elements.

Environmental contamination constrains microbial communities impacting diversity and total metabolic activity. The former S-3 Ponds contamination site at Oak Ridge Reservation (ORR), TN, has elevated concentrations of nitric acid and multiple metals from decades of processing nuclear material. To determine the nature of the metal contamination in the sediment, a three-step sequential chemical extraction (BCR) was performed on sediment segments from a core located upgradient (EB271, non-contaminated) and one downgradient (EB106, contaminated) of the S-3 Ponds. The resulting exchangeable, reducing, and oxidizing fractions were analyzed for 18 different elements. Comparison of the two cores revealed changes in operational speciation for several elements caused by the contamination. Those present from the S-3 Ponds, including Al, U, Co, Cu, Ni, and Cd, were not only elevated in concentration in the EB106 core but were also operationally more available with increased mobility in the acidic environment. Other elements, including Mg, Ca, P, V, As, and Mo, were less operationally available in EB106 having decreased concentrations in the exchangeable fraction. The bioavailability of essential macro nutrients Mg, Ca, and P from the two types of sediment was determined using three metal-tolerant bacteria previously isolated from ORR. Mg and Ca were available from both sediments for all three strains; however, P was not bioavailable from either sediment for any strain. The decreased operational speciation of P in contaminated ORR sediment may increase the dependence of the microbial community on other pools of P or select for microorganisms with increased P scavenging capabilities. Hence, the microbial community at the former S-3 Ponds contamination site may be constrained not only by increased toxic metal concentrations but also by the availability of essential elements, including P.

Geochemical constraints on bacteriophage infectivity in terrestrial environments.

Lytic phages can be potent and selective inhibitors of microbial growth and can have profound impacts on microbiome composition and function. However, there is uncertainty about the biogeochemical conditions under which phage predation modulates microbial ecosystem function, particularly in terrestrial systems. Ionic strength is critical for infection of bacteria by many phages, but quantitative data is limited on the ion thresholds for phage infection that can be compared with environmental ion concentrations. Similarly, while carbon composition varies in the environment, we do not know how this variability influences the impact of phage predation on microbiome function. Here, we measured the half-maximal effective concentrations (EC50) of 80 different inorganic ions for the infection of E. coli with two canonical dsDNA and ssRNA phages, T4 and MS2, respectively. Many alkaline earth metals and alkali metals enabled lytic infection but the ionic strength thresholds varied for different ions between phages. Additionally, using a freshwater nitrate-reducing microbiome, we found that the ability of lytic phages to influence nitrate reduction end-products depended upon the carbon source as well as ionic strength. For all phage:host pairs, the ion EC50s for phage infection exceeded the ion concentrations found in many terrestrial freshwater systems. Thus, our findings support a model where phages most influence terrestrial microbial functional ecology in hot spots and hot moments such as metazoan guts, drought influenced soils, or biofilms where ion concentration is locally or transiently elevated and nutrients are available to support the growth of specific phage hosts.