Issue |
Aquat. Living Resour.
Volume 32, 2019
|
|
---|---|---|
Article Number | 23 | |
Number of page(s) | 10 | |
DOI | https://doi.org/10.1051/alr/2019021 | |
Published online | 10 October 2019 |
Research Article
Genetic analysis of populations of brown trout (Salmo trutta L.) from the Romanian Carpathians
1
Department of Biochemistry and Molecular Biology, Faculty of Biology, University of Bucharest, 91-95 Splaiul Independenței, Bucharest, Romania
2
Applied Ecology Research Center, Faculty of Sciences, Lucian Blaga University of Sibiu, 5-7 Ion Rațiu Street, Sibiu, Romania
3
Department of Environmental Sciences, Faculty of Sciences, Lucian Blaga University of Sibiu, 5-7 Ion Rațiu Street, Sibiu, Romania
4
Department of Biology, Ecology and Environmental Protection, Faculty of Sciences, Vasile Alecsandri University of Bacău, 158 Calea Mărășești, Bacău, Romania
* Corresponding author: georgescu_se@yahoo.com
Handling Associate Editor: Carlos Saavedra
Received:
3
June
2019
Accepted:
4
September
2019
The Carpathian Mountains are one of the most complex orogenetic areas of Europe, with unique fauna, including the brown trout (Salmo trutta). In this study we performed population genetic analysis of 12 different S. trutta populations using two types of molecular markers: nine microsatellites and mitochondrial D-loop sequences. The following working hypothesis was considered: the Romanian Carpathians and their surrounding lowlands can be key relief units based on which the S. trutta genetic diversity, spread, distribution, connectivity, relative isolation and genetic divergence can be at least partially explained. The phylogenetic analysis revealed that the majority of sequences were grouped in the Danubian clade. The high haplotype diversity of the 12 analyzed brown trout populations can be explained by the high nucleotide diversity. The microsatellite analysis revealed an inbreeding event for all the loci and for the populations analyzed. The Romanian Carpathians' shape and geographic orientation play a zoogeographical key role driving force in respect to the S. trutta populations.
Key words: Brown trout / Carpathians Mountains / genetic diversity / microsatellites / mitochondrial DNA
© EDP Sciences 2019
1 Introduction
The brown trout (Salmo trutta, Linnaeus, 1758) is a member of the Salmonidae family; the only family currently placed in the order Salmoniformes (Behnke, 2002; Klemetsen et al., 2003). Members of the Salmonidae are fish species that rely on a multitude of habitats (resting, sheltering, feeding, spawning, nursery, etc.) since in this family are included both migratory and resident fish species. Moreover, these habitats, often differ from one age class to another (Elliott, 1994; Crisp, 2000). For the conservation of this fish species, the connection between the spawning and maturation habitats, alongside the habitats themselves must be preserved (Schwartz et al., 2007). In consequence, there is a need for regional management strategies that rely on supporting genetic diversity, and therefore on the adaptability potential of the populations (Reed and Frankham, 2003). Increased genetic diversity is also associated with greater resilience in the face of exploitation (Hilborn et al., 2003; Schindler et al., 2010).
S. trutta is a common fish species in the Romanian Carpathians and the most important salmonid species for angling, being found in many streams including in our area of interest. During the reproduction period the brown trout migrates upstream in the rivers and their tributaries. In general, the genetic diversity studies of brown trout populations show a decreasing trend of within population genetic diversity from Western to Eastern Europe (Kohout et al., 2013), which is clear in Romania when the few studies available are being considered.
All developmental stages of salmonids can be severely affected by environmental changes like low or high flow conditions and different alterations due to anthropic interventions (Warren et al., 2015). Therefore, genetic diversity is needed by the populations for adequate adaptive/evolutionary potential (Frankel and Soule, 1981). One of the ways to preserve genetic diversity is admixture, which takes place when spawners migrate haphazardly from one small tributary to another (Ostergren and Nilsson, 2012). This scenario occurs more often when the tributaries are spatially close, thus the genetic differences of the local populations are correlated with the geographical distance within the river system (Carlsson and Nilsson, 2000). Furthermore, anthropogenic activities and disastrous natural events can cause a reduction in population size, which increases inbreeding and genetic drift (Cunjak and Power, 1986). A series of studies have reported dams or hydropower plants constructions on some Romanian mountain streams and rivers that could have negative impacts on the aquatic and riverine habitats and fish communities (Bănăduc, 1999; Burghelea et al., 2013).
Non-indigenous Atlantic brown trout can hybridise with local populations, a phenomenon which took place between the Atlantic and Danubian lineages in the Upper/Middle Danube River basin (Hansen, 2002; Duftner et al., 2003; Sanz et al., 2006; Simonović et al., 2017). There are studies that highlight translocation of non-indigenous brown trout individuals from former Czechoslovakia to Bulgaria, and sporadic stocking in Serbia (Marić et al., 2006; Kohout et al., 2013). Attention must be drawn to the fact that the genetic variability between the brown trout populations is lost in central European parts, of the North, Black and Baltic Sea basins due to various restocking activities with individuals of unknown genetic diversity (Wenne, 2001; Wlodarczyk and Wenne, 2001; Was and Wenne, 2003; Kohout et al., 2012; Schenekar et al., 2014; Wenne et al., 2016; Osz et al., 2018).
The Carpathians represent an area where frequent chaotic stocking and restocking of S. trutta. These activities were carried out without considering the status of the population and the origin (lineage) of the individuals used for restocking. In this context, the genetic studies are more than necessary in respect to understanding this fish species' past zoogeographical paths and mechanisms, the present ecological status trends and the support for an appropriate management strategy (Curtean-Bănăduc et al., 2015). The necessity for research in this field is necessary, especially in the context of no public reports stating previous stocking activities (Didenko et al., 2011, 2014). Only few genetics studies were conducted on brown trout populations from the Romanian Carpathians (Popa et al., 2013, 2016; Nechifor et al., 2017) and a more complex analysis is necessary.
Given this context, our study aims to analyse the genetic diversity of some Romanian Carpathians brown trout populations by using nuclear and mitochondrial molecular markers. The study presented here was conducted in aleatory chosen rivers from the Romanian Carpathian Mountains, on their East to West axis. This paper's working hypothesis is that the Carpathians and their surrounding lowlands can be key relief units based on which the S. trutta genetic diversity, spread, distribution, connectivity, relative isolation and genetic divergence can be at least partially explained.
2 Materials and methods
2.1 Sampling and laboratory analyses
A total of 362 samples of brown trout was collected between 2012 and 2017. The individuals originated from 12 river drainages of the Danubian basin (Fig. 1). Fin clips were preserved in 96% ethanol and stored at 4 °C. The genomic DNA was extracted using the standard method with phenol/chloroform (Taggart et al., 1991).
A fragment of the mtDNA control region (partial D-loop region) of 1003 bp was amplified using the PST-FST primer pair (Cortey and Garcia-Marin, 2002) with the following PCR conditions: 95 °C for 10 min, 40 cycles at 95 °C for 30 s, 52 °C for 30 s and 72 °C for 1 min, followed by final extension at 72 °C for 10 min. Amplified fragments were sequenced on ABI Prism 3130 Genetic Analyzer (Applied Biosystems). Sequences were revised using BIOEDIT (Hall, 1999) and aligned using MAFFT online v7 (https://mafft.cbrc.jp/alignment/software/). For microsatellite analyses, we used nine primer pairs (Tab. S1) grouped in two 3-plexes (I and II), one duplex (III) and one monoplex (IV). The PCR conditions for the microsatellite amplification were: 95 °C for 10 min, 35 cycles at 95 °C for 30 s, group specific hybridisation temperature (Tab. S1) for 30 s and 72 °C for 1 min, followed by final extension at 72 °C for 60 min. The amplification reactions were done using a Veriti Thermal Cycler (Applied Biosystems) in a final volume of 25 µl with 1X PCR Buffer, 1.5 mM of MgCl2, 0.8 mM of dNTPs, 20 pmol of each primer, 1 unit of AmpliTaq Gold DNA polymerase, nuclease free water and 50 ng of DNA template. Amplified fragments were separated on an ABI Prism 310 Genetic Analyzer and the allele size was determined relative to the LIZ-500 Size Standard (Applied Biosystems) size standard using GENEMAPPER 4.0 (Applied Biosystems).
Fig. 1 Map of locations. (A) Sampling sites: Uz − Uz River, Otz − Oituz River, Uc − Ucea River, Por–Porumbacu River, Avg − Avrig River, Crt − Cârțișoara River, Top − Topolog River, BM − Bistra Mărului River, Sb − Sebeşel River, Cr − Craiului River, Bz − Bârzava River, Car − Caraş River. Rivers Porumbacu, Avrig, Cârțișoara, Ucea, Sebeşel discharge to the Northern/inner slope, while Toplog River discharge to the Southern/outer slope of the mountain ridge. (B) Map of region showing the locations of 12 sampling sites. |
2.2 Data analyses
For mtDNA, the number of haplotypes and the haplotype diversity indices were computed using DNASPv5 (Librado and Rozas, 2009). To identify new haplotypes and to reveal their phylogenetic relationships, all sequences of the brown trout control region of appropriate length available from GenBank were included in the analysis. These sequences were belonging to the following evolutionary lineages: Atlantic − At, Danubian-Da, Mediterranean − Me, and Adriatic − Ad (Bernatchez, 2001). Sequences from a separate lineage designated as marmoratus − Ma corresponding to marble trout (S. marmoratus) were also included in the data set despite the controversy about its taxonomic status (Pustovrh et al., 2014). The marble trout is considered either a separate species (Berrebi et al., 2000; Fumagalli et al., 2002; Splendiani et al., 2006) or a member of the S. trutta complex (Bernatchez, 2001; Meraner et al., 2007).
The AF133701–Salmo salar mitochondrial region 15662–16669 bp was used as outgroup (Tab. S2). The evolutionary history was inferred using the Maximum Likelihood method based on the Tamura 3-parameter model with a set of 141 D-loop sequences (the unique haplotypes identified in the 12 analysed populations along with sequences downloaded from GenBank, including the outgroup S. salar) with MEGA6 (Tamura et al., 2013) and the statistical support for branching patterns was estimated by 1000 bootstrap replications. A discrete Gamma distribution was used to model evolutionary rate differences among sites (5 categories, +G parameter =0.2809). The tree was drawn to scale, with branch lengths measured in the number of substitutions per site. There were a total of 946 positions in the final dataset. For an accurate phylogenetic analysis of the data set the extremely variable homopolimeric T region of D-loop was cut out as it might be lead to a phylogenetic pattern that not reflect the true phylogenetic relationships. The phylogenetic tree was visualized with FIGTREE v1 (http://tree.bio.ed.ac.uk/software/figtree/).
For microsatellites, we tested all loci with Micro-Checker (Van Oosterhout et al., 2004). The Polymorphic Information Content (PIC) values for each locus, alongside tests for deviation from Hardy-Weinberg proportions were computed with POWERMARKERv3.25 (Liu and Muse, 2005). Allele frequencies, FST values between pairs of populations and values of FIS, along with the Nm parameter values were computed with ARLEQUIN (Excoffier and Lischer, 2010) and tested for significance using 1023 nonparametric permutations. GENALEXv6.5 (Peakall and Smouse, 2006) was used to find private alleles in each population and in the 12 rivers. The same software was used to determine the PCoA based on Nei's genetic distance (D) between the individuals, along with the observed (HO) and expected (HE) heterozygosity. To explore indications of recent bottleneck events, the Garza–Williamson index was calculated across loci with the Arlequin software. This index is a ratio calculated by using the number of alleles and the range in allele size, based on the assumption that the number of alleles declines faster than the range in allele size during a bottleneck, while a value M lower than 0.68 can mean that a population has gone through a recent reduction in size (Garza and Williamson, 2001). The Bayesian-based clustering method in STRUCTUREv2.3.4 (Pritchard et al., 2000) was applied to infer the population structure and to reveal potential admixture between populations, without a priori assigned individuals to populations. The most probable number of genetic clusters (K) was estimated based on the posterior probability of the data for a given K and clarified using STRUCTURE HARVESTER (Earl and Von Holdt, 2012). For the estimation, genotypes were assigned into one to 20 groups and 100 iterations with 50.000 burn-in and 100.000 repetitions were applied for each K.
3 Results
3.1 Mitochondrial DNA
Out of the 362 samples processed for sequencing, we have obtained 316 usable D-loop sequences (Tab. 1, Tab. S3), this probably due to the various polymorphisms in this mitochondrial control region (Consuegra et al., 2015).
Among the sequences obtained for the 12 sampling sites, 80 haplotypes were revealed. The number of haplotypes per site varied between 1 for Caraș River and 14 for Bistra Mărului River (Tab. 1). The maximum number of sequences (106) belonged to one haplotype group which was composed of individuals that came from Avrig, Oituz, Uz, Topolog, Cârţişoara, and Porumbacu rivers populations. The haplotype diversity within the populations varied from 0.22 ± 0.10 (Bârzava River) to 0.94 ± 0.03 (Bistra Mărului River), with the most numerous polymorphic sites (46) being found in the Sebeşel population. The most numerous polymorphic sites were found in the first 500 bp of the D-loop sequence. The phylogenetic analysis revealed that the majority of sequences were grouped together in the Danubian clade, with the exception of Por14 (Porumbacu River), included in the Atlantic clade. Some sequences from Uz River (Uz1.1, Uz1.11, Uz2.6, Uz3.7 and Uz4.11) cluster separately. Also, the sequences Uz2.10 and Uz4.10 were grouped separately in a distinct clade together with a sequence from GenBank classified as Danubian haplotype (Fig. 2).
The number of D-loop sequences along with the number of haplotypes and genotypes obtained for each population. The abbreviation for each population name is in the brackets.
Fig. 2 Phylogenetic tree of D-loop sequences obtained with the Maximum Likelihood method based on the Tamura 3-parameter model with 1000 bootstrap replicates. The bootstrap values higher than 50 are displayed. Green circles − Danubian clade, yellow triangles − Atlantic clade, pink squares − Mediterranean clade, blue diamonds − Adriatic clade, red triangles − S. marmoratus. |
3.2 Microsatellites
For the nine microsatellites, a minimum number of alleles were identified for the inner Carpathian arch northern river basins (Cârţişoara, Porumbacu, Ucea, and Avrig) populations (Tab. S4). According to the results of Micro-checker test, we did not detect evidence for genotype inferring errors due to stuttering, neither for large allele dropout, nor for a high frequency of null alleles. A maximum number of alleles were identified for OmyFGT1 locus (24 alleles) in the Uz population (Tab. S4). Most of the loci were not in Hardy-Weinberg equilibrium (p < 0.05), and the PIC parameter varied from 0.03 for Str15, Str60 and Str73 (for Bistra Mărului, Ucea or Bârzava populations) to 0.90. The Garza–Williamson index had values from 0.20 ± 0.13 (Oituz population) to 0.36 ± 0.21 (Porumbacu population), while the observed and expected heterozygosity varied from 0.03 to 0.90 and 0.20 to 0.93, respectively (Tab. S4). Private alleles were recorded in all populations, for all loci, and the mean number of private alleles varied from 0.11 to 2.78 (Tab. S4). The mean coefficient of inbreeding FIS within populations was significant in Cârţişoara, Bistra Mărului, Bârzava, Oituz, Uz, Craiului, Caraş, and Topolog populations, reflecting a homozygote excess (Tab. 2).
The FST values across all Romanian brown trout populations varied from 0.08 (population pairs Porumbacu–Avrig, Topolog–Craiului, Topolog–Caraş) to 0.44 (population pair Ucea–Sebeşel), all values being statistically significant (p < 0.05). A high degree of genetic differentiation (FST > 0.25; p < 0.05) was found for the following population pairs: Ucea-all populations; Cârţişoara–Bistra Mărului; Caraş-Sebeşel; Porumbacu-Caraş and Sebeşel; Avrig–Bistra Mărului and Sebeşel; Bistra Mărului-Bârzava and Oituz; Bârzava-Sebeşel; Oituz-Sebeşel, and Uz–Sebeşel (Tab. 3).
The Bayesian analysis in STRUCTURE revealed that the 12 populations can be genetically grouped into four clusters (K = 4), due to maximum L(K), based on the nine analyzed microsatellites (Fig. 3, Fig. S1). The first cluster is made from all Ucea populations members; the second one represents the Cârţişoara, Porumbacu, Avrig and Bârzava individuals, with one individual genetically closer to the first cluster. The third group is composed of the individuals coming from the Uz and Oituz rivers, while the fourth one grouped the Bistra Mărului, Craiului, Caraş, Sebeşel and Topolog individuals.
Another type of analysis used to investigate the genetic structure of the populations, based on the genetic distance between the individuals, is the PCoA (Fig. 4). The individuals that are genetically close tend to cluster together, and the results from the Structure analysis support this representation: Ucea individuals are separated from the others, followed by a group of populations made of individuals from Cârţişoara and a part of the individuals from Porumbacu, Avrig, Bârzava, Oituz and Uz, and, finally a group formed by individuals from Sebeşel, Topolog, Caraş, Craiului and Bistra Mărului.
The inbreeding coefficient (FIS) values for the nine loci. *p < 0.05, tested with 1000 bootstrap replications. Abbreviated are used for the sampled populations.
FST values (above the diagonal) and Nm values (under the diagonal) per population. *p < 0.05, tested with 1023 permutations. 1–Ucea, 2–Cârțișoara, 3–Porumbacu, 4–Avrig, 5–Bistra Mărului, 6–Bârzava, 7–Oituz, 8–Uz, 9–Craiului, 10–Caraș, 11–Topolog, 12–Sebeșel.
Fig. 3 Bayesian Structure results showing the individual membership of the samples from 12 Romanian rivers populations (K = 4). Each individual is represented by a vertical line partitioned into segments according to the proportion of the genome assigned to each of the identified clusters. Vertical line-the membership coefficient (Q), horizontal line-individuals. The populations are grouped from east to west as follows: (1) − Uz, (2) − Oituz, (3) − Ucea, (4) − Cârțișoara, (5) − Porumbacu, (6) − Avrig, (7) − Topolog, (8) − Sebeșel, (9) − Bistra Mărului, (10) − Craiului, (11) − Bârzava, (12) − Caraș. |
Fig. 4 PCoA (Principal Coordinate Analysis) for the 12 brown trout populations based on Nei's genetic distance (D) for the nine microsatellites. The two axes (Coord. 1 and Coord. 2) explain 11.72% and 6.28% of the variation, respectively. |
4 Discussion
4.1 Mitochondrial DNA analysis
The high haplotype diversity of the 12 analysed brown trout populations can be explained by the high nucleotide diversity within the D-loop mitochondrial region. Analysing the nucleotide position at which polymorphic sites were found, we observed that the most polymorphic region is at the 5′ end of the D-loop, but also close to the homopolimeric T region, which is situated after the first 400 bp of the D-loop as it has been reported in literature (Dunner et al., 2000).
The presence of the Danubian lineage might reflect the existence of distinct groups of ancient populations in this part of Europe (Kohout et al., 2013). Past cyclic glacial events and vast changes in the interconnectivity of the Black Sea with Caspian and Aral basins have also impacted the populations (Bernatchez, 2001). The easternmost (Uz and Oituz Rivers) and westernmost (Bârzava and Caraş Rivers) populations have the highest genetic diversity. This observation could be explained by the climate and relief variability on this Eastern and Western limits of the studied Romanian Carpathians (Posea, 2006). Furthermore, the populations belonging to the northern oriented mountainous slopes basins (except Cârţişoara population), have shown a lower genetic variability. The shorter lengths of the inner Carpathians Rivers are favourable lotic sectors for S. trutta and the significant human impact pressure in the downstream sectors (Bănăduc, 1999) can offer a second category of arguments in explaining the upper described special case of the inner Carpathian arch populations.
In the Romanian Carpathians there still are pure Danubian S. trutta populations. The post-glacial colonization (Hamilton et al., 1989) and peri-alpine dispersion of relict populations of pure Danubian lineage S. trutta suggest that they settled headwater streams post-glacially earlier than the appearance and spread of the Atlantic lineage (Lerceteau-Kohler et al., 2013). The presence of the Atlantic individuals can be a result of some chaotic stocking and restocking with S. trutta in the Romanian Carpathian streams and rivers in the 19th and 20th centuries situation which is continuing in the 21st century too. Few individuals from the Atlantic clade from stocked rivers were found in the southern Danube Serbian mountain areas, but not from the Mediterranean clade (Marić et al., 2006). We should take into consideration also that single sequences retrieved from GenBank can be assigned erroneous to a taxon, so this aspect should be carefully treated.
For the populations with exclusively Danubian lineage haplotypes, a significant overall genetic diversity was highlighted. This might be the effect of the relatively low human impact for the whole area of alpine headwater rivers or at least in the upstream lotic sectors. So, this area can be considered in this context as a regional genetic pool of remnant populations of the Danubian clade of S. trutta. The exceptions represented by the Porumbacu river populations can be an effect of more accentuated and more chaotic stockings and restocking activities in the past and present. Furthermore, in the north of the Romanian Carpathians, in Ukraine, there are no data about S. trutta individuals of other clades, except for the Danubian one. In addition, the natural spread of fish from Upper Danube in the Lower Danube basin may be improbable since 1972, the year when the Iron Gate dam on the Danube was finished (Bănăduc et al., 2014).
4.2 Microsatellites analysis
It is easily observed that the genetic clusters resulted from the analysis of the 12 populations are not entirely related to their geographic closeness. It was observed that the groups were formed by individuals coming from streams of the Meridional and Western Carpathians (Cârţişoara–Porumbacu–Avrig–Bârzava populations group, followed by the Bistra Mărului–Craiului–Caraş–Topolog–Sebeşel), with the exception of the streams from the Eastern Carpathians (the Uz–Oituz populations group) which showed a more homogenous genetic structure. Thus, the very complex genetic structure of these clusters may be a result of the existence of ancient refugees on the Romanian Carpathians that allowed later propagation with original individuals which had a genetic structure adapted to the specific environmental conditions.
The microsatellites data reveals signals of bottleneck and inbreeding events. Therefore, the patterns of reduced genetic variability in some of the analysed populations like Cartisoara, Porumbacu, Ucea and Avrig represent probably an artefact of the described past bottleneck effect. However, other factors could be contributing to patterns of genetic diversity and this includes the effect of gene flow. Connected populations are expected to show higher genetic diversity than isolated populations. Likewise, effective population size can also explain patterns in genetic diversity. Populations with low effective population sizes are expected to have lower genetic diversity due to stronger genetic drift (Freeland et al., 2011). Furthermore, lower diversity may suggest a population was recently formed by a small number of individuals (Allendorf et al., 2012).
The maximum number of alleles identified in the extreme East Uz River population (Trotuş River basin) can be an effect of the possible stockings and restocking activities in that river. Ucea, Porumbacu and Avrig form again a separate group of inner Carpathians arch rivers (Cârţişoara is again an exception from this group) this time due to their extremely low coefficient of inbreeding, their close geographical locations allowed spawners to haphazardly migrate from one stream to another. The base for this high degree of genetic differentiation among different S. trutta populations can be based on the geographical induced lack of connectivity. These resulted clusters of intra- and interrelationships suggest that the zoogeographic paths were eased by the relief and hydrography evolution along time.
Regarding the number of private alleles for each population, based on the total number of alleles, we observed that the highest number of private alleles is found in the Uz population, followed by the population of Bistra Mărului, while the populations of Porumbacu and Avrig show the lowest number of private alleles (Figs. 3 and 4). This result suggests that isolated populations might be important as genetic reservoirs of metapopulations (Linløkken et al., 2014), as suggested by the high genetic diversity of brown trout from several small sized streams. Concerning the average FIS index, all populations show positive values (p < 0.05 for 10 out of the 12 populations), suggesting an inbreeding event for the nine nuclear loci analysed in this study. However, by analysing the same values per locus and per populations, we observed that for the Str60 locus 10 out of the 12 populations showed a heterozygote excess, as also suggested by the negative not statistically significant FIS values.
Furthermore, the Nm parameter (number of migrants) has values >1 for the following population pairs (Tab. 3). This suggests that there are migrants between the analysed populations, and that they are not genetically isolated (Frankham et al., 2009). So, the individuals may have been part of an ancient, isolated metapopulation, from which further migration into large rivers and into their tributaries took place, or sporadic restocking activities with similar genetic structure individuals were done, or the current structure is an effect of local adaptations to similar habitats. Also, the highly FST values might be a consequence of the low genetic within population diversity.
However, a profound understanding of the phenomenon of genetic dispersion, implicitly the existence of a gene flow, is more difficult since its environmental implications are observed on a broader time scale compared to a short-term period in which preservative management measures of natural populations are taken (Palumbi, 2003). Thus, measuring the gene flow and the immigrants flow through the Nm parameter, closely related to the FST parameter, may prove difficult for species with high dispersal rates, as it might be the case for the brown trout.
At the same time, data on the process of colonization with brown trout are not detailed for Romania, which makes it difficult to understand the current structure and its causes of the natural trout populations. More information is known about the distribution of the Atlantic lineage that appears to have existed in the post-glacial era, since the northern part of the Atlantic was covered by ice (Apostolidis et al., 1996). However, the existence of trout populations in regions not covered by ice is due to the existence of isolated populations in refuges (Ferguson and Fleming, 1983). Studies focused on allozymes and DNA analysis have shown that there may have been several stages of European colonization with brown trout: from a Mediterranean-Caspian refuge, from an Iberian one or from a refuge in the proximity of the English Channel (Ferguson and Fleming, 1983; Simonović et al., 2017). The influence of the Black Sea Basin on this hypothesis is not excluded, but this is still a controversy about this complex colonization process (Weiss et al., 2000) through one of the most complex historical ichthyofauna areas of convergence in the Lower Danube-Danube Delta-North West Black Sea (Bănăduc et al., 2016).
5 Conclusions
To summarize, the shape and geographical orientation of the inner Carpathians can play a “protective/conservative” role in this situation. The easternmost and westernmost populations have a higher genetic variability by comparing with the populations belonging to the northern oriented mountainous slopes basins. The results suggests that the complex structure of the analysed brown trout populations could be explained by an “inner/outer” Carpathians pattern alongside the anthropic intervention and past events of bottleneck and inbreeding.
Supplementary Material
Table S1 The characteristics of the primers used for the amplification of nine microsatellite loci from S. trutta.
Table S2 The GenBank sequences used in the phylogenetic analysis.
Table S3 Polymorphic sites distribution within D-loop analyzed region. PoS − polymorphic sites, SS − singleton sites, PaS − parsimonic sites.
Table S4 The nuclear markers analyzed for the 12 brown trout populations. Na − number of alleles, HO − observed heterozygosity, HE − expected heterozygosity, PIC − polymorphic information content, HW − exact test for testing the deviation from Hardy-Weinberg proportions, *p < 0.05, NS − non significant, number of private alleles per population, M − mean Garza-Williamson index, SD − standard deviation.
Figures S1. The Delta K parameter values (vertical axis) for each K (horizontal axis).
Access hereAcknowledgement
This work was supported by the UEFISCDI, project 16PED/2017.
References
- Allendorf F, Luikart G, Aitken, S. 2012. Conservation and the Genetics of Populations. 2nd edn. West Sussex: Wiley-Blackwell Publishing. [Google Scholar]
- Apostolidis A, Karakousis Y, Triantaphyllidis C. 1996. Genetic differentiation and phylogenetic relationships among Greek Salmo trutta L. (brown trout) populations as revealed by RFLP analysis of PCR amplified mitochondrial DNA segments. Heredity 77: 608–618. [CrossRef] [PubMed] [Google Scholar]
- Bănăduc D. 1999. Data concerning the human impact on the ichthyofauna of the upper and middle sectors of the Olt River. Transylv Rev Syst Ecol Res 1: 157–164. [Google Scholar]
- Bănăduc D, Curtean-Bănăduc A, Lenhardt M, Guti G. 2014. Porţile de Fier/Iron Gates Gorges area (Danube) fish fauna. Transylv Rev Syst Ecol Res 16: 171–196. [Google Scholar]
- Bănăduc D, Rey S, Trichkova T, Lenhardt M, Curtean-Bănăduc A. 2016. The Lower Danube River-Danube Delta-North West Black Sea: a pivotal area of major interest for the past, present and future of its fish fauna − a short review. Sci Total Environ 545–546: 137–151. [PubMed] [Google Scholar]
- Behnke R. 2002. Trout and Salmon of North America. New York: Free Press Publishing. [Google Scholar]
- Berrebi P, Povž M, Jesenšek D, Cattaneo-Berrebi G, Crivelli AJ. 2000, The genetic diversity of native, stocked and hybrid populations of marble trout in the Soča river, Slovenia. Heredity 85: 277–287. [CrossRef] [PubMed] [Google Scholar]
- Bernatchez L. 2001. The evolutionary history of brown trout (Salmo trutta L.) inferred from phylogeographic, nested clade, and mismatch analyses of mitochondrial DNA variation. Evolution 55: 351–379. [CrossRef] [PubMed] [Google Scholar]
- Burghelea B, Bănăduc D, Curtean-Bănăduc A. 2013. The Timiş River basin (Banat, Romania) natural and antropogenic elements. a study case − management challenges. Transylv Rev Syst Ecol Res 15: 173–206. [Google Scholar]
- Carlsson J, Nilsson J. 2000. Population genetic structure of brown trout (Salmo trutta L.) within a northern boreal forest stream. Hereditas 132: 173–181. [CrossRef] [PubMed] [Google Scholar]
- Consuegra S, Elgan J, Verspoor E, Garcia de Leaniz C. 2015. Patterns of natural selection acting on the mitochondrial genome of a locally adapted fish species. Genet Sel Evol 47: 58. [Google Scholar]
- Cortey M, Garcia-Marin J. 2002. Evidence for phylogeographically informative sequence variation in the mitochondrial control region of Atlantic brown trout. J Fish Biol 60: 1058–1063. [Google Scholar]
- Crisp D. 2000. Trout and Salmon − Ecology, Conservation and Rehabilitation. London: Blackwell Science. [CrossRef] [Google Scholar]
- Cunjak R, Power G. 1986. Winter habitat utilization by stream resident brook trout (Salvelinus fontinalis) and brown trout (Salmo trutta). Can J Fish Aquat Sci 43: 1970–1981. [Google Scholar]
- Curtean-Bănăduc A, Pauli S, Bănăduc D, Didenko A, Sender J, Marić S, Del Monte P, Khoshnood Z, Zakeyuddin S. 2015. Environmental aspects of implementation of micro hydro power plants − a short review. Transylv Rev Syst Ecol Res 17: 179–198. [Google Scholar]
- Didenko A, Velykopolsky I, Buzevich I. 2011. Illegal fishing in the Tisza River drainage within Ukraine: a threat for local fish stocks? Arch Pol Fisheries 19: 249–257. [Google Scholar]
- Didenko A, Velykopolsky I, Chuklin A. 2014. Use of poachers' catches for studying fish fauna in the water bodies of the Transcarpathian region (Ukraine). Transylv Rev Syst Ecol Res 16: 87–101. [Google Scholar]
- Duftner N, Weiss S, Medgyesy N, Sturmbauer C. 2003. Enhanced phylogeographic information about Austrian brown trout populations derived from complete mitochondrial control region sequences. J Fish Biol 62: 427–435. [Google Scholar]
- Dunner S, Rayo L, Cañon J. 2000. Genetic structure in Atlantic brown trout (Salmo trutta L.) populations in the Iberian Peninsula: evidence from mitochondrial and nuclear DNA analysis. J Anim Breed Genet 117: 105–120. [Google Scholar]
- Earl D, Von Holdt B. 2012. Structure Harvester: a website and program for visualizing Structure output and implementing the Evanno method. Conserv Genet Resour 4: 359–361. [Google Scholar]
- Elliott J. 1994. Quantitative Ecology and the Brown Trout, In: R. May P. Harvey (Eds.), Oxford Series in Ecology and Evolution. New York: Oxford University Press, pp. 1–298. [Google Scholar]
- Excoffier L, Lischer H. 2010. Arlequin suite v3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour 10: 564–567. [Google Scholar]
- Ferguson A, Fleming C. 1983. Evolutionary and taxonomic significance of protein variation in the brown trout (Salmo trutta L.). In: G. Oxford D. Rollinson (Eds.), Polymorphism: Adaptive and Taxonomic Significance. London: Academic Press, pp. 86–99. [Google Scholar]
- Frankel O, Soule M. 1981. Conservation and Evolution. New York: Cambridge University Press. [Google Scholar]
- Frankham A, Ballou J, Briscoe D. 2009. Introduction to Conservation Genetics 2nd edn. Cambridge: Cambridge University Press. [Google Scholar]
- Freeland JR, Petersen SD, Kirk H. 2011. Molecular Ecology, 2nd edn. West Sussex: Wiley-Blackwell Publishing. [Google Scholar]
- Fumagalli L, Snoj A, Jesenšek D, Balloux F, Jug T, Duron O, Brossier F, Crivelli AJ, Berrebi P. 2002. Extreme genetic differentiation among the remnant populations of marble trout (Salmo marmoratus) in Slovenia. Mol Ecol 11: 2711–2716. [CrossRef] [PubMed] [Google Scholar]
- Garza J, Williamson E. 2001. Detection of reduction in population size using data from microsatellite loci. Mol Ecol 10: 305–318. [CrossRef] [PubMed] [Google Scholar]
- Hall T. 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser 41: 95–98. [Google Scholar]
- Hamilton K, Ferguson A, Taggart J, Tómasson T, Walker A, Fahy E. 1989. Post-glacial colonization of brown trout, Salmo trutta L.: Ldh-5 as a phylogeographic marker locus. J Fish Biol 35: 651–664. [Google Scholar]
- Hansen M. 2002. Estimating the long-term effects of stocking domesticated trout into wild brown trout (Salmo trutta) populations: an approach using microsatellite DNA analysis of historical and contemporary samples. Mol Ecol 11: 1003–1015. [CrossRef] [PubMed] [Google Scholar]
- Hilborn R, Quinn T, Schindler D, Rogers D. 2003. Biocomplexity and Fisheries Sustainability. Proc Natl Acad Sci 100: 6564–6568. [CrossRef] [Google Scholar]
- Klemetsen A, Amundsen P, Dempson J, Jonsson B, Jonsson N, O'Connell M, Mortensen E. 2003. Atlantic salmon Salmo salar L., brown trout Salmo trutta L. and Arctic charr Salvelinus alpinus (L.): a review of aspects of their life histories. Ecol Freshw Fish 12: 1–59. [Google Scholar]
- Kohout J, Jašková I, Papoušek I, Šedivá A, Šlechta V. 2012. Effects of stocking on the genetic structure of brown trout, Salmo trutta, in Central Europe inferred from mitochondrial and nuclear DNA markers. Fish Manag Ecol 19: 252–263. [Google Scholar]
- Kohout J, Šedivá A, Apostolou A, Stefanov T, Maric S, Gaffaroglu M, Šlechta V. 2013. Genetic diversity and phylogenetic origin of brown trout Salmo trutta populations in eastern Balkans. Biologia 68: 1229–1237. [Google Scholar]
- Lerceteau-Kohler E, Schliewen U, Kopun T, Weiss S. 2013. Genetic variation in brown trout Salmo trutta across the Danube, Rhine, and Elbe headwaters: a failure of the phylogeographic paradigm? BMC Evol Biol 13: 176. [CrossRef] [PubMed] [Google Scholar]
- Librado P, Rozas J. 2009. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25: 1451–1452. [CrossRef] [PubMed] [Google Scholar]
- Linløkken A, Johansen W, Wilson R. 2014. Genetic structure of brown trout, Salmo trutta, populations from differently sized tributaries of Lake Mjøsa in south-east Norway. Fish Manag Ecol 21: 515–525. [CrossRef] [Google Scholar]
- Liu K, Muse S. 2005. PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics 21: 2128–2129. [CrossRef] [PubMed] [Google Scholar]
- Marić S, Sušnik S, Simonovic P, Snoj A. 2006. Phylogeographic study of brown trout from Serbia, based on mitochondrial DNA control region analysis. Genet Sel Evol 38: 411–430. [Google Scholar]
- Meraner A, Baric S, Pelster B, Dalla Via J. 2007. Trout (Salmo trutta) mitochondrial DNA polymorphism in the centre of the marble trout distribution area. Hydrobiologia 579: 337–349. [Google Scholar]
- Nechifor R, Popa GO, Samu M, Dudu A, Banaduc D, Costache M, Georgescu SE. 2017. The genetic profiles of two salmonid populations from Romania obtained through nuclear marker analysis. Scientific Papers. Animal Sci Biotechnol 50: 74–78. [Google Scholar]
- Ostergren J, Nilsson J. 2012. Importance of life-history and landscape characteristics for genetic structure and genetic diversity of brown trout (Salmo trutta L.). Ecol Freshw Fish 21, 119–133. [Google Scholar]
- Osz A, Horvath A, Hoitsy G, Sipos DK, Keszte S, Safrany AJ, Maric S, Palko C, Toth B, Urbanyi B, Kovacs B. 2018. The genetic status of the Hungarian brown trout populations: exploration of a blind spot on the European map of Salmo trutta studies. PEERJ 6: e5152. [CrossRef] [PubMed] [Google Scholar]
- Palumbi R. 2003. Population genetics, demographic connectivity, and the design of marine reserves. Ecol Appl 13: S146–S158. [Google Scholar]
- Peakall R, Smouse P. 2006. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes 6: 288–295. [Google Scholar]
- Popa GO, Khalaf M, Dudu A, Curtean-Banaduc A, Banaduc D, Georgescu SE, Costache M. 2013. Genetic diversity of brown trout populations using mithocondrial markers in relatively similar geographical and ecological conditions − a Carpathian study. Transylv Rev Syst Ecol Res 15: 125–132. [Google Scholar]
- Popa GO, Curtean-Banaduc A, Banaduc D, Florescu IE, Burcea A, Dudu A, Georgescu SE, Costache M. 2016. Molecular markers reveal reduced genetic diversity in Romanian populations of Brown Trout, Salmo trutta L., 1758 (Salmonidae). Acta Zool Bulg 68: 399–406. [Google Scholar]
- Posea G. 2006. Geografia fizică a României. Bucharest: Fundaţia România de mâine Publishing House (in Romanian). [Google Scholar]
- Pritchard J, Stephens M, Donnelly P. 2000. Inference of population structure using multilocus genotype data. Genetics 155: 945–959. [PubMed] [Google Scholar]
- Pustovrh G, Snoj A, Susnik Bajec S. 2014. Molecular phylogeny of Salmo of the western Balkans, based upon multiple nuclear loci. Genet Sel Evol 46: 1–12. [Google Scholar]
- Reed D, Frankham R. 2003. Correlation between population fitness and genetic diversity. Conserv Biol 17: 230–237. [Google Scholar]
- Sanz N, Cortey M, Pla C, García-Marín J. 2006. Hatchery introgression blurs ancient hybridization between brown trout (Salmo trutta) lineages as indicated by complementary allozymes and mtDNA markers. Biol Conserv 130: 278–289. [Google Scholar]
- Schenekar T, Lerceteau-Kohler E, Weiss S. 2014. Fine-scale phylogeographic contact zone in Austrian brown trout Salmo trutta reveals multiple waves of post-glacial colonization and a pre-dominance of natural versus anthropogenic admixture. Conserv Genet 15: 561–572. [Google Scholar]
- Schindler D, Hilborn R, Chasco B, Boatright CP, Quinn T, Rogers L, Webster M. 2010. Population diversity and the portfolio effect in an exploited species. Nature 465: 609–612. [CrossRef] [PubMed] [Google Scholar]
- Schwartz M, Luikart G, Waples R. 2007. Genetic monitoring as a promising tool for conservation and management. Trends Ecol Evol 22: 25–33. [CrossRef] [EDP Sciences] [PubMed] [Google Scholar]
- Simonović P, Tošić A, Dubravka S, Nikolić V, Piria M, Tomljanović T, Sprem N, Mrdak D, Milošević D, Beciraj A, Dekic R, Povz M. 2017. Diversity of Brown trout Salmo trutta (L.) in the River Danube basin of Western Balkans as assessed from the structure of their mitochondrial control region haplotypes. J Ichthyol 57: 603–616. [CrossRef] [Google Scholar]
- Splendiani A, Giovannotti M, Cerioni PN, Caniglia ML, Caputo V. 2006. Phylogeographic inferences on the native brown trout mtDNA variation in central Italy. Ital J Zool 73: 179–189. [CrossRef] [Google Scholar]
- Taggart J, Hynes R, Prodohl P, Ferguson A. 1991. A simplified protocol for routine total DNA isolation from salmonid fishes. J Fish Biol 4: 963–965. [Google Scholar]
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. 2013. MEGA6: Molecular Evolutionary Genetics Analysis, version 6.0. Mol Biol Evol 30: 2725–2729. [CrossRef] [PubMed] [Google Scholar]
- Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. 2004. MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes 4: 535–538. [Google Scholar]
- Warren M, Dunbar M, Smith C. 2015. River flow as a determinant of salmonid distribution and abundance: a review. Environ Biol Fishes 98: 1695–1717. [Google Scholar]
- Was A, Wenne R. 2003. Microsatellite DNA polymorphism in intensely enhanced populations of sea trout (Salmo trutta) in the Southern Baltic. Mar Biotechnol 5: 234–243. [Google Scholar]
- Weiss S, Antunes A, Schlötterer C, Alexandrino P. 2000. Mitochondrial haplotype diversity among Portuguese brown trout Salmo trutta L. populations: relevance to the post-Pleistocene recolonization of northern Europe. Mol Ecol 9: 691–698. [CrossRef] [PubMed] [Google Scholar]
- Wenne R. 2001. Mitochondrial DNA variation in sea trout from coastal rivers in the southern Baltic region. ICES J Mar Sci 58: 230–237. [Google Scholar]
- Wenne R, Bernas R, Pocwierz-Kotus A, Drywa A, Was A. 2016. Recent genetic changes in enhanced populations of sea trout (Salmo trutta m. trutta) in the southern Baltic rivers revealed with SNP analysis. Aquat Living Resour 29: 103. [CrossRef] [EDP Sciences] [Google Scholar]
- Wlodarczyk E, Wenne R. 2001. Mitochondrial DNA variation in sea trout from coastal rivers in the southern Baltic region. ICES J Mar Sci 58: 230–237. [Google Scholar]
Internet resources
- Allendorf F, Luikart G, Aitken, S. 2012. Conservation and the Genetics of Populations. 2nd edn. West Sussex: Wiley-Blackwell Publishing. [Google Scholar]
- Apostolidis A, Karakousis Y, Triantaphyllidis C. 1996. Genetic differentiation and phylogenetic relationships among Greek Salmo trutta L. (brown trout) populations as revealed by RFLP analysis of PCR amplified mitochondrial DNA segments. Heredity 77: 608–618. [CrossRef] [PubMed] [Google Scholar]
- Bănăduc D. 1999. Data concerning the human impact on the ichthyofauna of the upper and middle sectors of the Olt River. Transylv Rev Syst Ecol Res 1: 157–164. [Google Scholar]
- Bănăduc D, Curtean-Bănăduc A, Lenhardt M, Guti G. 2014. Porţile de Fier/Iron Gates Gorges area (Danube) fish fauna. Transylv Rev Syst Ecol Res 16: 171–196. [Google Scholar]
- Bănăduc D, Rey S, Trichkova T, Lenhardt M, Curtean-Bănăduc A. 2016. The Lower Danube River-Danube Delta-North West Black Sea: a pivotal area of major interest for the past, present and future of its fish fauna − a short review. Sci Total Environ 545–546: 137–151. [PubMed] [Google Scholar]
- Behnke R. 2002. Trout and Salmon of North America. New York: Free Press Publishing. [Google Scholar]
- Berrebi P, Povž M, Jesenšek D, Cattaneo-Berrebi G, Crivelli AJ. 2000, The genetic diversity of native, stocked and hybrid populations of marble trout in the Soča river, Slovenia. Heredity 85: 277–287. [CrossRef] [PubMed] [Google Scholar]
- Bernatchez L. 2001. The evolutionary history of brown trout (Salmo trutta L.) inferred from phylogeographic, nested clade, and mismatch analyses of mitochondrial DNA variation. Evolution 55: 351–379. [CrossRef] [PubMed] [Google Scholar]
- Burghelea B, Bănăduc D, Curtean-Bănăduc A. 2013. The Timiş River basin (Banat, Romania) natural and antropogenic elements. a study case − management challenges. Transylv Rev Syst Ecol Res 15: 173–206. [Google Scholar]
- Carlsson J, Nilsson J. 2000. Population genetic structure of brown trout (Salmo trutta L.) within a northern boreal forest stream. Hereditas 132: 173–181. [CrossRef] [PubMed] [Google Scholar]
- Consuegra S, Elgan J, Verspoor E, Garcia de Leaniz C. 2015. Patterns of natural selection acting on the mitochondrial genome of a locally adapted fish species. Genet Sel Evol 47: 58. [Google Scholar]
- Cortey M, Garcia-Marin J. 2002. Evidence for phylogeographically informative sequence variation in the mitochondrial control region of Atlantic brown trout. J Fish Biol 60: 1058–1063. [Google Scholar]
- Crisp D. 2000. Trout and Salmon − Ecology, Conservation and Rehabilitation. London: Blackwell Science. [CrossRef] [Google Scholar]
- Cunjak R, Power G. 1986. Winter habitat utilization by stream resident brook trout (Salvelinus fontinalis) and brown trout (Salmo trutta). Can J Fish Aquat Sci 43: 1970–1981. [Google Scholar]
- Curtean-Bănăduc A, Pauli S, Bănăduc D, Didenko A, Sender J, Marić S, Del Monte P, Khoshnood Z, Zakeyuddin S. 2015. Environmental aspects of implementation of micro hydro power plants − a short review. Transylv Rev Syst Ecol Res 17: 179–198. [Google Scholar]
- Didenko A, Velykopolsky I, Buzevich I. 2011. Illegal fishing in the Tisza River drainage within Ukraine: a threat for local fish stocks? Arch Pol Fisheries 19: 249–257. [Google Scholar]
- Didenko A, Velykopolsky I, Chuklin A. 2014. Use of poachers' catches for studying fish fauna in the water bodies of the Transcarpathian region (Ukraine). Transylv Rev Syst Ecol Res 16: 87–101. [Google Scholar]
- Duftner N, Weiss S, Medgyesy N, Sturmbauer C. 2003. Enhanced phylogeographic information about Austrian brown trout populations derived from complete mitochondrial control region sequences. J Fish Biol 62: 427–435. [Google Scholar]
- Dunner S, Rayo L, Cañon J. 2000. Genetic structure in Atlantic brown trout (Salmo trutta L.) populations in the Iberian Peninsula: evidence from mitochondrial and nuclear DNA analysis. J Anim Breed Genet 117: 105–120. [Google Scholar]
- Earl D, Von Holdt B. 2012. Structure Harvester: a website and program for visualizing Structure output and implementing the Evanno method. Conserv Genet Resour 4: 359–361. [Google Scholar]
- Elliott J. 1994. Quantitative Ecology and the Brown Trout, In: R. May P. Harvey (Eds.), Oxford Series in Ecology and Evolution. New York: Oxford University Press, pp. 1–298. [Google Scholar]
- Excoffier L, Lischer H. 2010. Arlequin suite v3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour 10: 564–567. [Google Scholar]
- Ferguson A, Fleming C. 1983. Evolutionary and taxonomic significance of protein variation in the brown trout (Salmo trutta L.). In: G. Oxford D. Rollinson (Eds.), Polymorphism: Adaptive and Taxonomic Significance. London: Academic Press, pp. 86–99. [Google Scholar]
- Frankel O, Soule M. 1981. Conservation and Evolution. New York: Cambridge University Press. [Google Scholar]
- Frankham A, Ballou J, Briscoe D. 2009. Introduction to Conservation Genetics 2nd edn. Cambridge: Cambridge University Press. [Google Scholar]
- Freeland JR, Petersen SD, Kirk H. 2011. Molecular Ecology, 2nd edn. West Sussex: Wiley-Blackwell Publishing. [Google Scholar]
- Fumagalli L, Snoj A, Jesenšek D, Balloux F, Jug T, Duron O, Brossier F, Crivelli AJ, Berrebi P. 2002. Extreme genetic differentiation among the remnant populations of marble trout (Salmo marmoratus) in Slovenia. Mol Ecol 11: 2711–2716. [CrossRef] [PubMed] [Google Scholar]
- Garza J, Williamson E. 2001. Detection of reduction in population size using data from microsatellite loci. Mol Ecol 10: 305–318. [CrossRef] [PubMed] [Google Scholar]
- Hall T. 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser 41: 95–98. [Google Scholar]
- Hamilton K, Ferguson A, Taggart J, Tómasson T, Walker A, Fahy E. 1989. Post-glacial colonization of brown trout, Salmo trutta L.: Ldh-5 as a phylogeographic marker locus. J Fish Biol 35: 651–664. [Google Scholar]
- Hansen M. 2002. Estimating the long-term effects of stocking domesticated trout into wild brown trout (Salmo trutta) populations: an approach using microsatellite DNA analysis of historical and contemporary samples. Mol Ecol 11: 1003–1015. [CrossRef] [PubMed] [Google Scholar]
- Hilborn R, Quinn T, Schindler D, Rogers D. 2003. Biocomplexity and Fisheries Sustainability. Proc Natl Acad Sci 100: 6564–6568. [CrossRef] [Google Scholar]
- Klemetsen A, Amundsen P, Dempson J, Jonsson B, Jonsson N, O'Connell M, Mortensen E. 2003. Atlantic salmon Salmo salar L., brown trout Salmo trutta L. and Arctic charr Salvelinus alpinus (L.): a review of aspects of their life histories. Ecol Freshw Fish 12: 1–59. [Google Scholar]
- Kohout J, Jašková I, Papoušek I, Šedivá A, Šlechta V. 2012. Effects of stocking on the genetic structure of brown trout, Salmo trutta, in Central Europe inferred from mitochondrial and nuclear DNA markers. Fish Manag Ecol 19: 252–263. [Google Scholar]
- Kohout J, Šedivá A, Apostolou A, Stefanov T, Maric S, Gaffaroglu M, Šlechta V. 2013. Genetic diversity and phylogenetic origin of brown trout Salmo trutta populations in eastern Balkans. Biologia 68: 1229–1237. [Google Scholar]
- Lerceteau-Kohler E, Schliewen U, Kopun T, Weiss S. 2013. Genetic variation in brown trout Salmo trutta across the Danube, Rhine, and Elbe headwaters: a failure of the phylogeographic paradigm? BMC Evol Biol 13: 176. [CrossRef] [PubMed] [Google Scholar]
- Librado P, Rozas J. 2009. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25: 1451–1452. [CrossRef] [PubMed] [Google Scholar]
- Linløkken A, Johansen W, Wilson R. 2014. Genetic structure of brown trout, Salmo trutta, populations from differently sized tributaries of Lake Mjøsa in south-east Norway. Fish Manag Ecol 21: 515–525. [CrossRef] [Google Scholar]
- Liu K, Muse S. 2005. PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics 21: 2128–2129. [CrossRef] [PubMed] [Google Scholar]
- Marić S, Sušnik S, Simonovic P, Snoj A. 2006. Phylogeographic study of brown trout from Serbia, based on mitochondrial DNA control region analysis. Genet Sel Evol 38: 411–430. [Google Scholar]
- Meraner A, Baric S, Pelster B, Dalla Via J. 2007. Trout (Salmo trutta) mitochondrial DNA polymorphism in the centre of the marble trout distribution area. Hydrobiologia 579: 337–349. [Google Scholar]
- Nechifor R, Popa GO, Samu M, Dudu A, Banaduc D, Costache M, Georgescu SE. 2017. The genetic profiles of two salmonid populations from Romania obtained through nuclear marker analysis. Scientific Papers. Animal Sci Biotechnol 50: 74–78. [Google Scholar]
- Ostergren J, Nilsson J. 2012. Importance of life-history and landscape characteristics for genetic structure and genetic diversity of brown trout (Salmo trutta L.). Ecol Freshw Fish 21, 119–133. [Google Scholar]
- Osz A, Horvath A, Hoitsy G, Sipos DK, Keszte S, Safrany AJ, Maric S, Palko C, Toth B, Urbanyi B, Kovacs B. 2018. The genetic status of the Hungarian brown trout populations: exploration of a blind spot on the European map of Salmo trutta studies. PEERJ 6: e5152. [CrossRef] [PubMed] [Google Scholar]
- Palumbi R. 2003. Population genetics, demographic connectivity, and the design of marine reserves. Ecol Appl 13: S146–S158. [Google Scholar]
- Peakall R, Smouse P. 2006. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes 6: 288–295. [Google Scholar]
- Popa GO, Khalaf M, Dudu A, Curtean-Banaduc A, Banaduc D, Georgescu SE, Costache M. 2013. Genetic diversity of brown trout populations using mithocondrial markers in relatively similar geographical and ecological conditions − a Carpathian study. Transylv Rev Syst Ecol Res 15: 125–132. [Google Scholar]
- Popa GO, Curtean-Banaduc A, Banaduc D, Florescu IE, Burcea A, Dudu A, Georgescu SE, Costache M. 2016. Molecular markers reveal reduced genetic diversity in Romanian populations of Brown Trout, Salmo trutta L., 1758 (Salmonidae). Acta Zool Bulg 68: 399–406. [Google Scholar]
- Posea G. 2006. Geografia fizică a României. Bucharest: Fundaţia România de mâine Publishing House (in Romanian). [Google Scholar]
- Pritchard J, Stephens M, Donnelly P. 2000. Inference of population structure using multilocus genotype data. Genetics 155: 945–959. [PubMed] [Google Scholar]
- Pustovrh G, Snoj A, Susnik Bajec S. 2014. Molecular phylogeny of Salmo of the western Balkans, based upon multiple nuclear loci. Genet Sel Evol 46: 1–12. [Google Scholar]
- Reed D, Frankham R. 2003. Correlation between population fitness and genetic diversity. Conserv Biol 17: 230–237. [Google Scholar]
- Sanz N, Cortey M, Pla C, García-Marín J. 2006. Hatchery introgression blurs ancient hybridization between brown trout (Salmo trutta) lineages as indicated by complementary allozymes and mtDNA markers. Biol Conserv 130: 278–289. [Google Scholar]
- Schenekar T, Lerceteau-Kohler E, Weiss S. 2014. Fine-scale phylogeographic contact zone in Austrian brown trout Salmo trutta reveals multiple waves of post-glacial colonization and a pre-dominance of natural versus anthropogenic admixture. Conserv Genet 15: 561–572. [Google Scholar]
- Schindler D, Hilborn R, Chasco B, Boatright CP, Quinn T, Rogers L, Webster M. 2010. Population diversity and the portfolio effect in an exploited species. Nature 465: 609–612. [CrossRef] [PubMed] [Google Scholar]
- Schwartz M, Luikart G, Waples R. 2007. Genetic monitoring as a promising tool for conservation and management. Trends Ecol Evol 22: 25–33. [CrossRef] [EDP Sciences] [PubMed] [Google Scholar]
- Simonović P, Tošić A, Dubravka S, Nikolić V, Piria M, Tomljanović T, Sprem N, Mrdak D, Milošević D, Beciraj A, Dekic R, Povz M. 2017. Diversity of Brown trout Salmo trutta (L.) in the River Danube basin of Western Balkans as assessed from the structure of their mitochondrial control region haplotypes. J Ichthyol 57: 603–616. [CrossRef] [Google Scholar]
- Splendiani A, Giovannotti M, Cerioni PN, Caniglia ML, Caputo V. 2006. Phylogeographic inferences on the native brown trout mtDNA variation in central Italy. Ital J Zool 73: 179–189. [CrossRef] [Google Scholar]
- Taggart J, Hynes R, Prodohl P, Ferguson A. 1991. A simplified protocol for routine total DNA isolation from salmonid fishes. J Fish Biol 4: 963–965. [Google Scholar]
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. 2013. MEGA6: Molecular Evolutionary Genetics Analysis, version 6.0. Mol Biol Evol 30: 2725–2729. [CrossRef] [PubMed] [Google Scholar]
- Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. 2004. MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes 4: 535–538. [Google Scholar]
- Warren M, Dunbar M, Smith C. 2015. River flow as a determinant of salmonid distribution and abundance: a review. Environ Biol Fishes 98: 1695–1717. [Google Scholar]
- Was A, Wenne R. 2003. Microsatellite DNA polymorphism in intensely enhanced populations of sea trout (Salmo trutta) in the Southern Baltic. Mar Biotechnol 5: 234–243. [Google Scholar]
- Weiss S, Antunes A, Schlötterer C, Alexandrino P. 2000. Mitochondrial haplotype diversity among Portuguese brown trout Salmo trutta L. populations: relevance to the post-Pleistocene recolonization of northern Europe. Mol Ecol 9: 691–698. [CrossRef] [PubMed] [Google Scholar]
- Wenne R. 2001. Mitochondrial DNA variation in sea trout from coastal rivers in the southern Baltic region. ICES J Mar Sci 58: 230–237. [Google Scholar]
- Wenne R, Bernas R, Pocwierz-Kotus A, Drywa A, Was A. 2016. Recent genetic changes in enhanced populations of sea trout (Salmo trutta m. trutta) in the southern Baltic rivers revealed with SNP analysis. Aquat Living Resour 29: 103. [CrossRef] [EDP Sciences] [Google Scholar]
- Wlodarczyk E, Wenne R. 2001. Mitochondrial DNA variation in sea trout from coastal rivers in the southern Baltic region. ICES J Mar Sci 58: 230–237. [Google Scholar]
Cite this article as: Popa G-O, Dudu A, Bănăduc D, Curtean-Bănăduc A, Burcea A, Ureche D, Nechifor R, Georgescu SE, Costache M. 2019. Genetic analysis of populations of brown trout (Salmo trutta L.) from the Romanian Carpathians. Aquat. Living Resour. 32: 23
All Tables
The number of D-loop sequences along with the number of haplotypes and genotypes obtained for each population. The abbreviation for each population name is in the brackets.
The inbreeding coefficient (FIS) values for the nine loci. *p < 0.05, tested with 1000 bootstrap replications. Abbreviated are used for the sampled populations.
FST values (above the diagonal) and Nm values (under the diagonal) per population. *p < 0.05, tested with 1023 permutations. 1–Ucea, 2–Cârțișoara, 3–Porumbacu, 4–Avrig, 5–Bistra Mărului, 6–Bârzava, 7–Oituz, 8–Uz, 9–Craiului, 10–Caraș, 11–Topolog, 12–Sebeșel.
All Figures
Fig. 1 Map of locations. (A) Sampling sites: Uz − Uz River, Otz − Oituz River, Uc − Ucea River, Por–Porumbacu River, Avg − Avrig River, Crt − Cârțișoara River, Top − Topolog River, BM − Bistra Mărului River, Sb − Sebeşel River, Cr − Craiului River, Bz − Bârzava River, Car − Caraş River. Rivers Porumbacu, Avrig, Cârțișoara, Ucea, Sebeşel discharge to the Northern/inner slope, while Toplog River discharge to the Southern/outer slope of the mountain ridge. (B) Map of region showing the locations of 12 sampling sites. |
|
In the text |
Fig. 2 Phylogenetic tree of D-loop sequences obtained with the Maximum Likelihood method based on the Tamura 3-parameter model with 1000 bootstrap replicates. The bootstrap values higher than 50 are displayed. Green circles − Danubian clade, yellow triangles − Atlantic clade, pink squares − Mediterranean clade, blue diamonds − Adriatic clade, red triangles − S. marmoratus. |
|
In the text |
Fig. 3 Bayesian Structure results showing the individual membership of the samples from 12 Romanian rivers populations (K = 4). Each individual is represented by a vertical line partitioned into segments according to the proportion of the genome assigned to each of the identified clusters. Vertical line-the membership coefficient (Q), horizontal line-individuals. The populations are grouped from east to west as follows: (1) − Uz, (2) − Oituz, (3) − Ucea, (4) − Cârțișoara, (5) − Porumbacu, (6) − Avrig, (7) − Topolog, (8) − Sebeșel, (9) − Bistra Mărului, (10) − Craiului, (11) − Bârzava, (12) − Caraș. |
|
In the text |
Fig. 4 PCoA (Principal Coordinate Analysis) for the 12 brown trout populations based on Nei's genetic distance (D) for the nine microsatellites. The two axes (Coord. 1 and Coord. 2) explain 11.72% and 6.28% of the variation, respectively. |
|
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.