Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Novelties in Hybrid Zones: Crossroads between Population Genomic and Ecological Approaches

  • Caroline Costedoat ,

    To whom correspondence should be addressed. E-mail: costedoa@up.univ-mrs.fr

    Affiliation Evolution Genome Environment (EA 3781 EGEE), Université de Provence, Marseille, France

  • Nicolas Pech,

    Affiliation Evolution Genome Environment (EA 3781 EGEE), Université de Provence, Marseille, France

  • Rémi Chappaz,

    Affiliation Evolution Genome Environment (EA 3781 EGEE), Université de Provence, Marseille, France

  • André Gilles

    Affiliation Evolution Genome Environment (EA 3781 EGEE), Université de Provence, Marseille, France

Abstract

Background

Interspecific hybridization is widespread, occurring in a taxonomically diverse array of species. The Cyprinidae family, which displays more than 30% hybridization, is a good candidate for studies of processes underlying isolation and speciation, such as genetic exchange between previously isolated lineages. This is particularly relevant in the case of recent hybridization between an invasive species, Chondrostoma nasus nasus (from Eastern Europe), and C. toxostoma toxostoma (a threatened species endemic to southern France), in which bidirectional introgressive hybridization has been demonstrated.

Methodology/Principal Findings

We studied 128 specimens from reference populations and 1495 hybrid zone specimens (two years of sampling and four stations), using five molecular markers (one mitochondrial gene, four nuclear introns), morphology (meristic and plastic characters) and life history traits (weight, size, coefficient of condition, sex, age, shoaling). We identified 65 hybrid combinations and visualized spatial and temporal changes in composition. The direction of mitochondrial introgression was density-dependent in favor of the rarer species and we demonstrate that the sexual selection hypothesis is a preponderant explanation in the asymmetry of introgression. Despite genomic evolution in the hybrid zone, convergence was observed for body shape and coefficient of condition, indicating changes in foraging behavior with respect to reference populations, reflecting strong environmental pressure.

Conclusions/Significance

The complex rules of hybrid zone dynamics are established very early in the contact zone. We propose “inheritance from the rare species” as a new evolutionary hypothesis for animal models. The endemic species was not assimilated by the invasive species. Survival rates for this species were highest in the middle of the river (the warmest part) due to a trade-off between food availability and fecundity. The environment-independent hybrid combination may result from nuclear-mitochondrial interactions involving the Tpi1b gene or a gene linked to this gene (Chromosome 16). This genomic region is also responsible for shoaling behavior in Danio rerio and is a promising zone for studies of changes in population dynamics and advances in integrated studies of hybrid zones.

Introduction

Hybridization between species and subsequent backcrossing may lead to the transfer of genes between species. This introgression of new genetic information provides new variability and is a major factor in subsequent evolution. However, the genetic assimilation of a rare species by more common congeners [e.g. in 1], [2], [3] has important implications for conservation. Genetic models have indicated that this process may occur extremely quickly, and many examples have been described in both plants [4], [5] and animals [3], [6]. Hybridization is therefore now considered a serious threat to endangered species. Wilson [7] reported that 38% of the 292 North American fishes listed as rare, threatened, or vulnerable in the Red Data Book published by the International Union for the Conservation of Nature and Natural Resources were threatened by hybridization. Scribner et al.[8] reviewed 158 articles covering 19 fish families from North America, Europe and Asia. In “natural” conditions (i.e. excluding fish farming conditions), species from the Cyprinidae included more hybrid pairs than any other family, with more than 30% hybrids, whereas hybridization rates were 8% in the Salmonidae and only 3% in the Cichlidae. The preponderance of this phenomenon in this family, its worldwide distribution and the inclusion in this family of a model species (Danio rerio) render this family particularly interesting for studies of hybridization phenomena in an evolutionary context.

Chondrostoma t. toxostoma is a threatened, protected endemic cyprinid species from southern France. Part of its distribution range has been colonized by the invasive Chondrostoma n. nasus, from Eastern Europe. This colonization has made hybridization between these two species possible, particularly in the Rhone Basin. In previous studies [9], [10], we reported bidirectional introgression between the two Chondrostoma species in the River Durance. The structuring of the hybrid zone was strongly linked to progressive urbanization and an increase in human activities along this river, with the construction of various dams preventing fish from swimming upstream or downstream. No continuum was found between the species inhabiting the tributaries (C. t. toxostoma) upstream and the species inhabiting the river (C. n. nasus) downstream. There was therefore no hybrid gradient based on the ecological preferences of the species, as would have been expected based on the route followed by C. n. nasus during the colonization process. These findings conflict with those for two other well known cyprinid models: a gradient hybrid zone has been reported between Barbus barbus and Barbus meridionalis in the Lergue (Europe) [11] and asymmetric introgression in two Luxilus species (America) [12] has been accounted for by historic effects (refuge populations and dispersal route).

Dams delimited four main stations displaying independent evolution as it was almost impossible for the fish to move from one station to another.

We evaluated the contributions of genome and environment (Eco-Geno context) to spatial and temporal changes in introgressive hybridization, to enable us to decipher hybridization phenomena and the rules governing them. We used species-specific molecular diagnosis profiles to quantify levels of hybridization and introgression in the two species. Changes in size may be associated with changes in shape, and may have functional significance [13]. However, few studies have considered changes in ecological characters or morphological analyses of hybridization, and most such studies have focused on a limited number of hybrid classes [e.g. in 14]. We studied each of the two Chondrostoma species in allopatric and parapatric populations (reference values) and then focused on the hybrid zone. The River Durance was divided into four sampling stations, from upstream to downstream, which were sampled in 2001 and in 2002. All 1623 specimens were identified on the basis of molecular markers (mtDNA cytochrome b and four nuclear introns), and life history traits, such as sex, age, weight, size, coefficient of condition, shoaling behavior and morphology (meristic and plastic characters), were analyzed.

Results

Specimen identification and dynamics of the hybrid zone

All 1623 specimens were identified based on the five molecular markers (cytochrome b plus the four introns). The allele dynamics for each marker are presented in Dataset S1, figures S1, S2, S3, S4; Tables S1, S2. We developed an identification code based on the results obtained for each marker (figure 1), defining 162 combinations: the two parental combinations plus 160 hybrid categories. As expected, all specimens from allopatric and parapatric populations were identified as “pure” specimens with our markers. In the hybrid zone, we observed 67 of the 162 possible combinations: the two parental combinations (H5 and T5) and 65 hybrid combinations.

thumbnail
Figure 1. IDENTIFICATION CODE.

Some examples of genomic combinations based on the alleles for each of the five molecular markers (1 mitochondrial gene and 4 nuclear introns).

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

We tested the null hypothesis of Hardy-Weinberg equilibrium, using chi-squared goodness-of-fit tests, with Benjamini-Hochberg correction for multiple comparisons. We classified the genomic combinations into three major groups: combinations overrepresented under the null hypothesis of Hardy-Weinberg equilibrium, those underrepresented and those observed in the expected proportions.

