Embracing Environmental Genomics and Machine Learning for Routine Biomonitoring

Using DNA metabarcoding to characterize national scale diatom-environment relationships and to develop indicators in streams and rivers of the United States

Recent advancements in DNA techniques, metabarcoding, and bioinformatics could help expand the use of benthic diatoms in monitoring and assessment programs by providing relatively quick and increasingly cost-effective ways to quantify diatom diversity in environmental samples. However, such applications of DNA-based approaches are relatively new, and in the United States, unknowns regarding their applications at large scales exist because only a few small-scale studies have been done. Here, we present results from the first nationwide survey to use DNA metabarcoding (rbcL) of benthic diatoms, which were collected from 1788 streams and rivers across nine ecoregions spanning the conterminous USA. At the national scale, we found that diatom assemblage structure (1) was strongly associated with total phosphorus and total nitrogen concentrations, conductivity, and pH and (2) had clear patterns that corresponded with differences in these variables among the nine ecoregions. These four variables were strong predictors of diatom assemblage structure in ecoregion-specific analyses, but our results also showed that diatom-environment relationships, the importance of environmental variables, and the ranges of these variables within which assemblage changes occurred differed among ecoregions. To further examine how assemblage data could be used for biomonitoring purposes, we used indicator species analysis to identify ecoregion-specific taxa that decreased or increased along each environmental gradient, and we used their relative abundances of gene reads in samples as metrics. These metrics were strongly correlated with their corresponding variable of interest (e.g., low phosphorus diatoms with total phosphorus concentrations), and generalized additive models showed how their relationships compared among ecoregions. These large-scale national patterns and nine sets of ecoregional results demonstrated that diatom DNA metabarcoding is a robust approach that could be useful to monitoring and assessment programs spanning the variety of conditions that exist throughout the conterminous United States.

Taxonomic accuracy and complementarity between bulk and eDNA metabarcoding provides an alternative to morphology for biological assessment of freshwater macroinvertebrates

Determining biological status of freshwater ecosystems is critical for ensuring ecosystem health and maintaining associated services to such ecosystems. Freshwater macroinvertebrates respond predictably to environmental disturbances and are widely used in biomonitoring programs. However, many freshwater species are difficult to capture and sort from debris or substrate and morphological identification is challenging, especially larval stages, damaged specimens, or hyperdiverse groups such as Diptera. The advent of high throughput sequencing technologies has enhanced DNA barcoding tools to automatise species identification for whole communities, as metabarcoding is increasingly used to monitor biodiversity. However, recent comparisons have revealed little congruence between morphological and molecular-based identifications. Using broad range universal primers for DNA barcode marker cox1, we compare community composition captured between morphological and molecular-based approaches from different sources - tissue-based (bulk benthic and bulk drift samples) and environmental DNA (eDNA, filtered water) metabarcoding - for samples collected along a gradient of anthropogenic disturbances. For comparability, metabarcoding taxonomic assignments were filtered by taxa included in the standardised national biological metric IBMWP. At the family level, bulk benthic metabarcoding showed the highest congruence with morphology, and the most abundant taxa were captured by all techniques. Richness captured by morphology and bulk benthic metabarcoding decreased along the gradient, whereas richness recorded by eDNA remained constant and increased downstream when sequencing bulk drift. Estimates of biological metrics were higher using molecular than morphological identification. At species level, diversity captured by bulk benthic samples were higher than the other techniques. Importantly, bulk benthic and eDNA metabarcoding captured different and complementary portions of the community - benthic versus water column, respectively - and their combined use is recommended. While bulk benthic metabarcoding can likely replace morphology using similar benthic biological indices, water eDNA will require new metrics because this technique sequences a different portion of the community.

Assessing the potential of nematode metabarcoding for benthic monitoring of offshore oil platforms

Environmental DNA metabarcoding is gaining momentum as a time and cost-effective tool for biomonitoring and environmental impact assessment. Yet, its use as a replacement for the conventional marine benthic monitoring based on morphological analysis of macrofauna is still challenging. Here we propose to study the meiofauna, which is much better represented in sediment DNA samples. We focus on nematodes, which are the most numerous and diverse group of meiofauna. Our aim is to assess the potential of nematode metabarcoding to monitor impacts associated with offshore oil platform activities. To achieve this goal, we used nematode-optimized marker (18S V1V2-Nema) and universal eukaryotic marker (18S V9) region to analyse 252 sediment DNA samples collected near three offshore oil platforms in the North Sea. For both markers, we analysed changes in alpha and beta diversity in relation to distance from the platforms and environmental variables. We also defined three impact classes based on selected environmental variables that are associated with oil extraction activities and used random forest classifiers to compare the predictive performance of both datasets. Our results show that alpha- and beta-diversity of nematodes varies with the increasing distance from the platforms. The variables directly related to platform activity, such as Ba and THC, strongly influence the nematode community. The nematode metabarcoding data provide more robust predictive models than eukaryotic data. Furthermore, the nematode community appears more stable in time and space, as illustrated by the overlap of nematode datasets obtained from the same platform three years apart. A significative negative correlation between distance and Shannon diversity also advocates for higher performance of the V1V2-Nema over the V9. Overall, these results suggest that the sensitivity of nematodes is higher compared to the eukaryotic community. Hence, nematode metabarcoding has the potential to become an effective tool for benthic monitoring in marine environment.

The geomicrobiology of limestone, sulfuric acid speleogenetic, and volcanic caves: basic concepts and future perspectives

Caves are ubiquitous subterranean voids, accounting for a still largely unexplored surface of the Earth underground. Due to the absence of sunlight and physical segregation, caves are naturally colonized by microorganisms that have developed distinctive capabilities to thrive under extreme conditions of darkness and oligotrophy. Here, the microbiomes colonizing three frequently studied cave types, i.e., limestone, sulfuric acid speleogenetic (SAS), and lava tubes among volcanic caves, have comparatively been reviewed. Geological configurations, nutrient availability, and energy flows in caves are key ecological drivers shaping cave microbiomes through photic, twilight, transient, and deep cave zones. Chemoheterotrophic microbial communities, whose sustenance depends on nutrients supplied from outside, are prevalent in limestone and volcanic caves, while elevated inorganic chemical energy is available in SAS caves, enabling primary production through chemolithoautotrophy. The 16S rRNA-based metataxonomic profiles of cave microbiomes were retrieved from previous studies employing the Illumina platform for sequencing the prokaryotic V3-V4 hypervariable region to compare the microbial community structures from different cave systems and environmental samples. Limestone caves and lava tubes are colonized by largely overlapping bacterial phyla, with the prevalence of Pseudomonadota and Actinomycetota, whereas the co-dominance of Pseudomonadota and Campylobacterota members characterizes SAS caves. Most of the metataxonomic profiling data have so far been collected from the twilight and transient zones, while deep cave zones remain elusive, deserving further exploration. Integrative approaches for future geomicrobiology studies are suggested to gain comprehensive insights into the different cave types and zones. This review also poses novel research questions for unveiling the metabolic and genomic capabilities of cave microorganisms, paving the way for their potential biotechnological applications.

How do biomarkers dance? Specific moves of defense and damage biomarkers for biological interpretation of dose-response model trends

Omics studies are currently increasingly used in ecotoxicology to highlight the induction of known or novel biomarkers when organisms are exposed to contaminants. Although it is virtually impossible to identify all biomarkers from all organisms, biomarkers can be grouped as defense or damage biomarkers, exhibiting a limited number of response trends. Our working hypothesis is that defense and damage biomarkers follow different dose-response patterns. A meta-analysis of 156 articles and 2595 observations of dose-response curves of defense and damage biomarkers was carried out in order to characterize the response trends of these biological parameters in a large panel of living organisms (18 phyla) exposed to inorganic or organic contaminants (176 in total). Using multinomial logistic regression models, defense biomarkers were found to describe biphasic responses (bell- and U-shaped) to a greater extent (2.5 times) than damage biomarkers. In contrast, damage biomarkers varied mainly monotonically (decreasing or increasing), representing 85% of the observations. Neither the nature of the contaminant nor the type of organisms belonging to 4 kingdoms, influence these specific responses. This result suggests that cellular defense and damage mechanisms are not specific to stressors and are conserved throughout life. Trend analysis of dose-response models as a biological interpretation of biomarkers could thus be a valuable way to exploit large omics datasets.

