Microcoleus (Cyanobacteria) form watershed-wide populations without strong gradients in population structure.

The relative importance of separation by distance and by environment to population genetic diversity can be conveniently tested in river networks, where these two drivers are often independently distributed over space. To evaluate the importance of dispersal and environmental conditions in shaping microbial population structures, we performed genome-resolved metagenomic analyses of benthic Microcoleus-dominated cyanobacterial mats collected in the Eel and Russian River networks (California, USA). The 64 Microcoleus genomes were clustered into three species that shared >96.5% average nucleotide identity (ANI). Most mats were dominated by one strain, but minor alleles within mats were often shared, even over large spatial distances (>300 km). Within the most common Microcoleus species, the ANI between the dominant strains within mats decreased with increasing spatial separation. However, over shorter spatial distances (tens of kilometres), mats from different subwatersheds had lower ANI than mats from the same subwatershed, suggesting that at shorter spatial distances environmental differences between subwatersheds in factors like canopy cover, conductivity, and mean annual temperature decreases ANI. Since mats in smaller creeks had similar levels of nucleotide diversity (π) as mats in larger downstream subwatersheds, within-mat genetic diversity does not appear to depend on the downstream accumulation of upstream-derived strains. The four-gamete test and sequence length bias suggest recombination occurs between almost all strains within each species, even between populations separated by large distances or living in different habitats. Overall, our results show that, despite some isolation by distance and environmental conditions, sufficient gene-flow occurs among cyanobacterial strains to prevent either driver from producing distinctive population structures across the watershed.

Genome Streamlining, Plasticity, and Metabolic Versatility Distinguish Co-occurring Toxic and Nontoxic Cyanobacterial Strains of Microcoleus.

Harmful cyanobacterial bloom occurrences have increased worldwide due to climate change and eutrophication, causing nuisance and animal deaths. Species from the benthic cyanobacterial genus Microcoleus are ubiquitous and form thick mats in freshwater systems, such as rivers, that are sometimes toxic due to the production of potent neurotoxins (anatoxins). Anatoxin-producing (toxic) strains typically coexist with non-anatoxin-producing (nontoxic) strains in mats, although the reason for this is unclear. To determine the genetic mechanisms differentiating toxic and nontoxic Microcoleus, we sequenced and assembled genomes from 11 cultures and compared these to another 31 Microcoleus genomes. Average nucleotide identities (ANI) indicate that toxic and nontoxic strains are distinct species (ANI, <95%), and only 6% of genes are shared across all 42 genomes, suggesting a high level of genetic divergence among Microcoleus strains. Comparative genomics showed substantial genome streamlining in toxic strains and a potential dependency on external sources for thiamine and sucrose. Toxic and nontoxic strains are further differentiated by an additional set of putative nitrate transporter (nitrogen uptake) and cyanophycin (carbon and nitrogen storage) genes, respectively. These genes likely confer distinct competitive advantages based on nutrient availability and suggest nontoxic strains are more robust to nutrient fluctuations. Nontoxic strains also possess twice as many transposable elements, potentially facilitating greater genetic adaptation to environmental changes. Our results offer insights into the divergent evolution of Microcoleus strains and the potential for cooperative and competitive interactions that contribute to the co-occurrence of toxic and nontoxic species within mats. IMPORTANCE Microcoleus autumnalis, and closely related Microcoleus species, compose a geographically widespread group of freshwater benthic cyanobacteria. Canine deaths due to anatoxin-a poisoning, following exposure to toxic proliferations, have been reported globally. While Microcoleus proliferations are on the rise, the mechanisms underpinning competition between, or coexistence of, toxic and nontoxic strains are unknown. This study identifies substantial genetic differences between anatoxin-producing and non-anatoxin-producing strains, pointing to reduced metabolic flexibility in toxic strains, and potential dependence on cohabiting nontoxic strains. Results provide insights into the metabolic and evolutionary differences between toxic and nontoxic Microcoleus, which may assist in predicting and managing aquatic proliferations.

Blue Waters, Green Bottoms: Benthic Filamentous Algal Blooms Are an Emerging Threat to Clear Lakes Worldwide.

Nearshore (littoral) habitats of clear lakes with high water quality are increasingly experiencing unexplained proliferations of filamentous algae that grow on submerged surfaces. These filamentous algal blooms (FABs) are sometimes associated with nutrient pollution in groundwater, but complex changes in climate, nutrient transport, lake hydrodynamics, and food web structure may also facilitate this emerging threat to clear lakes. A coordinated effort among members of the public, managers, and scientists is needed to document the occurrence of FABs, to standardize methods for measuring their severity, to adapt existing data collection networks to include nearshore habitats, and to mitigate and reverse this profound structural change in lake ecosystems. Current models of lake eutrophication do not explain this littoral greening. However, a cohesive response to it is essential for protecting some of the world's most valued lakes and the flora, fauna, and ecosystem services they sustain.