Forty-one of the 65 hybrid combinations observed (considering the entire data set) presented nasus mtDNA and 24, toxostoma mtDNA (significant difference: chi-square = 6.73; df = 1; p = 0.009). Both types of F1 hybrid were present, but at different frequencies (Table 1): 31 Hi4 specimens versus 5 Ti4 specimens (chi-square = 18.78; df = 1; p = 1.5×10−5). We then investigated the temporal and spatial dynamics of evolution in the hybrid zone, by studying the data sets for each year (2001 and 2002) and sampling station (Buech, Manosque, Pertuis and Cavaillon) separately (Table 1).

thumbnail
Table 1. OBSERVED AND EXPECTED EFFECTIVE SIZE OF EACH GENOMIC COMBINATIONS PRESENT IN THE HYBRID ZONE

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

Alternate asymmetric mtDNA introgression

We investigated the distribution of the hybrid combinations as a function of the abundance of the parental species. The mtDNA of the hybrid combinations (particularly for the overrepresented combinations) originated in most cases from the rarer of the two “parental” species. Indeed, for the total data set, more C. t. toxostoma individuals (677) than C. n. nasus individuals (251) were observed, whereas the hybrid combinations were more likely to present C. n. nasus mtDNA (188 vs 126; p<10−15) (Table 2). The same trend was observed if the data for each sampling year and for each sampling station were considered separately (Table 2). We call this phenomenon “inheritance from the rare species”.

thumbnail
Table 2. DISTRIBUTION OF THE OR TOXOSTOMA MTDNA MODALITY.

https://doi.org/10.1371/journal.pone.0000357.t002

Morphology and hybridization

We calculated the coordinates of each hybrid zone specimen on the discriminant axis defined by the reference populations (figure 2). The morphological ranges of the H5 and T5 classes were larger than those of the two species sampled in allopatry and parapatry and, furthermore, these ranges overlapped. The morphological range of the Hi4 specimens (F1 with nasus mtDNA) was intermediate between those of the H5 and T5 specimens (figure 2). See Dataset S2, figures S5 and S6.

thumbnail
Figure 2. VISUALIZATION OF THE VARIANCE FOR MORPHOLOGY.

Discriminant analysis based on meristic characters (x-axis) as a function of genetic class (y-axis). -1 = underrepresented combinations; 0 = combinations present in the expected proportions; 99 = overrepresented combinations other than T5, H5,T4i, Hi4 and H4T, which are represented separately (see the text for more detail). A) Meristic data set; B) Plastic data set on the global data set.

https://doi.org/10.1371/journal.pone.0000357.g002

Based on plastic morphology, H5 and T5 specimens in the hybrid zone were mostly classified in hybrid classes, rather than with their respective reference populations (Table S3), illustrating the difference in morphology between individuals in the hybrid zone and reference allopatry/parapatry populations (not observed with meristic characters). We therefore assessed the correlation between gene dilution (i.e. the proportion of alleles from each species present in each genomic combination) and morphology. The coefficient of determination was significant for both meristic and plastic characters (R2 = 0.804 for meristic characters; and R2 = 0.802 for plastic characters, figure 3).

thumbnail
Figure 3. CORRELATION BETWEEN MORPHOLOGY AND GENE DILUTION.

Correlation between morphology (meristic at the top and plastic at the bottom) and the degree of introgression (percentage of genes belonging to each of the two species) or dilution.

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

All genomic categories displayed a shift in body shape from those of the reference populations (figure 2). This shift corresponded to a movement towards the nasus type (figure 2). A similar shift was observed for individuals of all sizes and ages and from all stations (data not shown). We visualized this change, using deformation grids for different comparisons (figure 4). The main differences between the mean reference nasus individual and the mean H5 individual concerned snout shape and body height between the dorsal fin and the pelvic fin, which was smaller in the H5 specimen (figure 4). Surprisingly, the deformation from a reference mean toxostoma to a mean T5 was similar. Both the mean T5 specimen and the mean H5 specimen were more rounded than the corresponding mean reference specimens. All hybrid zone specimens presented similar body shape deformations, regardless of their genomic class. These changes resulted from an effect of the Durance river environment rather than a genomic effect (Table S4).

thumbnail
Figure 4. DEFORMATION GRIDS.

Deformation grids for various “mean specimen” comparisons (based on the 21 landmarks). These diagrams are schematic representations of the morphology encountered in allopatric/parapatric zones or in hybrid zones. The points indicate the main deformations observed on the corresponding grid.

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

Life history traits

In reference populations, the mean coefficient of condition value for C. n. nasus (KC.n.nasus = 12.50) differed significantly from that for C. t. toxostoma (KC.t.toxostoma = 9.84; p<10−15). However, no species effect was found in the hybrid zone (p = 0.16), with KH5 = 11.01; KT5 = 11.22. We can summarize this relationship as follows: KC.n.nasus>KH5 = KT5>KC.t.toxostoma

The greater K value for the reference C. n. nasus specimen than for H5 seems to result principally from the larger size of the H5 specimens (figure 5). Conversely, T5 specimens were smaller than C. t. toxostoma specimens and the difference in weight between the two was smaller than that for H5 and C. n. nasus. (More details in Dataset S2, figure S7, Table S5).

thumbnail
Figure 5. COEFFICIENT OF CONDITION AND CONFIDENCE SURFACES.

Confidence surfaces for mean K-size (a) and mean K-weight (b) plots, based on non-parametric bootstrap (1000 replicates). The circles correspond to 95% of bootstrapped values. H and T correspond to C. n. nasus and C. t. toxostoma, respectively, in the reference populations. The other classes are those of the hybrid zone (T5; H5; -1 = underrepresented hybrid combinations; 0 = hybrid combinations found in the expected proportions; 99 = overrepresented hybrid combinations; T4i = the only hybrid combination for which sufficient specimens were obtained for comparison alone, without pooling). The variance of K is very significantly different (Bartlett B = 63.7; p = 10−10) between classes.

https://doi.org/10.1371/journal.pone.0000357.g005

The oldest specimens (10 years old) were of the C. n. nasus species and the sex ratio in the hybrid zone was about 1 (Tables S6, S7).

Interactions between variables

We tested potential interactions between variables, using chi-square tests (second-order interactions) and log-linear models (third-order interactions) for qualitative variables. The third-order interaction of the saturated log linear model (genomic class×station×sampling year) was significant (dev. = 72.134; p = 5×10−10). The main interaction results are reported Table S8.

