Structural and functional analysis of the finished genome of the recently isolated toxic Anabaena sp. WA102

Metatranscriptomics reveals gene expression dynamics during an anatoxin-a producing Dolichospermum bloom in a western coastal lake

Cyanobacteria harmful algal blooms in lakes are primarily driven by nutrient and temperature conditions, yet the interplay of these abiotic factors with microbial community dynamics during bloom events is complex and challenging to unravel. Despite advances through deep sequencing approaches, the underlying transcriptomic changes occurring within blooming and non-blooming taxa remains an actively expanding area of study. In this work, we examined a spring-summer bloom event in Anderson Lake, WA, which has experienced recurring annual blooms dominated by the filamentous, anatoxin-a producing, diazotroph: Dolichospermum sp. WA102. Our data reveal the overall transcriptional dominance by Dolichospermum sp. WA102 during the bloom, initiated with increasing temperature and light intensity under high available phosphorus but low nitrogen conditions. We find that heterocyst differentiation was already transcriptionally initiated prior to the bloom, facilitating downstream gene cascades necessary for rapid nitrogen fixation and metabolism. As the bloom progresses, phosphorus becomes depleted, necessitating the expression of Pho regulon components in Dolichospermum sp. WA102 and possibly curtailing the bloom itself. We dissect toxin production and the transcriptional subtleties of the anatoxin-a synthesis locus. Additionally, co-occurring taxa exhibited distinct gene expression profiles, with competition for nutrients, light, and potential allelopathic interactions acting as drivers. Overall, our data provide a unique transcriptomic perspective on a single-taxa-driven, anatoxin-producing bloom, highlighting its competitive adaptation to nutrient acquisition and favorable conditions. This deeper understanding of the genetic mechanisms underlying algal bloom events may aid in predicting and preventing future blooms.

Cuspidothrix Is the First Genetically Proved Anatoxin A Producer in Bulgarian Lakes and Reservoirs

The paper presents the first proof of planktonic cyanoprokaryote genus Cuspidothrix as an anatoxin A (ATX) producer in Bulgarian wetlands. The results from polymerase chain reaction (PCR) obtained from two summer sampling campaigns in 26 selected lakes and reservoirs demonstrated presence of the anaC gene, responsible for ATX production in 21 strains of the genus. They were found in three waterbodies sampled in 2018 (coastal lake Vaya, coastal reservoir Poroy, inland reservoir Sinyata Reka) and in four waterbodies sampled in 2019 (inland reservoirs Duvanli, Koprinka, Plachidol 2, Sinyata Reka). The detected genetic diversity generally corresponds to the observations conducted by conventional light microscopy, by which we distinguished three species of Cuspidothrix (Cuspidothrix issatschenkoi, Cuspidothrix elenkinii and Cuspidothrix tropicalis, the latter considered alien in the country). Eleven strains showed high similarity to two sequences of C. issatschenkoi available from the National Centre for Biotechnology Information (NCBI). Ten other strains assembled in a group, which-in lack of available from NCBI genetic sequences-were presumed related to C. tropicalis and C. elenkinii after comparison with the results from light microscopy. Cuspidothrix strains found in Bulgarian waterbodies showed high genetic similarity to those isolated and sequenced from Asia (Japan, China) and Northern Europe (Norway, Finland).

7-epi-cylindrospermopsin and microcystin producers among diverse Anabaena/Dolichospermum/Aphanizomenon CyanoHABs in Oregon, USA

Several genomes of Nostocales ADA clade members from the US Pacific Northwest were recently sequenced. Biosynthetic genes for microcystin, cylindrospermopsin or anatoxin-a were present in 7 of the 15 Dolichospermum/Anabaena strains and none of the 5 Aphanizomenon flos-aquae (AFA) strains. Toxin analyses (ELISA and LC-MS/MS) were conducted to quantitate and identify microcystin (MC) and cylindrospermopsin (CYN) congeners/analogs in samples dominated by Dolichospermum spp. of known genome sequence. MC-LR was the main congener produced by Dolichospermum spp. from Junipers Reservoir, Lake Billy Chinook and Odell Lake, while a congener provisionally identified as [Dha⁷]MC-HtyR was produced by a Dolichospermum sp. in Detroit Reservoir. A second Dolichospermum sp. from Detroit Reservoir was found to produce 7-epi-CYN, with 7-deoxy-CYN also present, but no CYN. The monitoring history of each of these lakes indicates the capacity for high levels of cyanotoxins during periods when Dolichospermum spp. are the dominant cyanobacteria. The diversity of ADA strains found in the US Pacific NW emphasizes the importance of these cyanobacteria as potentially toxic HAB formers in this temperate climatic region. Our results linking congener and genetic identity add data points that will help guide development of improved tools for predicting congener specificity from cyanotoxin gene sequences.

