Michael Stat

Advances in high‐throughput sequencing (HTS) are revolutionizing monitoring in marine environments by enabling rapid, accurate and holistic detection of species within complex biological samples. Research institutions worldwide increasingly employ HTS methods for biodiversity assessments. However, variance in laboratory procedures, analytical workflows and bioinformatic pipelines impede the transferability and comparability of results across research groups. An international experiment was conducted to assess the consistency of metabarcoding results derived from identical samples and primer sets using varying laboratory procedures. Homogenized biofouling samples collected from four coastal locations (Australia, Canada, New Zealand and the USA) were distributed to 12 independent laboratories. Participants were asked to follow one of two HTS library preparation workflows. While DNA extraction, primers and bioinformatic analyses were purposefully standardized to allow comparison, many other technical variables were allowed to vary among laboratories (amplification protocols, type of instrument used, etc.). Despite substantial variation observed in raw results, the primary signal in the data was consistent, with the samples grouping strongly by geographical origin for all data sets. Simplepost hocdata clean‐up by removing low‐quality samples gave the best improvement in sample classification for nuclear 18S rRNA gene data, with an overall 92.81% correct group attribution. For mitochondrialCOIgene data, the best classification result (95.58%) was achieved after correction for contamination errors. The identified critical methodological factors that introduced the greatest variability (preservation buffer, sample defrosting, template concentration, DNA polymerase, PCR enhancer) should be of great assistance in standardizing future biodiversity studies using metabarcoding.

Aim: Environmental DNA (eDNA) metabarcoding has demonstrated its applicability as a highly sensitive biomonitoring tool across small spatial and temporal scales in marine ecosystems. However, it has rarely been tested across large spatial scales or biogeographical barriers. Here, we scale up marine eDNA metabarcoding, test its ability to detect a major marine biogeographic break and evaluate its use as a regional biomonitoring tool in Australia.Location: North‐western Australia (NWA).Methods: We applied metabarcoding assays targeting the mitochondrial 16S rRNA and CO1 genes to 284 surface seawater eDNA samples collected from 71 mid‐shelf, inshore, coastal and nearshore estuarine sites over 700 km of the NWA coastline.Results: Metabarcoding detected a wide range of bony fish (404 taxa), elasmobranchs (44) and aquatic reptiles (5). We detected bioregional and depth differentiation within inshore bony fish communities. These findings support the presence of a marine biogeographic break, which is purported to occur in the vicinity of Cape Leveque, demarcating the border between the Kimberley and Canning bioregions. Inshore bony fish and elasmobranch communities, as well as coastal bony fish assemblages, were additionally found to differ between the South and North Kimberley regions suggesting previously unrecognized subregional differentiation amongst these taxa. The overall compositional data have been used to update distribution information for a number of endangered, elusive and data‐deficient taxa, including sawfish (family: Pristidae), northern river shark (Glyphis garricki) and wedgefish (genus: Rhynchobatus).Main conclusions: eDNA metabarcoding demonstrated a high level of sensitivity that was able to discern fine‐scale patterns across the large‐scale, remote and oceanographically complex region of North‐western Australia. Importantly, this study highlights the potential of integrating broad‐scale eDNA metabarcoding alongside other baseline surveys and long‐term monitoring approaches, which are crucial for the sustainable management and conservation of marine biodiversity in this unique marine region.

Loss of biodiversity from lower to upper trophic levels reduces overall productivity and stability of coastal ecosystems in our oceans, but rarely are these changes documented across both time and space. The characterisation of environmental DNA (eDNA) from sediment and seawater using metabarcoding offers a powerful molecular lens to observe marine biota and provides a series of ‘snapshots’ across a broad spectrum of eukaryotic organisms. Using these next-generation tools and downstream analytical innovations including machine learning sequence assignment algorithms and co-occurrence network analyses, we examined how anthropogenic pressures may have impacted marine biodiversity on subtropical coral reefs in Okinawa, Japan. Based on 18 S ribosomal RNA, but not ITS2 sequence data due to inconsistent amplification for this marker, as well as proxies for anthropogenic disturbance, we show that eukaryotic richness at the family level significantly increases with medium and high levels of disturbance. This change in richness coincides with compositional changes, a decrease in connectedness among taxa, an increase in fragmentation of taxon co-occurrence networks, and a shift in indicator taxa. Taken together, these findings demonstrate the ability of eDNA to act as a barometer of disturbance and provide an exemplar of how biotic networks and coral reefs may be impacted by anthropogenic activities.