Positive interactions were found between H5 and the two extreme stations (the Buech station upstream and the Cavaillon station downstream) and between T5 and the two intermediate stations (Manosque and Pertuis). This is consistent with the distribution of the different genomic groups not following an upstream-downstream gradient. Some hybrid combinations were positively associated with a particular station (e.g. T4i and the Manosque station), and were not present at others. These combinations may be considered to be environment-dependent. In contrast, H4T (positively associated with the Buech station) was overrepresented at three of the four stations and present in the expected proportions at the fourth. This combination may therefore be considered environment-independent.

Discussion

From deciphering Chondrostoma hybrid zone dynamics to the emergence of an evolutionary scenario (figure 6)

thumbnail
Figure 6. GENOMIC EVOLUTION THROUGH THE HYBRID ZONE.

Global summary of the dynamics of the hybrid zone and main variable interactions (chi-square test and log-linear model) with respect to the various sampling years and stations (see the text for more detail). The percentages of the different categories are shown—from T5 in red to H5 in green, and “hybrids” in yellow—separately for each station and each sampling year.

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

Our study confirmed the bidirectionality of introgression between the Chrondrostoma species, via the two F1 hybrid pools identified as Hi4 and Ti4. The hybrid specimens displayed no signs of sterility, with hybrid males producing sperm during the reproductive period and hybrid females presenting normal development, with a micropyle (Stolzenberg personal communication). We identified environment-dependent and environment-independent combinations and describe here the first vertebrate hybrid zone following the “evolutionary novelty” model defined by Arnold [15], with two types of genotype variation encountered within the same species complex and in the same river. However, only one (H4T) of the 65 hybrid combinations observed was clearly environment-independent. This combination consisted of C. t. toxostoma alleles for the Tpi 1b intron and C. n. nasus alleles for all the other molecular markers. Recent studies have highlighted the importance of interactions between the mitochondrial and nuclear compartments, particularly for fitness (amount of available energy [16], [17]), in all biological models. The Tpi 1b gene product is involved in glycolysis. We therefore searched the Danio rerio chromosome map for other genes encoding proteins involved in glycolysis. G3PDH, like Tpi 1b, was found to be present on chromosome 16 (region displayed: 58.99-59.01 cM). We identified no other candidate genes from the Danio rerio map in this chromosome region, but a more precise estimate of the Chondrostoma genome is required to demonstrate such a direct interaction.

The hybrid zone presented a mosaic of genomic combinations evolving over space and time (figure 6). The tendency of some genomic combinations (H5, T5 and some hybrid ones) to be overrepresented may reflect assortative mating, leading to a tendency to display intraspecific shoaling in both species (as suggested by the frequent sampling of specimens of similar genomic combination and age together). Recent studies have highlighted the role of intraspecific and interspecific shoaling in the dynamics of fish populations and its consequences for individual fitness [18], [19], [20]. Such shoaling behavior cannot itself account for the pattern of evolution, but Persat (personal communication) has reported peaks in the numbers of C. nasus nasus in the River Rhone every four or five years that could account for dynamic changes in the hybrid zone (which tended to be confirmed by the distribution of age classes) in the two sampling years. Indeed, the structure of fish communities has great potential for rapid change in response to habitat loss and human activities modifying habitat. This is particularly true of group-spawning species (e.g. Cyprinidae) that focus their reproductive activities on specific spawning substrates, often within very restricted areas.

Our results show that environment may play a key role in determining the directionality of hybridization and the maintenance of particular hybrid combinations and we provide the first evidence of reversible directionality (or more precisely intensity) of gene flow between two species in the same river. The introgression observed was bidirectional, but the direction of mitochondrial introgression was density-dependent, favoring the rare species (figure 6). This “inheritance from the rare species” may be due to sexual selection [21], with all females fertilized but rare males having to compete directly with the males of the common species. Avise and Saunders [22] came close to formulating this hypothesis. However, they did not explain why the same argument did not apply to males and simply assumed differences in the abundances of the sexes in Lepomis species. Our field study (with a sex ratio of about 1) provides the first demonstration of asymmetric mtDNA introgression due to sexual selection. Indeed, we were able to test the hypothesis of “inheritance from the rare species” because we observed an alternation of gene flow between the two species in the same river, and because we knew that hybridization was contemporary, making it possible to “calibrate” each of the two species in allopatric zones. The station-dependence of the proportion of the rare species and of inheritance from the rare species indicates that this phenomenon is dependent on population dynamics.

Asymmetric mtDNA introgression and “inheritance from the rare species”: a ubiquitous pattern

Some studies have reported cases of unequal gene exchange between species [23], [24] or of cytonuclear disequilibrium and mitochondrial capture ([25], and more recently [26]), in which the mtDNA genotypes characteristic of one species occur in a predominant nuclear background from another species. The “capture” phenomenon is observed in both plants and animals, from fish to mice [25]. Most of these studies provisionally eliminated lineage sorting from a polymorphic ancestor as an alternative explanation and documented the widespread occurrence of historical gene exchange between species. They concluded that this capture phenomenon was due either to ancient introgressive hybridization between the related species with the “mother species” outnumbered in the past, or biases in species sampling [e.g. in], [27], ,[28; 29; 30]. We list examples of hybridization in the animal kingdom in Table 3, focusing particularly on the mitochondrial DNA present in the hybrids (similar conclusions could be inferred from Y chromosome data but the fish studied here has no sexual chromosome. We therefore limited our interpretation to mitochondrial DNA and the maternal lineage). In 65% of all reported cases, the hybrids contained the mitochondrial DNA of only one of the two parental species. Not all authors have tried to explain the absence of one of the two possible mtDNA types in the hybrids observed. Others have proposed various “prezygotic” and “postzygotic” hypotheses (cf. “authors conclusion” in Table 3). Rather than focusing on these hypotheses, we considered our own biological model and the characteristics of our hybrid zone. We concluded that the “most extreme case” presented by Wirtz [21] (with some crosses of two given species giving different results in different areas) was probably not an extreme scenario and was probably more common than initially thought and should be considered in other fish and metazoan models. Indeed, recent studies have shown that there are more hybrid taxa than expected in the animal kingdom [31], and comparisons of the results obtained in different studies on different biological models [e.g. 32], [33] may lead to the identification of general processes. For example, flycatcher hybrids are less fit than their parents [34]. However, the presence of these hybrids was accounted for by particular demographic situations in which it was better for females to mate with heterospecific males rather than not to mate at all. Arnold found this conclusion “intriguing” [35]. An “intriguing” feature of our work is that our model fish, which is very different in terms of reproduction and mate choice strategies (group spawner versus strong mating preference and partner choice), seems to display similar general biological behavior. Our “inheritance of rare species” hypothesis may thus be more general, possibly accounting for the bidirectional mtDNA introgression observed by Thulin et al. in 1997 and 2003 [36], [37] in hare species (Table 3). However, more detailed studies are required to confirm this.

Hybridization and life history traits (figures 7 and 8)

thumbnail
Figure 7. LIFE HISTORY TRAITS EVOLUTION THROUGH THE HYBRID ZONE.

