Rapid evolution in salmon life history induced by direct and indirect effects of fishing.

Understanding the drivers of evolution is a fundamental aim in biology. However, identifying the evolutionary impacts of human activities is challenging because of a lack of temporal data and limited knowledge of the genetic basis of most traits. Here, we identify the drivers of evolution toward maturity at an earlier age in Atlantic salmon through two types of fisheries-induced evolution acting in opposing directions: an indirect effect linked with harvest of a salmon prey species (capelin) at sea (selection against late maturation) and a direct effect due to net fishing in rivers (selection against early maturation). Because capelin are harvested as an aquaculture feed protein source, we hereby determine an indirect path by which salmon aquaculture may influence wild salmon populations.

Diversity and linkage disequilibrium in farmed Tasmanian Atlantic salmon.

Farmed Atlantic salmon (Salmo salar) is a globally important production species, including in Australia where breeding and selection has been in progress since the 1960s. The recent development of SNP genotyping platforms means genome-wide association and genomic prediction can now be implemented to speed genetic gain. As a precursor, this study collected genotypes at 218 132 SNPs in 777 fish from a Tasmanian breeding population to assess levels of genetic diversity, the strength of linkage disequilibrium (LD) and imputation accuracy. Genetic diversity in Tasmanian Atlantic salmon was lower than observed within European populations when compared using four diversity metrics. The distribution of allele frequencies also showed a clear difference, with the Tasmanian animals carrying an excess of low minor allele frequency variants. The strength of observed LD was high at short distances (<25 kb) and remained above background for marker pairs separated by large chromosomal distances (hundreds of kb), in sharp contrast to the European Atlantic salmon tested. Genotypes were used to evaluate the accuracy of imputation from low density (0.5 to 5 K) up to increased density SNP sets (78 K). This revealed high imputation accuracies (0.89-0.97), suggesting that the use of low density SNP sets will be a successful approach for genomic prediction in this population. The long-range LD, comparatively low genetic diversity and high imputation accuracy in Tasmanian salmon is consistent with known aspects of their population history, which involved a small founding population and an absence of subsequent introgression. The findings of this study represent an important first step towards the design of methods to apply genomics in this economically important population.

From population genomics to conservation and management: a workflow for targeted analysis of markers identified using genome-wide approaches in Atlantic salmon Salmo salar.

A genotyping assay for the Ion Torrent Ion PGM platform was developed for fast and cost-effective targeted genotyping of key single nucleotide polymorphisms (SNPs) earlier identified using a genome-wide SNP array in Atlantic salmon Salmo salar. The method comprised a simple primer design step for multiplex-polymerase chain reaction (PCR), followed by two rounds of Ion Torrent Ion PGM sequencing to empirically evaluate marker efficiency in large multiplexes and to optimise or exclude them when necessary. Of 282 primer pairs initially tested, 217 were successfully amplified, indicating good amplification success (>75%). These markers included the sdy partial gene product to determine genetic sex, as well as three additional modules comprising SNPs for assessing neutral genetic variation (NSNP = 150), examining functional genetic variation associated with sea age at maturity (NSNP = 5), and for performing genetic subpopulation assignment (NSNP = 61). The assay was primarily developed to monitor long-term genetic changes in S. salar from the Teno River, but modules are likely suitable for application in a wide range of S. salar populations. Furthermore, the fast and versatile assay development pipeline offers a strategy for developing targeted sequencing assays in any species.

Generation of a neutral FST baseline for testing local adaptation on gill raker number within and between European whitefish ecotypes in the Baltic Sea basin.