From biomarkers to community composition: Negative effects of UV/chlorine-treated reclaimed urban wastewater on freshwater biota

The use of urban wastewater reclaimed water has recently increased across the globe to restore stream environmental flows and mitigate the effects of water scarcity. Reclaimed water is disinfected using different treatments, but their effects into the receiving rivers are little studied. Physiological bioassays and biomarkers can detect sub-lethal effects on target species, but do not provide information on changes in community structure. In contrast, official monitoring programs use community structure information but often at coarse taxonomic resolution level that may fail to detect species level impacts. Here, we combined commonly used biomonitoring approaches from organism physiology to community species composition to scan a broad range of effects of disinfection of reclaimed water by UV-light only and both UV/chlorine on the biota. We (1) performed bioassays in one laboratory species (water flea Daphnia magna) and measured biomarkers in two wild species (caddisfly Hydropsyche exocellata and the barbel Luciobarbus graellsii), (2) calculated standard indices of biotic quality (IBQ) for diatoms, benthic macroinvertebrates, and fishes, and (3) analysed community species composition of eukaryotes determined by Cytochrome Oxidase C subunit I (cox1) metabarcoding. Only the UV/chlorine treatment caused significant changes in feeding rates of D. magna and reduced antioxidant defenses, increased anaerobic metabolism and altered the levels of lipid peroxidiation in H. exocellata. However, inputs of reclaimed water were significantly associated with a greater proportion of circulating neutrophils and LG-PAS cells in L. graellsii. Despite IBQ did not discriminate between the two water treatments, metabarcoding data detected community composition changes upon exposure to UV/chlorine reclaimed water. Overall, despite the effects of UV/chlorine-treated water were transient, our study suggests that UV-light treated is less harmful for freshwater biota than UV/chlorine-treated reclaimed water, but those effects depend of the organizational level.

Persistence of location-specific microbial signatures on boats during voyages

Objects collect microorganisms from their surroundings and develop a microbial "fingerprint" that may be useful for determining an object's past location (provenance). It may be possible to use ubiquitous microorganisms for forensics or as environmental sensors. Here, we use microbial communities in the Chesapeake Bay region to demonstrate the use of natural microorganisms as biological sensors to determine the past location of boats. The microbiomes of two boats and of the open water were sampled as these vessels traveled from the Port of Baltimore to the Port of Norfolk, and back to Baltimore. 16S rRNA sequencing was performed to identify microorganisms. Differential abundance and machine learning analyses were utilized to identify microbial signatures and predicted probabilities which were used to determine the vessel's previous location. The work presented here provides a better understanding of how microbes in aquatic systems can be leveraged as utility for object biosensors.

An IoT framework for quality analysis of aquatic water data using time-series convolutional neural network

Water quality monitoring and analysis in fish farms are of paramount importance for the aquaculture sector; however, traditional methods can pose difficulties. To address this challenge, this study proposes an IoT-based deep learning model using a time-series convolution neural network (TMS-CNN) for monitoring and analyzing water quality in fish farms. The proposed TMS-CNN model can handle spatial-temporal data effectively by considering temporal and spatial dependencies between data points, which allows it to capture patterns and trends that would not be possible with traditional models. The model calculates the water quality index (WQI) using correlation analysis and assigns class labels to the data based on the WQI. Then, the TMS-CNN model analyzed the time-series data. It produces high accuracy of 96.2% in analysis of water quality parameters for fish growth and mortality conditions. The proposed model accuracy is higher than the current best model MANN, which has only had an accuracy of 91%.

Supervised machine learning improves general applicability of eDNA metabarcoding for reservoir health monitoring

Effective and standardized monitoring methodologies are vital for successful reservoir restoration and management. Environmental DNA (eDNA) metabarcoding sequencing offers a promising alternative for biomonitoring and can overcome many limitations of traditional morphological bioassessment. Recent attempts have even shown that supervised machine learning (SML) can directly infer biotic indices (BI) from eDNA metabarcoding data, bypassing the cumbersome calculation process of BI regardless of the taxonomic assignment of eDNA sequences. However, questions surrounding the general applicability of this taxonomy-free approach to monitoring reservoir health remain unclear, including model stability, feature selection, algorithm choice, and multi-season biomonitoring. Here, we firstly developed a novel biological integrity index (Me-IBI) that integrates multitrophic interactions and environmental information, based on taxonomy-assigned eDNA metabarcoding data. The Me-IBI can better distinguish the actual health status of the Three Gorges Reservoir (TGR) than physicochemical assessments and have a clear response to human activity. Then, taking this reliable Me-IBI as a supervised label, we compared the impact of selecting different numbers of features and SML algorithms on the stability and predictive performance of the model for predicting ecological conditions in multiple seasons using taxonomy-free eDNA metabarcoding data. We discovered that even with a small number of features, different SML algorithms can establish a stable model and obtain excellent predictive performance. Finally, we proposed a four-step strategy for standardized routine biomonitoring using SML tools. Our study firstly explores the general applicability problem of the taxonomy-free eDNA-SML approach and establishes a solid foundation for the large-scale and standardized biomonitoring application.

Microbial biomarkers of tree water status for next-generation biomonitoring of forest ecosystems

Next-generation biomonitoring proposes to combine machine-learning algorithms with environmental DNA data to automate the monitoring of the Earth's major ecosystems. In the present study, we searched for molecular biomarkers of tree water status to develop next-generation biomonitoring of forest ecosystems. Because phyllosphere microbial communities respond to both tree physiology and climate change, we investigated whether environmental DNA data from tree phyllosphere could be used as molecular biomarkers of tree water status in forest ecosystems. Using an amplicon sequencing approach, we analysed phyllosphere microbial communities of four tree species (Quercus ilex, Quercus robur, Pinus pinaster and Betula pendula) in a forest experiment composed of irrigated and non-irrigated plots. We used these microbial community data to train a machine-learning algorithm (Random Forest) to classify irrigated and non-irrigated trees. The Random Forest algorithm detected tree water status from phyllosphere microbial community composition with more than 90% accuracy for oak species, and more than 75% for pine and birch. Phyllosphere fungal communities were more informative than phyllosphere bacterial communities in all tree species. Seven fungal amplicon sequence variants were identified as candidates for the development of molecular biomarkers of water status in oak trees. Altogether, our results show that microbial community data from tree phyllosphere provides information on tree water status in forest ecosystems and could be included in next-generation biomonitoring programmes that would use in situ, real-time sequencing of environmental DNA to help monitor the health of European temperate forest ecosystems.

Environmental DNA metabarcoding serves as a promising method for aquatic species monitoring and management: A review focused on its workflow, applications, challenges and prospects

Marine and freshwater biodiversity is under threat from both natural and manmade causes. Biological monitoring is currently a top priority for biodiversity protection. Given present limitations, traditional biological monitoring methods may not achieve the proposed monitoring aims. Environmental DNA metabarcoding technology reflects species information by capturing and extracting DNA from environmental samples, using molecular biology techniques to sequence and analyze the DNA, and comparing the obtained information with existing reference libraries to obtain species identification. However, its practical application has highlighted several limitations. This paper summarizes the main steps in the environmental application of eDNA metabarcoding technology in aquatic ecosystems, including the discovery of unknown species, the detection of invasive species, and evaluations of biodiversity. At present, with the rapid development of big data and artificial intelligence, certain advanced technologies and devices can be combined with environmental DNA metabarcoding technology to promote further development of aquatic species monitoring and management.

A combination of machine-learning and eDNA reveals the genetic signature of environmental change at the landscape levels

