Genomics in marine monitoring: new opportunities for assessing marine health status

Forever young: The successful story of a marine biotic index

In 2000, the AZTI's Marine Biotic Index (AMBI) was published and was one of a number of marine benthic indices development to assess the ecological status of soft-bottom macroinvertebrates. This index, and its derivatives, has been very successful in its application to different geographical areas, across the world, as well as to different environments, from the intertidal to the abyssal, or from tidal freshwater to offshore habitats. In this review, we explain the story behind the AMBI development, and look for an explanation of the index's success. For doing that, we comment on the current practicalities of the index, we present the new AMBI species list, with 9251 taxa, we dismantle six myths around this index, and examine the past and the future of the index. We show that the solid ecological roots of the index make it a robust tool to assess the status of marine benthic communities, at any time and environment. Hence, we think that it will be 'forever young' helping managers in taking informed decisions to improve benthic marine systems.

Harmful phytoplankton diversity and dynamics in an upwelling region (Sagres, SW Portugal) revealed by ribosomal RNA microarray combined with microscopy

The study region in Sagres, SW Portugal, is subject to natural eutrophication of coastal waters by wind-driven upwelling, which stimulates high primary productivity facilitating the recent economic expansion of bivalve aquaculture in the region. However, this economic activity is threatened by harmful algal blooms (HAB) caused by the diatoms Pseudo-nitzschia spp., Dinophysis spp. and other HAB dinoflagellates, all of which can produce toxins, that can induce Amnesic Shellfish Poisoning (ASP), Diarrhetic Shellfish Poisoning (DSP) and Paralytic Shellfish Poisoning (PSP). This study couples traditional microscopy with 18S/28S rRNA microarray to improve the detection of HAB species and investigates the relation between HAB and the specific oceanographic conditions in the region. Good agreement was obtained between microscopy and microarray data for diatoms of genus Pseudo-nitzschia and dinoflagellates Dinophysis spp., Gymnodinium catenatum and raphidophyte Heterosigma akashiwo, with less effective results for Prorocentrum. Microarray provided detection of flagellates Prymnesium spp., Pseudochattonella spp., Chloromorum toxicum and the important HAB dinoflagellates of the genera Alexandrium and Azadinium, with the latter being one of the first records from the study region. Seasonality and upwelling induced by northerly winds were found to be the driving forces of HAB development, with Pseudo-nitzschia spp. causing the risk of ASP during spring and summer upwelling season, and dinoflagellates causing the risk of DSP and PSP during upwelling relaxation, mainly in summer and autumn. The findings were in agreement with the results from toxicity monitoring of shellfish by the Portuguese Institute for Sea and Atmosphere and confirm the suitability of the RNA microarray method for HABs detection and aquaculture management applications.

Towards an eDNA metabarcode-based performance indicator for full-scale municipal wastewater treatment plants

Knowledge-driven management for wastewater treatment plant (WWTP) performance becomes increasingly important considering the globally growing production of wastewater and the rising demand of clean water supply. Even though the potential of microbial organisms (bacteria and protists) as bioindicators for WWTP performance is well known, it is far from being fully exploited for routine monitoring programs. Therefore, we here used massive sequencing of environmental (e)DNA metabarcodes from bacterial (V3-V4 region of the SSU rRNA gene) and eukaryote (V9 region of the SSU rRNA gene) communities in 21 activated sludge samples obtained from full-scale municipal WWTPs in Germany. Microbial community patterns were compared to standard WWTP operating parameters and two traditionally used WWTP performance indicators (Sludge Biotic Index and Sludge Index). Both indices showed low concordance and hardly correlated with chemical WWTP performance parameters nor did they correlate with microbial community structures. In contrast, microbial community profiles significantly correlated with WWTP performance parameters and operating conditions of the plants under study. Therefore, eDNA metabarcode profiles of whole microbial communities indicate the performance of WWTP and can provide useful information for management strategies. We here suggest a strategy for the development of an eDNA metabarcode based bioindicator system, which can be implemented in future standard monitoring programs for WWTP performance and effluent quality.

Toxicology tailored low density oligonucleotide microarray for the thicklip grey mullets (Chelon labrosus): Biomarker gene transcription profile after caging in a polluted harbour

In aquatic organisms inhabiting polluted waters genes are activated to build an adaptive/compensatory defence against the possible effects of pollutants. Such responses can be used as biomarkers of exposure to chemical compounds, outlining the molecular mechanisms activated under specific pollution scenarios. With the aim of exploiting such approach in environmental health assessment, toxicologically relevant gene fragments were sequenced in the thicklip grey mullet (Chelon labrosus) and a toxicologically tailored low-density (160 genes) oligonucleotide microarray was customised. The tool was validated comparing organ/sex specific gene expression profiles and characterising responses under laboratory exposure to model chemicals. Finally, juvenile mullets were caged in a polluted harbour and hepatic gene expression profiles analysed after 5 and 21 days of deployment. Cages were deployed in the inner (IH) and outer (OH) Pasaia harbour, Bay of Biscay. Mussels (Mytilus galloprovincialis) were also caged as biological matrix for chemical bioaccumulation analysis and stress biomarkers measurements. Slightly higher concentrations of chemicals (metals, tributyltin, PAHs, phthalates) were quantified in IH than in OH, fish bile metabolites also revealing higher availability of PAHs in IH. Lysosome membrane stability in mussels was reduced, indicating stress condition in both sites. The developed microarray discriminated mullets showing distinctive expression profiles depending on site and deployment time. Genes related to immune and hypoxia responses were regulated comparing IH and OH at day 5. Phase I and II biotransformation genes, such as cyp2, cyp3 and ugt, were up-regulated in IH, together with the aryl hydrocarbon receptor 2 (ahr2) and the ahr repressor. Similarly, TBT-binding proteins and genes involved in lipid metabolism (pparγ, cyp7) were up-regulated with deployment time. Even if nowadays higher throughput approaches for gene expression analyses are available, the developed mullet tool constitutes a comprehensive tool to assess molecular responses of mullets exposed to pollutants, although it remains to be explored whether it can be applied to assess pollutant exposure in active pollution monitorings and in environmental health assessment.

