Temporal allele frequency changes in large-effect loci reveal potential fishing impacts on salmon life-history diversity.

Fishing has the potential to influence the life-history traits of exploited populations. However, our understanding of how fisheries can induce evolutionary genetic changes remains incomplete. The discovery of large-effect loci linked with ecologically important life-history traits, such as age at maturity in Atlantic salmon (Salmo salar), provides an opportunity to study the impacts of temporally varying fishing pressures on these traits. A 93-year archive of fish scales from wild Atlantic salmon catches from the northern Baltic Sea region allowed us to monitor variation in adaptive genetic diversity linked with age at maturity of wild Atlantic salmon populations. The dataset consisted of samples from both commercial and recreational fisheries that target salmon on their spawning migration. Using a genotyping-by-sequencing approach (GT-seq), we discovered strong within-season allele frequency changes at the vgll3 locus linked with Atlantic salmon age at maturity: fishing in the early season preferentially targeted the vgll3 variant linked with older maturation. We also found within-season temporal variation in catch proportions of different wild Atlantic salmon subpopulations. Therefore, selective pressures of harvesting may vary depending on the seasonal timing of fishing, which has the potential to cause evolutionary changes in key life-history traits and their diversity. This knowledge can be used to guide fisheries management to reduce the effects of fishing practices on salmon life-history diversity. Thus, this study provides a tangible example of using genomic approaches to infer, monitor and help mitigate human impacts on adaptively important genetic variation in nature.

Life-history genotype explains variation in migration activity in Atlantic salmon (Salmo salar).

One of the most well-known life-history continuums is the fast-slow axis, where 'fast' individuals mature earlier than 'slow' individuals. 'Fast' individuals are predicted to be more active than 'slow' individuals because high activity is required to maintain a fast life-history strategy. Recent meta-analyses revealed mixed evidence for such integration. Here, we test whether known life-history genotypes differ in activity expression by using Atlantic salmon (Salmo salar) as a model. In salmon, variation in Vgll3, a transcription cofactor, explains approximately 40% of variation in maturation timing. We predicted that the allele related to early maturation (vgll3*E) would be associated with higher activity. We used an automated surveillance system to follow approximately 1900 juveniles including both migrants and non-migrants (i.e. smolt and parr fish, respectively) in semi-natural conditions over 31 days (approx. 580 000 activity measurements). In migrants, but not in non-migrants, vgll3 explained variation in activity according to our prediction in a sex-dependent manner. Specifically, in females the vgll3*E allele was related to increasing activity, whereas in males the vgll3*L allele (later maturation allele) was related to increasing activity. These sex-dependent effects might be a mechanism maintaining within-population genetic life-history variation.

Population genomics reveals repeated signals of adaptive divergence in the Atlantic salmon of north-eastern Europe.

Our ability to examine genetic variation across entire genomes has enabled many studies searching for the genetic basis of local adaptation. These studies have identified numerous loci as candidates for differential local selection; however, relatively few have examined the overlap among candidate loci identified from independent studies of the same species in different geographic areas or evolutionary lineages. We used an allelotyping approach with a 220K SNP array to characterize the population genetic structure of Atlantic salmon in north-eastern Europe and ask whether the same genomic segments emerged as outliers among populations in different geographic regions. Genome-wide data recapitulated the phylogeographic structure previously inferred from mtDNA and microsatellite markers. Independent analyses of three genetically and geographically distinct groups of populations repeatedly inferred the same 17 haploblocks to contain loci under differential local selection. The most strongly supported of these replicated haploblocks had known strong associations with life-history variation or immune response in Atlantic salmon. Our results are consistent with these genomic segments harbouring large-effect loci which have a major role in Atlantic salmon diversification and are ideal targets for validation studies.

Genomic signatures of fine-scale local selection in Atlantic salmon suggest involvement of sexual maturation, energy homeostasis and immune defence-related genes.

Elucidating the genetic basis of adaptation to the local environment can improve our understanding of how the diversity of life has evolved. In this study, we used a dense SNP array to identify candidate loci potentially underlying fine-scale local adaptation within a large Atlantic salmon (Salmo salar) population. By combining outlier, gene-environment association and haplotype homozygosity analyses, we identified multiple regions of the genome with strong evidence for diversifying selection. Several of these candidate regions had previously been identified in other studies, demonstrating that the same loci could be adaptively important in Atlantic salmon at subdrainage, regional and continental scales. Notably, we identified signals consistent with local selection around genes associated with variation in sexual maturation, energy homeostasis and immune defence. These included the large-effect age-at-maturity gene vgll3, the known obesity gene mc4r, and major histocompatibility complex II. Most strikingly, we confirmed a genomic region on Ssa09 that was extremely differentiated among subpopulations and that is also a candidate for local selection over the global range of Atlantic salmon. This region colocalized with a haplotype strongly associated with spawning ecotype in sockeye salmon (Oncorhynchus nerka), with circumstantial evidence that the same gene (six6) may be the selective target in both cases. The phenotypic effect of this region in Atlantic salmon remains cryptic, although allelic variation is related to upstream catchment area and covaries with timing of the return spawning migration. Our results further inform management of Atlantic salmon and open multiple avenues for future research.