The current advances of environmental DNA (eDNA) bring profound changes to ecological monitoring and provide unique insights on the biological diversity of ecosystems. The very nature of eDNA data is challenging yet also revolutionizing how biological monitoring information is analysed. In particular, new metrics and approaches should take full advantage of the extent and detail of molecular data produced by genetic methods. In this perspective, machine learning algorithms are particularly promising as they can capture complex relationships between the multiple environmental pressures and the diversity of biological communities. We investigated the potential of a new generation of biomonitoring tools that implement machine-learning techniques to fully exploit eDNA datasets. We trained a machine learning model to discriminate between reference and impacted communities of freshwater macroinvertebrates and assessed its performances using a large eDNA dataset collected at 64 standard federal monitoring sites across Switzerland. We show that a model trained on eDNA is significantly better than a naive model and performs similarly to a model trained on traditional data. Our proof-of-concept shows that such a combination of eDNA and machine learning approaches has the potential to complement or even replace traditional environmental monitoring, and could be scaled along temporal or spatial dimensions.

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.

Current Status of Omics in Biological Quality Elements for Freshwater Biomonitoring

Freshwater ecosystems have been experiencing various forms of threats, mainly since the last century. The severity of this adverse scenario presents unprecedented challenges to human health, water supply, agriculture, forestry, ecological systems, and biodiversity, among other areas. Despite the progress made in various biomonitoring techniques tailored to specific countries and biotic communities, significant constraints exist, particularly in assessing and quantifying biodiversity and its interplay with detrimental factors. Incorporating modern techniques into biomonitoring methodologies presents a challenging topic with multiple perspectives and assertions. This review aims to present a comprehensive overview of the contemporary advancements in freshwater biomonitoring, specifically by utilizing omics methodologies such as genomics, metagenomics, transcriptomics, proteomics, metabolomics, and multi-omics. The present study aims to elucidate the rationale behind the imperative need for modernization in this field. This will be achieved by presenting case studies, examining the diverse range of organisms that have been studied, and evaluating the potential benefits and drawbacks associated with the utilization of these methodologies. The utilization of advanced high-throughput bioinformatics techniques represents a sophisticated approach that necessitates a significant departure from the conventional practices of contemporary freshwater biomonitoring. The significant contributions of omics techniques in the context of biological quality elements (BQEs) and their interpretations in ecological problems are crucial for biomonitoring programs. Such contributions are primarily attributed to the previously overlooked identification of interactions between different levels of biological organization and their responses, isolated and combined, to specific critical conditions.

Comparing quantile regression spline analyses and supervised machine learning for environmental quality assessment at coastal marine aquaculture installations

Organic enrichment associated with marine finfish aquaculture is a local stressor of marine coastal ecosystems. To maintain ecosystem services, the implementation of biomonitoring programs focusing on benthic diversity is required. Traditionally, impact-indices are determined by extracting and identifying benthic macroinvertebrates from samples. However, this is a time-consuming and expensive method with low upscaling potential. A more rapid, inexpensive, and robust method to infer the environmental quality of marine environments is eDNA metabarcoding of bacterial communities. To infer the environmental quality of coastal habitats from metabarcoding data, two taxonomy-free approaches have been successfully applied for different geographical regions and monitoring goals, namely quantile regression splines (QRS) and supervised machine learning (SML). However, their comparative performance remains untested for monitoring the impact of organic enrichment introduced by aquaculture on marine coastal environments. We compared the performance of QRS and SML using bacterial metabarcoding data to infer the environmental quality of 230 aquaculture samples collected from seven farms in Norway and seven farms in Scotland along an organic enrichment gradient. As a measure of environmental quality, we used the Infaunal Quality Index (IQI) calculated from benthic macrofauna data (reference index). The QRS analysis plotted the abundance of amplicon sequence variants (ASVs) as a function to the IQI from which the ASVs with a defined abundance peak were assigned to eco-groups and a molecular IQI was subsequently calculated. In contrast, the SML approach built a random forest model to directly predict the macrofauna-based IQI. Our results show that both QRS and SML perform well in inferring the environmental quality with 89% and 90% accuracy, respectively. For both geographic regions, there was high correspondence between the reference IQI and both the inferred molecular IQIs (p < 0.001), with the SML model showing a higher coefficient of determination compared to QRS. Among the 20 most important ASVs identified by the SML approach, 15 were congruent with the good quality spline ASV indicators identified via QRS for both Norwegian and Scottish salmon farms. More research on the response of the ASVs to organic enrichment and the co-influence of other environmental parameters is necessary to eventually select the most powerful stressor-specific indicators. Even though both approaches are promising to infer environmental quality based on metabarcoding data, SML showed to be more powerful in handling the natural variability. For the improvement of the SML model, addition of new samples is still required, as background noise introduced by high spatio-temporal variability can be reduced. Overall, we recommend the development of a powerful SML approach that will be onwards applied for monitoring the impact of aquaculture on marine ecosystems based on eDNA metabarcoding data.

Compact and automated eDNA sampler for in situ monitoring of marine environments

Using environmental DNA (eDNA) to monitor biodiversity in aquatic environments is becoming an efficient and cost-effective alternative to other methods such as visual and acoustic identification. Until recently, eDNA sampling was accomplished primarily through manual sampling methods; however, with technological advances, automated samplers are being developed to make sampling easier and more accessible. This paper describes a new eDNA sampler capable of self-cleaning and multi-sample capture and preservation, all within a single unit capable of being deployed by a single person. The first in-field test of this sampler took place in the Bedford Basin, Nova Scotia, Canada alongside parallel samples taken using the typical Niskin bottle collection and post-collection filtration method. Both methods were able to capture the same aquatic microbial community and counts of representative DNA sequences were well correlated between methods with R[Formula: see text] values ranging from 0.71-0.93. The two collection methods returned the same top 10 families in near identical relative abundance, demonstrating that the sampler was able to capture the same community composition of common microbes as the Niskin. The presented eDNA sampler provides a robust alternative to manual sampling methods, is amenable to autonomous vehicle payload constraints, and will facilitate persistent monitoring of remote and inaccessible sites.

Using Machine Learning to make nanomaterials sustainable

Sustainable development is a key challenge for contemporary human societies; failure to achieve sustainability could threaten human survival. In this review article, we illustrate how Machine Learning (ML) could support more sustainable development, covering the basics of data gathering through each step of the Environmental Risk Assessment (ERA). The literature provides several examples showing how ML can be employed in most steps of a typical ERA.A key observation is that there are currently no clear guidance for using such autonomous technologies in ERAs or which standards/checks are required. Steering thus seems to be the most important task for supporting the use of ML in the ERA of nano- and smart-materials. Resources should be devoted to developing a strategy for implementing ML in ERA with a strong emphasis on data foundations, methodologies, and the related sensitivities/uncertainties. We should recognise historical errors and biases (e.g., in data) to avoid embedding them during ML programming.

Shared and contrasting associations in the dynamic nano- and picoplankton communities of two close but contrasting sites from the Bay of Biscay

Pico- and nanoplankton are key players in the marine ecosystems due to their implication in the biogeochemical cycles, nutrient recycling and the pelagic food webs. However, the specific dynamics and niches of most bacterial, archaeal and eukaryotic plankton remain unknown, as well as the interactions between them. Better characterization of these is critical for understanding and predicting ecosystem functioning under anthropogenic pressures. We used environmental DNA metabarcoding across a 6-year time series to explore the structure and seasonality of pico- and nanoplankton communities in two sites of the Bay of Biscay, one coastal and one offshore, and construct association networks to reveal potential keystone and connector taxa. Temporal trends in alpha diversity were similar between the two sites, and concurrent communities more similar than within the same site at different times. However, we found differences between the network topologies of the two sites, with both shared and site-specific keystones and connectors. For example, Micromonas, with lower abundance in the offshore site is a keystone here, indicating a stronger effect of associations such as resource competition. This study provides an example of how time series and association network analysis can reveal how similar communities may function differently despite being geographically close.