A comparative analysis of metabarcoding and morphology-based identification of benthic communities across different regional seas

In a world of declining biodiversity, monitoring is becoming crucial. Molecular methods, such as metabarcoding, have the potential to rapidly expand our knowledge of biodiversity, supporting assessment, management, and conservation. In the marine environment, where hard substrata are more difficult to access than soft bottoms for quantitative ecological studies, Artificial Substrate Units (ASUs) allow for standardized sampling. We deployed ASUs within five regional seas (Baltic Sea, Northeast Atlantic Ocean, Mediterranean Sea, Black Sea, and Red Sea) for 12-26 months to measure the diversity and community composition of macroinvertebrates. We identified invertebrates using a traditional approach based on morphological characters, and by metabarcoding of the mitochondrial cytochrome oxidase I (COI) gene. We compared community composition and diversity metrics obtained using the two methods. Diversity was significantly correlated between data types. Metabarcoding of ASUs allowed for robust comparisons of community composition and diversity, but not all groups were successfully sequenced. All locations were significantly different in taxonomic composition as measured with both kinds of data. We recovered previously known regional biogeographical patterns in both datasets (e.g., low species diversity in the Black and Baltic Seas, affinity between the Bay of Biscay and the Mediterranean). We conclude that the two approaches provide complementary information and that metabarcoding shows great promise for marine monitoring. However, until its pitfalls are addressed, the use of metabarcoding in monitoring of rocky benthic assemblages should be used in addition to classical approaches rather than instead of them.

Multi-marker metabarcoding approach to study mesozooplankton at basin scale

Zooplankton plays a pivotal role in marine ecosystems and the characterisation of its biodiversity still represents a challenge for marine ecologists. In this study, mesozooplankton composition from 46 samples collected in summer along the western Adriatic Sea, was retrieved by DNA metabarcoding analysis. For the first time, the highly variable fragments of the mtDNA COI and the V9 region of 18S rRNA genes were used in a combined matrix to compile an inventory of mesozooplankton at basin scale. The number of sequences retrieved after quality filtering were 824,148 and 223,273 for COI and 18S (V9), respectively. The taxonomical assignment against reference sequences, using 95% (for COI) and 97% (for 18S) similarity thresholds, recovered 234 taxa. NMDS plots and cluster analysis divided coastal from offshore samples and the most representative species of these clusters were distributed according to the dominant surface current pattern of the Adriatic for the summer period. For selected sampling sites, mesozooplankton species were also identified under a stereo microscope providing insights on the strength and weakness of the two approaches. In addition, DNA metabarcoding was shown to be helpful for the monitoring of non-indigenous marine metazoans and spawning areas of commercial fish species. We defined pros and cons of applying this approach at basin scale and the benefits of combining the datasets from two genetic markers.

Implementation options for DNA-based identification into ecological status assessment under the European Water Framework Directive

Assessment of ecological status for the European Water Framework Directive (WFD) is based on "Biological Quality Elements" (BQEs), namely phytoplankton, benthic flora, benthic invertebrates and fish. Morphological identification of these organisms is a time-consuming and expensive procedure. Here, we assess the options for complementing and, perhaps, replacing morphological identification with procedures using eDNA, metabarcoding or similar approaches. We rate the applicability of DNA-based identification for the individual BQEs and water categories (rivers, lakes, transitional and coastal waters) against eleven criteria, summarised under the headlines representativeness (for example suitability of current sampling methods for DNA-based identification, errors from DNA-based species detection), sensitivity (for example capability to detect sensitive taxa, unassigned reads), precision of DNA-based identification (knowledge about uncertainty), comparability with conventional approaches (for example sensitivity of metrics to differences in DNA-based identification), cost effectiveness and environmental impact. Overall, suitability of DNA-based identification is particularly high for fish, as eDNA is a well-suited sampling approach which can replace expensive and potentially harmful methods such as gill-netting, trawling or electrofishing. Furthermore, there are attempts to replace absolute by relative abundance in metric calculations. For invertebrates and phytobenthos, the main challenges include the modification of indices and completing barcode libraries. For phytoplankton, the barcode libraries are even more problematic, due to the high taxonomic diversity in plankton samples. If current assessment concepts are kept, DNA-based identification is least appropriate for macrophytes (rivers, lakes) and angiosperms/macroalgae (transitional and coastal waters), which are surveyed rather than sampled. We discuss general implications of implementing DNA-based identification into standard ecological assessment, in particular considering any adaptations to the WFD that may be required to facilitate the transition to molecular data.

Evaluation of biomarkers in Mytilus galloprovincialis as an integrated measure of biofilm-membrane bioreactor (BF-MBR) system efficiency in mitigating the impact of oily wastewater discharge to marine environment: a microcosm approach

The large volumes of oily wastewater discharged to marine environment cause heavy impacts on the coastal marine ecosystem. The selection of an appropriate technology to reduce these impacts should be based on the respect of the discharge limits and on the effective assessment and monitoring of its effects on biological organism preservation. To this aim, we set up a controlled microcosm-scale system to compare the effects of a treated and untreated oily wastewater discharge in which the restore process is performed through a Membrane Bio-Reactor. The system is completed by other three microcosms to control and isolate any possible concurrent effect on the Mytilus galloprovincialis, used as sentinel organism. Mytilus galloprovincialis have been kept in all these microcosms, and then mRNA expression and morphology were evaluated on gills and digestive gland. The genes considered in this work are Heat Shock Protein 70 and Metallothionein 10, involved in response to physicochemical sublethal stressors, Superoxide dismutase 1, Catalase, and Cytochrome P450 involved in oxidative stress response. Our results evidenced a significant overexpression, both in gills and digestive gland, of HSP70 in samples maintained in the microcosm receiving the untreated effluent, and of MT10 in those animals kept in microcosm where the effluent was treated. Even though the mRNA modifications are considered "primary" and transient responses which do not always correspond to protein content, the study of these modifications can help to gain insights into the mechanisms of action of xenobiotic exposure. Morphological analysis suggested that, although different, depending on the microcosm, the most serious damages were found in the gill epithelium accompanied with severe haemocyte infiltration, whilst in digestive gland the tissue architecture alterations and the haemocyte infiltration were less pronounced. These observations suggest that the immune system was activated as a general response to stressful stimuli such as the presence of toxic compounds. Moreover, the results indicate that the treatment process is useful. In fact, samples derived from the microcosm receiving the treated effluent, even though presenting some signs of stress, seemed to partially recover the normal structure, although their mRNA expression indicated some cellular suffering.

