Mitochondrial haplotype variation and phylogeography of Iberian brown trout populations

Regional environmental and climatic concerns on preserving native gene pools of a least concern species: Brown trout lineages in Mediterranean streams

The European brown trout, Salmo trutta, is a cold-adapted fish reported as a Least Concern species in the IUCN Red List. This species colonized new territories from southern refuges during the last glacial melting, but during the 20th century suffered from anthropic impacts on its habitats. The long-time survival of the species relies on the genetic diversity within and among populations. Brown trout is among the genetically most diverse vertebrate species; however, native populations in Mediterranean rivers have dramatically suffered of introgressive hybridization from extensive releases of evolutionary distant non-native Atlantic stocks. In addition, in Mediterranean rivers climate change will result in unsuitable conditions for the species during the 21st century. Using brown trout populations at the headstreams of a Pyrenean river as a model, this paper revised how hatchery releases have affected the native gene pools and how environmental and climatic variables controlled the amount of local introgression at intra-basin level. Introgressive hybridization was detected in all studied sites. Ten times larger divergence was observed among populations at tributaries than among populations along the main stem. A highly impacted population distributed in a long transect in the main stem suggested that hatchery fish move towards the main stem wherever released. From already highly impacted populations and despite the cessation of hatchery releases, warmer temperatures and lower precipitation expected from climate change will extend the introgressive hybridization along the basin, contributing to the extinction of the native gene pools. Based on available morphological distinction of native, hatchery and hybrid brown trout, we advocate the involvement of regional social groups (e.g. riverside dwellers, anglers, conservationists, hikers) in citizen science programs to detect the spread of non-native phenotypes along the rivers. These are cheap and fast methods to collaborate with fishery managers in the preservation and recovery of the regional native populations.

Genomic Hatchery Introgression in Brown Trout (Salmo trutta L.): Development of a Diagnostic SNP Panel for Monitoring the Impacted Mediterranean Rivers

Brown trout (Salmo trutta L.) populations have been restocked during recent decades to satisfy angling demand and counterbalance the decline of wild populations. Millions of fertile brown trout individuals were released into Mediterranean and Atlantic rivers from hatcheries with homogeneous central European stocks. Consequently, many native gene pools have become endangered by introgressive hybridization with those hatchery stocks. Different genetic tools have been used to identify and evaluate the degree of introgression starting from pure native and restocking reference populations (e.g., LDH-C* locus, microsatellites). However, due to the high genetic structuring of brown trout, the definition of the "native pool" is hard to achieve. Additionally, although the LDH-C* locus is useful for determining the introgression degree at the population level, its consistency at individual level is far from being accurate, especially after several generations were since releases. Accordingly, the development of a more powerful and cost-effective tool is essential for an appropriate monitoring to recover brown-trout-native gene pools. Here, we used the 2b restriction site-associated DNA sequencing (2b-RADseq) and Stacks 2 with a reference genome to identify single-nucleotide polymorphisms (SNPs) diagnostic for hatchery-native fish discrimination in the Atlantic and Mediterranean drainages of the Iberian Peninsula. A final set of 20 SNPs was validated in a MassARRAY® System genotyping by contrasting data with the whole SNP dataset using samples with different degree of introgression from those previously recorded. Heterogeneous introgression impact was confirmed among and within river basins, and was the highest in the Mediterranean Slope. The SNP tool reported here should be assessed in a broader sample scenario in Southern Europe considering its potential for monitoring recovery plans.

A genomic perspective on an old question: Salmo trouts or Salmo trutta (Teleostei: Salmonidae)?

