# Sea lamprey nests promote the diversity of benthic macroinvertebrate assemblages

Marius Dhamelincourt, Jacques Rives, Marie Pons, Aitor Larrañaga, Cédric Tentelier, Arturo Elosegi

<figure><img src="image_1.png" /></figure>

## To cite this version:

Marius Dhamelincourt, Jacques Rives, Marie Pons, Aitor Larrañaga, Cédric Tentelier, et al. Sea lamprey nests promote the diversity of benthic macroinvertebrate assemblages. PLoS ONE, 2022, 17 (12), pp.e0274719. 10.1371/journal.pone.0274719. hal-04028896

HAL Id: hal-04028896

https://hal.inrae.fr/hal-04028896

Submitted on 14 Mar 2023

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Sea lamprey nests promote the diversity of benthic macroinvertebrate assemblages

Marius Dhamelincourt<sup>1,2*</sup>, Jacques Rives<sup>1,2</sup>, Marie Pons<sup>3</sup>, Aitor Larrañaga<sup>4</sup>, Cédric Tentelier<sup>1,2</sup>, Arturo Elosegi<sup>4</sup>

1 Université de Pau et des Pays de l'Adour, E2S UPPA, INRAE, ECOBIOP, Saint-Pée-sur-Nivelle, France,
2 Pôle Gestion des Migrateurs Amphihalins dans leur Environnement, OFB, INRAE, Agrocampus Ouest, Univ Pau & Pays Adour/ E2S UPPA, Pau, France, 3 Bureau d'études Aquabio, Route de Créon, Saint-Germain-du-Puch, France, 4 Department of Plant Biology and Ecology, University of the Basque Country (UPV/EHU), Bilbao, Spain

* marius.dhamelincourt@inrae.fr

<figure><img src="image_3.png" /></figure>

## Abstract

The habitat heterogeneity hypothesis states that increased habitat heterogeneity promotes species diversity through increased availability of ecological niches. We aimed at describing the local-scale (i.e. nest and adjacent substrate) effects of nests of the sea lamprey (*Petromyzon marinus L.*) as ecosystem engineer on macroinvertebrate assemblages. We hypothesized that increased streambed physical heterogeneity caused by sea lamprey spawning would modify invertebrate assemblages and specific biologic traits and promote reach-scale diversity. We sampled thirty lamprey nests of the Nive River, a river of the south western France with a length of 79.3 km and tributary of the Adour River, in three zones: the unmodified riverbed (upstream) and zones corresponding to the nest: the area excavated (pit) and the downstream accumulation of pebbles and cobbles (mound). The increased habitat heterogeneity created by lamprey was accompanied by biological heterogeneity with a reduced density of invertebrates (3777 ± 1332 individuals per $m^2$ in upstream, 2649 ± 1386 individuals per $m^2$ in pit and 3833 ± 1052 individuals per $m^2$ in mound) and number of taxa (23.5 ± 3.9 taxa for upstream, 18.6 ± 3.9 taxa in pit and 21.2 ± 4.5 taxa for mound) in the pit compared to other zones. However the overall taxa diversity in nest increased with 82 ± 14 taxa compared to the 69 ± 8 taxa estimated in upstream zone. Diversity indices were consistent with the previous results indicating a loss of $\alpha$ diversity in pit but a higher $\beta$ diversity between a pit and a mound than between two upstream zones, especially considering Morisita index accounting for taxa abundance. Trait analysis showed high functional diversity within zones with a reduced proportion of collectors, scrapers, shredders, litter/mud preference and small invertebrates in mound, while the proportion of “slabs, blocks, stones and pebbles” preference and largest invertebrates increased. Pit presented the opposite trend, while upstream had globally intermediate trait proportions. Our results highlight important effects on species and functional diversity due to habitat heterogeneity created by a nest-building species, what can ultimately influence food webs and nutrient processes in river ecosystems.

<figure><img src="image_4.png" /></figure>

## OPEN ACCESS

**Citation:** Dhamelincourt M, Rives J, Pons M, Larrañaga A, Tentelier C, Elosegi A (2022) Sea lamprey nests promote the diversity of benthic macroinvertebrate assemblages. PLoS ONE 17(12): e0274719. https://doi.org/10.1371/journal.pone.0274719

**Editor:** Michael A Chadwick, King's College London, UNITED KINGDOM

**Received:** September 2, 2022

**Accepted:** November 16, 2022

**Published:** December 15, 2022

**Peer Review History:** PLOS recognizes the benefits of transparency in the peer review process; therefore, we enable the publication of all of the content of peer review and author responses alongside final, published articles. The editorial history of this article is available here: https://doi.org/10.1371/journal.pone.0274719

