Are Oligotypes Meaningful Ecological and Phylogenetic Units? A Case Study of Microcystis in Freshwater Lakes

Comparative genomic analysis of Microcystis strain diversity using conserved marker genes

Microcystis-dominated cyanobacterial harmful algal blooms (cyanoHABs) have a global impact on freshwater environments, affecting both wildlife and human health. Microcystis diversity and function in field samples and laboratory cultures can be determined by sequencing whole genomes of cultured isolates or natural populations, but these methods remain computationally and financially expensive. Amplicon sequencing of marker genes is a lower cost and higher throughput alternative to characterize strain composition and diversity in mixed samples. However, the selection of appropriate marker gene region(s) and primers requires prior understanding of the relationship between single gene genotype, whole genome content, and phenotype. To identify phylogenetic markers of Microcystis strain diversity, we compared phylogenetic trees built from each of 2,351 individual core genes to an established phylogeny and assessed the ability of these core genes to predict whole genome content and bioactive compound genotypes. We identified single-copy core genes better able to resolve Microcystis phylogenies than previously identified marker genes. We developed primers suitable for current Illumina-based amplicon sequencing with near-complete coverage of available Microcystis genomes and demonstrate that they outperform existing options for assessing Microcystis strain composition. Results showed that genetic markers can be used to infer Microcystis gene content and phenotypes such as potential production of bioactive compounds , although marker performance varies by bioactive compound gene and sequence similarity. Finally, we demonstrate that these markers can be used to characterize the Microcystis strain composition of laboratory or field samples like those collected for surveillance and modeling of Microcystis-dominated cyanobacterial harmful algal blooms.

Spatiotemporal diversity and community structure of cyanobacteria and associated bacteria in the large shallow subtropical Lake Okeechobee (Florida, United States)

Lake Okeechobee is a large eutrophic, shallow, subtropical lake in south Florida, United States. Due to decades of nutrient loading and phosphorus rich sediments, the lake is eutrophic and frequently experiences cyanobacterial harmful algal blooms (cyanoHABs). In the past, surveys of the phytoplankton community structure in the lake have been conducted by morphological studies, whereas molecular based studies have been seldom employed. With increased frequency of cyanoHABs in Lake Okeechobee (e.g., 2016 and 2018 Microcystis-dominated blooms), it is imperative to determine the diversity of cyanobacterial taxa that exist within the lake and the limnological parameters that drive bloom-forming genera. A spatiotemporal study of the lake was conducted over the course of 1 year to characterize the (cyano)bacterial community structure, using 16S rRNA metabarcoding, with coincident collection of limnological parameters (e.g., nutrients, water temperature, major ions), and cyanotoxins. The objectives of this study were to elucidate spatiotemporal trends of community structure, identify drivers of community structure, and examine cyanobacteria-bacterial relationships within the lake. Results indicated that cyanobacterial communities within the lake were significantly different between the wet and dry season, but not between periods of nitrogen limitation and co-nutrient limitation. Throughout the year, the lake was primarily dominated by the picocyanobacterium Cyanobium. The bloom-forming genera Cuspidothrix, Dolichospermum, Microcystis, and Raphidiopsis were highly abundant throughout the lake and had disparate nutrient requirements and niches within the lake. Anatoxin-a, microcystins, and nodularins were detected throughout the lake across both seasons. There were no correlated (cyano)bacteria shared between the common bloom-forming cyanobacteria Dolichospermum, Microcystis, and Raphidiopsis. This study is the first of its kind to use molecular based methods to assess the cyanobacterial community structure within the lake. These data greatly improve our understanding of the cyanobacterial community structure within the lake and the physiochemical parameters which may drive the bloom-forming taxa within Lake Okeechobee.

Impacts of an invasive filter-feeder on bacterial biodiversity are context dependent

Bacteria represent most of the biodiversity and play key roles in virtually every ecosystem. In doing so, bacteria act as part of complex communities shaped by interactions across all domains of life. Here, we report on direct interactions between bacteria and dreissenid mussels, a group of invasive filter-feeders threatening global aquatic systems due to high filtration rates. Previous studies showed that dreissenids can impact bacterial community structure by changing trait distributions and abundances of specific taxa. However, studies on bacterial community effects were conducted using water from Lake Michigan (an oligotrophic lake) only, and it is unknown whether similar patterns are observed in systems with differing nutrient regimes. We conducted ten short-term dreissenid grazing experiments in 2019 using water from two eutrophic lake regions-the western basin of Lake Erie and Saginaw Bay in Lake Huron. Predation by dreissenids led to decline in overall bacterial abundance and diversity in both lakes. However, feeding on bacteria was not observed during every experiment. We also found that traits related to feeding resistance are less phylogenetically conserved than previously thought. Our results highlight the role of temporal, spatial, and genomic heterogeneity in bacterial response dynamics to a globally important invasive filter feeder.

Uptake of Phytoplankton-Derived Carbon and Cobalamins by Novel Acidobacteria Genera in Microcystis Blooms Inferred from Metagenomic and Metatranscriptomic Evidence

Interactions between bacteria and phytoplankton can influence primary production, community composition, and algal bloom development. However, these interactions are poorly described for many consortia, particularly for freshwater bloom-forming cyanobacteria. Here, we assessed the gene content and expression of two uncultivated Acidobacteria from Lake Erie Microcystis blooms. These organisms were targeted because they were previously identified as important catalase producers in Microcystis blooms, suggesting that they protect Microcystis from H₂O₂. Metatranscriptomics revealed that both Acidobacteria transcribed genes for uptake of organic compounds that are known cyanobacterial products and exudates, including lactate, glycolate, amino acids, peptides, and cobalamins. Expressed genes for amino acid metabolism and peptide transport and degradation suggest that use of amino acids and peptides by Acidobacteria may regenerate nitrogen for cyanobacteria and other organisms. The Acidobacteria genomes lacked genes for biosynthesis of cobalamins but expressed genes for its transport and remodeling. This indicates that the Acidobacteria obtained cobalamins externally, potentially from Microcystis, which has a complete gene repertoire for pseudocobalamin biosynthesis; expressed them in field samples; and produced pseudocobalamin in axenic culture. Both Acidobacteria were detected in Microcystis blooms worldwide. Together, the data support the hypotheses that uncultured and previously unidentified Acidobacteria taxa exchange metabolites with phytoplankton during harmful cyanobacterial blooms and influence nitrogen available to phytoplankton. Thus, novel Acidobacteria may play a role in cyanobacterial physiology and bloom development. IMPORTANCE Interactions between heterotrophic bacteria and phytoplankton influence competition and successions between phytoplankton taxa, thereby influencing ecosystem-wide processes such as carbon cycling and algal bloom development. The cyanobacterium Microcystis forms harmful blooms in freshwaters worldwide and grows in buoyant colonies that harbor other bacteria in their phycospheres. Bacteria in the phycosphere and in the surrounding community likely influence Microcystis physiology and ecology and thus the development of freshwater harmful cyanobacterial blooms. However, the impacts and mechanisms of interaction between bacteria and Microcystis are not fully understood. This study explores the mechanisms of interaction between Microcystis and uncultured members of its phycosphere in situ with population genome resolution to investigate the cooccurrence of Microcystis and freshwater Acidobacteria in blooms worldwide.