Depth-dependent influence of biochar application on the abundance and community structure of diazotrophic under sugarcane growth

Despite progress in understanding diazotrophic distribution in surface soils, few studies have investigated the distribution of diazotrophic bacteria in deeper soil layers. Here, we leveraged high-throughput sequencing (HTS) of nifH genes obtained to assess the influence of biochar amended soil (BC) and control (CK), and soil depths (0–20, 20–40 and 40–60 cm) on diazotrophic abundance and community structures, soil enzyme activities and physio-chemical properties. Multivariate ANOVA analysis revealed that soil depth had profound impact on majority of the soil parameters measured than fertilization. Although soil physio-chemical properties, enzymes activities, diazotrophic genera and enriched operational taxonomic units (OTUs) were significantly influenced across the entire soil profiles, we also observed that BC amended soil significantly increased cane stalk height and weight, nitrate (NO₃^-), ammonium (NH₄⁺), organic matter (OM), total carbon (TC) and available potassium (AK), and enhanced diazotrophic genera in soil depth 0–20 cm compared to CK treatment. Soil TC, total nitrogen (TN), OM and NH₄⁺ were the major impact factors shifting diazotrophic community structures in soil depth 0–20 cm. Overall, these results were more pronounced in 0–20 cm soil depth in BC than CK treatment.

Comparative genomics of the ADA clade within the Nostocales

The ADA clade of Nostocales cyanobacteria, a group that is prominent in current harmful algal bloom events, now includes over 40 genome sequences with the recent addition of sixteen novel sequenced genomes (Dreher et al., Harmful Algae, 2021). Fourteen genomes are complete (closed), enabling highly detailed assessments of gene content and genome architecture. ADA genomes contain 5 rRNA operons, genes expected to support a photoautotrophic and diazotrophic lifestyle, and a varied array of genes for the synthesis of bioactive secondary metabolites. Genes for the production of the taste-and-odor compound geosmin and the four major classes of cyanotoxins - anatoxin-a, cylindrospermopsin, microcystin and saxitoxin - are represented in members of the ADA clade. Notably, the gene array for the synthesis of cylindrospermopsin by Dolichospermum sp. DET69 was located on a plasmid, raising the possibility of facile horizontal transmission. However, genes supporting independent conjugative transfer of this plasmid are lacking. Further, analysis of genomic loci containing this and other cyanotoxin gene arrays shows evidence that these arrays have long-term stability and do not appear to be genomic islands easily capable of horizontal transmission to other cells. There is considerable diversity in the gene complements of individual ADA genomes, including the variable presence of physiologically important genes: genomes in three species-level subclades lack the gas vesicle genes that facilitate a planktonic lifestyle, and, surprisingly, the genome of Cuspidothrix issatschenkoi CHARLIE-1, a reported diazotroph, lacks the genes for nitrogen fixation. Notably, phylogenetically related genomes possess limited synteny, indicating a prominent role for chromosome rearrangements during ADA strain evolution. The genomes contain abundant insertion sequences and repetitive transposase genes, which could be the main drivers of genome rearrangement through active transposition and homologous recombination. No prophages were found, and no evidence of viral infection was observed in the bloom population samples from which the genomes discussed here were derived. Phages thus seem to have a limited influence on ADA evolution.

Complete genomes derived by directly sequencing freshwater bloom populations emphasize the significance of the genus level ADA clade within the Nostocales