Analysis of bacterial metagenomes from the Southwestern Gulf of Mexico for pathogens detection

Little is known about the diversity of bacteria in the Southwestern Gulf of Mexico. The aim of the study illustrated in this perspective was to search for the presence of bacterial pathogens in this ecosystem, using metagenomic data recently generated by the Mexican research group known as the Gulf of Mexico Research Consortium. Several genera of bacteria annotated as pathogens were detected in water and sediment marine samples. As expected, native and ubiquitous pathogenic bacteria genera such as Burkolderia, Halomonas, Pseudomonas, Shewanella and Vibrio were highly represented. Surprisingly, non-native genera of public health concern were also detected, including Borrelia, Ehrlichia, Leptospira, Mycobacterium, Mycoplasma, Salmonella, Staphylococcus, Streptococcus and Treponema. While there are no previous metagenomics studies of this environment, the potential influences of natural, anthropogenic and ecological factors on the diversity of putative pathogenic bacteria found in it are reviewed. The taxonomic annotation herein reported provides a starting point for an improved understanding of bacterial biodiversity in the Southwestern Gulf of Mexico. It also represents a useful tool in public health as it may help identify infectious diseases associated with exposure to marine water and ingestion of fish or shellfish, and thus may be useful in predicting and preventing waterborne disease outbreaks.

A cross-taxa study using environmental DNA/RNA metabarcoding to measure biological impacts of offshore oil and gas drilling and production operations

Standardized ecosystem-based monitoring surveys are critical for providing information on marine ecosystem health. Environmental DNA/RNA (eDNA/eRNA) metabarcoding may facilitate such surveys by quickly and effectively characterizing multi-trophic levels. In this study, we assessed the suitability of eDNA/eRNA metabarcoding to evaluate changes in benthic assemblages of bacteria, Foraminifera and other eukaryotes along transects at three offshore oil and gas (O&G) drilling and production sites, and compared these to morphologically characterized macro-faunal assemblages. Bacterial communities were the most responsive to O&G activities, followed by Foraminifera, and macro-fauna (the latter assessed by morphology). The molecular approach enabled detection of hydrocarbon degrading taxa such as the bacteria Alcanivorax and Microbulbifer at petroleum impacted stations. Most identified indicator taxa, notably among macro-fauna, were highly specific to site conditions. Based on our results we suggest that eDNA/eRNA metabarcoding can be used as a stand-alone method for biodiversity assessment or as a complement to morphology-based monitoring approaches.

Early detection monitoring for aquatic non-indigenous species: Optimizing surveillance, incorporating advanced technologies, and identifying research needs

Following decades of ecologic and economic impacts from a growing list of nonindigenous and invasive species, government and management entities are committing to systematic early- detection monitoring (EDM). This has reinvigorated investment in the science underpinning such monitoring, as well as the need to convey that science in practical terms to those tasked with EDM implementation. Using the context of nonindigenous species in the North American Great Lakes, this article summarizes the current scientific tools and knowledge - including limitations, research needs, and likely future developments - relevant to various aspects of planning and conducting comprehensive EDM. We begin with the scope of the effort, contrasting target-species with broad-spectrum monitoring, reviewing information to support prioritization based on species and locations, and exploring the challenge of moving beyond individual surveys towards a coordinated monitoring network. Next, we discuss survey design, including effort to expend and its allocation over space and time. A section on sample collection and analysis overviews the merits of collecting actual organisms versus shed DNA, reviews the capabilities and limitations of identification by morphology, DNA target markers, or DNA barcoding, and examines best practices for sample handling and data verification. We end with a section addressing the analysis of monitoring data, including methods to evaluate survey performance and characterize and communicate uncertainty. Although the body of science supporting EDM implementation is already substantial, research and information needs (many already actively being addressed) include: better data to support risk assessments that guide choice of taxa and locations to monitor; improved understanding of spatiotemporal scales for sample collection; further development of DNA target markers, reference barcodes, genomic workflows, and synergies between DNA-based and morphology-based taxonomy; and tools and information management systems for better evaluating and communicating survey outcomes and uncertainty.

Indigenous species barcode database improves the identification of zooplankton

Incompleteness and inaccuracy of DNA barcode databases is considered an important hindrance to the use of metabarcoding in biodiversity analysis of zooplankton at the species-level. Species barcoding by Sanger sequencing is inefficient for organisms with small body sizes, such as zooplankton. Here mitochondrial cytochrome c oxidase I (COI) fragment barcodes from 910 freshwater zooplankton specimens (87 morphospecies) were recovered by a high-throughput sequencing platform, Ion Torrent PGM. Intraspecific divergence of most zooplanktons was < 5%, except Branchionus leydign (Rotifer, 14.3%), Trichocerca elongate (Rotifer, 11.5%), Lecane bulla (Rotifer, 15.9%), Synchaeta oblonga (Rotifer, 5.95%) and Schmackeria forbesi (Copepod, 6.5%). Metabarcoding data of 28 environmental samples from Lake Tai were annotated by both an indigenous database and NCBI Genbank database. The indigenous database improved the taxonomic assignment of metabarcoding of zooplankton. Most zooplankton (81%) with barcode sequences in the indigenous database were identified by metabarcoding monitoring. Furthermore, the frequency and distribution of zooplankton were also consistent between metabarcoding and morphology identification. Overall, the indigenous database improved the taxonomic assignment of zooplankton.