Metagenomic and Metatranscriptomic Insights into Population Diversity of Microcystis Blooms: Spatial and Temporal Dynamics of mcy Genotypes, Including a Partial Operon That Can Be Abundant and Expressed

Cyanobacterial harmful algal blooms (cyanoHABs) degrade freshwater ecosystems globally. Microcystis aeruginosa often dominates cyanoHABs and produces microcystin (MC), a class of hepatotoxins that poses threats to human and animal health. Microcystin toxicity is influenced by distinct structural elements across a diversity of related molecules encoded by variant mcy operons. However, the composition and distribution of mcy operon variants in natural blooms remain poorly understood. Here, we characterized the variant composition of mcy genes in western Lake Erie Microcystis blooms from 2014 and 2018. Sampling was conducted across several spatial and temporal scales, including different bloom phases within 2014, extensive spatial coverage on the same day (2018), and frequent, autonomous sampling over a 2-week period (2018). Mapping of metagenomic and metatranscriptomic sequences to reference sequences revealed three Microcystis mcy genotypes: complete (all genes present [mcyA-J]), partial (truncated mcyA, complete mcyBC, and missing mcyD-J), and absent (no mcy genes). We also detected two different variants of mcyB that may influence the production of microcystin congeners. The relative abundance of these genotypes was correlated with pH and nitrate concentrations. Metatranscriptomic analysis revealed that partial operons were, at times, the most abundant genotype and expressed in situ, suggesting the potential biosynthesis of truncated products. Quantification of genetic divergence between genotypes suggests that the observed strains are the result of preexisting heterogeneity rather than de novo mutation during the sampling period. Overall, our results show that natural Microcystis populations contain several cooccurring mcy genotypes that dynamically shift in abundance spatiotemporally via strain succession and likely influence the observed diversity of the produced congeners. IMPORTANCE Cyanobacteria are responsible for producing microcystins (MCs), a class of potent and structurally diverse toxins, in freshwater systems around the world. While microcystins have been studied for over 50 years, the diversity of their chemical forms and how this variation is encoded at the genetic level remain poorly understood, especially within natural populations of cyanobacterial harmful algal blooms (cyanoHABs). Here, we leverage community DNA and RNA sequences to track shifts in mcy genes responsible for producing microcystin, uncovering the relative abundance, expression, and variation of these genes. We studied this phenomenon in western Lake Erie, which suffers annually from cyanoHAB events, with impacts on drinking water, recreation, tourism, and commercial fishing.

Comparing microscopy and DNA metabarcoding techniques for identifying cyanobacteria assemblages across hundreds of lakes

Accurately identifying the species present in an ecosystem is vital to lake managers and successful bioassessment programs. This is particularly important when monitoring cyanobacteria, as numerous taxa produce toxins and can have major negative impacts on aquatic ecosystems. Increasingly, DNA-based techniques such as metabarcoding are being used for measuring aquatic biodiversity, as they could accelerate processing time, decrease costs and reduce some of the biases associated with traditional light microscopy. Despite the continuing use of traditional microscopy and the growing use of DNA metabarcoding to identify cyanobacteria assemblages, methodological comparisons between the two approaches have rarely been reported from a wide suite of lake types. Here, we compare planktonic cyanobacteria assemblages generated by inverted light microscopy and DNA metabarcoding from a 379-lake dataset spanning a longitudinal and trophic gradient. We found moderate levels of congruence between methods at the broadest taxonomic levels (i.e., Order, RV=0.40, p < 0.0001). This comparison revealed distinct cyanobacteria communities from lakes of different trophic states, with Microcystis, Aphanizomenon and Dolichospermum dominating with both methods in eutrophic and hypereutrophic sites. This finding supports the use of either method when monitoring eutrophication in lake surface waters. The biggest difference between the two methods was the detection of picocyanobacteria, which are typically underestimated by light microscopy. This reveals that the communities generated by each method currently are complementary as opposed to identical and promotes a combined-method strategy when monitoring a range of trophic systems. For example, microscopy can provide measures of cyanobacteria biomass, which are critical data in managing lakes. Going forward, we believe that molecular genetic methods will be increasingly adopted as reference databases are routinely updated with more representative sequences and will improve as cyanobacteria taxonomy is resolved with the increase in available genetic information.

Bacterioplankton Zonation Does Exist in High Elevation, Polymictic Lakes

The assessment of distribution patterns or zonation of planktonic microbes along the water column is a crucial step to interpret their function in the ecosystem. In lakes without seasonal thermal stratification or polymictic systems such as high elevation tropical lakes, planktonic bacterial taxa are probably homogeneously distributed in the water column in contrast to what is known for thermally stratified lakes. However, we know little about bacterial distribution patterns in polymictic lakes and their relation to environmental gradients other than temperature. Here we assessed the diversity, microdiversity, and bacterial community composition at different discrete depths in three high elevation lakes (4,400-4,550 m above sea level) from the Andean plateau to test whether bacterial zonation patterns exist along the water column. For this objective, we analyzed bulk DNA and the putatively active fraction (cDNA) of the 16S rRNA gene. Although a clear gradient of temperature and oxygen was not detected along the water column, a significant vertical spatial zonation of the bacterial communities was present in two out of the three lakes, with microdiversity contributing to such pattern. Our results provide a reference for understanding how changing environmental conditions could affect high elevation aquatic ecosystems, particularly when warming is amplified with elevation, accelerating changes in hydrological regimes and biodiversity. Finally, our results highlight the importance of incorporating the whole water column in ecological studies of aquatic ecosystems lacking temporal or permanent thermal stratification.

Patterns of Relative Bacterial Richness and Community Composition in Seawater and Marine Sediment Are Robust for Both Operational Taxonomic Units and Amplicon Sequence Variants

To understand the relative influences of operational taxonomic units (OTUs) and amplicon sequence variants (ASVs) on patterns of marine microbial diversity and community composition, we examined bacterial diversity and community composition of seawater from 12 sites in the North Atlantic Ocean and Canadian Arctic and sediment from two sites in the North Atlantic. For the seawater analyses, we included samples from three to six zones in the water column of each site. For the sediment analyses, we included over 20 sediment horizons at each of two sites. For all samples, we amplified the V4-V5 hypervariable region of the 16S ribosomal RNA (rRNA) gene. We analyzed each sample in two different ways: (i) by clustering its reads into 97%-similar OTUs and (ii) by assigning sequences to unique ASVs. OTU richness is much higher than ASV richness for every sample, but both OTUs and ASVs exhibit similar vertical patterns of relative diversity in both the water column and the sediment. Bacterial richness is highest just below the photic zone in the water column and at the seafloor in the sediment. For both OTUs and ASVs, richness estimates depend on the number of sequences analyzed. Both methods yield broadly similar community compositions for each sample at the taxonomic levels of phyla to families. While the two methods yield different richness values, broad-scale patterns of relative richness and community composition are similar with both methods.

Deep sequencing across multiple host species tests pine-endophyte specificity

