Globally occurring pelagiphage infections create ribosome-deprived cells.

Phages play an essential role in controlling bacterial populations. Those infecting Pelagibacterales (SAR11), the dominant bacteria in surface oceans, have been studied in silico and by cultivation attempts. However, little is known about the quantity of phage-infected cells in the environment. Using fluorescence in situ hybridization techniques, we here show pelagiphage-infected SAR11 cells across multiple global ecosystems and present evidence for tight community control of pelagiphages on the SAR11 hosts in a case study. Up to 19% of SAR11 cells were phage-infected during a phytoplankton bloom, coinciding with a ~90% reduction in SAR11 cell abundance within 5 days. Frequently, a fraction of the infected SAR11 cells were devoid of detectable ribosomes, which appear to be a yet undescribed possible stage during pelagiphage infection. We dubbed such cells zombies and propose, among other possible explanations, a mechanism in which ribosomal RNA is used as a resource for the synthesis of new phage genomes. On a global scale, we detected phage-infected SAR11 and zombie cells in the Atlantic, Pacific, and Southern Oceans. Our findings illuminate the important impact of pelagiphages on SAR11 populations and unveil the presence of ribosome-deprived zombie cells as part of the infection cycle.

Taxonomic and functional stability overrules seasonality in polar benthic microbiomes.

Coastal shelf sediments are hot spots of organic matter mineralization. They receive up to 50% of primary production, which, in higher latitudes, is strongly seasonal. Polar and temperate benthic bacterial communities, however, show a stable composition based on comparative 16S rRNA gene sequencing despite different microbial activity levels. Here, we aimed to resolve this contradiction by identifying seasonal changes at the functional level, in particular with respect to algal polysaccharide degradation genes, by combining metagenomics, metatranscriptomics, and glycan analysis in sandy surface sediments from Isfjorden, Svalbard. Gene expressions of diverse carbohydrate-active enzymes changed between winter and spring. For example, ß-1,3-glucosidases (e.g. GH30, GH17, GH16) degrading laminarin, an energy storage molecule of algae, were elevated in spring, while enzymes related to α-glucan degradation were expressed in both seasons with maxima in winter (e.g. GH63, GH13_18, and GH15). Also, the expression of GH23 involved in peptidoglycan degradation was prevalent, which is in line with recycling of bacterial biomass. Sugar extractions from bulk sediments were low in concentrations during winter but higher in spring samples, with glucose constituting the largest fraction of measured monosaccharides (84% ± 14%). In porewater, glycan concentrations were ~18-fold higher than in overlying seawater (1107 ± 484 vs. 62 ± 101 µg C l-1) and were depleted in glucose. Our data indicate that microbial communities in sandy sediments digest and transform labile parts of photosynthesis-derived particulate organic matter and likely release more stable, glucose-depleted residual glycans of unknown structures, quantities, and residence times into the ocean, thus modulating the glycan composition of marine coastal waters.

Possible Role of CHAD Proteins in Copper Resistance.

Conserved Histidine Alpha-helical Domain (CHAD) proteins attached to the surface of polyphosphate (PolyP) have been studied in some bacteria and one archaeon. However, the activity of CHAD proteins is unknown beyond their interaction with PolyP granules. By using bioinformatic analysis, we report that several species of the biomining acidophilic bacteria contain orthologs of CHAD proteins with high sequence identity. Furthermore, the gene coding for the CHAD protein is in the same genetic context of the enzyme polyphosphate kinase (PPK), which is in charge of PolyP synthesis. Particularly, the group of ppk and CHAD genes is highly conserved. Metallosphaera sedula and other acidophilic archaea used in biomining also contain CHAD proteins. These archaea show high levels of identity in genes coding for a cluster having the same organization. Amongst these genes are chad and ppx. In general, both biomining bacteria and archaea contain high PolyP levels and are highly resistant to heavy metals. Therefore, the presence of this conserved genetic organization suggests a high relevance for their metabolism. It has been formerly reported that a crystallized CHAD protein contains a copper-binding site. Based on this previous knowledge, in the present report, it was determined that all analyzed CHAD proteins are very conserved at their structural level. In addition, it was found that the lack of YgiF, an Escherichia coli CHAD-containing protein, decreases copper resistance in this bacterium. This phenotype was not only complemented by transforming E. coli with YgiF but also by expressing CHAD from Acidithiobacillus ferrooxidans in it. Interestingly, the strains in which the possible copper-binding sites were mutated were also more metal sensitive. Based on these results, we propose that CHAD proteins are involved in copper resistance in microorganisms. These findings are very interesting and may eventually improve biomining operations in the future.