Global summary of changes in ecological and morphological traits in the hybrid zone with respect to the reference populations (see the text for more detail). Arrows indicate the overall trend for overall changes in snout orientation, “?” indicates that we could not strictly determine the diet of the Chondrostoma specimens in the hybrid zone, * indicates a significant difference, the other tests could not be done because the sample was too small, but we extrapolated a trend. “ΣD12” indicates the sum of degree-days above 12°C for each station (mean established for an 8-year data set) Chappaz (personal communication).

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

thumbnail
Figure 8. MORPHOLOGY AND GENOME DILUTION

. The graphs show the distribution of individuals for each genome dilution class (x-axis) as a function of i) their probability based on numbers of specimens; y-axis (first graph) ii) their probability based on numbers of combinations; y-axis (second graph). The color gradient (observed values) extends from red (T5) to green (H5). The intermediate morphologies and genomic combinations are shown in yellow-orange, with theoretical values shown in gray (see the text for more detail). The fish morphologies shown are a simplistic representation of the morphological gradient (not all intermediate morphologies are represented).

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

We observed a significant change in the size of individuals, with, in particular, T5 specimens being smaller than the C. t. toxostoma specimens sampled from reference populations. This change in size may have limited competition with the invasive C. n. nasus species. However, the coefficients of condition (species effect in the allopatric zones but no genomic combination effect in the hybrid zone) obtained could be interpreted as showing convergence in the feeding of the two species in the Durance or homogeneous food assimilation. No such convergence was assumed in published data obtained for allopatric zones [38, 39 for C. t. toxostoma in the Rhone and Garonne Basins, respectively; 40 for C. n. nasus in the Rhone Basin]. Moreover, our hypotheses seem to be confirmed by the first results obtained in stable isotope (13C and 15N) studies on 168 Durance Chondrostoma specimens for which we can identify no species, station, or year effect.

T5 individuals reached sexual maturity at a significantly younger age than reference C. t. toxostoma individuals (2–3 years old, rather than 3–4 years old). A similar trend was observed for H5 and C. n. nasus (4–5 years old rather than 5–6 years old), but the sample size was too small for confirmation. Hybrid zone females also had significantly fewer eggs than reference population females (figure 7), consistent with previous findings [41], [42]. Changes in temperature may be responsible for changes in age at sexual maturity and the number of eggs per female [as in 43]. Thus, interspecies competition is not the only factor responsible for ecological displacement; environmental pressure also plays a key role in the evolution of sympatric species.

We observed a morphological gradient, passing from an upper limit of about 80% of alleles from one species, through an intermediate morphology to the morphology of the other species (figure 8). However, we identified changes in the general body shape of hybrid zone individuals, with a shift in mean morphology for all genomic combinations of the hybrid zone, regardless of sampling station or the age or sex of specimens. Surprisingly, this shift was towards the nasus morphology rather than that of toxostoma, which is presumably better adapted to its niche. This shift in morphology may be associated with the Durance River rather than a direct effect of hybridization, because the level of hybridization differed between stations, making it highly unlikely that the direction of change would be the same at all stations. This result tends to confirm the influence of environment on the Chondrostoma complex, as suggested by Peacor et al. [44]. So, is the regular flow of water from upstream to downstream (via the dams) a major factor in the global unidirectional deformation of body shape regardless of genomic combination, and of the size and age of specimens? Studies of other cyprinid hybrid zones would be required to test these hypotheses and to evaluate the role of ecological conditions in speciation [45].

Crossroads between population genomic and ecological approaches

Our results indicate that the interspecies population dynamics of hybrid zones are governed by a genomic component and an ecological component, which may be difficult to dissociate. It is widely thought that shoal members benefit most from being in phenotypically uniform shoals [20]. Could this account for the unidirectional effect of the Durance River on the morphology of the two Chondrostoma species?

A previous study [46] demonstrated that shoaling behavior has a genetic component and provided a preliminary estimate of its heritability in zebrafish (Danio rerio). More recently [47], QTL for shoaling and boldness have been identified on chromosome 16 of Danio rerio. This is particularly interesting because the Tpi 1b gene is also present in this linkage group. This result supports our biological conclusions and confirms the need to focus on these aspects for future “eco-geno” studies of hybrid zones.

Materials and Methods

Collection of samples and area studied

For each of the two species, a population in allopatry (38 individuals from the Garonne for C. t. toxostoma; 33 individuals from the Flet for C. n. nasus) and a population in parapatry (30 individuals from the Doubs for C. t. toxostoma; 27 individuals from the Allier for C. n. nasus) were analyzed and treated as reference populations (i.e. “pure” representative specimens of the two species; figure 9). The choice of these reference populations was dictated by a previous molecular study [9], in which these populations displayed the largest range of variability. The River Durance (hybrid zone) is the last major affluent of the Rhone on the eastern side. This river has many dams and weirs isolating different parts of the river and preventing the migration of freshwater fauna. Construction of the major dams began around 1950, after C. n. nasus began to colonize the Durance.

thumbnail
Figure 9. GEOGRAPHICAL AREA STUDIED.

Allopatric zone: Garonne (C. t. toxostoma); Flet (C. n. nasus); parapatric zones: Doubs (C. t. toxostoma); Allier (C. n. nasus) and hybrid zone: the River Durance and the four sampling stations: Buech, Manosque, Pertuis and Cavaillon.

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

In the hybrid zone, 642 individuals were sampled in 2001 and 853 in 2002, at four independent stations (figure 9). Each station was sampled repeatedly between March and the end of June (Appendix S1). An identifier was attributed to each fish, in order of sampling.

Molecular marker

Total genomic DNA was extracted from muscle tissue or from a piece of caudal fin, using a standard phenol chloroform protocol [48].

mtDNA.

A 610 bp fragment of the cytochrome b mtDNA sequence was amplified using primers and an amplification protocol described elsewhere [10]. The single strand was used as a template for the sequencing of 532 bp with an automated DNA sequencer (Macrogen Inc.) (Genbank accession number EF363338-EF363452).

Nuclear introns.

We first tested 13 introns (Gilles and Costedoat unpublished data). Each was used to search the Genbank database for potential paralogs representing a major risk in our study of hybridization between species of the Cyprinidae family, due to the presence of 3R replication events in teleosts in particular [49], [50], [51]. The position of the gene was systematically checked, based on the chromosome map of Danio rerio (Teleostei; Cyprinidae). We finally chose to study four introns (ribosomal protein gene S7 intron 1; triose phosphate isomerase 1b, Tpi1b intron 4; α-tropomyosin, α-Trop intron 5; and recombination activating gene 1, RAG 1 intron 2) because they presented a sufficiently high level of polymorphism and alleles capable of discriminating between the two Chondrostoma species. Sequence variation within and between populations and species was characterized by analyzing single-stranded conformational polymorphisms (SSCP) [52], [53], [54], using [α32P]-dATP (IRSN Corp.)-labeled DNA fragments. The amplification protocols and primers are presented in the Dataset S3. As the SSCP method requires the migration of relatively small fragments, the sequence lengths of the intron amplicons were 392 bp for RAG-1 (Genbank accession number EF363453-EF363458), 186 bp for S7 (Genbank accession number EF363459-EF363466), 327 bp for Tpi1b (Genbank accession number EF363467-EF363472) and 224 bp for α-Trop (Genbank accession number EF363473-EF363479).

