Consistent diversity loss but divergent functional responses of freshwater microbial communities to salinity and nutrient enrichment across climatic regions
Data files
Jul 17, 2026 version files 3.31 GB
-
16S_patterns_taxonomy.csv
6.84 MB
-
Functional_data.csv
9.63 KB
-
pond_data_for_16S.csv
7.70 KB
-
README.md
5.32 KB
-
Sequences.zip
3.30 GB
Aug 13, 2026 version files 3.35 GB
-
16S_patterns_taxonomy.csv
6.84 MB
-
ASV_factors.csv
5.43 KB
-
Fuggle_et_al_2026_Proc_B_code.Rmd
98.38 KB
-
Functional_data.csv
10.06 KB
-
pond_data_for_16S.csv
7.70 KB
-
README.md
6.87 KB
-
Sequences.zip
3.30 GB
-
Shannons_Data.csv
3.67 KB
-
TAXA_16S_GTDB.csv
36.71 MB
Abstract
Freshwater salinisation and nutrient enrichment are intensifying globally. Despite their widespread co-occurrence the extent to which their combined ecological effects are generalisable across distinct climatic regions remains unclear, preventing appropriate scales for policy implementation to be identified. Using outdoor 1,000 L open freshwater ponds established > 10 years ago, replicated across inland Mediterranean, coastal Mediterranean and mountain regions, we examined the independent and interactive effects of salinity and nutrient enrichment across regions using benthic biofilms as a model system. We tested the effects of elevated salinity (0, +5, +15 g NaCl L -1) and nutrient enrichment (10:1 N:P additions) on freshwater benthic bacterial community structure (richness, diversity, taxonomic structure) and whole-community biofilm functions after 30 and 60 days. Across all regions, salinity was the dominant driver of bacterial community structure and function with effects strengthening through time. Nutrient enrichment had no significant main or interactive effects, possibly due to the strong nutrient limitations in the systems. After 30 days, the structure of high and ambient salinity bacterial communities differed significantly at all regions, and by day 60 all salinity levels differed significantly from each other. Bacterial richness and diversity (Shannon’s index) also declined significantly with increasing salinity at all regions. Functional responses to salinity, however, varied by region. In the inland Mediterranean region, biofilm biomass, gross primary productivity (GPP) and net primary productivity (NPP) were significantly higher under high salinity compared to moderate and ambient salinity; in the coastal Mediterranean region, GPP and NPP were lower under high salinity compared to moderate and ambient salinity but biomass did not differ; while in the mountain region salinity had no effect on the measured functional variables. Photosynthetic efficiency was unaffected by treatment in any region. Policy implications: In sum, salinisation consistently decreased freshwater bacterial diversity in biofilms across three distinct climatic regions, whereas functional outcomes were regionally contingent. This suggests management policymakers can set cross-regional salinity thresholds to maintain bacterial diversity, but functional responses to salinisation must be assessed regionally to predict and manage the impacts of freshwater salinisation under accelerating global change.
Principle Investigator Contact Information
Name: Rose Fuggle
Institution: University of Sydney
Email: rose.fuggle@sydney.edu.au
Alternate Contact Information
Name: Ziggy Marzinelli
Institution: University of Sydney
Email: e.marzinelli@sydney.edu.au
Dataset Overview
This dataset contains the data in Fuggle et al (in review), testing the hypothesis that the structure and function of freshwater benthic microbial communities are affected by salinity and nutrient enrichment consistently across three regions. Data was obtained from benthic microbial communities grown on tiles in outdoor mesocosms in three regions: Jaca (mountain), Porto (coastal Mediterranean) and Evora (inland Mediterranean). The dataset includes 1) biofilm in situ functional metrics, including total chlorophyll a concentration in biofilms, photosynthetic yield, gross primary productivity (GPP), net primary productivity (NPP), community respiration (CR) 2) 16S rRNA ASV sequences from swabbed biofilms, 3) pond water data including oxygen concentrations, chlorophyll a, pH, conductivity, nitrate and phosphate.
Dates of Data Collection
Data was collected in May - July 2023.
Data Spatial Scope
Data was collected from an experimental field sites in Jaca (Spain), Porto (Portugal), and Evora (Portugal) as part of the Iberian Pond Network infrastructure.
Funding
This study has received funding from the AQUACOSM-plus project of the European Union's Horizon 2020 research and innovation programme under grant agreement No. 871081 and Transnational Access project SalNutFun. Additional funds were received from The University of Sydney WG Murrell funds.
Ethics Approval
N/A
Files and Folders
Folder:Sequences.zip
This is a zip folder of all raw fasta files for each sample. Samples are named as [Site][Treatment][Pond][Time][Read].fasta. Treatment group abbreviations differ from those in the manuscript and are: LS = low salinity, MS = moderate salinity, HS = high salinity, LSHN = low salinity, high nutrients, MSHN = moderate salinity, high nutrients, HSHN = high salinity, high nutrients.
File:pond_data_for_16S.csv
This is a csv file with the pond water data. Columns are:
- Site: [JA = Jaca, PO = Porto, EV = Evora] this is the site where data were collected.
- Pond_Number: [Site][PondID] This is the unique pond identifer
- Treatment: Pond salinity and nutrient treatment. Acroynyms match those in the manuscript and are: HSAN = high salinity, ambient nutrients. HSEN = high salinity, enriched nutrients. MSAN = moderate salinity, ambient nutrients. MSEN = moderate salinity, enriched nutrients. ASAN = Ambient salinity, ambient nutrients. ASEN = ambient salinity, enriched nutrients.
- Pond: Pond number that corresponds to 16S Pond_ID
- Sample_Point: Experiment day data were collected (T0 = day 30 and T2 = day 60)
- Conductivity: Conductivity in water in mS/cm
- Chlorophyll_a: Chlorophyll a concentration in water in in ug/L
- pH: pH in water
- Oxygen_Concentration: Dissolved oxygen concentration in water in mg O2/L
- Nitrite: Nitrite concentration in water in in mg/L. Missing data on Day 60 is because nutrient data was not collected for this timepoint.
- Nitrate: Nitrate concentration in water in in mg/L. Missing data on Day 60 is because nutrient data was not collected for this timepoint.
- Phosphate: Phosphate concentration in water in in mg/L. Missing data on Day 60 is because nutrient data was not collected for this timepoint.
File: Functional_data.csv
This is a csv file with the univariate data collected in situ from the biofilms. Columns are:
- Sample_Name: [Site][Treatment][Pond-replicate]_[Time]
- Location: I = Inland, M = Mountain, C = Coastal.
- Pond_Replicate: Pond replicates (P1-P4) within each location
- Sample_Point: T1 = Day 30, T2 = Day 60
- PAM_Yield: A unitless value representing photosynthetic efficiency
- Concentration_Green_Algae: Concentration of green algae in biofilms in the nits ug chl-a/cm2
- Concentration_Cyanobacteria: Concentration of cyanobacteria in biofilms in the units ug chl-a/cm2
- Concentration_Diatoms: Concentration of diatoms in biofilms in the units ug chl-a/cm2
- CR: Community Respiration rates in mg O2/L/hr. Missing values for samples where Location = M, Sample_Point = T2 and Treatment = ASAN, MSAN and MSEN are due to equipment failure.
- NPP: Net Primary Productivity rates in mg O2/L/hr. Missing values for samples where Location = M, Sample_Point = T2 and Treatment = ASAN, MSAN and MSEN are due to equipment failure.
- GPP: Gross Primary Productivity rates in mg O2/L/hr. Missing values for samples where Location = M, Sample_Point = T2 and Treatment = ASAN, MSAN and MSEN are due to equipment failure.
File: 16S_patterns_taxonomy.csv
This is a csv file with the ASV counts for each sample merged with their taxonomic identifiers from the Genome Taxonomy Database (GTDB). Columns are:
- ASV: The genetic sequence of the amplicon
- Columns B through to EL: Unique biofilm samples named in the format [Site][Treatment]_ [Pond]\[Time]
- Columns EM through to ES: Taxonomic identifiers for each ASV from the levels Kingdom, Phylum, Class, Order, Family, Genus, Species.
Note: blanks are used to represent absent or missing data
File: Shannons_Data.csv
This is a csv file that contains the Shannon's diversity index for each sample. Columns are:
- Sample_Name: Unique biofilm samples named in the format [Site]_[Treatme**nt]_Pond][Time]
- H: Shannon's diversity index.
File: TAXA_16S_GTDB.csv
This is a csv file that assigns taxonomy to each ASV in the complete dataset. Columns are:
- ASV: Unique amplicon sequence variant
- Columns B through to H: Taxonomic identifiers for each ASV from the levels Kingdom, Phylum, Class, Order, Family, Genus, Species.
- Columns I through to O: Bootstraps for the relevant taxonomic identifier.
File: ASV_factors.csv
This is a csv file that assigns the factors to each sample for the ASV analyses. Columns are:
- Sample: Unique biofilm samples named in the format [Site]_[Treatment]_Pond][Time]
- Site: EV = Evora, JA = Jaca, PO = Porto
- Treatment: As referenced in the manuscript, the salinity and nutrient treatment
- Pond: Unique pond environment labeled [Site][Treatment][Pond_replicate]
- Time: Sample point T0 = 30 day, T2 = day 60
- Sal: Salinity treatment
- Nut: Nutrient treatment
File: Fuggle_et_al_2026_Proc_B_code.Rmd
This is an R Markdown file containing the code for all analyses and figure generation. To run this script, download the files contained within this Dryad repository and move them into an R project containing this R markdown script.
Changes after Jul 17, 2026:
Addition of R Markdown file containing code for analyses and figures.
Addition of datasets to enable complete running of R Markdown script.
