Synthetic and natural rubber associated chemicals drive functional and structural changes as well as adaptations to antibiotics in in vitro marine microbiomes.

The leaching of additives from plastics and elastomers (rubbers) has raised concerns due to their potential negative impacts on the environment and the development of antibiotic resistance. In this study, we investigated the effects of chemicals extracted from two types of rubber on microbiomes derived from a benthic sea urchin and two pelagic fish species. Additionally, we examined whether bacterial communities preconditioned with rubber-associated chemicals displayed adaptations to antibiotics. At the highest tested concentrations of chemicals, we observed reduced maximum growth rates and yields, prolonged lag phases, and increased alpha diversity. While the effects on alpha and beta diversity were not always conclusive, several bacterial genera were significantly influenced by chemicals from the two rubber sources. Subsequent exposure of sea urchin microbiomes preconditioned with rubber chemicals to the antibiotic ciprofloxacin resulted in decreased maximum growth rates. This indicates a more sensitive microbiome to ciprofloxacin when preconditioned with rubber chemicals. Although no significant interaction effects between rubber chemicals and ciprofloxacin exposure were observed in bacterial alpha and beta diversity, we observed log-fold changes in two bacterial genera in response to ciprofloxacin exposure. These findings highlight the structural and functional alterations in microbiomes originating from various marine species when exposed to rubber-associated chemicals and underscore the potential risks posed to marine life.

Bacterial bioindicators enable biological status classification along the continental Danube river.

Despite the importance of bacteria in aquatic ecosystems and their predictable diversity patterns across space and time, biomonitoring tools for status assessment relying on these organisms are widely lacking. This is partly due to insufficient data and models to identify reliable microbial predictors. Here, we show metabarcoding in combination with multivariate statistics and machine learning allows to identify bacterial bioindicators for existing biological status classification systems. Bacterial beta-diversity dynamics follow environmental gradients and the observed associations highlight potential bioindicators for ecological outcomes. Spatio-temporal links spanning the microbial communities along the river allow accurate prediction of downstream biological status from upstream information. Network analysis on amplicon sequence veariants identify as good indicators genera Fluviicola, Acinetobacter, Flavobacterium, and Rhodoluna, and reveal informational redundancy among taxa, which coincides with taxonomic relatedness. The redundancy among bacterial bioindicators reveals mutually exclusive taxa, which allow accurate biological status modeling using as few as 2-3 amplicon sequence variants. As such our models show that using a few bacterial amplicon sequence variants from globally distributed genera allows for biological status assessment along river systems.

Occupational exposure to inhalable pathogenic microorganisms in waste sorting.

This study assessed microorganisms in personal inhalable work air samples aiming to identify potential human pathogens, and correlate exposure to adverse health outcomes in waste workers. Full-shift personal exposure was measured in six different waste sorting plants. Microbial concentrations in inhalable air samples were analysed using MALDI-TOF MS for cultivable, and next generation sequencing (NGS) for non-cultivable microorganisms. Concentrations of bacterial and fungal CFUs varied substantially within and between waste sorting plants, ranging from no identifiable organisms to a maximum concentration in the order of 105 CFU/m3. Bacillus and Staphylococcus were among the most abundant bacterial genera, whilst fungal genera were dominated by Aspergillus and Penicillium. Approximately 15% of all identified species were human pathogens classified in risk group 2, whereas 7% belonged to risk group 1. Furthermore, significant correlations between concentrations of fungi in risk group 1 and self-reported adverse symptoms, such as wheezing were identified in exposed workers. The combination of culture-based methods and NGS facilitated the investigation of infectious microbial species with potential pathophysiological properties as well as non-infectious biological agents in inhalable work air samples and thereby contributed to the risk assessment of occupational exposure in waste sorting.

Trajectories of freshwater microbial genomics and greenhouse gas saturation upon glacial retreat.

Due to climate warming, ice sheets around the world are losing mass, contributing to changes across terrestrial landscapes on decadal time spans. However, landscape repercussions on climate are poorly constrained mostly due to limited knowledge on microbial responses to deglaciation. Here, we reveal the genomic succession from chemolithotrophy to photo- and heterotrophy and increases in methane supersaturation in freshwater lakes upon glacial retreat. Arctic lakes at Svalbard also revealed strong microbial signatures form nutrient fertilization by birds. Although methanotrophs were present and increased along lake chronosequences, methane consumption rates were low even in supersaturated systems. Nitrous oxide oversaturation and genomic information suggest active nitrogen cycling across the entire deglaciated landscape, and in the high Arctic, increasing bird populations serve as major modulators at many sites. Our findings show diverse microbial succession patterns, and trajectories in carbon and nitrogen cycle processes representing a positive feedback loop of deglaciation on climate warming.