PREMISE

Foliar fungal endophytes vary in their distributions across landscapes or plant host taxa, indicative of specialized ecologies and host specific adaptations. Accounts of specialization, however, depend on the taxonomic breadth and geographic range of the host plants included in each study. A broad region-scale study or deep sampling of diverse potential host species still remains relatively rare but is becoming increasingly possible with high-throughput sequencing.

METHODS

Amplicon sequencing was used to rapidly identify the fungal endophytic community among six pine (Pinus, Pinaceae) species co-occurring across northeastern United States and to test for site and host specialization. We focused on the endophytic genus Lophodermium (Rhytismataceae), whose species members are thought to specialize on different pine species, to test if amplicon sequencing could rapidly verify previously implied or discover new patterns of host specificity.

RESULTS

While amplicon sequencing could analyze more samples at greater depths and recover greater numbers of unique Lophodermium taxa than when endophyte communities were surveyed with traditional culturing methods, patterns of specialization were not better supported. This may be because amplicon sequencing can indiscriminately capture non-host specific organisms found incidentally from plant tissues or because we have overestimated host-specificity in the past with biased culturing techniques.

CONCLUSIONS

Amplicon sequencing can quickly identify patterns of host specificity by allowing large-scale surveys but has limitations in quantifying the level of intimacy of these relationships.

Impact of Electrolyzed Water on the Microbial Spoilage Profile of Piedmontese Steak Tartare

A low initial contamination level of the meat surface is the sine qua non to extend the subsequent shelf life of ground beef for as long as possible. Therefore, the short- and long-term effects of a pregrinding treatment with electrolyzed water (EW) on the microbiological and physicochemical features of Piedmontese steak tartare were here assessed on site, by following two production runs through storage under vacuum packaging conditions at 4°C. The immersion of muscle meat in EW solution at 100 ppm of free active chlorine for 90 s produced an initial surface decontamination with no side effects or compositional modifications, except for an external color change that was subsequently masked by the grinding step. However, the initially measured decontamination was no longer detectable in ground beef, perhaps due to a quick recovery by bacteria during the grinding step from the transient oxidative stress induced by the EW. We observed different RNA-based metataxonomic profiles and metabolomic biomarkers (volatile organic compounds [VOCs], free amino acids [FAA], and biogenic amines [BA]) between production runs. Interestingly, the potentially active microbiota of the meat from each production run, investigated through operational taxonomic unit (OTU)-, oligotyping-, and amplicon sequence variant (ASV)-based bioinformatic pipelines, differed as soon as the early stages of storage, whereas microbial counts and biomarker dynamics were significantly distinguishable only after the expiration date. Higher diversity, richness, and abundance of Streptococcus organisms were identified as the main indicators of the faster spoilage observed in one of the two production runs, while Lactococcus piscium development was the main marker of shelf life end in both production runs. IMPORTANCE Treatment with EW prior to grinding did not result in an effective intervention to prolong the shelf life of Piedmontese steak tartare. Our RNA-based approach clearly highlighted a microbiota that changed markedly between production runs but little during the first shelf life stages. Under these conditions, an early metataxonomic profiling might provide the best prediction of the microbiological fate of each batch of the product.

Genotyping and multivariate regression trees reveal ecological diversification within the Microcystis aeruginosa complex along a wide environmental gradient

Addressing the ecological and evolutionary processes underlying biodiversity patterns is essential to identify the mechanisms shaping community structure and function. In bacteria, the formation of new ecologically distinct populations (ecotypes) is proposed as one of the main drivers of diversification. New ecotypes arise when mutations in key functional genes or acquisition of new metabolic pathways by horizontal gene transfer allow the population to exploit new resources, permitting their coexistence with the parental population. We previously reported the presence of microcystin-producing organisms of the Microcystis aeruginosa complex (toxic MAC) through an 800 km environmental gradient ranging from freshwater to estuarine-marine waters in South America. We hypothesize that the success of toxic MAC in such a gradient is due to the existence of very closely related populations that are ecologically distinct (ecotypes), each specialized to a specific arrangement of environmental variables. Here, we analyzed toxic MAC genetic diversity through qPCR and high-resolution melting analysis (HRMA) of a functional gene (mcyJ, microcystin synthetase cluster). We explored the variability of the mcyJ gene along the environmental gradient by multivariate classification and regression trees (mCART). Six groups of mcyJ genotypes were distinguished and associated with different combinations of water temperature, conductivity and turbidity. We propose that each mcyJ variant associated to a defined environmental condition is an ecotype (or species) whose relative abundances vary according to their fitness in the local environment. This mechanism would explain the success of toxic MAC in such a wide array of environmental conditions. Importance Organisms of the Microcystis aeruginosa Complex form harmful algal blooms (HABs) in nutrient-rich water bodies worldwide. MAC HABs are difficult to manage owing to the production of potent toxins (microcystins) that resist water treatment. Besides, the role of microcystins in the ecology of MAC organisms is still elusive, meaning that the environmental conditions driving the toxicity of the bloom are not clear. Furthermore, the lack of coherence between morphology-based and genomic-based species classification makes it difficult to draw sound conclusions about when and where each member species of the MAC will dominate the bloom. Here, we propose that the diversification process and success of toxic MAC in a wide range of waterbodies involves the generation of ecotypes, each specialized in a particular niche, whose relative abundance varies according to its fitness in the local environment. This knowledge can improve the generation of accurate prediction models of MAC growth and toxicity, helping to prevent human and animal intoxication.

Oligotyping and metagenomics reveal distinct Candidatus Accumulibacter communities in side-stream versus conventional full-scale enhanced biological phosphorus removal (EBPR) systems

Candidatus Accumulibacter phosphatis (CAP) and its clade-level micro-diversity has been associated with and implicated in functional differences in phosphorus removal performance in enhanced biological phosphorus removal (EBPR) systems. Side-stream EBPR (S2EBPR) is an emerging process that has been shown to present a suite of advantages over the conventional EBPR design, however, large knowledge gaps remain in terms of its underlying ecological mechanisms. Here, we compared and revealed the higher-resolution differences in microbial ecology of CAP between a full-scale side-stream EBPR configuration and a conventional A2O EBPR process that were operated in parallel and with the same influent feed. Even though the relative abundance of CAP, revealed by 16S rRNA gene amplicon sequencing, was similar in both treatment trains, a clade-level analysis, using combined 16S rRNA-gene based amplicon sequencing and oligotyping analysis and metagenomics analysis, revealed the distinct CAP microdiversity between the S2EBPR and A2O configurations that likely attributed to the improved performance in S2EBPR in comparison to conventional EBPR. Furthermore, genome-resolved metagenomics enabled extraction of three metagenome-assembled genomes (MAGs) belonging to CAP clades IIB (RCAB4-2), IIC (RC14) and II (RC18), from full-scale EBPR sludge for the first time, including a distinct Ca. Accumulibacter clade that is dominant and associated only with the S2EBPR configuration. The results also revealed the temporally increasing predominance of RC14, which belonged to Clade IIC, during the implementation of the S2EBPR configuration. Finally, we also show the existence of previously uncharacterized diversity of clades of CAP, namely the clades IIB and as yet unidentified clade of type II, in full-scale EBPR communities, highlighting the unknown diversity of CAP communities in full-scale EBPR systems.

