Genetic assessment of meiobenthic community composition and spatial distribution in coastal sediments along northern Gulf of Mexico

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.

Transition of an estuarine benthic meiofauna assemblage 1.7 and 2.8 years after a mining disaster

Background

Estuaries are transitional coastal ecosystems that are threatened by multiple sources of human pollution. In 2015, mining tailings from an upstream dam failure caused massive metal contamination that impacted benthic assemblages on the Brazilian Rio Doce estuary.

Methods

In this study, we investigate and compare meiofaunal assemblages with eDNA metabarcoding 1.7 years (2017) and 2.8 years (2018) after the initial contamination by mine tailings in order to evaluate the continued impact of sediment mine tailing contaminants on the structure of benthic assemblages after the disaster.

Results

The community was dominated by Arthropoda and Nematoda 1.7 yr after the impacts (42 and 29% of meiofaunal sequence reads, respectively) but after 2.8 years Arthropoda (64.8% of meiofaunal sequence reads) and Rotifera (11.8%) were the most common taxa. This continued impact on meiofaunal assemblage revealed a lower phylogenetic diversity (7.8-fold) in 2018, despite overall decrease in metal concentration (Al, Ba, Cr, As, Fe, Zn, Mn, Pb, Cd, Co) in sediments. Our data suggests that differences in benthic assemblages and loss of diversity may be influenced by contaminants in sediments of this estuary, and indicate that broad eDNA assessments are greatly useful to understand the full range of biodiversity changes in dynamic estuarine ecosystems.

Monitoring of benthic eukaryotic communities in two tropical coastal lagoons through eDNA metabarcoding: a spatial and temporal approximation

Tropical coastal lagoons are important ecosystems that support high levels of biodiversity and provide several goods and services. Monitoring of benthic biodiversity and detection of harmful or invasive species is crucial, particularly in relation to seasonal and spatial variation of environmental conditions. In this study, eDNA metabarcoding was used in two tropical coastal lagoons, Chacahua (CH) and Corralero (C) (Southern Mexican Pacific), to describe the benthic biodiversity and its spatial–temporal dynamics. The distribution of benthic diversity within the lagoons showed a very particular pattern evidencing a transition from freshwater to seawater. Although the two lagoon systems are similar in terms of the species composition of metazoans and microeukaryotes, our findings indicate that they are different in taxa richness and structure, resulting in regional partitioning of the diversity with salinity as the driving factor of community composition in CH. Harmful, invasive, non-indigenous species, bioindicators and species of commercial importance were detected, demonstrating the reach of this technique for biodiversity monitoring along with the continued efforts of building species reference libraries.

River Flow Impacts Bacterial and Archaeal Community Structure in Surface Sediments in the Northern Gulf of Mexico

Meiobenthic community structure in the northern Gulf of Mexico has been shown to be driven by geographical differences due to inshore-offshore gradients and location relative to river discharge. Samples collected along three transects spanning Mobile Bay, Alabama, showed significant differences in meiobenthic communities east of the bay compared to those sampled from the west. In contrast, analysis of bacterial and archaeal communities from the same sediment samples shows that the inshore-offshore gradient has minimal impact on their community structure. Significant differences in community structure were observed for Bacteria and Archaea between the east and west samples, but there was no difference in richness or diversity. Grouped by sediment type, higher richness was observed in silty samples compared to sandy samples. Significant differences were also observed among sediment types for community structure with bacteria communities in silty samples having more anaerobic sulfate reducers compared to aerobic heterotrophs, which had higher abundances in sandy sediments. This is likely due to increased organic matter in the silty sediments from the overlying river leading to low oxygen habitats. Most archaeal sequences represented poorly characterized high-level taxa, limiting interpretation of their distributions. Overlap between groups based on transect and sediment characteristics made determining which factor is more important in structuring bacterial and archaeal communities difficult. However, both factors are driven by discharge from the Mobile River. Although inshore-offshore gradients do not affect Bacteria or Archaea to the same extent as the meiobenthic communities, all three groups are strongly affected by sediment characteristics.

The use of metabarcoding for meiofauna ecological patterns assessment

Marine meiofauna comprises up to 22 phyla. Its morphological identification requires time and taxonomists' expertise, and molecular tools can make this task faster. We aim to disentangle meiofaunal diversity patterns at Araçá Bay by applying a model selection approach and estimating the effectiveness of metabarcoding (18S rDNA) and morphological methods for estimating the response of meiofauna diversity in small-scale interactions with environmental variables. A rarefaction curve indicated that ten samples were sufficient for estimating the total number of meiofauna OTUs in a tidal flat. In both approaches, richness was predicted by mean sand percentage, sediment sorting, and bacteria concentration. Nematode genera composition differed significantly between approaches, the result of taxonomic mismatch in the genetic database. The similarity between the model selected for diversity descriptors, the richness of nematode genera and meiofauna composition emphasized the utility of predictive models for metabarcoding estimates to detect small-scale interactions of these organisms.

Spatial and temporal variation of intertidal nematodes in the northern Gulf of Mexico after the Deepwater Horizon oil spill

Nematodes are an abundant and diverse interstitial component of sedimentary habitats that have been reported to serve as important bioindicators. Though the 2010 Deepwater Horizon (DWH) disaster occurred 60 km offshore in the Gulf of Mexico (GOM) at a depth of 1525 m, oil rose to the surface and washed ashore, subjecting large segments of coastline in the northern GOM to contamination. Previous metabarcoding work shows intertidal nematode communities were negatively affected by the oil spill. Here we examine the subsequent recovery of nematode community structure at five sites along the Alabama coast over a two-year period. The latter part of the study (July 2011-July 2012) also included an examination of nematode vertical distribution in intertidal sediments. Results showed nematode composition within this region was more influenced by sample locality than time and depth. The five sampling sites were characterized by distinct nematode assemblages that varied by sampling dates. Nematode diversity decreased four months after the oil spill but increased after one year, returning to previous levels at all sites except Bayfront Park (BP). There was no significant difference among nematode assemblages in reference to vertical distribution. Although the composition of nematode assemblages changed, the feeding guilds they represented were not significantly different even though some variation was noted. Data from morphological observations integrated with metabarcoding data indicated similar spatial variation in nematode distribution patterns, indicating the potential of using these faster approaches to examine overall disturbance impact trends within communities. Heterogeneity of microhabitats in the intertidal zone indicates that future sampling and fine-scale studies of nematodes are needed to examine such anthropogenic effects.

Spatio-temporal monitoring of deep-sea communities using metabarcoding of sediment DNA and RNA

We assessed spatio-temporal patterns of diversity in deep-sea sediment communities using metabarcoding. We chose a recently developed eukaryotic marker based on the v7 region of the 18S rRNA gene. Our study was performed in a submarine canyon and its adjacent slope in the Northwestern Mediterranean Sea, sampled along a depth gradient at two different seasons. We found a total of 5,569 molecular operational taxonomic units (MOTUs), dominated by Metazoa, Alveolata and Rhizaria. Among metazoans, Nematoda, Arthropoda and Annelida were the most diverse. We found a marked heterogeneity at all scales, with important differences between layers of sediment and significant changes in community composition with zone (canyon vs slope), depth, and season. We compared the information obtained from metabarcoding DNA and RNA and found more total MOTUs and more MOTUs per sample with DNA (ca. 20% and 40% increase, respectively). Both datasets showed overall similar spatial trends, but most groups had higher MOTU richness with the DNA template, while others, such as nematodes, were more diverse in the RNA dataset. We provide metabarcoding protocols and guidelines for biomonitoring of these key communities in order to generate information applicable to management efforts.