An ANI gap within bacterial species that advances the definitions of intra-species units.

IMPORTANCE: Bacterial strains and clonal complexes are two cornerstone concepts for microbiology that remain loosely defined, which confuses communication and research. Here we identify a natural gap in genome sequence comparisons among isolate genomes of all well-sequenced species that has gone unnoticed so far and could be used to more accurately and precisely define these and related concepts compared to current methods. These findings advance the molecular toolbox for accurately delineating and following the important units of diversity within prokaryotic species and thus should greatly facilitate future epidemiological and micro-diversity studies across clinical and environmental settings.

Diversity, Taxonomic Novelty, and Encoded Functions of Salar de Ascotán Microbiota, as Revealed by Metagenome-Assembled Genomes.

Salar de Ascotán is a high-altitude arsenic-rich salt flat exposed to high ultraviolet radiation in the Atacama Desert, Chile. It hosts unique endemic flora and fauna and is an essential habitat for migratory birds, making it an important site for conservation and protection. However, there is limited information on the resident microbiota's diversity, genomic features, metabolic potential, and molecular mechanisms that enable it to thrive in this extreme environment. We used long- and short-read metagenomics to investigate the microbial communities in Ascotán's water, sediment, and soil. Bacteria predominated, mainly Pseudomonadota, Acidobacteriota, and Bacteroidota, with a remarkable diversity of archaea in the soil. Following hybrid assembly, we recovered high-quality bacterial (101) and archaeal (6) metagenome-assembled genomes (MAGs), including representatives of two putative novel families of Patescibacteria and Pseudomonadota and two novel orders from the archaeal classes Halobacteriota and Thermoplasmata. We found different metabolic capabilities across distinct lineages and a widespread presence of genes related to stress response, DNA repair, and resistance to arsenic and other metals. These results highlight the remarkable diversity and taxonomic novelty of the Salar de Ascotán microbiota and its rich functional repertoire, making it able to resist different harsh conditions. The highly complete MAGs described here could serve future studies and bioprospection efforts focused on salt flat extremophiles, and contribute to enriching databases with microbial genome data from underrepresented regions of our planet.

In situ cell division and mortality rates of SAR11, SAR86, Bacteroidetes, and Aurantivirga during phytoplankton blooms reveal differences in population controls.

Net growth of microbial populations, that is, changes in abundances over time, can be studied using 16S rRNA fluorescence in situ hybridization (FISH). However, this approach does not differentiate between mortality and cell division rates. We used FISH-based image cytometry in combination with dilution culture experiments to study net growth, cell division, and mortality rates of four bacterial taxa over two distinct phytoplankton blooms: the oligotrophs SAR11 and SAR86, and the copiotrophic phylum Bacteroidetes, and its genus Aurantivirga. Cell volumes, ribosome content, and frequency of dividing cells (FDC) co-varied over time. Among the three, FDC was the most suitable predictor to calculate cell division rates for the selected taxa. The FDC-derived cell division rates for SAR86 of up to 0.8/day and Aurantivirga of up to 1.9/day differed, as expected for oligotrophs and copiotrophs. Surprisingly, SAR11 also reached high cell division rates of up to 1.9/day, even before the onset of phytoplankton blooms. For all four taxonomic groups, the abundance-derived net growth (-0.6 to 0.5/day) was about an order of magnitude lower than the cell division rates. Consequently, mortality rates were comparably high to cell division rates, indicating that about 90% of bacterial production is recycled without apparent time lag within 1 day. Our study shows that determining taxon-specific cell division rates complements omics-based tools and provides unprecedented clues on individual bacterial growth strategies including bottom-up and top-down controls. IMPORTANCE The growth of a microbial population is often calculated from their numerical abundance over time. However, this does not take cell division and mortality rates into account, which are important for deriving ecological processes like bottom-up and top-down control. In this study, we determined growth by numerical abundance and calibrated microscopy-based methods to determine the frequency of dividing cells and subsequently calculate taxon-specific cell division rates in situ. The cell division and mortality rates of two oligotrophic (SAR11 and SAR86) and two copiotrophic (Bacteroidetes and Aurantivirga) taxa during two spring phytoplankton blooms showed a tight coupling for all four taxa throughout the blooms without any temporal offset. Unexpectedly, SAR11 showed high cell division rates days before the bloom while cell abundances remained constant, which is indicative of strong top-down control. Microscopy remains the method of choice to understand ecological processes like top-down and bottom-up control on a cellular level.