Mitochondrial and nuclear allele affiliations and phylogenetic reconstructions are presented in the Dataset S4.

Specimen identification and Hardy-Weinberg equilibrium

Each specimen of the hybrid zone was identified by combining the alleles observed for each of the five molecular markers (cf.Dataset S4). In “pure” C. n. nasus specimens, all five markers (cytochrome b plus the four introns) originated from C. n. nasus (the reciprocal applied to “pure” C. t. toxostoma individuals). We could potentially identify 162 genomic combinations (the two “pure” or parental combinations plus 160 hybrid or introgressed combinations). We used a previously described approach ([55] and [56]) to investigate non randomness in genotype frequencies (departure from Hardy-Weinberg expectations and from linkage equilibrium). We used chi-square tests to check for non randomness in genotype frequencies, with Benjamini-Hochberg correction for multiple comparisons. Indeed, as we carried out as many tests as there were combinations for the Hardy-Weinberg test of equilibrium, (i.e. 162), we had to control for false positives. We used the false discovery rate (FDR) and its majoration proposed by Benjamini and Hochberg [57]. Let (Pi)i be the ordered P-value for the 162 tests,3Then we have FDR (t) ≤ α

In practical terms, each test giving a P-value Pi such that was declared significant.

With this approach, we were able to assign the different genomic combinations to three major groups: combinations that were overrepresented, combinations that were underrepresented and combinations present in the expected proportions. This approach also permitted robust use of the statistical models.

Morphological markers

We ensured that the morphological ranges of the two species were as large as possible by also adding two other allopatric populations—32 C. t. toxostoma specimens from the Berre (a Mediterranean coastal river) and 24 C. n. nasus specimens from Romania.

Morphology was analyzed by studying four meristic characters (numbers of dorsal, anal and pelvic fins, rays and number of scales along the lateral line) and through a morphometric approach, based on the landmarks method [58]. We defined 21 homologous landmarks on the body (figure S5). All landmarks were digitized, using TpsDig software [59]. Landmark-based geometric morphometric methods were used to capture information about shape, by obtaining the x and y coordinates of homologous landmarks. Differences in the sets of coordinates between specimens due to scaling, rotation and translation were eliminated by a typical geometric morphometric approach [60], [61] in which the specimens were placed in a procruste superimposition on the iteratively estimated mean reference form, using the generalized procruste analysis (GPA) procedure. Points representing landmark configurations were then projected into Euclidean tangent space approximating the curved shape space.

Morphological variables (both meristic and plastic variables) were first analyzed for the allopatric and parapatric populations, and then for the hybrid zone. For each analysis, a MANCOVA was performed on size, age, sex, station and species (or genetic class) effects. Size was calculated as centroid size (i.e. the squared root of the sum of squared distances between each landmark and the centroid (centre of gravity) of a landmark configuration). For significant effects, discriminant analysis was then performed, with deformation grids used for plastic variables, using Bookstein's [62] thin plate spline algorithm. All statistical analyses were performed using R [63]; deformation grids were visualized with the tpsgrid function of the R shape package [60].

Life history traits

Weight, size, age and sex were determined for each specimen.

The relative robustness, or degree of well-being, can be expressed as the “coefficient of condition” [64]. Variations in the coefficient of condition primarily reflect the state of sexual maturity and nourishment. The formula most often used is: K = W/L3 (W = the weight in grams; L = the standard length in decimeters). We carried out an analysis of variance to determine whether the mean coefficient of condition differed significantly i) between the two Chondrostoma species ii) between the different populations of each of the Chondrostoma species (allopatry and parapatry) and iii) between the reference populations and the hybrid zone populations.

As both species are protected, the individuals sampled in 2002 were analyzed directly in the field and then released. However, the Chondrostoma species do not display sexual dimorphism and have no sexual chromosome. Sex was therefore determined by the canulation method for this sampling year (Stolzenberg personal communication). As a result, it was very difficult to determine the sex of the animals in some sampling campaigns and for some specimens.

Variable interactions

We analyzed interactions between the various qualitative variables, using chi-square tests or log-linear models [10]. For quantitative variables (coefficient of condition and morphology), interactions were analyzed by variance-based methods—ANOVA and MANCOVA, respectively.

Supporting Information

Appendix S1.

SAMPLING INFORMATION. Sampling campaigns and sample size for each Durance station.

https://doi.org/10.1371/journal.pone.0000357.s001

(0.16 MB DOC)

Dataset S1.

Molecular marker history and allele distribution.

https://doi.org/10.1371/journal.pone.0000357.s002

(0.06 MB DOC)

Dataset S2.

Hybrid zone and life history traits

https://doi.org/10.1371/journal.pone.0000357.s003

(0.04 MB DOC)

Dataset S3.

Nuclear amplification protocols and SSCP method

https://doi.org/10.1371/journal.pone.0000357.s004

(0.04 MB DOC)

Figure S1.

MITOCHONDRIAL C.T. TOXOSTOMA NETWORK. C.t. toxostoma mtDNA haplotype network.

https://doi.org/10.1371/journal.pone.0000357.s006

(2.44 MB TIF)

Figure S2.

MITOCHONDRIAL C.N. NASUS NETWORK. C.n. nasus mtDNA haplotype network.

https://doi.org/10.1371/journal.pone.0000357.s007

(1.12 MB TIF)

Figure S3.

PHYLOGENETIC RECONSTRUCTION. Neighbor-Joining (with Kimura-2-parameter distance model) based on the sequences obtained for each of the four introns. All the topologies were tested by non-parametric bootstrap resampling (500 replicates).

https://doi.org/10.1371/journal.pone.0000357.s008

(0.86 MB TIF)

Figure S4.

EVOLUTIONARY HISTORY OF THE NUCLEAR INTRONS. Left side: summary of the phylogentic relationships between haplotypes of the two species (C. nasus and C. toxostoma) obtained with NJ and parsimony methods. Right side: haplotype relationships under haplotype network representation. C. toxostoma haplotypes are in red, C. nasus haplotypes are in green. For network reconstruction, the relationships between the two species haplotypes required a decrease in the parsimony limit (90%) obtained with TCS software. The arrow corresponds to the root as observed on the NJ trees. For the S7 intron, the Rh2 haplotype corresponds to a Romanian C. n. nasus haplotype.

https://doi.org/10.1371/journal.pone.0000357.s009

