Species-specific responses to landscape fragmentation: implications for management strategies

Habitat fragmentation affects the integrity of many species, but little is known about species-specific sensitivity to fragmentation. Here, we compared the genetic structure of four freshwater fish species differing in their body size (Leuciscus cephalus; Leuciscus leuciscus; Gobio gobio and Phoxinus phoxinus) between a fragmented and a continuous landscape. We tested if, overall, fragmentation affected the genetic structure of these fish species, and if these species differed in their sensitivity to fragmentation. Fragmentation negatively affected the genetic structure of these species. Indeed, irrespective of the species identity, allelic richness and heterozygosity were lower, and population divergence was higher in the fragmented than in the continuous landscape. This response to fragmentation was highly species-specific, with the smallest fish species (P. phoxinus) being slightly affected by fragmentation. On the contrary, fish species of intermediate body size (L. leuciscus and G. gobio) were highly affected, whereas the largest fish species (L. cephalus) was intermediately affected by fragmentation. We discuss the relative role of dispersal ability and effective population size on the responses to fragmentation we report here. The weirs studied here are of considerable historical importance. We therefore conclude that restoration programmes will need to consider both this societal context and the biological characteristics of the species sharing this ecosystem.


Introduction
Habitat fragmentation is probably the most pervasive effect humans impose on wild species (Vitousek et al. 1997;Lawler et al. 2006). It is also among the most highly studied phenomenon of the process of global change, and many studies have improved our knowledge on the ecological and evolutionary outcomes of fragmented populations (Fahrig 2003;Ewers and Didham 2006;Lawler et al. 2006). Typically, at the meta-population scale, habitat fragmentation causes habitat patches to be reduced in size and to be isolated from one another, hence decreasing gene flow between patches (Fahrig 2003). For most species, this spatial rearrangement also decreases the effective population size as well as genetic diversity at the patch scale, through the processes of genetic drift and inbreeding (Frankham 1998;Couvet 2002;DiBattista 2008). Ultimately, in the better cases, changes in the genetic structure of connected populations modify their evolutionary trajectory (by changes in life-history characteristics for instance Johansson et al. 2007;Waples et al. most other species) that are used to define the appropriate ways landscapes must be managed for preserving diversity (Lambeck 1997). Indeed, several phenotypic traits can make a species more or less susceptible to habitat fragmentation. Body size is one such trait because it correlates with several characteristics of individuals or populations (e.g. dispersal, effective population size or trophic status) that are known to affect species sensibility to fragmentation (Davies et al. 2000;Henle et al. 2004;Ewers and Didham 2006). For instance, all things being equal, theory predicts that large-bodied species should have a higher dispersal capability enabling them to rescue or re-colonize distant patches and hence to support patch isolation, even in a highly fragmented habitat (reviewed in Ewers and Didham 2006;Henle et al. 2004). Alternatively, large-bodied species often have a small population size (Cotgreave 1993;Blackburn and Gaston 1997), and rare species are expected to be more sensitive to fragmentation, probably because of their low effective population size (Ne) (Davies et al. 2000). However, empirical studies testing the hypothesis that body size underlies species-specific responses to fragmentation remain scant because most studies target single species (but see Davies et al. 2000).
Fragmentation caused by humans is particularly prevalent in freshwater ecosystems (Nilsson et al. 2005). Indeed, for a very long time, humans have built dams for managing irrigation, producing energy or for recreational purposes. Dams are so widespread that the water flow of over half the large river systems in the world is affected by them (Nilsson et al. 2005). The ecological and evolutionary impacts of dams are obvious and range from the loss of biodiversity to the malfunctioning of the whole ecosystem (Loot et al. 2007;Poulet 2007;Grenouillet et al. 2008;Maloney et al. 2008). At the genetic level, the negative effects of dams have been unambiguously demonstrated for many fish species (Alo and Turner 2005;Wofford et al. 2005;Raeymaekers et al. 2008). Based on these results, a straightforward restoration tool consists in 'un-damming' rivers (Bednarek 2001;Palmer et al. 2005Palmer et al. , 2008. Although environmentally beneficial, this solution can however prove undesirable when dams are also part of the historical and cultural heritage of a region or when they provide benefits for society. For example, many small to medium European rivers are scattered with water mills associated to small weirs (2-3 m high) that, at least partially, impede the movement of fish (Ovidio and Philippart 2002;Raeymaekers et al. 2008). Many of these man-made structures date back to the 15th century, and are thus part of the cultural heritage of local populations (see Fig. 1 and Raeymaekers et al. 2009). In such cases, adopting management practices requires consideration of both the ecological and social dimensions of these structures.
In this study, we present a comparative study aimed at evaluating the effects of anthropogenic river fragmentation on the genetic integrity of four widespread fish species (the European chub, Leuciscus cephalus; the rostrum dace, Leuciscus leuciscus; the gudgeon, Gobio gobio and the European minnow, Phoxinus phoxinus). These fish species differ mostly with respect to their maximal body size and we hypothesized that they may display species-specific responses to river fragmentation (Davies et al. 2000). To test this hypothesis, we sampled and genotyped (using microsatellite loci) these species along two nearby river basins (i.e. landscapes); one being highly fragmented (the river Viaur) and the other being a continuous landscape (the river Célé). We first compared the spatial patterns of genetic diversity between the landscapes and among the fish species using allelic richness (AR) and heterozygosity as two measures of genetic diversity. Irrespective of the species, we expected to find a lower genetic diversity in the fragmented landscape because of the stronger influence of patch isolation on genetic drift and inbreeding (Frankham 1998;Couvet 2002;Keller and Waller 2002). We also expected the smallest fish species to exhibit the greatest response to fragmentation if dispersal ability is the most important characteristic explaining species sensitivity to fragmentation (Ewers and Didham 2006). The reverse is expected (i.e. the largest species are more sensitive) if, alternatively, effective population size is the best explanation for sensitivity to fragmentation (Davies et al. 2000). By combining the use of F-statistics (Wright 1951) and the isolation-by-distance (IBD) framework (hereafter IBD Hutchison and Templeton 1999;Slatkin 1993), we then compared the population genetic structure of each species between the two landscapes. We expect a higher level of genetic structure (i.e. higher F st values) for all species in the fragmented landscape because weirs might partially impede dispersal between patches (Raeymaekers Figure 1 Traditional mill weir of the river Viaur (the fragmented landscape). They are made of rocks and most date back to the Middle Ages (12th-15th centuries). They are hence part of the local and cultural heritage, but they also impede the dispersal of several fish species. More than 50 such weirs are scattered along the river Viaur. Copyright: Loïc Tudesque. et al. 2008). Particularly, the smallest or the largest species should be more strongly affected if dispersal ability or effective population size respectively explain sensitivity to fragmentation. We also expect species-specific responses on IBD patterns. Specific responses would depend upon the initial equilibrium between genetic drift and gene flow of each species in the absence of fragmentation (Hutchison and Templeton 1999;Templeton et al. 2001). For all species, IBD in the fragmented river basin should tend towards a pattern whereby genetic drift is higher and gene flow is lower. Finally, we discuss our results in the light of the cultural and historical status of the particular human fragmentation studied here (Fig. 1).

Study area
The area studied is located in Southwestern France and encompasses two river basins: the River Célé and the River Viaur (Fig. 2). We focused on these two river basins because (i) they are geographically very close to one another (their sources are roughly 50 km apart) and (ii) they both belong to the Adour-Garonne river drainage, which suggest a shared geological, colonization and evolutionary histories (Costedoat et al. 2006). These two rivers are very similar in terms of hydrological and chemical characteristics (G. Loot and S. Blanchet, unpublished data). Respectively, the River Célé and the River Viaur are 136 and 168 km long, their drainage areas cover 1350 and 1530 km 2 and their annual mean flow ranges between 7-30 m 3 /s and 8-25 m 3 /s. However, these two rivers differ with regards to the level of anthropogenic fragmentation. The river Viaur is highly fragmented and will be hereafter referred to as the 'fragmented landscape'. More than 50 small mill weirs (2-3 m high; see Fig. 1) are scattered along its main channel (approximately one weir every 2-3 km; Grenouillet et al. 2008;Poulet 2007). Most of these structures were built during the 15th century (it seems that weir building even began during the 12th century, see also Raeymaekers et al. 2009) and very few of them have been equipped with fish ladders. Based on non-standardized interviews, these weirs (even if often non-functional) are often viewed as parts of the local culture and heritage (G. Loot and S. Blanchet, unpublished data). In the first half of the 20th century, two hydroelectric dams (30 m high) were built; one located 30 km from the source and the other located 80 km from the source of the river (Fig. 2). The River Célé is weakly fragmented and will be hereafter referred to as the 'continuous landscape'. Ten small weirs (2-3 m high) were constructed on the River Célé during the 20th century but most were equipped with fish ladders, enabling nearly all of the fish species to disperse between river segments (Fig. 2).

Sampling design
We sampled 11 sites on the river Célé (from C1 to C11; Fig. 2) and 10 sites on the river Viaur (from V1 to V10; Fig. 2). Because the genetic structure of freshwater fish is often spatially structured along the upstream-downstream gradient (Hanfling and Weetman 2006;Raeymaekers et al. 2008), the sites were directly on the main channel of each river and were chosen at regular distances along the gradient to cover this whole upstream-downstream gradient (Fig. 2). The riparian distance between each sampling site and the river source was calculated using a geo-referenced map (ArcView v.9ª, ESRI, Redlands, CA, USA). A maximum of 24 adults per species and per site were caught by electrofishing (Table 1). For each individual, a small piece of pelvic fins was collected and preserved in 70% ethanol. Two fish species (L. leuciscus and L. cephalus) were absent from some sampling sites (i.e. the upper sampling sites, Table 1) while the other two species (G. gobio and P. phoxinus) were present all along the sampled gradient. At each site we also estimated the density of each species following a standardized two-pass electrofishing method (Bohlin et al. 1989).

Figure 2
Map showing the geographical position of the river Cé lé (the continuous landscape) and the river Viaur (the fragmented landscape). Both rivers belong to the Adour-Garonne basin drainage. This map also shows the sampling sites (black dots) within each river. Weirs (small lines) and hydro-electric dams (large lines in the river Viaur) are also shown. However, note that not all weirs (over 50) are represented on the river Viaur.

Biological models
All four fish species are widely distributed in Europe and belong to the Cyprinids family, which is the most diverse family in Europe ). These species varied strikingly according to their maximal body length (Keith and Allardi 2001). The largest is L. cephalus which can reach a maximum body length of 600 mm, followed by L. leuciscus (400 mm), G. gobio (200 mm) and P. phoxinus (140 mm) (Keith and Allardi 2001). The average body lengths of the sampled fish were highly correlated with the maximal body lengths reported in the literature (r 2 = 0.96, n = 4; see also Table 1). All these species are mainly insectivorous, with some feeding in the water column (L. cephalus and P. phoxinus) and the other preferentially feeding on the bottom (L. leuciscus and G. gobio) (Keith and Allardi 2001). By focusing on species belonging to the same trophic level and the same family, we should limit potential comparative biases inherent to phylogenetic constraints and to phenotypic attributes other than fish body size (Henle et al. 2004). Finally, by focusing on fish species that are not targeted by anglers, we limit the possibility of stocking effects and uninformed translocation between river drainages.
In fish, as in most animals, body length is positively correlated to dispersal ability at the interspecific level, and body size can hence be viewed as a good proxy of dispersal capability (Hugueny 1989;Jenkins et al. 2007). Although we lack of precise data on the dispersal ability of each of these four species, there are evidences that L. cephalus and L. leuciscus are good dispersers as they are both able to perform long-distance upstream movements (up to 10 km) and to migrate into section upstream weirs (Clough and Beaumont 1998;Lucas 2000;Bolland et al. 2008;De Leeuw and Winter 2008). On the contrary, the two other species (G. gobio and P. phoxinus) seems to have lower dispersal ability (Holthe et al. 2005;Knaepkens et al. 2007), which is not surprising given their relatively small body sizes. Moreover, a negative relationship between species body size and population density is expected at the interspecific level (Cotgreave 1993;Blackburn and Gaston 1997). Accordingly, the two largest species (L. cephalus and L. leuciscus) had lower densities than the two smallest species (G. gobio and P. Phoxinus) (see details in Table 1), which suggested a lower effective population size for these two largest species. Effective population sizes were estimated using microsatellites data (see below the section Descriptive genetic analyses). Both dispersal ability and effective population size are considered as important predictors of vulnerability to fragmentation; the relative role of both characteristics will be considered further in the discussion. This table shows, for each landscape (fragmented or continuous), the proportion of sampling sites that was occupied, the mean sample size, the mean population density, the mean effective population size, and the mean body length (total body length in mm) for each of the four sampled species (95% CI, 95% confidence interval).