There are particular challenges in defining the taxonomic status of recently radiated groups due to the low level of phylogenetic signal. Members of the Salmo trutta species-complex, which mostly evolved during and following the Pleistocene, show high morphological and ecological diversity that, along with their very wide geographic distribution, have led to morphological description of 47 extant nominal species. However, many of these species have not been supported by previous phylogenetic studies, which could be partly due to lack of significant genetic differences among them, the limited resolution offered by molecular methods previously used, as well as the often local scale of these studies. The development of next-generation sequencing (NGS) and related analytical tools have enhanced our ability to address such challenging questions. In this study, Genotyping-by-Sequencing (GBS) of 15,169 filtered SNPs and mitochondrial DNA (mtDNA) D-loop sequences were combined to assess the phylogenetic relationships among 166 brown trouts representing 21 described species and three undescribed groups collected from 84 localities throughout their natural distribution in Europe, west Asia, and North Africa. The data were analysed using different clustering algorithms (admixture analysis and discriminant analysis of principal components-DAPC), a Bayes Factor Delimitation (BFD) test, species tree reconstruction, gene flow tests (three- and four-population tests), and Rogue taxa identification tests. Genomic contributions of the Atlantic lineage brown trout were found in all major sea basins excluding the North African and Aral Sea basins, suggesting introgressive hybridization of native brown trouts driven by stocking using strains of the Atlantic lineage. After removing the phylogenetic noise caused by the Atlantic brown trout, admixture clusters and DAPC clustering based on GBS data, respectively, resolved 11 and 13 clusters among the previously described brown trout species, which were also supported by BFD test results. Our results suggest that natural hybridization between different brown trout lineages has probably played an important role in the origin of several of the putative species, including S. marmoratus, S. carpio, S. farioides, S. pellegrini, S. caspius (in the Kura River drainage) and Salmo sp. in the Danube River basin. Overall, our results support a multi-species taxonomy for brown trouts. They also resolve some species in the Adriatic-Mediterranean and Black Sea drainages as members of very closely related genomic clusters that may need taxonomic revision. However, any final conclusions pertaining to the taxonomy of the brown trout complex should be based on an integrative approach combining genomic, morphological, and ecological data. To avoid challenges in taxonomy and conservation of species complexes like brown trouts, it is suggested to describe species based on genomic clusters of populations instead of describing species based only on morphologically differentiated single type populations.

Salmo macrostigma (Teleostei, Salmonidae): Nothing more than a brown trout (S. trutta) lineage?

We examined specimens of the macrostigma trout Salmo macrostigma, which refers to big black spots on the flanks, to assess whether it is an example of taxonomic inflation within the brown trout Salmo trutta complex. Using new specimens, publicly available data and a mitogenomic protocol to amplify the control and cytochrome b regions of the mitochondrial genome from degraded museum samples, including one syntype specimen, the present study shows that the macrostigma trout is not a valid species. Our results suggest the occurrence of a distinct evolutionary lineage of S. trutta in North Africa and Sicily. The name of the North African lineage is proposed for this lineage, which was found to be sister to the Atlantic lineage of brown trout, S. trutta.

Molecular phylogeny of Salmo of the western Balkans, based upon multiple nuclear loci

Background

Classification of species within the genus Salmo is still a matter of discussion due to their high level of diversity and to the low power of resolution of mitochondrial (mt)DNA-based phylogeny analyses that have been traditionally used in evolutionary studies of the genus. We apply a new marker system based on nuclear (n)DNA loci to present a novel view of the phylogeny of Salmo representatives and we compare it with the mtDNA-based phylogeny.

Methods

Twenty-two nDNA loci were sequenced for 76 individuals of the brown trout complex: Salmo trutta (Danubian, Atlantic, Adriatic, Mediterranean and Duero mtDNA lineages), Salmo marmoratus (marble trout), Salmo obtusirostris (softmouth trout), and Salmo ohridanus (Ohrid belvica or belushka). Sequences were phylogenetically analyzed using maximum-likelihood and Bayesian Inference methods. The divergence time of the major clades was estimated using the program BEAST.

Results