Patterns of Gene Content and Co-occurrence Constrain the Evolutionary Path toward Animal Association in Candidate Phyla Radiation Bacteria.

Candidate Phyla Radiation (CPR) bacteria are small, likely episymbiotic organisms found across Earth's ecosystems. Despite their prevalence, the distribution of CPR lineages across habitats and the genomic signatures of transitions among these habitats remain unclear. Here, we expand the genome inventory for Absconditabacteria (SR1), Gracilibacteria, and Saccharibacteria (TM7), CPR bacteria known to occur in both animal-associated and environmental microbiomes, and investigate variation in gene content with habitat of origin. By overlaying phylogeny with habitat information, we show that bacteria from these three lineages have undergone multiple transitions from environmental habitats into animal microbiomes. Based on co-occurrence analyses of hundreds of metagenomes, we extend the prior suggestion that certain Saccharibacteria have broad bacterial host ranges and constrain possible host relationships for Absconditabacteria and Gracilibacteria. Full-proteome analyses show that animal-associated Saccharibacteria have smaller gene repertoires than their environmental counterparts and are enriched in numerous protein families, including those likely functioning in amino acid metabolism, phage defense, and detoxification of peroxide. In contrast, some freshwater Saccharibacteria encode a putative rhodopsin. For protein families exhibiting the clearest patterns of differential habitat distribution, we compared protein and species phylogenies to estimate the incidence of lateral gene transfer and genomic loss occurring over the species tree. These analyses suggest that habitat transitions were likely not accompanied by large transfer or loss events but rather were associated with continuous proteome remodeling. Thus, we speculate that CPR habitat transitions were driven largely by availability of suitable host taxa and were reinforced by acquisition and loss of some capacities. IMPORTANCE Studying the genetic differences between related microorganisms from different environment types can indicate factors associated with their movement among habitats. This is particularly interesting for bacteria from the Candidate Phyla Radiation because their minimal metabolic capabilities require associations with microbial hosts. We found that shifts of Absconditabacteria, Gracilibacteria, and Saccharibacteria between environmental ecosystems and mammalian mouths/guts probably did not involve major episodes of gene gain and loss; rather, gradual genomic change likely followed habitat migration. The results inform our understanding of how little-known microorganisms establish in the human microbiota where they may ultimately impact health.

inStrain profiles population microdiversity from metagenomic data and sensitively detects shared microbial strains.

Coexisting microbial cells of the same species often exhibit genetic variation that can affect phenotypes ranging from nutrient preference to pathogenicity. Here we present inStrain, a program that uses metagenomic paired reads to profile intra-population genetic diversity (microdiversity) across whole genomes and compares microbial populations in a microdiversity-aware manner, greatly increasing the accuracy of genomic comparisons when benchmarked against existing methods. We use inStrain to profile >1,000 fecal metagenomes from newborn premature infants and find that siblings share significantly more strains than unrelated infants, although identical twins share no more strains than fraternal siblings. Infants born by cesarean section harbor Klebsiella with significantly higher nucleotide diversity than infants delivered vaginally, potentially reflecting acquisition from hospital rather than maternal microbiomes. Genomic loci that show diversity in individual infants include variants found between other infants, possibly reflecting inoculation from diverse hospital-associated sources. inStrain can be applied to any metagenomic dataset for microdiversity analysis and rigorous strain comparison.

Soil bacterial populations are shaped by recombination and gene-specific selection across a grassland meadow.

Soil microbial diversity is often studied from the perspective of community composition, but less is known about genetic heterogeneity within species. The relative impacts of clonal interference, gene-specific selection, and recombination in many abundant but rarely cultivated soil microbes remain unknown. Here we track genome-wide population genetic variation for 19 highly abundant bacterial species sampled from across a grassland meadow. Genomic inferences about population structure are made using the millions of sequencing reads that are assembled de novo into consensus genomes from metagenomes, as each read pair describes a short genomic sequence from a cell in each population. Genomic nucleotide identity of assembled genomes was significantly associated with local geography for over half of the populations studied, and for a majority of populations within-sample nucleotide diversity could often be as high as meadow-wide nucleotide diversity. Genes involved in metabolite biosynthesis and extracellular transport were characterized by elevated nucleotide diversity in multiple species. Microbial populations displayed varying degrees of homologous recombination and recombinant variants were often detected at 7-36% of loci genome-wide. Within multiple populations we identified genes with unusually high spatial differentiation of alleles, fewer recombinant events, elevated ratios of nonsynonymous to synonymous variants, and lower nucleotide diversity, suggesting recent selective sweeps for gene variants. Taken together, these results indicate that recombination and gene-specific selection commonly shape genetic variation in several understudied soil bacterial lineages.