Genetic analyses
Genomic DNA was extracted from the pelvic fins using a salt-extraction protocol (Aljanabi and Martinez 1997). We employed a cross-amplification approach for selecting the loci for each species. Using an initial set of 40 markers developed for target and sister species (see Table S1), we cross-amplified the 40 markers for each of the four species and we conserved only the markers that displayed highly readable and repeatable profiles. These crossamplified markers were then tested for the presence of null-alleles using the software micro-checker 2.3 (Van Oosterhout et al. 2004). The loci that showed evidence of the presence of null-alleles were dropped out from our final choice. This approach resulted in the selection of a set of 8-15 loci according to the species [L. cephalus (n = 10); L. leuciscus (n = 15); G. gobio (n = 8), P. phoxinus (n = 8); Table S1]. The selected loci were co-amplified using the QIAGEN Ò Multiplex PCR Kit (Qiagen, Valencia, CA, USA). Polymerase chain reaction (PCR) reactions were carried out in a 10 lL final volume containing 5-20 ng of genomic DNA, 5 lL of 2· QIAGEN Multiplex PCR Master Mix, and locus-specific optimized combination of primers (detailed recipes are available upon request). PCR amplifications were performed in a Mastercycler PCR apparatus (Eppendorf Ò , Hauppauge, NY, USA) under the following conditions: 15 min at 95°C followed by 30 cycles of 1 min at 94°C, 1 min at 60°C and 1 min at 72°C and finally followed by a 60 min final elongation step at 72°C. Amplified fragments were then separated on an ABI PRISMÔ 3730 automated capillary sequencer (Applied Biosystems). Allelic sizes were then scored using GENEMAPPERÔ v.4.0 (Applied Biosystems, Foster City, CA, USA).