Rapid Cyanobacteria Species Identification with High Sensitivity Using Native Mass Spectrometry

Cyanobacteria have evolved over billions of years to adapt and survive in diverse climates. Environmentally, this presents a huge challenge because cyanobacteria can now rapidly form algae blooms that are detrimental to aquatic life. In addition, many cyanobacteria produce toxins, making them hazardous to animals and humans that they encounter. Rapid identification of cyanobacteria is essential to monitor and prevent toxic algae blooms. Here, we show for the first time how native mass spectrometry can quickly and precisely identify cyanobacteria from diverse aquatic environments. By monitoring phycobiliproteins, abundant protein complexes within cyanobacteria, simple, easy-to-understand mass spectral "fingerprints" were created that were unique to each species. Moreover, our method is 10-fold more sensitive than the current MALDI-TOF mass spectrometric methods, meaning that cyanobacteria can be monitored using this technology prior to bloom formation. Together, the data show great promise for the simultaneous detection and identification of co-existing cyanobacteria in situ.

Individual Microcystis colonies harbour distinct bacterial communities that differ by Microcystis oligotype and with time

Interactions between bacteria and phytoplankton in the phycosphere have impacts at the scale of whole ecosystems, including the development of harmful algal blooms. The cyanobacterium Microcystis causes toxic blooms that threaten freshwater ecosystems and human health globally. Microcystis grows in colonies that harbour dense assemblages of other bacteria, yet the taxonomic composition of these phycosphere communities and the nature of their interactions with Microcystis are not well characterized. To identify the taxa and compositional variance within Microcystis phycosphere communities, we performed 16S rRNA V4 region amplicon sequencing on individual Microcystis colonies collected biweekly via high-throughput droplet encapsulation during a western Lake Erie cyanobacterial bloom. The Microcystis phycosphere communities were distinct from microbial communities in whole water and bulk phytoplankton seston in western Lake Erie but lacked 'core' taxa found across all colonies. However, dissimilarity in phycosphere community composition correlated with sampling date and the Microcystis 16S rRNA oligotype. Several taxa in the phycosphere were specific to and conserved with Microcystis of a single oligotype or sampling date. Together, this suggests that physiological differences between Microcystis strains, temporal changes in strain phenotypes, and the composition of seeding communities may impact community composition of the Microcystis phycosphere.

Metagenomics Approaches for the Detection and Surveillance of Emerging and Recurrent Plant Pathogens

Globalization has a dramatic effect on the trade and movement of seeds, fruits and vegetables, with a corresponding increase in economic losses caused by the introduction of transboundary plant pathogens. Current diagnostic techniques provide a useful and precise tool to enact surveillance protocols regarding specific organisms, but this approach is strictly targeted, while metabarcoding and shotgun metagenomics could be used to simultaneously detect all known pathogens and potentially new ones. This review aims to present the current status of high-throughput sequencing (HTS) diagnostics of fungal and bacterial plant pathogens, discuss the challenges that need to be addressed, and provide direction for the development of methods for the detection of a restricted number of related taxa (specific surveillance) or all of the microorganisms present in a sample (general surveillance). HTS techniques, particularly metabarcoding, could be useful for the surveillance of soilborne, seedborne and airborne pathogens, as well as for identifying new pathogens and determining the origin of outbreaks. Metabarcoding and shotgun metagenomics still suffer from low precision, but this issue can be limited by carefully choosing primers and bioinformatic algorithms. Advances in bioinformatics will greatly accelerate the use of metagenomics to address critical aspects related to the detection and surveillance of plant pathogens in plant material and foodstuffs.

Physical drivers facilitating a toxigenic cyanobacterial bloom in a major Great Lakes tributary

The Maumee River is the primary source for nutrients fueling seasonal Microcystis-dominated blooms in western Lake Erie's open waters though such blooms in the river are infrequent. The river also serves as source water for multiple public water systems and a large food services facility in northwest Ohio, USA. On 20 September 2017, an unprecedented bloom was reported in the Maumee River estuary within the Toledo metropolitan area, which triggered a recreational water advisory. Here we (1) explore physical drivers likely contributing to the bloom's occurrence, and (2) describe the toxin concentration and bacterioplankton taxonomic composition. A historical analysis using ten-years of seasonal river discharge, water level, and local wind data identified two instances when high-retention conditions occurred over ≥10 days in the Maumee River estuary: in 2016 and during the 2017 bloom. Observation by remote sensing imagery supported the advection of cyanobacterial cells into the estuary from the lake during 2017 and the lack of an estuary bloom in 2016 due to a weak cyanobacterial bloom in the lake. A rapid-response survey during the 2017 bloom determined levels of the cyanotoxins, specifically microcystins, in excess of recreational contact limits at sites within the lower 20 km of the river while amplicon sequencing found these sites were dominated by Microcystis. These results highlight the need to broaden our understanding of physical drivers of cyanobacterial blooms within the interface between riverine and lacustrine systems, particularly as such blooms are expected to become more prominent in response to a changing climate.

Spatiotemporal Distribution of the Environmental Microbiota in Food Processing Plants as Impacted by Cleaning and Sanitizing Procedures: the Case of Slaughterhouses and Gaseous Ozone

Microbial complexity and contamination levels in food processing plants heavily impact the final product fate and are mainly controlled by proper environmental cleaning and sanitizing. Among the emerging disinfection technologies, ozonation is considered an effective strategy to improve the ordinary cleaning and sanitizing of slaughterhouses. However, its effects on contamination levels and environmental microbiota still need to be understood. For this purpose, we monitored the changes in microbiota composition in different slaughterhouse environments during the phases of cleaning/sanitizing and ozonation at 40, 20, or 4 ppm. Overall, the meat processing plant microbiota differed significantly between secondary processing rooms and deboning rooms, with a greater presence of psychrotrophic taxa in secondary processing rooms because of their lower temperatures. Cleaning/sanitizing procedures significantly reduced the contamination levels and in parallel increased the number of detectable operational taxonomic units (OTUs), by removing the masking effect of the most abundant human/animal-derived OTUs, which belonged to the phylum Firmicutes Subsequently, ozonation at 40 or 20 ppm effectively decreased the remaining viable bacterial populations. However, we could observe selective ozone-mediated inactivation of psychrotrophic bacteria only in the secondary processing rooms. There, the Brochothrix and Pseudomonas abundances and their viable counts were significantly affected by 40 or 20 ppm of ozone, while more ubiquitous genera like Staphylococcus showed a remarkable resistance to the same treatments. This study showed the effectiveness of highly concentrated gaseous ozone as an adjunct sanitizing method that can minimize cross-contamination and so extend the meat shelf life.IMPORTANCE Our in situ survey demonstrates that RNA-based sequencing of 16S rRNA amplicons is a reliable approach to qualitatively probe, at high taxonomic resolution, the changes triggered by new and existing cleaning/sanitizing strategies in the environmental microbiota in human-built environments. This approach could soon represent a fast tool to clearly define which routine sanitizing interventions are more suitable for a specific food processing environment, thus limiting the costs of special cleaning interventions and potential product loss.