Interactive effects of nutrition, reproductive state and pollution on molecular stress responses of mussels, Mytilus galloprovincialis Lamarck, 1819

Marine bivalves including mussels Mytilus galloprovincialis are commonly used as sentinels for pollution monitoring and ecosystem health assessment in the coastal zones. Use of biomarkers to assess the pollution effects assumes that the effects of pollutants on the biomarkers exceed the natural background variability; yet this assumption has rarely been tested. We exposed mussels at different reproductive stages and nutritive states to two concentrations of a polycyclic aromatic hydrocarbon (fluoranthene, 3 and 60 μg L^-1) for three weeks. Expression levels of the molecular biomarkers related to the detoxification and general stress response [cytochrome P450 oxidase (CYP450), glutathione S-transferases (GST-α; GST-S1; GST-S2), the multixenobiotic resistance protein P-glycoprotein (PgP), metallothioneins (MT10 and MT20), heat shock proteins (HSP22, HSP70-2; HSP70-3; HSP70-4), as well as mRNA expression of two reproduction-related genes, vitellogenin (Vitel) and vitelline coat lysin M7 (VCLM7)] were measured. The mussels' nutrition and reproductive state affected the baseline mRNA levels of molecular biomarkers and modulated the transcriptional responses of biomarker genes to the pollutant exposure. Thus, mussel physiological state could act as a confounding factor in the evaluation of the response of pollution through molecular biomarkers. The biomarker baseline levels must be determined across a range of physiological states to enable the use of biomarkers in monitoring programs.

Persistence of marine fish environmental DNA and the influence of sunlight

Harnessing information encoded in environmental DNA (eDNA) in marine waters has the potential to revolutionize marine biomonitoring. Whether using organism-specific quantitative PCR assays or metabarcoding in conjunction with amplicon sequencing, scientists have illustrated that realistic organism censuses can be inferred from eDNA. The next step is establishing ways to link information obtained from eDNA analyses to actual organism abundance. This is only possible by understanding the processes that control eDNA concentrations. The present study uses mesocosm experiments to study the persistence of eDNA in marine waters and explore the role of sunlight in modulating eDNA persistence. We seeded solute-permeable dialysis bags with water containing indigenous eDNA and suspended them in a large tank containing seawater. Bags were subjected to two treatments: half the bags were suspended near the water surface where they received high doses of sunlight, and half at depth where they received lower doses of sunlight. Bags were destructively sampled over the course of 87 hours. eDNA was extracted from water samples and used as template for a Scomber japonicus qPCR assay and a marine fish-specific 12S rRNA PCR assay. The latter was subsequently sequenced using a metabarcoding approach. S. japonicus eDNA, as measured by qPCR, exhibited first order decay with a rate constant ~0.01 hr ^-1 with no difference in decay rate constants between the two experimental treatments. eDNA metabarcoding identified 190 organizational taxonomic units (OTUs) assigned to varying taxonomic ranks. There was no difference in marine fish communities as measured by eDNA metabarcoding between the two experimental treatments, but there was an effect of time. Given the differences in UVA and UVB fluence received by the two experimental treatments, we conclude that sunlight is not the main driver of fish eDNA decay in the experiments. However, there are clearly temporal effects that need to be considered when interpreting information obtained using eDNA approaches.

Multiple stressors in sediments impact adjacent hard substrate habitats and across biological domains

Coastal systems are increasingly impacted by human activities. While the direct effects of individual contaminants have been investigated, the potential for multiple contaminants to impact adjacent hard substrate habitats is poorly understood. Sediment-bound contaminants pose a risk to water column organisms through resuspension and the fluxing of dissolved nutrients and metals. This study experimentally manipulated contaminated coastal sediments in mesocosms with additions of a common fertiliser to investigate the impact on both bacterial biofilms and macrofouling communities on nearby hard substrates. Field mesocosms were deployed sub-tidally for two weeks in a fully crossed design with two levels of metal contamination (ambient or high) and three levels of organic enrichment (ambient, low and high). Developing biofilm and macrofaunal communities were collected on acetate settlement sheets above the mesocosm sediments and censused with a combination of high-throughput sequencing (biofilm) and microscopy (macrofauna). Organic enrichment of sediments induced compositional shifts in biofilm communities, reducing their diversity, evenness and richness. Furthermore, co-occurrence networks built from microbial assemblages exposed to contaminated sediments displayed reduced connectivity compared to controls, suggesting a more stochastic assembly dynamic, where microbial interactions are reduced. Macrofouling community composition shifted in response to increased enrichment with separate and interactive effects of metals also observed for individual taxa. Specifically, antagonistic stressor interactions were observed for colonial ascidians and arborescent bryozoans; metal contamination decreased abundances of these taxa, except under high enrichment conditions. Together these micro- and macrofaunal responses indicate selection for depauperate, but contaminant-tolerant, communities and a potential breakdown in biotic connectivity through multiple stressor impacts across habitat boundaries.

DNA extraction replicates improve diversity and compositional dissimilarity in metabarcoding of eukaryotes in marine sediments

Human impact on marine benthic communities has traditionally been assessed using visible morphological traits and has focused on the macrobenthos, whereas the ecologically important organisms of the meio- and microbenthos have received less attention. DNA metabarcoding offers an alternative to this approach and enables a larger fraction of the biodiversity in marine sediments to be monitored in a cost-efficient manner. Although this methodology remains poorly standardised and challenged by biases inherent to rRNA copy number variation, DNA extraction, PCR, and limitations related to taxonomic identification, it has been shown to be semi-quantitative and useful for comparing taxon abundances between samples. Here, we evaluate the effect of replicating genomic DNA extraction in order to counteract small scale spatial heterogeneity and improve diversity and community structure estimates in metabarcoding-based monitoring. For this purpose, we used ten technical replicates from three different marine sediment samples. The effect of sequence depth was also assessed, and in silico pooling of DNA extraction replicates carried out in order to maintain the number of reads constant. Our analyses demonstrated that both sequencing depth and DNA extraction replicates could improve diversity estimates as well as the ability to separate samples with different characteristics. We could not identify a “sufficient” replicate number or sequence depth, where further improvements had a less significant effect. Based on these results, we consider replication an attractive alternative to directly increasing the amount of sample used for DNA extraction and strongly recommend it for future metabarcoding studies and routine assessments of sediment biodiversity.