Descriptive analyses
For each species, locus-by-locus heterozygosity (observed and expected) as well as F is estimates were calculated using GENETIX version 4.05.2 (Belkhir et al. 2002). Departure from the Hardy-Weinberg equilibrium was tested using GENEPOP version 3.4 (Raymond and Rousset 1995). Linkage disequilibrium between all pairs of loci was tested in the program Fstat version 2.9.3.2 (Goudet 1995). We used the software NeEstimator 1.3 (Queensland Government, Brisbane, Australia) to estimate the effective population size for each species. We used a point estimation method based on linkage/gametic disequilibrium (Hill 1981). This method was preferred over the 'heterozygote excess method' (Pudovkin et al. 1996) because the later was less effective for converging. It was applied on each sampling site independently and we calculated a mean effective population size over all sampling sites. Point estimation methods must be interpreted cautiously and more sophisticated methods involving computationally intensive algorithms and/or multiple cohorts are often preferred to estimate properly effective population size (Waples 1989). The results presented here will therefore be used only for inter-species comparison but not for inter-landscapes comparison.
In addition, we used unpaired t-test (two-tails) to test if, within-species, density, effective population size and mean body size differed between the two landscapes.

Effect of fragmentation on genetic diversity
Two measures of genetic diversity were computed for each sampling site and for each species. AR corrected for the minimum sampling size was computed using Fstat version 2.9.3.2 (Goudet 1995), and observed heterozygosity (Hobs) was computed using GENETIX version 4.05.2 (Belkhir et al. 2002). For each species we computed generalized linear models (GLM) with either AR or Hobs as the dependent variables, landscape status (fragmented or continuous landscapes) as a categorical predictor and the position of each sampling site from the river source as a continuous predictor. As observed elsewhere (e.g. Raeymaekers et al. 2008) we expected both AR and Hobs to be positively correlated to the distance from the river's source. We also fitted the two-way interaction to test for a significant effect of the landscape status on the slope of the relationship between either AR or Hobs and the distance from the source. For all models we assumed a Gaussian error-term distribution and the sampling site was the replicate unit.
In addition to testing for a significant effect, we also used a meta-analytic framework to investigate how the magnitude of the effect of fragmentation (i.e. the effect size; Nakagawa and Cuthill 2007) differed between species. From mean and standard deviation data, we calculated the Hedges' d-value for each species and each measure of genetic diversity (see Rosenberg et al. 2000 for formulas). The Hedges' d-value is a standardized effect size that is comprised between )¥ and +¥, with a negative value representing effects where the control group (here the continuous landscape) attains a greater value than the experimental group (here the fragmented landscape) (Rosenberg et al. 2000). We then calculated the cumulative effect size ( E) and the 95% confidence interval for all pooled species. This procedure allows the calculation of a mean effect size weighted for differential sampling sizes between species (Rosenberg et al. 2000). We then verified whether the effect size calculated for each species fell within the 95% confidence interval calculated for all species (i.e. did an individual effect significantly differ from the mean effect?). Finally, we evaluated if the set of effect sizes calculated for each species was homogeneous or not. The total heterogeneity of a sample (Q t ) was calculated as described in Rosenberg et al. (2000) and its significance was tested using chi-square statistics.