Opioid agonist and antagonist use and the gut microbiota: associations among people in addiction treatment

Murine models suggest that opioids alter the gut microbiota, which may impact opioid tolerance and psychopathology. We examined how gut microbiota characteristics related to use of opioid agonists and antagonists among people receiving outpatient addiction treatment. Patients (n = 46) collected stool samples and were grouped by use of opioid agonists (heroin, prescription opioids), antagonists (naltrexone), agonist-antagonist combinations (buprenorphine-naloxone), or neither agonists nor antagonists within the month before enrollment. We sequenced the V4 region of the 16S rRNA gene using Illumina MiSeq to examine how alpha diversity, enterotypes, and relative abundance of bacterial genera varied by opioid agonist and antagonist exposures. Compared to 31 participants who used neither agonists nor antagonists, 5 participants who used opioid agonists (without antagonists) had lower microbiota diversity, Bacteroides enterotypes, and lower relative abundance of Roseburia, a butyrate producing genus, and Bilophila, a bile acid metabolizing genus. There were no differences in gut microbiota features between those using agonist + antagonists (n = 4), antagonists only (n = 6), and neither agonists nor antagonists. Similar to murine morphine exposure models, opioid agonist use was associated with lower microbiota diversity. Lower abundance of Roseburia and Bilophila may relate to the gut inflammation/permeability and dysregulated bile acid metabolism observed in opioid-exposed mice.

Dissolved Microcystin Release Coincident with Lysis of a Bloom Dominated by Microcystis spp. in Western Lake Erie Attributed to a Novel Cyanophage

Western Lake Erie (Laurentian Great Lakes) is prone to annual cyanobacterial harmful algal blooms (cHABs) dominated by Microcystis spp. that often yield microcystin toxin concentrations exceeding the federal EPA recreational contact advisory of 8 μg liter-1 In August 2014, microcystin levels were detected in finished drinking water above the World Health Organization 1.0 μg liter-1 threshold for consumption, leading to a 2-day disruption in the supply of drinking water for >400,000 residents of Toledo, Ohio (USA). Subsequent metatranscriptomic analysis of the 2014 bloom event provided evidence that release of toxin into the water supply was likely caused by cyanophage lysis that transformed a portion of the intracellular microcystin pool into the dissolved fraction, rendering it more difficult to eliminate during treatment. In August 2019, a similar increase in dissolved microcystins at the Toledo water intake was coincident with a viral lytic event caused by a phage consortium different in composition from what was detected following the 2014 Toledo water crisis. The most abundant viral sequence in metagenomic data sets was a scaffold from a putative member of the Siphoviridae, distinct from the Ma-LMM01-like Myoviridae that are typically documented to occur in western Lake Erie. This study provides further evidence that viral activity in western Lake Erie plays a significant role in transformation of microcystins from the particulate to the dissolved fraction and therefore requires monitoring efforts from local water treatment plants. Additionally, identification of multiple lytic cyanophages will enable the development of a quantitative PCR toolbox to assess viral activity during cHABs.IMPORTANCE Viral attack on cHABs may contribute to changes in community composition during blooms, as well as bloom decline, yet loss of bloom biomass does not eliminate the threat of cHAB toxicity. Rather, it may increase risks to the public by delivering a pool of dissolved toxin directly into water treatment utilities when the dominating Microcystis spp. are capable of producing microcystins. Detecting, characterizing, and quantifying the major cyanophages involved in lytic events will assist water treatment plant operators in making rapid decisions regarding the pool of microcystins entering the plant and the corresponding best practices to neutralize the toxin.

Microbial communities of the Laurentian Great Lakes reflect connectivity and local biogeochemistry

The Laurentian Great Lakes are a vast, interconnected freshwater system spanning strong physicochemical gradients, thus constituting a powerful natural laboratory for addressing fundamental questions about microbial ecology and evolution. We present a comparative analysis of pelagic microbial communities across all five Laurentian Great Lakes, focusing on Bacterial and Archaeal picoplankton characterized via 16S rRNA amplicon sequencing. We collected samples throughout the water column from the major basins of each lake in spring and summer over 2 years. Two oligotypes, classified as LD12 (Alphaproteobacteria) and acI‐B1 (Actinobacteria), were among the most abundant in every sample. At the same time, microbial communities showed distinct patterns with depth during summer stratification. Deep hypolimnion samples were frequently dominated by a Chloroflexi oligotype that reached up to 19% relative abundance. Stratified surface communities differed between the colder, less productive upper lakes (Superior, Michigan, Huron) and warmer, more productive lower lakes (Erie, Ontario), in part due to an Actinobacteria oligotype (acI‐C2) that averaged 7.7% of sequences in the lower lakes but <0.2% in the upper lakes. Together, our findings suggest that both hydrologic connectivity and local selective pressures shape microbial communities in the Great Lakes and establish a framework for future investigations.

Effects of Habitat Partitioning on the Distribution of Bacterioplankton in Deep Lakes

In deep lakes, many investigations highlighted the existence of exclusive groups of bacteria adapted to deep oxygenated and hypoxic and anoxic hypolimnia. Nevertheless, the extent of bacterial strain diversity has been much less scrutinized. This aspect is essential for an unbiased estimation of genetic variation, biodiversity, and population structure, which are essential for studying important research questions such as biogeographical patterns, temporal and spatial variability and the environmental factors affecting this variability. This study investigated the bacterioplankton community in the epilimnetic layers and in the oxygenated and hypoxic/anoxic hypolimnia of five large and deep lakes located at the southern border of the Alps using high throughput sequencing (HTS) analyses (16S rDNA) and identification of amplicon sequence variants (ASVs) resolving reads differing by as little as one nucleotide. The study sites, which included two oligomictic (Garda and Como) and three meromictic lakes (Iseo, Lugano, and Idro) with maximum depths spanning from 124 to 410 m, were chosen among large lakes to represent an oxic-hypoxic gradient. The analyses showed the existence of several unique ASVs in the three layers of the five lakes. In the case of cyanobacteria, this confirmed previous analyses made at the level of strains or based on oligotyping methods. As expected, the communities in the hypoxic/anoxic monimolimnia showed a strong differentiation from the oxygenated layer, with the exclusive presence in single lakes of several unique ASVs. In the meromictic lakes, results supported the hypothesis that the formation of isolated monimolimnia sustained the development of highly diversified bacterial communities through ecological selection, leading to the establishment of distinctive biodiversity zones. The genera identified in these layers are well-known to activate a wide range of redox reactions at low O₂ conditions. As inferred from 16S rDNA data, the highly diversified and coupled processes sustained by the monimolimnetic microbiota are essential ecosystem services that enhance mineralization of organic matter and formation of reduced compounds, and also abatement of undesirable greenhouse gasses.

Microbial Similarity between Students in a Common Dormitory Environment Reveals the Forensic Potential of Individual Microbial Signatures