(1.31 MB TIF)

Figure S5.

LANDMARKS. The 21 chosen landmarks used for morphometric analysis.

https://doi.org/10.1371/journal.pone.0000357.s010

(0.51 MB TIF)

Figure S6.

DISCRIMINANT ANALYSIS REALIZED ON PLASTIC MORPHOLOGICAL CHARACTERS. Illustration of the Durance tendency. T = C.t.toxostoma in allopatry/paraptry; H = C.n.nasus in allopatry/paraptry; H5 = C.n.nasus in the hybrid zone; T5 = C.t.toxostoma in the hybrid zone and the different hybrid groups. The DA was carried out with all the groups, but we illustrate only the variance of the T (in blue), T5 (in red), H (in green) and H5 (purple) groups.

https://doi.org/10.1371/journal.pone.0000357.s011

(1.19 MB TIF)

Figure S7.

VARIANCE OF THE COEFFICIENT OF CONDITION. Graphic representation of variance analyses of the coefficient of condition in sympatry. A) by year; B) as a function of sex F = female; M = male; I = immature, C) by station, 17 = Buech; 18 = Manosque; 19 = Pertuis; 20 = Cavaillon, D) by age and E) by genomic combination, H5; H4T; Hi4; T4i; T5; -1 = under-represented combinations; 0 = combinations in the expected proportions; 99 = over-represented combinations; *significant value.

https://doi.org/10.1371/journal.pone.0000357.s012

(1.59 MB TIF)

Table S1.

MTDNA HAPLOTYPE DISTRIBUTION. Letter A = C. toxostoma; B = C. nasus and y = hybrid zone haplotypes. Haplotype diversity (H) and nucleotide diversity (Pi) were calculated for each population.

https://doi.org/10.1371/journal.pone.0000357.s013

(0.09 MB DOC)

Table S2.

NUCLEAR ALLELES DISTRIBUTION. Genotype distribution of the 4 introns (S7; &#x03B1;-Trop; Rag; Tpi) as a function of the population sampled (reference populations and hybrid zone populations).

https://doi.org/10.1371/journal.pone.0000357.s014

(0.14 MB DOC)

Table S3.

CROSS VALIDATION PROCEDURE FOR MORPHOLOGICAL CHARACTERS. C.n.n = C.n.nasus in reference zones; C.t.t = C.t. toxosotoma in reference zones. -1 = underrepresented combinations; 0 = combinations in expected proportion; 1 = overrepresented combinations; H5 = C.n.nasus from hybrid zone; T5 = C.t. toxosotoma from hybrid zone; H4T; Hi4; T4i = hybrid combinations (cf text for more explanations).

https://doi.org/10.1371/journal.pone.0000357.s015

(0.11 MB DOC)

Table S4.

MANCOVA ON MORPHOLOGICAL DATA SET. Order-3 analysis. Cond: Coefficient of condition, Gen. Group: Genomic group (cf. text for more details).

https://doi.org/10.1371/journal.pone.0000357.s016

(0.13 MB DOC)

Table S5.

ANOVA ON COEFFICIENT OF CONDITION DATA SET.

https://doi.org/10.1371/journal.pone.0000357.s017

(0.03 MB DOC)

Table S6.

AGE DISTRIBUTION. Distribution of hybrid zone combinations as a function of genome dilution and age (data are in percentages).

https://doi.org/10.1371/journal.pone.0000357.s018

(0.06 MB DOC)

Table S7.

SEX RATIO IN THE HYBRID ZONE. References did not indicate a sex ratio significantly different from 1 for the two species studied in allopatric populations [1], [2]. F = female, M = male and I = indeterminate. 1) Gozlan RE (1998) Environmental biology and morphodynamics of the sofie Chondrostoma toxostoma (Cyprinidae), with emphasis on early development. PhD, Université Toulouse, France. 196p. 2) Nelva-Pasqual A (1985) Biogéographie, démographie et écologie de Chondrostoma nasus nasus (L., 1758) Hotu (Poisson, Téléostéen, Cyprinidé). PhD Université Claude-Bernard-Lyon I. 340p.

https://doi.org/10.1371/journal.pone.0000357.s019

(1.51 MB DOC)

Table S8.

VARIABLES INTERACTIONS. A) order-two interactions (chi2); B) order-three interactions (log-linear model) and main associated results.

https://doi.org/10.1371/journal.pone.0000357.s020

(0.12 MB DOC)

Acknowledgments

We would like to thank the members of the “Hydrobiology” laboratory from their help with field work. We particularly thank Nicolas Stolzenberg, Adrian Ionascu and Laurent Cavalli for determining the age and fecundity of specimens, and Jonathan Grey for the stable isotope results. We thank the officers of the “Conseil Supérieur de la pêche” for their help and interest in our work. Our English usage was corrected by Alex Edelman and Associates.

Author Contributions

Conceived and designed the experiments: AG. Performed the experiments: CC. Analyzed the data: AG CC NP. Contributed reagents/materials/analysis tools: AG CC NP RC. Wrote the paper: CC. Other: Manuscript revision and field work: AG. Field work: CC. Manuscript revision and field work: NP. Ecological and financial support, field work: RC.

