Data and code from: When giant amphibians rule the river: Community trophic structure of Japanese giant salamander ecosystems
Data files
Jun 02, 2026 version files 97.93 KB
-
data_trophic_ecology_community_giant_salamanders.zip
92.07 KB
-
README.md
5.86 KB
Abstract
Trophic interactions among species determine the structure and functioning of food webs and ultimately govern ecosystem organization. Top predators, at the highest trophic levels, play a pivotal role in shaping food webs. Yet, the ecological consequences of body-size–related trophic shifts in large predators are poorly explored in freshwater ecosystems, despite being well documented in marine and terrestrial systems. This gap is particularly evident in river ecosystems dominated by giant salamanders. We investigated the trophic ecology of a freshwater community dominated by the Japanese giant salamander (Andrias japonicus) to assess how extreme ontogenetic body-size variation modifies intra- and interspecific interactions and contributes to overall community structure. We also aimed to identify the dominant carbon sources fueling the community by estimating the relative contributions of terrestrial versus aquatic energy pathways. Using stable isotope analysis, we quantified trophic positions and isotopic niche metrics for all major consumer species. Our results reveal clear trophic shifts in A. japonicus: smaller individuals overlapped with mesopredators (teleost fish, prawns, turtles), suggesting potential competition, while trophic position increased progressively with size, placing large individuals at the top of the food chain. At the community level, the presence of large salamanders broadened the δ¹⁵N range, increased trophic evenness, and expanded isotopic space, enhancing vertical complexity. Basal source tracing indicated complementary reliance on aquatic and terrestrial carbon across taxa, revealing multiple energy pathways sustaining the food web. These findings show that gigantism in A. japonicus reshapes trophic hierarchies by reducing redundancy and strengthening trophic stratification. More broadly, they underscore the ecological importance of ontogenetic niche shifts and body size in structuring freshwater food webs, highlighting how such shifts cascade to influence ecosystem functioning. By integrating isotopic niche metrics with basal source contributions, our study provides a generalizable framework for linking predator ontogeny, energy pathways, and food web structure, with direct relevance for conserving vulnerable freshwater ecosystems. Moreover, it highlights that endangered top predators such as A. japonicus may contribute to food web resilience and ecological integrity in vulnerable freshwater ecosystems.
All analyses were performed in R version 4.3.3. Required packages are listed and loaded at the top of each script.
Main folder: data_trophic_ecology_community_giant_salamanders.zip
Folder 1: Clustering of Japanese giant salamander data
In this folder, data and scripts are provided to perform the clustering of Japanese giant salamanders into 4 size-based groups.
The file "cluster_giant_salamanders.csv" contains the following data:
- column 'SVL': body size of each individual (cm)
- column 'd13C': δ13C values (‰)
- column 'd15N': δ15N values (‰)
Folder 2: PLot isotopic ellipses 4 giant salamanders groups
In this folder, data and scripts are provided to compute and plot the isotopic niches of the 4 giant salamander groups obtained from the clustering procedure.
The file "SIBER_4groups_salamanders.csv" contains the following data:
- column 'iso1': δ13C values (‰)
- column 'iso2': δ15N values (‰)
- column 'group': size group of Japanese giant salamander (resulting from the clustering procedure)
- column 'community': community in which the group is located (here only 1 community)
Folder 3: Plot isotopic ellipses all groups
In this folder, data and scripts are provided to compute and plot the isotopic niches of all consumers of the community. There are 2 parts: (1) the full community with all fish pooled as a single group “fish” (Figure 5a), and (2) a zoom on the subset of species showing niche overlap, with fish considered separately (Figure 5b).
There are 4 data files provided:
- "SIBER_community_global.csv": file containing isotopic data of all consumers to build isotopic ellipses (Figure 5a)
- column 'iso1': δ13C values (‰)
- column 'iso2': δ15N values (‰)
- column 'group': consumers (see next file for correspondence between group numbers and species)
- column 'community': community in which is group is located (here only 1 community)
- "key_community.csv": file containing the identity of each of the 12 groups in the previous file
- column 'Group': group number
- column 'Name': name of the species or group of species
- "SIBER_community_zoom.csv": file containing isotopic data to build ellipses for the overlapping species, considering each fish species separately (Figure 5b)
- column 'iso1': δ13C values (‰)
- column 'iso2': δ15N values (‰)
- column 'group': consumers (see next file for correspondence between group numbers and species)
- column 'community': community in which is group is located (here only 1 community)
- "key_community_zoom.csv": file containing the identity of each of the 12 groups in the previous file
- column 'Group': group number
- column 'Name': name of the species or group of species
Folder 4: Isotopic niche width and overlap
In this folder, data and scripts are provided to calculate standard ellipse area (SEA) and the niche overlap between the different groups.
There are 2 data files provided:
- "SIBER_global.csv": file containing isotopic data of all consumers to calculate and plot SEA for each group (Figure 3)
- column 'iso1': δ13C values (‰)
- column 'iso2': δ15N values (‰)
- column 'group': consumers (see next file for correspondence between group numbers and species)
- column 'community': community in which is group is located (here only 1 community)
- "key_global.csv": file containing the identity of each of the 17 groups in the previous file
- column 'Group': group number
- column 'Name': name of the species or group of species
Folder 5: Layman's metrics for each species group through boostrapping
In this folder, data and scripts are provided to calculate Layman's metrics for each species group through bootstrapping and to create the corresponding figure containing 3 Layman's metrics (Figure 4).
The file "data_boot.csv" contains isotopic data of all consumers to calculate and plot the 3 Layman's metrics:
- column 'iso1': δ13C values (‰)
- column 'iso2': δ15N values (‰)
- column 'group': consumers (see the file named "key_global.csv" in the previous part for the correpondance number/species)
- column 'ID': identification number of each individual of each species group
Folder 6: Layman's metrics at the community level
In this folder, data and scripts are provided to calculate Layman's metrics at the community level under four scenarios of Japanese giant salamander inclusion (Figure 6).
The files used here are the same than in Folder 4.
Folder 7: Trophic position and proportional contribution of terrestrial carbon sources (α)
In this folder, data and scripts are provided to calculate and plot trophic position and proportional contribution of terrestrial carbon sources (α) of each group.
There are 2 data files provided:
- "data_trophic_position.csv" contains the following data:
- column 'community' corresponds to the species group
- column 'Species': within each species group indicates whether the row corresponds to a consumer (e.g., “aj1” for Japanese giant salamander group 1), “lit” (leaf litter; terrestrial baseline), or “peri” (periphyton; aquatic baseline)
- column 'FG' standing for functional group: "cons" = consumer; "aquatic_BL" = aquatic baselines; "terrestrial_BL" = terrestrial baselines
- column 'd13C': δ13C values
- column 'd15N': δ15N values
- "TP_Alpha_results.csv" contains the results of trophic position and alpha to produce the plot
- column 'species': species group
- column 'metric': TP (trophic position) or α (alpha)
- column 'mode': value of the estimated mode for each group
- columns 'X25', 'X125', 'X250', 'X750', 'X875', and 'X975': credible interval bounds (respectively 2.5%, 12.5%, 25%, 75%, 87.5%, 97.5%)