Humans leave behind a microbial trail, regardless of intention. This may allow for the identification of individuals based on the “microbial signatures” they shed in built environments. In a shared living environment, these trails intersect, and through interaction with common surfaces may become homogenized, potentially confounding our ability to link individuals to their associated microbiota. We sought to understand the factors that influence the mixing of individual signatures and how best to process sequencing data to best tease apart these signatures.

The microbiota of the built environment is an amalgamation of both human and environmental sources. While human sources have been examined within single-family households or in public environments, it is unclear what effect a large number of cohabitating people have on the microbial communities of their shared environment. We sampled the public and private spaces of a college dormitory, disentangling individual microbial signatures and their impact on the microbiota of common spaces. We compared multiple methods for marker gene sequence clustering and found that minimum entropy decomposition (MED) was best able to distinguish between the microbial signatures of different individuals and was able to uncover more discriminative taxa across all taxonomic groups. Further, weighted UniFrac- and random forest-based graph analyses uncovered two distinct spheres of hand- or shoe-associated samples. Using graph-based clustering, we identified spheres of interaction and found that connection between these clusters was enriched for hands, implicating them as a primary means of transmission. In contrast, shoe-associated samples were found to be freely interacting, with individual shoes more connected to each other than to the floors they interact with. Individual interactions were highly dynamic, with groups of samples originating from individuals clustering freely with samples from other individuals, while all floor and shoe samples consistently clustered together.

DAIRYdb: a manually curated reference database for improved taxonomy annotation of 16S rRNA gene sequences from dairy products

Background

Reads assignment to taxonomic units is a key step in microbiome analysis pipelines. To date, accurate taxonomy annotation of 16S reads, particularly at species rank, is still challenging due to the short size of read sequences and differently curated classification databases. The close phylogenetic relationship between species encountered in dairy products, however, makes it crucial to annotate species accurately to achieve sufficient phylogenetic resolution for further downstream ecological studies or for food diagnostics. Curated databases dedicated to the environment of interest are expected to improve the accuracy and resolution of taxonomy annotation.

Results

We provide a manually curated database composed of 10’290 full-length 16S rRNA gene sequences from prokaryotes tailored for dairy products analysis (https://github.com/marcomeola/DAIRYdb).

The performance of the DAIRYdb was compared with the universal databases Silva, LTP, RDP and Greengenes. The DAIRYdb significantly outperformed all other databases independently of the classification algorithm by enabling higher accurate taxonomy annotation down to the species rank. The DAIRYdb accurately annotates over 90% of the sequences of either single or paired hypervariable regions automatically.

The manually curated DAIRYdb strongly improves taxonomic annotation accuracy for microbiome studies in dairy environments. The DAIRYdb is a practical solution that enables automatization of this key step, thus facilitating the routine application of NGS microbiome analyses for microbial ecology studies and diagnostics in dairy products.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5914-8) contains supplementary material, which is available to authorized users.

Understanding and overcoming the pitfalls and biases of next-generation sequencing (NGS) methods for use in the routine clinical microbiological diagnostic laboratory

Recent advancements in next-generation sequencing (NGS) have provided the foundation for modern studies into the composition of microbial communities. The use of these NGS methods allows for the detection and identification of (‘difficult-to-culture’) microorganisms using a culture-independent strategy. In the field of routine clinical diagnostics however, the application of NGS is currently limited to microbial strain typing for epidemiological purposes only, even though the implementation of NGS for microbial community analysis may yield clinically important information. This lack of NGS implementation is due to many different factors, including issues relating to NGS method standardization and result reproducibility. In this review article, the authors provide a general introduction to the most widely used NGS methods currently available (i.e., targeted amplicon sequencing and shotgun metagenomics) and the strengths and weaknesses of each method is discussed. The focus of the publication then shifts toward 16S rRNA gene NGS methods, which are currently the most cost-effective and widely used NGS methods for research purposes, and are therefore more likely to be successfully implemented into routine clinical diagnostics in the short term. In this respect, the experimental pitfalls and biases created at each step of the 16S rRNA gene NGS workflow are explained, as well as their potential solutions. Finally, a novel diagnostic microbiota profiling platform (‘MYcrobiota’) is introduced, which was developed by the authors by taking into consideration the pitfalls, biases, and solutions explained in this article. The development of the MYcrobiota, and future NGS methodologies, will help pave the way toward the successful implementation of NGS methodologies into routine clinical diagnostics.

The respiratory microbiome and susceptibility to influenza virus infection

Influenza is a major cause of morbidity and mortality worldwide. However, vaccine effectiveness has been low to moderate in recent years and vaccine coverage remains low, especially in low- and middle-income countries. Supplementary methods of prevention should be explored to reduce the high burden of influenza. A potential target is the respiratory tract microbiome, complex microbial communities which envelop the respiratory epithelium and play an important role in shaping host immunity. Using a household transmission study, we examined whether the nose/throat microbiota was associated with influenza susceptibility among participants exposed to influenza virus in the household. Further, we characterized changes in the nose/throat microbiota to explore whether community stability was influenced by influenza virus infection. Using a generalized linear mixed effects model, we found a nasal/oropharyngeal community state type (CST) associated with decreased susceptibility to influenza. The CST was rare and transitory among young children but a prevalent and stable CST among adults. Using boosting and linear mixed effects models, we found associations between the nose/throat microbiota and influenza also existed at the taxa level, specifically with the relative abundance of Alloprevotella, Prevotella, and Bacteroides oligotypes. We found high rates of change between bacterial community states among both secondary cases and household contacts who were not infected during follow up. Further work is needed to separate the effect of influenza virus infection from the considerable short-term changes that occur even in the absence of virus. Lastly, age was strongly associated with susceptibility to influenza and the nose/throat bacterial community structure. Although additional studies are needed to determine causality, our results suggest the nose/throat microbiome may be a potential target for reducing the burden of influenza.

Spatial-temporal survey of Microcystis oligopeptide chemotypes in reservoirs with dissimilar waterbody features and their relation to genetic variation

The ability of cyanobacteria to produce toxins and other secondary metabolites is patchily distributed in natural populations, enabling the use of cellular oligopeptide compositions as markers to classify strains into ecologically-relevant chemotypical subpopulations. The composition and spatiotemporal distribution of Microcystis chemotypes within and among waterbodies was studied at different time scales by analyzing (i) Microcystis strains isolated between 1998 and 2007 from different Spanish reservoirs and (ii) individual Microcystis aeruginosa colonies collected from pelagic and littoral habitats in Valmayor reservoir (Spain) during a bloom. No agreement between chemotypes and both morphotypes and genotypes (based on cpcBA-IGS, 16S-23S rRNA ITS and mcyB genes) was found, suggesting that oligopeptide profiles in individual strains evolve independently across morphospecies and phylogenetic genotypes, and that the diversity of microcystin variants produced cannot be explained by mcyB gene variations alone. The presence of identical chemotypes in spatially-distant reservoirs with dissimilar trophic state, lithology or depth indicate that waterbody characteristics and geographical boundaries weakly affect chemotype composition and distribution. At smaller spatiotemporal scales (i.e. during bloom), M. aeruginosa populations showed high number of chemotypes, as well as marked differences in chemotype composition and relative abundance among the littoral and pelagic habitats. This indicates that the factors influencing chemotype composition, relative abundance and dynamics operate at short spatial and temporal scales, and supports emerging hypotheses about interactions with antagonistic microorganisms as possible drivers for widespread chemical polymorphisms in cyanobacteria.