Divergent selection at ecologically important traits is thought to be a major factor driving phenotypic differentiation between populations. To elucidate the role of different evolutionary processes shaping the variation in gill raker number of European whitefish (Coregonus lavaretus sensu lato) in the Baltic Sea basin, we assessed the relationships between genetic and phenotypic variation among and within three whitefish ecotypes (sea spawners, river spawners and lake spawners). To generate expected neutral distribution of FST and to evaluate whether highly variable microsatellite loci resulted in deflated FST estimates compared to less variable markers, we performed population genetic simulations under finite island and hierarchical island models. The genetic divergence observed among (FCT = 0.010) and within (FST = 0.014-0.041) ecotypes was rather low. The divergence in gill raker number, however, was substantially higher between sea and river spawners compared to observed microsatellite data and simulated neutral baseline (PCT > FCT ). This suggests that the differences in gill raker number between sea and river spawners are likely driven by divergent natural selection. We also found strong support for divergent selection on gill raker number among different populations of sea spawners (PST > FST ), most likely caused by highly variable habitat use and diverse diet. The putative role of divergent selection within lake spawners initially inferred from empirical microsatellite data was not supported by simulated FST distributions. This work provides a first formal test of divergent selection on gill raker number in Baltic whitefish, and demonstrates the usefulness of population genetic simulations to generate informative neutral baselines for PST -FST analyses helping to disentangle the effects of stochastic evolutionary processes from natural selection.

Sympatric divergence and clinal variation in multiple coloration traits of Ficedula flycatchers.

Geographic variation in phenotypes plays a key role in fundamental evolutionary processes such as local adaptation, population differentiation and speciation, but the selective forces behind it are rarely known. We found support for the hypothesis that geographic variation in plumage traits of the pied flycatcher Ficedula hypoleuca is explained by character displacement with the collared flycatcher Ficedula albicollis in the contact zone. The plumage traits of the pied flycatcher differed strongly from the more conspicuous collared flycatcher in a sympatric area but increased in conspicuousness with increasing distance to there. Phenotypic differentiation (PST ) was higher than that in neutral genetic markers (FST ), and the effect of geographic distance remained when statistically controlling for neutral genetic differentiation. This suggests that a cline created by character displacement and gene flow explains phenotypic variation across the distribution of this species. The different plumage traits of the pied flycatcher are strongly to moderately correlated, indicating that they evolve non-independently from each other. The flycatchers provide an example of plumage patterns diverging in two species that differ in several aspects of appearance. The divergence in sympatry and convergence in allopatry in these birds provide a possibility to study the evolutionary mechanisms behind the highly divergent avian plumage patterns.

Genetic mixed-stock analysis of lake-run brown trout Salmo trutta fishery catches in the Inari Basin, northern Finland: implications for conservation and management.

Genetic mixed-stock analysis (MSA) of wild lake-run brown trout Salmo trutta fishery catches (n = 665) from the Inari Basin (northern Finland) between 2006 and 2008 was carried out using a previously characterized baseline with 30 populations (n = 813) and 13 microsatellite loci. Altogether, 12 populations contributed significantly to mixed-stock fisheries, with the Ivalojoki system being the major contributor (70%) to the total catch. When catches were analysed regionally, geographically nearby populations were the main contributors to the local catches, indicating that a large proportion of S. trutta occupy lacustrine areas near the natal river mouth rather than dispersing throughout the lake. Similarly, far upstream populations contributed insignificantly to catches. These findings have important implications for the conservation and sustainable fishery management of the Inari system.

Annotated genes and nonannotated genomes: cross-species use of Gene Ontology in ecology and evolution research.

Recent advances in molecular technologies have opened up unprecedented opportunities for molecular ecologists to better understand the molecular basis of traits of ecological and evolutionary importance in almost any organism. Nevertheless, reliable and systematic inference of functionally relevant information from these masses of data remains challenging. The aim of this review is to highlight how the Gene Ontology (GO) database can be of use in resolving this challenge. The GO provides a largely species-neutral source of information on the molecular function, biological role and cellular location of tens of thousands of gene products. As it is designed to be species-neutral, the GO is well suited for cross-species use, meaning that, functional annotation derived from model organisms can be transferred to inferred orthologues in newly sequenced species. In other words, the GO can provide gene annotation information for species with nonannotated genomes. In this review, we describe the GO database, how functional information is linked with genes/gene products in model organisms, and how molecular ecologists can utilize this information to annotate their own data. Then, we outline various applications of GO for enhancing the understanding of molecular basis of traits in ecologically relevant species. We also highlight potential pitfalls, provide step-by-step recommendations for conducting a sound study in nonmodel organisms, suggest avenues for future research and outline a strategy for maximizing the benefits of a more ecological and evolutionary genomics-oriented ontology by ensuring its compatibility with the GO.