Comparing genomes recovered from time-series metagenomes using long- and short-read sequencing technologies.

BACKGROUND: Over the past years, sequencing technologies have expanded our ability to examine novel microbial metabolisms and diversity previously obscured by isolation approaches. Long-read sequencing promises to revolutionize the metagenomic field and recover less fragmented genomes from environmental samples. Nonetheless, how to best benefit from long-read sequencing and whether long-read sequencing can provide recovered genomes of similar characteristics as short-read approaches remains unclear. RESULTS: We recovered metagenome-assembled genomes (MAGs) from the free-living fraction at four-time points during a spring bloom in the North Sea. The taxonomic composition of all MAGs recovered was comparable between technologies. However, differences consisted of higher sequencing depth for contigs and higher genome population diversity in short-read compared to long-read metagenomes. When pairing population genomes recovered from both sequencing approaches that shared ≥ 99% average nucleotide identity, long-read MAGs were composed of fewer contigs, a higher N50, and a higher number of predicted genes when compared to short-read MAGs. Moreover, 88% of the total long-read MAGs carried a 16S rRNA gene compared to only 23% of MAGs recovered from short-read metagenomes. Relative abundances for population genomes recovered using both technologies were similar, although disagreements were observed for high and low GC content MAGs. CONCLUSIONS: Our results highlight that short-read technologies recovered more MAGs and a higher number of species than long-read due to an overall higher sequencing depth. Long-read samples produced higher quality MAGs and similar species composition compared to short-read sequencing. Differences in the GC content recovered by each sequencing technology resulted in divergences in the diversity recovered and relative abundance of MAGs within the GC content boundaries.

ROCker Models for Reliable Detection and Typing of Short-Read Sequences Carrying ß-Lactamase Genes.

Identification of genes encoding ß-lactamases (BLs) from short-read sequences remains challenging due to the high frequency of shared amino acid functional domains and motifs in proteins encoded by BL genes and related non-BL gene sequences. Divergent BL homologs can be frequently missed during similarity searches, which has important practical consequences for monitoring antibiotic resistance. To address this limitation, we built ROCker models that targeted broad classes (e.g., class A, B, C, and D) and individual families (e.g., TEM) of BLs and challenged them with mock 150-bp- and 250-bp-read data sets of known composition. ROCker identifies most-discriminant bit score thresholds in sliding windows along the sequence of the target protein sequence and hence can account for nondiscriminative domains shared by unrelated proteins. BL ROCker models showed a 0% false-positive rate (FPR), a 0% to 4% false-negative rate (FNR), and an up-to-50-fold-higher F1 score [2 × precision × recall/(precision + recall)] compared to alternative methods, such as similarity searches using BLASTx with various e-value thresholds and BL hidden Markov models, or tools like DeepARG, ShortBRED, and AMRFinder. The ROCker models and the underlying protein sequence reference data sets and phylogenetic trees for read placement are freely available through http://enve-omics.ce.gatech.edu/data/rocker-bla. Application of these BL ROCker models to metagenomics, metatranscriptomics, and high-throughput PCR gene amplicon data should facilitate the reliable detection and quantification of BL variants encoded by environmental or clinical isolates and microbiomes and more accurate assessment of the associated public health risk, compared to the current practice. IMPORTANCE Resistance genes encoding ß-lactamases (BLs) confer resistance to the widely prescribed antibiotic class ß-lactams. Therefore, it is important to assess the prevalence of BL genes in clinical or environmental samples for monitoring the spreading of these genes into pathogens and estimating public health risk. However, detecting BLs in short-read sequence data is technically challenging. Our ROCker model-based bioinformatics approach showcases the reliable detection and typing of BLs in complex data sets and thus contributes toward solving an important problem in antibiotic resistance surveillance. The ROCker models developed substantially expand the toolbox for monitoring antibiotic resistance in clinical or environmental settings.