Drivers and variability of CO2:O2 saturation along a gradient from boreal to Arctic lakes.

Lakes are significant players for the global climate since they sequester terrestrially derived dissolved organic carbon (DOC), and emit greenhouse gases like CO2 to the atmosphere. However, the differences in environmental drivers of CO2 concentrations are not well constrained along latitudinal and thus climate gradients. Our aim here is to provide a better understanding of net heterotrophy and gas balance at the catchment scale in a set of boreal, sub-Arctic and high-Arctic lakes. We assessed water chemistry and concentrations of dissolved O2 and CO2, as well as the CO2:O2 ratio in three groups of lakes separated by steps of approximately 10 degrees latitude in South-Eastern Norway (near 60° N), sub-Arctic lakes in the northernmost part of the Norwegian mainland (near 70° N) and high-Arctic lakes on Svalbard (near 80° N). Across all regions, CO2 saturation levels varied more (6-1374%) than O2 saturation levels (85-148%) and hence CO2 saturation governed the CO2:O2 ratio. The boreal lakes were generally undersaturated with O2, while the sub-Arctic and high-Arctic lakes ranged from O2 saturated to oversaturated. Regardless of location, the majority of the lakes were CO2 supersaturated. In the boreal lakes the CO2:O2 ratio was mainly related to DOC concentration, in contrast to the sub-Arctic and high-Arctic localities, where conductivity was the major statistical determinant. While the southern part is dominated by granitic and metamorphic bedrock, the sub-Arctic sites are scattered across a range of granitic to sedimentary bed rocks, and the majority of the high-Arctic lakes are situated on limestone, resulting in contrasting lake alkalinities between the regions. DOC dependency of the CO2:O2 ratio in the boreal region together with low alkalinity suggests that in-lake heterotrophic respiration was a major source of lake CO2. Contrastingly, the conductivity dependency indicates that CO2 saturation in the sub-Arctic and high-Arctic lakes was to a large part explained by DIC input from catchment respiration and carbonate weathering.

Detecting aquatic pathogens with field-compatible dried qPCR assays.

Field-ready qPCR assays with extended shelf-life support monitoring programs for emerging aquatic pathogens and enable quick conservation and management decisions. Here, we developed, validated, and tested the shelf-life of qPCR assays targeting Gyrodactylus salaris and Aphanomyces astaci with lyophilization and air-drying.

Disturbance-based management of ecosystem services and disservices in partial nitritation-anammox biofilms.

The resistance and resilience provided by functional redundancy, a common feature of microbial communities, is not always advantageous. An example is nitrite oxidation in partial nitritation-anammox (PNA) reactors designed for nitrogen removal in wastewater treatment, where suppression of nitrite oxidizers like Nitrospira is sought. In these ecosystems, biofilms provide microhabitats with oxygen gradients, allowing the coexistence of aerobic and anaerobic bacteria. We designed a disturbance experiment where PNA biofilms, treating water from a high-rate activated sludge process, were constantly or intermittently exposed to anaerobic sidestream wastewater, which has been proposed to inhibit nitrite oxidizers. With increasing sidestream exposure we observed decreased abundance, alpha-diversity, functional versatility, and hence functional redundancy, among Nitrospira in the PNA biofilms, while the opposite patterns were observed for anammox bacteria within Brocadia. At the same time, species turnover was observed for aerobic ammonia-oxidizing Nitrosomonas populations. The different exposure regimens were associated with metagenomic assembled genomes of Nitrosomonas, Nitrospira, and Brocadia, encoding genes related to N-cycling, substrate usage, and osmotic stress response, possibly explaining the three different patterns by niche differentiation. These findings imply that disturbances can be used to manage the functional redundancy of biofilm microbiomes in a desirable direction, which should be considered when designing operational strategies for wastewater treatment.

Two different approaches of microbial community structure characterization in riverine epilithic biofilms under multiple stressors conditions: Developing molecular indicators.