The existence of five genetic units i.e. S. salar, S. ohridanus, S. obtusirostris, S. marmoratus and the S. trutta complex, including its major phylogenetic lineages was confirmed. Contrary to previous observations, S. obtusirostris was found to be sister to the S. trutta complex and the S. marmoratus clade rather than to the S. ohridanus clade. Reticulate evolution of S. obtusirostris was confirmed and a time for its pre-glacial origin suggested. S. marmoratus was found to be a separate species as S. trutta and S. obtusirostris. Relationships among lineages within the S. trutta complex were weakly supported and remain largely unresolved.

Conclusions

Nuclear DNA-based results showed a fairly good match with the phylogeny of Salmo inferred from mtDNA analyses. The comparison of nDNA and mtDNA data revealed at least four cases of mitochondrial–nuclear DNA discordance observed that were all confined to the Adriatic basin of the Western Balkans. Together with the well-known extensive morphological and genetic variability of Balkan trouts, this observation highlights an interesting and variegated evolutionary history of Salmo in this area.

Restricted gene flow at the micro- and macro-geographical scale in marble trout based on mtDNA and microsatellite polymorphism

BACKGROUND

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.

RESULTS

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.

CONCLUSION

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.

Across the great divide: genetic forensics reveals misidentification of endangered cutthroat trout populations

Accurate assessment of species identity is fundamental for conservation biology. Using molecular markers from the mitochondrial and nuclear genomes, we discovered that many putatively native populations of greenback cutthroat trout (Oncorhynchus clarkii stomias) comprised another subspecies of cutthroat trout, Colorado River cutthroat trout (Oncorhynchus clarkii pleuriticus). The error can be explained by the introduction of Colorado River cutthroat trout throughout the native range of greenback cutthroat trout in the late 19th and early 20th centuries by fish stocking activities. Our results suggest greenback cutthroat trout within its native range is at a higher risk of extinction than ever before despite conservation activities spanning more than two decades.

An evaluation of the effects of conservation and fishery enhancement hatcheries on wild populations of salmon

The historical, political and scientific aspects of salmon hatchery programmes designed to enhance fishery production, or to recover endangered populations, are reviewed. We start by pointing out that the establishment of hatcheries has been a political response to societal demands for harvest and conservation; given this social context, we then critically examined the levels of activity, the biological risks, and the economic analysis associated with salmon hatchery programmes. A rigorous analysis of the impacts of hatchery programmes was hindered by the lack of standardized data on release sizes and survival rates at all ecological scales, and since hatchery programme objectives are rarely defined, it was also difficult to measure their effectiveness at meeting release objectives. Debates on the genetic effects of hatchery programmes on wild fish have been dominated by whether correct management practices can reduce negative outcomes, but we noted that there has been an absence of programmatic research approaches addressing this important issue. Competitive interactions between hatchery and wild fish were observed to be complex, but studies researching approaches to reduce these interactions at all ecological scales during the entire salmon life history have been rare, and thus are not typically considered in hatchery management. Harvesting of salmon released from fishery enhancement hatcheries likely impacts vulnerable wild populations; managers have responded to this problem by mass marking hatchery fish, so that fishing effort can be directed towards hatchery populations. However, we noted that the effectiveness of this approach is dependant on accurate marking and production of hatchery fish with high survival rates, and it is not yet clear whether selective fishing will prevent overharvest of wild populations. Finally, research demonstrating disease transmission from hatchery fish to wild populations was observed to be equivocal; evidence in this area has been constrained by the lack of effective approaches to studying the fate of pathogens in the wild. We then reviewed several approaches to studying the economic consequences of hatchery activities intended to inform the social decisions surrounding programmes, but recognized that placing monetary value on conservation efforts or on hatcheries that mitigate cultural groups' loss of historical harvest opportunities may complicate these analyses. We noted that economic issues have rarely been included in decision making on hatchery programmes. We end by identifying existing major knowledge gaps, which, if filled, could contribute towards a fuller understanding of the role that hatchery programmes could play in meeting divergent goals. However, we also recognized that many management recommendations arising from such research may involve trade-offs between different risks, and that decisions about these trade-offs must occur within a social context. Hatcheries have played an important role in sustaining some highly endangered populations, and it is possible that reform of practices will lead to an increase in the number of successful programmes. However, a serious appraisal of the role of hatcheries in meeting broader needs is urgently warranted and should take place at the scientific, but more effectively, at the societal level.

