Data and code from: Estuarine benthic diversity declines with temperature across a tropical-temperate gradient
Data files
Aug 06, 2026 version files 117.75 MB
-
ASV_Taxa_sequence_18S_DNA_Latitude_edit.csv
117.21 MB
-
estuary_lat_long.csv
851 B
-
gradistat_latitude.csv
67.51 KB
-
LatitudeDNA_1.Data.R
61.48 KB
-
LatitudeDNA_2.ABND_DivAnalysis.R
115.76 KB
-
LatitudeDNA_2.RARAlphaDivAnalysis.R
107.94 KB
-
LatitudeDNA_3.CommAnalysis2.R
34.49 KB
-
LatitudeDNA_4.Supps.R
35.08 KB
-
Metadata_Latitude_ECOG.csv
116.82 KB
-
README.md
4.97 KB
Abstract
Estuarine benthic sediments are important habitats that support ecosystem health and productivity, but the drivers of large-scale patterns in benthic diversity are not well documented. Using 18S amplicon sequencing of 425 sediment DNA samples from estuaries spanning 20o of latitude along Australia’s east coast, we assessed tropical–temperate patterns in benthic eukaryotic richness, community composition, and environmental structuring. Contrary to classical expectations, ASV richness declined towards the tropics, with higher richness in temperate estuaries. This inverse gradient was evident for total eukaryotes, microeukaryotes, macrofaunal metazoans, and multicellular algae, while microbial metazoans and unicellular algae showed weaker or no latitudinal trends. Temperature was one of the strongest predictors of richness and community composition, although its high correlation with latitude means that thermal effects cannot be separated from other latitudinally structured processes. Sediment properties were also important, particularly sediment heterogeneity, which was consistently positively associated with richness. Nutrients and water-column physicochemical variables, including salinity, dissolved oxygen, and pH, had more group-specific effects on richness but contributed more clearly to shifts in community composition with latitude and taxon-level abundances. Over the latitudinal gradient, differences in habitat conditions may have contributed to these observed patterns. Compared with higher-latitude temperate estuaries, lower-latitude tropical estuaries tended to have more heterogeneous sediments but lower nutrient levels, lower dissolved oxygen, and lower pH, suggesting that warmer environments were associated with less favourable benthic habitat. These findings show that estuarine benthic eukaryotic diversity does not necessarily follow classical tropical-richness expectations, but instead reflects interactions between large-scale gradients, local sediment habitat conditions, and taxon-specific environmental responses. More broadly, our results suggest that warming and changes in sediment dynamics could reduce benthic biodiversity and compromise the habitat quality and functioning of estuarine ecosystems.
Dataset DOI: 10.5061/dryad.x3ffbg80j
Description of the data and file structure
This repository contains code and data to reproduce the results of the manuscript “Estuarine benthic diversity declines with temperature across a tropical-temperate gradient”. Data includes 18S amplicon sequencing of 425 sediment samples collected from 17 estuaries along Australia's east coast, covering both tropical and temperate locations.
Files and variables
File: Metadata_Latitude_ECOG.csv
Description: Field and laboratory metadata for estuarine sediment samples used to analyse environmental drivers of 18S rRNA ASV richness and community composition across ~20° of latitude in eastern Australia; each row corresponds to an individual sediment grab (unique Sample_ID), spanning 17 estuaries with multiple sites (1–5 per estuary) and technical replicates (n=5), and is intended to link environmental context to ASV richness, beta-diversity, and compositional change.
Variables
- Sample_ID: Sample identifier
- Year: Year of collection
- Study: Study (L=latitude)
- Code: Abbreviated estuary name
- Site: Site number (1-5)
- Rep: Rep number (1-5)
- Estuary: Estuary name
- Estuary_lat: Estuary latitude
- Estuary_long: Estuary longitude
- Site_lat: Site latitude
- Site_long: Site longitude
- Site_lat_dec: Site decimal latitude
- Site_long_dec: Site decimal longitude
- UTM_x: Universal Transverse Mercator x coordinates
- UTM_y: Universal Transverse Mercator y coordinates
- Date: Date of collection
- Time: Time of collection
- Depth_m: Depth of grab
- Temperature: Water temperature (oC)
- Salinity_ppt: Water salinity (ppt)
- Salinity_mS: Water salinity (mS)
- pH: pH
- Turbidity: Turbidity
- DO_pc: Dissolved oxygen (%)
- DO_mg: Dissolved oxygen (mg)
- OM_pc: Organic matter content (%)
- gravel_pc: Gravel content (%)
- sand_pc: Sand content (%)
- mud_pc: Silt content (%)
- TP: Total phosphorus (mg/kg)
- TN: Total nitrogen (mg/kg)
- TOC: Total organic carbon (%)
- TC: Total carbon (%)
- Al:Zn: Metal concentrations (mg/kg)
File: ASV_Taxa_sequence_18S_DNA_Latitude_edit.csv
Description: Amplicon sequence variant (ASV) table for 18S rRNA metabarcoding of estuarine sediments. Each row represents a unique ASV with its nucleotide sequence, taxonomic assignment, and per-sample read counts. The file supports community analyses of eukaryotic diversity and biogeographic patterns across a ~20° latitudinal gradient.
File: estuary_lat_long.csv
Description: This file provides latitude and longitude information for each sampled estuary, required as part of the map in Fig. 1.
Variables
- STATE: Australian state
- Estuary: Estuary name
- DLAT: Estuary latitude
- DLONG: Estuary longitude
File: gradistat_latitude.csv
Description: Sediment grain size analysis results from the program GRADISTAT (Blott & Pye, 2001)
Variables
- Sample_ID: Sample identifier
- Year: Year of collection
- Code: Estuary abbreviated code
- Site: Site number (1-5)
- Rep: Rep number (1-3)
- Gradistat_sample_no.: Gradistat identifier
- FW_um_mean: Folk and Ward (FW) mean grain size (µm)
- FW_um_sort: Sorting (µm)
- FW_um_skew: Skewness (µm)
- FW_um_kurt: Kurtosis (µm)
- FW_phi_mean: Mean grain size (phi units)
- FW_phi_sort: Sorting (phi units)
- FW_phi_skew: Skewness (phi units)
- FW_phi_kurt: Kurtosis (phi units)
- FW_desc_mean: Mean grain size (descriptive)
- FW_desc_sort: Sorting (descriptive)
- FW_desc_skew: Skewness (descriptive)
- FW_desc_kurt: Kurtosis (descriptive)
Code/software
All diversity calculations, statistical analyses, and figures were done using R version 4.4.0. The following packages are required to run each script:
- vegan
- tidyverse
- rnaturalearth
- ggspatial
- patchwork
- sf
- osmdata
- lme4
- lmerTest
- DHARMa
- piecewiseSEM
- dismo
- pheatmap
- cowplot
- gridExtra
- RColorBrewer
- gllvm
- psych
- Hmisc
LatitudeDNA_1.Data.R
Description: Must be run first (L11-L1704) before other scripts. Loads data and does basic data manipulation and calculates relative abundance data required for further analysis.
LatitudeDNA_2.RARAlphaDivAnalysis.R
Description: This script runs rarefied richness analyses and boosted regression tree analysis and produces figures 1-4.
LatitudeDNA_2.ABND_DivAnalysis.R
Description: This script runs relative abundance-based richness analyses and boosted regression tree analysis presented in the Supplementary material Appendix 3.
LatitudeDNA_3.CommAnalysis2.R
Description: This script runs beta diversity and multivariate gllvm analysis and produces figures 5 and 6.
LatitudeDNA_4.Supps.R
Description: This script runs analysis and produces figures S1-S8 as part of the Supplemental Material B.