Microbial communities are major players in the biogeochemical processes and ecosystem functioning of river networks. Despite their importance in the ecosystem, biomonitoring tools relying on prokaryotes are still lacking. Only a few studies have employed both metabarcoding and quantitative techniques such as catalysed reported deposition fluorescence in situ hybridization (CARD-FISH) to analyse prokaryotic communities of epilithic biofilms in river ecosystems. We intended to investigate the efficacy of both techniques in detecting changes in microbial community structure associated with environmental drivers. We report a significant correlation between the prokaryotic community composition and pH in rivers from two different geographical areas in Norway. Both CARD-FISH and metabarcoding data were following the pattern of the environmental variables, but the main feature distinguishing the community composition was the regional difference itself. Beta-dispersion analyses on both CARD-FISH abundance and metabarcoding data revealed higher accuracy of metabarcoding to differentiate regions and river systems. The CARD-FISH results showed high variability, even for samples within the same river, probably due to some unmeasured microscale ecological variability which we could not resolve. We also present a statistical method, which uses variation coefficient and overall prevalence of taxonomic groups, to detect possible biological indicators among prokaryotes using metabarcoding data. The development of new prokaryotic bioindicators would benefit from both techniques used in this study, but metabarcoding seems to be faster and more reliable than CARD-FISH for large scale bio-assessment.

Multiple thresholds and trajectories of microbial biodiversity predicted across browning gradients by neural networks and decision tree learning.

Ecological association studies often assume monotonicity such as between biodiversity and environmental properties although there is growing evidence that nonmonotonic relations dominate in nature. Here, we apply machine-learning algorithms to reveal the nonmonotonic association between microbial diversity and an anthropogenic-induced large-scale change, the browning of freshwaters, along a longitudinal gradient covering 70 boreal lakes in Scandinavia. Measures of bacterial richness and evenness (alpha-diversity) showed nonmonotonic trends in relation to environmental gradients, peaking at intermediate levels of browning. Depending on the statistical methods, variables indicative for browning could explain 5% of the variance in bacterial community composition (beta-diversity) when applying standard methods assuming monotonic relations and up to 45% with machine-learning methods taking non-monotonicity into account. This non-monotonicity observed at the community level was explained by the complex interchangeable nature of individual taxa responses as shown by a high degree of nonmonotonic responses of individual bacterial sequence variants to browning. Furthermore, the nonmonotonic models provide the position of thresholds and predict alternative bacterial diversity trajectories in boreal freshwater as a result of ongoing climate and land-use changes, which in turn will affect entire ecosystem metabolism and likely greenhouse gas production.

Corrigendum: Phosphorus Availability Promotes Bacterial DOC-Mineralization, but Not Cumulative CO2-Production.

[This corrects the article DOI: 10.3389/fmicb.2020.569879.].

Phosphorus Availability Promotes Bacterial DOC-Mineralization, but Not Cumulative CO2-Production.

The current trend of increasing input of terrestrially derived dissolved organic carbon (DOC) to boreal freshwater systems is causing increased levels of carbon dioxide (CO2) supersaturation and degassing. Phosphorus (P) is often the most limiting nutrient for bacterial growth and would thus be expected to increase overall mineralization rates and CO2 production. However, high carbon (C) to P ratios of terrestrially derived DOC could also cause elevated cell-specific respiration of the excess C in heterotrophic bacteria. Using data from a survey of 75 Scandinavian lakes along an ecosystem gradient of DOC, we estimated in situ CO2 production rates. These rates showed a unimodal response with DOC-specific CO2 production negatively related to DOC:total phosphorus (TP) ratio, and a turning point at 5 mg C L-1, indicating higher DOC turnover rates in productive than in unproductive lakes. To further assess the dependency of bacterial respiration (BR) on DOC and P, we monitored CO2 production in incubations of water with a gradient of DOC crossed with two levels of inorganic P. Finally, we crossed DOC and P with a temperature gradient to test the temperature dependency of respiration rates [as oxygen (O2) consumption]. While total CO2 production seemed to be unaffected by P additions, respiration rates, and growth yields, as estimated by ribosomal gene copy numbers, suggest increased bacterial growth and decreased cell-specific respiration under non-limited P conditions. Respiration rates showed a sigmoid response to increasing DOC availability reaching a plateau at about 20 mg C L-1 of initial DOC concentrations. In addition to these P and DOC level effects, respiration rates responded in a non-monotonic fashion to temperature with an increase in respiration rates by a factor of 2.6 (±0.2) from 15 to 25°C and a decrease above 30°C. The combined results from the survey and experiments highlight DOC as the major determinant of CO2 production in boreal lakes, with P and temperature as significant modulators of respiration kinetics.

