Genetic diversity affects ecosystem functions across trophic levels as much as species diversity, but in an opposite direction

  1. Laura Fargeot  Is a corresponding author
  2. Camille Poesy
  3. Maxim Lefort
  4. Jerome G Prunier
  5. Madoka Krick
  6. Rik Verdonck
  7. Charlotte Veyssiere
  8. Murielle Richard
  9. Delphine Legrand
  10. Geraldine Loot
  11. Blanchet Simon  Is a corresponding author
  1. Centre National de la Recherche Scientifique (CNRS), Station d’Ecologie Théorique et Expérimentale, France
  2. Université Paul Sabatier (UPS) Toulouse III, France
  3. Center for Environmental Sciences, Environmental Biology, Hasselt University, Belgium
5 figures, 2 tables and 3 additional files

Figures

Location of sampling sites and illustration of the trophic chain.

(a) Distribution of the 52 sampling sites (black dots) spanning an east-west gradient at the foothills of the Pyrenees Mountains (France). Each site is denoted by a six-letter code, with the first three letters indicating the river name and the last three letters indicating the closest city or village. (b) Our study focused on a tri-trophic food chain commonly found in mountain rivers, consisting of riparian trees, macroinvertebrate shredders, and fishes (from bottom to top). Within each trophic level, we measured two facets of biodiversity: genetic diversity in a single target species within each trophic level (specifically P. dragarum, Gammarus sp., and A. glutinosa), and species diversity from communities.

Example of one of the seven causal models used to quantify the relationships between (species and genetic) diversity and ecosystem functions.

We focused on seven ecosystem functions associated with genetic and species diversity at three trophic levels (green for primary producer, orange for primary consumer, and blue for secondary consumer). Each relationship between biodiversity and ecosystem functions (n=6 values per function, but for some functions for which irrelevant links were not considered, see the text, n=32 values in total) was measured at the same trophic level (triangles) or at another trophic level (dots).

Details of the seven causal models linking abiotic parameters, species and genetic diversity, and ecosystem functions.

Each causal graph (a–g) represents a simplified illustration of the relationships between the two principal component analysis (PCA) axes synthesizing the environmental parameters of each sampling site (Environmental axis 1 and 2), the species diversity estimated at each trophic level (boxes ‘Fish’, ‘Shredders’, and ‘Trees’), the genomic diversity estimated from each focal species at each trophic level (boxes ‘Phoxinus’, ‘Gammarus’, and ‘Alnus’), and each ecosystem function (one model per function). Only the relationships for which the p-value was inferior to 0.20 are indicated for visual simplification. Full arrows indicated positive effects, whereas dotted arrows indicated negative effects. The width of the arrows is proportional to the size of their effects. The percentage of variance explained by environmental and biodiversity effects on ecosystem functions (r2) is indicated for each function.

General description of individual effect sizes measured between biodiversity estimates and ecosystem functions (BEFs) in a riverine trophic chain.

(a) The magnitude and direction of individual effect sizes (Zr) of biodiversity is shown for each ecosystem functions as a biplot between Zr associated with genetic diversity (y-axis, genetic BEFs) measured for one of three target species (A. glutinosa, Gammarus sp., and P. dragarum) and Zr associated with species diversity (x-axis, species BEFs) measured for one of three trophic levels (trees, invertebrates, and fish). For each ecosystem function (but the biomass of trees and of A. glutinosa), a total of six Zr are depicted in the biplot; four of them are associated with biodiversity measured at another trophic level than the one of the target functions (red symbols, e.g. effect of fish diversity on invertebrate biomass) and two of them are associated with biodiversity measured at the same trophic level than the one of the target functions (blue symbols, e.g. effect of fish diversity on fish biomass). The arrows indicate significant Zr (95% confidence intervals excluded 0, see Supplementary file 2); vertical arrows are for significant genetic BEFs, horizontal arrows are for species BEFs. White quadrats stand for situation in which genetic and species BEFs are in the same direction, whereas gray quadrats indicate situation in which genetic and species BEFs are in the opposite direction. Within each quadrat, sub-quadrats indicate the relative magnitude of BEFs, i.e., whether genetic BEFs are stronger, weaker, or equal in magnitude than species BEFs. (b) Density plots displaying the distribution of individual Zr for species and genetic BEFs (dotted and full lines respectively).