Environmental DNA-based profiling of benthic bacterial and eukaryote communities along a crude oil spill gradient in a coral reef in the Persian Gulf

Coral reef ecosystems in the Persian Gulf are frequently exposed to crude oil spills. We investigated benthic bacterial and eukaryote community structures at such coral reef sites subjected to different degrees of polycyclic aromatic hydrocarbon (PAH) pollution using environmental DNA (eDNA) metabarcoding. Both bacterial and eukaryote communities responded with pronounced shifts to crude oil pollution and distinguished control sites, moderately and heavily impacted sites with significant confidentiality. The observed community patterns were predominantly driven by Alphaproteobacteria and metazoans. Among these, we identified individual genera that were previously linked to oil spill stress, but also taxa, for which a link to hydrocarbon still remains to be established. Considering the lack of an early-warning system for the environmental status of coral reef ecosystems exposed to frequent crude-oil spills, our results encourage further research towards the development of an eDNA-based biomonitoring tool that exploits benthic bacterial and eukaryote communities as bioindicators.

Supervised Machine Learning Enables Geospatial Microbial Provenance

The recent increase in publicly available metagenomic datasets with geospatial metadata has made it possible to determine location-specific, microbial fingerprints from around the world. Such fingerprints can be useful for comparing microbial niches for environmental research, as well as for applications within forensic science and public health. To determine the regional specificity for environmental metagenomes, we examined 4305 shotgun-sequenced samples from the MetaSUB Consortium dataset-the most extensive public collection of urban microbiomes, spanning 60 different cities, 30 countries, and 6 continents. We were able to identify city-specific microbial fingerprints using supervised machine learning (SML) on the taxonomic classifications, and we also compared the performance of ten SML classifiers. We then further evaluated the five algorithms with the highest accuracy, with the city and continental accuracy ranging from 85-89% to 90-94%, respectively. Thereafter, we used these results to develop Cassandra, a random-forest-based classifier that identifies bioindicator species to aid in fingerprinting and can infer higher-order microbial interactions at each site. We further tested the Cassandra algorithm on the Tara Oceans dataset, the largest collection of marine-based microbial genomes, where it classified the oceanic sample locations with 83% accuracy. These results and code show the utility of SML methods and Cassandra to identify bioindicator species across both oceanic and urban environments, which can help guide ongoing efforts in biotracing, environmental monitoring, and microbial forensics (MF).

Metagenomics versus total RNA sequencing: most accurate data-processing tools, microbial identification accuracy and perspectives for ecological assessments

Metagenomics and total RNA sequencing (total RNA-Seq) have the potential to improve the taxonomic identification of diverse microbial communities, which could allow for the incorporation of microbes into routine ecological assessments. However, these target-PCR-free techniques require more testing and optimization. In this study, we processed metagenomics and total RNA-Seq data from a commercially available microbial mock community using 672 data-processing workflows, identified the most accurate data-processing tools, and compared their microbial identification accuracy at equal and increasing sequencing depths. The accuracy of data-processing tools substantially varied among replicates. Total RNA-Seq was more accurate than metagenomics at equal sequencing depths and even at sequencing depths almost one order of magnitude lower than those of metagenomics. We show that while data-processing tools require further exploration, total RNA-Seq might be a favorable alternative to metagenomics for target-PCR-free taxonomic identifications of microbial communities and might enable a substantial reduction in sequencing costs while maintaining accuracy. This could be particularly an advantage for routine ecological assessments, which require cost-effective yet accurate methods, and might allow for the incorporation of microbes into ecological assessments.

Developing Indicators of Nutrient Pollution in Streams Using 16S rRNA Gene Metabarcoding of Periphyton-Associated Bacteria

Indicators based on nutrient-biota relationships in streams can inform water quality restoration and protection programs. Bacterial assemblages could be particularly useful indicators of nutrient effects because they are species-rich, important contributors to ecosystem processes in streams, and responsive to rapidly changing conditions. Here, we sampled 25 streams weekly (12-14 times each) and used 16S rRNA gene metabarcoding of periphyton-associated bacteria to quantify the effects of total phosphorus (TP) and total nitrogen (TN). Threshold indicator taxa analysis identified assemblage-level changes and amplicon sequence variants (ASVs) that increased or decreased with increasing TP and TN concentrations (i.e., low P, high P, low N, and high N ASVs). Boosted regression trees confirmed that relative abundances of gene sequence reads for these four indicator groups were associated with nutrient concentrations. Gradient forest analysis complemented these results by using multiple predictors and random forest models for each ASV to identify portions of TP and TN gradients at which the greatest changes in assemblage structure occurred. Synthesized statistical results showed bacterial assemblage structure began changing at 24 μg TP/L with the greatest changes occurring from 110 to 195 μg/L. Changes in the bacterial assemblages associated with TN gradually occurred from 275 to 855 μg/L. Taxonomic and phylogenetic analyses showed that low nutrient ASVs were commonly Firmicutes, Verrucomicrobiota, Flavobacteriales, and Caulobacterales, Pseudomonadales, and Rhodobacterales of Proteobacteria, whereas other groups, such as Chitinophagales of Bacteroidota, and Burkholderiales, Rhizobiales, Sphingomonadales, and Steroidobacterales of Proteobacteria comprised the high nutrient ASVs. Overall, the responses of bacterial ASV indicators in this study highlight the utility of metabarcoding periphyton-associated bacteria for quantifying biotic responses to nutrient inputs in streams.

Surface microbiota of Mediterranean loggerhead sea turtles unraveled by 16S and 18S amplicon sequencing

The loggerhead sea turtle is considered a keystone species with a major ecological role in Mediterranean marine environment. As is the case with other wild reptiles, their outer microbiome is rarely studied. Although there are several studies on sea turtle’s macro-epibionts and endo-microbiota, there has been little research on epibiotic microbiota associated with turtle skin and carapace. Therefore we aimed to provide the identification of combined epibiotic eukaryotic, bacterial and archaeal microbiota on Mediterranean loggerhead sea turtles. In this study, we sampled skins and carapaces of 26 loggerheads from the Mediterranean Sea during 2018 and 2019. To investigate the overall microbial diversity and composition, amplicon sequencing of 16S and 18S rRNA genes was performed. We found that the Mediterranean loggerhead sea turtle epibiotic microbiota is a reservoir of a vast variety of microbial species. Microbial communities mostly varied by different locations and seas, while within bacterial communities’ significant difference was observed between sampled body sites (carapace vs. skin). In terms of relative abundance, Proteobacteria and Bacteroidota were the most represented phyla within prokaryotes, while Alveolata and Stramenopiles thrived among eukaryotes. This study, besides providing a first survey of microbial eukaryotes on loggerheads via metabarcoding, identifies fine differences within both bacterial and eukaryotic microbial communities that seem to reflect the host anatomy and habitat. Multi-domain epi-microbiome surveys provide additional layers of information that are complementary with previous morphological studies and enable better understanding of the biology and ecology of these vulnerable marine reptiles.

Characterizing temporal variability in streams supports nutrient indicator development using diatom and bacterial DNA metabarcoding

Interest in developing periphytic diatom and bacterial indicators of nutrient effects continues to grow in support of the assessment and management of stream ecosystems and their watersheds. However, temporal variability could confound relationships between indicators and nutrients, subsequently affecting assessment outcomes. To document how temporal variability affects measures of diatom and bacterial assemblages obtained from DNA metabarcoding, we conducted weekly periphyton and nutrient sampling from July to October 2016 in 25 streams in a 1293 km² mixed land use watershed. Measures of both diatom and bacterial assemblages were strongly associated with the percent agriculture in upstream watersheds and total phosphorus (TP) and total nitrogen (TN) concentrations. Temporal variability in TP and TN concentrations increased with greater amounts of agriculture in watersheds, but overall diatom and bacterial assemblage variability within sites-measured as mean distance among samples to corresponding site centroids in ordination space-remained consistent. This consistency was due in part to offsets between decreasing variability in relative abundances of taxa typical of low nutrient conditions and increasing variability in those typical of high nutrient conditions as mean concentrations of TP and TN increased within sites. Weekly low and high nutrient diatom and bacterial metrics were more strongly correlated with site mean nutrient concentrations over the sampling period than with same day measurements and more strongly correlated with TP than with TN. Correlations with TP concentrations were consistently strong throughout the study except briefly following two major precipitation events. Following these events, biotic relationships with TP reestablished within one to three weeks. Collectively, these results can strengthen interpretations of survey results and inform monitoring strategies and decision making. These findings have direct applications for improving the use of diatoms and bacteria, and the use of DNA metabarcoding, in monitoring programs and stream site assessments.