The genome sequences of 16 Nostocales cyanobacteria have been determined. Most of them are complete or near-complete genome sequences derived by long-read metagenome sequencing of recent harmful algal blooms (HABs) in freshwater lakes without the potential bias of culture isolation. The genomes are all members of the recently recognized ADA clade (Driscoll et al., Harmful Algae, 77:93, 2018), which we argue represents a genus. We identify 10 putative species-level branches within the clade, on the basis of 91-gene phylogenomic and average nucleotide identity analyses. The assembled genomes each correspond to a single morphotype in the original sample, but distinct genomes from different HABs in some cases correspond to similar morphotypes. We present data indicating that the ADA clade is a highly significant component of current cyanobacterial HABs, including members assigned to the prevalent Dolichospermum and Aphanizomenon genera, as well as Cuspidothrix and Anabaena. In general, currently used genus and species names within the ADA clade are not monophyletic. We infer that the morphological characters routinely used in taxonomic assignments are not reliable for discriminating species within the ADA clade. Taxonomic revisions will be needed to create a genus with a single name (we recommend Anabaena) and to adopt species names that do not depend on morphological traits that lack sufficient discrimination and specificity, while recognizing the utility of some easily observable and distinct morphologies.

Exploring syntenic conservation across genomes for phylogenetic studies of organisms subjected to horizontal gene transfers: A case study with Cyanobacteria and cyanolichens

Understanding the evolutionary history of symbiotic Cyanobacteria at a fine scale is essential to unveil patterns of associations with their hosts and factors driving their spatiotemporal interactions. As for bacteria in general, Horizontal Gene Transfers (HGT) are expected to be rampant throughout their evolution, which justified the use of single-locus phylogenies in macroevolutionary studies of these photoautotrophic bacteria. Genomic approaches have greatly increased the amount of molecular data available, but the selection of orthologous, congruent genes that are more likely to reflect bacterial macroevolutionary histories remains problematic. In this study, we developed a synteny-based approach and searched for Collinear Orthologous Regions (COR), under the assumption that genes that are present in the same order and orientation across a wide monophyletic clade are less likely to have undergone HGT. We searched sixteen reference Nostocales genomes and identified 99 genes, part of 28 COR comprising three to eight genes each. We then developed a bioinformatic pipeline, designed to minimize inter-genome contamination and processed twelve Nostoc-associated lichen metagenomes. This reduced our original dataset to 90 genes representing 25 COR, which were used to infer phylogenetic relationships within Nostocales and among lichenized Cyanobacteria. This dataset was narrowed down further to 71 genes representing 22 COR by selecting only genes part of one (largest) operon per COR. We found a relatively high level of congruence among trees derived from the 90-gene dataset, but congruence was only slightly higher among genes within a COR compared to genes across COR. However, topological congruence was significantly higher among the 71 genes part of one operon per COR. Nostocales phylogenies resulting from concatenation and species tree approaches based on the 90- and 71-gene datasets were highly congruent, but the most highly supported result was obtained when using synteny, collinearity, and operon information (i.e., 71-gene dataset) as gene selection criteria, which outperformed larger datasets with more genes.

Using a lake sediment record to infer the long-term history of cyanobacteria and the recent rise of an anatoxin producing Dolichospermum sp

Lakes that experience recurrent toxic cyanobacterial harmful algae blooms (cyanoHABS) are often subject to cultural eutrophication, where landscape development and upland activities increase the nutrient inputs to the water column and fuel cyanoHABS. Few studies have focused on the response of a lake to nutrient inputs for which the natural geomorphic setting predisposes a nutrient-rich water column to already support abundant cyanobacteria. Here, we present a sediment core record from a lake surrounded by parkland that experiences recurrent cyanoHABs which produce dangerous levels of the neurotoxin, anatoxin-a, impacting the recreational use of the lake and park. Using photoautotrophic pigments in the sediment record, we establish cyanobacteria have long been part of the diverse and abundant phytoplankton community within the lake. Despite this long record, shotgun metagenome and other DNA analyses of the sediment record suggest that the current anatoxin-a producer Dolichospermum sp. WA102 only emerged to dominate the cyanobacterial community in the mid-1990s. A period of lakeshore farming that finished in the 1950s-1960s and possibly the stocking of rainbow trout fry (1970-2016) coincide with a progressive shift in primary production, together with a change in bacterial communities. Based on the history of the lake and contemporary ecology of Dolichospermum, we propose that the legacy of nutrient inputs and changes in nutrient cycling within the lake has encouraged the development of an ecosystem where the toxin producing Dolichospermum sp WA102 is highly competitive. Understanding the historical presence of cyanobacteria in the lake provides a context for current-day management strategies of cyanoHABs.