Effect of fragmentation on genetic structure
We first investigated the population genetic structure by estimating the F st (Weir and Cockerham 1984). We calculated global F st values for each locus and each species using the program Fstat version 2.9.3.2 (Goudet 1995). For each species, we computed a GLM with F st as the dependent variable and landscape status (fragmented or continuous landscapes) as a categorical predictor. For all models. we assumed a Gaussian error-term distribution and the locus was the replicate unit. The meta-analytic framework we described above was used to compare the magnitude of the effect of fragmentation on the F st values for each species.
We then investigated the spatial genetic structure using the IBD framework (Hutchison and Templeton 1999;Templeton et al. 2001). Pairwise F st between all pairs of populations were computed for each landscape and each species separately using GENETIX version 4.05.2 (Belkhir et al. 2002). Pairwise riparian distances between sites were computed using the Geographical Information System. We first graphically inspected the linear relationship between F st /(1 ) F st ) and the distance between pairs of sampling sites as proposed by Hutchison and Templeton (1999) to infer the migration drift equilibrium for each species in each landscape. Because of the nonindependence of the replicate unit (i.e. a pair of sampling sites), we assessed the strength (as the Pearson coefficient of correlation) and the significance of each linear relationship using a re-sampling procedure (i.e. procedure analogous to the Mantel test; Manly 1997). We then investigated the effect of fragmentation on the slope of the linear relationship between F st /(1 ) F st ) and the distance between a pair of sampling sites using another re-sampling procedure in which the slope of the fragmented landscape for a given species was compared with the slope of the continuous landscape for this species (Manly 1997; see also Epps et al. 2005). If gene flow is reduced between sub-populations, we would expect a significantly higher slope in the fragmented than in the continuous landscape.