Verrucomicrobiota are specialist consumers of sulfated methyl pentoses during diatom blooms.

Marine algae annually sequester petagrams of carbon dioxide into polysaccharides, which are a central metabolic fuel for marine carbon cycling. Diatom microalgae produce sulfated polysaccharides containing methyl pentoses that are challenging to degrade for bacteria compared to other monomers, implicating these sugars as a potential carbon sink. Free-living bacteria occurring in phytoplankton blooms that specialise on consuming microalgal sugars, containing fucose and rhamnose remain unknown. Here, genomic and proteomic data indicate that small, coccoid, free-living Verrucomicrobiota specialise in fucose and rhamnose consumption during spring algal blooms in the North Sea. Verrucomicrobiota cell abundance was coupled with the algae bloom onset and accounted for up to 8% of the bacterioplankton. Glycoside hydrolases, sulfatases, and bacterial microcompartments, critical proteins for the consumption of fucosylated and sulfated polysaccharides, were actively expressed during consecutive spring bloom events. These specialised pathways were assigned to novel and discrete candidate species of the Akkermansiaceae and Puniceicoccaceae families, which we here describe as Candidatus Mariakkermansia forsetii and Candidatus Fucivorax forsetii. Moreover, our results suggest specialised metabolic pathways could determine the fate of complex polysaccharides consumed during algae blooms. Thus the sequestration of phytoplankton organic matter via methyl pentose sugars likely depend on the activity of specialised Verrucomicrobiota populations.

Transcriptomic and rRNA:rDNA Signatures of Environmental versus Enteric Enterococcus faecalis Isolates under Oligotrophic Freshwater Conditions.

The use of enterococci as a fecal indicator bacterial group for public health risk assessment has been brought into question by recent studies showing that "naturalized" populations of Enterococcus faecalis exist in the extraenteric environment. The extent to which these naturalized E. faecalis organisms can confound water quality monitoring is unclear. To determine if strains isolated from different habitats display different survival strategies and responses, we compared the decay patterns of three E. faecalis isolates from the natural environment (environmental strains) against three human gut isolates (enteric strains) in laboratory mesocosms that simulate an oligotrophic, aerobic freshwater environment. Our results showed similar overall decay rates between enteric and environmental isolates based on viable plate and quantitative PCR (qPCR) counts. However, the enteric isolates exhibited a spike in copy number ratios of 16S rRNA gene transcripts to 16S rRNA gene DNA copies (rRNA:rDNA ratios) between days 1 and 3 of the mesocosm incubations that was not observed in environmental isolates, which could indicate a different stress response. Nevertheless, there was no strong evidence of differential gene expression between environmental and enteric isolates related to habitat adaptation in the accompanying mesocosm metatranscriptomes. Overall, our results provide novel information on how rRNA levels may vary over different growth conditions (e.g., standard lab versus oligotrophic) for this important indicator bacteria. We also observed some evidence for habitat adaptation in E. faecalis; however, this adaptation may not be substantial or consistent enough for integration in water quality monitoring. IMPORTANCE Enterococci are commonly used worldwide to monitor environmental fecal contamination and public health risk for waterborne diseases. However, closely related enterococci strains adapted to living in the extraenteric environment may represent a lower public health risk and confound water quality estimates. We developed an rRNA:rDNA viability assay for E. faecalis (a predominant species within this fecal group) and tested it against both enteric and environmental isolates in freshwater mesocosms to assess whether this approach can serve as a more sensitive water quality monitoring tool. We were unable to reliably distinguish the different isolate types using this assay under the conditions tested; thus, environmental strains should continue to be counted during routine water monitoring. However, this assay could be useful for distinguishing more recent (i.e., higher-risk) fecal pollution because rRNA levels significantly decreased after 1 week in all isolates.

The influence of alfalfa-switchgrass intercropping on microbial community structure and function.