Clades of huge phages from across Earth's ecosystems.

Bacteriophages typically have small genomes1 and depend on their bacterial hosts for replication2. Here we sequenced DNA from diverse ecosystems and found hundreds of phage genomes with lengths of more than 200 kilobases (kb), including a genome of 735 kb, which is-to our knowledge-the largest phage genome to be described to date. Thirty-five genomes were manually curated to completion (circular and no gaps). Expanded genetic repertoires include diverse and previously undescribed CRISPR-Cas systems, transfer RNAs (tRNAs), tRNA synthetases, tRNA-modification enzymes, translation-initiation and elongation factors, and ribosomal proteins. The CRISPR-Cas systems of phages have the capacity to silence host transcription factors and translational genes, potentially as part of a larger interaction network that intercepts translation to redirect biosynthesis to phage-encoded functions. In addition, some phages may repurpose bacterial CRISPR-Cas systems to eliminate competing phages. We phylogenetically define the major clades of huge phages from human and other animal microbiomes, as well as from oceans, lakes, sediments, soils and the built environment. We conclude that the large gene inventories of huge phages reflect a conserved biological strategy, and that the phages are distributed across a broad bacterial host range and across Earth's ecosystems.

Toxic benthic freshwater cyanobacterial proliferations: Challenges and solutions for enhancing knowledge and improving monitoring and mitigation.

1. This review summarises knowledge on the ecology, toxin production, and impacts of toxic freshwater benthic cyanobacterial proliferations. It documents monitoring, management, and sampling strategies, and explores mitigation options. 2. Toxic proliferations of freshwater benthic cyanobacteria (taxa that grow attached to substrates) occur in streams, rivers, lakes, and thermal and meltwater ponds, and have been reported in 19 countries. Anatoxin- and microcystin-containing mats are most commonly reported (eight and 10 countries, respectively). 3. Studies exploring factors that promote toxic benthic cyanobacterial proliferations are limited to a few species and habitats. There is a hierarchy of importance in environmental and biological factors that regulate proliferations with variables such as flow (rivers), fine sediment deposition, nutrients, associated microbes, and grazing identified as key drivers. Regulating factors differ among colonisation, expansion, and dispersal phases. 4. New -omics-based approaches are providing novel insights into the physiological attributes of benthic cyanobacteria and the role of associated microorganisms in facilitating their proliferation. 5. Proliferations are commonly comprised of both toxic and non-toxic strains, and the relative proportion of these is the key factor contributing to the overall toxin content of each mat. 6. While these events are becoming more commonly reported globally, we currently lack standardised approaches to detect, monitor, and manage this emerging health issue. To solve these critical gaps, global collaborations are needed to facilitate the rapid transfer of knowledge and promote the development of standardised techniques that can be applied to diverse habitats and species, and ultimately lead to improved management.

Multiple cyanotoxin congeners produced by sub-dominant cyanobacterial taxa in riverine cyanobacterial and algal mats.

Benthic cyanobacterial proliferations in rivers are have been reported with increasing frequency worldwide. In the Eel and Russian rivers of California, more than a dozen dog deaths have been attributed to cyanotoxin toxicosis since 2000. Periphyton proliferations in these rivers comprise multiple cyanobacterial taxa capable of cyanotoxin production, hence there is uncertainty regarding which taxa are producing toxins. In this study, periphyton samples dominated by the cyanobacterial genera Anabaena spp. and Microcoleus spp. and the green alga Cladophora glomerata were collected from four sites in the Eel River catchment and one site in the Russian River. Samples were analysed for potential cyanotoxin producers using polymerase chain reaction (PCR) in concert with Sanger sequencing. Cyanotoxin concentrations were measured using liquid chromatography tandem-mass spectrometry, and anatoxin quota (the amount of cyanobacterial anatoxins per toxigenic cell) determined using droplet digital PCR. Sequencing indicated Microcoleus sp. and Nodularia sp. were the putative producers of cyanobacterial anatoxins and nodularins, respectively, regardless of the dominant taxa in the mat. Anatoxin concentrations in the mat samples varied from 0.1 to 18.6 µg g-1 and were significantly different among sites (p < 0.01, Wilcoxon test); however, anatoxin quotas were less variable (< 5-fold). Dihydroanatoxin-a was generally the most abundant variant in samples comprising 38% to 71% of the total anatoxins measured. Mats dominated by the green alga C. glomerata contained both anatoxins and nodularin-R at concentrations similar to those of cyanobacteria-dominated mats. This highlights that even when cyanobacteria are not the dominant taxa in periphyton, these mats may still pose a serious health risk and indicates that more widespread monitoring of all mats in a river are necessary.

Impacts of microbial assemblage and environmental conditions on the distribution of anatoxin-a producing cyanobacteria within a river network.