Metabarcoding monitoring analysis: the pros and cons of using co-extracted environmental DNA and RNA data to assess offshore oil production impacts on benthic communities

Sequencing environmental DNA (eDNA) is increasingly being used as an alternative to traditional morphological-based identification to characterize biological assemblages and monitor anthropogenic impacts in marine environments. Most studies only assess eDNA which, compared to eRNA, can persist longer in the environment after cell death. Therefore, eRNA may provide a more immediate census of the environment due to its relatively weaker stability, leading some researchers to advocate for the use of eRNA as an additional, or perhaps superior proxy for portraying ecological changes. A variety of pre-treatment techniques for screening eDNA and eRNA derived operational taxonomic units (OTUs) have been employed prior to statistical analyses, including removing singleton taxa (i.e., OTUs found only once) and discarding those not present in both eDNA and eRNA datasets. In this study, we used bacterial (16S ribosomal RNA gene) and eukaryotic (18S ribosomal RNA gene) eDNA- and eRNA-derived data from benthic communities collected at increasing distances along a transect from an oil production platform (Taranaki, New Zealand). Macro-infauna (visual classification of benthic invertebrates) and physico-chemical data were analyzed in parallel. We tested the effect of removing singleton taxa, and removing taxa not present in the eDNA and eRNA libraries from the same environmental sample (trimmed by shared OTUs), by comparing the impact of the oil production platform on alpha- and beta-diversity of the eDNA/eRNA-based biological assemblages, and by correlating these to the morphologically identified macro-faunal communities and the physico-chemical data. When trimmed by singletons, presence/absence information from eRNA data represented the best proxy to detect changes on species diversity for both bacteria and eukaryotes. However, assessment of quantitative beta-diversity from read abundance information of bacteria eRNA did not, contrary to eDNA, reveal any impact from the oil production activity. Overall, the data appeared more robust when trimmed by shared OTUs, showing a greater effect of the platform on alpha- and beta-diversity. Trimming by shared OTUs likely removes taxa derived from legacy DNA and technical artefacts introduced through reverse transcriptase, polymerase-chain-reaction and sequencing. Findings from our scoping study suggest that metabarcoding-based biomonitoring surveys should, if funds, time and expertise allow, be assessed using both eDNA and eRNA products.

Reconstructing a herbivore's diet using a novel rbcL DNA mini-barcode for plants

Next Generation Sequencing and the application of metagenomic analyses can be used to answer questions about animal diet choice and study the consequences of selective foraging by herbivores. The quantification of herbivore diet choice with respect to native versus exotic plant species is particularly relevant given concerns of invasive species establishment and their effects on ecosystems. While increased abundance of white-tailed deer (Odocoileus virginianus) appears to correlate with increased incidence of invasive plant species, data supporting a causal link is scarce. We used a metabarcoding approach (PCR amplicons of the plant rbcL gene) to survey the diet of white-tailed deer (fecal samples), from a forested site in Warren County, Virginia with a comprehensive plant species inventory and corresponding reference collection of plant barcode and chloroplast sequences. We sampled fecal pellet piles and extracted DNA from 12 individual deer in October 2014. These samples were compared to a reference DNA library of plant species collected within the study area. For 72 % of the amplicons, we were able to assign taxonomy at the species level, which provides for the first time-sufficient taxonomic resolution to quantify the relative frequency at which native and exotic plant species are being consumed by white-tailed deer. For each of the 12 individual deer we collected three subsamples from the same fecal sample, resulting in sequencing 36 total samples. Using Qiime, we quantified the plant DNA found in all 36 samples, and found that variance within samples was less than variance between samples (F = 1.73, P = 0.004), indicating additional subsamples may not be necessary. Species level diversity ranged from 60 to 93 OTUs per individual and nearly 70 % of all plant sequences recovered were from native plant species. The number of species detected did reduce significantly (range 4-12) when we excluded species whose OTU composed <1 % of each sample's total. When compared to the abundance of native and non-natives plants inventoried in the local community, our results support the observation that white-tailed deer have strong foraging preferences, but these preferences were not consistent for species in either class. Deer forage behaviour may favour some exotic species, but not all.

Aquatic environmental DNA detects seasonal fish abundance and habitat preference in an urban estuary

The difficulty of censusing marine animal populations hampers effective ocean management. Analyzing water for DNA traces shed by organisms may aid assessment. Here we tested aquatic environmental DNA (eDNA) as an indicator of fish presence in the lower Hudson River estuary. A checklist of local marine fish and their relative abundance was prepared by compiling 12 traditional surveys conducted between 1988–2015. To improve eDNA identification success, 31 specimens representing 18 marine fish species were sequenced for two mitochondrial gene regions, boosting coverage of the 12S eDNA target sequence to 80% of local taxa. We collected 76 one-liter shoreline surface water samples at two contrasting estuary locations over six months beginning in January 2016. eDNA was amplified with vertebrate-specific 12S primers. Bioinformatic analysis of amplified DNA, using a reference library of GenBank and our newly generated 12S sequences, detected most (81%) locally abundant or common species and relatively few (23%) uncommon taxa, and corresponded to seasonal presence and habitat preference as determined by traditional surveys. Approximately 2% of fish reads were commonly consumed species that are rare or absent in local waters, consistent with wastewater input. Freshwater species were rarely detected despite Hudson River inflow. These results support further exploration and suggest eDNA will facilitate fine-scale geographic and temporal mapping of marine fish populations at relatively low cost.

Integrated chemical and biological assessment of contaminant impacts in selected European coastal and offshore marine areas