Viral Eco-Genomic Tools: Development and Implementation for Aquatic Biomonitoring

Enteric viruses (EVs) occurrence within aquatic environments varies and leads to significant risk on public health of humans, animals, and diversity of aquatic taxa. Early and efficacious recognition of cultivable and fastidious EVs in aquatic systems are important to ensure the sanitary level of aquatic water and implement required treatment strategies. Herein, we provided a comprehensive overview of the conventional and up-to-date eco-genomic tools for aquatic biomonitoring of EVs, aiming to develop better water pollution monitoring tools. In combination with bioinformatics techniques, genetic tools including cloning sequencing analysis, DNA microarray, next-generation sequencing (NGS), and metagenomic sequencing technologies are implemented to make informed decisions about the global burden of waterborne EVs-associated diseases. The data presented in this review are helpful to recommend that: (1) Each viral pollution detection method has its own merits and demerits; therefore, it would be advantageous for viral pollution evaluation to be integrated as a complementary platform. (2) The total viral genome pool extracted from aquatic environmental samples is a real reflection of pollution status of the aquatic eco-systems; therefore, it is recommended to conduct regular sampling through the year to establish an updated monitoring system for EVs, and quantify viral peak concentrations, viral typing, and genotyping. (3) Despite that conventional detection methods are cheaper, it is highly recommended to implement molecular-based technologies to complement aquatic ecosystems biomonitoring due to numerous advantages including high-throughput capability. (4) Continuous implementation of the eco-genetic detection tools for monitoring the EVs in aquatic ecosystems is recommended.

Applying convolutional neural networks to speed up environmental DNA annotation in a highly diverse ecosystem

High-throughput DNA sequencing is becoming an increasingly important tool to monitor and better understand biodiversity responses to environmental changes in a standardized and reproducible way. Environmental DNA (eDNA) from organisms can be captured in ecosystem samples and sequenced using metabarcoding, but processing large volumes of eDNA data and annotating sequences to recognized taxa remains computationally expensive. Speed and accuracy are two major bottlenecks in this critical step. Here, we evaluated the ability of convolutional neural networks (CNNs) to process short eDNA sequences and associate them with taxonomic labels. Using a unique eDNA data set collected in highly diverse Tropical South America, we compared the speed and accuracy of CNNs with that of a well-known bioinformatic pipeline (OBITools) in processing a small region (60 bp) of the 12S ribosomal DNA targeting freshwater fishes. We found that the taxonomic labels from the CNNs were comparable to those from OBITools, with high correlation levels for the composition of the regional fish fauna. The CNNs enabled the processing of raw fastq files at a rate of approximately 1 million sequences per minute, which was about 150 times faster than with OBITools. Given the good performance of CNNs in the highly diverse ecosystem considered here, the development of more elaborate CNNs promises fast deployment for future biodiversity inventories using eDNA.

Scientific and budgetary trade-offs between morphological and molecular methods for deep-sea biodiversity assessment

Deep-sea biodiversity, a source of critical ecological functions and ecosystem services, is increasingly subject to the threat of disturbance from existing practices (e.g., fishing, waste disposal, oil and gas extraction) as well as emerging industries such as deep-seabed mining. Current scientific tools may not be adequate for monitoring and assessing subsequent changes to biodiversity. In this paper, we evaluate the scientific and budgetary trade-offs associated with morphology-based taxonomy and metabarcoding approaches to biodiversity surveys in the context of nascent deep-seabed mining for polymetallic nodules in the Clarion-Clipperton Zone, the area of most intense interest. For the dominant taxa of benthic meiofauna, we discuss the types of information produced by these methods and use cost-effectiveness analysis to compare their abilities to yield biological and ecological data for use in environmental assessment and management. On the basis of our evaluation, morphology-based taxonomy is less cost-effective than metabarcoding but offers scientific advantages, such as the generation of density, biomass, and size structure data. Approaches that combine the two methods during the environmental assessment phase of commercial activities may facilitate future biodiversity monitoring and assessment for deep-seabed mining and for other activities in remote deep-sea habitats, for which taxonomic data and expertise are limited. Integr Environ Assess Manag 2022;18:655-663. © 2021 SETAC.

The Effects of Water Level Fluctuation on Zooplankton Communities in Shahu Lake Based on DNA Metabarcoding and Morphological Methods

Background: The water level of Poyang Lake (China) fluctuates seasonally. Shahu Lake, a smaller body of water connected to Poyang Lake during the wet season, is separated in the dry season. Due to a special fishing method termed ‘lake enclosed in autumn’, the water level is lowered and reaches its lowest point in January, which is <0.5 m deep in the middle of the lake. Our research investigated the effect of water level changes on the zooplankton community composition in Shahu Lake. Methods: We used both DNA metabarcoding method (MBC) (18S rRNA gene V4 region) and morphological method (MOI) to track the zooplankton community structure over four seasons in Shahu Lake (China). Results: Totals of 90 and 98 species of zooplankton were detected by MOI and MBC, respectively, with rotifers being the main zooplankton component. The α-diversity index of both methods increased from spring to summer and decreased from summer to autumn, reaching the lowest value in winter. NMDS and a cluster analysis showed that all zooplankton communities detected by MOI and MBC were significantly separated by season. The zooplankton community in winter was separated from that of the other three seasons, but the summer and autumn communities were more similar. Conclusions: Changes in the water level had significant effects on the zooplankton community composition. We found that MBC was more able to detect the differences in the zooplankton composition than MOI. MBC also had more advantages in copepod recognition. In our study, 37 species of copepods were detected by MBC, but only 11 species were detected by MOI. We concluded that MBC should be used to research the seasonal variations of zooplankton.

Multifeature superposition analysis of the effects of microplastics on microbial communities in realistic environments

Microplastic (MP) contamination has become an increasingly serious environmental problem. However, the risks of MP contamination in complex global climatic and geographic scenarios remain unclear. We established a multifeature superposition analysis boosting (MFAB) machine learning (ML) approach to address the above knowledge gap. MFAB-ML identified and predicted the importance, interaction networks and superposition effects of multiple features, including 34 characteristic variables (e.g., MP contamination and climatic and geographic variables), from 1354 samples distributed globally. MFAB-ML analysis achieved realistic and significant results, in some cases even opposite to those obtained using a single or a few features, revealing the importance of considering complicated scenarios. We found that the microbial diversity in East Asian seas will continually decrease due to the superposition effects of MPs with ocean warming; for example, the Chao1 index will decrease by 10.32% by 2065. The present work provides a powerful approach to identify and predict the multifeature superposition effects of pollutants on realistic environments in complicated climatic and geographic scenarios, overcoming the bias from general studies.

Advancement of Metatranscriptomics towards Productive Agriculture and Sustainable Environment: A Review