Blooms of planktonic cyanobacteria have long been of concern in lakes, but more recently, harmful impacts of riverine benthic cyanobacterial mats been recognized. As yet, we know little about how various benthic cyanobacteria are distributed in river networks, or how environmental conditions or other associated microbes in their consortia affect their biosynthetic capacities. We performed metagenomic sequencing for 22 Oscillatoriales-dominated (Cyanobacteria) microbial mats collected across the Eel River network in Northern California and investigated factors associated with anatoxin-a producing cyanobacteria. All microbial communities were dominated by one or two cyanobacterial species, so the key mat metabolisms involve oxygenic photosynthesis and carbon oxidation. Only a few metabolisms fueled the growth of the mat communities, with little evidence for anaerobic metabolic pathways. We genomically defined four cyanobacterial species, all which shared <96% average nucleotide identity with reference Oscillatoriales genomes and are potentially novel species in the genus Microcoleus. One of the Microcoleus species contained the anatoxin-a biosynthesis genes, and we describe the first anatoxin-a gene cluster from the Microcoleus clade within Oscillatoriales. Occurrence of these four Microcoleus species in the watershed was correlated with total dissolved nitrogen and phosphorus concentrations, and the species that contains the anatoxin-a gene cluster was found in sites with higher nitrogen concentrations. Microbial assemblages in mat samples with the anatoxin-a gene cluster consistently had a lower abundance of Burkholderiales (Betaproteobacteria) species than did mats without the anatoxin-producing genes. The associations of water nutrient concentrations and certain co-occurring microbes with anatoxin-a producing Microcoleus motivate further exploration for their roles as potential controls on the distributions of toxigenic benthic cyanobacteria in river networks.

Widespread anatoxin-a detection in benthic cyanobacterial mats throughout a river network.

Benthic algae fuel summer food webs in many sunlit rivers, and are hotspots for primary and secondary production and biogeochemical cycling. Concerningly, riverine benthic algal assemblages can become dominated by toxic cyanobacteria, threatening water quality and public health. In the Eel River in Northern California, over a dozen dog deaths have been attributed to cyanotoxin poisonings since 2000. During the summers of 2013-2015, we documented spatial and temporal patterns of cyanotoxin concentrations in the watershed, showing widespread distribution of anatoxin-a in benthic cyanobacterial mats. Solid phase adsorption toxin tracking (SPATT) samplers were deployed weekly to record dissolved microcystin and anatoxin-a levels at 10 sites throughout the watershed, and 187 Anabaena-dominated or Phormidium-dominated cyanobacterial mat samples were collected from 27 locations to measure intracellular anatoxin-a (ATX) and microcystins (MCY). Anatoxin-a levels were higher than microcystin for both SPATT (mean MCY = 0.8 and ATX = 4.8 ng g resin-1 day-1) and cyanobacterial mat samples (mean MCY = 0.074 and ATX = 1.89 µg g-1 DW). Of the benthic mats sampled, 58.9% had detectable anatoxin-a (max = 70.93 µg g-1 DW), while 37.6% had detectable microcystins (max = 2.29 µg g-1 DW). SPATT cyanotoxin levels peaked in mid-summer in warm mainstem reaches of the watershed. This is one of the first documentations of widespread anatoxin-a occurrence in benthic cyanobacterial mats in a North American watershed.

Rise and fall of toxic benthic freshwater cyanobacteria (Anabaena spp.) in the Eel river: Buoyancy and dispersal.

Benthic cyanobacteria in rivers produce cyanotoxins and affect aquatic food webs, but knowledge of their ecology lags behind planktonic cyanobacteria. The buoyancy of benthic Anabaena spp. mats was studied to understand implications for Anabaena dispersal in the Eel River, California. Field experiments were used to investigate the effects of oxygen bubble production and dissolution on the buoyancy of Anabaena dominated benthic mats in response to light exposure. Samples of Anabaena dominated mats were harvested from the South Fork Eel River and placed in settling columns to measure floating and sinking velocities, or deployed into in situ ambient and low light treatments to measure the effect of light on flotation. Floating and sinking occurred within minutes and were driven by oxygen bubbles produced during photosynthesis, rather than intracellular changes in carbohydrates or gas vesicles. Light experiment results showed that in a natural ambient light regime, mats remained floating for at least 4days, while in low light mats begin to sink in <24h. Floating Anabaena samples were collected from five sites in the watershed and found to contain the cyanotoxins anatoxin-a and microcystin, with higher concentrations of anatoxin-a (median 560, max 30,693ng/gDW) than microcystin (median 30, max 37ng/gDW). The ability of Anabaena mats to maintain their buoyancy will markedly increase their downstream dispersal distances. Increased buoyancy also allows toxin-containing mats to collect along shorelines, increasing threats to human and animal public health.