Descriptive analyses
We did not find evidence for any significant heterozygous deficit (after Bonferroni correction) for any of the loci that were considered (Table S1). No pair of loci displayed significant linkage disequilibrium (results not shown). As expected, population density correlates quite strongly to effective population size (r 2 = 0.766, n = 8, P < 0.005; see Table 1). Roughly, the two largest bodied-species (L. cephalus and L. leuciscus) tended to have lower effective population size than the two smallest-bodied species (G. gobio and P. phoxinus) (Table 1). However, within-species, we detected no significant differences between landscapes in term of mean body length, population density and effective population size (t-test, all P > 0.05, see Table 1).

Effect of fragmentation on genetic diversity
In case of the AR, we found significant differences between the fragmented versus the continuous landscape for all the species except P. phoxinus (Table 2; Fig. 3A). In all significant cases, the AR was lower in the fragmented landscape (Fig. 3A). It is noteworthy that the AR (irrespective of the landscape) increases as the maximum body size of the considered species decreases (Fig. 3A, see also Frankham 1996). Over all species, the magnitude of the effect imposed by fragmentation was significantly different from 0 (see the 95% confidence interval in Fig. 4A), but was significantly heterogeneous among species (Q t = 24.44, d.f. = 3, P < 0.001). This means that not all species respond equally to river fragmentation. Specifically, the effect sizes for L. leuciscus and G. gobio were greater than the mean whereas the effect size was lower than the mean for P. phoxinus (Fig. 4A). For one species (G. gobio), the AR was not affected by the distance of the sampling site from the river's source (Table 2D). For L. leuciscus and P. phoxinus, there was a significant positive correlation between AR and distance from the river's source (Table 2D; see also Fig. S1C,G). We also found a significant interaction between 'fragmented status' and 'distance from the river's source' for L. cephalus (Table 2A). This indicated that the relationship between AR and distance from the source was positive and significant for the continuous landscape (r = 0.75, n = 8, P = 0.033) but not for the fragmented landscape (r = 0.03, n = 8, P = 0.936; see Fig. S1A).
Concerning observed heterozygosity, we found significant differences between the fragmented versus the continuous landscape for all species (Table 2; Fig. 3B). In all cases, the observed heterozygosity was lower in the fragmented landscape (Fig. 3B). Over all species, the magnitude of the effect imposed by fragmentation was significantly different from 0 (Fig. 4B), and did not differ significantly among species (Q t = 1.58, d.f. = 3, P = 0.661). The effect sizes for P. phoxinus tended to be slightly lower than for the other species (Fig. 4B). The effect of the distance of the sampling site from the river's source on Hobs was significant only for P. phoxinus (Table 2D). For this species, there was also a significant interaction between 'fragmented status' and 'distance from the river's source' (Table 2D). This showed that the relationship between AR and distance from the source was positive for both the continuous and the fragmented landscape, but the slope of the relationship was steeper in the continuous landscape (Fig. S1H).