Transcription and redox enzyme activities: comparison of equilibrium and disequilibrium levels in the three-spined stickleback.

Evolutionary and acclimatory responses require functional variability, but in contrast with mRNA and protein abundance data, most physiological measurements cannot be obtained in a high-throughput manner. Consequently, one must either rely on high-throughput transcriptomic or proteomic data with only predicted functional information, or accept the limitation that most physiological measurements can give fewer data than those provided by transcriptomics or proteomics. We evaluated how transcriptional and redox enzyme activity data agreed with regard to population differentiation (i.e. a system in steady state in which any time lag between transcription, translation and post-translational effects would be irrelevant) and in response to an acute 6°C increase in temperature (i.e. a disequilibrium state wherein translation could not have caught up with transcription) in the three-spined stickleback (Gasterosteus aculeatus). Transcriptional and enzyme activity data corresponded well with regard to population differentiation, but less so with regard to acute temperature increase. The data thus suggest that transcriptional and functional measurements can lead to similar conclusions when a biological system is in a steady state. The responses to acute changes must, as has been demonstrated earlier, be based on changes in cellular conditions or properties of existing proteins without significant de novo synthesis of new gene products.

Phylogenetic status of brown trout Salmo trutta populations in five rivers from the southern Caspian Sea and two inland lake basins, Iran: a morphogenetic approach.

Interrelationships, origin and phylogenetic affinities of brown trout Salmo trutta populations from the southern Caspian Sea basin, Orumieh and Namak Lake basins in Iran were analysed from complete mtDNA control region sequences, 12 microsatellite loci and morphological characters. Among 129 specimens from six populations, seven haplotypes were observed. Based on mtDNA haplotype data, the Orumieh and southern Caspian populations did not differ significantly, but the Namak basin-Karaj population presented a unique haplotype closely related to the haplotypes of the other populations (0·1% Kimura two-parameter, K2P divergence). All Iranian haplotypes clustered as a distinct group within the Danube phylogenetic grouping, with an average K2P distance of 0·41% relative to other Danubian haplotypes. The Karaj haplotype in the Namak basin was related to a haplotype (Da26) formerly identified in the Tigris basin in Turkey, to a Salmo trutta oxianus haplotype from the Aral Sea basin, and to haplotype Da1a with two mutational steps, as well as to other Iranian haplotypes with one to two mutational steps, which may indicate a centre of origin in the Caspian basin. In contrast to results of the mtDNA analysis, more pronounced differentiation was observed among the populations studied in the morphological and microsatellite DNA data, except for the two populations from the Orumieh basin, which were similar, possibly due to anthropogenic causes.

Genetic differentiation between two sympatric morphs of the blind Iran cave barb Iranocypris typhlops.

The phylogenetic relationship between two sympatric morphotypes of the Iran cave barb Iranocypris typhlops, and Garra rufa, was investigated by sequencing the cytochrome c oxidase I (coI) region (788 bp) providing the first molecular evidence of their phylogeny. Consistent with their morphological differences, the mean genetic distance between the two forms of I. typhlops was significantly higher than generally reported for intraspecific divergence in freshwater fishes. They were phylogenetically closer to G. rufa than to any other species.

Candidate genes for colour and vision exhibit signals of selection across the pied flycatcher (Ficedula hypoleuca) breeding range.

