ORIGINAL RESEARCH article

Front. Mar. Sci., 02 February 2026

Sec. Ocean Observation

Volume 12 - 2025 | https://doi.org/10.3389/fmars.2025.1674917

Autonomous reef monitoring structures as a tool for assessing UK marine benthic biodiversity

  • 1. The Georgina Mace Centre for the Living Planet, Imperial College London, Ascot, United Kingdom

  • 2. Zoological Society London, Institute of Zoology, London, United Kingdom

  • 3. School of Geography, Queen Mary University of London, London, United Kingdom

  • 4. Institute of Agricultural and Environmental Sciences, Estonian University of Life Sciences, Tartu, Estonia

Abstract

In benthic marine ecosystems, small organisms that dwell within the habitat matrix comprise the majority of species richness but are inherently difficult to sample. This limits our ability to document the biodiversity of these ecosystems. Autonomous Reef Monitoring Structures (ARMS), artificial reefs that mimic the three-dimensional complexity of benthic habitats, can alleviate this challenge. However, ARMS have been applied infrequently in many temperate locations, including the United Kingdom (UK). To showcase the applicability of ARMS to the UK, this paper applies standardized Smithsonian ARMS protocols for image analysis as well as DNA barcoding and metabarcoding of the cytochrome c oxidase subunit I (COI) gene to 3 ARMS units deployed in the Sussex Inshore Fisheries and Conservation Authority (IFCA) Nearshore Trawling Exclusion zone, a conservation area created in 2021 to facilitate the recovery of seabed habitats. We document 176 genera across 22 eukaryotic phyla as well as 1,920 > 2 mm motile specimens with a biomass of 312 g. We identify 15 notable species, including the 4th UK record of Cephalothrix simula, a non-native poisonous nemertean with the potential to enter the food supply. We also assess the complementarity of image analysis and DNA metabarcoding in describing sessile communities, finding that the two methods produce meaningfully different estimates of relative abundance for some phyla, particularly Cnidaria (77x difference) and Bryozoa (72x difference). As our oceans come increasingly under threat from global change, it is imperative that we can accurately describe their biodiversity. We advocate for ARMS as a critical tool for measuring UK marine benthic biodiversity and discuss the use of ARMS in closing gaps in reference databases as well as in assessing ecosystem function and environmental disturbance in benthic habitats.

Introduction

Natural ecosystems are facing a range of increasing threats that are cumulatively contributing to global biodiversity loss, species range shifts, and ecosystem degradation (). In the UK, many ecosystems have suffered significant historical biodiversity losses and indicators of species abundance show continued declines, resulting in the UK being considered “one of the most nature-depleted countries on Earth” (). To arrest this decline in biodiversity, the UK government has enacted a variety of legal frameworks, such as the 2021 Environment Act, which established legally binding biodiversity targets (), and the Environmental Improvement Plan 2023, which prioritized “thriving plants and wildlife” as its apex goal ().

To assess the impact of management decisions aiming to enhance and protect biodiversity (e.g. the establishment of 377 Marine Protected Areas in the UK; Noble-James et al., 2023), and the impact of continued human development on biodiversity (e.g. the development of offshore wind farms; Watson et al., 2024), rigorous and comprehensive monitoring of biodiversity is imperative. The UK Government has set out to develop and implement a monitoring framework to achieve this, the Outcome Indicator Framework (OIF) (), but to date many key indicators remain in development and significant increases in biodiversity have yet to be achieved (Office for Environmental Protection, 2024).

The lack of biodiversity monitoring data suitable for evaluating long term trends is a gap that is particularly acute for marine benthic (seafloor) ecosystems. Life on the seafloor in the UK is highly diverse, with more than 10,000 species, however it is challenging to obtain data about these organisms due to where they live (). While the OIF provides long term tracking of the relative abundance of 1,177 species under “Wildlife” (Theme D), the only benthic marine species included are a small number of fish, none of which are given priority status (). Under “Seas and Estuaries” (Theme C) benthic habitats are assessed directly (indicator C4, “Diverse seas: condition of seafloor habitats”) but this indicator includes just 4 benthic habitat types, covers a narrow taxonomic range, and doesn’t provide species level data (). The lack of information on the state of the marine environment means that it is not possible to determine progress towards legal targets and the current infrequency of monitoring makes it difficult to provide timely course correction (Office for Environmental Protection, 2024).

This is a major concern because marine biota constitute billions of pounds of public value to the UK in the form of fisheries, carbon sequestration, recreation, and protecting against coastal erosion (Office for Natural Statistics, 2021). Benthic habitats act as refugia and nesting grounds for a broad array of taxa and so play a significant role in supporting overall marine biodiversity (). For example, over 1800 species have been documented in UK kelp forests, with a single kelp holdfast containing up to 70 macrofaunal species (Smale et al., 2013). Seagrass (Zostera sp.) beds, European oyster (Ostrea edulis) reefs, and rocky reefs have also been shown to exhibit high levels of macrofaunal biodiversity (Tu Do et al., 2012; ; ) as well as provide ecosystem services such as carbon sequestration, water purification, and wave and flood attenuation (; ). Therefore, losing these habitats also means losing much of the diverse “cryptobiome” of small organisms that comprises the majority of diversity in these three-dimensionally complex benthic habitats. These organisms include diverse invertebrate (e.g. sponges, ascidians, and crabs), algal, and microbial communities, which are major regenerators of organic (; ) and inorganic (; Rasheed et al., 2002) nutrients.