Life on the edge: the long-term persistence and contrasting spatial genetic structure of distinct brown trout life histories at their ecological limits

The Atlantic Iberian brown trout is at the southwestern limit of its distribution. At this ecological edge, which was once a glacial refugia, anadromy becomes less common as increased water temperatures restricted populations closer to the headwaters. We examined 847 individuals from 20 populations from throughout this region and assessed spatial genetic structure using 11 protein and four microsatellite loci. The higher levels of heterozygosity and allelic diversity north of the southernmost limit of anadromy (SLA), as well as an isolation-by-distance model of population structure, likely influenced by the anadromous forms, suggest that more stable demographic conditions existed over time in this region. Populations south of the SLA were highly differentiated given the restricted size of the area (protein F(ST) = 0.16 in the north and 0.63 in the south of the SLA; microsatellite F(ST) = 0.18 in the north and 0.70 in the south of the SLA). The low levels of heterozygosity and the pattern of southward allele depletion in resident populations is indicative of fragmentation, caused by stressful ecological conditions that reduced the anadromy (restricting gene flow) and the effective population sizes (higher genetic drift), which, in combination, dramatically decreased within-population genetic variation and increased among-population genetic variation. The higher gene diversity north of the SLA does not reflect ancestry but rather the signature of a population size expansion, as evidence suggest the persistence of older populations (with several private alleles) south of the SLA. These data support a scenario that demonstrates how contemporary events (critical ecological conditions) can moderate historical influences, suggesting that careful interpretation of the evolutionary history of glacial refugia is necessary, especially where populations persisted for a long time but not always with optimal ecological conditions. These peripheral populations are of high conservation value and should be preserved to help conserve the future potential of the species.

Influence of hydrogeographic history and hybridization on the distribution of genetic variation in the pupfishes Cyprinodon atrorus and C. bifasciatus

The evolutionary importance of hybridization in animals has been subject of much debate. In this study, we examined the influence of hydrogeographic history and hybridization on the present distribution of nuclear and mitochondrial DNA variation in two pupfish species, Cyprinodon atrorus and Cyprinodon bifasciatus. Results presented here indicate that there has been limited introgression of nuclear genes; however, mtDNA introgression has been substantial, with complete replacement of the C. bifasciatus mitochondrial genome by that of C. atrorus. Subsequent to this replacement, there has been diversification of mitochondrial haplotypes along major geographic regions in the basin. Evidence was also found that mitochondrial replacement follows a predictable, cyclical pattern in this system, with isolation and diversification followed by re-contact and replacement of C. bifasciatus mitochondrial haplotypes by those of C. atrorus. This pattern is best explained by a combination of a numeric bias towards C. atrorus and mating site selection rather than selection for C. atrorus mitochondrial genome. These results demonstrate the important role hybridization can play in evolution.

Historical biogeography of Mediterranean trout

Complete sequencing of the mitochondrial control region was used to describe phylogenetic relationships of brown trout populations (Salmo trutta) in the Mediterranean river basins of Iberia and to review the historical biogeography of trout from the Mediterranean regions. Phylogenetic relationships among trout lineages suggested that the Danubian one is the most ancestral, in accordance with the eastern origin of most of the European freshwater fish species. Nested-clade and mismatch analyses suggested that the present distribution of haplotypes of the Adriatic and Mediterranean lineages resulted from population expansions originated, respectively, from central and western Europe, which favoured extensive secondary contacts between lineages. Reduced diversity detected within 50% of the analysed populations and large intrabasin differentiation indicated restricted gene flow in post-glacial periods.