Environmental DNA (eDNA) metabarcoding, a technique for retrieving multispecies DNA from environmental samples, can detect a diverse array of marine species from filtered seawater samples. There is a growing potential to integrate eDNA alongside existing monitoring methods in order to establish or improve the assessment of species diversity. Remote island reefs are increasingly vulnerable to climate‐related threats and as such there is a pressing need for cost‐effective whole‐ecosystem surveying to baseline biodiversity, study assemblage changes and ultimately develop sustainable management plans. We investigated the utility of eDNA metabarcoding as a high‐resolution, multitrophic biomonitoring tool at the Cocos (Keeling) Islands, Australia (CKI)—a remote tropical coral reef atoll situated within the eastern Indian Ocean. Metabarcoding assays targeting the mitochondrial 16S rRNA and CO1 genes, as well as the 18S rRNA nuclear gene, were applied to 252 surface seawater samples collected from 42 sites within a 140 km² area. Our assays successfully detected a wide range of bony fish and elasmobranchs (244 taxa), crustaceans (88), molluscs (37) and echinoderms (7). Assemblage composition varied significantly between sites, reflecting habitat partitioning across the island ecosystem and demonstrating the localisation of eDNA signals, despite extensive tidal and oceanic movements. In addition, we document putative new occurrence records for 46 taxa and compare the efficiency of our eDNA approach to visual survey techniques at CKI. Our study demonstrates the utility of a multimarker metabarcoding approach in capturing multitrophic biodiversity across an entire coral reef atoll and sets an important baseline for ongoing monitoring and management.

Protein and carbohydrate have profound impacts on the ecology of gut microbiota, but disentangling the single and interactive effects of different dietary constituents is challenging. Here we used a multidimensional approach, the Geometric Framework, to study the interactions between nutrition and bacterial abundances with respect to protein and carbohydrate intakes in field cricket, Teleogryllus oceanicus. Our study revealed that species richness decreased as crickets consumed more macronutrients, and species evenness peaked at high intake of protein rich diets. Sex and protein:carbohydrate (P:C) ratios in diets were the primary factors influencing the gut bacterial community, but most of the microbial operational taxonomic units (OTUs) that were significantly different between males and females were present in low abundance. In contrast, protein intake had a greater influence than carbohydrate consumption on the relative abundances of the core bacterial taxa, as an increase in dietary protein availability could remove the growth constraint imposed by limited nitrogen. Taken together, the use of the Geometric Framework provides a deeper insight into how nutritional intakes influence the relative abundances of gut microbes, and could be a useful tool to integrate the study of gut microbiome and fitness traits in a host.

Marine ecosystems are changing rapidly as the oceans warm and become more acidic. The physical factors and the changes to ocean chemistry that they drive can all be measured with great precision. Changes in the biological composition of communities in different ocean regions are far more challenging to measure because most biological monitoring methods focus on a limited taxonomic or size range. Environmental DNA (eDNA) analysis has the potential to solve this problem in biological oceanography, as it is capable of identifying a huge phylogenetic range of organisms to species level. Here we develop and apply a novel multi-gene molecular toolkit to eDNA isolated from bulk plankton samples collected over a five-year period from a single site. This temporal scale and level of detail is unprecedented in eDNA studies. We identified consistent seasonal assemblages of zooplankton species, which demonstrates the ability of our toolkit to audit community composition. We were also able to detect clear departures from the regular seasonal patterns that occurred during an extreme marine heatwave. The integration of eDNA analyses with existing biotic and abiotic surveys delivers a powerful new long-term approach to monitoring the health of our world’s oceans in the context of a rapidly changing climate. All environments contain genetic remnants of the life they contain and support. For example, samples collected from the ocean contain biological material such as microscopic organisms, shed cells, excrement and saliva—the DNA from which reveals the surrounding marine biodiversity. Environmental DNA (eDNA) approaches have the ability to identify marine species that are notoriously difficult to identify using morphology alone. Here we develop, and apply, a DNA ‘toolkit’ to five years of samples collected from a single site in the Indian ocean. It is rare to find a temporal series of samples of this duration that are also suitable for DNA analysis. We show that eDNA techniques have the capacity to monitor ocean biology in fine detail. We demonstrate how the biological communities of plankton and fish respond to normal seasonal changes and, more importantly, to that of an uncharacteristic heatwave. The methods embodied in this paper are applicable to marine studies across the globe and, as such, pave the way for the design of long-term monitoring programs using eDNA.