Regulatory Architecture of Gene Expression Variation in the Threespine Stickleback Gasterosteus aculeatus.

Much adaptive evolutionary change is underlain by mutational variation in regions of the genome that regulate gene expression rather than in the coding regions of the genes themselves. An understanding of the role of gene expression variation in facilitating local adaptation will be aided by an understanding of underlying regulatory networks. Here, we characterize the genetic architecture of gene expression variation in the threespine stickleback (Gasterosteus aculeatus), an important model in the study of adaptive evolution. We collected transcriptomic and genomic data from 60 half-sib families using an expression microarray and genotyping-by-sequencing, and located expression quantitative trait loci (eQTL) underlying the variation in gene expression in liver tissue using an interval mapping approach. We identified eQTL for several thousand expression traits. Expression was influenced by polymorphism in both cis- and trans-regulatory regions. Trans-eQTL clustered into hotspots. We did not identify master transcriptional regulators in hotspot locations: rather, the presence of hotspots may be driven by complex interactions between multiple transcription factors. One observed hotspot colocated with a QTL recently found to underlie salinity tolerance in the threespine stickleback. However, most other observed hotspots did not colocate with regions of the genome known to be involved in adaptive divergence between marine and freshwater habitats.

Single nucleotide polymorphisms to discriminate different classes of hybrid between wild Atlantic salmon and aquaculture escapees.

Many wild Atlantic salmon (Salmo salar) populations are threatened by introgressive hybridization from domesticated fish that have escaped from aquaculture facilities. A detailed understanding of the hybridization dynamics between wild salmon and aquaculture escapees requires discrimination of different hybrid classes; however, markers currently available to discriminate the two types of parental genome have limited power to do this. Using a high-density Atlantic salmon single nucleotide polymorphism (SNP) array, in combination with pooled-sample allelotyping and an Fst outlier approach, we identified 200 SNPs that differentiated an important Atlantic salmon stock from the escapees potentially hybridizing with it. By simulating multiple generations of wild-escapee hybridization, involving wild populations in two major phylogeographic lineages and a genetically diverse set of escapees, we showed that both the complete set of SNPs and smaller subsets could reliably assign individuals to different hybrid classes up to the third hybrid (F3) generation. This set of markers will be a useful tool for investigating the genetic interactions between native wild fish and aquaculture escapees in many Atlantic salmon populations.

Integrative Taxonomy Recognizes Evolutionary Units Despite Widespread Mitonuclear Discordance: Evidence from a Rotifer Cryptic Species Complex.

Mitonuclear discordance across taxa is increasingly recognized as posing a major challenge to species delimitation based on DNA sequence data. Integrative taxonomy has been proposed as a promising framework to help address this problem. However, we still lack compelling empirical evidence scrutinizing the efficacy of integrative taxonomy in relation to, for instance, complex introgression scenarios involving many species. Here, we report remarkably widespread mitonuclear discordance between about 15 mitochondrial and 4 nuclear Brachionus calyciflorus groups identified using different species delimitation approaches. Using coalescent-, Bayesian admixture-, and allele sharing-based methods with DNA sequence or microsatellite data, we provide strong evidence in support of hybridization as a driver of the observed discordance. We then describe our combined molecular, morphological, and ecological approaches to resolving phylogenetic conflict and inferring species boundaries. Species delimitations based on the ITS1 and 28S nuclear DNA markers proved a more reliable predictor of morphological variation than delimitations using the mitochondrial COI gene. A short-term competition experiment further revealed systematic differences in the competitive ability between two of the nuclear-delimited species under six different growth conditions, independent of COI delimitations; hybrids were also observed. In light of these findings, we discuss the failure of the COI marker to estimate morphological stasis and morphological plasticity in the B. calyciflorus complex. By using B. calyciflorus as a representative case, we demonstrate the potential of integrative taxonomy to guide species delimitation in the presence of mitonuclear phylogenetic conflicts.

Genome-wide SNP analysis reveals a genetic basis for sea-age variation in a wild population of Atlantic salmon (Salmo salar).