Phylogenomic Analysis of Secondary Metabolism in the Toxic Cyanobacterial Genera Anabaena, Dolichospermum and Aphanizomenon

Cyanobacteria produce an array of toxins that pose serious health risks to humans and animals. The closely related diazotrophic genera, Anabaena, Dolichospermum, and Aphanizomenon, frequently form poisonous blooms in lakes and brackish waters around the world. These genera form a complex now termed the Anabaena, Dolichospermum, and Aphanizomenon (ADA) clade and produce a greater array of toxins than any other cyanobacteria group. However, taxonomic confusion masks the distribution of toxin biosynthetic pathways in cyanobacteria. Here we obtained 11 new draft genomes to improve the understanding of toxin production in these genera. Comparison of secondary metabolite pathways in all available 31 genomes for these three genera suggests that the ability to produce microcystin, anatoxin-a, and saxitoxin is associated with specific subgroups. Each toxin gene cluster was concentrated or even limited to a certain subgroup within the ADA clade. Our results indicate that members of the ADA clade encode a variety of secondary metabolites following the phylogenetic clustering of constituent species. The newly sequenced members of the ADA clade show that phylogenetic separation of planktonic Dolichospermum and benthic Anabaena is not complete. This underscores the importance of taxonomic revision of Anabaena, Dolichospermum, and Aphanizomenon genera to reflect current phylogenomic understanding.

Biosynthesis of Anatoxins in Cyanobacteria: Identification of the Carboxy-anatoxins as the Penultimate Biosynthetic Intermediates

Anatoxin-a, homoanatoxin-a, and dihydroanatoxin-a are potent cyanobacterial neurotoxins. They are biosynthesized in cyanobacteria from proline and acetate by a pathway involving three polyketide synthases. We report the identification of carboxy-anatoxin-a, carboxy-homoanatoxin-a, and carboxy-dihydroanatoxin-a in acidic extracts of Cuspidothrix issatschenkoi CHARLIE-1, Oscillatoria sp. PCC 6506, and Cylindrospermum stagnale PCC 7417, respectively, using liquid chromatography coupled to mass spectrometry. The structure of these carboxy derivatives was confirmed by mass spectrometry and by isotopic incorporation experiments using labeled proline and acetate. Each of these three cyanobacteria only produce one carboxy-anatoxin, suggesting that these metabolites are the product of the hydrolysis by AnaA, the type II thioesterase, of the thioesters bound to AnaG, the last polyketide synthase of the pathway. By measuring the rate of isotopic incorporation of labeled proline into carboxy-homoanatoxin-a and homoanatoxin-a produced by Oscillatoria sp. PCC 6506, we show that carboxy-homoanatoxin-a is the intracellular precursor of homoanatoxin-a, and that homoanatoxin-a is then excreted into the extracellular medium. The transformation of carboxy-homoanatoxin-a into homoanatoxin-a is a very slow two-step process, with accumulation of carboxy-homoanatoxin-a, suggesting that the decarboxylation is spontaneous and not enzymatically catalyzed. However, an unidentified and extracellular catalyst accelerates the decarboxylation when the cell extracts are prepared at neutral pH.

Comprehensive Analyses of Cytochrome P450 Monooxygenases and Secondary Metabolite Biosynthetic Gene Clusters in Cyanobacteria

The prokaryotic phylum Cyanobacteria are some of the oldest known photosynthetic organisms responsible for the oxygenation of the earth. Cyanobacterial species have been recognised as a prosperous source of bioactive secondary metabolites with antibacterial, antiviral, antifungal and/or anticancer activities. Cytochrome P450 monooxygenases (CYPs/P450s) contribute to the production and diversity of various secondary metabolites. To better understand the metabolic potential of cyanobacterial species, we have carried out comprehensive analyses of P450s, predicted secondary metabolite biosynthetic gene clusters (BGCs), and P450s located in secondary metabolite BGCs. Analysis of the genomes of 114 cyanobacterial species identified 341 P450s in 88 species, belonging to 36 families and 79 subfamilies. In total, 770 secondary metabolite BGCs were found in 103 cyanobacterial species. Only 8% of P450s were found to be part of BGCs. Comparative analyses with other bacteria Bacillus, Streptomyces and mycobacterial species have revealed a lower number of P450s and BGCs and a percentage of P450s forming part of BGCs in cyanobacterial species. A mathematical formula presented in this study revealed that cyanobacterial species have the highest gene-cluster diversity percentage compared to Bacillus and mycobacterial species, indicating that these diverse gene clusters are destined to produce different types of secondary metabolites. The study provides fundamental knowledge of P450s and those associated with secondary metabolism in cyanobacterial species, which may illuminate their value for the pharmaceutical and cosmetics industries.

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.