While chemical fertilisers and pesticides indeed enhance agricultural productivity, their excessive usage has been detrimental to environmental health. In addressing this matter, the use of environmental microbiomes has been greatly favoured as a ‘greener’ alternative to these inorganic chemicals’ application. Challenged by a significant proportion of unidentified microbiomes with unknown ecological functions, advanced high throughput metatranscriptomics is prudent to overcome the technological limitations in unfolding the previously undiscovered functional profiles of the beneficial microbiomes. Under this context, this review begins by summarising (1) the evolution of next-generation sequencing and metatranscriptomics in leveraging the microbiome transcriptome profiles through whole gene expression profiling. Next, the current environmental metatranscriptomics studies are reviewed, with the discussion centred on (2) the emerging application of the beneficial microbiomes in developing fertile soils and (3) the development of disease-suppressive soils as greener alternatives against biotic stress. As sustainable agriculture focuses not only on crop productivity but also long-term environmental sustainability, the second half of the review highlights the metatranscriptomics’ contribution in (4) revolutionising the pollution monitoring systems via specific bioindicators. Overall, growing knowledge on the complex microbiome functional profiles is imperative to unlock the unlimited potential of agricultural microbiome-based practices, which we believe hold the key to productive agriculture and sustainable environment.

Tracking the impacts of nutrient inputs on estuary ecosystem function

Estuaries are one of the most impacted coastal environments globally, subjected to multiple stressors from urban, industry and coastal development. With increasing anthropogenic activity surrounding estuarine systems, sewage inputs have become a common concern. Stable isotope analysis provides a well-established tool to investigate the incorporation of nitrogen into marine organisms and identify major nutrient sources. Benthic macroinvertebrate communities are often used as bioindicators in ecological studies as they typically display predictable responses to anthropogenic pressures, however have a suite of limitations and costs associated with their use. 16S rDNA amplicon sequencing techniques allow for investigation of the microbial communities inhabiting complex environmental samples, with potential as a tool in the ecological assessment of pollution. These communities have not yet been adequately considered for ecological studies and biomonitoring, with a need to better understand interactions with environmental stressors and implications for ecosystem function. This study used a combination of stable isotope analysis to trace the uptake of anthropogenic nitrogen in biota, traditional assessment of benthic macroinvertebrate communities, and 16S rDNA genotyping of benthic microbial communities. Stable isotope analysis of seagrass and epiphytes identified multiple treated and untreated sewage inputs, ranges of 5.2-7.2‰ and 1.9-4.0‰ for δ¹⁵N respectively, as the dominant nitrogen source at specific locations. The benthic macroinvertebrate community reflected these inputs with shifts in dominant taxa and high abundances of polychaetes at some sites. Microbial communities provided a sensitive indication of impact with a breadth of information not available using traditional techniques. Composition and predicted function reflected sewage inputs, particularly within sediments, with the relative abundance of specific taxa and putative pathogens linked to these inputs. This research supports the growing body of evidence that benthic microbial communities respond rapidly to anthropogenic stressors and have potential as a monitoring tool in urban estuarine systems.

Alkaline mine drainage drives stream sediment microbial community structure and function

With advances in eDNA metabarcoding, environmental microbiomes are increasingly used as cost-effective tools for monitoring ecosystem health. Stream ecosystems in Central Appalachia, heavily impacted by alkaline drainage from mountaintop coal mining, present ideal opportunities for biomonitoring using stream microbiomes, but the structural and functional responses of microbial communities in different environmental compartments are not well understood. We investigated sediment microbiomes in mining impacted streams to determine how community composition and function respond to mining and to look for potential microbial bioindicators. Using 16s rRNA gene amplicon sequencing, we found that mining leads to shifts in microbial community structure, with the phylum Planctomycetes enriched by 1-6% at mined sites. We observed ~51% increase in species richness in bulk sediments. In contrast, of the 31 predicted metabolic pathways that changed significantly with mining, 23 responded negatively. Mining explained 15-18% of the variance in community structure and S, Se, %C and %N were the main drivers of community and functional pathway composition. We identified 12 microbial indicators prevalent in the ecosystem and sensitive to mining. Overall, alkaline mountaintop mining drainage causes a restructuration of the sediment microbiome, and our study identified promising microbial indicators for the long-term monitoring of these impacted streams.

Microbial diversity alteration reveals biomarkers of contamination in soil-river-lake continuum

Microbial communities inhabiting soil-water-sediment continuum in coastal areas provide important ecosystem services. Their adaptation in response to environmental stressors, particularly mitigating the impact of pollutants discharged from human activities, has been considered for the development of microbial biomonitoring tools, but their use is still in the infancy. Here, chemical and molecular (16S rRNA gene metabarcoding) approaches were combined in order to determine the impact of pollutants on microbial assemblages inhabiting the aquatic network of a soil-water-sediment continuum around the Ichkeul Lake (Tunisia), an area highly impacted by human activities. Samples were collected within the soil-river-lake continuum at three stations in dry (summer) and wet (winter) seasons. The contaminant pressure index (PI), which integrates Polycyclic aromatic hydrocarbons (PAHs), alkanes, Organochlorine pesticides (OCPs) and metal contents, and the microbial pressure index microgAMBI, based on bacterial community structure, showed significant correlation with contamination level and differences between seasons. The comparison of prokaryotic communities further revealed specific assemblages for soil, river and lake sediments. Correlation analyses identified potential "specialist" genera for the different compartments, whose abundances were correlated with the pollutant type found. Additionally, PICRUSt analysis revealed the metabolic potential for pollutant transformation or degradation of the identified "specialist" species, providing information to estimate the recovery capacity of the ecosystem. Such findings offer the possibility to define a relevant set of microbial indicators for assessing the effects of human activities on aquatic ecosystems. Microbial indicators, including the detection of "specialist" and sensitive taxa, and their functional capacity, might be useful, in combination with integrative microbial indices, to constitute accurate biomonitoring tools for the management and restoration of complex coastal aquatic systems.

Consideration of Multitrophic Biodiversity and Ecosystem Functions Improves Indices on River Ecological Status

Biological quality elements have been developed worldwide to assess whether a water body is in a good status or not. However, current studies mainly focus on a single taxonomic group or a small set of species, often limited by methods of morphological identification, and lack further aspects of biodiversity (e.g., across taxa and multiple attributes) and ecosystem functions. Here, we advance a framework for assessing the river's ecological status based on complete biodiversity data measured by environmental DNA (eDNA) metabarcoding and measurements of ecosystem functions in addition to physicochemical elements across a large riverine system in China. We identified 40 indicators of biodiversity and ecosystem functions, covering five taxonomic groups from bacteria to invertebrates, and associated with multiple attributes of biodiversity and ecosystem functions. Our data show that human impact on ecosystems could be accurately predicted by these eDNA-based indicators and ecosystem functions, using cross-validation with a known stressor gradient. Moreover, indices based on these indicators of biodiversity and ecosystem functions not only distinguish the physicochemical characteristics of the sites but also improve the assessment accuracy of 20-30% for the river's ecological status. Overall, by incorporating eDNA-based biodiversity with physicochemical and ecosystem functional elements, the multidimensional perspectives of ecosystem states provide additional information to protect and maintain a good ecological status of rivers.

Revision and annotation of DNA barcode records for marine invertebrates: report of the 8th iBOL conference hackathon

The accuracy of specimen identification through DNA barcoding and metabarcoding relies on reference libraries containing records with reliable taxonomy and sequence quality. The considerable growth in barcode data requires stringent data curation, especially in taxonomically difficult groups such as marine invertebrates. A major effort in curating marine barcode data in the Barcode of Life Data Systems (BOLD) was undertaken during the 8th International Barcode of Life Conference (Trondheim, Norway, 2019). Major taxonomic groups (crustaceans, echinoderms, molluscs, and polychaetes) were reviewed to identify those which had disagreement between Linnaean names and Barcode Index Numbers (BINs). The records with disagreement were annotated with four tags: a) MIS-ID (misidentified, mislabeled, or contaminated records), b) AMBIG (ambiguous records unresolved with the existing data), c) COMPLEX (species names occurring in multiple BINs), and d) SHARE (barcodes shared between species). A total of 83,712 specimen records corresponding to 7,576 species were reviewed and 39% of the species were tagged (7% MIS-ID, 17% AMBIG, 14% COMPLEX, and 1% SHARE). High percentages (&gt;50%) of AMBIG tags were recorded in gastropods, whereas COMPLEX tags dominated in crustaceans and polychaetes. The high proportion of tagged species reflects either flaws in the barcoding workflow (e.g., misidentification, cross-contamination) or taxonomic difficulties (e.g., synonyms, undescribed species). Although data curation is essential for barcode applications, such manual attempts to examine large datasets are unsustainable and automated solutions are extremely desirable.