Effect of fragmentation on genetic structure
Concerning F st , we found significant differences between the continuous and fragmented landscape for the two largest species [L. cephalus (F 1,18 = 9.828, P = 0.006) and L. leuciscus (F 1,28 = 10.844, P = 0.003)] but not for the two smallest species [G. gobio (F 1,14 = 1.964, P = 0.182) and P. phoxinus (F 1,14 = 0.159, P = 0.696)]. In both significant cases, the level of population differentiation (i.e. F st ) was higher in the fragmented landscape (Fig. 3C). Over all species, the magnitude of the effect was significantly different from zero (Fig. 4C), and did not differ significantly between species (Q t = 3.518, d.f. = 3, P = 0.320). There was a positive tendency between the maximum body size of each species and the magnitude of the effect (Fig. 4C).
Concerning IBD, there were two species (L. cephalus and P. phoxinus) for which there were significant relationships between genetic differentiation and geographical distance for both the continuous and the fragmented landscapes (Table 3; Fig. 5A,D). In both cases, we detected no significant differences between the slope of the continuous landscape and the slope of the fragmented landscape (re-sampling tests, L. cephalus, P = 0.711; P. phoxinus, P = 0.118; Fig. 5A,D). For the other two species (L. leuciscus and G. gobio), we found a significant relationship between genetic differentiation and geographical distance only for the fragmented landscape (Table 3; Fig. 5B,C). The pattern was not significant for the continuous landscape (Table 3; Fig. 5B,C) indicating that distance between sampling sites had no influence on the genetic structure of these species in the continuous landscape. In both cases, there was a significant difference between the slope of the continuous landscape and the slope of the fragmented landscape (re-sampling tests, L. leuciscus, P = 0.014; G. gobio, P = 0.002; Fig. 5B,C).

Discussion
Using a meta-analytic framework, we highlight significant differences between fragmented and continuous landscapes, both for genetic diversity and the genetic structure of the set of species studied. In streams and rivers, gene flow is thought to be quasi-unidirectional as it is mainly driven by the direction of the water flow (Fraser et al. 2004;Hanfling andWeetman 2006, Crispo et al. 2006). This downstream-biased gene flow can impede our ability to disentangle the effects of historical colonization, man-made migration barriers and downstream-biased gene flow (Raeymaekers et al. 2008). We ought to bypass this difficulty by using a synchronic comparison, i.e. a comparison between a fragmented and a continuous landscape that share similar evolutionary and geological histories (Hoehn et al. 2007;Hendry et al. 2008).
Although we can not totally discard that the pattern observed here are a hangover from the past or downstream-biased gene flow (Hanfling and Weetman 2006), such a comparison makes it likely that the betweenlandscape differences reported here are (at least in part) the result of anthropogenic fragmentation. Thus, our results support numerous studies on many animal species showing that human fragmentation modifies the genetic integrity of natural populations, notably in decreasing gene flow among populations and genetic diversity within the remnant population patch (e.g. Van den Bussche et al. 2003;Alo and Turner 2005;Trizio et al. 2005;Hoehn et al. 2007;Schiffer et al. 2007). It is worth noting that our preliminary results on the ecological characteristics of these fish populations (see Table 1) indicate that, from an ecological point of view, these populations have not suffered intensively from the changes in genetic diversity and structure we report here. Therefore, it is possible that these species have undergone evolutionary changes that allow them to cope with fragmentation (and the associated genetic changes), and/or that genetic consequences reported here are not strong enough to influence significantly the biological fitness of populations. More interestingly, we detected obvious species-specific responses to fragmentation. This important result is consistent with previous papers showing that not all species are equally affected by fragmentation (Davies et al. 2000;Henle et al. 2004;Ewers and Didham 2006;Hoehn et al. 2007), and hence suggest that the effect of fragmentation in streams and rivers can not be generalized from one fish species to another. Contrary to certain expectations (Ewers and Didham 2006), we failed to verify the hypothesis that large-bodied species are less sensitive to habitat fragmentation. On the contrary, all genetic indices clearly indicated that the smallest-bodied species (P. phoxinus; Table 1) was the least affected by fragmentation. We further provide evidence that species with intermediate body size might in fact be the most affected by fragmentation, at least at the genetic level. Species differences in the responses to fragmentation were obvious and significant for two discriminant genetic indices; the AR and IBD patterns. For these two genetic indices, L. leuciscus and G. gobio (which are the two intermediate-bodied species; Table 1) were most affected by the presence of weirs. Indeed, in the case of AR, these two species had an effect size that was significantly lower than the mean (Fig. 4A). With regard to IBD, these two species showed patterns that significantly differed between landscapes: a pattern evocative of high gene flow and low genetic drift in the continuous landscape, and a pattern evocative of a Figure 3 Bar plots showing the mean allelic richness (A), the mean observed heterozygosity (B) and the mean F st (C) for each of the four species studied independently and each landscape (continuous landscape, black bars; fragmented landscape, grey bars). The four species are ranked from the largest to the smallest from left to right (see also Table 1). The stars (*P < 0.05; **P < 0.01; ***P < 0.001) indicated significant differences (generalized linear models; see Table 2) between the two landscapes. n.s., non-significant differences. Error bars are standard errors.