**Copyright:** © 2022 Dhamelincourt et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** Data and code to run the analysis are available in the INRAE dataverse repository (https://doi.org/10.15454/BEHOOK) under the name "Sea lamprey nests promote the
diversity of benthic macroinvertebrate assemblages".

**Funding:** Functioning was financed by Pôle Gestion des Migrateurs Amphihalins dans leur Environnement (https://www6.rennes.inrae.fr/u3e/PRESENTATION/Organisation/Pole-MIAME). M.D. PhDs was financed by University of Pau and Pays de l'Adour (https://www.univ-pau.fr/fr/index.html) and UPV/EHU (https://www.ehu.eus/es/home). Field work used resources from the IE ECP Experimental Facility of the UMR Ecobiop (https://www6.bordeaux-aquitaine.inrae.fr/ie-ecp-ecobiop). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Grant no. IT1471-22 from the Basque Government funded this study.

**Competing interests:** The authors have declared that no competing interests exist.

Introduction

Habitat heterogeneity has been widely recognized as an important factor structuring biological assemblages [1, 2]. The habitat heterogeneity hypothesis [3] states that increased habitat heterogeneity promotes species diversity thanks to a greater number of available ecological niches. Li and Reynolds [4] defined heterogeneity in an operational way as "the complexity and/or variability of a system property in space and/or time". Heterogeneity can be considered at different scales, and depending on the scale considered, it has contrasted effects on ecological processes [5]. In the case of stream invertebrates, substrate heterogeneity can affect food distribution [6] and thus condition their local-scale (i.e. nest-scale) feeding and movements [7], whereas movements at larger scales can respond to drift if there are no available interstices or substrate for attachment.

In aquatic ecosystems, heterogeneity analyzes are often used to explain macroinvertebrate assemblages. Beisel et al. [8] found that the number of macroinvertebrate taxa in a given area (considering 0.5 to 4 m radius zones of substrate in a fourth-order stream) was higher in a heterogeneous environment with more substrate types and an elevated patchiness. Furthermore, the same study found that homogeneous environments are more prone to the dominance of few taxa, probably due to the lack of competition with taxa from neighbouring patches. At a very fine scale, Boyero [9] found a positive relationship between substrate heterogeneity and taxonomic diversity across 225 cm² samples collected within the same square meter of riffler. Substrate complexity appears to be a key driver of macroinvertebrate assemblage structure followed by current velocity and depth [10]. However, increased heterogeneity does not necessarily enhance macroinvertebrate diversity, at least in stream restoration projects [11].

Ecosystem engineering [12, 13] is defined as the change of availability of resources to other species due to physical changes in biotic or abiotic materials, directly or indirectly mediated by organisms. During the process they modify the habitats or create new ones. Effects of engineer species [14], alongside with the interaction between hydraulics and the substrate, are abiotic and biotic processes influencing heterogeneity. Together, they affect habitat heterogeneity, which in turn can affect ecosystem functioning [6]. Substrate disturbance by foraging fish is another cause of substrate heterogeneity, which may destabilize the substrate and cause an increased bedload transport during subsequent flood events [15].

Fish nesting, especially in species such as salmonids, which move large volumes of sediment and alter local bed morphology, is a biotic process that can increase riverbed substrate heterogeneity [16]. For migratory species with high nest density in the spawning grounds, the ecological effects can be important and cascade to affect biological assemblages and ecosystem processes. For instance, Moore and Schindler [17, 18] found that bed disturbance by salmon in spawning grounds caused severe seasonal declines in periphyton and benthic invertebrate abundance, with no recovery within the same season. However, other species such as chubs of the genus *Nocomis* build nests that promote macroinvertebrate density and even allow some taxa to persist in degraded streams [19]. Migratory species likely to affect substrate heterogeneity and the assemblages include the sea lamprey (*Petromyzon marinus* L.), which builds 40–220 cm-wide nests [20–22] in late spring ([23], Fig 1). Both male and female lamprey remove cobbles with their mouth and release them downstream so that the resulting nest consists of a pit carpeted with sand or gravel, followed by a mound of pebbles and cobbles. This structural heterogeneity might have functional consequences, as the density of some taxa such as Hydropsychidae, Philopotamidae and Heptageniidae can increase up to tenfolds in mounds [24] in an oligotrophic stream. Another study [25] found an increased Simuliidae density in areas with disturbed substrate and added sea lamprey carcasses, possibly resulting from the association of the cleared substrate with suitable characteristics and food supply from dead lampreys.
<figure><img src="image_5.png" /><figcaption>**Fig 1.** Photo of two sea lamprey nests in the Nive River. The mound of cobbles, immediately downstream of the pit, is clearly visible. ©INRAE-GLISE.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0274719.g001

These results suggest an effect of sea lamprey nesting activity on macroinvertebrate assemblages, with potential functional consequences.

Here, we aimed at describing the local-scale (i.e. nest and adjacent substrate) effects of sea lamprey as ecosystem engineer on macroinvertebrate assemblages at the nest. We expect the taxonomic and functional diversity on the nest as a whole (pit+mound) to be higher than on undisturbed substrate. We hypothesized that increased streambed physical heterogeneity caused by sea lamprey nesting would promote invertebrate diversity and specific biologic traits. Pits, with deeper and more lentic conditions than unmodified riverbed (upstream) would attract lentophilic and fine sediment fauna, whereas the shallower and more lotic mounds would attract taxa with opposite preferences. In addition, nests with more heterogeneous current and depth may host more diverse communities. To test our hypotheses, we first assessed the variability of depth and current velocity on pit, mound and upstream as these variables are among the most likely to explain differences in macroinvertebrate assemblages [10]. Then, we compared the density, taxa richness and diversity of macroinvertebrates across nest zones. Finally, we studied ecological traits describing food type, mode of alimentation, substrate preference and body size, to infer the effects of nests on invertebrate functional heterogeneity.

## Methods

No permit was required to sample invertebrates in lamprey nests, because the sampling was performed after the larvae left the nests, thereby preventing disturbance to lamprey lifecycle.
## Study site and experimental period

The study took place in the Nive River, a 79 km long river situated in Northern Basque Country, France, and draining a basin of 1030 km² (river width: 18 m; bankfull width: 25 m; mean river depth: 52 cm; average discharge: 21 m³/s). The selected section corresponds to a reach of 400 m long, located in Saint-Martin-d'Arrossa (43° 14' 34.926" N, 1° 18' 27.305" W) and bypassed by a hydropower scheme. It is mainly composed of riffles and runs suitable for sea lamprey nesting, as shown by the 46 nests found on a previous survey [26]. To avoid disturbing sea lamprey spawning, macroinvertebrate samples were collected on July 7 2020, at the end of the reproductive season, after two weeks of not detecting any lamprey near the nests. We considered these two weeks a period sufficient to allow for stability of local invertebrate assemblages [27], whereas reducing the risk of scouring by floods, which are frequent at the study site during late spring and early summer.

## Macroinvertebrate samples and nest characteristics

We studied 30 completed nests where no lamprey was observed (Fig 2). For each nest, we defined three zones. The **upstream zone**, 20 to 50 cm upstream from the upstream verge of the nest pit, indicated the initial conditions of the habitat chosen by the lamprey to build their nest, being both close to the nest and less likely to be disturbed by digging than the left, right

<figure><img src="image_6.png" /><figcaption>Fig 2. Location of the nests sampled in the selected section of the Niver river in Saint-Martin-d'Arrossa, bypassed by a hydropower scheme.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0274719.g002
or immediate downstream of the nest. The second zone was the **pit**, delimited as the area excavated by the lamprey, where water was deeper and the substrate finer than upstream. The third zone was the **mound**, corresponding to the downstream accumulation of pebbles and cobbles. For all three zones current velocity (± 1 cm/s) was measured at 5 cm from the bottom (using a magnetic flowmeter “FLO-MATE 2000”) and water depth (± 1 cm) was measured with a Vernier gauge. Transversal (perpendicular to the current) and longitudinal (parallel to the current) diameters of the pit were measured using a measuring tape (± 1 cm).

Macroinvertebrates were sampled using a Surber net (1/20 m², 0,5 mm mesh) by stirring the entire top layer of the sediment covered by the net frame (no more than 50 mm). Three samples, one per zone, were collected on each nest. After collection the samples were stored in 70° ethanol up to the determination. Most invertebrates were identified following AFNOR standard XP T 90-388 at the genus level except for Oligochaeta (subclass), Diptera (family) and some individuals for whom genus determination was not certain (determined at the family level). Macroinvertebrate traits were determined following Tachet et al. [28] and expressed as the frequency of individuals belonging to a certain trait over the total number of individuals in the sample. We selected food type, mode of alimentation, substrate preference and body size as these traits are likely to vary according to the physical characteristics of the nests. Food type and mode of alimentation directly reflect the resources provided by a given environment. With different current and substrate characteristics those traits should vary between zones. As substrate will differ, substrate preferences may reflect the sorting made by spawners between zones, if the scale of the nest has an impact on macroinvertebrate spatial location. Finally, body size emphasizes the carrying capacity for smaller or bigger invertebrates, which completes the substrate preference.

## Invertebrate diversity

All analyses were performed using R software version 4.1.0 [29] and a significance level of 0.05.

We used several α and β indices to measure the effects of lamprey nests on the diversity of macroinvertebrate assemblages. For α indices we calculated the indices in each zone to compare the pit and the mound with the upstream zone. First, we measured the species richness with the *specpool* function of the *vegan* package [30, version 2.5-7] using a Chao1 equation [31, Eq 1]. Then, Shannon α diversity index [32, Eq 2], calculated using the *diversity* function of the *vegan* package [30, version 2.5-7], indicated how diverse the taxa in each zone are, increasing with the increase of evenness and taxa richness in a sample. Shannon index was completed with Pielou index [33, Eq 3] to specifically analyze the evenness of the assemblages. These indices were transformed into Log Response Ratio (LRR), corresponding to the effect size of the diversity differences between nest (pit/mound) zones and upstream zones (Eq 4 and Fig 6) where the “index modified zone” corresponds either to the index of a pit or to the index of a mound. Indices were calculated for each pit and each mound and compared to the mean of Shannon and Pielou indices obtained for the upstream zones (“index upstream zone”). When the LRR standard deviation does not cross the value “0”, it means that the effect is significantly positive (above 0) or negative (below 0; [34]). Behind this analysis, we compared how a given index (Shannon or Pielou) differed on a pit or on a mound—thus reflecting a disturbed habitat—compared to the average value of these indices in the section considered.

$$Chao 1 = S + \frac{a^2}{2b} \qquad (1)$$
Where S is the number of taxa in a sample, a the number of taxa solely represented by one individual in the sample, and b the number of taxa represented by exactly two individuals in the sample.

$$ Shannon\ index(H') = -\sum_{i=1}^{S} p_i \ln p_i \qquad (2) $$

Where $p_i$ is the proportion of a taxon i on the total number of individuals N and S is the total number of taxa in the sample.

$$ Pielou\ evenness\ index(J') = \frac{H'}{H'_{\max}} \qquad (3) $$

Where $H'$ is the Shannon diversity index and $H'_{\max}$ is the maximum possible value of $H'$ (if every species was equally likely): $H'_{\max} = -\sum_{i=1}^{S} \frac{1}{S} \ln \frac{1}{S}$ and S being the total number of taxa in the sample.

$$ Log\ Response\ Ratio = \ln \left( \frac{index\ modified\ zone}{index\ upstream\ zone} \right) \qquad (4) $$

Following this $\alpha$ diversity overview we studied several $\beta$ diversity aspects to compare nest and upstream. For each index we calculated all the pairwise differences between pit and mound and between two upstream zones. Firstly, we calculated Sorensen dissimilarity index [35, Eq 5] which includes both replacement and nestedness and corresponds to a global indicator of taxa diversity. Then we calculated Simpson $\beta$ diversity index to study taxa replacement [35, Eq 6]. The second component of the Sorensen index, the taxa nestedness, was obtained by computing the difference between Sorensen and Simpson indices [36]. Finally, to test the diversity considering taxa abundance and not solely their occurrence, we calculated the Morisita dissimilarity index [37, Eq 7]. As for the $\alpha$ diversity indices the $\beta$ indices were transformed into Log Response Ratio (Eq 4 and Fig 7) where the “index modified zone” corresponds to a pairwise index between a pit and a mound and the “index upstream zone” to a pairwise index between two upstream zones. This LRR indicates if the diversity difference between two locations of a nest (always a pit versus a mound) differs from the same comparison made between two upstream zones, thus reflecting a difference in taxa heterogeneity. To avoid a nest effect we did not compare the taxa of a pit versus the mound of the same nest and therefore we randomly selected the pit and the mound to be compared in each index value. Indeed, pit and mound of the same nest may not be independent, and this analysis aimed at describing the heterogeneity of communities between pits and mounds in general, without the influence of similarities existing between adjacent samples.

$$ Sorensen\ \beta\ diversity\ index = \frac{(b+c)}{(2a+b+c)} \qquad (5) $$

Where a is the number of common taxa between two sites, b the number of taxa exclusive to a site, and c the number of taxa exclusive to the other site.

$$ Simpson\ \beta\ diversity\ index = \frac{\min(b, c)}{\min(b, c) + a} \qquad (6) $$

Where $\min(b, c)$ is the lowest number of exclusive taxa between two sites, b being the number of taxa exclusive to a site, and c the number of taxa exclusive to the other site; a is the number of
common taxa between the two sites.

$$Morisita overlap index(C_D) = \frac{2 \sum_{i=1}^{S} x_i y_i}{(D_x + D_y)XY} \qquad (7)$$

Where $x_i$ is the number of times taxon i is represented in the total X from one site; $y_i$ is the number of times taxon i is represented in the total Y from another site; $D_x$ and $D_y$ are the Simpson's $\alpha$ diversity index for the x and y sites respectively:

$$D_i = \sum_{i=1}^{S} p_i^2 \qquad (8)$$

With S the total number of taxa in the sample and $p_i$ the proportion of a taxon i on the total number of individuals N.

## Data analysis

The significance of differences between nest zones for depth, current, density, diversity, Shannon index and Pielou index was determined using linear mixed model and the *lmer* function of the *lme4* package [38] with a random effect of nest identity.

For density and taxa diversity (Fig 4A and 4B), a pairwise Tukey test was used to check the differences between each zone rather than globally.

To test the effect of current and depth on abundance and taxa diversity, a General Linear Model was used with current and depth as covariates and following a Gaussian distribution.

A two-sample t-test was used to test the significance of Chao1 diversity index calculated for nest (pit + mound) and upstream zone.

Trait differences between zones were determined using binomial mixed models with upstream or mound as reference condition and nest identity as a random effect (function *glmer* of the *lme4* package; Bates et al. [38]).

The links between trait and nest characteristics were studied using Redundancy analysis (RDA), a multivariate approach of linear models applied in *vegan* package [30, version 2.5-7] with *rda* function and setting traits as response variables. Traits were implemented after Hellinger transformation to give low weight to traits with low counts and many zeros. Explanatory variables were the five variables describing nest characteristics: transversal and longitudinal diameter, depth, current and depth difference between pit and mound. The latter was calculated assuming nests with higher depth difference could be more heterogeneous with a more important volume and habitat complexity.

## Results

### Nests characteristics

Lamprey nests differed greatly in their dimensions, transversal diameter of the pit averaging 129.2 ± 37.2 cm (range, 70–210 cm) and longitudinal diameter 117.4 ± 45.7 cm (range, 50–260 cm, Table 1). The pit was the deepest zone (59.2 ± 9.9 cm) followed by the upstream (51.8 ± 10.1 cm) and the mound (38.3 ± 10.1 cm), and differences were statistically significant considering mixed models p-values (Table 2). The coefficient of variation for depth was similar for all three zones (Table 2). Current velocity ranged from 0 to 65 cm/s (Table 1) and varied among zones inversely to depth (Fig 3), although the coefficient of variation was larger. Again, differences were statistically significant considering mixed models results.
<table><caption>Table 1. Summary of the nests characteristics measured in the study; mean ± sd (cv).</caption><thead><tr><th></th><th>Pit</th><th>Mound</th><th>Upstream</th></tr></thead><tbody><tr><td>Transversal diameter (cm)</td><td>129.2 ± 37.2 (0.3)</td><td>-</td><td>-</td></tr><tr><td>Longitudinal diameter (cm)</td><td>117.4 ± 45.7 (0.4)</td><td>-</td><td>-</td></tr><tr><td>Depth (cm)</td><td>59.2 ± 9.9 (0.2)</td><td>38.3 ± 10.1 (0.3)</td><td>51.8 ± 10.1 (0.2)</td></tr><tr><td>Current (cm/s)</td><td>8.8 ± 7.8 (0.9)</td><td>43.8 ± 11.6 (0.3)</td><td>19.8 ± 11.1 (0.6)</td></tr><tr><td>Depth difference (cm)</td><td>20.9 ± 4.8 (0.2)</td><td>-</td><td>-</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0274719.t001

## Macroinvertebrate density and diversity

The density of macroinvertebrates (Fig 4A) was 2649 ± 1386 individuals per $m^2$ in pit, 3833 ± 1052 individuals per $m^2$ in mound and 3777 ± 1332 individuals per $m^2$ in upstream. It was significantly lower in pit than in upstream and mound (Tukey test: $Df = 58, t.ratio = 3.529, p-value = 0.0023$ and $Df = 58, t.ratio = -3.705, p-value = 0.0014$ respectively) but there were no differences between upstream and mound (Tukey test: $Df = 58, t.ratio = -0.175, p-value = 0.9832$). The taxa richness (Fig 4B) also did not vary significantly between upstream (23.5 ± 3.9 taxa) and mound (21.2 ± 4.5 taxa; Tukey test: $Df = 58, t.ratio = 2.293, p-value = 0.0647$) but was higher in these two zones than in pit (18.6 ± 3.9 taxa; Tukey test: $Df = 58, t.ratio = 5.066, p-value \le 0.0001$ for pit compared to upstream and Tukey test: $Df = 58, t.ratio = -2.772, p-value = 0.0201$ for pit compared to mound).

General Linear Models indicated that current affected neither the number of taxa nor the macroinvertebrate density ($tvalue = 0.391; Pr(> |t|) = 0.697$ and $tvalue = 1.470; Pr(> |t|) = 0.1452$ respectively). Depth influenced negatively the macroinvertebrate density ($tvalue = -2.026; Pr(> |t|) = 0.0458$) but not the number of taxa ($tvalue = -0.821; Pr(> |t|) = 0.414$).

The Venn diagram (Fig 5), describing exclusivity or sharing of taxa, indicates an important taxa overlap with 59% of the taxa found on all zones. The significant pattern of reduced taxa richness found at the nest scale was not found by summing all the different taxa identified in the thirty samples of each zone and visualizing it in the Venn diagram.

## Diversity

The estimate of Chao1 diversity index was 82 ± 14 taxa for nest (pit + mound) and 69 ± 8 taxa for the upstream zone. In spite of the slight overlap of the standard error of those estimates the species richness was higher in the nest than upstream, as supported by results of two-sample t-test ($Df = 44.739, t = 4.5017, p-value \le 0.0001$). Log response ratios of the $\alpha$ diversity indices (Fig 6) showed an overall trend of reduced diversity and equitability for both Shannon and Pielou indices in pit and mound, compared to the average values observed in the upstream zone.

<table><caption>Table 2. Mixed model results for analyses of differences between zones for nest characteristics, density and diversity of macroinvertebrates and $\alpha$ diversity indices, with nest identity as a random effect. With $*P \le 0.05, **P \le 0.01, ***P \le 0.001$.</caption><thead><tr><th>Variables</th><th>Df</th><th>&chi;<sup>2</sup></th><th>P</th></tr></thead><tbody><tr><td>Depth</td><td>2</td><td>503.5</td><td>&le; 0.001***</td></tr><tr><td>Current</td><td>2</td><td>278.41</td><td>&le; 0.001***</td></tr><tr><td>Density</td><td>2</td><td>17.475</td><td>&le; 0.001***</td></tr><tr><td>Diversity</td><td>2</td><td>25.736</td><td>&le; 0.001***</td></tr><tr><td>Shannon index (upstream/pit/mound)</td><td>2</td><td>12.739</td><td>0.0017 **</td></tr><tr><td>Pielou index (upstream/pit/mound)</td><td>2</td><td>6.6937</td><td>0.035 *</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0274719.t002
<figure><img src="image_7.png" /><figcaption>Fig 3. Relationship between depth and current velocity for each nest at each zone. Ellipses correspond to multivariate t-distribution.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0274719.g003

However, only the log response ratio of the Shannon index of the pit is strictly below 0 and indicates a significantly reduced diversity in pit compared to the upstream zone (-10.5 ± 10.2%). The log response ratios of β diversity indices (Fig 7) showed a higher overall β diversity between a pit and a mound than between two upstream zones. However the log response ratios were highly variable and none were significantly different from 0. Nevertheless the higher Morisita index ratio compared to other indices seems to indicate a differentiation made by taxa abundance rather than by replacement or nestedness.

## Trait analysis

The traits studied tended to be more similar in mound and upstream than between these two locations and the pit (Table 3 and Figs 8–10). Considering alimentation traits (Fig 8), the proportion of collectors and scrapers was lower in the mound than in the pit and the proportion of filterers was lower in pit than either in mound or upstream. Predators were more present in mound than in pit and upstream, whereas shredders were less abundant in mound than in pit and upstream. However, food traits, related to alimentation aspects, did not highlight significant differences. Substrate preference analysis (Fig 9) showed a similar pattern for litter and mud with higher proportion in pit than in mound. Preference for gravel and sand was more important in pit, followed by upstream and then mound. Finally, preference for slabs, blocks, stones and pebbles was highest in mound and lowest in pit. The smallest size category (Fig 10)
<figure><img src="image_8.png" /><figcaption>Fig 4. Density of macroinvertebrates (ind/m²) (A) and number of taxa per sample (B) for the three zones studied. With ns P > 0.05, *P ≤ 0.05, **P ≤ 0.01, ****P ≤ 0.0001.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0274719.g004

<figure><img src="image_9.png" /><figcaption>Fig 5. Venn diagram with the distribution of taxa among zones. The area of each part is proportional to the number of taxa indicated in absolute number and percentage of total taxa richness.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0274719.g005
<figure><img src="image_10.png" /><figcaption>**Fig 6.** Log response ratios of Shannon and Pielou α diversity indices for each pit and each mound compared to the mean of Shannon and Pielou indices obtained for the upstream zones. Dots correspond to the mean and error bars correspond to the standard deviation. When the LRR standard deviation does not cross the value "0", it means that the effect is significantly positive (above 0) or negative (below 0).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0274719.g006

was more represented inside the pit, followed by upstream and mound, whereas larger size categories followed the inverse pattern.

**Relationship between traits and nest characteristics**

Redundancy analysis (Fig 11) yielded a significant model ($F = 10.876, Pr(> F) = 0.001$), although the only significant term was the current ($F = 46.8357, Pr(> F) = 0.001$), depth and depth difference being marginally insignificant ($F = 2.4535, Pr(> F) = 0.077$ and $F = 2.6966, Pr(> F) = 0.077$ respectively). The ellipses indicated a separation of pit from mound and upstream structured by current and "< 5 mm", "scrapers", "sand" and "gravel" trait categories belonging to the pit ellipse. Upstream and mound were not clearly separated but "5–10 mm", "slabs, blocks, stones and pebbles" and "filterers" belong to the mound ellipse. Those results are consistent with the specific trait analyzes (Figs 8–10). Depth, if not the most discriminant variable, explained the separation of pit from other zones.

**Discussion**

Firstly, our results globally showed that habitat diversity and structural heterogeneity due to lamprey nesting activities modify several aspects of biodiversity. The clear distinction between
<figure><img src="image_11.png" /><figcaption>Fig 7. Log response ratios of *β* diversity indices of all mound and pit pairwise indices, each compared to a pairwise upstream index among all possible upstream pairwise comparisons (upstream zone randomly selected, not directly upstream from the pit and mound considered). Dots correspond to the mean and error bars correspond to the standard deviation. When the LRR standard deviation does not cross the value "0", it means that the effect is significantly positive (above 0) or negative (below 0).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0274719.g007

nest zones regarding depth and riverbed current velocity implied biological heterogeneity, a non trivial result as previously highlighted [11]. The pit had a lower density of invertebrates and a lower number of taxa than the other zones, which did not differ between them, although the number of taxa tended to be higher upstream. Log response ratio of *α* Shannon index showed a similar trend concerning the diversity inside the pit. Those results indicate that nest building reduced the local abundance and diversity of macroinvertebrates in the pit, at least
<table><caption>Table 3. Binomial mixed models results comparing frequency of selected traits between zones with nest identity as a random effect. With *P* ≤ 0.05, **P** ≤ 0.001.</caption><thead><tr><th></th><th>Zone</th><th>Df</th><th>Z value</th><th>P</th></tr></thead><tbody><tr><td rowspan="3">Food</td><td>Mound (reference: upstream)</td><td>2</td><td>-1.016</td><td>0.31</td></tr><tr><td>Pit (reference: upstream)</td><td>2</td><td>-0.614</td><td>0.539</td></tr><tr><td>Mound (reference: pit)</td><td>2</td><td>-0.401</td><td>0.688</td></tr><tr><td rowspan="3">Alimentation</td><td>Mound (reference: upstream)</td><td>2</td><td>1.877</td><td>6.05e-2</td></tr><tr><td>Pit (reference: upstream)</td><td>2</td><td>-2.201</td><td>2.77e-2 *</td></tr><tr><td>Mound (reference: pit)</td><td>2</td><td>4.063</td><td>4.84e-5 ***</td></tr><tr><td rowspan="3">Substrate</td><td>Mound (reference: upstream)</td><td>2</td><td>-6.918</td><td>4.57e-12 ***</td></tr><tr><td>Pit (reference: upstream)</td><td>2</td><td>10.779</td><td>&lt;2e-16 ***</td></tr><tr><td>Mound (reference: pit)</td><td>2</td><td>-17.391</td><td>&lt;2e-16 ***</td></tr><tr><td rowspan="3">Size</td><td>Mound (reference: upstream)</td><td>2</td><td>5.682</td><td>1.33e-08 ***</td></tr><tr><td>Pit (reference: upstream)</td><td>2</td><td>-11.63</td><td>&lt;2e-16 ***</td></tr><tr><td>Mound (reference: pit)</td><td>2</td><td>16.71</td><td>&lt;2e-16 ***</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0274719.t003

during the first weeks after nest digging, which are the most significant from the point of view of the lamprey, as larvae abandon the nest shortly after hatching [39]. Patches of fine substrate such as the pit tend to be less diverse than zones with coarser grain size [40]. Reduced invertebrate abundance and richness have been also reported in nests of the largemouth bass (*Micropterus salmoides*) [41] and the pink salmon (*Oncorhynchus gorbuscha*) [42] during the spawning season, which was attributed to spawning-related disturbance. On the other hand, Hogg et al. [24] reported that the density of invertebrates in sea lamprey nest mounds were twice that found in pits and 77% higher than in upstream zone. Their results of reduced macroinvertebrate pit abundance are consistent with ours but we did not observe an increased mound density. An explanation for the important density in mounds could be the dominance of the Chironomidae within the river studied, not found for the Nive River (S1 Fig). The difference between our results and theirs seems not to derive from differences in macroinvertebrate taxa composition, as it was relatively similar between both studies (with Chironomidae, Hydropsychidae, Heptageniidae and Ephemerellidae dominant in [24] and abundant in the Nive River: S1 Fig and S1 Table. The significantly higher number of taxa in nest (mound and pit combined) than upstream, highlighted by the Chao1 index, indicates that sea lamprey nests create heterogeneity but also increase local species diversity.

In a local scale, i.e. the zone, we found for each pit and mound several taxa not present upstream (see S2 Fig). Indeed, whereas upstream is more complex and offers more substrate types with its unsorted grain size (due to the absence of sorting by lamprey during the nest building compared to the nest)—likely to provide suitable habitat for more taxa— [43], mound and pit provide larger areas of some specific characteristics, including substrate. This wider area may potentially allow the establishment of more taxa sharing similar ecological preferences, when these would be less represented and therefore potentially absent or very rare upstream. Considering the nest scale, rare taxa should be found more easily, either in mound or in pit considering their preferences than upstream (considering a same sample size).

In addition to increasing the structural and biological heterogeneity, sea lamprey nests seem to shape invertebrate assemblages. Indeed, three traits over the four studied showed differences between nest and the upstream zones, although the direction of change did not always follow our expectations. In particular, we expected herbivores to increase in mounds, as epilithic algae tend to be favored by shallow water, coarse sediments and moderately fast flow [44, 45], but no significant difference was found globally for alimentation traits, perhaps because of
<figure><img src="image_12.png" /><figcaption>Fig 8. Median percentage of macroinvertebrate alimentation traits for the three zones studied.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0274719.g008

lack of time for algae to grow on overturned cobbles. The increased proportion of collectors in pits was expected, as these are a preferential place for organic matter deposits. Similarly, collectors tend to be more abundant in fine substrate, which collected more detritus [46]. The higher concentration of organic matter might also explain the higher abundance of shredders and scrapers in pit than in mound. Also, according to our expectations, filterers were more represented in the mound, as they are favoured by fast-flowing areas [47], particularly during low flow conditions [48]. Those results reflect the different dynamics of food provided within the nest zones.
<figure><img src="image_13.png" /><figcaption>Fig 9. Median percentage of macroinvertebrate substrate preferences for the three zones studied.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0274719.g009

Not surprisingly given the differences in substrate among zones, invertebrates also differed in their preferred substrate traits. Macroinvertebrates associated with coarse substrate were especially abundant in mounds, whereas taxa associated with fine substrate were more abundant in pits. It must be noted that we visually noticed an absence of sand in the substrate surface of the study reach, apart from lamprey nests and some marginal areas, which suggests that nesting lampreys create patches of habitat for lentophylic species in river stretches where these species could not dwell otherwise.

The macroinvertebrate size traits showed a greater proportion of small invertebrates (< 5 mm) within the pit, including *Esolus sp.*, the most abundant taxon in pit (S1 Fig) and
<figure><img src="image_14.png" /><figcaption>Fig 10. Median percentage of macroinvertebrate size traits for the three zones studied.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0274719.g010

*Psychomyia pusilla*. This result could be explained by the carrying capacity of fine substrate. A fine substrate has less shelters than a coarser one. Shelters of fine substrate are more suitable for small size range of macroinvertebrates, whereas bigger macroinvertebrates are more prone to find shelter in the coarser substrate of the mound or upstream. Bêche et al. [49] suggested a better exploitation of refugia for small sized macroinvertebrates. Data on the grain size distribution of nest zones may help to confirm this hypothesis. Furthermore, small body size could be a possible resistance trait to disturbances such as riverbed movement [50].

Our results demonstrate that sea lamprey creates physical heterogeneity, which then enhances biological heterogeneity both in assemblage composition and function. The
<figure><img src="image_15.png" /><figcaption>**Fig 11. Redundancy analysis between nest characteristics and macroinvertebrate traits.** Slabs = slabs, blocks, stones and pebbles.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0274719.g011

heterogeneity demonstrated in this study was created at the nest scale, but local-scale processes can have a major impact at higher scales such as the reach scale [51]. In the case of lamprey nests in the Nive River, pits average 1.15 m² and so the total modified streambed area averages 34.5 m², representing 4% of the river streambed in our studied zone. With an average reduction of 30% of macroinvertebrate density in pit (compared to upstream) a total reduction of 1.2% of density is expected in the studied zone. In a river reacting in a similar way to the Nive, this surface with decreased macroinvertebrate density does not seem to be compensated by a higher density in the mound and so, sea lamprey nests globally decrease the macroinvertebrate density. To complete this finding, samples collected on a substrate further from the nests than our upstream samples may help to verify than they were not negatively affected by the downstream digging activity of lampreys.

Measuring the effect of lamprey nests on reach-scale macroinvertebrate density would require a BACI design [52] including periods and reaches with and without nests. If invertebrates affected by lamprey nesting responded with small-scale movements, the result would be an increase in macroinvertebrate density in upstream samples of reaches with nests compared to upstream samples of reaches without nests. If the macroinvertebrates responded drifting, the results would be more complex to interpret and depend on the distribution of nests of sea lamprey. An increase of macroinvertebrate densities could be observed if spawning sites occur upstream in the river. Macroinvertebrates drifting from those nests could colonize and increase the density of “upstream” sites situated downstream. When a spawning ground is located on the upper limit of the area colonized by lampreys, drift will not affect the adjacent area of this spawning ground but may affect downstream sites.

Our results highlight the need for a better understanding of the effects of nest-building species on the river ecosystem. Unlike salmon, whose effects are relatively well-studied [53–55], the consequences of sea lamprey spawning are still poorly known. Whereas lamprey nests only represent 4% of the river streambed in our study, the occupied area may be much higher, especially for invasive populations of the Great Lakes [56] or below impassable dams, likely to
dramatically increase the number of spawners in immediately downstream spawning grounds [26, 57]. Studies exist on the effects of carcasses [25, 58] or substrate modification [24] on a local scale, but not at the river or at the global scales. However, macroinvertebrates have a basal to intermediate position on the food chain [59]. Such position implies an important relationship with ecosystem components. Firstly, macroinvertebrate assemblages can have an important influence on numerous processes [59], such as nutrient cycles [60–62], primary productivity [63], decomposition [64–66], and translocation of materials [67, 68]. Then, they are an important resource for organisms belonging to higher trophic levels such as fish species. Macroinvertebrates are considered as the most important source of food, widespread in all freshwater ecosystems [69].

In our study some trait aspects can provide clues about the effects of sea lamprey on general ecosystem functioning in the river. Increased proportion of shredders in the pit suggests they find suitable conditions such as increased plant debris. On a more global scale a certain amount of material (depending on the nest density and count) could be retained by nests in spite of being carried away by the current, and then consumed by macroinvertebrates. Seeing this impact in late spring when litter input is reduced also indicates that the proportion of shredders should be more important in autumn, where these inputs can represent 73% of annual allochtonous inputs (in temperate deciduous forests; Abelho and Graça [70]). But this retention can occur only if the lamprey nests remain in autumn, without being filled by floods. Changes in streambed complexity created by nests persist during the autumn [24] but a change in hydrological conditions is likely to affect this process. Filterers are another trait category possibly affecting organic matter processing, being one pathway through which carbon and nutrients are transferred from the water column to the sediments [71]. Filterers being more present in proportion in the mound than upstream suggests that mounds are hotspots of nutrient processing. Therefore, sea lamprey nests could have an important effect on nutrient spiraling in streams, but the persistence of their influence on ecosystem functioning until the next spawning season remains to be assessed. Although nests may persist considering their physical structure, it is likely that the subsidies provided by eggs affect macroinvertebrate colonization dynamics and the distribution of functional groups, as predation was depicted in other species [72, 73]. An assessment of the effects of eggs deposited in the pit on functionality should be done to determine what is their potential contribution to the observed differences.

All these possible different effects imply that the modifications of macroinvertebrate assemblages by sea lampreys highlighted here and in previous studies [24, 25] are likely to modify the general functioning of the rivers. However, it is difficult to determine these consequences precisely due to the complexity inside macroinvertebrate assemblages. Such assessment requires specific studies on the components previously described. We advocate for studies focusing on these effects, important in a context of endangered native populations, likely to disappear in some rivers and consequently modifying their dynamics.

**Supporting information**

S1 Fig. Median of the abundances for the 10 most abundant taxa within each zone (Upstream, Pit, Mound).
(TIFF)

S2 Fig. Number of exclusive taxa in mound or pit compared to the upstream. With
**P** ≤ 0.01.
(TIFF)
**S1 Table. List of taxa found in each zone of each nest.** With, for column names: 1 Leuctra; 2 Perlodes; 3 Perla; 4 Rhyacophila; 5 Glossosomatidae; 6 Agapetus; 7 Hydroptila; 8 Hydropsychidae; 9 Hydropsyche; 10 Polycentropodidae; 11 Polycentropus; 12 Psychomyiidae; 13 Psychomyia; 14 Oligoplectrum; 15 Micrasema; 16 Goeridae; 17 Lepidostoma; 18 Athripsodes; 19 Ceraclea; 20 Sericostomatidae; 21 Baetis; 22 Oligoneuriella; 23 Heptageniidae; 24 Epeorus; 25 Rhithrogena; 26 Ecdyonurus; 27 Ephemerella ignita; 28 Caenis; 29 Ephemera; 30 Potamanthus; 31 Hydraena; 32 Stenelmis; 33 Elmis; 34 Esolus; 35 Oulimnius; 36 Limnius; 37 Normandia; 38 Micronecta; 39 Aphelocheirus; 40 Blephariceridae; 41 Limoniidae; 42 Simuliidae; 43 Tanypodinae; 44 Ceratopogonidae; 45 Empididae; 46 Athericidae; 47 Gammaridae; 48 Echinogammarus; 49 Hydracarina; 50 Piscicolidae; 51 Erpobdellidae; 52 Oligochaeta; 53 Theodoxus; 54 Potamopyrgus; 55 Radix; 56 Planorbidae; 57 Ancylus; 58 Dugesiidae; 59 Nematoda; 60 Hydroporinae; 61 Hydrophilinae; 62 Prostoma; 63 Hydrozoa; 64 Ostracoda; 65 Copepoda; 66 Agapetinae; 67 Physella; 68 Leuctra geniculata; 69 Chironomidae excl. Tanypodinae.
(PDF)

**Acknowledgments**

Field work used resources from the IE ECP Experimental Facility of the UMR Ecobiop [74].

**Author Contributions**

**Conceptualization:** Marius Dhamelincourt, Cédric Tentelier, Arturo Elosegi.

**Data curation:** Marius Dhamelincourt, Jacques Rives, Marie Pons, Cédric Tentelier, Arturo Elosegi.

**Formal analysis:** Marius Dhamelincourt, Marie Pons, Aitor Larrañaga.

**Funding acquisition:** Cédric Tentelier.

**Investigation:** Marius Dhamelincourt, Jacques Rives, Cédric Tentelier, Arturo Elosegi.

**Methodology:** Marius Dhamelincourt, Jacques Rives, Cédric Tentelier, Arturo Elosegi.

**Supervision:** Cédric Tentelier, Arturo Elosegi.

**Validation:** Cédric Tentelier, Arturo Elosegi.

**Writing – original draft:** Marius Dhamelincourt.

**Writing – review & editing:** Marius Dhamelincourt, Aitor Larrañaga, Cédric Tentelier, Arturo Elosegi.

**References**

1. Michael J. Cramer and Michael R. Willig. Habitat heterogeneity, species diversity and null models. *Oikos*, 108(2):209-218, 2005. ISSN 1600-0706. https://doi.org/10.1111/j.0030-1299.2005.12944.x

2. Tews J., Brose U., Grimm V., Tielbörger K., Wichmann M. C., Schwager M., et al. Animal species diversity driven by habitat heterogeneity/diversity: The importance of keystone structures. *Journal of Biogeography*, 31(1):79-92, 2004. ISSN 1365-2699. https://doi.org/10.1046/j.0305-0270.2003.00994.x

3. MacArthur Robert H. and MacArthur John W. On Bird Species Diversity. *Ecology*, 42(3):594-598, 1961. ISSN 1939-9170. https://doi.org/10.2307/1932254

4. Li H. and Reynolds J. F. On Definition and Quantification of Heterogeneity. *Oikos*, 73(2):280-284, 1995. ISSN 0030-1299. https://doi.org/10.2307/3545921

5. González-Megías Adela, María Gómez José, and Sánchez-Piñero Francisco. Diversity-habitat heterogeneity relationship at different spatial and temporal scales. *Ecography*, 30(1):31-41, 2007. ISSN 1600-0587. https://doi.org/10.1111/j.0906-7590.2007.04867.x
6. Elosegi Arturo, Flores Lorea, and Díez Joserra. The importance of local processes on river habitat characteristics: a Basque stream case study. *Limnética*, 30(2):0183-196, 2011. ISSN 0213-8409.

7. Flores Lorea, Ramón Díez Jose, Larrañaga Aitor, Pascoal Cláudia, and Elosegi Arturo. Effects of retention site on breakdown of organic matter in a mountain stream. *Freshwater Biology*, 58(6):1267-1278, 2013. ISSN 1365-2427. https://doi.org/10.1111/fwb.12125

8. Beisel Jean-Nicolas, Usseglio-Polatera Philippe, and Moreteau Jean-Claude. The Spatial Heterogeneity of a River Bottom: A Key Factor Determining Macroinvertebrate Communities. *Hydrobiologia*, 422-423:163-171, April 2000. https://doi.org/10.1023/A:1017094606335

9. Boyero Luz. The quantification of local substrate heterogeneity in streams and its significance for macroinvertebrate assemblages. *Hydrobiologia*, 499(1):161-168, June 2003. ISSN 1573-5117. https://doi.org/10.1023/A:1026321331092

10. Beisel Jean-Nicolas, Usseglio-Polatera Philippe, Thomas Sandra, and Moreteau Jean-Claude. Stream community structure in relation to spatial variation: The influence of mesohabitat characteristics. *Hydro-biologia*, 389(1):73-88, December 1998. ISSN 1573-5117. https://doi.org/10.1023/A:1003519429979

11. Palmer Margaret A., Menninger Holly L., and Bernhardt Emily. River restoration, habitat heterogeneity and biodiversity: A failure of theory or practice? *Freshwater Biology*, 55(s1):205-222, 2010. ISSN 1365-2427. https://doi.org/10.1111/j.1365-2427.2009.02372.x

12. Jones Clive G., Lawton John H., and Shachak Moshe. Organisms as Ecosystem Engineers. *Oikos*, 69 (3):373-386, 1994. ISSN 0030-1299. https://doi.org/10.2307/3545850

13. Butler David R. *Zoogeomorphology: Animals as Geomorphic Agents*. Cambridge University Press, Cambridge, 1995. ISBN 978-0-521-43343-3.

14. Wright Justin P., Jones Clive G., and Flecker Alexander S. An ecosystem engineer, the beaver, increases species richness at the landscape scale. *Oecologia*, 132(1):96-101, June 2002. ISSN 0029-8549, 1432-1939. https://doi.org/10.1007/s00442-002-0929-1

15. Pledger A. G., Rice S. P., and Millett J. Bed disturbance via foraging fish increases bedload transport during subsequent high flows and is controlled by fish size and species. *Geomorphology*, 253:83-93., January 2016. ISSN 0169-555X.

16. Gottesfeld Allen, Hassan Marwan, and Tunnicliffe Jon. Salmon Bioturbation and Stream Process. *American Fisheries Society Symposium*, 65, January 2008.

17. Moore J. W. and Schindler D. E. Impacts of Bioturbation by Spawning Salmon on the Community Dynamics and Ecosystem Processes of Alaskan Streams. *AGU Spring Meeting Abstracts*, 21:NB21B-06, May 2005.

18. Moore Jonathan W. and Schindler Daniel E. Biotic disturbance and benthic community dynamics in salmon-bearing streams. *Journal of Animal Ecology*, 77(2):275-284, 2008. ISSN 1365-2656. https://doi.org/10.1111/j.1365-2656.2007.01336.x PMID: 18081781

19. Swartwout Meredith C., Keating Fiona, and Frimpong Emmanuel A. A survey of macroinvertebrates colonizing bluehead chub nests in a Virginia stream. *Journal of Freshwater Ecology*, 31(1):147-152, January 2016. ISSN 0270-5060. https://doi.org/10.1080/02705060.2015.1036943

20. Gardner Cory, Coghlan Stephen M. Jr, and Zydlewski Joseph. Distribution and Abundance of Anadromous Sea Lamprey Spawners in a Fragmented Stream: Current Status and Potential Range Expansion Following Barrier Removal. *Northeastern Naturalist*, 19(1):99-110, March 2012. ISSN 1092-6194, 1938-5307. https://doi.org/10.1656/045.019.0108

21. Sousa R., Araújo M. J., and Antunes C. Habitat modifications by sea lampreys (*Petromyzon marinus*) during the spawning season: Effects on sediments. *Journal of Applied Ichthyology*, 28(5):766-771, October 2012. ISSN 1439-0426. https://doi.org/10.1111/j.1439-0426.2012.02025.x

22. Marius Dhamelincourt. Individual and group characteristics affecting nest building in sea lamprey, 2021b.

23. Patrick Grellier. La biologie de la lamproie marine (*Petromyzon marinus* Linne 1758) de la côte atlantique française. *Archimer*, page 18, 1996.

24. Hogg Robert S., Coghlan Stephen M., Zydlewski Joseph, and Simon Kevin S. Anadromous sea lampreys (*Petromyzon marinus*) are ecosystem engineers in a spawning tributary. *Freshwater Biology*, 59 (6):1294-1307, 2014. ISSN 1365-2427. https://doi.org/10.1111/fwb.12349

25. Weaver Daniel M., Coghlan Stephen M. Jr, and Zydlewski Joseph. Effects of sea lamprey substrate modification and carcass nutrients on macroinvertebrate assemblages in a small Atlantic coastal stream. *Journal of Freshwater Ecology*, 33(1):19-30, January 2018. ISSN 0270-5060. https://doi.org/10.1080/02705060.2017.1417168

26. Migradour. Suivi de la reproduction de la Lamproie marine sur le bassin de l'Adour—Tranche 1/3, gaves et nives. Technical report, 2010.
27. Khalaf G. and Tachet H. La dynamique de colonisation des substrats artificiels par les macroinvertebres d'un cours d'eau. *Annales de Limnologie*, 13(2):169-190, 1977. ISSN 0003-4088, 2100-000X. https://doi.org/10.1002/9781118445112.stat03432.pub2

28. Henri Tachet, Philippe Richoux, Michel Bournaud, and Philippe Usseglio-Polatera. *Invertébrés d'eau Douce*. CNRS éditions édition, May 2000. ISBN 978-2-271-05745-7.

29. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, 2020.

30. Jari Oksanen, Roeland Kindt, Pierre Legendre, Bob Hara, Gavin Simpson, Peter Solymos, et al. The vegan Package. January 2009.

31. Anne Chao and Chun-Huo Chiu. Species Richness: Estimation and Comparison. In *Wiley StatsRef: Statistics Reference Online*, pages 1-26. American Cancer Society, 2016. ISBN 978-1-118-44511-2.

32. Shannon C. E. A Mathematical Theory of Communication. *Bell System Technical Journal*, 27(3):379-423, 1948. ISSN 1538-7305. https://doi.org/10.1002/j.1538-7305.1948.tb01338.x

33. Pielou E. C. *Ecological Diversity*. Wiley, New York, 1975. ISBN 978-0-471-68925-6.

34. Fazlioglu Fatih and Chen Luzhen. Introduced non-native mangroves express better growth performance than co-occurring native mangroves. *Scientific Reports*, 10:1-11, March 2020. https://doi.org/10.1038/s41598-020-60454-z PMID: 32123225

35. Koleff Patricia, Gaston Kevin J., and Lennon Jack J. Measuring beta diversity for presence-absence data. *Journal of Animal Ecology*, 72(3):367-382, 2003. ISSN 1365-2656. https://doi.org/10.1046/j.1365-2656.2003.00710.x

36. Baselga Andrés. Partitioning the turnover and nestedness components of beta diversity. *Global Ecology and Biogeography*, 19(1):134-143, 2010. ISSN 1466-8238. https://doi.org/10.1111/j.1466-8238.2009.00490.x

37. Wolda Henk. Similarity indices, sample size and diversity. *Oecologia*, 50(3):296-302, September 1981. ISSN 1432-1939. https://doi.org/10.1007/BF00344966 PMID: 28309044

38. Bates Douglas, Mächler Martin, Bolker Ben, and Walker Steve. Fitting linear mixed-effects models using lme4. *Journal of Statistical Software*, 67(1):1-48, October 2015. ISSN 1548-7660. https://doi.org/10.18637/jss.v067.i01

39. Derosier Amy L., Jones Michael L., and Scribner Kim T. Dispersal of sea lamprey larvae during early life: Relevance for recruitment dynamics. *Environmental Biology of Fishes*, 78(3):271-284, March 2007. ISSN 1573-5133. https://doi.org/10.1007/s10641-006-9095-3

40. Duan Xuehua, Wang Zhaoyin, Xu Mengzhen, and Zhang Kang. Effect of streambed sediment on benthic ecology. *International Journal of Sediment Research*, 24(3):325-338, September 2009. ISSN 1001-6279. https://doi.org/10.1016/S1001-6279(10)60007-8

41. Kupriyanova E. K. and Bailey C. F. The Disturbance Effect of Largemouth Bass Nesting on a Benthic Macroinvertebrate Community. *Journal of Freshwater Ecology*, 13(3):333-341, September 1998. ISSN 0270-5060. https://doi.org/10.1080/02705060.1998.9663626

42. Campbell Emily Y., Merritt Richard W., Cummins Kenneth W., and Eric Benbow M. Spatial and Temporal Variability of Macroinvertebrates in Spawning and Non-Spawning Habitats during a Salmon Run in Southeast Alaska. *PLOS ONE*, 7(6):e39254, June 2012. ISSN 1932-6203. https://doi.org/10.1371/journal.pone.0039254 PMID: 22745724

43. Barnes Jacqueline B., Vaughan Ian P., and Ormerod Steve J. Reappraising the effects of habitat structure on river macroinvertebrates. *Freshwater Biology*, 58(10):2154-2167, 2013. ISSN 1365-2427. https://doi.org/10.1111/fwb.12198

44. Lau Y. L. and Liu D. Effect of flow rate on biofilm accumulation in open channels. *Water Research*, 27(3):355-360, March 1993. ISSN 0043-1354. https://doi.org/10.1016/0043-1354(93)90034-F

45. Horner R. R. and Welch E. B. Stream Periphyton Development in Relation to Current Velocity and Nutrients. *Canadian Journal of Fisheries and Aquatic Sciences*, April 2011. https://doi.org/10.1139/f81-062

46. Parker Michael S. Effect of Substrate Composition on Detritus Accumulation and Macroinvertebrate Distribution in a Southern Nevada Desert Stream. *The Southwestern Naturalist*, 34(2):181-187, 1989. ISSN 0038-4909. https://doi.org/10.2307/3671727

47. Habdija Ivan, Primc Habdija Biserka, Matoničkin Renata, Kučinić Mladen, Radanović Ines, Miliša Marko, et al. Current velocity and food supply as factors affecting the composition of macroinvertebrates in bryophyte habitats in karst running water. *Biologia (Bratislava)*, 59(5):577, 2004.

48. Brooks Andrew J. and Haeusler Tim. Invertebrate responses to flow: Trait-velocity relationships during low and moderate flows. *Hydrobiologia*, 773(1):23-34, June 2016. ISSN 0018-8158, 1573-5117. https://doi.org/10.1007/s10750-016-2676-z
49. Bêche Leah A., Mcelravy Eric P., and Resh Vincent H. Long-term seasonal variation in the biological traits of benthic-macroinvertebrates in two Mediterranean-climate streams in California, U.S.A. *Freshwater Biology*, 51(1):56–75, 2006. ISSN 1365-2427. https://doi.org/10.1111/j.1365-2427.2005.01473.x

50. Townsend Colin, Dolédec Sylvain, and Scarsbrook Mike. Species traits in relation to temporal and spatial heterogeneity in streams: A test of habitat templet theory. *Freshwater Biology*, 37(2):367–387, 1997. ISSN 1365-2427. https://doi.org/10.1046/j.1365-2427.1997.00166.x

51. Thrush Simon F., Hewitt Judi E., Cummings Vonda J., Green Malcolm O., Funnell Greig A., and Wilkinson Michelle R. The Generality of Field Experiments: Interactions Between Local and Broad-Scale Processes. *Ecology*, 81(2):399–415, 2000. ISSN 1939-9170. https://doi.org/10.1890/0012-9658(2000)081%5B0399:TGOFEI%5D2.0.CO;2

52. Underwood A. J. Beyond BACI: The detection of environmental impacts on populations in the real, but variable, world. *Journal of Experimental Marine Biology and Ecology*, 161(2):145–178, October 1992. ISSN 0022-0981. https://doi.org/10.1016/0022-0981(76)90138-6

53. Bilby Robert E, Fransen Brian R, and Bisson Peter A. Incorporation of nitrogen and carbon from spawning coho salmon into the trophic system of small streams: Evidence from stable isotopes. *Canadian Journal of Fisheries and Aquatic Sciences*, 53(1):164–173, January 1996. ISSN 0706-652X. https://doi.org/10.1139/f95-159

54. Janetski David J., Chaloner Dominic T., Tiegs Scott D., and Lamberti Gary A. Pacific salmon effects on stream ecosystems: A quantitative synthesis. *Oecologia*, 159(3):583–595, March 2009. ISSN 1432-1939. https://doi.org/10.1007/s00442-008-1249-x

55. Verspoor Jan J., Braun Douglas C., Stubbs Morgan M., and Reynolds John D. Persistent ecological effects of a salmon-derived nutrient pulse on stream invertebrate communities. *Ecosphere*, 2(2):art18, 2011. ISSN 2150-8925. https://doi.org/10.1890/ES10-00011.1

56. Bernard R. Smith. Sea lampreys in the Great Lakes of North America. pages 207–247, 1971.

57. Vernon C. Applegate. Natural history of the sea lamprey, *Petromyzon marinus*, in Michigan. Federal Government Series 55, U.S. Fish and Wildlife Service, Ann Arbor, Michigan: University of Michigan Library, 1950.

58. Weaver Daniel, Coghlan Stephen M., and Zydlewski Joseph D. Sea Lamprey Carcasses Exert Local and Variable Food Web Effects in a Nutrient-limited Atlantic Coastal Stream. *Canadian Journal of Fisheries and Aquatic Sciences*, 73, May 2016. https://doi.org/10.1139/cjfas-2015-0506

59. Bruce Wallace J. and Webster Jackson R. The Role of Macroinvertebrates in Stream Ecosystem Function. *Annual Review of Entomology*, 41(1):115–139, 1996. https://doi.org/10.1146/annurev.en.41.010196.000555

60. Bruce Wallace J. and Gurtz Martin E. Response of Baetis Mayflies (Ephemeroptera) to Catchment Logging. *The American Midland Naturalist*, 115(1):25–41, 1986. ISSN 0003-0031. https://doi.org/10.2307/2425834

61. Grimm Nancy B. Role of Macroinvertebrates in Nitrogen Dynamics of a Desert Stream. *Ecology*, 69(6):1884–1893, 1988. ISSN 1939-9170. https://doi.org/10.2307/1941165

62. Hunter Mark D., Adl Sina, Pringle Catherine M., and Coleman David C. Relative effects of macroinvertebrates and habitat on the chemistry of litter during decomposition. *Pedobiologia*, 47(2):101–115, January 2003. ISSN 0031-4056. https://doi.org/10.1078/0031-4056-00174

63. Hill Walter R., Boston Harry L., and Steinman Alan D. Grazers and Nutrients Simultaneously Limit Lotic Primary Productivity. *Canadian Journal of Fisheries and Aquatic Sciences*, 49(3):504–512, March 1992. ISSN 0706-652X. https://doi.org/10.1139/f92-059

64. Ferreira Verónica, Chauvet Eric, and Canhoto Cristina. Effects of experimental warming, litter species, and presence of macroinvertebrates on litter decomposition and associated decomposers in a temperate mountain stream. *Canadian Journal of Fisheries and Aquatic Sciences*, 72(2):206–216, February 2015. ISSN 0706-652X. https://doi.org/10.1139/cjfas-2014-0119

65. Mutshekwa Thendo, Cuthbert Ross N., Wasserman Ryan J., Murungweni Florence M., and Dalu Tatenda. Macroinvertebrate colonisation associated with native and invasive leaf litter decomposition. *Knowledge & Management of Aquatic Ecosystems*, (421): 32, 2020. ISSN 1961-9502. https://doi.org/10.1051/kmae/2020025

66. Anderson N H and Sedell J R. Detritus Processing by Macroinvertebrates in Stream Ecosystems. *Annual Review of Entomology*, 24(1):351–377, 1979. https://doi.org/10.1146/annurev.en.24.010179.002031

67. Windsor Fredric M., Tilley Rosie M., Tyler Charles R., and Ormerod Steve J. Microplastic ingestion by riverine macroinvertebrates. *Science of The Total Environment*, 646:68–74, January 2019. ISSN 0048-9697. https://doi.org/10.1016/j.scitotenv.2018.07.271 PMID: 30048870
68. Mermillod-Blondin Florian, Gaudet Jean-Paul, Gérino Magali, Desrosiers Gaston, and Creuzé des Châ-telliers Michel. Influence of macroinvertebrates on physico-chemical and microbial processes in hypor-heic sediments. *Hydrological Processes*, 17(4):779–794, 2003. ISSN 1099-1085. https://doi.org/10.1002/hyp.1165

69. Cummins Kenneth W. Hynes H. B. N. The ecology of running waters. *Limnology and Oceanography*, 16(3):593–594, 1971. ISSN 1939-5590. https://doi.org/10.4319/lo.1971.16.3.0593

70. Abelho Manuela and Graça M. A. S. Litter in a first-order stream of a temperate deciduous forest (Margarça Forest, central Portugal). *Hydrobiologia*, 386(1):147–152, November 1998. ISSN 1573-5117. https://doi.org/10.1023/A:1003532921432

71. Monaghan Michael, Thomas Steven, Minshall G, Newbold Denis, and Cushing Colbert. The influence of filter-feeding benthic macroinvertebrates on the transport and deposition of particulate organic matter and diatoms in two streams. *Limnology and oceanography, 46*:1091–1099, July 2001. https://doi.org/10.4319/lo.2001.46.5.1091

72. Brown Andrew F. and Diamond Mark. The consumption of rainbow trout (*Salmo gairdneri Richardson*) eggs by macroinvertebrates in the field. *Freshwater Biology*, 14(2):211–215, 1984. ISSN 1365-2427. https://doi.org/10.1111/j.1365-2427.1984.tb00036.x

73. Mills Chris A. Egg population dynamics of naturally spawning dace, *Leuciscus leuciscus* (L.). *Environmental Biology of Fishes*, 6(2):151–158, May 1981. ISSN 1573-5133. https://doi.org/10.1007/BF00002778

74. ECP. Ecology and Fish Population Biology Facility. INRAE, 2018.