Consistent declines in aquatic biodiversity across diverse domains of life in rivers impacted by surface coal mining

The rivers of Appalachia (United States) are among the most biologically diverse freshwater ecosystems in the temperate zone and are home to numerous endemic aquatic organisms. Throughout the Central Appalachian ecoregion, extensive surface coal mines generate alkaline mine drainage that raises the pH, salinity, and trace element concentrations in downstream waters. Previous regional assessments have found significant declines in stream macroinvertebrate and fish communities after draining these mined areas. Here, we expand these assessments with a more comprehensive evaluation across a broad range of organisms (bacteria, algae, macroinvertebrates, all eukaryotes, and fish) using high-throughput amplicon sequencing of environmental DNA (eDNA). We collected water samples from 93 streams in Central Appalachia (West Virginia, United States) spanning a gradient of mountaintop coal mining intensity and legacy to assess how this land use alters downstream water chemistry and affects aquatic biodiversity. For each group of organisms, we identified the sensitive and tolerant taxa along the gradient and calculated stream specific conductivity thresholds in which large synchronous declines in diversity were observed. Streams below mining operations had steep declines in diversity (-18 to -41%) and substantial shifts in community composition that were consistent across multiple taxonomic groups. Overall, large synchronous declines in bacterial, algal, and macroinvertebrate communities occurred even at low levels of mining impact at stream specific conductivity thresholds of 150-200 µS/cm that are substantially below the current U.S. Environmental Protection Agency aquatic life benchmark of 300 µS/cm for Central Appalachian streams. We show that extensive coal surface mining activities led to the extirpation of 40% of biodiversity from impacted rivers throughout the region and that current water quality criteria are likely not protective for many groups of aquatic organisms.

Ecosystems monitoring powered by environmental genomics: A review of current strategies with an implementation roadmap

A decade after environmental scientists integrated high-throughput sequencing technologies in their toolbox, the genomics-based monitoring of anthropogenic impacts on the biodiversity and functioning of ecosystems is yet to be implemented by regulatory frameworks. Despite the broadly acknowledged potential of environmental genomics to this end, technical limitations and conceptual issues still stand in the way of its broad application by end-users. In addition, the multiplicity of potential implementation strategies may contribute to a perception that the routine application of this methodology is premature or "in development", hence restraining regulators from binding these tools into legal frameworks. Here, we review recent implementations of environmental genomics-based methods, applied to the biomonitoring of ecosystems. By taking a general overview, without narrowing our perspective to particular habitats or groups of organisms, this paper aims to compare, review and discuss the strengths and limitations of four general implementation strategies of environmental genomics for monitoring: (a) Taxonomy-based analyses focused on identification of known bioindicators or described taxa; (b) De novo bioindicator analyses; (c) Structural community metrics including inferred ecological networks; and (d) Functional community metrics (metagenomics or metatranscriptomics). We emphasise the utility of the three latter strategies to integrate meiofauna and microorganisms that are not traditionally utilised in biomonitoring because of difficult taxonomic identification. Finally, we propose a roadmap for the implementation of environmental genomics into routine monitoring programmes that leverage recent analytical advancements, while pointing out current limitations and future research needs.

A pilot experiment to assess the efficiency of pharmaceutical plant wastewater treatment and the decreasing effluent toxicity to periphytic biofilms

Pharmaceutical industry effluents are complex and highly variable in time. Assessing the efficiency of a pharmaceutical industry wastewater treatment plant (WWTP) and the resulting decrease in effluent toxicity and ecological risk is thus not straightforward. We set up an original in situ pilot directly connected to a pharmaceutical WWTP to monitor the chronic toxicity of successive effluents using natural periphytic biofilms. Their structural and functional responses to effluent exposure were assessed by combining (i) a molecular approach to characterize the bacterial and diatom diversity and (ii) functional measurements of photosynthetic and enzyme activities. Effluent contamination by pharmaceuticals strongly decreased after the quaternary treatment (activated carbon). Most of the structural biological characteristics improved with cumulative WWTP treatment (bacterial diversity, microbial genetic structure, and biological diatom index), showing community recovery along the treatment process. However, functional parameters did not show clear links with treatment steps, suggesting that microbial activities were not solely driven by pharmaceuticals produced during the experimental period. Operationally, this type of pilot system offers a useful tool for biomonitoring approaches and offers new approaches for industrial managers to assess the ecological risk of production effluents in receiving water.

KAUST Metagenomic Analysis Platform (KMAP), enabling access to massive analytics of re-annotated metagenomic data

Exponential rise of metagenomics sequencing is delivering massive functional environmental genomics data. However, this also generates a procedural bottleneck for on-going re-analysis as reference databases grow and methods improve, and analyses need be updated for consistency, which require acceess to increasingly demanding bioinformatic and computational resources. Here, we present the KAUST Metagenomic Analysis Platform (KMAP), a new integrated open web-based tool for the comprehensive exploration of shotgun metagenomic data. We illustrate the capacities KMAP provides through the re-assembly of ~ 27,000 public metagenomic samples captured in ~ 450 studies sampled across ~ 77 diverse habitats. A small subset of these metagenomic assemblies is used in this pilot study grouped into 36 new habitat-specific gene catalogs, all based on full-length (complete) genes. Extensive taxonomic and gene annotations are stored in Gene Information Tables (GITs), a simple tractable data integration format useful for analysis through command line or for database management. KMAP pilot study provides the exploration and comparison of microbial GITs across different habitats with over 275 million genes. KMAP access to data and analyses is available at https://www.cbrc.kaust.edu.sa/aamg/kmap.start .

Identifying the minimum amplicon sequence depth to adequately predict classes in eDNA-based marine biomonitoring using supervised machine learning

Environmental DNA metabarcoding is a powerful approach for use in biomonitoring and impact assessments. Amplicon-based eDNA sequence data are characteristically highly divergent in sequencing depth (total reads per sample) as influenced inter alia by the number of samples simultaneously analyzed per sequencing run. The random forest (RF) machine learning algorithm has been successfully employed to accurately classify unknown samples into monitoring categories. To employ RF to eDNA data, and avoid sequencing-depth artifacts, sequence data across samples are normalized using rarefaction, a process that inherently loses information. The aim of this study was to inform future sampling designs in terms of the relationship between sampling depth and RF accuracy. We analyzed three published and one new bacterial amplicon datasets, using a RF, based initially on the maximal rarefied data available (minimum mean of > 30,000 reads across all datasets) to give our baseline performance. We then evaluated the RF classification success based on increasingly rarefied datasets. We found that extreme to moderate rarefaction (50-5000 sequences per sample) was sufficient to achieve prediction performance commensurate to the full data, depending on the classification task. We did not find that the number of classification classes, data balance across classes, or the total number of sequences or samples, were associated with predictive accuracy. We identified the ability of the training data to adequately characterize the classes being mapped as the most important criterion and discuss how this finding can inform future sampling design for eDNA based biomonitoring to reduce costs and computation time.

Quantification of the covariation of lake microbiomes and environmental variables using a machine learning-based framework

It is known that microorganisms are essential for the functioning of ecosystems, but the extent to which microorganisms respond to different environmental variables in their natural habitats is not clear. In the current study, we present a methodological framework to quantify the covariation of the microbial community of a habitat and environmental variables of this habitat. It is built on theoretical considerations of systems ecology, makes use of state-of-the-art machine learning techniques and can be used to identify bioindicators. We apply the framework to a data set containing operational taxonomic units (OTUs) as well as more than twenty physicochemical and geographic variables measured in a large-scale survey of European lakes. While a large part of variation (up to 61%) in many environmental variables can be explained by microbial community composition, some variables do not show significant covariation with the microbial lake community. Moreover, we have identified OTUs that act as "multitask" bioindicators, i.e., that are indicative for multiple environmental variables, and thus could be candidates for lake water monitoring schemes. Our results represent, for the first time, a quantification of the covariation of the lake microbiome and a wide array of environmental variables for lake ecosystems. Building on the results and methodology presented here, it will be possible to identify microbial taxa and processes that are essential for functioning and stability of lake ecosystems.