Many of these ecosystems have shown significant declines due to physical disturbance from fishing, sand extraction, and offshore construction (), as well as from sedimentation (), excess nitrogen, low light levels (), and over-harvesting (). Four canopy-forming kelp species have decreased in abundance in some parts of their range (Smale et al., 2013; Yesson et al., 2015), severe population declines have been shown for seagrass (since the 1869, peaking in the 1970’s; ), European oysters have declined across their biogeographic range (), and declines in rock-dependent taxa such as seaweeds imply increased pressure on rocky reefs (). Simultaneously, the UK marine environment has seen increased numbers of invasive species () and artificial structures such as wind farms (), changing the composition and habitat structure of marine benthic communities (Smale et al., 2013; Yesson et al., 2015). Restoration efforts are underway to conserve and restore many of these natural coastal ecosystems (e.g. the Sussex Kelp Recovery Project), but measuring the impact of these efforts remains an important task.

Here we evaluate the value that Autonomous Reef Monitoring Structures (ARMS) (Zimmerman and Martin, 2004) can provide in measuring biodiversity gains for such projects and, more generally, for monitoring benthic biodiversity in the UK. ARMS consist of a stack of PVC plates that mimic the three-dimensional complexity of benthic habitats such as coral, oyster and chalk reefs, kelp forests, and human-made structures. After being deployed and colonized by the diverse and hard to sample “cryptobiome” found in these habitats, ARMS are retrieved with little disturbance to natural substrates and analyzed using standardized protocols for image and molecular analysis (Leray and Knowlton, 2015; Ransome et al., 2017). This provides a comprehensive insight into diversity that goes beyond conspicuous macrofauna. Originally designed for use on tropical reefs, very few studies using ARMS have been carried out in temperate regions (; ; Pagnier et al., 2025), and ARMS have yet to be used to assess UK benthic biodiversity.

Materials and methods

ARMS deployment & collection

ARMS (n = 3) were deployed on July 15, 2022 at a site off the coast of Sussex (50.784452, -0.372378; Figure 1A) in an area targeted for seabed recovery after the loss of an extensive kelp habitat (Williams and Davies, 2019). Depth at deployment was 7m (at chart datum).To promote recovery of essential fish and kelp habitat in this area, a nearshore trawling byelaw was implemented in 2021 by the Sussex Inshore Fisheries and Conservation Authority (IFCA) to protect the seabed from the impact of demersal trawls and enable habitat restoration (Sussex Inshore Fisheries and Conservation Authority (IFCA), 2019). The seabed at the site was consistent with EUNIS habitat classification 5.5 (sublittoral macrophyte-dominated sediment). After 15 months, on October 30, 2023, ARMS were collected and processed according to Smithsonian Global ARMS protocols 2–7 and 9 (The Smithsonian Institution, 2024). See Figure 1B for an image of one of the ARMS post-colonization.

Figure 1

ARMS plate image annotation

Coverage of the sessile biota on ARMS plate faces (e.g. Figure 1B) was assessed according to Smithsonian Global ARMS protocol 15 (The Smithsonian Institution, 2024). Images were analyzed in CoralNet () using 225 annotation points per image, generated using a 15 x 15 cell grid, with one point randomly placed within each cell. Each annotation point was assigned 1 of 32 labels following the National Oceanic and Atmospheric Administration’s “Coral Reef Ecosystem Division/Program for standardizing ARMS analysis” protocol (), with some additional groupings (see Supplementary Table 1), with all identifications performed manually without CoralNet’s machine learning engine. Labels covered 21 sessile functional groups, 1 motile fauna group, 7 substrate types, unknown taxa, no recruitment or unavailable recruitment space (see Supplementary Table 1 for full descriptions).

Motile specimen morphosorting and DNA barcoding

Large (> 2 mm) motile specimens were sorted into morphotypes based on physical characteristics (e.g., number and shape of appendages) cross referencing vs . Individual specimens were weighed to within 0.01 g and up to 5 specimens from each morphotype were selected as type specimens. Some arthropods and annelids were fragmented during morphotyping, so counts are approximate for these taxa. Photographs of specimens were taken using either a Celestron microscope or a Nikon D810 Digital SLR camera with a 105 mm lens. All specimens were preserved in 95% ethanol at -20°C prior to and following morphosorting.