Diazotroph Genomes and Their Seasonal Dynamics in a Stratified Humic Bog Lake.

Aquatic N-fixation is generally associated with the growth and mass development of Cyanobacteria in nitrogen-deprived photic zones. However, sequenced genomes and environmental surveys suggest active aquatic N-fixation also by many non-cyanobacterial groups. Here, we revealed the seasonal variation and genomic diversity of potential N-fixers in a humic bog lake using metagenomic data and nif gene clusters analysis. Groups with diazotrophic operons were functionally divergent and included Cholorobi, Geobacter, Desulfobacterales, Methylococcales, and Acidobacteria. In addition to nifH (a gene that encodes the dinitrogenase reductase component of the molybdenum nitrogenase), we also identified sequences corresponding to vanadium and iron-only nitrogenase genes. Within the Chlorobi population, the nitrogenase (nifH) cluster was included in a well-structured retrotransposon. Furthermore, the presence of light-harvesting photosynthesis genes implies that anoxygenic photosynthesis may fuel nitrogen fixation under the prevailing low-irradiance conditions. The presence of rnf genes (related to the expression of H+/Na+-translocating ferredoxin: NAD+ oxidoreductase) in Methylococcales and Desulfobacterales suggests that other energy-generating processes may drive the costly N-fixation in the absence of photosynthesis. The highly reducing environment of the anoxic bottom layer of Trout Bog Lake may thus also provide a suitable niche for active N-fixers and primary producers. While future studies on the activity of these potential N-fixers are needed to clarify their role in freshwater nitrogen cycling, the metagenomic data presented here enabled an initial characterization of previously overlooked diazotrophs in freshwater biomes.

Archaea in boreal Swedish lakes are diverse, dominated by Woesearchaeota and follow deterministic community assembly.

Despite their key role in biogeochemical processes, particularly the methane cycle, archaea are widely underrepresented in molecular surveys because of their lower abundance compared with bacteria and eukaryotes. Here, we use parallel high-resolution small subunit rRNA gene sequencing to explore archaeal diversity in 109 Swedish lakes and correlate archaeal community assembly mechanisms to large-scale latitudinal, climatic (nemoral to arctic) and nutrient (oligotrophic to eutrophic) gradients. Sequencing with universal primers showed the contribution of archaea was on average 0.8% but increased up to 1.5% of the three domains in forest lakes. Archaea-specific sequencing revealed that freshwater archaeal diversity could be partly explained by lake variables associated with nutrient status. Combined with deterministic co-occurrence patterns this finding suggests that ecological drift is overridden by environmental sorting, as well as other deterministic processes such as biogeographic and evolutionary history, leading to lake-specific archaeal biodiversity. Acetoclastic, hydrogenotrophic and methylotrophic methanogens as well as ammonia-oxidizing archaea were frequently detected across the lakes. Archaea-specific sequencing also revealed representatives of Woesearchaeota and other phyla of the DPANN superphylum. This study adds to our understanding of the ecological range of key archaea in freshwaters and links these taxa to hypotheses about processes governing biogeochemical cycles in lakes.

Novel Autotrophic Organisms Contribute Significantly to the Internal Carbon Cycling Potential of a Boreal Lake.

Oxygen-stratified lakes are typical for the boreal zone and also a major source of greenhouse gas emissions in the region. Due to shallow light penetration, restricting the growth of phototrophic organisms, and large allochthonous organic carbon inputs from the catchment area, the lake metabolism is expected to be dominated by heterotrophic organisms. In this study, we test this assumption and show that the potential for autotrophic carbon fixation and internal carbon cycling is high throughout the water column. Further, we show that during the summer stratification carbon fixation can exceed respiration in a boreal lake even below the euphotic zone. Metagenome-assembled genomes and 16S profiling of a vertical transect of the lake revealed multiple organisms in an oxygen-depleted compartment belonging to novel or poorly characterized phyla. Many of these organisms were chemolithotrophic, potentially deriving their energy from reactions related to sulfur, iron, and nitrogen transformations. The community, as well as the functions, was stratified along the redox gradient. The autotrophic potential in the lake metagenome below the oxygenic zone was high, pointing toward a need for revising our concepts of internal carbon cycling in boreal lakes. Further, the importance of chemolithoautotrophy for the internal carbon cycling suggests that many predicted climate change-associated fluctuations in the physical properties of the lake, such as altered mixing patterns, likely have consequences for the whole-lake metabolism even beyond the impact to the phototrophic community.IMPORTANCE Autotrophic organisms at the base of the food web are the only life form capable of turning inorganic carbon into the organic form, facilitating the survival of all other organisms. In certain environments, the autotrophic production is limited by environmental conditions and the food web is supported by external carbon inputs. One such environment is stratified boreal lakes, which are one of the biggest natural sources of greenhouse gas emissions in the boreal region. Thus, carbon cycling in these habitats is of utmost importance for the future climate. Here, we demonstrate a high potential for internal carbon cycling via phototrophic and novel chemolithotrophic organisms in the anoxic, poorly illuminated layers of a boreal lake. Our results significantly increase our knowledge on the microbial communities and their metabolic potential in oxygen-depleted freshwaters and help to understand and predict how climate change-induced alterations could impact the lake carbon dynamics.