The use of nitrogen fertilizer on bioenergy crops such as switchgrass results in increased costs, nitrogen leaching and emissions of N2 O, a potent greenhouse gas. Intercropping with nitrogen-fixing alfalfa has been proposed as an environmentally sustainable alternative, but the effects of synthetic fertilizer versus intercropping on soil microbial community functionality remain uncharacterized. We analysed 24 metagenomes from the upper soil layer of agricultural fields from Prosser, WA over two growing seasons and representing three agricultural practices: unfertilized switchgrass (control), fertilized switchgrass and switchgrass intercropped with alfalfa. The synthetic fertilization and intercropping did not result in major shifts of microbial community taxonomic and functional composition compared with the control plots, but a few significant changes were noted. Most notably, mycorrhizal fungi, ammonia-oxidizing archaea and bacteria increased in abundance with intercropping and fertilization. However, only betaproteobacterial ammonia-oxidizing bacteria abundance in fertilized plots significantly correlated to N2 O emission and companion qPCR data. Collectively, a short period of intercropping elicits minor but significant changes in the soil microbial community toward nitrogen preservation and that intercropping may be a viable alternative to synthetic fertilization.

Microbial metagenome-assembled genomes of the Fram Strait from short and long read sequencing platforms.

The impacts of climate change on the Arctic Ocean are manifesting throughout the ecosystem at an unprecedented rate. Of global importance are the impacts on heat and freshwater exchange between the Arctic and North Atlantic Oceans. An expanding Atlantic influence in the Arctic has accelerated sea-ice decline, weakened water column stability and supported the northward shift of temperate species. The only deep-water gateway connecting the Arctic and North Atlantic and thus, fundamental for these exchange processes is the Fram Strait. Previous research in this region is extensive, however, data on the ecology of microbial communities is limited, reflecting the wider bias towards temperate and tropical latitudes. Therefore, we present 14 metagenomes, 11 short-read from Illumina and three long-read from PacBio Sequel II, of the 0.2-3 µm fraction to help alleviate such biases and support future analyses on changing ecological patterns. Additionally, we provide 136 species-representative, manually refined metagenome-assembled genomes which can be used for comparative genomics analyses and addressing questions regarding functionality or distribution of taxa.

Diversity of microbial communities and genes involved in nitrous oxide emissions in Antarctic soils impacted by marine animals as revealed by metagenomics and 100 metagenome-assembled genomes.

Antarctic soils generally have low temperatures and limited availability of liquid water and nutrients. However, animals can increase the nutrient availability of ice-free areas by transferring nutrients from marine to terrestrial ecosystems, mainly through their excreta. In this study, we employed shotgun metagenomics and population genome binning techniques to study the diversity of microbial communities in Antarctic soils impacted by marine pinnipeds and birds relative to soils with no evident animal presence. We obtained ~285,000 16S rRNA gene-carrying metagenomic reads representing ~60 phyla and 100 metagenome-assembled genomes (MAGs) representing eight phyla. Only nine of these 100 MAGs represented previously described species, revealing that these soils harbor extensive novel diversity. Proteobacteria, Actinobacteria, and Bacteroidetes were the most abundant phyla in all samples, with Rhodanobacter being one of the most abundant genera in the bird-impacted soils. Further, the relative abundance of genes related to denitrification was at least double in soils impacted by birds than soils without animal influence. These results advance our understanding of the microbial populations and their genes involved in nitrous oxide emissions in ice-free coastal Antarctic soils impacted by marine animals and reveal novel microbial diversity associated with these ecosystems.

Metagenomic Characterization of Soil Microbial Communities in the Luquillo Experimental Forest (Puerto Rico) and Implications for Nitrogen Cycling.