Monitoring of microbial dynamics in a drinking water distribution system using the culture-free, user-friendly, MYcrobiota platform

Drinking water utilities currently rely on a range of microbiological detection techniques to evaluate the quality of their drinking water (DW). However, microbiota profiling using culture-free 16S rRNA gene next-generation sequencing (NGS) provides an opportunity for improved monitoring of the microbial ecology and quality of DW. Here, we evaluated the utility of a previously validated microbiota profiling platform (MYcrobiota) to investigate the microbial dynamics of a full-scale, non-chlorinated DW distribution system (DWDS). In contrast to conventional methods, we observed spatial and temporal bacterial genus changes (expressed as operational taxonomic units - OTUs) within the DWDS. Further, a small subset of bacterial OTUs dominated with abundances that shifted across the length of the DWDS, and were particularly affected by a post-disinfection step. We also found seasonal variation in OTUs within the DWDS and that many OTUs could not be identified, even though MYcrobiota is specifically designed to reduce potential PCR sequencing artefacts. This suggests that our current knowledge about the microbial ecology of DW communities is limited. Our findings demonstrate that the user-friendly MYcrobiota platform facilitates culture-free, standardized microbial dynamics monitoring and has the capacity to facilitate the introduction of microbiota profiling into the management of drinking water quality.

Temporal and spatial distribution of Microcystis biomass and genotype in bloom areas of Lake Taihu

Cyanobacterial blooms as a global environmental issue are of public health concern. In this study, we investigated the spatial (10 sites) and temporal (June, August and October) variations in: 1) their biomass based on chlorophyll-a (chl-a) concentration, 2) their toxic genotype based on gene copy ratio of mcyJ to cpcBA, and 3) their cpcBA genotype composition of Microcystis during cyanobacterial bloom in Lake Taihu. While spatial-temporal variations were found in chl-a and mcyJ/cpcBA ratio, only spatial variation was observed in cpcBA genotype composition. Samples from northwestern part had a higher chl-a, but mcyJ/cpcBA ratio didn't vary among the sites. High chl-a was observed in August, while mcyJ/cpcBA ratio and genotypic richness increased with time. The spatial variations in chl-a and mcyJ/cpcBA ratio and temporal variation in cpcBA genotype were correlated negatively with dissolved N and positively with dissolved P. Spatial distribution of Microcystis biomass was positively correlated with nitrite and P excluding October, but no correlation was found for spatial distribution of mcyJ/cpcBA ratio and cpcBA genotype. Spatial distribution of toxic and cpcBA genotypes may result from horizontal transport of Microcystis colonies, while spatial variation in Microcystis biomass was probably controlled by both nutrient-mediated growth and horizontal transport of Microcystis. The temporal variation in Microcystis biomass, toxic genotype and cpcBA genotype composition were related to nutrient levels, but cause-and-effect relationships require further study.

Diversity and Cyclical Seasonal Transitions in the Bacterial Community in a Large and Deep Perialpine Lake

High-throughput sequencing (HTS) was used to analyze the seasonal variations in the bacterioplankton community composition (BCC) in the euphotic layer of a large and deep lake south of the Alps (Lake Garda). The BCC was analyzed throughout two annual cycles by monthly samplings using the amplification and sequencing of the V3-V4 hypervariable region of the 16S rRNA gene by the MiSeq Illumina platform. The dominant and most diverse bacterioplankton phyla were among the more frequently reported in freshwater ecosystems, including the Proteobacteria, Cyanobacteria, Bacteroidetes, Verrucomicrobia, Actinobacteria, and Planctomycetes. As a distinctive feature, the development of the BCC showed a cyclical temporal pattern in the two analyzed years and throughout the euphotic layer. The recurring temporal development was controlled by the strong seasonality in water temperature and thermal stratification, and by cyclical temporal changes in nutrients and, possibly, by the remarkable annual cyclical development of cyanobacteria and eukaryotic phytoplankton hosting bacterioplankton that characterizes Lake Garda. Further downstream analyses of operational taxonomic units associated to cyanobacteria allowed confirming the presence of the most abundant taxa previously identified by microscopy and/or phylogenetic analyses, as well as the presence of other small Synechococcales/Chroococcales and rare Nostocales never identified so far in the deep lakes south of the Alps. The implications of the high diversity and strong seasonality are relevant, opening perspectives for the definition of common and discriminating patterns characterizing the temporal and spatial distribution in the BCC, and for the application of the new sequencing technologies in the monitoring of water quality in large and deep lakes.

Characterization of Microcystis (Cyanobacteria) Genotypes Based on the Internal Transcribed Spacer Region of rRNA by Next-Generation Sequencing

Microcystis is one of the most common and dominant bloom-forming cyanobacteria in freshwater worldwide. The method for genotype detection based on traditional molecular cloning is expensive and time consuming and generates a limited number of sequences. In this study, a high-throughput sequencing (HTS) method was developed to detect the internal transcribed spacer (ITS) regions between 16S and 23S rRNA region of Microcystis populations along a typical water system in Yuqiao Reservoir-Haihe River in Tianjin, northern China. A total of 629,341 reads were obtained and clustered into 2005 operational taxonomic units (OTUs). Analysis of alpha diversity indices showed that the Haihe River is more diverse than Yuqiao Reservoir. In general, the two water areas exhibit a clear differentiation pattern in OTU abundance, sharing genotypes from a small part of Yuqiao Reservoir with those in the Haihe River. Phylogenetic analysis further indicated the possible flexible evolution of Microcystis genotypes occurring in the research areas. This study provides the first exhaustive description of HTS method for detection of ITS region to evaluate Microcystis intra-species diversity and relationship.

Identification and Removal of Contaminant Sequences From Ribosomal Gene Databases: Lessons From the Census of Deep Life

Earth’s subsurface environment is one of the largest, yet least studied, biomes on Earth, and many questions remain regarding what microorganisms are indigenous to the subsurface. Through the activity of the Census of Deep Life (CoDL) and the Deep Carbon Observatory, an open access 16S ribosomal RNA gene sequence database from diverse subsurface environments has been compiled. However, due to low quantities of biomass in the deep subsurface, the potential for incorporation of contaminants from reagents used during sample collection, processing, and/or sequencing is high. Thus, to understand the ecology of subsurface microorganisms (i.e., the distribution, richness, or survival), it is necessary to minimize, identify, and remove contaminant sequences that will skew the relative abundances of all taxa in the sample. In this meta-analysis, we identify putative contaminants associated with the CoDL dataset, recommend best practices for removing contaminants from samples, and propose a series of best practices for subsurface microbiology sampling. The most abundant putative contaminant genera observed, independent of evenness across samples, were Propionibacterium, Aquabacterium, Ralstonia, and Acinetobacter. While the top five most frequently observed genera were Pseudomonas, Propionibacterium, Acinetobacter, Ralstonia, and Sphingomonas. The majority of the most frequently observed genera (high evenness) were associated with reagent or potential human contamination. Additionally, in DNA extraction blanks, we observed potential archaeal contaminants, including methanogens, which have not been discussed in previous contamination studies. Such contaminants would directly affect the interpretation of subsurface molecular studies, as methanogenesis is an important subsurface biogeochemical process. Utilizing previously identified contaminant genera, we found that ∼27% of the total dataset were identified as contaminant sequences that likely originate from DNA extraction and DNA cleanup methods. Thus, controls must be taken at every step of the collection and processing procedure when working with low biomass environments such as, but not limited to, portions of Earth’s deep subsurface. Taken together, we stress that the CoDL dataset is an incredible resource for the broader research community interested in subsurface life, and steps to remove contamination derived sequences must be taken prior to using this dataset.