Variability in the anatoxin gene clusters of Cuspidothrix issatschenkoi from Germany, New Zealand, China and Japan

Anatoxin-a and homoanatoxin-a are neurotoxic cyanotoxins produced by benthic and planktonic cyanobacteria worldwide. These toxins are produced by the cyanobacterial genera Dolichospermum, Cuspidothrix, Phormidium, Oscillatoria, Tychonema and Cylindrospermum. In the present study the ana gene clusters (anaA-anaG; c. 21.1 kilobases) of two anatoxin producing Cuspidothrix issatschenkoi strains from Germany: (NIVA-CYA 711) and New Zealand (CAWBG02) were sequenced and compared with the ana gene clusters of two C. issatschenkoi strains from Japan (RM-6 and LBRI48) and one from China (CHABD3). All five ana gene clusters are characterized by the same gene order for anaA-anaG. Similarities were highest (99.56–99.57%) between German (NIVA-CYA 711), New Zealand (CAWBG02) and Chinese (CHABD3) strains. Similarities were lower (91.40–91.67%) when compared to the Japanese strains (RM-6 and LBRI48). Collectively, 2,037 variable sites (328 single nucleotide polymorphisms and 9 insertions/deletions, comprising 1,709 nucleotides) were found in the ana gene clusters of the German, New Zealand and Japanese strains compared to the Chinese strain (CHABD3). The ana gene clusters of the German (NIVA-CYA 711), New Zealand (CAWBG02) and Japanese (RM-6 and LBRI48) strains were characterized by 83, 84, 255 and 231 SNP’s compared to the Chinese strain (CHABD3), respectively. The anaE and anaF genes showed the highest variability in all five strains and are recommended as the best genetic markers for further phylogenetic studies of the ana gene cluster from C. issatschenkoi.

A closely-related clade of globally distributed bloom-forming cyanobacteria within the Nostocales

In order to better understand the relationships among current Nostocales cyanobacterial blooms, eight genomes were sequenced from cultured isolates or from environmental metagenomes of recent planktonic Nostocales blooms. Phylogenomic analysis of publicly available sequences placed the new genomes among a group of 15 genomes from four continents in a distinct ADA clade (Anabaena/Dolichospermum/Aphanizomenon) within the Nostocales. This clade contains four species-level groups, two of which include members with both Anabaena-like and Aphanizomenon flos-aquae-like morphology. The genomes contain many repetitive genetic elements and a sizable pangenome, in which ABC-type transporters are highly represented. Alongside common core genes for photosynthesis, the differentiation of N₂-fixing heterocysts, and the uptake and incorporation of the major nutrients P, N and S, we identified several gene pathways in the pangenome that may contribute to niche partitioning. Genes for problematic secondary metabolites-cyanotoxins and taste-and-odor compounds-were sporadically present, as were other polyketide synthase (PKS) and nonribosomal peptide synthetase (NRPS) gene clusters. By contrast, genes predicted to encode the ribosomally generated bacteriocin peptides were found in all genomes.

Optofluidic Single-Cell Genome Amplification of Sub-micron Bacteria in the Ocean Subsurface