Background: Effective biodiversity monitoring is fundamental in tracking changes in ecosystems as it relates to commercial, recreational, and conservation interests. Current approaches to survey coral reef ecosystems center on the use of indicator species and repeat surveying at specific sites. However, such approaches are often limited by the narrow snapshot of total marine biodiversity that they describe and are thus hindered in their ability to contribute to holistic ecosystem-based monitoring. In tandem, environmental DNA (eDNA) and next-generation sequencing metabarcoding methods provide a new opportunity to rapidly assess the presence of a broad spectrum of eukaryotic organisms within our oceans, ranging from microbes to macrofauna. Methods: We here investigate the potential for rapid universal metabarcoding surveys (RUMS) of eDNA in sediment samples to provide snapshots of eukaryotic subtropical biodiversity along a depth gradient at two coral reefs in Okinawa, Japan based on 18S rRNA. Results: Using 18S rRNA metabarcoding, we found that there were significant separations in eukaryotic community assemblages (at the family level) detected in sediments when compared across different depths ranging from 10 to 40 m (p = 0.001). Significant depth zonation was observed across operational taxonomic units assigned to the class Demospongiae (sponges), the most diverse class (contributing 81% of species) within the phylum Porifera; the oldest metazoan phylum on the planet. However, zonation was not observed across the class Anthozoa (i.e., anemones, stony corals, soft corals, and octocorals), suggesting that the former may serve as a better source of indicator species based on sampling over fine spatial scales and using this universal assay. Furthermore, despite their abundance on the examined coral reefs, we did not detect any octocoral DNA, which may be due to low cellular shedding rates, assay sensitivities, or primer biases. Discussion: Overall, our pilot study demonstrates the importance of exploring depth effects in eDNA and suggest that RUMS may be applied to provide a baseline of information on eukaryotic marine taxa at coastal sites of economic and conservation importance.

Diet studies provide base understanding of trophic structure and are a valuable initial step for many fields of marine ecology, including conservation and fisheries biology. Considerable complexity in marine trophic structure can exist due to the presence of highly mobile species with long life spans. Mobula rays are highly mobile, large, planktivorous elasmobranchs that are frequently caught either directly or as bycatch in fisheries, which, combined with their conservative life history strategy, makes their populations susceptible to decline in intensely fished regions. Effective management of these iconic and vulnerable species requires an understanding of the diets that sustain them, which can be difficult to determine using conventional sampling methods. We use three DNA metabarcode assays to identify 44 distinct taxa from the stomachs (n = 101) of four sympatric Mobula ray species (Mobula birostris, Mobula tarapacana, Mobula japanica, and Mobula thurstoni) caught over 3 years (2013–2015) in a direct fishery off Bohol in the Philippines. The diversity and incidence of bony fishes observed in ray diets were unprecedented. Nevertheless, rays showed dietary overlap, with krill (Euphausia) dominating their diet. Our results provide a more detailed assessment of sympatric ray diets than was previously described and reveal the complexity that can exist in food webs at critical foraging habitats.

While in recent years environmental DNA (eDNA) metabarcoding surveys have shown great promise as an alternative monitoring method, the integration into existing marine monitoring programs may be confounded by the dispersal of the eDNA signal. Currents and tidal influences could transport eDNA over great distances, inducing false‐positive species detection, leading to inaccurate biodiversity assessments and, ultimately, mismanagement of marine environments. In this study, we determined the ability of eDNA metabarcoding surveys to distinguish localized signals obtained from four marine habitats within a small spatial scale (<5 km) subject to significant tidal and along‐shore water flow. Our eDNA metabarcoding survey detected 86 genera, within 77 families and across 11 phyla using three established metabarcoding assays targeting fish (16S rRNA gene), crustacean (16S rRNA gene) and eukaryotic (cytochrome oxidase subunit 1) diversity. Ordination and cluster analyses for both taxonomic and OTU data sets show distinct eDNA signals between the sampled habitats, suggesting dispersal of eDNA among habitats was limited. Individual taxa with strong habitat preferences displayed localized eDNA signals in accordance with their respective habitat, whereas taxa known to be less habitat‐specific generated more ubiquitous signals. Our data add to evidence that eDNA metabarcoding surveys in marine environments detect a broad range of taxa that are spatially discrete. Our work also highlights that refinement of assay choice is essential to realize the full potential of eDNA metabarcoding surveys in marine biodiversity monitoring programs.