Genetic structure of northwestern Spanish brown trout (Salmo trutta L.) populations, differences between microsatellite and allozyme loci

Genetic variation in nine wild brown trout (Salmo trutta L.) populations was studied by means of allozyme and microsatellite markers. All brown trout populations were clearly separated into two clusters that represented the Sil and Duero basins. Although both markers revealed a strong genetic differentiation between basins, microsatellite loci resulted much more accurate when population structure at the intrabasin level was analysed. Also pairwise multilocus FST estimates and assignment tests of individual fish to the set of sampled populations demonstrated a much higher efficiency of microsatellites compared to allozymes. The analysis of both markers provides new insights in defining the conservation units at this local area and confirms the existence of a recognized sub-lineage in the Duero basin. The management implications of these findings are discussed and changes in trout release activity are recommended to avoid mixing of trout gene pools mainly in the Sil basin.

Mitochondrial phylogeography of the Woodmouse (Apodemus sylvaticus) in the Western Palearctic region

We sequenced 965 base pairs of the mitochondrial DNA cytochrome b from 102 woodmice (Apodemus sylvaticus) collected from 40 European localities. The aims of the study were to answer the following questions. (i) Did the Mediterranean peninsulas play a role as refuge for woodmice? (ii) Is genetic variability of A. sylvaticus higher in the Mediterranean region compared with northern Europe? (iii) Are the patterns of the postglacial colonization of Europe by woodmice similar to those presently recognized for other European species? The results provide a clear picture of the impact of the Quaternary glaciations on the genetic and geographical structure of the woodmouse. Our analyses indicate a higher genetic variability of woodmice in the Mediterranean peninsulas compared to northern Europe, suggesting a role of the former as refuge regions for this small mammal. An original pattern of postglacial colonization is proposed where the Iberian and southern France refuge populations colonized almost all European regions. The Sicilian population appears to be very differentiated and highly variable. This emphasizes the importance of this island as a 'hot spot' for the intraspecific genetic diversity of the woodmouse. Finally, woodmice in North Africa originated from southwestern Europe, most probably as a result of a recent anthropogenic introduction.

Phylogeographic congruence between mtDNA and rDNA ITS markers in brown trout

Variation in the internal transcribed spacer (ITS) of rDNA was examined throughout the range of the brown trout (Salmo trutta) to analyze the usefulness of this molecular marker for phylogeographic analysis. The results were compared with those previously obtained with mtDNA, a region exhaustively analyzed along the brown trout distribution. ITS2 was essentially conserved at all populations sampled, no informative characters being detected across the main lineages described in this species. Conversely, ITS1 showed a greater homogenization than other genetic markers at a microgeographic scale, with variation partitioning into several major phylogenetic groups. Phylogeographic patterns were partially congruent between both ITS1 and mtDNA. The main discrepancies were the detection of intra-individual variation and putative recombinant ITS1 sequences in hybridization areas between genetically different, yet historically overlapping, assemblages. Also, the existence of an ancient ITS1 sequence in the Mediterranean-southeastern area (rMEDA), not revealed by mtDNA analysis, was evidenced after rDNA ITS1 analysis.

The role of nuclear genes in intraspecific evolutionary inference: genealogy of the transferrin gene in the brown trout