The phylogenetic and functional diversities of microbial communities in tropical rainforests and how these differ from those of temperate communities remain poorly described but are directly related to the increased fluxes of greenhouse gases such as nitrous oxide (N2O) from the tropics. Toward closing these knowledge gaps, we analyzed replicated shotgun metagenomes representing distinct life zones and an elevation gradient from four locations in the Luquillo Experimental Forest (LEF), Puerto Rico. These soils had a distinct microbial community composition and lower species diversity compared to those of temperate grasslands or agricultural soils. In contrast to the overall distinct community composition, the relative abundances and nucleotide sequences of N2O reductases (nosZ) were highly similar between tropical forest and temperate soils. However, respiratory NO reductase (norB) was 2-fold more abundant in the tropical soils, which might be relatable to their greater N2O emissions. Nitrogen fixation (nifH) also showed higher relative abundance in rainforest than in temperate soils, i.e., 20% versus 0.1 to 0.3% of bacterial genomes in each soil type harbored the gene, respectively. Finally, unlike temperate soils, LEF soils showed little stratification with depth in the first 0 to 30 cm, with ∼45% of community composition differences explained solely by location. Collectively, these results advance our understanding of spatial diversity and metabolic repertoire of tropical rainforest soil communities and should facilitate future ecological studies of these ecosystems. IMPORTANCE Tropical rainforests are the largest terrestrial sinks of atmospheric CO2 and the largest natural source of N2O emissions, two greenhouse gases that are critical for the climate. The microbial communities of rainforest soils that directly or indirectly, through affecting plant growth, contribute to these fluxes remain poorly described by cultured-independent methods. To close this knowledge gap, the present study applied shotgun metagenomics to samples selected from three distinct life zones within the Puerto Rico rainforest. The results advance our understanding of microbial community diversity in rainforest soils and should facilitate future studies of natural or manipulated perturbations of these critical ecosystems.

Distinct ecotypes within a natural haloarchaeal population enable adaptation to changing environmental conditions without causing population sweeps.

Microbial communities thriving in hypersaline brines of solar salterns are highly resistant and resilient to environmental changes, and salinity is a major factor that deterministically influences community structure. Here, we demonstrate that this resilience occurs even after rapid osmotic shocks caused by a threefold change in salinity (a reduction from 34 to 12% salts) leading to massive amounts of archaeal cell lysis. Specifically, our temporal metagenomic datasets identified two co-occurring ecotypes within the most dominant archaeal population of the brines Haloquadratum walsbyi that exhibited different salt concentration preferences. The dominant ecotype was generally more abundant and occurred in high-salt conditions (34%); the low abundance ecotype always co-occurred but was enriched at salinities around 20% or lower and carried unique gene content related to solute transport and gene regulation. Despite their apparent distinct ecological preferences, the ecotypes did not outcompete each other presumably due to weak functional differentiation between them. Further, the osmotic shock selected for a temporal increase in taxonomic and functional diversity at both the Hqr. walsbyi population and whole-community levels supporting the specialization-disturbance hypothesis, that is, the expectation that disturbance favors generalists. Altogether, our results provide new insights into how intraspecies diversity is maintained in light of substantial gene-content differences and major environmental perturbations.

Candidatus Abditibacter, a novel genus within the Cryomorphaceae, thriving in the North Sea.

Coastal phytoplankton blooms are frequently followed by successive blooms of heterotrophic bacterial clades. The class Flavobacteriia within the Bacteroidetes has been shown to play an important role in the degradation of high molecular weight substrates that become available in the later stages of such blooms. One of the flavobacterial clades repeatedly observed over the course of several years during phytoplankton blooms off the coast of Helgoland, North Sea, is Vis6. This genus-level clade belongs to the family Cryomorphaceae and has been resistant to cultivation to date. Based on metagenome assembled genomes, comparative 16S rRNA gene sequence analyses and fluorescence in situ hybridization, we here propose a novel candidate genus Abditibacter, comprising three novel species Candidatus Abditibacter vernus, Candidatus Abditibacter forsetii and Candidatus Abditibacter autumni. While the small genomes of the three novel photoheterotrophic species encode highly similar gene repertoires, including genes for degradation of proteins and algal storage polysaccharides such as laminarin, two of them - Ca. A. vernus and Ca. A. forsetii - seem to have a preference for spring blooms, while Ca. A. autumni almost exclusively occurs in late summer and autumn.

Metagenomics as a Public Health Risk Assessment Tool in a Study of Natural Creek Sediments Influenced by Agricultural and Livestock Runoff: Potential and Limitations.