Model Communities Hint at Promiscuous Metabolic Linkages between Ubiquitous Free-Living Freshwater Bacteria.

Genome streamlining is frequently observed in free-living aquatic microorganisms and results in physiological dependencies between microorganisms. However, we know little about the specificity of these microbial associations. In order to examine the specificity and extent of these associations, we established mixed cultures from three different freshwater environments and analyzed the cooccurrence of organisms using a metagenomic time series. Free-living microorganisms with streamlined genomes lacking multiple biosynthetic pathways showed no clear recurring pattern in their interaction partners. Free-living freshwater bacteria form promiscuous cooperative associations. This notion contrasts with the well-documented high specificities of interaction partners in host-associated bacteria. Considering all data together, we suggest that highly abundant free-living bacterial lineages are functionally versatile in their interactions despite their distinct streamlining tendencies at the single-cell level. This metabolic versatility facilitates interactions with a variable set of community members.

Environmental DNA (eDNA) detects the pool frog (Pelophylax lessonae) at times when traditional monitoring methods are insensitive.

Detection of endangered species is invaluable for conservation efforts, yet many traditional sampling techniques are ineffective at low population abundances or during certain periods of the year. Here, we compared results from a newly developed eDNA approach and the traditional observational method for the endangered pool frog (Pelophylax lessonae). Analysis using an occupancy-modeling framework indicated that the probability of pools being occupied using eDNA (0.93) was higher than for the traditional method of counting calling males and silent observed individuals (0.72). Detailed analysis revealed complementarity among the methods. That is, the traditional method gave a high rate of observation in June, whereas eDNA gave at least as many or more observations during other parts of the year. Discrepancies among the methods depended on the dominant lifecycle stage, and eDNA concentrations were higher when juveniles were present than at times when spawning occurred. eDNA concentrations were also positively related to P. lessonae observations. Our study demonstrates that an eDNA protocol for monitoring of endangered amphibian species can be particularly valuable during periods when individuals are hard to detect by observational methods, and provides guidance to sampling efforts for research and monitoring programs in other regions and systems.

Effects of predation stress and food ration on perch gut microbiota.

BACKGROUND: Gut microbiota provide functions of importance to influence hosts' food digestion, metabolism, and protection against pathogens. Factors that affect the composition and functions of gut microbial communities are well studied in humans and other animals; however, we have limited knowledge of how natural food web factors such as stress from predators and food resource rations could affect hosts' gut microbiota and how it interacts with host sex. In this study, we designed a two-factorial experiment exposing perch (Perca fluviatilis) to a predator (pike, Esox lucius), and different food ratios, to examine the compositional and functional changes of perch gut microbiota based on 16S rRNA amplicon sequencing. We also investigated if those changes are host sex dependent. RESULTS: We showed that overall gut microbiota composition among individual perch significantly responded to food ration and predator presence. We found that species richness decreased with predator presence, and we identified 23 taxa from a diverse set of phyla that were over-represented when a predator was present. For example, Fusobacteria increased both at the lowest food ration and at predation stress conditions, suggesting that Fusobacteria are favored by stressful situations for the host. In concordance, both food ration and predation stress seemed to influence the metabolic repertoire of the gut microbiota, such as biosynthesis of other secondary metabolites, metabolism of cofactors, and vitamins. In addition, the identified interaction between food ration and sex emphasizes sex-specific responses to diet quantity in gut microbiota. CONCLUSIONS: Collectively, our findings emphasize an alternative state in gut microbiota with responses to changes in natural food webs depending on host sex. The obtained knowledge from this study provided us with an important perspective on gut microbiota in a food web context.