References

  1. 1. Rieseberg LH (1998) Molecular ecology of hybridization. In: Carvalho GR, editor. Advances in Molecular Ecology. IOS press. pp. 243–265.
  2. 2. Rieseberg LH, Gerber D (1995) Hybridization in a Catalina mahogany: RAPD evidence. Cons Biol 9: 199–203.
  3. 3. Rhymer JM, Simberloff D (1996) Extinction by hybridization and introgression. Annu Rev Syst 27: 83–109.
  4. 4. Rieseberg LH (1991) Homoploid reticulate evolution in Helianthus: evidence from ribosomal genes. Am J Bot 78: 1218–1237.
  5. 5. Ellstrand NC, Elam DR (1993) Population genetic consequences of small population size: implications for plant conservation. Annu Rev Ecol Syst 24: 217–242.
  6. 6. Dowling TE, Secor C (1997) The role of hybridization and introgression in the diversification of animals. Annu Rev Ecol Syst 28: 593–619.
  7. 7. Wilson EO (1992) The Diversity of Life. Massachusetts: Harvard University Press.
  8. 8. Scribner , Page KS, Bartron ML (2001) Hybridization in freshwater fishes: a review of case studies and cytonuclear methods of biological inference. Fish Biol&Fish 10: 293–323.
  9. 9. Costedoat C, Pech N, Chappaz R, Salducci D, Lim P, Gilles A (2004) Study of introgressive hybridization between Chondrostoma t. toxostoma and Chondrostoma n. nasus (Teleostei, Cyprinidae) using multiple approaches. Cybium 28: 51–61.
  10. 10. Costedoat C, Pech N, Salducci MD, Chappaz R, Gilles A (2005) Evolution of mosaic hybrid zone between invasive and endemic species of Cyprinidae through space and time. Biol J Linn Soc 85: 135–155.
  11. 11. Crespin L, Berrebi P, Lebreton JD (2002) Spatially varying natural selection in a fish hybrid zone. J Fish Biol 3: 696–711.
  12. 12. Dowling TE, Broughton RE, DeMarais BD (1997) Significant role for historical effects in the evolution of reproductive isolation: Evidence from patterns of introgression between the cyprinid fishes, Luxilus cornutus and Luxilus chrysocephalus. Evolution 51: 1574–1583.
  13. 13. Dayan T, Simberloff D (2005) Ecological and community-wide character displacement: the next generation. Ecol Lett 8: 875–894.
  14. 14. Valentin A, Sevigny JM, Chanut JP (2002) Geometric morphometrics reveals body shape differences between sympatric redfish Sebastes mentella, Sebastes fasciatus and their hybrids in the Gulf of St Lawrence. J Fish Biol 60: 857–875.
  15. 15. Arnold ML (1997) Natural hybridization and evolution. Oxford University press.
  16. 16. Das J (2006) The role of mitochondrial respiration in physiological and evolutionary adaptation. Bioessays 28: 890–901.
  17. 17. Willett CS (2006) Deleterious epistatic interactions between electron transport system protein-coding loci in the copepod Tigriopus californicus. Genetics 173: 1465–1477.
  18. 18. Pritchar VL, Lawrence J, Butlin RK, Krause J (2001) Shoal choice in zebrafish, Danio rerio: the influence of shoal size and activity. Anim Behav 62: 1085–1088.
  19. 19. Ward AJW, Krause J (2001) Body length assortative shoaling in the European minnow, Phoxinus phoxinus. Anim Behav 62: 617–621.
  20. 20. Ward AJW, Axford S, Krause J (2002) Mixed-species shoaling in fish: the sensory mechanisms and costs of shoal choice. Behav Ecol Sociobiol 52: 182–187.
  21. 21. Wirtz P (1999) Mother species-father species: unidirectional hybridization in animals with female choice. Anim Behav 58: 1–12.
  22. 22. Avise JC, Saunders NC (1984) Hybridization and introgression among species of sunfish (Lepomis): analysis by mitochondrial DNA and allozyme markers. Genetics, 108: 237–255.
  23. 23. McDonald DB, Clay RP, Brumfield RT, Braun MJ (2001) Sexual selection on plumage and behavior in an avian hybrid zone: Experimental tests of male-male interactions. Evolution 57: 1443–1451.
  24. 24. Thulin CG, Tegelstrom H (2002) Biased geographical distribution of mitochondrial DNA that passed the species barrier from mountain hares to brown hares (genus Lepus): an effect of genetic incompatibility and mating behaviour? J Zool 258: 299–306.
  25. 25. Harrison RG (1989) Animal mitochondrial-DNA as a genetic–marker in population and evolutionary biology. Trends Ecol Evol 4: 6–11.
  26. 26. Avise JC (2004) Molecular markers, Natural History and Evolution, 2 edn. New York: Chapman & Hall..
  27. 27. Dowling TE, Smith GR, Brown WM (1989) Reproductive isolation and introgression between Notropis cornutus and Notropis chrysocephalus (Family Cyprinidae): comparison of morphology, allozymes and mitochondrial DNA. Evolution 43: 620–634.
  28. 28. Taylor DJ, Hebert PDN (1993) Cryptic intercontinental hybridization in Daphnia (Crustacea)-the ghost of introductions past. Proc Roy Soc London III 254: 163–168.
  29. 29. Wayne RK (1993) Molecular evolution of the dog family. Trends Genet 9: 218–224.
  30. 30. Carson EW, Dowling TE (2006) Influence of hydrogeographic history and hybridization on the distribution of genetic variation in the pupfishes Cyprinodon atrorus and C. bifasciatus. Mol Ecol 15: 667–679.
  31. 31. Mallet J (2005) Hybridization as an invasion of the genome. Trends Ecol Evol 20: 229–237.
  32. 32. Joron M, Papa R, Beltran M, Chamberlain N, Mavarez J, et al. (2006) A conserved supergene locus controls colour pattern diversity in Heliconius butterlies. Plos Biol 4: 1831–1840.
  33. 33. Terai Y, Seehausen O, Sasaki T, Takahashi K, Mizoiri S, et al. (2006) Divergent selection on opsins drives incipient speciation in Lake Victoria cichlids. Plos Biol 4: 2244–2251.
  34. 34. Veen T, Borge T, Griffith SC, Saetre GP, Bures S, et al. (2001) Hybridization and adaptive mate choice in flycatchers. Nature 411: 45–50.
  35. 35. Arnold ML (2006) Evolution through genetic exchange. Oxford University Press..
  36. 36. Thulin CG, Jaarola M, Tegelstrom H (1997) The occurrence of mountain hare mitochondrial DNA in wild brown hares. Mol Ecol 6: 463–467.
  37. 37. Thulin CG, Tegelstrom H, Fredga K (2003) Haplotype-diversity of mountain hare mtDNA among native mountain-hares and introduced brown hares in Scandinavia. Ann Zool Fenn 40: 45–52.
  38. 38. Chappaz R, Brun G, Olivari G (1989) New data about biology and ecology of a little studied fish–Chondrostoma toxostoma (Vallot, 1836) comparison with nase, Chondrostoma nasus (L 1766). CR Acad Sciences series III 5: 181–186.
  39. 39. Gozlan RE (1998) Environmental biology and morphodynamics of the sofie Chondrostoma toxostoma (Cyprinidae), with emphasis on early development. France196: PhD, Université Toulouse,.
  40. 40. Nelva-Pasqual A (1985) Biogéographie, démographie et écologie de Chondrostoma nasus nasus (L., 1758) Hotu (Poisson, Téléostéen, Cyprinidé). PhD Université Claude-Bernard-Lyon I.
  41. 41. Angilletta MJ, Wilson RS, Navas CA, James RS (2003) Tradeoffs and the evolution of thermal reaction norms. Trends Ecol Evol 18: 234–240.
  42. 42. Angilletta MJ, Steury TD, Sears MW (2004) Temperature, growth rate, and body size in ecotherms: fitting pieces of a life-history puzzle. Integr Comp Biol 44: 498–509.
  43. 43. Sandstrom O, Abrahamsson I, Andersson J, Vetemaa M (1997) Temperature effects on spawning and egg development in Eurasian perch. J Fish Biol 51: 1015–1024.
  44. 44. Peacor SD, Allesina S, Riolo RL, Pascual M (2006) Phenotypic plasticity opposes species invasions by altering fitness surface. Plos Biol 4: 2112–2120.
  45. 45. McKinnon JS, Mori S, Blackman BK, David L, Kingsley DL, et al. (2004) Evidence for ecology's role in speciation. Nature 429: 294–298.
  46. 46. Wright D, Rimmer LB, Pritchard VL, Krause J, Butlin RK (2003) Inter and intra-population variation in shoaling and boldness in the zebrafish (Danio rerio). Naturwissenschaften 90: 374–377.
  47. 47. Wright D, Nakamichi R, Krause J, Butlin RK (2006) QTL analysis of behavioral and morphological differentiation between wild and laboratory zebrafish (Danio rerio). Behav Genet 36: 271–284.
  48. 48. Sambrook J, Fritsch EF, Maniatis T (1989) Molecular cloning: A Laboratory Manual, 2dn edition. New York: Cold Spring Harbor Laboratory Press..
  49. 49. Vandepoele K, De Vos W, Taylor JS, Meyer A, Van de Peer Y (2004) Major events in the genome evolution of vertebrates: Paranome age and size differ considerably between ray-finned fishes and land vertebrates. Proc Natl Acad Sci USA 101: 1638–1643.
  50. 50. Meyer A, Van de Peer Y (2005) From 2R to 3R: evidence for a fish-specific genome duplication (FSGD). Bioessays 27: 937–945.
  51. 51. Steinke D, Hoegg S, Brinkmann H, Meyer A (2006) Three rounds (1R/2R/3R) of genome duplications and the evolution of the glycolytic pathway in vertebrates. BMC Biol 4: 1–7.
  52. 52. Hayashi K (1991) PCR-SSCP: a simple and sensitive method for detection of mutations in the genomic DNA. PCR Methods Appl 1: 34–38.
  53. 53. Orti G, Hare MP, Avise JC (1997) Detection and isolation of nuclear haplotypes by PCR-SSCP. Mol Ecol 6: 575–580.
  54. 54. Sunnucks P, Wilson AC, Beheregaray LB, et al. (2000) SSCP is not so difficult: the application and utility of single-stranded conformation polymorphism in evolutionary biology and molecular ecology. Mol Ecol 9: 1699–1710.
  55. 55. Orive ME, Barton NH (2002) Associations between cytoplasmic and nuclear loci in hybridizing populations. Genetics 162: 1469–1485.
  56. 56. Taylor EB, Boughman JW, Groenenboom M, Sniatynsky M, Schluter D, et al. (2006) Speciation in reverse: morphological and genetic evidence of the collapse of a three-spined stickleback (Gasterosteus aculeatus) species pair. Mol Ecol 2: 343–355.
  57. 57. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate-a practical and powerful approach to multiple testing. J Roy Stat Soc B 57: 289–300.
  58. 58. Adams DC, Rohlf FJ, Slice DE (2004) Geometric morphometrics: ten years of progress following the ‘revolution’. It J Zool 71: 5–16.
  59. 59. Rohlf FJ (2004) tpsDig, version 1.40. Department of Ecology and Evolution, State University of New York at Stony Brook.
  60. 60. Dryden I, Mardia K (1998) Statistical shape analysis. New-York: John Wiley&Sons, USA..
  61. 61. Rohlf FJ (1999) Shape statistics: procustes superimpositions and tangent spaces. J. Class. 16: 197–223.
  62. 62. Bookstein FL (1991) Morphometric tools for landmark data: geometry and biology. Cambridge: Cambridge University Press..
  63. 63. R Development Core Team (2005) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, URL http://www.R-project.org.
  64. 64. Gray MA, Curry AR, Munkittrick KR (2002) Non lethal sampling methods for assessing environmental impacts using small bodied sentinel fish species. Water Qual. Res. J. Canada. 37: 195–211.
  65. 65. Lamb T, Avise JC (1986) Directional introgression of mitochondrial-DNA in a hybrid population of tree frogs-the influence of mating-behavior. Proc Natl Acad Sci USA 83: 2526–2530.
  66. 66. Spolsky C, Uzzell T (1986) Evolutionary history of the hybridogenetic hybrid frog Rana esculenta as deduced from mtDNA analyses. Mol Biol Evol 3: 44–56.
  67. 67. Gill FB (1997) Local cytonuclear extinction of the golden-winged warbler. Evolution 51: 519–525.
  68. 68. Andersson (1999) Hybridization and skua phylogeny. Proc Roy Soc London III 266: 1579–1585.
  69. 69. Edwards CA, Skibinski DOF (1987) Genetic-variation of mitochondrial-DNA in mussel (Mytilus edulis and Mytilus galloprovincialis) populations from south-west England and South Wales. Marine Biol 94: 547–556.
  70. 70. Glemet H, Blier P, Bernatchez L (1998) Geographical extent of Arctic char (Salvelinus alpinus) mtDNA introgression in brook char populations (S. fontinalis) from eastern Quebec, Canada. Mol Ecol 7: 1655–1662.
  71. 71. Wilson CC, Bernatchez L (1998) The ghost of hybrids past: fixation of arctic charr (Salvelinus alpinus) mitochondrial DNA in an introgressed population of lake trout (S. namaycush). Mol Ecol 7: 127–132.
  72. 72. Avise JC, Pierce PC, VandenAvyle MJ, Smith MH, Nelson WS, Asmussen MA (1997) Cytonuclear introgressive swamping and species turnover of bass after an introduction. J Hered 88: 14–20.
  73. 73. Dowling TE, Smith GR, Brown WM (1989) Reproductive isolation and introgression between Notropis cornutus and Notropis chrysocephalus (family Cyprinidae)-comparison of morphology, allozymes, and mitochondrial DNA. Evolution 43: 620–634.
  74. 74. Rognon X, Guyomard R (2003) Large extent of mitochondrial DNA transfer from Oreochromis aureus to O. niloticus in West Africa. Mol Ecol 12: 435–445.
  75. 75. Jiggins CD, McMillan WO, King P, Mallet J (1997) The maintenance of species differences across a Heliconius hybrid zone. Heredity 79: 495–505.
  76. 76. Harrison RG, Rand DM, Wheeler WC (1987) Mitochondrial-DNA variation in field crickets across a narrow hybrid zone. Mol Biol Evol 4: 144–158.
  77. 77. Desalle R, Giddings LV, Kaneshiro KY (1986) Mitochondrial-DNA variability in natural-populations of Hawaiian Drosophila. 2. Genetic and phylogenetic relationships of natural populations of Drosophila silvestris and Drosophila heteroneura. Heredity 56: 87–96.
  78. 78. Ferris SD, Sage RD, Huang CM, Nielsen JT, Ritte U, Wilson AC (1983) Flow of mitochondrial DNA across a species boundary. Proc Natl Acad Sci USA 80: 2290–2294.
  79. 79. Good JM, Demboski JR, Nagorsen DW, Sullivan J (2003) Phylogeography and introgressive hybridization: Chipmunks (genus Tamias) in the northern Rocky Mountains. Evolution 57: 1900–1916.