Modelling spatial patterns in host-associated microbial communities

Microbial communities exhibit spatial structure at different scales, due to constant interactions with their environment and dispersal limitation. While this spatial structure is often considered in studies focusing on free-living environmental communities, it has received less attention in the context of host-associated microbial communities or microbiota. The wider adoption of methods accounting for spatial variation in these communities will help to address open questions in basic microbial ecology as well as realize the full potential of microbiome-aided medicine. Here, we first overview known factors affecting the composition of microbiota across diverse host types and at different scales, with a focus on the human gut as one of the most actively studied microbiota. We outline a number of topical open questions in the field related to spatial variation and patterns. We then review the existing methodology for the spatial modelling of microbiota. We suggest that methodology from related fields, such as systems biology and macro-organismal ecology, could be adapted to obtain more accurate models of spatial structure. We further posit that methodological developments in the spatial modelling and analysis of microbiota could in turn broadly benefit theoretical and applied ecology and contribute to the development of novel industrial and clinical applications.

Debar: A sequence-by-sequence denoiser for COI-5P DNA barcode data

DNA barcoding and metabarcoding are now widely used to advance species discovery and biodiversity assessments. High-throughput sequencing (HTS) has expanded the volume and scope of these analyses, but elevated error rates introduce noise into sequence records that can inflate estimates of biodiversity. Denoising -the separation of biological signal from instrument (technical) noise-of barcode and metabarcode data currently employs abundance-based methods which do not capitalize on the highly conserved structure of the cytochrome c oxidase subunit I (COI) region employed as the animal barcode. This manuscript introduces debar, an R package that utilizes a profile hidden Markov model to denoise indel errors in COI sequences introduced by instrument error. In silico studies demonstrated that debar recognized 95% of artificially introduced indels in COI sequences. When applied to real-world data, debar reduced indel errors in circular consensus sequences obtained with the Sequel platform by 75%, and those generated on the Ion Torrent S5 by 94%. The false correction rate was less than 0.1%, indicating that debar is receptive to the majority of true COI variation in the animal kingdom. In conclusion, the debar package improves DNA barcode and metabarcode workflows by aiding the generation of more accurate sequences aiding the characterization of species diversity.

Microbiome Composition and Function in Aquatic Vertebrates: Small Organisms Making Big Impacts on Aquatic Animal Health

Aquatic ecosystems are under increasing stress from global anthropogenic and natural changes, including climate change, eutrophication, ocean acidification, and pollution. In this critical review, we synthesize research on the microbiota of aquatic vertebrates and discuss the impact of emerging stressors on aquatic microbial communities using two case studies, that of toxic cyanobacteria and microplastics. Most studies to date are focused on host-associated microbiomes of individual organisms, however, few studies take an integrative approach to examine aquatic vertebrate microbiomes by considering both host-associated and free-living microbiota within an ecosystem. We highlight what is known about microbiota in aquatic ecosystems, with a focus on the interface between water, fish, and marine mammals. Though microbiomes in water vary with geography, temperature, depth, and other factors, core microbial functions such as primary production, nitrogen cycling, and nutrient metabolism are often conserved across aquatic environments. We outline knowledge on the composition and function of tissue-specific microbiomes in fish and marine mammals and discuss the environmental factors influencing their structure. The microbiota of aquatic mammals and fish are highly unique to species and a delicate balance between respiratory, skin, and gastrointestinal microbiota exists within the host. In aquatic vertebrates, water conditions and ecological niche are driving factors behind microbial composition and function. We also generate a comprehensive catalog of marine mammal and fish microbial genera, revealing commonalities in composition and function among aquatic species, and discuss the potential use of microbiomes as indicators of health and ecological status of aquatic ecosystems. We also discuss the importance of a focus on the functional relevance of microbial communities in relation to organism physiology and their ability to overcome stressors related to global change. Understanding the dynamic relationship between aquatic microbiota and the animals they colonize is critical for monitoring water quality and population health.

Expanding ecological assessment by integrating microorganisms into routine freshwater biomonitoring

Bioindication has become an indispensable part of water quality monitoring in most countries of the world, with the presence and abundance of bioindicator taxa, mostly multicellular eukaryotes, used for biotic indices. In contrast, microbes (bacteria, archaea and protists) are seldom used as bioindicators in routine assessments, although they have been recognized for their importance in environmental processes. Recently, the use of molecular methods has revealed unexpected diversity within known functional groups and novel metabolic pathways that are particularly important in energy and nutrient cycling. In various habitats, microbial communities respond to eutrophication, metals, and natural or anthropogenic organic pollutants through changes in diversity and function. In this review, we evaluated the common trends in these changes, documenting that they have value as bioindicators and can be used not only for monitoring but also for improving our understanding of the major processes in lotic and lentic environments. Current knowledge provides a solid foundation for exploiting microbial taxa, community structures and diversity, as well as functional genes, in novel monitoring programs. These microbial community measures can also be combined into biotic indices, improving the resolution of individual bioindicators. Here, we assess particular molecular approaches complemented by advanced bioinformatic analysis, as these are the most promising with respect to detailed bioindication value. We conclude that microbial community dynamics are a missing link important for our understanding of rapid changes in the structure and function of aquatic ecosystems, and should be addressed in the future environmental monitoring of freshwater ecosystems.

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.

The Seagrass Holobiont: What We Know and What We Still Need to Disclose for Its Possible Use as an Ecological Indicator

Microbes and seagrass establish symbiotic relationships constituting a functional unit called the holobiont that reacts as a whole to environmental changes. Recent studies have shown that the seagrass microbial associated community varies according to host species, environmental conditions and the host’s health status, suggesting that the microbial communities respond rapidly to environmental disturbances and changes. These changes, dynamics of which are still far from being clear, could represent a sensitive monitoring tool and ecological indicator to detect early stages of seagrass stress. In this review, the state of art on seagrass holobiont is discussed in this perspective, with the aim of disentangling the influence of different factors in shaping it. As an example, we expand on the widely studied Halophila stipulacea’s associated microbial community, highlighting the changing and the constant components of the associated microbes, in different environmental conditions. These studies represent a pivotal contribution to understanding the holobiont’s dynamics and variability pattern, and to the potential development of ecological/ecotoxicological indices. The influences of the host’s physiological and environmental status in changing the seagrass holobiont, alongside the bioinformatic tools for data analysis, are key topics that need to be deepened, in order to use the seagrass-microbial interactions as a source of ecological information.

Machine learning applications in microbial ecology, human microbiome studies, and environmental monitoring

Advances in nucleic acid sequencing technology have enabled expansion of our ability to profile microbial diversity. These large datasets of taxonomic and functional diversity are key to better understanding microbial ecology. Machine learning has proven to be a useful approach for analyzing microbial community data and making predictions about outcomes including human and environmental health. Machine learning applied to microbial community profiles has been used to predict disease states in human health, environmental quality and presence of contamination in the environment, and as trace evidence in forensics. Machine learning has appeal as a powerful tool that can provide deep insights into microbial communities and identify patterns in microbial community data. However, often machine learning models can be used as black boxes to predict a specific outcome, with little understanding of how the models arrived at predictions. Complex machine learning algorithms often may value higher accuracy and performance at the sacrifice of interpretability. In order to leverage machine learning into more translational research related to the microbiome and strengthen our ability to extract meaningful biological information, it is important for models to be interpretable. Here we review current trends in machine learning applications in microbial ecology as well as some of the important challenges and opportunities for more broad application of machine learning to understanding microbial communities.