Depth-structured sponge assemblages offer limited evidence for mesophotic refuge potential under the deep reef refuge hypothesis
Data files
Feb 02, 2026 version files 136.98 KB
-
PK_all_sites_sponges.xlsx
55.43 KB
-
README.md
3.58 KB
-
Sponge_assemblage_stats_and_plots_ecography.R
77.97 KB
Abstract
Dataset DOI: 10.5061/dryad.7m0cfxq9b
Description of the data and file structure
This repository contains the data and R code used to analyse sponge assemblage structure across depth gradients at the Poor Knights Islands Marine Reserve, New Zealand. The analyses assess sponge assemblage composition, diversity, and depth-related abundance patterns in shallow and mesophotic reef zones, with the aim of evaluating the potential for mesophotic ecosystems to act as ecological refuges for shallow sponge species.
The dataset accompanies the manuscript:
Depth-structured sponge assemblages offer limited evidence for mesophotic refuge potential under the deep reef refuge hypothesis
Files and variables
File: PK_all_sites_sponges.xlsx
Description: Excel file containing sponge percentage cover data derived from benthic image analysis.
Sheets:
NAW– Northern Arch WestNAE– Northern Arch EastIP– Imagination Point
Variables:
Sample– Unique identifier for each image/quadratSite– Sampling siteDepth– Sampling depth (m)Zone– Depth zone (Shallow<30 m,Mesophotic>30 m)- Columns 5 onward – Sponge species or operational taxonomic units (OTUs), expressed as percentage cover per image
Notes:
- Percentage cover values represent proportional cover of the benthos per image.
- OTUs include both identified species and morphospecies where taxonomic resolution was not possible from imagery.
File: Sponge_assemblage_stats_and_plots_ecography.R
Description: R script containing all data processing, statistical analyses, and figure generation used in the manuscript.
Main analyses include:
- Data import and formatting
- Non-metric multidimensional scaling (nMDS) using Bray–Curtis dissimilarities
- PERMANOVA and PERMDISP analyses testing for differences among depths and zones
- Beta diversity partitioning (turnover vs. nestedness) using Jaccard dissimilarity
- Linear regression analyses of total sponge cover and depth-generalist sponge abundance
- Figure generation for assemblage ordinations, dispersion plots, and diversity summaries
- Export of PERMANOVA and PERMDISP results to Excel files
Reproducibility notes:
- A fixed random seed (
set.seed(123)) is used for permutation-based analyses. - All analyses were conducted in R version 4.4.0.
- File paths in the script may need to be updated to match the user’s local directory structure.
Code/software
Key R packages used:
tidyverse, vegan, pairwiseAdonis, betapart, ggplot2, patchwork, lme4, emmeans, among others.
Relationship Between Data, Code, and Manuscript
- The Excel file (
PK_all_sites_sponges.xlsx) contains the raw percentage cover data used for all analyses. - The R script (
sponge assemblage stats and plots ecography.R) imports these data and performs all statistical analyses reported in the manuscript. - Figures and statistical outputs generated by the script correspond to figures and tables in the main text and supplementary materials.
How to Run the Analyses
- Download all files to a single working directory.
- Open
sponge assemblage stats and plots ecography.Rin R or RStudio. - Update the
setwd()paths if necessary. - Ensure all required R packages are installed.
- Run the script sequentially from top to bottom.