This paper reports a full assessment of results from ICON, an international workshop on marine integrated contaminant monitoring, encompassing different matrices (sediment, fish, mussels, gastropods), areas (Iceland, North Sea, Baltic, Wadden Sea, Seine estuary and the western Mediterranean) and endpoints (chemical analyses, biological effects). ICON has demonstrated the use of a framework for integrated contaminant assessment on European coastal and offshore areas. The assessment showed that chemical contamination did not always correspond with biological effects, indicating that both are required. The framework can be used to develop assessments for EU directives. If a 95% target were to be used as a regional indicator of MSFD GES, Iceland and offshore North Sea would achieve the target using the ICON dataset, but inshore North Sea, Baltic and Spanish Mediterranean regions would fail.

The best of both worlds: A combined approach for analyzing microalgal diversity via metabarcoding and morphology-based methods

An increasing number of studies use next generation sequencing (NGS) to analyze complex communities, but is the method sensitive enough when it comes to identification and quantification of species? We compared NGS with morphology-based identification methods in an analysis of microalgal (periphyton) communities. We conducted a mesocosm experiment in which we allowed two benthic grazer species to feed upon benthic biofilms, which resulted in altered periphyton communities. Morphology-based identification and 454 (Roche) pyrosequencing of the V4 region in the small ribosomal unit (18S) rDNA gene were used to investigate the community change caused by grazing. Both the NGS-based data and the morphology-based method detected a marked shift in the biofilm composition, though the two methods varied strongly in their abilities to detect and quantify specific taxa, and neither method was able to detect all species in the biofilms. For quantitative analysis, we therefore recommend using both metabarcoding and microscopic identification when assessing the community composition of eukaryotic microorganisms.

Hepatic gene transcription profiles in turbot (Scophthalmus maximus) experimentally exposed to heavy fuel oil nº 6 and to styrene

Oil and chemical spills in the marine environment, although sporadic, are highly dangerous to biota inhabiting coastal and estuarine areas. Effects of spilled compounds in exposed organisms occur at different biological organization levels: from molecular, cellular or tissue levels to the physiological one. The present study aims to determine the specific hepatic gene transcription profiles observed in turbot juveniles under exposure to fuel oil n °6 and styrene vs controls using an immune enriched turbot (Scophthalmus maximus) oligo-microarray containing 2716 specific gene probes. After 3 days of exposure, fuel oil specifically induced aryl hydrocarbon receptor mediated transcriptional response through up-regulation of genes, such as ahrr and cyp1a1. More gene transcripts were regulated after 14 days of exposure involved in ribosomal biosynthesis, immune modulation, and oxidative response among the most significantly regulated functional pathways. On the contrary, gene transcription alterations caused by styrene did not highlight any significantly regulated molecular or metabolic pathway. This was also previously reported at cell and tissue level where no apparent responses were distinguishable. For the fuel oil experiment, obtained specific gene profiles could be related to changes in cell-tissue organization in the same individuals, such as increased hepatocyte vacuolization, decrease in melano-macrophage centers and the regulation of leukocyte numbers. In conclusion, the mode of action reflected by gene transcription profiles analyzed hereby in turbot livers could be linked with the responses previously reported at higher biological organization levels. Molecular alterations described hereby could be preceding observed alterations at cell and tissue levels.

Exploring the Impacts of Anthropogenic Disturbance on Seawater and Sediment Microbial Communities in Korean Coastal Waters Using Metagenomics Analysis

The coastal ecosystems are considered as one of the most dynamic and vulnerable environments under various anthropogenic developments and the effects of climate change. Variations in the composition and diversity of microbial communities may be a good indicator for determining whether the marine ecosystems are affected by complex forcing stressors. DNA sequence-based metagenomics has recently emerged as a promising tool for analyzing the structure and diversity of microbial communities based on environmental DNA (eDNA). However, few studies have so far been performed using this approach to assess the impacts of human activities on the microbial communities in marine systems. In this study, using metagenomic DNA sequencing (16S ribosomal RNA gene), we analyzed and compared seawater and sediment communities between sand mining and control (natural) sites in southern coastal waters of Korea to assess whether anthropogenic activities have significantly affected the microbial communities. The sand mining sites harbored considerably lower levels of microbial diversities in the surface seawater community during spring compared with control sites. Moreover, the sand mining areas had distinct microbial taxonomic group compositions, particularly during spring season. The microbial groups detected solely in the sediment load/dredging areas (e.g., Marinobacter, Alcanivorax, Novosphingobium) are known to be involved in degradation of toxic chemicals such as hydrocarbon, oil, and aromatic compounds, and they also contain potential pathogens. This study highlights the versatility of metagenomics in monitoring and diagnosing the impacts of human disturbance on the environmental health of marine ecosystems from eDNA.

De novo transcriptome assembly of the marine gastropod Reishia clavigera for supporting toxic mechanism studies

The intertidal whelk Reishia clavigera is commonly used as a biomonitor of chemical contamination in the marine environment along Western Pacific region, and as a model for mechanistic studies of organotin-mediated imposex development. However, limited genomic resources of R. clavigera have restricted its role for the investigation of molecular mechanisms of such endocrine disruptions. This study, therefore, aimed to establish tissue-specific transcriptomes of the digestive gland, gonad, head ganglia, penis and the remaining body part of the male and female R. clavigera. By combining the results, a global transcriptome was obtained. A total of 578,134,720 high-quality filtered reads were obtained using Illumina sequencing. The R. clavigera transcriptome comprised of 38,466 transcripts and 32,798 unigenes with predicted open reading frames. The average length of transcripts was 1,709bp with N50 of 2,236bp. Based on sequence similarity searches against public databases, 28,657 transcripts and 24,403 unigenes had at least one BLAST hit. There were 17,530 transcripts and 14,897 unigenes annotated with at least one Gene Ontology (GO) term. Moreover, 5,776 transcripts and 5,137 unigenes were associated with 333 Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathways. The numbers of unigenes were similar among the five target tissues and between sexes, but tissue-specific expression profiles were revealed by multivariate analyses. Based on the functional annotation, putative steroid hormone-associated unigenes were identified. In particular, we highlighted the presence of steroid hormone receptor homologues that could be the targets for mechanistic studies of the organotin-mediated imposex development in marine gastropods. This newly generated transcriptome assembly of R. clavigera provides a valuable molecular resource for ecotoxicological and environmental genomic studies.