Delaying sexual maturation can lead to larger body size and higher reproductive success, but carries an increased risk of death before reproducing. Classical life history theory predicts that trade-offs between reproductive success and survival should lead to the evolution of an optimal strategy in a given population. However, variation in mating strategies generally persists, and in general, there remains a poor understanding of genetic and physiological mechanisms underlying this variation. One extreme case of this is in the Atlantic salmon (Salmo salar), which can show variation in the age at which they return from their marine migration to spawn (i.e. their 'sea age'). This results in large size differences between strategies, with direct implications for individual fitness. Here, we used an Illumina Infinium SNP array to identify regions of the genome associated with variation in sea age in a large population of Atlantic salmon in Northern Europe, implementing individual-based genome-wide association studies (GWAS) and population-based FST outlier analyses. We identified several regions of the genome which vary in association with phenotype and/or selection between sea ages, with nearby genes having functions related to muscle development, metabolism, immune response and mate choice. In addition, we found that individuals of different sea ages belong to different, yet sympatric populations in this system, indicating that reproductive isolation may be driven by divergence between stable strategies. Overall, this study demonstrates how genome-wide methodologies can be integrated with samples collected from wild, structured populations to understand their ecology and evolution in a natural context.

Chromosome rearrangements, recombination suppression, and limited segregation distortion in hybrids between Yellowstone cutthroat trout (Oncorhynchus clarkii bouvieri) and rainbow trout (O. mykiss).

BACKGROUND: Introgressive hybridization is an important evolutionary process that can lead to the creation of novel genome structures and thus potentially new genetic variation for selection to act upon. On the other hand, hybridization with introduced species can threaten native species, such as cutthroat trout (Oncorhynchus clarkii) following the introduction of rainbow trout (O. mykiss). Neither the evolutionary consequences nor conservation implications of rainbow trout introgression in cutthroat trout is well understood. Therefore, we generated a genetic linkage map for rainbow-Yellowstone cutthroat trout (O. clarkii bouvieri) hybrids to evaluate genome processes that may help explain how introgression affects hybrid genome evolution. RESULTS: The hybrid map closely aligned with the rainbow trout map (a cutthroat trout map does not exist), sharing all but one linkage group. This linkage group (RYHyb20) represented a fusion between an acrocentric (Omy28) and a metacentric chromosome (Omy20) in rainbow trout. Additional mapping in Yellowstone cutthroat trout indicated the two rainbow trout homologues were fused in the Yellowstone genome. Variation in the number of hybrid linkage groups (28 or 29) likely depended on a Robertsonian rearrangement polymorphism within the rainbow trout stock. Comparison between the female-merged F1 map and a female consensus rainbow trout map revealed that introgression suppressed recombination across large genomic regions in 5 hybrid linkage groups. Two of these linkage groups (RYHyb20 and RYHyb25_29) contained confirmed chromosome rearrangements between rainbow and Yellowstone cutthroat trout indicating that rearrangements may suppress recombination. The frequency of allelic and genotypic segregation distortion varied among parents and families, suggesting few incompatibilities exist between rainbow and Yellowstone cutthroat trout genomes. CONCLUSIONS: Chromosome rearrangements suppressed recombination in the hybrids. This result supports several previous findings demonstrating that recombination suppression restricts gene flow between chromosomes that differ by arrangement. Conservation of synteny and map order between the hybrid and rainbow trout maps and minimal segregation distortion in the hybrids suggest rainbow and Yellowstone cutthroat trout genomes freely introgress across chromosomes with similar arrangement. Taken together, these results suggest that rearrangements impede introgression. Recombination suppression across rearrangements could enable large portions of non-recombined chromosomes to persist within admixed populations.

The genomic trajectory of hybrid swarms: outcomes of repeated crosses between populations of Tigriopus californicus.

Introgressive hybridization between genetically divergent populations is an important evolutionary process. The degree to which repeated hybridization events between the same parental taxa lead to similar genomic outcomes is unknown. This study addressed this question by following genomic trajectories of replicate hybrid swarms of the copepod Tigriopus californicus over many generations of free mating. Swarm composition was determined both by differential reproductive success of founder individuals and subsequent selection on hybrid genotypes. For one cross, between two populations showing differential fitness in the laboratory and no hybrid breakdown, the genetic trajectory was highly repeatable: replicates rapidly became dominated by alleles from the fitter parent. In a second cross, between two populations showing similar fitness and significant F2 hybrid breakdown, alleles from alternative populations dominated different replicates. Swarms exhibited a general temporal trend of decreasing cytonuclear mismatch. Some patterns of differential introgression across the genome were strikingly congruent amongst swarm replicates, both within and between cross types, and reflected patterns of segregation distortion previously observed within controlled crosses between the same parental populations. Differences in heterozygosity between the sexes, and evidence for a previously suspected sex-distortion locus, suggest that complex interactions between sex and genotype influence hybrid swarm outcome.