Technical and biological hurdles have precluded the retrieval of nuclear gene genealogies within most species. Among these obstacles, the possibility of intragenic recombination is one of the most demanding challenges. We studied the utility of nuclear genes for intraspecific evolutionary inferences by selecting a nuclear gene that exhibits patterns of considerable geographic differentiation in the brown trout (Salmo trutta) species complex. Haplotype variation from a nucleotide sequence of approximately 3.7 kb encompassing a portion of the transferrin (TF) gene was surveyed in 31 brown trout individuals collected across the native Eurasian range. Statistically significant recombination and gene conversion events were detected. However, we showed that the substantial cladistic structure was not disrupted by recombination or gene conversion events and the additional structure was estimated to have emerged after those events. Because loci with unusually high levels of variation might indicate the presence of selection, we tested the hypothesis of neutrality and found some evidence for directional selection. The strong geographic signal observed in the TF genealogy, coupled with the current spatial distribution of electromorphs, gave us the ability to draw empirical phylogeographic inferences. We delineated the composition of current brown trout populations on the basis of 3,625 individuals electrophoretically scored for the TF locus. We hypothesized scenarios of historical radiation and dispersal events, thus providing new insights refining previous allozyme and mtDNA inferences. We infer that the most ancestral brown trout populations inhabited tributaries from the Black, Caspian, and Aral Sea drainages. An early radiation of the species occurred throughout the Mediterranean, followed by independent dispersal events from the Adriatic to the southernmost Iberian Atlantic and, more recently, a rapid expansion throughout most of the Atlantic drainages.

Evolution of the mitochondrial control region in Palaearctic brown trout (Salmo trutta) populations: the biogeographical role of the Iberian Peninsula

In order to extend present knowledge of brown trout phylogeography in the Palaearctic, we analysed the complete mitochondrial D-loop sequence (1025-1027 bp) of all mitochondrial haplotypes of Salmo trutta found in the Iberian Peninsula and one North African haplotype. Sequence variation in the mitochondrial control region serves to identify four major haplotype groups within the Iberian Peninsula, i.e. Atlantic, Duero, Mediterranean and Andalusian. Including the Iberian haplotypes, the five main European groups previously established were increased to six: (i) an Atlantic group including two different clusters, South European and North Atlantic; (ii) a group representing an endemism restricted to the Duero basin in the Iberian Peninsula; (iii) an Adriatic-Andalusian group found in two vicariant areas including Adriatic-Ionian populations in the Mediterranean and the Andalusian basins of the southern Iberian Peninsula; (iv) a Mediterranean group with a distribution range that extends from the southwestern basins of the Iberian Peninsula to the Ionian basins of the Greek Peninsula; (v) a Danube group of wide distribution in the Black, Aral and Caspian basins; and (vi) a group comprising the S. t. marmoratus subspecies confined to the Adriatic Sea. The Iberian Peninsula appears to have acted as a physical boundary between haplotypes corresponding to Atlantic- and Mediterranean-draining rivers. Owing to its geographical position, this area has played a major role in present Palaearctic species distribution, as illustrated by its haplotype diversity.

Allozymic evidence of parapatric differentiation of brown trout (Salmo trutta L.) within an Atlantic river basin of the Iberian Peninsula

The genetic variation of brown trout from Duero, one of the main Atlantic Iberian river basins, was assessed at 34 enzymatic loci in 62 native populations. A strong intrabasin differentiation was detected (G(ST) = 0.46; range D: 0-0.066), mainly attributable to the existence of two divergent groups of populations within Duero: southern and northern groups. This divergence was mainly a consequence of the unequal distribution of *75 and *100 alleles at sMDH-B1,2* isoloci, which were correlated with substantial differences in genetic diversity among regions. The Lower Course region (nearly fixed for the *100 allele) and Pisuerga River (nearly fixed for the *75 allele) showed lower heterozygosities (H approximately 0.8%) in contrast with adjacent areas, which evidenced intermediate frequencies for both alleles and higher heterozygosities (H: 2.2-3.1%). Vicariance appeared as the more probable explanation for the significant positive correlation detected between genetic and geographical distances in Duero Basin. Genetic relationships with adjacent Iberian drainages indicate a close similarity between the southern group and Cantabric trout, whereas the northern group constitutes an ancient form from this basin. This study confirmed complex genetic relationships in brown trout from northwest Iberia, reasserting the existence of clines at several loci and for genetic diversity. The interaction between Cantabric and Duero trout, as well as the location of the limit of the anadromous form around the 42 degrees N parallel, are both required to understand the genetic characteristics of brown trout from this area.