Species-specific responses to fragmentation
Blanchet et al.
regional equilibrium between gene flow and drift in the fragmented landscape (Hutchison and Templeton 1999). As expected, this latter result suggests that barriers limit gene flow between populations and favour genetic drift within them. Very interestingly, this demonstrates that human fragmentation has disturbed the initial balance between genetic drift and gene flow, and has constrained all these fish species in a similar population structure whereby an equilibrium between genetic drift and gene flow is attained (Raeymaekers et al. 2008). We can hence conclude that (i) not all species react similarly to fragmentation, (ii) the smallest sized species (P. phoxinus) is the least affected by fragmentation and (iii) there is a tendency for the intermediate-bodied species (L. leuciscus and G. gobio) to be the more affected.
Many traits have been proposed as potential predictors of species sensitivity to fragmentation (reviewed in Ewers and Didham 2006;Henle et al. 2004). Predictions regarding body size depend upon the level of organization being considered. For instance, it is often claimed that largebodied species are at greater risk of extinction than smallbodied species because large-bodied species are often at higher trophic levels and hence, have lower population size and more unstable population dynamics (Gaston and Blackburn 1996;Henle et al. 2004). However, we control for this possibility by selecting species of the same trophic level, which orientates our prediction towards the hypothesis that large-bodied species are better dispersers and are therefore less sensitive to fragmentation (see Fig. 2Q,R in Ewers and Didham 2006). As stated above, this hypothesis was not supported, as the smallest-bodied species was the least sensitive to fragmentation. In fish, small-bodied species generally have a narrower habitat range than large-bodied species (Pyron 1999;Rosenfield 2002). Therefore, we can speculate that small-bodied species are not affected by fragmentation because their home range is smaller than the distance separating two barriers. A nonmutually exclusive speculation relates to the possibility that small-bodied species also have large effective population size (which is verified in our dataset; Table 1), which might blur the effect of genetic drift within remnant population patches (Lawton et al. 1994;Davies et al. 2000). Because of the strong covariance between body size and several other traits such as home range, population size and dispersal ability, it is currently difficult to ascertain  The first line of each row reports the coefficient of correlation (Pearson r, 1000 iterations) of the relationships and of its standard deviation (in brackets). The second line of each row is the P-value (italic) associated to such coefficients. Bold P-values are significant (P < 0.05). Results are shown for each landscape (fragmented or continuous) and each of the four species [from (A) to (D), in decreasing order of body length; see Table 1] separately.
which mechanisms underlie the observed pattern for the smallest-bodied species of our dataset (Henle et al. 2004;Theodorou et al. 2009). Nevertheless, we claim that effective population size cannot explain why the two intermediate-bodied species are the most impacted by human fragmentation. Indeed, because these two species varied greatly in their respective abundance and effective population size (Table 1) while being equally affected by fragmentation, it is more likely that this common pattern arises from other traits such as dispersal ability or behaviour. The idea that species with intermediate dispersal ability might be more affected by fragmentation than species with either high or low dispersal has already been predicted theoretically (Fahrig 1998). Empirically, this has been demonstrated in a butterfly community in which species with intermediate mobility were more likely to decline in abundance following habitat fragmentation than butterflies with either low or high mobility (Thomas 2000). Although this remains to be tested, empirical data, however, support the idea that L. leuciscus and G. gobio display different dispersal ability (see the Biological models section; Knaepkens et al. 2007;Lucas 2000). It is therefore highly probable that movement behaviour and life-history strategies linked to feeding habits for instance (L. leuciscus and G. gobio feed on the bottom while the other two species feed on the water column) interact with dispersal ability to explain species sensitivity to fragmentation. For instance, L. cephalus and L. leuciscus are known to have similar dispersal ability (Lucas 2000;De Leeuw and Winter 2008;Ovidio et al. 2008), but are however differentially affected by weirs (our study). We have shown that in absence of weirs, L. cephalus population display a significant IBD pattern in the unfragmented landscape, which might indicate that they are more sedentary than L. leuciscus. This means that, in that case, dispersal behaviour rather than dispersal ability per se could explain species discrepancy in sensitivity to weirs. Given the descriptive nature of our study, further experimental work is needed to sort out the relative role of dispersal ability, dispersal behaviour and effective population size as drivers of species sensitivity to fragmentation. However, this means that certain phenotypic traits make species more or less sensitive to fragmentation, and that evolutionary processes (e.g. local adaptation by natural selection) can themselves differentially affect species.