The role of natural selection in shaping adaptive trait differentiation in natural populations has long been recognized. Determining its molecular basis, however, remains a challenge. Here, we search for signals of selection in candidate genes for colour and its perception in a passerine bird. Pied flycatcher plumage varies geographically in both its structural and pigment-based properties. Both characteristics appear to be shaped by selection. A single-locus outlier test revealed 2 of 14 loci to show significantly elevated signals of divergence. The first of these, the follistatin gene, is expressed in the developing feather bud and is found in pathways with genes that determine the structure of feathers and may thus be important in generating variation in structural colouration. The second is a gene potentially underlying the ability to detect this variation: SWS1 opsin. These two loci were most differentiated in two Spanish pied flycatcher populations, which are also among the populations that have the highest UV reflectance. The follistatin and SWS1 opsin genes thus provide strong candidates for future investigations on the molecular basis of adaptively significant traits and their co-evolution.

Microsatellite standardization and evaluation of genotyping error in a large multi-partner research programme for conservation of Atlantic salmon (Salmo salar L.).

Microsatellite genotyping is a common DNA characterization technique in population, ecological and evolutionary genetics research. Since different alleles are sized relative to internal size-standards, different laboratories must calibrate and standardize allelic designations when exchanging data. This interchange of microsatellite data can often prove problematic. Here, 16 microsatellite loci were calibrated and standardized for the Atlantic salmon, Salmo salar, across 12 laboratories. Although inconsistencies were observed, particularly due to differences between migration of DNA fragments and actual allelic size ('size shifts'), inter-laboratory calibration was successful. Standardization also allowed an assessment of the degree and partitioning of genotyping error. Notably, the global allelic error rate was reduced from 0.05 ± 0.01 prior to calibration to 0.01 ± 0.002 post-calibration. Most errors were found to occur during analysis (i.e. when size-calling alleles; the mean proportion of all errors that were analytical errors across loci was 0.58 after calibration). No evidence was found of an association between the degree of error and allelic size range of a locus, number of alleles, nor repeat type, nor was there evidence that genotyping errors were more prevalent when a laboratory analyzed samples outside of the usual geographic area they encounter. The microsatellite calibration between laboratories presented here will be especially important for genetic assignment of marine-caught Atlantic salmon, enabling analysis of marine mortality, a major factor in the observed declines of this highly valued species.

Strong gene flow and lack of stable population structure in the face of rapid adaptation to local temperature in a spring-spawning salmonid, the European grayling (Thymallus thymallus).

Gene flow has the potential to both constrain and facilitate adaptation to local environmental conditions. The early stages of population divergence can be unstable because of fluctuating levels of gene flow. Investigating temporal variation in gene flow during the initial stages of population divergence can therefore provide insights to the role of gene flow in adaptive evolution. Since the recent colonization of Lake Lesjaskogsvatnet in Norway by European grayling (Thymallus thymallus), local populations have been established in over 20 tributaries. Multiple founder events appear to have resulted in reduced neutral variation. Nevertheless, there is evidence for local adaptation in early life-history traits to different temperature regimes. In this study, microsatellite data from almost a decade of sampling were assessed to infer population structuring and its temporal stability. Several alternative analyses indicated that spatial variation explained 2-3 times more of the divergence in the system than temporal variation. Over all samples and years, there was a significant correlation between genetic and geographic distance. However, decomposed pairwise regression analysis revealed differing patterns of genetic structure among local populations and indicated that migration outweighs genetic drift in the majority of populations. In addition, isolation by distance was observable in only three of the six years, and signals of population bottlenecks were observed in the majority of samples. Combined, the results suggest that habitat-specific adaptation in this system has preceded the development of consistent population substructuring in the face of high levels of gene flow from divergent environments.

High level of population genetic structuring in lake-run brown trout, Salmo trutta, of the Inari Basin, northern Finland.

Rivers draining into (Lake) Inarijärvi, northern Finland, sustain a number of lake-run brown trout, Salmo trutta, populations but, as with most lake-run S. trutta systems, the level of population genetic structuring among populations is unknown. To address this and to assist fish stock management in the region, the population genetic structure of S. trutta collected from 28 sampling sites in rivers flowing into Inarijärvi was studied using 13 microsatellite loci. Populations were clustered into three separate groups, largely corresponding to geographic regions, with between-region F(ST) values ranging from 0·11 to 0·16. The significant differentiation observed between most populations within each region also implies that individual populations should be recognized as separate management units and actions to improve, and subsequently maintain, conditions for natural spawning should be prioritized. The results of this study further indicate that the trout from each of these regions may have different biological characteristics, such as local-lake feeding behaviour among the western populations and strong isolation among the northern stocks. As a consequence, further research is warranted to better understand the level of ecological uniqueness of lake-run S. trutta populations.