Influence of pulsed and continuous substrate inputs on freshwater bacterial community composition and functioning in bioreactors.

Aquatic environments are typically not homogenous, but characterized by changing substrate concentration gradients and nutrient patches. This heterogeneity in substrate availability creates a multitude of niches allowing bacteria with different substrate utilization strategies to hypothetically coexist even when competing for the same substrate. To study the impact of heterogeneous distribution of organic substrates on bacterioplankton, bioreactors with freshwater bacterial communities were fed artificial freshwater medium with acetate supplied either continuously or in pulses. After a month-long incubation, bacterial biomass and community-level substrate uptake rates were twice as high in the pulsed treatment compared to the continuously fed reactors even if the same total amount of acetate was supplied to both treatments. The composition of the bacterial communities emerging in the two treatments differed significantly with specific taxa overrepresented in the respective treatments. The higher estimated growth yield in cultures that received pulsed substrate inputs, imply that such conditions enable bacteria to use resources more efficiently for biomass production. This finding agrees with established concepts of basal maintenance energy requirements and high energetic costs to assimilate substrates at low concentration. Our results further imply that degradation of organic matter is influenced by temporal and spatial heterogeneity in substrate availability.

Vitamin and Amino Acid Auxotrophy in Anaerobic Consortia Operating under Methanogenic Conditions.

Syntrophy among Archaea and Bacteria facilitates the anaerobic degradation of organic compounds to CH4 and CO2. Particularly during aliphatic and aromatic hydrocarbon mineralization, as in the case of crude oil reservoirs and petroleum-contaminated sediments, metabolic interactions between obligate mutualistic microbial partners are of central importance. Using micromanipulation combined with shotgun metagenomic approaches, we describe the genomes of complex consortia within short-chain alkane-degrading cultures operating under methanogenic conditions. Metabolic reconstruction revealed that only a small fraction of genes in the metagenome-assembled genomes encode the capacity for fermentation of alkanes facilitated by energy conservation linked to H2 metabolism. Instead, the presence of inferred lifestyles based on scavenging anabolic products and intermediate fermentation products derived from detrital biomass was a common feature. Additionally, inferred auxotrophy for vitamins and amino acids suggests that the hydrocarbon-degrading microbial assemblages are structured and maintained by multiple interactions beyond the canonical H2-producing and syntrophic alkane degrader-methanogen partnership. Compared to previous work, our report points to a higher order of complexity in microbial consortia engaged in anaerobic hydrocarbon transformation. IMPORTANCE Microbial interactions between Archaea and Bacteria mediate many important chemical transformations in the biosphere from degrading abundant polymers to synthesis of toxic compounds. Two of the most pressing issues in microbial interactions are how consortia are established and how we can modulate these microbial communities to express desirable functions. Here, we propose that public goods (i.e., metabolites of high energy demand in biosynthesis) facilitate energy conservation for life under energy-limited conditions and determine the assembly and function of the consortia. Our report suggests that an understanding of public good dynamics could result in new ways to improve microbial pollutant degradation in anaerobic systems.

Comparison of four DNA extraction methods for comprehensive assessment of 16S rRNA bacterial diversity in marine biofilms using high-throughput sequencing.

High-throughput DNA sequencing technologies are increasingly used for the metagenomic characterisation of microbial biodiversity. However, basic issues, such as the choice of an appropriate DNA extraction method, are still not resolved for non-model microbial communities. This study evaluates four commonly used DNA extraction methods for marine periphyton biofilms in terms of DNA yield, efficiency, purity, integrity and resulting 16S rRNA bacterial diversity. Among the tested methods, the Plant DNAzol® Reagent (PlantDNAzol) and the FastDNA® SPIN Kit for Soil (FastDNA Soil) methods were best suited to extract high quantities of DNA (77-130 µg g wet wt-1). Lower amounts of DNA were obtained (<37 µg g wet wt-1) with the Power Plant® Pro DNA Isolation Kit (PowerPlant) and the Power Biofilm® DNA Isolation Kit (PowerBiofilm) methods, but integrity and purity of the extracted DNA were higher. Results from 16S rRNA amplicon sequencing demonstrate that the choice of a DNA extraction method significantly influences the bacterial community profiles generated. A higher number of bacterial OTUs were detected when DNA was extracted with the PowerBiofilm and the PlantDNAzol methods. Overall, this study demonstrates the potential bias in metagenomic diversity estimates associated with different DNA extraction methods.