Little is known about the public health risks associated with natural creek sediments that are affected by runoff and fecal pollution from agricultural and livestock practices. For instance, the persistence of foodborne pathogens such as Shiga toxin-producing Escherichia coli (STEC) originating from these practices remains poorly quantified. Towards closing these knowledge gaps, the water-sediment interface of two creeks in the Salinas River Valley of California was sampled over a 9-month period using metagenomics and traditional culture-based tests for STEC. Our results revealed that these sediment communities are extremely diverse and have functional and taxonomic diversity comparable to that observed in soils. With our sequencing effort (∼4 Gbp per library), we were unable to detect any pathogenic E. coli in the metagenomes of 11 samples that had tested positive using culture-based methods, apparently due to relatively low abundance. Furthermore, there were no significant differences in the abundance of human- or cow-specific gut microbiome sequences in the downstream impacted sites compared to that in upstream more pristine (control) sites, indicating natural dilution of anthropogenic inputs. Notably, the high number of metagenomic reads carrying antibiotic resistance genes (ARGs) found in all samples was significantly higher than ARG reads in other available freshwater and soil metagenomes, suggesting that these communities may be natural reservoirs of ARGs. The work presented here should serve as a guide for sampling volumes, amount of sequencing to apply, and what bioinformatics analyses to perform when using metagenomics for public health risk studies of environmental samples such as sediments.IMPORTANCE Current agricultural and livestock practices contribute to fecal contamination in the environment and the spread of food- and waterborne disease and antibiotic resistance genes (ARGs). Traditionally, the level of pollution and risk to public health are assessed by culture-based tests for the intestinal bacterium Escherichia coli However, the accuracy of these traditional methods (e.g., low accuracy in quantification, and false-positive signal when PCR based) and their suitability for sediments remain unclear. We collected sediments for a time series metagenomics study from one of the most highly productive agricultural regions in the United States in order to assess how agricultural runoff affects the native microbial communities and if the presence of Shiga toxin-producing Escherichia coli (STEC) in sediment samples can be detected directly by sequencing. Our study provided important information on the potential for using metagenomics as a tool for assessment of public health risk in natural environments.

Comparing DNA, RNA and protein levels for measuring microbial dynamics in soil microcosms amended with nitrogen fertilizer.

To what extent multi-omic techniques could reflect in situ microbial process rates remains unclear, especially for highly diverse habitats like soils. Here, we performed microcosm incubations using sandy soil from an agricultural site in Midwest USA. Microcosms amended with isotopically labeled ammonium and urea to simulate a fertilization event showed nitrification (up to 4.1 ± 0.87 µg N-NO3- g-1 dry soil d-1) and accumulation of N2O after 192 hours of incubation. Nitrification activity (NH4+ → NH2OH → NO → NO2- → NO3-) was accompanied by a 6-fold increase in relative expression of the 16S rRNA gene (RNA/DNA) between 10 and 192 hours of incubation for ammonia-oxidizing bacteria Nitrosomonas and Nitrosospira, unlike archaea and comammox bacteria, which showed stable gene expression. A strong relationship between nitrification activity and betaproteobacterial ammonia monooxygenase and nitrite oxidoreductase transcript abundances revealed that mRNA quantitatively reflected measured activity and was generally more sensitive than DNA under these conditions. Although peptides related to housekeeping proteins from nitrite-oxidizing microorganisms were detected, their abundance was not significantly correlated with activity, revealing that meta-proteomics provided only a qualitative assessment of activity. Altogether, these findings underscore the strengths and limitations of multi-omic approaches for assessing diverse microbial communities in soils and provide new insights into nitrification.

Niche differentiation among annually recurrent coastal Marine Group II Euryarchaeota.

Since the discovery of archaeoplankton in 1992, the euryarchaeotal Marine Group II (MGII) remains uncultured and less understood than other planktonic archaea. We characterized the seasonal dynamics of MGII populations in the southern North Sea on a genomic and microscopic level over the course of four years. We recovered 34 metagenome-assembled genomes (MAGs) of MGIIa and MGIIb that corroborated proteorhodopsin-based photoheterotrophic lifestyles. However, MGIIa and MGIIb MAG genome sizes differed considerably (~1.9 vs. ~1.4 Mbp), as did their transporter, peptidase, flagella and sulfate assimilation gene repertoires. MGIIb populations were characteristic of winter samples, whereas MGIIa accounted for up to 23% of the community at the beginning of summer. Both clades consisted of annually recurring, sequence-discrete populations with low intra-population sequence diversity. Oligotyping of filtered cell-size fractions and microscopy consistently suggested that MGII cells were predominantly free-living. Cells were coccoid and ~0.7 µm in diameter, likely resulting in grazing avoidance. Based on multiple lines of evidence, we propose distinct niche adaptations of MGIIa and MGIIb Euryarchaeota populations that are characteristic of summer and winter conditions in the coastal North Sea.