Niche Separation Increases With Genetic Distance Among Bloom-Forming Cyanobacteria

Bacterial communities are composed of distinct groups of potentially interacting lineages, each thought to occupy a distinct ecological niche. It remains unclear, however, how quickly niche preference evolves and whether more closely related lineages are more likely to share ecological niches. We addressed these questions by following the dynamics of two bloom-forming cyanobacterial genera over an 8-year time-course in Lake Champlain, Canada, using 16S amplicon sequencing and measurements of several environmental parameters. The two genera, Microcystis (M) and Dolichospermum (D), are frequently observed simultaneously during bloom events and thus have partially overlapping niches. However, the extent of their niche overlap is debated, and it is also unclear to what extent niche partitioning occurs among strains within each genus. To identify strains within each genus, we applied minimum entropy decomposition (MED) to 16S rRNA gene sequences. We confirmed that at a genus level, M and D have different preferences for nitrogen and phosphorus concentrations. Within each genus, we also identified strains differentially associated with temperature, precipitation, and concentrations of nutrients and toxins. In general, niche similarity between strains (as measured by co-occurrence over time) declined with genetic distance. This pattern is consistent with habitat filtering - in which closely related taxa are ecologically similar, and therefore tend to co-occur under similar environmental conditions. In contrast with this general pattern, similarity in certain niche dimensions (notably particulate nitrogen and phosphorus) did not decline linearly with genetic distance, and instead showed a complex polynomial relationship. This observation suggests the importance of processes other than habitat filtering - such as competition between closely related taxa, or convergent trait evolution in distantly related taxa - in shaping particular traits in microbial communities.

The oral microbiome

The human microbiome is receiving a great deal of attention as its role in health and disease becomes ever more apparent. The oral microbiome, perhaps due to the ease with which we can obtain samples, is arguably the most well-studied human microbiome to date. It is obvious, however, that we have only just begun to scratch the surface of the complex bacterial and bacterial–host interactions within this complex community. Here, we describe the factors which are known to influence the development of the seemingly globally conserved, core, oral microbiome and those which are likely to be responsible for the observed differences at the individual level. We discuss the paradoxical situation of maintaining a stable core microbiome which is at the same time incredibly resilient and adaptable to many different stresses encountered in the open environment of the oral cavity. Finally, we explore the interactions of the oral microbiome with the host and discuss the interactions underlying human health and disease.

Ubiquity and quantitative significance of bacterioplankton lineages inhabiting the oxygenated hypolimnion of deep freshwater lakes

The oxygenated hypolimnion accounts for a volumetrically significant part of the global freshwater systems. Previous studies have proposed the presence of hypolimnion-specific bacterioplankton lineages that are distinct from those inhabiting the epilimnion. To date, however, no consensus exists regarding their ubiquity and abundance, which is necessary to evaluate their ecological importance. The present study investigated the bacterioplankton community in the oxygenated hypolimnia of 10 deep freshwater lakes. Despite the broad geochemical characteristics of the lakes, 16S rRNA gene sequencing demonstrated that the communities in the oxygenated hypolimnia were distinct from those in the epilimnia and identified several predominant lineages inhabiting multiple lakes. Catalyzed reporter deposition fluorescence in situ hybridization revealed that abundant hypolimnion-specific lineages, CL500-11 (Chloroflexi), CL500-3, CL500-37, CL500-15 (Planctomycetes) and Marine Group I (Thaumarchaeota), together accounted for 1.5-32.9% of all bacterioplankton in the hypolimnion of the lakes. Furthermore, an analysis of single-nucleotide variation in the partial 16S rRNA gene sequence (oligotyping) suggested the presence of different sub-populations between lakes and water layers among the lineages occurring in the entire water layer (for example, acI-B1 and acI-A7). Collectively, these results provide the first comprehensive overview of the bacterioplankton community in the oxygenated hypolimnion of deep freshwater lakes.

Variation of Oxygenation Conditions on a Hydrocarbonoclastic Microbial Community Reveals Alcanivorax and Cycloclasticus Ecotypes

Deciphering the ecology of marine obligate hydrocarbonoclastic bacteria (MOHCB) is of crucial importance for understanding their success in occupying distinct niches in hydrocarbon-contaminated marine environments after oil spills. In marine coastal sediments, MOHCB are particularly subjected to extreme fluctuating conditions due to redox oscillations several times a day as a result of mechanical (tide, waves and currents) and biological (bioturbation) reworking of the sediment. The adaptation of MOHCB to the redox oscillations was investigated by an experimental ecology approach, subjecting a hydrocarbon-degrading microbial community to contrasting oxygenation regimes including permanent anoxic conditions, anoxic/oxic oscillations and permanent oxic conditions. The most ubiquitous MOHCB, Alcanivorax and Cycloclasticus, showed different behaviors, especially under anoxic/oxic oscillation conditions, which were more favorable for Alcanivorax than for Cycloclasticus. The micro-diversity of 16S rRNA gene transcripts from these genera revealed specific ecotypes for different oxygenation conditions and their dynamics. It is likely that such ecotypes allow the colonization of distinct ecological niches that may explain the success of Alcanivorax and Cycloclasticus in hydrocarbon-contaminated coastal sediments during oil-spills.

Exact sequence variants should replace operational taxonomic units in marker-gene data analysis

Recent advances have made it possible to analyze high-throughput marker-gene sequencing data without resorting to the customary construction of molecular operational taxonomic units (OTUs): clusters of sequencing reads that differ by less than a fixed dissimilarity threshold. New methods control errors sufficiently such that amplicon sequence variants (ASVs) can be resolved exactly, down to the level of single-nucleotide differences over the sequenced gene region. The benefits of finer resolution are immediately apparent, and arguments for ASV methods have focused on their improved resolution. Less obvious, but we believe more important, are the broad benefits that derive from the status of ASVs as consistent labels with intrinsic biological meaning identified independently from a reference database. Here we discuss how these features grant ASVs the combined advantages of closed-reference OTUs—including computational costs that scale linearly with study size, simple merging between independently processed data sets, and forward prediction—and of de novo OTUs—including accurate measurement of diversity and applicability to communities lacking deep coverage in reference databases. We argue that the improvements in reusability, reproducibility and comprehensiveness are sufficiently great that ASVs should replace OTUs as the standard unit of marker-gene analysis and reporting.