First evaluation of foraminiferal metabarcoding for monitoring environmental impact from an offshore oil drilling site

At present, environmental impacts from offshore oil and gas activities are partly determined by measuring changes in macrofauna diversity. Morphological identification of macrofauna is time-consuming, expensive and dependent on taxonomic expertise. In this study, we evaluated the applicability of using foraminiferal-specific metabarcoding for routine monitoring. Sediment samples were collected along distance gradients from two oil platforms off Taranaki (New Zealand) and their physico-chemical properties, foraminiferal environmental DNA/RNA, and macrofaunal composition analyzed. Macrofaunal and foraminiferal assemblages showed similar shifts along impact gradients, but responded differently to environmental perturbations. Macrofauna were affected by hypoxia, whereas sediment grain size appeared to drive shifts in foraminifera. We identified eight foraminiferal molecular operational taxonomic units that have potential to be used as bioindicator taxa. Our results show that metabarcoding represents an effective tool for assessing foraminiferal communities near offshore oil and gas platforms, and that it can be used to complement current monitoring techniques.

Monitoring and managing microbes in aquaculture - Towards a sustainable industry

Microorganisms are of great importance to aquaculture where they occur naturally, and can be added artificially, fulfilling different roles. They recycle nutrients, degrade organic matter and, occasionally, they infect and kill the fish, their larvae or the live feed. Also, some microorganisms may protect fish and larvae against disease. Hence, monitoring and manipulating the microbial communities in aquaculture environments hold great potential; both in terms of assessing and improving water quality, but also in terms of controlling the development of microbial infections. Using microbial communities to monitor water quality and to efficiently carry out ecosystem services within the aquaculture systems may only be a few years away. Initially, however, we need to thoroughly understand the microbiomes of both healthy and diseased aquaculture systems, and we need to determine how to successfully manipulate and engineer these microbiomes. Similarly, we can reduce the need to apply antibiotics in aquaculture through manipulation of the microbiome, i.e. by the use of probiotic bacteria. Recent studies have demonstrated that fish pathogenic bacteria in live feed can be controlled by probiotics and that mortality of infected fish larvae can be reduced significantly by probiotic bacteria. However, the successful management of the aquaculture microbiota is currently hampered by our lack of knowledge of relevant microbial interactions and the overall ecology of these systems.

Deep COI sequencing of standardized benthic samples unveils overlooked diversity of Jordanian coral reefs in the northern Red Sea

High-throughput sequencing (HTS) of DNA barcodes (metabarcoding), particularly when combined with standardized sampling protocols, is one of the most promising approaches for censusing overlooked cryptic invertebrate communities. We present biodiversity estimates based on sequencing of the cytochrome c oxidase subunit 1 (COI) gene for coral reefs of the Gulf of Aqaba, a semi-enclosed system in the northern Red Sea. Samples were obtained from standardized sampling devices (Autonomous Reef Monitoring Structures (ARMS)) deployed for 18 months. DNA barcoding of non-sessile specimens >2 mm revealed 83 OTUs in six phyla, of which only 25% matched a reference sequence in public databases. Metabarcoding of the 2 mm - 500 μm and sessile bulk fractions revealed 1197 OTUs in 15 animal phyla, of which only 4.9% matched reference barcodes. These results highlight the scarcity of COI data for cryptobenthic organisms of the Red Sea. Compared with data obtained using similar methods, our results suggest that Gulf of Aqaba reefs are less diverse than two Pacific coral reefs but much more diverse than an Atlantic oyster reef at a similar latitude. The standardized approaches used here show promise for establishing baseline data on biodiversity, monitoring the impacts of environmental change, and quantifying patterns of diversity at regional and global scales.

Phylogenetics and population genetics of Plotosus canius (Siluriformes: Plotosidae) from Malaysian coastal waters

Plotosus canius (Hamilton, 1822) is a significant marine species in Malaysia from nutritional and commercial perspectives. Despite numerous fundamental research on biological characteristics of P. canius, there are various concerns on the level of population differentiation, genomic structure, and the level of genetic variability among their populations due to deficiency of genetic-based studies. Deficiency on basic contexts such as stock identification, phylogenetic relationship and population genetic structure would negatively impact their sustainable conservation. Hence, this study was conducted to characterize the genetic structure of P. canius for the first time through the application of mitochondrial Cytochrome Oxidase I (COI) gene, cross amplification of Tandanus tandanus microsatellites, and a total of 117 collected specimens across five selected populations of Malaysia. The experimental results of the mitochondrial analysis revealed that the haplotype diversity and nucleotide diversity varied from 0.395–0.771 and 0.033–0.65 respectively. Moreover, the statistical analysis of microsatellites addressed a considerable heterozygote insufficiency in all populations, with average observed heterozygosity (H_o) value of 0.2168, which was lower than the standard heterozygosity in marine populations (H_o = 0.79). This alongside the high Fis values estimation, high pairwise differentiation among populations and low within population variations are supposed to be associated with small sample size, and inbreeding system. Besides, the significant finding of this study was the sharing of common haplotype KR086940, which reflects a historical genetic connectivity between Peninsular Malaysia and Borneo populations due to the geological history of Southeast Asia during Pleistocene era. Demographic analyses showed that all populations were in an equilibrium state with no significant evidence of population expansion. To put it briefly, the current study has managed to provide an initial genomic database toward understanding of the genetic characterization, phylogenetic, molecular diversification and population structure in P. canius, and should be necessary highlighted for appropriate management and conservation of species. Further studies must be carried out involving more geographical and sampling sites, larger population size per site, and utilization of species specific microsatellites loci.

Omics of the marine medaka (Oryzias melastigma) and its relevance to marine environmental research