Implications for management
That fragmentation differentially affects fish species means that not all populations are concerned by fragmentation, and hence that programmes of restoration can be prioritized. For instance, programmes aimed at avoiding the genetic erosion of species (e.g. hatchery and stocking programmes) do not need to be undertaken for all species at the same time, and it is here clear that the first measures that are needed must concern L. leuciscus and G. gobio. Being able to prioritize conservation efforts is very important, notably when funding is limited. Continuous lines indicate significant univariate relationships (P < 0.05; see Table 3), while dotted lines are not significant.
In case of river fragmentation by dams or weirs, four practical management solutions can be undertaken to restore the genetic integrity of populations: (i) the stocking of fish produced in a hatchery, (ii) the translocation of individuals from site to site, (iii) the construction of fish passages and (iv) the un-damming of rivers (Bednarek 2001;Palmer et al. 2005Palmer et al. , 2008Raeymaekers et al. 2009). Stocking is often associated with other evolutionary problems that make this solution undesirable (Blanchet et al. 2008;Frankham 2008). Translocation is too time-consuming to be efficient. For weirs such as those in the river Viaur (the fragmented landscape), un-damming is not a solution because these weirs date back to the Middle Ages (see also Raeymaekers et al. 2009) and are thus part of the local heritage and culture. However, there is some evidence that building fish passages can be an efficient restoration tool to preserve both the genetic integrity of fish species and the authenticity of the weirs (Raeymaekers et al. 2009). In the river Viaur, there are more than 70 weirs that have been identified. Because of the high financial cost of constructing fish passages on all the weirs, a priority must be given to weirs that have a high impact on the fish genetic structure, and the passages must be efficient (i.e. all fish can benefit from the passage). Prioritization in weirs restoration can be evaluated qualitatively and quantitatively by using the same procedure as the one used in Raeymaekers et al. (2009). According to the 'umbrella species' principle (Lambeck 1997), L. leuciscus and G. gobio appear as good candidates species for investigating prioritization in weirs restoration in the near future. The construction of fish passages should, however, follow the biological requirements of most species inhabiting such river systems. Following the recommendations made by Raeymaekers et al. (2009) in a similar system, we can expect that priority will be given to the highest weirs. However, the construction of fish passages should follow the biological requirements of most species inhabiting such river systems. Following the recommendations made by Raeymaekers et al. (2009) in a similar system, we propose that priority will be given to the highest weirs.
It is worth remembering that the river Viaur also contains two recently built hydro-electrics dams that are more than 30-m high. The effect of these dams (relative to the effect of weirs) on the evolution and ecology of fish populations is currently not known, and it should be the focus of our future studies. However, the flow rate of the river Viaur has been considerably altered by these two dams, being severely decreased downstream of the dams (i.e. a major part of the flow is diverted for agricultural purposes), which has increased sedimentation above small weirs (S. Blanchet and G. Loot personal observation). If fish ladders are installed in the river Viaur, they will need to be built according to the flow rate imposed by dams to ensure that these ladders function well (Larinier et al. 2002).
The approach used here for evaluating the effect of fragmentation on fish populations differs from previous studies in that (i) it focuses on four fish species and (ii) it uses a synchronic sampling design (Hendry et al. 2008). These characteristics make this study highly valuable from a conservation perspective, as it allows inferring the effect of fragmentation at a 'community-wide' level. Furthermore, it allows detecting subtle effects that can be blurred when more traditional approaches are used (i.e. when only a single-fragmented landscape is considered). However, this approach implicitly assumes that differences observed between landscapes are mainly the result of human fragmentation. This assumption is hard to ascertain and as a consequence, our study should be viewed as a basis for future studies considering simultaneously the effects of historical and contemporary histories on the genetic structure of wild fish populations. The use of a larger spatial scale sampling design and/or different genetic markers could be use for reaching this goal. This improvement should get a more complete picture of the relative role of humans in shaping the nowadays genetic structure of wild fish populations, and hence should improve management policies.

Conclusion
To conclude, our study demonstrates the importance of considering several species when investigating the effect of fragmentation mediated by humans. We have shown that a biological trait such as body size is important to understand the influence of fragmentation on the genetic integrity of fish populations. This means that the evolutionary outcomes of fragmentation will depend upon the traits that characterize species. Moreover, this implies that restoration programmes will need to be prioritized and targeted according to the specificity of each species. Future experimental studies are needed to ascertain the traits that affect the genetic sensitivity to fragmentation, and to test restoration programmes that can benefit a spectrum of species that are functionally dissimilar.

Supporting Information
Additional Supporting Information may be found in the online version of this article: Figure S1. Relationships between distance from the source and genetic diversity (AR and Obs H) for the four species studied here. Table S1. Description of the microsatellites loci used in this study.
Please note: Wiley-Blackwell are not responsible for the content or functionality of any supporting materials supplied by the authors. Any queries (other than missing material) should be directed to the corresponding author for the article.