Magnitude and direction of the mean effects sizes estimated from the relationships between biodiversity and ecosystem functions (BEFs) measured in a riverine trophic chain.

(a) The magnitude and direction of BEFs are expressed as effect sizes (Zr) and are displayed according to the facet used to measure biodiversity (genetic or species diversity, light gray and white boxplots respectively) and to the type of BEFs (within-trophic level BEFs or between-trophic level BEFs, triangles and dots respectively). Red color and stars indicate global effect sizes that are significantly different from zero (p-value<0.05). Large symbols are mean ± 95%IC estimated as marginal effects from the meta-regressions. Small symbols are raw estimates. (b) Same representation as (a) but with details at each trophic level (mean ± 95% IC estimated as marginal effects from the meta-regressions, green for primary producer, orange for primary consumer, and blue for secondary consumer). The trophic level at which BEFs are measured is coherent across all trophic levels.

Tables

Table 1
Characteristics of the first two principal components identified by the principal component analysis (PCA) ran on the 13 environmental variables.

The part of the total environmental variance (%) and the contribution of each variable on each component are shown. The variables that contributed significantly to the axis are highlighted in bold.

Component 1Component 2
Part of total variance (%)21.6616.37
River width0.5960.320
Connectivity–0.1550.646
Altitude0.648–0.385
Distance from outlet0.528–0.331
East-west gradient0.105–0.795
Oxygen concentration0.7380.343
Oxygen saturation0.266–0.509
Water temperature–0.594–0.012
Specific conductivity–0.463–0.045
pH0.5730.398
Concentration in NO3+NO2–0.369–0.286
Concentration in NH4+–0.019–0.207
Concentration in PO43-–0.2790.236
Global characteristicLow altitude, poorly oxygenated site – high altitude, highly oxygenated sitesPoorly connected east site – highly connected west site
Table 2
ANOVA table for the linear mixed model testing whether the relationships between biodiversity and ecosystem functions (BEFs) measured in a riverine trophic chain differ between the biodiversity facets (species or genetic diversity) and the types of BEF (within- or between-trophic levels).

A Wald chi-square test is used to test the significance of each fixed effect.

Degree of freedomChisq-valuep-Value
(Intercept)10.2870.595
Biodiversity facet10.2320.630
Type of BEF15.3930.020
Biodiversity facet*Type of BEF15.5670.018

Additional files

Supplementary file 1

ANOVA table for the linear mixed model testing whether the relationships between biodiversity and ecosystem functions (BEFs) measured in a riverine trophic chain differ between the biodiversity facets (species or genetic diversity), the types of BEF (within- or between-trophic levels), and the trophic levels at which BEFs are estimated (primary producers, primary consumers, or secondary consumers).

A Wald chi-square test is used to test the significance of each fixed effect.

https://cdn.elifesciences.org/articles/100041/elife-100041-supp1-v1.docx
Supplementary file 2

Estimates of individual effect sizes of biodiversity and ecosystem functions (BEFs) (Zr, n=34) for each ecosystem function, each biodiversity facet (genetic or species diversity), and each type of BEF (within- or between-trophic levels).

95% confidence intervals are provided together with the estimate of each BEF. BEFs are considered as significant when the 95% CI does not overlap 0. p-Values estimated from t-test are also provided.

https://cdn.elifesciences.org/articles/100041/elife-100041-supp2-v1.docx
MDAR checklist
https://cdn.elifesciences.org/articles/100041/elife-100041-mdarchecklist1-v1.docx

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Laura Fargeot
  2. Camille Poesy
  3. Maxim Lefort
  4. Jerome G Prunier
  5. Madoka Krick
  6. Rik Verdonck
  7. Charlotte Veyssiere
  8. Murielle Richard
  9. Delphine Legrand
  10. Geraldine Loot
  11. Blanchet Simon
(2025)
Genetic diversity affects ecosystem functions across trophic levels as much as species diversity, but in an opposite direction
eLife 13:RP100041.
https://doi.org/10.7554/eLife.100041.4