Restricted gene flow at the micro- and macro-geographical scale in marble trout based on mtDNA and microsatellite polymorphism
© Pujolar et al; licensee BioMed Central Ltd. 2011
Received: 12 November 2010
Accepted: 14 April 2011
Published: 14 April 2011
The genetic structure of the marble trout Salmo trutta marmoratus, an endemic salmonid of northern Italy and the Balkan peninsula, was explored at the macro- and micro-scale level using a combination of mitochondrial DNA (mtDNA) and microsatellite data.
Sequence variation in the mitochondrial control region showed the presence of nonindigenous haplotypes indicative of introgression from brown trout into marble trout. This was confirmed using microsatellite markers, which showed a higher introgression at nuclear level. Microsatellite loci revealed a strong genetic differentiation across the geographical range of marble trout, which suggests restricted gene flow both at the micro-geographic (within rivers) and macro-geographic (among river systems) scale. A pattern of Isolation-by-Distance was found, in which genetic samples were correlated with hydrographic distances. A general West-to-East partition of the microsatellite polymorphism was observed, which was supported by the geographic distribution of mitochondrial haplotypes.
While introgression at both mitochondrial and nuclear level is unlikely to result from natural migration and might be the consequence of current restocking practices, the pattern of genetic substructuring found at microsatellites has been likely shaped by historical colonization patterns determined by the geological evolution of the hydrographic networks.
Freshwater fish species show a greater average degree of genetic differentiation among locations than marine species, often resulting from the isolation of fish populations among drainages . The marked zoogeography produced by the historical patterns of isolation among drainages is subjected to the continuous remodeling of the river drainage and to climatic fluctuations of which the glacial/interglacial periods of the Pleistocene played a crucial role . Species respond actively to fluctuations in their natural range, with hydrographic networks being used for colonization and for retreat throughout geological times. In the case of salmonids, anadromy should increase gene flow among river basins . However, Bernatchez et al.  reported significant heterogeneity at the mitochondrial DNA level among populations of the pan-European brown trout Salmo trutta, which lead to the definition of five distinct evolutionary lineages: four of them, the Adriatic (AD), the Atlantic (AT), the Danube (DA) and the Mediterranean (ME) lineage correspond to the major drainage basins, whereas the fifth lineage is represented by the marble trout (MA). Examination of the evolutionary relationships among the five trout lineages suggested a scenario in which the most ancient phylogeographic separation resulted from allopatric isolation in the Pleistocene between the three major European drainages, Atlantic (AT lineage), Danube (DA lineage) and Mediterranean (AD, MA and ME lineages), followed by later divergence within the Mediterranean basin linked to the existence of glacial refugia that originated the AD, MA and ME lineages . More recently, a sixth lineage has been described, restricted to the Duero river in the Iberian Peninsula . Complex genetic structures exist in many trout species, resulting from multiple colonization episodes by distinct lineages and secondary intergradiations, with restricted gene flow and extreme genetic drift producing deep genetic divergences [7–11].
The marble trout is characterized by a particular marbled colour pattern and different biological and ecological traits such as a very large body size (up to 1 m in length and 25 kg in weight) and the complete lack of black spots on the lateral side of the body. It has a restricted native geographic distribution within the Adriatic basin, found on the alpine part of the river Po in Northern Italy, former Yugoslavia and Albania . The taxonomic status of marble trout is controversial, being recognized either as a separate species Salmo marmoratus[13–16] or as a sub-species Salmo trutta marmoratus within the Salmo trutta complex [17, 18].
The existence of marble trout has been seriously threatened by the introduction of nonindigenous brown trout throughout its natural range by means of massive restocking practices . Hybridization between marble and brown trout has been so extensive that hybrids now dominate most rivers, and typical marble trout forms are rare and limited to secluded headwater streams. Molecular data confirm a high level of introgression in the western part of the Po river in Northern Italy [17, 18], the Soca river system in Slovenia [14–16], and recently, the Adige river system in South Tyrol [20, 21], although all studies also reported pure populations of marble trout in headwaters of all river systems.
Up to now, molecular genetic studies on marble trout have focused on species identification using mitochondrial DNA haplotypes and assessing the rate of introgression among lineages in the Salmo trutta complex, while examination of the population genetic structure of marble trout has been mostly neglected. Using allozymes, Giuffra et al.  found no apparent differences among marble trout samples collected along tributaries of the Po river system in Northern Italy, which contrasted with the fixed or nearly fixed differences for alternate alleles at 8 loci found between marble and brown trout (Nei's genetic distance of 0.16-0.18). However, extreme genetic heterogeneity has been found among remnant marble trout populations from the Soca river system in Slovenia using microsatellites, with a global FST = 0.66 (p < 0.001) and all pairwise FST ranging from 0.31 to 0.88 .
Frequencies of marble and brown trout mitochondrial haplotypes across samples.
>10% Brown trout
Introgression rate (%)
Summary of genetic variability at microsatellite loci across marble trout samples including number of individuals (N), observed (H o ) and expected (H e ) heterozygosities and allelic richness (AR) considering a) all individuals and b) only non-introgressed individuals.
In order to estimate the native patterns of population genetic structure and gene diversity, data was re-analyzed removing all individuals with >10% expected introgression. As a consequence, Adda (reduced to only two individuals) and Koritnica (all individuals introgressed) were excluded from the analysis. The effect of removing introgressed individuals was a expected drop in genetic diversity (Table 2). Heterozygosities were similar across samples except for Predelica, which showed a 2-4-fold lower heterozygosity (He = 0.080) in comparison with the rest of samples (He = 0.164-0.388). Following Bonferroni correction, only one out of 45 tests departed significantly from HWE (locus Str85 at Piave). The software MICROCHECKER showed no evidence for scoring errors due to stuttering or large alle dropout. An homozygote excess was found at locus Str85, which might be attributable to the presence of null alleles, at Adige, Isonzo and Piave, but not at the rest of samples. No linkage disequilibrium was observed between any pair of loci after Bonferroni correction.
Since our study was carried out with a limited number of microsatellite markers, statistical power to detect genetic divergence was assessed with POWSIM prior to the genetic differentiation analysis. Simulations using our empirical microsatellite data, taking into account the sample size for each population, and a wide range of predefined FST values (0.001, 0.005, 0.01, 0.05, 0.1) showed that our markers have enough statistical power to detect FST values > 0.005 (p = 1.000). This suggests that our markers can provide a reasonably accurate picture of the genetic population structure of marble trout.
Pairwise FST estimates (above diagonal) and genetic distances (DCE, below diagonal) using Cavalli-Sforza and Edwards (1967) among all locations (labelled as in Figure 1).
Genetic differentiation (FST) using Hierarchical AMOVA includes FCT (among groups) and FSC (within groups).
2 (Western/Central-Italy; Eastern-Italy/Slovenia)
0.244 (p = 0.000)***
0.013 (p = 0.430)
0.231 (p = 0.000)***
3 (Western/Central-Italy; Eastern-Italy; Slovenia)
0.237 (p = 0.000)***
0.003 (p = 0.321)
0.234 (p = 0.000)***
3 (Western-Italy; Central Italy; Eastern Italy/Slovenia)
0.231 (p = 0.000)***
0.078 (p = 0.053)
0.167 (p = 0.000)***
4 (Western-Italy; Central Italy; Eastern Italy; Slovenia)
0.252 (p = 0.000)***
0.098 (p = 0.025)*
0.171 (p = 0.000)***
Results of the bottleneck analyses using BOTTLENECK (Piry et al. 1999) and M_P_VAL (Garza and Williamson, 2001), and estimates of effective population size (Ne) using DIY-ABC (Cornuet et al. 2008) including 95% confidence interval.
196 - 10700
341 - 16100
232 - 10600
604 - 25400
156 - 9360
538 - 28600
308 - 13900
184 - 7600
131 - 2790
Introgression at the mitochondrial and microsatellite level
Mitochondrial data confirmed the existence of genetic introgression from nonindigenous brown trout into marble trout, suggested by previous mitochondrial DNA polymorphism studies [17, 20], likely consequence of restocking practices. Presence of non-marble trout haplotypes was highly variable across samples, and included (i) pure marble trout populations (Pellice, Predelica), (ii) low introgression of one or two individuals per population in the North-central locations (Prissiano, Adige, Brenta, Piave and Tagliamento), and (iii) moderate to high introgression (from 11% to 55%) in the remaining populations (Ticino, Adda, Isonzo, Koritnica). By comparison, introgression at nuclear (microsatellite) level was much higher and occurred in all locations except Predelica, reaching extreme levels of introgression in Adda and especially Koritnica, where all individuals were introgressed. A general correspondence was observed between mitochondrial-nuclear introgression, so that locations with no mitochondrial introgression showed <3.5% nuclear introgression, locations with a few introgressed individuals at the mitochondrial introgression showed moderate nuclear introgression (up to 20%), while locations with high introgression at the mitochondrial level (>45%) showed all or practically all individuals introgressed at the nuclear level.
When comparing results with other studies using mitochondrial data, the study of Giuffra et al.  with partial sampling overlap with our study (Pellice, Ticino, Brenta) reported a smaller amount of introgression after observing 6 out of 8 pure marble trout populations, which might be due to the lower sampling sizes (N = 8-12). While both studies agree with populations from the western side of the Po river at Pellice being pure marble trout, our study suggested a smaller introgression (7.5%) at Brenta in comparison with the 25% reported in the study of Giuffra et al. . Furthermore, the latter study found no signatures of introgression at Ticino, while our study found 30% exogenous haplotypes. In addition, our sampling overlaps partially with the study of Meraner et al.  focusing on the Adige region in South Tyrol (Northern Italy), which reported a higher introgression with all locations showing >20% exogenous haplotypes. By contrast, the two locations examined in our study from the Adige river system only showed 1 out of 30 (Prissiano) and 1 out of 72 (Adige) exogenous haplotypes. The differences across studies might be due to the different tributaries or river sections sampled, and/or to change of stocking practices along time.
Our study also identified for the first time the presence of introgressed marble trout populations in the river Soca drainage off Slovenia both at the mitochondrial and nuclear level. While our data confirmed the Predelica sample to be pure marble trout in accordance with Snoj et al.  and Fumagalli et al. , introgression was evident at Koritnica (about only 15 km from Predelica) and Isonzo (on the Italian side of the Soca river).
The majority of exogenous mitochondrial haplotypes observed in our marble trout data corresponded to the allochthonous Atlantic haplotype AT-S1. The occurrence of Atlantic haplotypes is unlikely to result from natural migration and might be the consequence of restocking practices. The AT-S1 haplotype group is found at high frequencies in hatchery populations and is thought to represent the predominant haplotype in commercial hatchery strains throughout southern Europe [22, 23]. Massive restocking has been conducted in the Po river and its tributaries, to the extent that hybrids with intermediate morphometric traits between brown and marble trout are common [17, 18, 22]. Similarly, the presence in high frequencies of the nonindigenous Danubian haplotype DA-S1 in the samples of the Ticino region (23-45%) might be consequence of the documented translocations of individuals from Soca river (Slovenia) to Ticino and Adda rivers for restocking purposes . This might also explain (i) the presence of the marble trout MA-S1 haplotype in Ticino, which is typical of the Slovenia region and absent west off Piave, and (ii) the presence in Ticino and Adda of the exclusive marble trout haplotype MA-S4, which is more similar to the Slovenian haplotype (MA-S1) than to the West Italian haplotype (MA-S2). Alternatively, a natural colonisation of Danubian individuals carrying the DA-S1 haplotype due to temporary connections between basins during the Pleistocene cannot be ruled out. It is known that during the Lago Mare period following the Messinian salinity crisis many exchanges and dispersions of brackish and freshwater fish fauna occurred . A similar scenario has been suggested by Meraner et al.  in order to explain the presence of DA-S1 haplotypes in the Adige river, with a possible connection via the Drave river drainage allowing migration between Danubian and Adriatic neighbouring tributaries.
Genetic Differentiation at micro- and macro-geographic scale
Our study reveals strong genetic differentiation among native populations across the geographical range of marble trout using microsatellite data, which suggests restricted gene flow both at the micro-geographic (within rivers) and macro-geographic (among river systems) scale.
Genetic heterogeneity was found not only between geographically distant rivers but also within the same river basin. Samples collected no more than 15 km apart in tributaries of the same river, for instance Adige and Prissiano, resulted significantly differentiated, indicative of limited genetic exchange. However, genetic differentiation in our study was much lower in comparison with the study of Fumagalli et al.  on Slovenian marble trout. The authors found extreme FST values between samples using microsatellite loci, with an average pairwise FST value of 0.79 (about three times higher than the average pairwise FST value found in our study of 0.276), indicative of complete absence of genetic exchanges. They suggested the remnant marble trout populations from Slovenia to be completely isolated from each other and from the main river system as immigration from downstream is precluded by impassable barriers. Furthermore, in the populations analyzed by Fumagalli et al. , microsatellite allele numbers (1.2-2.0) and heterozygosities (He = 0.04- 0.22) were remarkably low, pointing to possible past or present bottleneck episodes. By comparison, most populations in our study showed expected heterozygosities up to 0.388 and allelic richness ranging between 2.3 and 3.4, with the only exception being Predelica (He = 0.08; AR = 1.4), a location in the Soca river in Slovenia shared with the study of Fumagalli et al. . Concordantly, an Approximately Bayesian Calculation approach suggested the smallest effective population size for Predelica (Ne = 858). The fact that the rest of locations in our study showed (i) moderate to high genetic diversity, (ii) no signatures of a genetic bottleneck, (iii) and an inferred effective population size over 2000 individuals in all cases suggests that marble trout population in Northern Italy are not extremely isolated and small, unlike the secluded headwater streams of the Soca river.
At the macro-geographic level microsatellite data reveal a highly significant genetic differentiation consistent with a pattern of Isolation-by-Distance. A significant correlation was found between genetic and hydrographic distances that suggests a reduction in gene flow among populations directly related to the increase of waterway distances among sampling locations. According to the Isolation-by-Distance effect, positioning of samples on the horizontal axis of the Multi-Dimentional Scaling plots (Figure 3) clearly reflected the West-to-East geographic distribution of the samples. This pattern was supported by the distribution of mitochondrial haplotypes, with haplotype MA-s2 being fixed in all tributaries from Eastern and Central Italy and gradually decreasing in frequency moving toward Slovenian rivers, where it is replaced by haplotype MA-s1.
It is generally assumed that the intra-specific genetic structure largely reflects the timing and the pattern of habitat colonisation across the entire range distribution of the species. This holds true especially for freshwater species, for which colonisation patterns have been strongly determined by the geological evolution of the hydrographic networks in the past. The Wurmian glacial period (with a maximum peak 18000 years ago) is thought to have played a fundamental role in determining habitat expansion and contraction of most of native freshwater species in the North-Italian and North-western Balkan regions . During the last glacial age, the river Po extended its natural palaeocourse beyond the "meso-Adriatic" sea pit, gathering in the same drainage all the main rivers of the north-east Italian (rivers Tagliamento, Piave, Sile, Brenta and Adige and their hydrographic networks) and Balkan (river Soca) districts. Such a scenario, with a single non-fragmented freshwater basin, suggest an initial common gene pool for all the North-Italian/Slovenian marble trout populations. At the beginning of the postglacial period, the receding of the coastline determined a progressive submerging of the relict outlet of the Po river. This process sequentially isolated the outlet of the rivers presently flowing into the North Adriatic Sea, likely determining a strong reduction of the gene flow among populations confined in different river basins. As a consequence of confinement, individuals tended to find mates from the same or nearby populations rather than distant populations, and as a result, populations that live near each other became genetically more similar than populations that live further apart, generating the pattern of Isolation-by-Distance found in our study.
Besides historical patterns of isolation among drainages, the high rates of introgression may have affected the native patterns of population genetic structure. While a similar East-to-West partition of the microsatellite polymorphism following an Isolation-by-Distance pattern and comparable FST values were found when including all individuals in the analysis and when including only native (non-introgressed) individuals, some introgression could remain considering that we excluded introgressed individuals using a 10% criterion. In addition to past biogeographic events, the higher introgression observed in Central Italy could contribute to the deep genetic break between Pellice and Ticino, and between Isonzo and Predelica, and the differences within river basins. In regards to the differences in genetic diversity, introgression might have also contributed to a higher native genetic diversity in Central Italy in comparison with the more extreme locations (Predelica, Pellice), but the fact that a similar native genetic diversity is observed in locations with high (Ticino, Brenta) and low introgression (Adige) suggests that introgression alone cannot account for those differences.
Preserving genetic diversity within a population or species sustains its evolutionary potential to cope with current and future environmental changes. In this sense, native species are specially vulnerable to human activities including introduction of nonindigenous species. Loss of genetic integrity as a consequence of introgression of exogenous gene pools is considered the most serious thread to native gene pools. Thus, detecting hybridization and introgression in wild populations is a primary concern in developing conservation and management strategies. The risk of genetic erosion caused by introgression cannot be ignored as it obscures native gene pools through alteration of genetic composition, reduction of fitness, diminution of population sizes, loss of genetic diversity and extinction of local populations.
Materials and methods
A total of 429 marble trout Salmo trutta marmoratus individuals were caught using electrofishing over the period 1995-1997. Collecting sites included most of the geographic distribution of the species, ranging from North-western Italy to Slovenia (Figure 1). Samples were obtained at the headwaters of eleven tributaries of the Po and Soca (known as Isonzo in Italy) rivers. During field work, fin clips were collected from anaesthetised individuals that were immediately released back to the water. Specimens were photographed and classified as marble trout on the basis of typical morphological features and skin colour patters according to the criteria described by Gandolfi et al. . An additional pure brown trout sample (confirmed at the mitochondrial and microsatellite level) was obtained at the river Chisone in North-West Italy for introgression studies.
DNA was extracted from ethanol-preserved finclips following a resin-based protocol (Chelex 100; Bio-rad) as described in Walsh et al. . A 352 bp segment of the mtDNA control region (D-loop) was amplified using primers LPro2 (5' AACTCCCACCACTAACTCCCAAAGC 3'; ) and H2 (5' CGTTGGTCGGTTCTTAC 3'; ) in a GeneAmp System 9700 Termocycler (Applied Biosystems). PCR products were used for SSCP (Single Strand Conformation Polymorphism) analysis . Electrophoresis was performed for 40 h in a 12% acrylamide gel at constant room temperature (4°C). DNA fragments were visualized by silver staining. All rare (low frequency) electrophoretic profiles plus a subset of 25 individuals from each common profile were sequenced using a 725 VISTRA automatic sequencer (Molecular Dynamics). Sequences were analyzed using the CHROMAS software (http://www.technelysium.com.au/chromas.html) and aligned using CLUSTALW .
All samples were scored for five microsatellite loci: four previously described in brown trout Salmo trutta, Str15 , Str85, Str543, Str591 ; and one previously described in Atlantic salmon Salmo salar: Ssa85 . PCR products were obtained in a GeneAmp PCR System 9700 Thermocycler (Applied Biosystems). PCR reactions consisted of 1 μl template DNA, 0.8 μl 10X reaction buffer (Promega), 0.4 μl 1 mM dNTPs, 0.8 μl 25 mM MgCl2, 0.3 μl 10 mM forward and reverse primers, 0.04 μl Taq polymerase (Promega, 5 u/μl) and water up to 8 μl. PCR products were run on a 725 VISTRA automatic sequencer (Molecular Dynamics) and analyzed using FRAGMENTOR (Molecular Dynamics).
Analysis of mtDNA data was performed using ARLEQUIN v. 3.5 , including computation of haplotype frequencies and exact test of population differentiation. Observed haplotypes were compared with those considered typical of marble trout and brown trout  in order to detect possible signatures of introgression. The presence of exogenous DNA material was also tested at microsatellite loci using the software STRUCTURE , a model based clustering algorithm that infers the most likely number of groups (species) in the data. The analysis included all marble trout samples from our study plus a sample of brown trout obtained in the river Chisone (North-western Italy, N = 23). The software organizes individuals into a predefined number of clusters (K) with a given likelihood, which may represent putative populations or species. The analysis was performed with 1 < K < 10 to account for population substructuring within species, using the admixture model. The most likely K was determined using the criterion of Evanno et al.  and then used to assign each individual. A burn in length of 104 iterations followed by 105 additional Markov Chain Monte Carlo (MCMC) iterations were performed as recommended by Pritchard et al. .
At microsatellite loci, within-sample genetic diversity statistics were assessed by observed (Ho) and expected (He) heterozygosities per locus using GENETIX version 4.05 , and allelic richness (AR) using FSTAT . Differences in genetic diversity among samples were tested by one-way ANOVA using STATISTICA version 6.0 (StatSoft Inc.). Deviations from Hardy-Weinberg Equilibrium (HWE), linkage disequilibrium, and differences in allele and genotype frequencies among samples were tested using GENEPOP version 3.4 . Presence of null alleles was tested using the program MICRO-CHECKER version 2.2.3 .
Prior to the population structure analysis, the statistical power of the markers employed in our study for investigating genetic differentiation was assessed with the program POWSIM . We tested a range of predefined levels of expected divergence (FST = 0.001, 0.005, 0.01, 0.05 and 0.1).
Population structure was studied using non-hierarchical and hierarchical F-Statistics  calculated using ARLEQUIN . Significance tests were assessed with 10,000 permutation tests. In all cases, significance levels were corrected for multiple comparisons using Bonferroni . Pairwise multilocus comparisons between samples were calculated by Cavalli-Sforza and Edwards chord distance  and a multivariate ordination was conducted by Multi-Dimensional Scaling (MDS) analysis using STATISTICA version 6.0 (StatSoft). Isolation-by-Distance (IBD) was tested using single and partial Mantel tests implemented in GENETIX, by correlating linearized genetic distance (FST/(1- FST)) vs. geographic distance (shortest waterway distance between pairs of samples).
The software DIYABC  was used to estimate the effective population size for all population samples. The software allows to make inference based on Approximate Bayesian Computation (ABC), in which scenarios can be customized by the user to fit many complex situations involving any number of populations and samples. We estimated effective population size for each population. For all DIYABC runs, prior uniform distributions were set for Ne (lower limit 100 and upper limit 100 000), mean μ (lower limit 1 × 10-5 and upper limit 1 × 10-3) and locus μ (lower limit 1 × 10-5 and upper limit 1 × 10-2). Microsatellite repeat motif and allele range were specified for each locus. The stepwise mutation model was used and 500,000 replicate runs were performed to generate the reference tables.
Finally, we tested for a genetic bottleneck episode in two different ways. First, we used the software BOTTLENECK , based on the principle that after a recent reduction of their effective population size, number of alleles (k) decreases faster than heterozygosity (He) at polymorphic loci. Thus, in a recently bottlenecked population, the observed gene diversity is higher than the expected equilibrium gene diversity (Heq) which is computed from the observed number of alleles (k), under the assumption of a constant-size (equilibrium) population . We used the Multiple-step Stepwise (TPM) model , which consists of mostly one-step mutations but a small percentage of multi-step changes, and is the recommended model for microsatellite data sets rather than the Infinite Alleles (IAM) or Single-step Stepwise (SMM) models . The proportion of singlestep mutation events was set to 90%, variance = 12%. Observed and expected heterozygosities were compared using a Wilcoxon sign-rank test as suggested by Piry et al. . Second, we calculated M, the mean ratio between number of alleles (k) and range in allele size (r), assuming that during a bottleneck episode k decreases faster than r (M_P_Val; ). Hence, the value of M decreases when a population is reduced in size. Average M was calculated across loci and compared with the critical value Mcrit estimated after 10,000 simulations and assuming the population to be at equilibrium. In all simulations, three different values of θ were used (5, 10 and 20). A range of mutation models were examined and conservative values were used for ps (frequency of one-step mutations) and Δg (average size of non one-step mutations), ps = 0.90 and Δg = 3.5.
We are grateful to the following colleagues who provided support in sampling: G. Forneris, F. Merati, A. Atz, P. Turin, R. Pria, G. Amirante, G. Sigalotti and M. Bertok. We thank L Bargeloni for his very useful technical and conceptual advise. We also thank Prof G. Gandolfi, L. Zane and L. Congiu for useful discussions. This work was supported by a doctoral fellowship to ANL and by a grant from MURST to TP.
- Ward RD, Woodward M, Skibinski DOF: A comparison of genetic diversity levels in marine, freshwater and anadromous fishes. J Fish Biol. 2004, 44: 213-227. 10.1111/j.1095-8649.1994.tb01200.x.View ArticleGoogle Scholar
- Volckaert FAM, Hanfling B, Hellemans B, Carvalho GR: Timing of the population dynamics of bullhead (Cottus gobio) during the Pleistocene. J Evol Biol. 2002, 15: 930-944. 10.1046/j.1420-9101.2002.00469.x.View ArticleGoogle Scholar
- Ryman N: Patterns of distribution of biochemical genetic variation in salmonids: differences between species. Aquaculture. 1983, 33: 247-262. 10.1016/0044-8486(83)90382-4.View ArticleGoogle Scholar
- Bernatchez L: The evolutionary history of brown trout (Salmo trutta) inferred from phylogeographic, nested clade, and mismatch analyses of mitochondrial DNA variation. Evolution. 2001, 55: 351-379.View ArticlePubMedGoogle Scholar
- Bernatchez L, Guyomard R, Bonhomme F: DNA sequence variation of the mitochondrial control region among geographicall yand morphologically remote European brown trout Salmo trutta populations. Mol Ecol. 1992, 1: 161-173. 10.1111/j.1365-294X.1992.tb00172.x.View ArticlePubMedGoogle Scholar
- Vera M, Cortey M, Sanz N, Garcia-Marin JL: Maintenance of an endemic lineage of brown trout (Salmo trutta) within the Duero river basin. J Zool Syst Evol Res. 2009, 48: 181-187. 10.1111/j.1439-0469.2009.00547.x.View ArticleGoogle Scholar
- Garcia-Marin JL, Pla C: Origins and relationships of native populations of brown trout (Salmo trutta) in Spain. Heredity. 1996, 77: 313-323. 10.1038/hdy.1996.145.View ArticleGoogle Scholar
- Machordom A, Suarez J, Almodovar A, Bautista JM: Mitochondrial haplotype variation and phylogeography of Iberian brown trout populations. Mol Ecol. 2000, 9: 1324-1338. 10.1046/j.1365-294x.2000.01015.x.View ArticlePubMedGoogle Scholar
- Sanz N, Garcia-Marin JL, Pla C: Divergence of brown trout (Salmo trutta) within glacial refugia. Can J Fish Aquat Sci. 2000, 57: 2201-2210. 10.1139/cjfas-57-11-2201.View ArticleGoogle Scholar
- Sanz N, Garcia-Marin JL, Pla C: Managing fish populations under mosaic relationships. The case of brown trout (Salmo trutta) in peripheral Mediterranean populations. Conserv Genet. 2002, 3: 385-400. 10.1023/A:1020527420654.View ArticleGoogle Scholar
- Cortey M, Garcia-Marin JL, Pla C: Historical biogeography of Mediterranean brown trout. Mol Phyl Evol. 2004, 57: 2201-2210.Google Scholar
- Crivelli AJ, Poizat G, Berrebi P, Jesensek D, Rubin JF: Conservation biology applied to fish: the example of a project for rehabilitating the marble trout in Slovenia. Cybium. 2000, 24: 211-230.Google Scholar
- Sommani E: Il Salmo marmoratus: sua origine e distribuzione nell'Italia settentrionale. Bollettino di Pesca, Piscicoltura e di Idrobiologia. 1960, 17: 243-248.Google Scholar
- Berrebi P, Povz M, Jesensek D, Cattaneo-Berrebi G, Crivelli : The genetic diversity of native, stocked and hybrid populations of marble trout in the Soca river, Slovenia. Heredity. 2000, 85: 277-287. 10.1046/j.1365-2540.2000.00753.x.View ArticlePubMedGoogle Scholar
- Snoj A, Jug T, Milkic E, Susnik S, Pohar J, Dovc P: Mitochondrial and microsatellite DNA analysis of marble trout in Slovenia. Quaderni ETP. 2000, 29: 5-11.Google Scholar
- Fumagalli L, Snoj A, Jesensek D, Balloux F, Jug T, Duron O, Brossier F, Crivelli AJ, Berrebi P: Extreme genetic differentiation among the remnant populations of marble trout Salmo marmoratus in Slovenia. Mol Ecol. 2002, 11: 2711-2716. 10.1046/j.1365-294X.2002.01648.x.View ArticlePubMedGoogle Scholar
- Giuffra E, Bernatchez L, Guyomard R: Mitochondrial control region and protein coding genes sequence variation among phenotypic forms of brown trout Salmo trutta from northern Italy. Mol Ecol. 1994, 3: 161-171. 10.1111/j.1365-294X.1994.tb00117.x.View ArticlePubMedGoogle Scholar
- Giuffra E, Guyomard R, Forneris G: Phylogenetic relationships and introgression patterns between incipient parapatric species of Italian brown trout Salmo trutta complex. Mol Ecol. 1996, 5: 207-220.View ArticleGoogle Scholar
- Crivelli AJ: Are fish introductions a threat to endemic freshwater fishes in the northern Mediterranean region. Biol Conserv. 1995, 72: 311-319. 10.1016/0006-3207(94)00092-5.View ArticleGoogle Scholar
- Meraner A, Baric S, Pelster B, Dalla Via J: Trout (Salmo trutta) mitochondrial DNA polymorphism in the centre of the marble trout distribution area. Hydrobiologia. 2007, 579: 337-349. 10.1007/s10750-006-0479-3.View ArticleGoogle Scholar
- Meraner A, Baric S, Pelster B, Dalla Via J: Mitochondrial DNA data point to extensive but incomplete admixture in a marble and brown trout hybridization zone. Conserv Genet. 2010, 11: 985-998. 10.1007/s10592-009-9942-9.View ArticleGoogle Scholar
- Argenton F, Bargelloni L, Patarnello T, Colombo L, Bortolussi M: Risk of introgressive hybridization between fario and marmoratus morphs of Salmo trutta in north-eastern Italy as evidenced by mitochondrial DNA analysis. Riv Italiana Acquacoltura. 1992, 27: 119-126.Google Scholar
- Weiss S, Schlotterer C, Waidbacher H, Jungwirth M: Haplotype (mtDNA) diversity of brown trout Salmo trutta in tributaries of the Austrian Danube: massive introgression of Atlantic basin fish - by man or nature?. Mol Ecol. 2000, 10: 1241-1246. 10.1046/j.1365-294X.2001.01261.x.View ArticleGoogle Scholar
- Economidis P, Miller P: Systematics and freshwater gobies from Greece. J Zool. 1990, 221: 125-170. 10.1111/j.1469-7998.1990.tb03781.x.View ArticleGoogle Scholar
- Climap : The surface of the ice-age Earth. Science. 1976, 191: 1131-1144. 10.1126/science.191.4232.1131.View ArticleGoogle Scholar
- Gandolfi G, Zerunian S, Torricelli P, Marconato A: I pesci delle acque interne italiane. 1991, Ministerio dell'Ambiente e Unione Zoologica Italiana, RomeGoogle Scholar
- Walsh PS, Metzger DA, Higuchi R: Chelex-100 as a medium for simple extraction of DNA for PCR-based typing from forensic material. Biotechniques. 1991, 10: 298-317.Google Scholar
- Astolfi L, Dupanloup I, Rossi R, Bisol PM, Faure E, Congiu L: Mitochondrial variability of sand smelt Atherina boyeri populations from North Mediterranean coastal lagoons. Mar Ecol Prog Ser. 2005, 297: 233-243. 10.3354/meps297233.View ArticleGoogle Scholar
- Bernatchez L, Danzmann G: Congruence in control region sequence and restriction site variation in mitochondrial DNA of brook charr (Salvelinus fontinalis). Mol Biol Evol. 1993, 10: 1002-1014.Google Scholar
- Orita M, Iwahana H, Kanazawa H, Hayashi K, Sekiya T: Detection of polymorphisms of human DNA by gel electrophoresis as single-strand conformation polymorphisms. Proc Natl Acad Sci USA. 1989, 86: 2766-2770. 10.1073/pnas.86.8.2766.PubMed CentralView ArticlePubMedGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucl Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralView ArticlePubMedGoogle Scholar
- Estoup A, Presa P, Krieg F, Valman F, Guyomard R: (CT)n and (GT)n microsatellites: a new class of genetic markers for Salmo trutta (brown trout). Heredity. 1993, 71: 488-496. 10.1038/hdy.1993.167.View ArticlePubMedGoogle Scholar
- Presa P, Guyomard R: Conservation of microsatellites in three species of salmonids. J Fish Biol. 1996, 49: 1326-1329.Google Scholar
- O'Reilly PT, Hamilton LC, McConnell SK, Wright JM: Rapid analysis of genetic variation in Atlantic salmon (Salmo salar) by PCR multiplexing of dinucleotide and tetranucleotide microsatellites. Can J Fish Aquat Sci. 1996, 53: 2292-2298.Google Scholar
- Schneider S, Roessli D, Excoffier L: ARLEQUIN: a software for population genetics data analysis ver 3.0. 2000, Genetics and Biometry Lab, Dept. of Anthropology, University of Geneva, SwitzerlandGoogle Scholar
- Pritchard JK, Stephens M, Donelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMed CentralPubMedGoogle Scholar
- Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.View ArticlePubMedGoogle Scholar
- Belkhir K, Borsa P, Goudet J, Bonhomme F: GENETIX 4.05: Logiciel sous Windows pour la genetique des populations. 2005, Laboratoire Genome & Population, CNRS-UPR/Universite de Montpellier II, Montpellier, FranceGoogle Scholar
- Goudet J: FSTAT version 184.108.40.206. 2002, [http://www2.unil.ch/popgen/softwares/fstat.htm]Google Scholar
- Raymond M, Rousset F: GENEPOP (version 1.2): a population genetics software for exact tests and ecumenicism. J Hered. 1995, 86: 248-249.Google Scholar
- Van Oosterhout CV, Hutchinson WF, Wills DPM, Shipley P: Micro-checker: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol. 2004, 4: 535-538. 10.1111/j.1471-8286.2004.00684.x.View ArticleGoogle Scholar
- Ryman N, Palm S: POWSIM: a computer program for assessing statistical power when testing genetic differentiation. Mol Ecol Notes. 2006, 6: 600-602. 10.1111/j.1471-8286.2006.01378.x.View ArticleGoogle Scholar
- Weir BS, Cockerman CC: Estimating F-Statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.View ArticleGoogle Scholar
- Rice WR: Analyzing tables and statistical tests. Evolution. 1989, 43: 223-225. 10.2307/2409177.View ArticleGoogle Scholar
- Cavalli-Sforza LL, Edwards AWF: Phylogenetic analysis models and estimation procedures. Evolution. 1967, 32: 550-570. 10.2307/2406616.View ArticleGoogle Scholar
- Cornuet JM, Santos F, Beaumont MA, Robert CP, Marin JM, Balding DJ, Guillemaud T, Estoup A: Inferring population history with DIYABC: a user-friendly approach to Approximate Bayesian Computation. Bioinformatics. 2008, 24: 2713-2719. 10.1093/bioinformatics/btn514.PubMed CentralView ArticlePubMedGoogle Scholar
- Piry SG, Luikart G, Cornuet JM: BOTTLENECK: a computer program for detecting recent reductions in the effective population size using allele frequency data. J Hered. 1999, 90: 502-503. 10.1093/jhered/90.4.502.View ArticleGoogle Scholar
- Cornuet JM, Luikart G: Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics. 1996, 144: 2001-2014.PubMed CentralPubMedGoogle Scholar
- Di Rienzo A, Peterson AC, Garza JC, Valdes AM, Slatkin M, Freimer NB: Mutational processes of single-sequence repeat loci in human populations. Proc Natl Acad Sci USA. 1994, 91: 3166-3170. 10.1073/pnas.91.8.3166.PubMed CentralView ArticlePubMedGoogle Scholar
- Luikart G, Sherwin WB, Steele BM, Allendorf FW: Usefulness of molecular markers for detecting population bottlenecks via monitoring genetic change. Mol Ecol. 1998, 7: 963-974. 10.1046/j.1365-294x.1998.00414.x.View ArticlePubMedGoogle Scholar
- Garza JC, Williamson EG: Detection of reduction in population size using data from microsatellite loci. Mol Ecol. 2001, 10: 305-318. 10.1046/j.1365-294x.2001.01190.x.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.