In recent years, the marine medaka (Oryzias melastigma), also known as the Indian medaka or brackish medaka, has been recognized as a model fish species for ecotoxicology and environmental research in the Asian region. O. melastigma has several promising features for research, which include a short generation period (3-4 months), daily spawning, small size (3-4 cm), transparent embryos, sexual dimorphism, and ease of mass culture in the laboratory. There have been extensive transcriptome and genome studies on the marine medaka in the past decade. Such omics data can be useful in understanding the signal transduction pathways of small teleosts in response to environmental stressors. An omics-integrated approach in the study of the marine medaka is important for strengthening its role as a small fish model for marine environmental studies. In this review, we present current omics information about the marine medaka and discuss its potential applications in the study of various molecular pathways that can be targets of marine environmental stressors, such as chemical pollutants. We believe that this review will encourage the use of this small fish as a model species in marine environmental research.

Molecular screening of microbial communities for candidate indicators of multiple metal impacts in marine sediments from northern Australia

Coastal sediments accumulate metals from anthropogenic sources and as a consequence industry is required to monitor sediment health. The total concentration of a metal does not necessarily reflect its potential toxicity or biological impact, so biological assessment tools are useful for monitoring. Rapid biological assessment tools sensitive enough to detect relatively small increases in metal concentrations would provide early warning of future ecosystem impact. The authors investigated in situ populations of Archaea and Bacteria as potential tools for rapid biological assessment in sediment at 4 northern Australian coastal locations over 2 yr, in both wet and dry seasons. The 1 M HCl-extractable concentrations of metals in sediment were measured, and Archaeal and Bacterial community profiles were obtained by next-generation sequencing of sediment deoxyribonucleic acid (DNA). Species response curves were used to identify several taxonomic groups with potential as biological indicators of metal impact. Spatial variation, sediment grain size, water depth, and dissolved oxygen also correlated with microbial population shifts. Seasonal variation was less important than geographic location. Metal-challenge culture trials supported the identification of metal-resistant and -sensitive taxa. In situ Archaea and Bacteria are potentially sensitive indicators for changes in bioavailable concentrations of metals; however, the complexity of the system suggests it is important to identify metal-specific functional genes that may be informed by these sequencing surveys, and thus provide a useful addition to identity-based assays.

Marine microbial communities of the Great Barrier Reef lagoon are influenced by riverine floodwaters and seasonal weather events

The role of microorganisms in maintaining coral reef health is increasingly recognized. Riverine floodwater containing herbicides and excess nutrients from fertilizers compromises water quality in the inshore Great Barrier Reef (GBR), with unknown consequences for planktonic marine microbial communities and thus coral reefs. In this baseline study, inshore GBR microbial communities were monitored along a 124 km long transect between 2011 and 2013 using 16S rRNA gene amplicon sequencing. Members of the bacterial orders Rickettsiales (e.g., Pelagibacteraceae) and Synechococcales (e.g., Prochlorococcus), and of the archaeal class Marine Group II were prevalent in all samples, exhibiting a clear seasonal dynamics. Microbial communities near the Tully river mouth included a mixture of taxa from offshore marine sites and from the river system. The environmental parameters collected could be summarized into four groups, represented by salinity, rainfall, temperature and water quality, that drove the composition of microbial communities. During the wet season, lower salinity and a lower water quality index resulting from higher river discharge corresponded to increases in riverine taxa at sites near the river mouth. Particularly large, transient changes in microbial community structure were seen during the extreme wet season 2010–11, and may be partially attributed to the effects of wind and waves, which resuspend sediments and homogenize the water column in shallow near-shore regions. This work shows that anthropogenic floodwaters and other environmental parameters work in conjunction to drive the spatial distribution of microorganisms in the GBR lagoon, as well as their seasonal and daily dynamics.

The Next-Generation PCR-Based Quantification Method for Ambient Waters: Digital PCR

Real-time quantitative PCR (qPCR) is increasingly being used for ambient water monitoring, but development of digital polymerase chain reaction (digital PCR) has the potential to further advance the use of molecular techniques in such applications. Digital PCR refines qPCR by partitioning the sample into thousands to millions of miniature reactions that are examined individually for binary endpoint results, with DNA density calculated from the fraction of positives using Poisson statistics. This direct quantification removes the need for standard curves, eliminating the labor and materials associated with creating and running standards with each batch, and removing biases associated with standard variability and mismatching amplification efficiency between standards and samples. Confining reactions and binary endpoint measurements to small partitions also leads to other performance advantages, including reduced susceptibility to inhibition, increased repeatability and reproducibility, and increased capacity to measure multiple targets in one analysis. As such, digital PCR is well suited for ambient water monitoring applications and is particularly advantageous as molecular methods move toward autonomous field application.

An evaluation of coral lophelia pertusa mucus as an analytical matrix for environmental monitoring: A preliminary proteomic study

For the environmental monitoring of coral, mucus appears to be an appropriate biological matrix due to its array of functions in coral biology and the non-intrusive manner in which it can be collected. The aim of the present study was to evaluate the feasibility of using mucus of the stony coral Lophelia pertusa (L. pertusa) as an analytical matrix for discovery of biomarkers used for environmental monitoring. More specifically, to assess whether a mass-spectrometry-based proteomic approach can be applied to characterize the protein composition of coral mucus and changes related to petroleum discharges at the seafloor. Surface-enhanced laser desorption/ionization-time of flight mass spectrometry (SELDI-TOF MS) screening analyses of orange and white L. pertusa showed that the mucosal protein composition varies significantly with color phenotype, a pattern not reported prior to this study. Hence, to reduce variability from phenotype difference, L. pertusa white individuals only were selected to characterize in more detail the basal protein composition in mucus using liquid chromatography, mass spectrometry, mass spectrometry (LC-MS/MS). In total, 297 proteins were identified in L. pertusa mucus of unexposed coral individuals. Individuals exposed to drill cuttings in the range 2 to 12 mg/L showed modifications in coral mucus protein composition compared to unexposed corals. Although the results were somewhat inconsistent between individuals and require further validation in both the lab and the field, this study demonstrated preliminary encouraging results for discovery of protein markers in coral mucus that might provide more comprehensive insight into potential consequences attributed to anthropogenic stressors and may be used in future monitoring of coral health.