Discovery and characterization of novel genetic markers for use in the management of Lahontan cutthroat trout (Oncorhynchus clarkii henshawi).

The Lahontan cutthroat trout (Oncorhynchus clarkii henshawi) is threatened by habitat destruction, over-harvest and hybridization with nonnative trout. Currently, three Geographic Management Units (GMUs) are recognized within the taxon. Here, we describe a suite of 68 single-nucleotide polymorphism (SNP) genetic markers for use in the study and management of Lahontan cutthroat trout and a closely related subspecies, the Paiute cutthroat trout (O. c. seleneris). These include markers variable within the two subspecies (n = 35), diagnostic for the two subspecies (n = 23) and diagnostic for Yellowstone cutthroat trout (O. c. bouvieri) and other closely related subspecies (n = 10). Sixty-three markers were discovered by Sanger sequencing of 171 EST loci in an ascertainment panel including Lahontan cutthroat trout from four populations representing all GMUs. Five markers were identified in a secondary sequencing effort with a single population of Lahontan cutthroat trout. TaqMan assays were validated on six Lahontan cutthroat trout populations and a diverse panel of other trout. Over 90% of the markers variable in Lahontan cutthroat trout were polymorphic in at least two populations, and 66% were variable within all three GMUs. All Lahontan diagnostic markers were also fixed for the Lahontan allele in Paiute cutthroat trout. Most of the Yellowstone diagnostic markers can also be used for this purpose in other cutthroat trout subspecies. This is the first set of SNP markers to be developed for Lahontan cutthroat trout, and will be an important tool for conservation and management.

Interpopulation hybridization results in widespread viability selection across the genome in Tigriopus californicus.

BACKGROUND: Genetic interactions within hybrids influence their overall fitness. Understanding the details of these interactions can improve our understanding of speciation. One experimental approach is to investigate deviations from Mendelian expectations (segregation distortion) in the inheritance of mapped genetic markers. In this study, we used the copepod Tigriopus californicus, a species which exhibits high genetic divergence between populations and a general pattern of reduced fitness in F2 interpopulation hybrids. Previous studies have implicated both nuclear-cytoplasmic and nuclear-nuclear interactions in causing this fitness reduction. We identified and mapped population-diagnostic single nucleotide polymorphisms (SNPs) and used these to examine segregation distortion across the genome within F2 hybrids. RESULTS: We generated a linkage map which included 45 newly elucidated SNPs and 8 population-diagnostic microsatellites used in previous studies. The map, the first available for the Copepoda, was estimated to cover 75% of the genome and included markers on all 12 T. californicus chromosomes. We observed little segregation distortion in newly hatched F2 hybrid larvae (fewer than 10% of markers at p < 0.05), but strikingly higher distortion in F2 hybrid adult males (45% of markers at p < 0.05). Hence, segregation distortion was primarily caused by selection against particular genetic combinations which acted between hatching and maturity. Distorted markers were not distributed randomly across the genome but clustered on particular chromosomes. In contrast to other studies in this species we found little evidence for cytonuclear coadaptation. Instead, different linkage groups exhibited markedly different patterns of distortion, which appear to have been influenced by nuclear-nuclear epistatic interactions and may also reflect genetic load carried within the parental lines. CONCLUSION: Adult male F2 hybrids between two populations of T. californius exhibit dramatic segregation distortion across the genome. Distorted loci are clustered within specific linkage groups, and the direction of distortion differs between chromosomes. This segregation distortion is due to selection acting between hatching and adulthood.

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.

Inter and intra-population variation in shoaling and boldness in the zebrafish (Danio rerio).

Population differences in anti-predator behaviour have been demonstrated in several species, although less is known about the genetic basis of these traits. To determine the extent of genetic differences in boldness (defined as exploration of a novel object) and shoaling within and between zebrafish (Danio rerio) populations, and to examine the genetic basis of shoaling behaviour in general, we carried out a study that involved laboratory-raised fish derived from four wild-caught populations. Controlling for differences in rearing environment, significant inter-population differences were found in boldness but not shoaling. A larger shoaling experiment was also performed using one of the populations as the basis of a North Carolina type II breeding design (174 fish in total) to estimate heritability of shoaling tendency. A narrow-sense heritability estimate of 0.40 was obtained, with no apparent dominance effects.