Beyond MHC: signals of elevated selection pressure on Atlantic salmon (Salmo salar) immune-relevant loci.

Using Atlantic salmon (Salmo salar) as a model system, we investigated whether 18 microsatellites tightly linked to immune-relevant genes have experienced different selection pressures than 76 loci with no obvious association with immune function. Immune-relevant loci were identified as outliers by two outlier tests significantly more often than nonimmune linked loci (22% vs. 1.6%). In addition, the allele frequencies of immune relevant markers were more often correlated with latitude and temperature. Combined, these results support the hypothesis that immune-relevant loci more frequently exhibit footprints of selection than other loci. They also indicate that the correlation between immune-relevant loci and latitude may be due to temperature-induced differences in pathogen-driven selection or some other environmental factor correlated with latitude.

High Gyrodactylus salaris infection rate in triploid Atlantic salmon Salmo salar.

We describe an unusually high infection rate of Gyrodactylus salaris Malmberg in juvenile Atlantic salmon Salmo salar L. of Baltic Sea origin, which are generally believed to be more resistant to G. salaris than East Atlantic salmon populations. Based on analyses of mitochondrial (complete cytochrome oxidase 1 [CO1] gene, 1548 bp) and nuclear (ADNAM1, 435 bp; internal transcribed spacer [ITS] rDNA region, 1232 bp) DNA fragments, the closest relatives of the characterized Estonian G. salaris strain were parasites found off the Swedish west coast and in Raasakka hatchery, Iijoki (Baltic Sea, Finland). Analyses of 14 microsatellite loci of the host S. salarrevealed that approximately 40% of studied fish were triploids. We subsequently identified triploid Atlantic salmon of Baltic origin as more susceptible to G. salaris infection than their diploid counterparts, possibly due to compromised complement-dependent immune pathways in triploid salmon. This is in accordance with earlier studies that have shown elevated susceptibility of triploids to various viral or bacterial pathogens, and represents one of the first reports of increased susceptibility of triploid salmonid fish to an ectoparasite. However, further experimental work is needed to determine whether triploid Atlantic salmon is generally more susceptible to G. salaris compared to their diploid counterparts, irrespective of the particular triploidization method and population of origin.

Seventy new microsatellites for the pied flycatcher, Ficedula hypoleuca and amplification in other passerine birds.

The pied flycatcher (Ficedula hypoleuca) is a small migratory passerine bird commonly distributed across Europe which has been the focus of considerable ecological and evolutionary research. Here, we present details of 70 microsatellite markers for the species adding to the six which are currently available. Sixty-six markers were also polymorphic in the closely related collared flycatcher (Ficedula albicollis), while 54 were polymorphic in another related passerine, the bluethroat (Luscinia svecica), and 12 were polymorphic in the more distantly related Siberian jay (Perisoreus infaustus).

Use of differential expression data for identification of novel immune relevant expressed sequence tag-linked microsatellite markers in Atlantic salmon (Salmo salar L.).

Despite the large number of genes contributing to the immune response, wildlife immunogenetic studies have tended to focus mostly on the major histocompatibility complex-related genes. Here, we utilized previously published microarray and competitive RNA hybridization information to identify 3750 immune relevant Atlantic salmon (Salmo salar) expressed sequence tags. We then identified those expressed sequence tags containing microsatellites and subsequently designed 48 primer pairs and tested them for polymorphism in Atlantic salmon. Altogether, 16 polymorphic markers were characterized, with allele numbers ranging from two to 18, and these 16 loci were further tested in five other salmonid species.