DNA was extracted from 10–20 mg of tissue (if available) from individual specimens covering all morphotypes (n = 46; see Supplementary Table 2 for photos and details) using a DNEasy Blood and Tissue Kit (Qiagen) following manufacturer’s instructions for the spin column protocol. Multiple morphologically distinct specimens were sequenced for more abundant morphotypes. DNA extracts were amplified through Polymerase Chain Reaction (PCR) targeting a 658 base pair region of the cytochrome c oxidase subunit I (COI) gene (primers dgLCO1490/dgHCO2198; Meyer, 2003). Because these protocols were applied to novel taxa, multiple rounds of PCR were run to establish optimal parameters for amplification for each specimen (see Supplementary Table 2). Results were imaged on a 1% agarose gel. All samples were run in 50 μL volumes, which included 10 μL HOTFIREPol Taq mix (Solis BioDyne), 2 μL forward primer at 10 mM, 2 μL reverse primer at 10 mM, 35 μL PCR-grade H2O, and 1 μL extracted DNA. A negative control was included in all runs. The PCR consisted of an initial denaturation at 95°C for 15 minutes, 37 cycles of 95°C denaturation for 40 seconds, annealing for 40 seconds at various temperatures, and 72°C extension for 60 seconds, followed by a final 72°C extension for 10 minutes. Amplicons were quantified using a Qubit BR Assay Kit (Thermo Fisher Scientific) following manufacturer’s instructions and sent to Eurofins Genomics for purification and Sanger sequencing of the forward read.

DNA metabarcoding of bulk fractions

The sessile, 106 µm-500 µm motile, and 500 µm-2 mm motile bulk fraction samples (n = 9) were prepared for DNA extraction according to Smithsonian ARMS protocol 10 (The Smithsonian Institution, 2024) with some alterations (see Supplementary Materials for details). Adhering to Smithsonian ARMS protocols, DMSO preserved samples were used for sessile fraction metabarcoding and ethanol preserved samples were used for motile fraction metabarcoding. DMSO better preserves DNA in the homogenized sessile fraction (Ransome et al., 2017), whereas it does not preserve the morphological features of some invertebrates (Ogiso-Tanaka et al., 2025), so is avoided for unhomogenised motile specimens.

DNA was extracted and cleaned using a DNEasy PowerMax Soil Kit (Qiagen) and a DNEasy PowerClean Pro Cleanup Kit (Qiagen), according to Smithsonian ARMS protocol 11 (The Smithsonian Institution, 2024) except where the protocol specified it followed the manufacturer’s instructions. DNA concentrations of cleaned samples were measured using a Qubit HS Assay Kit (Thermo Fisher Scientific) following manufacturer’s instructions. A 313 base pair region of the COI gene (Leray et al., 2013; The Smithsonian Institution, 2024) was amplified via PCR. PCRs were run in triplicates in 20 μL volumes and included 4 μL HOT FIREPol Taq Mix (Solis BioDyne), 1 μL of forward primer mlCOIint_F (including Illumina forward overhang adapter sequence) at 10 mM, 1 μL of reverse primer jgHCO2198 (including Illumina reverse overhang adapter sequence) at 10 mM, and 2 ng of DNA for sessile samples or 5 ng of DNA for motile samples. The PCR profile from Smithsonian protocol 11 was used with the initial denaturation step extended to 15 minutes and the final extension step extended to 10 minutes. A negative control was run for each PCR reaction and results were imaged on a 1% agarose gel.

Library preparation consisted of a custom version of the Illumina 16S Metagenomic Sequencing Library Preparation protocol. PCR triplicates were pooled and cleaned using AMPure XP beads (Beckman Coulter). An indexing PCR was run in 50 μL volumes and included 10 μL HOT FIREPol Taq Mix (Solis BioDyne), 5 μL of Index primer 1, 5 μL of Index primer 2, 25 μL PCR-grade H2O, and 5 μL cleaned first round PCR product. Indexes were chosen from a Nextera XT DNA Library Preparation Kit (Illumina) to provide a unique combination for each sample. The PCR profile used was as follows: an initial denaturation at 95°C for 15 minutes, then 8 cycles of 95°C denaturation for 30 seconds, 55°C annealing for 60 seconds, and 72°C extension for 60 seconds, followed by a final 72°C extension for 10 minutes. Indexed product was cleaned using AMPure XP beads (Beckman Coulter) and DNA concentration measured using a Qubit HS Assay Kit (Thermo Fisher Scientific) following manufacturer’s instructions. Samples were then normalized and pooled. Libraries were denatured, diluted, spiked with 10% PhiX, loaded onto a MiSeq Reagent Nano Kit v2 (500-cycles) (Illumina) and sequenced on a MiSeq platform.

Data analysis

Unless otherwise noted, analysis was performed using R (v4.4.2) in R Studio (v2024.04.2 + 764) and applying a common set of packages for data manipulation, summary, and plotting. Package details and scripts are available at https://github.com/MGJordan/SKRP/.

Image annotations were extracted from CoralNet and the percentage cover of individual sessile phyla on available recruitment space was calculated. Organisms were visually identified to species level where possible (Marine Life Information Network (MarLIN), 2016; ). Community composition across ARMS was visualised via non-metric multidimensional scaling (NMDS) using the vegan R package (Oksanen et al., 2024) with Bray-Curtis distances and a fourth-root transformation of count data to prevent highly abundant species from dominating the results. Using the negative binomial distribution and a log link function, the manyglm function from the mvabund R package (v.4.2.1) (Wang et al., 2012) was used to run a multivariate generalized linear model (GLM) on the effect of plate face (top or bottom ARMS plate surfaces, see Figure 1B) on sessile community composition. Univariate GLMs and subsequent analyses of variance (ANOVA) were then run to assess the effect of plate face on the abundance of individual sessile functional groupings with at least 30 annotation points.