Optofluidic single-cell genome amplification was used to obtain genome sequences from sub-micron cells collected from the euphotic and mesopelagic zones of the northwestern Sargasso Sea. Plankton cells were visually selected and manually sorted with an optical trap, yielding 20 partial genome sequences representing seven bacterial phyla. Two organisms, E01-9C-26 (Gammaproteobacteria), represented by four single cell genomes, and Opi.OSU.00C, an uncharacterized Verrucomicrobia, were the first of their types retrieved by single cell genome sequencing and were studied in detail. Metagenomic data showed that E01-9C-26 is found throughout the dark ocean, while Opi.OSU.00C was observed to bloom transiently in the nutrient-depleted euphotic zone of the late spring and early summer. The E01-9C-26 genomes had an estimated size of 4.76-5.05 Mbps, and contained "O" and "W"-type monooxygenase genes related to methane and ammonium monooxygenases that were previously reported from ocean metagenomes. Metabolic reconstruction indicated E01-9C-26 are likely versatile methylotrophs capable of scavenging C1 compounds, methylated compounds, reduced sulfur compounds, and a wide range of amines, including D-amino acids. The genome sequences identified E01-9C-26 as a source of "O" and "W"-type monooxygenase genes related to methane and ammonium monooxygenases that were previously reported from ocean metagenomes, but are of unknown function. In contrast, Opi.OSU.00C genomes encode genes for catabolizing carbohydrate compounds normally associated with eukaryotic phytoplankton. This exploration of optofluidics showed that it was effective for retrieving diverse single-cell bacterioplankton genomes and has potential advantages in microbiology applications that require working with small sample volumes or targeting cells by their morphology.

Benthic periphyton from Pennsylvania, USA is a source for both hepatotoxins (microcystins/nodularin) and neurotoxins (anatoxin-a/homoanatoxin-a)

In 2016, the Pennsylvania Department of Environmental Protection conducted a limited survey of streams in the Susquehanna River basin in Pennsylvania, USA, to screen for microcystins/nodularins, anatoxin-a (ATX) and homoanatoxin-a (HTX). Testing revealed the presence of HTX in samples collected from the Pine Creek basin, with ATX present at lower levels. Microcystins/nodularins (MCs/NODs) were also tested and found to be concomitant, with NOD-R confirmed present by LC-MS/MS.

Recovering Genomics Clusters of Secondary Metabolites from Lakes Using Genome-Resolved Metagenomics

Metagenomic approaches became increasingly popular in the past decades due to decreasing costs of DNA sequencing and bioinformatics development. So far, however, the recovery of long genes coding for secondary metabolites still represents a big challenge. Often, the quality of metagenome assemblies is poor, especially in environments with a high microbial diversity where sequence coverage is low and complexity of natural communities high. Recently, new and improved algorithms for binning environmental reads and contigs have been developed to overcome such limitations. Some of these algorithms use a similarity detection approach to classify the obtained reads into taxonomical units and to assemble draft genomes. This approach, however, is quite limited since it can classify exclusively sequences similar to those available (and well classified) in the databases. In this work, we used draft genomes from Lake Stechlin, north-eastern Germany, recovered by MetaBat, an efficient binning tool that integrates empirical probabilistic distances of genome abundance, and tetranucleotide frequency for accurate metagenome binning. These genomes were screened for secondary metabolism genes, such as polyketide synthases (PKS) and non-ribosomal peptide synthases (NRPS), using the Anti-SMASH and NAPDOS workflows. With this approach we were able to identify 243 secondary metabolite clusters from 121 genomes recovered from our lake samples. A total of 18 NRPS, 19 PKS, and 3 hybrid PKS/NRPS clusters were found. In addition, it was possible to predict the partial structure of several secondary metabolite clusters allowing for taxonomical classifications and phylogenetic inferences. Our approach revealed a high potential to recover and study secondary metabolites genes from any aquatic ecosystem.

SAR202 Genomes from the Dark Ocean Predict Pathways for the Oxidation of Recalcitrant Dissolved Organic Matter

Deep-ocean regions beyond the reach of sunlight contain an estimated 615 Pg of dissolved organic matter (DOM), much of which persists for thousands of years. It is thought that bacteria oxidize DOM until it is too dilute or refractory to support microbial activity. We analyzed five single-amplified genomes (SAGs) from the abundant SAR202 clade of dark-ocean bacterioplankton and found they encode multiple families of paralogous enzymes involved in carbon catabolism, including several families of oxidative enzymes that we hypothesize participate in the degradation of cyclic alkanes. The five partial genomes encoded 152 flavin mononucleotide/F420-dependent monooxygenases (FMNOs), many of which are predicted to be type II Baeyer-Villiger monooxygenases (BVMOs) that catalyze oxygen insertion into semilabile alicyclic alkanes. The large number of oxidative enzymes, as well as other families of enzymes that appear to play complementary roles in catabolic pathways, suggests that SAR202 might catalyze final steps in the biological oxidation of relatively recalcitrant organic compounds to refractory compounds that persist.IMPORTANCE Carbon in the ocean is massively sequestered in a complex mixture of biologically refractory molecules that accumulate as the chemical end member of biological oxidation and diagenetic change. However, few details are known about the biochemical machinery of carbon sequestration in the deep ocean. Reconstruction of the metabolism of a deep-ocean microbial clade, SAR202, led to postulation of new biochemical pathways that may be the penultimate stages of DOM oxidation to refractory forms that persist. These pathways are tied to a proliferation of oxidative enzymes. This research illuminates dark-ocean biochemistry that is broadly consequential for reconstructing the global carbon cycle.

Towards long-read metagenomics: complete assembly of three novel genomes from bacteria dependent on a diazotrophic cyanobacterium in a freshwater lake co-culture

Here we report three complete bacterial genome assemblies from a PacBio shotgun metagenome of a co-culture from Upper Klamath Lake, OR. Genome annotations and culture conditions indicate these bacteria are dependent on carbon and nitrogen fixation from the cyanobacterium Aphanizomenon flos-aquae, whose genome was assembled to draft-quality. Due to their taxonomic novelty relative to previously sequenced bacteria, we have temporarily designated these bacteria as incertae sedis Hyphomonadaceae strain UKL13-1 (3,501,508 bp and 56.12% GC), incertae sedis Betaproteobacterium strain UKL13-2 (3,387,087 bp and 54.98% GC), and incertae sedis Bacteroidetes strain UKL13-3 (3,236,529 bp and 37.33% GC). Each genome consists of a single circular chromosome with no identified plasmids. When compared with binned Illumina assemblies of the same three genomes, there was ~7% discrepancy in total genome length. Gaps where Illumina assemblies broke were often due to repetitive elements. Within these missing sequences were essential genes and genes associated with a variety of functional categories. Annotated gene content reveals that both Proteobacteria are aerobic anoxygenic phototrophs, with Betaproteobacterium UKL13-2 potentially capable of phototrophic oxidation of sulfur compounds. Both proteobacterial genomes contain transporters suggesting they are scavenging fixed nitrogen from A. flos-aquae in the form of ammonium. Bacteroidetes UKL13-3 has few completely annotated biosynthetic pathways, and has a comparatively higher proportion of unannotated genes. The genomes were detected in only a few other freshwater metagenomes, suggesting that these bacteria are not ubiquitous in freshwater systems. Our results indicate that long-read sequencing is a viable method for sequencing dominant members from low-diversity microbial communities, and should be considered for environmental metagenomics when conditions meet these requirements.

Moving towards adaptive management of cyanotoxin-impaired water bodies

The cyanobacteria are a phylum of bacteria that have played a key role in shaping the Earth's biosphere due to their pioneering ability to perform oxygenic photosynthesis. Throughout their history, cyanobacteria have experienced major biogeochemical changes accompanying Earth's geochemical evolution over the past 2.5+ billion years, including periods of extreme climatic change, hydrologic, nutrient and radiation stress. Today, they remain remarkably successful, exploiting human nutrient over‐enrichment as nuisance “blooms.” Cyanobacteria produce an array of unique metabolites, the functions and biotic ramifications of which are the subject of diverse ecophysiological studies. These metabolites are relevant from organismal and ecosystem function perspectives because some can be toxic and fatal to diverse biota, including zooplankton and fish consumers of algal biomass, and high‐level consumers of aquatic food sources and drinking water, including humans. Given the long history of environmental extremes and selection pressures that cyanobacteria have experienced, it is likely that that these toxins serve ecophysiological functions aimed at optimizing growth and fitness during periods of environmental stress. Here, we explore the molecular and ecophysiological mechanisms underlying cyanotoxin production, with emphasis on key environmental conditions potentially controlling toxin production. Based on this information, we offer potential management strategies for reducing cyanotoxin potentials in natural waters; for cyanotoxins with no clear drivers yet elucidated, we highlight the data gaps and research questions that are still lacking. We focus on the four major classes of toxins (anatoxins, cylindrospermopsins, microcystins and saxitoxins) that have thus far been identified as relevant from environmental health perspectives, but caution there may be other harmful metabolites waiting to be elucidated.