Barcoding sequences were trimmed on both the 5’ and 3’ ends using the Trim Ends function in Geneious Prime (v2024.07) with an Error Probability Limit of 0.05. The dgLCO1490 primer was trimmed allowing up to 5 mismatches. Trimmed sequences were queried using the Barcode of Life Datasystem (BOLDv5; using BOLDigger3 command line application (–db 3 –mode 3) with default similarity thresholds (species: 97%, genus: 95%, family: 90%, order: 85%) and the Genbank nucleotide collection (Geneious Prime (v2024.0.7); using Megablast search and default parameters). A taxonomic assignment from each database was produced for each sequence. For BOLD, the BOLDigger3 algorithm assigned taxonomy following . For Genbank, up to 100 hits were returned and the taxonomy with the greatest number of hits that were > 85% query coverage and > 97% pairwise identity was used (Ransome et al., 2017). Genbank results were checked for synonymity using the World Registry of Marine Species (WoRMS) and resolved to the accepted taxonomy (World Registry of Marine Species (WoRMS) Editorial Board, 2025). Sequences for which taxonomic assignments did not match specimen morphology (e.g., an arthropod specimen assignment to a sponge taxonomy) were left unassigned. Where BOLD and Genbank produced different assignments (n = 1), the assignment from the database with more hits above the threshold was used.

Metabarcoding sequences were processed using the SimpleMetaPipeline (v0.0.0.9000) command line application running in a miniconda (v24.5.0) environment (Williams et al., 2024). Briefly, this pipeline uses DADA2 (v1.26.0) to perform primer trimming, quality filtering, denoising, merging, and chimera removal, VSEARCH (v2.14.2) to cluster amplicon sequence variants (ASVs) into operational taxonomic units (OTUs), and LULU (v0.1.0) to curate OTUs by identifying erroneous OTUs and merging them with their parent (; Rognes et al., 2016; ). The IDTAXA algorithm (Murali et al., 2018) was used to classify sequences against a custom classifier designed for marine COI metabarcoding by curating MIDORI2 (v248) (Leray et al., 2022) and BIOCODE (Ransome et al., 2017) databases. The confidence of assignment threshold was 60%, corresponding to very high confidence (Murali et al., 2018), which is appropriate for identifying taxa at lower taxonomic levels. The pipeline also queried sequences against the same custom database using BLAST () and assignments above specific similarity threshold were retained depending on the taxonomic rank (species: 97%, genus: 95%, family: 90%, order: 85%, class: 80%, phylum: 0%). IDATAXA and BLAST results were imported into R Studio using the SimpleMetaPackage (v0.0.0.9000) R package (Williams et al., 2024) and curated OTUs for clustering. BLAST results were filtered to include only assignments with >= 85% query coverage and >= 97% percent identity. Because IDTAXA has been shown to produce lower rates of over-classification (Murali et al., 2018), final assignments were first granted to IDTAXA and then granted to BLAST if no IDTAXA assignment was made. Final assignments were cleaned and (to facilitate comparison between datasets) taxonomies missing a phylum were assigned a phylum by searching the next lowest taxon versus the Catalogue of Life (). The microViz (v0.12.5) R package was used to remove non-eukaryotic OTUs () to ensure non-target OTUs were removed from our dataset. Multiple rarefactions were performed (1,000 iterations) to the lowest read count among samples (30,295) (Mikryukov, 2023). Rarefied reads were calculated by taking the mean reads per OTU across iterations.

Negative controls were checked for OTUs. The sessile negative control contained no OTUs. The motile negative control contained 10 OTUs of which 7 were eukaryotic (Amphinema dinema, Balanidae sp., Crepidula fornicata, Demospongiae sp., Phyllodocidae sp., Polynoidae sp., Sabellaria spinulosa). These 7 OTUs were also present in all sessile samples, suggesting that their presence in the motile negative control was due to cross-contamination between samples. Because this does not impact comparisons reported below, these OTUs were not removed to preserve the fullest estimate of biodiversity at the site.

Species of special interest

The presence of species of conservation value, of commercial value, or previously undocumented in the UK was evaluated across image, barcoding, and metabarcoding datasets using WoRMS (World Registry of Marine Species (WoRMS) Editorial Board, 2025), MarLIN (Marine Life Information Network (MarLIN), 2016), the Global Biodiversity Information Facility (), and the National Biodiversity Network (NBN) Atlas (National Biodiversity Network (NBN) Trust, 2024). The presence of non-native and invasive species was evaluated using the World Register of Introduced Marine Species () versus records downloaded for the Atlantic Ocean (General Sea Area) including subareas and distributions from synonyms. The presence of threatened seaweed species was also evaluated based on the British Red List of seaweeds () and for 7 kelps with broad UK distribution (Smale et al., 2013). Due to uncertainty in potential species assignment, species of special interest with UK congenerics and with missing representation from reference libraries were not reported as special interest but remained in the results.

Comparison between image analysis and metabarcoding data

Phyla-level community composition recovered by image analysis and metabarcoding approaches was compared, using only phyla identified in image analysis. For image analysis, the relative abundance of each phylum was calculated as the number of annotation points assigned to that phylum out of the total number of points assigned to living sessile organisms. Any points assigned to empty or broken shells or bare substrate were not included in this calculation. For metabarcoding, the relative abundance of each phylum was calculated as the number of rarefied eukaryotic reads assigned to a phylum out of the total number of rarefied eukaryotic reads among phyla only identified in image analysis.

Results

A total of 176 genera, from 22 phyla were identified across three ARMS

A total of 176 genera of organisms, from 153 families (across 42 classes, 81 orders, and 22 phyla), were recovered from ~0.15m2 of planar benthic surface area (ARMS n = 3) at a single site off the Sussex coast using barcoding, metabarcoding and image analysis. The majority of the species retrieved were animals (82%) with some algae (15%) and a small proportion of other groups, including fungi (2%) and non-algal protists (2%) (see Supplementary Table 3 for a full list of identified taxa). The three phyla showing the greatest richness at Family level were the Arthropoda (35 families), Annelida (23 families), and Mollusca (21 families).

Image analysis alone identified 10 phyla across 18 sessile functional groups

Image analysis of ARMS plates identified 10 phyla across 13 sessile invertebrate and 5 algal functional groupings (see Supplementary Table 1 for details). A total of 16 species were identified (see Supplementary Table 3) including three species of interest: the invasive slipper limpet Crepidula fornicata (), the habitat forming honeycomb worm Sabellaria alveolata and sand mason worm Lanice conchilega (; ).

Among the 10,347 points annotated, 57.9% were assigned to sessile organisms (53.9% invertebrates, 4.1% algae), 21.2% to substrate (e.g. sediment), 15.8% to no recruitment, 4.5% to unknown, and 0.4% to motile organisms. Among all annotation points assigned to live sessile organisms, a majority (93%) were assigned to sessile invertebrates, in particular bryozoans (43%), annelids (33%), and porifera (10%) (Figure 2). ARMS communities separated into two distinct clusters in an NMDS, based on plate surface (top or bottom) (see Supplementary Figure 1). Plate face had a significant effect on the abundance of 7 out of 10 sessile invertebrate functional groupings (ANOVA, p<0.05; see Supplementary Table 4), with crustose coralline algae and serpulid worm sediment found in higher abundances on top plate surfaces, and bryozoans, calcareous worms, and sponges found in higher abundances on bottom plate surfaces.

Figure 2

Recovered large motile (> 2mm) fauna consisted of 1920 specimens and 312 grams of biota mass

One thousand nine hundred and twenty motile organisms, sorted into 35 morphospecies, were collected from 3 ARMS. Of the 46 specimens picked for barcoding, 28 specimens (19 morphotypes) were assigned a taxonomy, 27 to species-level. BOLD and Genbank taxonomic assignments were identical for all but 1 specimen. Taxonomic identification confirmed unique identities of each morphotype (Supplementary Table 2). These 19 morphotypes accounted for 91% of all large motile individuals collected and 75% of the total motile biomass (312 grams) (see Figure 3, Supplementary Table 6 for morphotype counts and weights). The remaining 18 specimens (15 morphotypes) were not assigned a taxonomy for a variety of reasons including lack of matches to taxonomic records despite average sequence quality, poor sequence quality, and non-specific PCR amplification (see Supplementary Table 2 for full details).

Figure 3

The > 2 mm motile assemblage was dominated by arthropods (63.5% of specimens, 37.1% of biomass) and molluscs (26% of specimens, 36.5% of biomass) (Figure 3), including Pisidia longicornis (long-clawed porcelain crab), Tritia incrassata (thick-lipped dogwhelk), Athanas nitescens (hooded shrimp), and Pilumnus hirtellus (bristly crab) (Figure 4), which collectively comprised 80% of all specimens and 57% of the total biomass.

Figure 4

Metabarcoding of sessile and small motile (<2mm) communities recovered 215 eukaryotic OTUs

A total of 1,051,408 reads were generated from metabarcoding 9 bulk ARMS samples, of which 889,457 passed quality control, denoising, merging, and chimera-removal steps. The mean cleaned reads per sample was 98,823. Reads were clustered into 1,641 OTUs of which 87% (41% of reads) were removed from further analysis as they were not identified as eukaryotic, leaving 215 eukaryotic OTUs comprising 525,029 reads for downstream analysis.

In terms of read abundance, the ARMS assemblage was dominated by annelids (61% of reads) including Pholoe inornata (scale worm, 23% of reads) and Sabellaria spinulosa (ross worm, 21% of reads) (Figure 5). Arthropods (15% of reads), sponges (10% of reads), cnidarians (5% of reads), and molluscs (4% of reads) were also abundant. The assemblage was much more evenly distributed when considering sequence diversity (OTU richness; Figure 5), with annelids as the most diverse phyla (21% of OTUs), followed by arthropods (19% of OTUs), molluscs (11% of OTUs), sponges (9% of OTUs), and cnidarians (7% of OTUs). Multiple OTUs were also recovered for the rhodophytes (8% of OTUs), bryozoa (5% of OTUs), chordates (5% of OTUs), and ochrophytes (5% of OTUs), despite relatively low read abundances for these phyla (Figure 5).

Figure 5

Taxonomic assignment at phylum and species-level was 99.6% and 86.4% of reads assigned respectively. Taxonomic assignment of OTUs was also high at the phylum level (99.5%), but considerably lower at the species-level (70.7%) (Figure 5).

Fifteen species of interest were identified, but no kelp was detected

Species of interest identified included 4 habitat forming species, 3 commercially fished species, 1 established non-native species, 2 recently introduced non-native species, 4 British Red List seaweeds, and 1 species with no prior recorded distribution in the UK (Table 1). Of these, 4 were identified by image analysis, none were identified by barcoding, and 11 were identified by metabarcoding, with no overlap between methods. No kelp were identified in metabarcoding libraries, although all 7 UK kelp species had representative sequences in our reference databases.

Table 1

SpeciesCommon nameStatus notesMethod
Lanice conchilegaSand Mason worm• Habitat forming
• OIF (C4)
IA
Ostrea edulisEuropean oyster• Habitat forming
• OSPAR
• OIF (C6)
MB
Sabellaria alveolataHoneycomb worm• Habitat forming
• OIF (C4)
IA
Sabellaria spinulosaRoss worm• Habitat forming
• OSPAR
• OIF (C6)
MB
Aequipecten opercularisQueen scallop• Commercially fished in UKIA
Raja clavataThornback skate• Commercially fished in UK
• IRL Near-Threatened
MB
Symphodus melopsCorkwing wrasse• Commercially fished in UKMB
Crepidula fornicataSlipper limpet• Long established non-native
• Invasive
MB
Cephalothrix simula–• Recently introduced non-native
• Invasive
MB
Perophora japonica–• Recently introduced non-nativeMB
Halidrys siliquosaSea oak• NRL Near-ThreatenedMB
Polyides rotundaDiscoid Forked Weed• NRL Near-ThreatenedIA
Porphyra umbilicalisLaver• NRL Critically EndangeredMB
Rhodymenia delicatula–• NRL VulnerableMB
Igernella notabilis–• Potential first UK recordMB

Notable species identified with notes on status and assignment method [metabarcoding (MB) or image analysis (IA)].

Status notes include the species’ conservation status (National Red List for British Seaweeds (NRL), ; IUCN Red List of Threatened Species (IRL), ) and whether the species is habitat forming (Oslo and Paris (Oslo and Paris Conventions (OSPAR) Commission, 2024a, 2024), OSPAR listed (Oslo and Paris Conventions (OSPAR) Commission, 2024a; 2024b), included in an OIF indicator (with which indicator) (), commercially fished in the UK (; ; ), not native (Martin et al., 2006; Nishikawa et al., 2000), invasive (), or a potential first UK record (World Registry of Marine Species (WoRMS) Editorial Board, 2025).

Taxonomic variability between image analysis and metabarcoding

Of the 10 phyla identified using image analysis, metabarcoding identified 9, omitting Foraminifera. A two-fold difference in relative abundance was observed between the two methods for 6 phyla (Figure 6). The largest differences were observed for Cnidaria (77x higher relative abundance in metabarcoding than in image analysis) and Bryozoa (72x higher relative abundance in image analysis than in metabarcoding) (Figure 6).

Figure 6

Discussion

Here we recovered a diverse benthic community in an important UK seabed recovery zone, highlighting the utility of ARMS for seabed biomonitoring. In total, 176 genera were identified across 22 eukaryotic phyla using image analysis, barcoding, and metabarcoding approaches, which included 15 species of interest. We observed that after one-and-a-half years, sessile organisms covered about half (57.9%) of available ARMS plate surfaces, providing habitat to 1,920 motile specimens with a biomass of 312 g across just 0.15 m2 of planar space. Additionally, metabarcoding detected the presence of 215 eukaryotic OTUs.

Using ARMS provides a more detailed visual analysis of sessile benthic taxonomic richness compared to previous efforts that are biased toward conspicuous taxa such as diver surveys (Seasearch, 2025), towed video cameras (Mallinson and Yesson, 2021), and baited remote underwater video cameras (). Because these methods rely on visual inspection of the habitat surface through the water column, they are prone to under-detection of the small or cryptic organisms (; Willis, 2001; ; Lowry et al., 2012) that comprise the majority of benthic biodiversity (Steyaert et al., 2022a). In contrast, image analysis of ARMS plate faces facilitates identification of small and cryptic organisms by exposing the habitat matrix interior and obviating visibility bias caused by factors such as water turbidity (), camera orientation (), and habitat complexity (). This is evident when comparing our results to Mallinson and Yesson (2021) towed video camera census in Sussex, which identified sessile animalia from 4 phyla (Annelida, Bryozoa, Cnidaria, Porifera); ARMS image analysis identified sessile animalia from 8 phyla (the same plus Arthropoda, Mollusca, Chordata, and Foraminifera), and included small organisms such as barnacles and tunicates that are unlikely to be evident in towed video images.

ARMS also provide the opportunity to collect benthic samples for molecular analysis in a non-destructive and standardized way, while still recovering diverse epibenthic communities. Previous UK environmental DNA biodiversity monitoring studies have typically sampled DNA from the demersal layer (water above the benthos; ) or sediment cores (Steyaert et al., 2020), and protocols have recently been developed for eDNA benthic sampling from boat grabs (e.g. Van Veen grabs or box cores) for the statutory monitoring of benthic species and habitats of conservation importance (NatureMetrics, 2022). However, while demersal samples are far easier to obtain, they fail to detect species that are embedded in the sediment (e.g. crustaceans and molluscs; Le Joncour et al., 2024), and in reef crevices (Nichols et al., 2022). They have also been shown to be a poor proxy for the analysis of benthic structure compared to direct sampling methods (), potentially underrepresenting groups such as polychaetes, shrimps, mollusks, brittle stars, and red algae (Nichols et al., 2022). In contrast, while eDNA sampling from sediment cores and grabs likely provide more pertinent information about benthic community structure, these methods are destructive to natural ecosystems and are likely difficult to sample in a standardized way when dominated by mixed or rocky substrate due to the complex structure of hard habitats (Rees, 2009). ARMS plug this gap for UK biodiversity monitoring by providing a standardized sample area and access to communities found in three-dimensional habitats, such as those found living on rocky reefs, seagrass and kelp.

In addition to providing more complete biodiversity estimates, these results demonstrate the ability of ARMS to non-destructively recover both species of conservation value and invasive species. Included in this list is Ostrea edulis, Rhodymenia delicatula, Crepidula fornicata, and Cephalothrix simula. O. edulis has seen massive declines with urgent need for active management (), making in situ observations valuable. Nationally vulnerable red seaweed R. delicatula is predominantly reported from SW England, Wales and West Scotland, with no recent observations in the eastern Channel (), thus expanding our knowledge of its current distribution is important. The observed invasive C. fornicata has been shown to alter benthic community composition and trophic structure (Martin et al., 2006; ), and the invasive C. simula is notable because its tissue contains Tetrodotoxin (TTX), a potent neurotoxin (Turner et al., 2018), with an individual worm potentially carrying lethal doses of TTX for humans (). Although the sample first reported in the UK was below the lethal dose, it is nonetheless important to understand the distribution of C. simula in the UK due to its co-occurrence with shellfish and potential to enter the food supply (Turner et al., 2018). With only two other confirmed records in the UK, both in Southern England (; ), our study presents a valuable data point about its distribution.

No kelp taxa (families Laminariaceae, Phyllariaceae) were identified in our study. Notably, the pseudo-kelp Saccorhiza polyschides is commonly observed near the deployment site (Ray Ward pers. com.) but was not recovered. However, testing the metabarcoding primers versus publicly accessible COI accessions for this species shows a 5 base pair mismatch for this species in the forward primer. This highlights the value of understanding primer applicability for local species of interest. To date this research has focused on particular taxa such as fish (e.g. Schenekar et al., 2020), rather than habitat forming benthic species.

Metabarcoding studies also depend on quality reference sequences for the identification of organisms (Weigand et al., 2019). However, there are currently large taxonomic gaps in marine barcoding databases, particularly for marine invertebrate COI sequences in the EU (Weigand et al., 2019), hampering taxonomic assignment, particularly to the species level (Wort et al., 2022). We show that species level identification rates vary greatly for UK marine benthic habitats. For example, 15% of Arthropoda OTUs and 47% of Porifera OTUs were left unidentified at the species level. Some less abundant phyla such as Nematoda had no OTUs assigned at the species level. This lack of species level information makes confirmation of range extensions difficult and makes it challenging to understand the impact of conservation measures on marine biodiversity. While direct sampling of biomass is best suited for building taxonomies (Nichols et al., 2022), it is challenging to obtain these samples from benthic grabs, where the community has been mixed and crushed. ARMS provide the opportunity to collect clean samples for barcoding, which would greatly improve our taxonomic coverage of UK marine benthic taxa in reference databases.

The combination of image and molecular analysis of ARMS communities can provide important insights for benthic biomonitoring. Few studies to date have analysed images from ARMS plates (Sembiring et al., 2023; Steyaert et al., 2022a, 2022; ), and even fewer in combination with metabarcoding (; ). Here, integrating ARMS image data resulted in the identification of four species of special interest that are not present in our metabarcoding libraries and provided a less costly alternative to metabarcoding. However, the resulting biodiversity data from both methods need to be interpreted with care. While our study found that metabarcoding and image analysis produced highly similar estimates for sponge abundance, for other phyla, different results were obtained across methods. Metabarcoding detected 9 out of the 10 phyla identified by image analysis. Foraminifera, the only phylum not identified by metabarcoding, requires taxon specific primers (Macher et al., 2022). Metabarcoding retrieved higher abundances of motile phyla (e.g. Arthropoda) and Cnidaria (e.g. hydroids), both of which are difficult to pick up from images. In contrast, image analysis retrieved significantly higher abundances of Bryozoa, which are one of the most abundant and diverse groups of invertebrates in coastal marine ecosystems, but are underrepresented in genetic databases due to the difficulty of sampling (). These results highlight strength in using both methods. However, the accumulating global dataset of ARMS images is a vastly untapped resource, and with the use of machine-learning approaches such as CoralNet’s neural networks, may provide expedited and important outputs when analysed alongside key environmental and anthropogenic variables.

A further untapped use of ARMS is their ability to provide standardized biomass-size spectrum data of benthic communities. Here we show that while morphotypes identified in this study were similar to the number retrieved in other temperate studies (e.g. in Virginia, USA 38 morphotypes were retrieved from 9 ARMS after a 6 month deployment; Leray and Knowlton, 2015), we retrieved significantly more large (>2mm) specimens (n = 1,920) than other studies. The number of large motile specimens identified here from 3 ARMS is almost 5x greater than the 498 specimens across 9 ARMS in Virgina. It is also ~3x greater than the 725 specimens identified from 12 ARMS deployed for 24 months in Singapore (). In size-based ecosystem modelling, the biomass-size spectrum is the most widely used method to express the size structure of a community () and refers to the distribution of living biomass across organism size range (Tung et al., 2023). It has been applied in many settings, including in benthic intertidal communities, showing that the physical environment affects communities by creating habitats for organisms of different size classes (Schwinghamer et al., 1986). Communities of different size spectrums may carve distinct niches by influencing biological processes (e.g. metabolism and carbon flux; ) and deviations in size spectrum may be used to identify environmental impacts (Tung et al., 2023). For example, as larger organisms have higher energetic requirements, slower growth, and lower reproductive output (), they are more vulnerable to environmental disturbances (). Given the history of bottom trawling as well as the concerns about pollution in Sussex waters (Williams and Davies, 2019), the large abundance of small organisms found here may be a direct response to the degraded environment or a result of species flux during recovery at our sites, as has been noted in other coastal ecosystems (). We suggest that biomass, size and count data is useful to collect and share across ARMS studies and could provide data about environmental disturbance and ecosystem function for UK biomonitoring.

Conclusion

Although ARMS were originally designed to mimic the matrix of coral reefs (Zimmerman and Martin, 2004), deployments in polar and temperate regions have demonstrated that they have broader applicability to benthic biodiversity monitoring (). Our results highlight the variety of applications of the multidimensional data collected under standardized ARMS protocols. This includes describing patterns of community composition and abundance, identifying species that are invasive or of conversation value, improving sequence reference libraries, and gaining information on ecosystem function and environmental disturbance in a non-destructive way. The processing of ARMS is time-intensive and so we do not recommend that these tools supersede less time-intensive methods, such as eDNA sampling, or established methods such as benthic grabs. Rather we advocate for the use of ARMS as a tool to ground truth eDNA methods, provide important benthic data that other methods cannot, and establish trends in sensitive habitats (e.g. MPAs) where benthic grabs are not warranted.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

Ethical approval was not required for the study involving animals in accordance with the local legislation and institutional requirements because research involved the collection of surplus material for DNA analysis and did not involve the removal of material from any natural habitat.

Author contributions

MJ: Data curation, Formal Analysis, Investigation, Methodology, Visualization, Writing – original draft. MS: Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Writing – review & editing. NN: Data curation, Formal Analysis, Investigation, Methodology, Visualization, Writing – review & editing. KH: Methodology, Supervision, Writing – review & editing. RW: Investigation, Resources, Writing – review & editing. CY: Conceptualization, Funding acquisition, Investigation, Project administration, Supervision, Writing – review & editing. ER: Conceptualization, Funding acquisition, Investigation, Project administration, Supervision, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. The Sussex Kelp Recovery Project provided funds for deployment and recovery of ARMS.

Conflict of interest

The authors declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2025.1674917/full#supplementary-material

Supplementary Figure 1

NMDS of Image Analysis Functional Grouping Community Composition x Plate Face.

Supplementary Figure 2

Boxplots of Percentage Cover Between Plate Faces.

References

Summary

Keywords

barcoding, benthic biodiversity, kelp, marine, metabarcoding

Citation

Jordan M, Steyaert M, Ng NYN, Hopkins K, Ward RD, Yesson C and Ransome E (2026) Autonomous reef monitoring structures as a tool for assessing UK marine benthic biodiversity. Front. Mar. Sci. 12:1674917. doi: 10.3389/fmars.2025.1674917

Received

28 July 2025

Revised

26 November 2025

Accepted

27 November 2025

Published

02 February 2026

Volume

12 - 2025

Edited by

Pierre Yves Le Traon, Mercator Ocean, France

Reviewed by

Cory Douglas Bishop, St. Francis Xavier University, Canada

Hyunsung Yang, Korea Institute of Ocean Science and Technology (KIOST), Republic of Korea

Updates

Copyright

*Correspondence: Emma Ransome,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics