Candidate loci reveal genetic differentiation between temporally divergent migratory runs of Chinook salmon (Oncorhynchus tshawytscha)

Rhythmic Clock Gene Expression in Atlantic Salmon Parr Brain

To better understand the complexity of clock genes in salmonids, a taxon with an additional whole genome duplication, an analysis was performed to identify and classify gene family members (clock, arntl, period, cryptochrome, nr1d, ror, and csnk1). The majority of clock genes, in zebrafish and Northern pike, appeared to be duplicated. In comparison to the 29 clock genes described in zebrafish, 48 clock genes were discovered in salmonid species. There was also evidence of species-specific reciprocal gene losses conserved to the Oncorhynchus sister clade. From the six period genes identified three were highly significantly rhythmic, and circadian in their expression patterns (per1a.1, per1a.2, per1b) and two was significantly rhythmically expressed (per2a, per2b). The transcriptomic study of juvenile Atlantic salmon (parr) brain tissues confirmed gene identification and revealed that there were 2,864 rhythmically expressed genes (p < 0.001), including 1,215 genes with a circadian expression pattern, of which 11 were clock genes. The majority of circadian expressed genes peaked 2 h before and after daylight. These findings provide a foundation for further research into the function of clock genes circadian rhythmicity and the role of an enriched number of clock genes relating to seasonal driven life history in salmonids.

Spatial and environmental influences on selection in a clock gene coding trinucleotide repeat in Canada lynx (Lynx canadensis)

Clock genes exhibit substantial control over gene expression and ultimately life-histories using external cues such as photoperiod, and are thus likely to be critical for adaptation to shifting seasonal conditions and novel environments as species redistribute their ranges under climate change. Coding trinucleotide repeats (cTNRs) are found within several clock genes, and may be interesting targets of selection due to their containment within exonic regions and elevated mutation rates. Here, we conduct inter-specific characterization of the NR1D1 cTNR between Canada lynx and bobcat, and intra-specific spatial and environmental association analyses of neutral microsatellites and our functional cTNR marker, to investigate the role of selection on this locus in Canada lynx. We report signatures of divergent selection between lynx and bobcat, with the potential for hybrid-mediated gene flow in the area of range overlap. We also provide evidence that this locus is under selection across Canada lynx in eastern Canada, with both spatial and environmental variables significantly contributing to the explained variation, after controlling for neutral population structure. These results suggest that cTNRs may play an important role in the generation of functional diversity within some mammal species, and allow for contemporary rates of adaptation in wild populations in response to environmental change. We encourage continued investment into the study of cTNR markers to better understand their broader relevance to the evolution and adaptation of mammals.

Estimating the contribution of Greenland Halibut (Reinhardtius hippoglossoides) stocks to nurseries by means of genotyping-by-sequencing: Sex and time matter

Identification of stocks and quantification of their relative contribution to recruitment are major objectives toward improving the management and conservation of marine exploited species. Next-generation sequencing allows for thousands of genomic markers to be analyzed, which provides the resolution needed to address these questions in marine species with weakly differentiated populations. Greenland Halibut (Reinhardtius hippoglossoides) is one of the most important exploited demersal species throughout the North Atlantic, in particular in the Gulf of St. Lawrence, Canada. There, two nurseries are known, the St. Lawrence Estuary and the northern Anticosti Island, but their contribution to the renewal of stocks remains unknown. The goals of this study were (a) to document the genetic structure and (b) to estimate the contribution of the different identified breeding stocks to nurseries. We sampled 100 juveniles per nursery and 50 adults from seven sites ranging from Saguenay Fjord to offshore Newfoundland, with some sites sampled over two consecutive years in order to evaluate the temporal stability of the contribution. Our results show that after removing sex-linked markers, the Estuary/Gulf of St. Lawrence represents a single stock which is genetically distinct from the Atlantic around Newfoundland (F ST = 0.00146, p-value = .001). Population assignment showed that recruitment in both nurseries is predominantly associated with the St. Lawrence stock. However, we found that the relative contribution of both stocks to the nurseries is temporally variable with 1% contribution of the Newfoundland stock one year but up to 33% for the second year, which may be caused by year-to-year variation in larval transport into the Gulf of St. Lawrence. This study serves as a model for the identification of stocks for fisheries resources in a context where few barriers to dispersal occur, in addition to demonstrating the importance of considering sex-linked markers and temporal replicates in studies of population genomics.

Temporal dynamics of migration-linked genetic variation are driven by streamflows and riverscape permeability

Landscape permeability is often explored spatially, but may also vary temporally. Landscape permeability, including partial barriers, influences migratory animals that move across the landscape. Partial barriers are common in rivers where barrier passage varies with streamflow. We explore the influence of partial barriers on the spatial and temporal distribution of migration‐linked genotypes of Oncorhynchus mykiss, a salmonid fish with co‐occurring resident and migratory forms, in tributaries to the South Fork Eel River, California, USA, Elder and Fox Creeks. We genotyped >4,000 individuals using RAD‐capture and classified individuals as resident, heterozygous or migratory genotypes using life history‐associated loci. Across four years of study (2014–2017), the permeability of partial barriers varied across dry and wet years. In Elder Creek, the largest waterfall was passable for adults migrating up‐river 4–39 days each year. In this stream, the overall spatial pattern, with fewer migratory genotypes above the waterfall, remained true across dry and wet years (67%–76% of migratory alleles were downstream of the waterfall). We also observed a strong relationship between distance upstream and proportion of migratory alleles. In Fox Creek, the primary barrier is at the mouth, and we found that the migratory allele frequency varied with the annual timing of high flow events. In years when rain events occurred during the peak breeding season, migratory allele frequency was high (60%–68%), but otherwise it was low (30% in two years). We highlight that partial barriers and landscape permeability can be temporally dynamic, and this effect can be observed through changing genotype frequencies in migratory animals.

Testing the expression of circadian clock genes in the tissues of Chinook salmon, Oncorhynchus tshawytscha

Animals have an endogenous circadian clock that temporally regulates 24 hour (h) oscillations in behavior and physiology. This highly conserved mechanism consists of two positive regulators, Bmal and Clock, and two negative regulators, Cry and Per, that run with a 24-h cycle that synchronizes itself with environmental changes in light, food, and temperature. We examined the circadian clock in Chinook salmon (Oncorhynchus tshawytscha), a non-model organism in which the function of the clock has not been studied. Recent studies indicate that clock genes in Chinook salmon play a role in its evolution of local adaptation, possibly by influencing migration timing. We designed real-time quantitative PCR (RT-qPCR) assays to quantify the transcription of components of the clock system, and validated these for PCR efficiency and specificity in detecting Chinook target genes. Chinook salmon tissue samples were collected in 3-h intervals, over the course of 24 h, from five different organs. Our data indicate that the circadian clock functions differently in each of these tissues. In the liver, positive and negative regulators exhibit anti-phasic peaking in the evening and morning, respectively. However, in the heart, these same regulators peak and trough with a different timing, indicating that the liver and heart are not synchronous. The digestive tract displays yet another difference: simultaneous phases in the expression of positive and negative clock regulators, and we do not observe significant rhythms in clock gene expression in the retina. Our data show that there is a functional clock in Chinook salmon tissues, but that this clock behaves in a tissue-specific manner, regardless of the whole animal being exposed to the same environmental cues. These results highlight the adaptive role of the clock in Chinook salmon and that it may have different positive and negative effects depending on tissue function.

Inter-generational resemblance of methylation levels at circadian genes and associations with phenology in the barn swallow

Regulation of gene expression can occur via epigenetic effects as mediated by DNA methylation. The potential for epigenetic effects to be transmitted across generations, thus modulating phenotypic variation and affecting ecological and evolutionary processes, is increasingly appreciated. However, the study of variation in epigenomes and inter-generational transmission of epigenetic alterations in wild populations is at its very infancy. We studied sex- and age-related variation in DNA methylation and parent-offspring resemblance in methylation profiles in the barn swallows. We focused on a class of highly conserved 'clock' genes (clock, cry1, per2, per3, timeless) relevant in the timing of activities of major ecological importance. In addition, we considerably expanded previous analyses on the relationship between methylation at clock genes and breeding date, a key fitness trait in barn swallows. We found positive assortative mating for methylation at one clock locus. Methylation varied between the nestling and the adult stage, and according to sex. Individuals with relatively high methylation as nestlings also had high methylation levels when adults. Extensive parent-nestling resemblance in methylation levels was observed. Occurrence of extra-pair fertilizations allowed to disclose evidence hinting at a prevalence of paternal germline or sperm quality effects over common environment effects in generating father-offspring resemblance in methylation. Finally, we found an association between methylation at the clock poly-Q region, but not at other loci, and breeding date. We thus provided evidence for sex-dependent variation and the first account of parent-offspring resemblance in methylation in any wild vertebrate. We also showed that epigenetics may influence phenotypic plasticity of timing of life cycle events, thus having a major impact on fitness.

The genetics and epigenetics of animal migration and orientation: birds, butterflies and beyond

Migration is a complex behavioural adaptation for survival that has evolved across the animal kingdom from invertebrates to mammals. In some taxa, closely related migratory species, or even populations of the same species, exhibit different migratory phenotypes, including timing and orientation of migration. In these species, a significant proportion of the phenotypic variance in migratory traits is genetic. In others, the migratory phenotype and direction is triggered by seasonal changes in the environment, suggesting an epigenetic control of their migration. The genes and epigenetic changes underpinning migratory behaviour remain largely unknown. The revolution in (epi)genomics and functional genomic tools holds great promise to rapidly move the field of migration genetics forward. Here, we review our current understanding of the genetic and epigenetic architecture of migratory traits, focusing on two emerging models: the European blackcap and the North American monarch butterfly. We also outline a vision of how technical advances and integrative approaches could be employed to identify and functionally validate candidate genes and cis-regulatory elements on these and other migratory species across both small and broad phylogenetic scales to significantly advance the field of genetics of animal migration.

Non-Pleiotropic Coupling of Daily and Seasonal Temporal Isolation in the European Corn Borer

Speciation often involves the coupling of multiple isolating barriers to produce reproductive isolation, but how coupling is generated among different premating barriers is unknown. We measure the degree of coupling between the daily mating time and seasonal mating time between strains of European corn borer (Ostrinia nubilalis) and evaluate the hypothesis that the coupling of different forms of allochrony is due to a shared genetic architecture, involving genes with pleiotropic effects on both timing phenotypes. We measure differences in gene expression at peak mating times and compare these genes to previously identified candidates that are associated with changes in seasonal mating time between the corn borer strains. We find that the E strain, which mates earlier in the season, also mates 2.7 h earlier in the night than the Z strain. Earlier daily mating is correlated with the differences in expression of the circadian clock genes cycle, slimb, and vrille. However, different circadian clock genes associate with daily and seasonal timing, suggesting that the coupling of timing traits is maintained by natural selection rather than pleiotropy. Juvenile hormone gene expression was associated with both types of timing, suggesting that circadian genes activate common downstream modules that may impose constraint on future evolution of these traits.

Examination of Clock and Adcyap1 gene variation in a neotropical migratory passerine

Complex behavioral traits, such as those making up a migratory phenotype, are regulated by multiple environmental factors and multiple genes. We investigated possible relationships between microsatellite variation at two candidate genes implicated in the control of migratory behavior, Clock and Adcyap1, and several aspects of migratory life-history and evolutionary divergence in the Painted Bunting (Passerina ciris), a species that shows wide variation in migratory and molting strategies across a disjunct distribution. We focused on Clock and Adcyap1 microsatellite variation across three Painted Bunting populations in Oklahoma, Louisiana, and North Carolina, and for the Oklahoma breeding population we used published migration tracking data on adult males to explore phenotypic variation in individual migratory behavior. We found no correlation between microsatellite allele size within either Clock and Adcyap1 relative to the initiation or duration of fall migration in adult males breeding in Oklahoma. We also show the lack of significant correlations with aspects of the migratory phenotype for the Louisiana population. Our research highlights the limitations of studying microsatellite allelic mutations that are of undetermined functional influence relative to complex behavioral phenotypes.

Divergence in expression of candidate genes for the smoltification process between juvenile resident rainbow and anadromous steelhead trout

Rainbow and steelhead trout (Oncorhynchus mykiss), among other salmonid fishes, exhibit tremendous life history diversity, foremost of which is variation in migratory propensity. While some individuals possess the ability to undertake an anadromous marine migration, others remain resident in freshwater throughout their life cycle. Those that will migrate undergo tremendous physiological, morphological, and behavioral transformations in a process called smoltification which transitions freshwater-adapted parr to marine-adapted smolts. While the behavior, ecology, and physiology of smoltification are well described, our understanding of the proximate genetic mechanisms that trigger the process are not well known. Quantitative genetic analyses have identified several genomic regions associated with smoltification and migration-related traits within this species. Here we investigate the divergence in gene expression of 18 functional and positional candidate genes for the smoltification process in the brain, gill, and liver tissues of migratory smolts, resident parr, and precocious mature male trout at the developmental stage of out-migration. Our analysis reveals several genes differentially expressed between life history classes and validates the candidate nature of several genes in the parr-smolt transformation including Clock1α, FSHβ, GR, GH2, GHR1, GHR2, NDK7, p53, SC6a7, Taldo1, THRα, THRβ, and Vdac2.

Time scale matters: genetic analysis does not support adaptation-by-time as the mechanism for adaptive seasonal declines in kokanee reproductive life span

Seasonal declines of fitness-related traits are often attributed to environmental effects or individual-level decisions about reproductive timing and effort, but genetic variation may also play a role. In populations of Pacific salmon (Oncorhynchus spp.), seasonal declines in reproductive life span have been attributed to adaptation-by-time, in which divergent selection for different traits occurs among reproductively isolated temporal components of a population. We evaluated this hypothesis in kokanee (freshwater obligate Oncorhynchus nerka) by testing for temporal genetic structure in neutral and circadian-linked loci. We detected no genetic differences in presumably neutral loci among kokanee with different arrival and maturation dates within a spawning season. Similarly, we detected no temporal genetic structure in OtsClock1b, Omy1009uw, or OmyFbxw11, candidate loci associated with circadian function. The genetic evidence from this study and others indicates a lack of support for adaptation-by-time as an important evolutionary mechanism underlying seasonal declines in reproductive life span and a need for greater consideration of other mechanisms such as time-dependent, adaptive adjustment of reproductive effort.

Testing advances in molecular discrimination among Chinook salmon life histories: evidence from a blind test

The application of DNA-based markers toward the task of discriminating among alternate salmon runs has evolved in accordance with ongoing genomic developments and increasingly has enabled resolution of which genetic markers associate with important life-history differences. Accurate and efficient identification of the most likely origin for salmon encountered during ocean fisheries, or at salvage from fresh water diversion and monitoring facilities, has far-reaching consequences for improving measures for management, restoration and conservation. Near-real-time provision of high-resolution identity information enables prompt response to changes in encounter rates. We thus continue to develop new tools to provide the greatest statistical power for run identification. As a proof of concept for genetic identification improvements, we conducted simulation and blind tests for 623 known-origin Chinook salmon (Oncorhynchus tshawytscha) to compare and contrast the accuracy of different population sampling baselines and microsatellite loci panels. This test included 35 microsatellite loci (1266 alleles), some known to be associated with specific coding regions of functional significance, such as the circadian rhythm cryptochrome genes, and others not known to be associated with any functional importance. The identification of fall run with unprecedented accuracy was demonstrated. Overall, the top performing panel and baseline (HMSC21) were predicted to have a success rate of 98%, but the blind-test success rate was 84%. Findings for bias or non-bias are discussed to target primary areas for further research and resolution.

Adaptive genetic markers discriminate migratory runs of Chinook salmon (Oncorhynchus tshawytscha) amid continued gene flow

Neutral genetic markers are routinely used to define distinct units within species that warrant discrete management. Human-induced changes to gene flow however may reduce the power of such an approach. We tested the efficiency of adaptive versus neutral genetic markers in differentiating temporally divergent migratory runs of Chinook salmon (Oncorhynchus tshawytscha) amid high gene flow owing to artificial propagation and habitat alteration. We compared seven putative migration timing genes to ten microsatellite loci in delineating three migratory groups of Chinook in the Feather River, CA: offspring of fall-run hatchery broodstock that returned as adults to freshwater in fall (fall run), spring-run offspring that returned in spring (spring run), and fall-run offspring that returned in spring (FRS). We found evidence for significant differentiation between the fall and federally listed threatened spring groups based on divergence at three circadian clock genes (OtsClock1b, OmyFbxw11, and Omy1009UW), but not neutral markers. We thus demonstrate the importance of genetic marker choice in resolving complex life history types. These findings directly impact conservation management strategies and add to previous evidence from Pacific and Atlantic salmon indicating that circadian clock genes influence migration timing.

Evaluating adaptive divergence between migratory and nonmigratory ecotypes of a salmonid fish, Oncorhynchus mykiss

Next-generation sequencing and the application of population genomic and association approaches have made it possible to detect selection and unravel the genetic basis to variable phenotypic traits. The use of these two approaches in parallel is especially attractive in nonmodel organisms that lack a sequenced and annotated genome, but only works well when population structure is not confounded with the phenotype of interest. Herein, we use population genomics in a nonmodel fish species, rainbow trout (Oncorhynchus mykiss), to better understand adaptive divergence between migratory and nonmigratory ecotypes and to further our understanding about the genetic basis of migration. Restriction site-associated DNA (RAD) tag sequencing was used to identify single-nucleotide polymorphisms (SNPs) in migrant and resident O. mykiss from two systems, one in Alaska and the other in Oregon. A total of 7920 and 6755 SNPs met filtering criteria in the Alaska and Oregon data sets, respectively. Population genetic tests determined that 1423 SNPs were candidates for selection when loci were compared between resident and migrant samples. Previous linkage mapping studies that used RAD DNA tag SNPs were available to determine the position of 1990 markers. Several significant SNPs are located in genome regions that contain quantitative trait loci for migratory-related traits, reinforcing the importance of these regions in the genetic basis of migration/residency. Annotation of genome regions linked to significant SNPs revealed genes involved in processes known to be important in migration (such as osmoregulatory function). This study adds to our growing knowledge on adaptive divergence between migratory and nonmigratory ecotypes of this species; across studies, this complex trait appears to be controlled by many loci of small effect, with some in common, but many loci not shared between populations studied.

Application of single nucleotide polymorphisms to non-model species: a technical review

Single nucleotide polymorphisms (SNPs) have gained wide use in humans and model species and are becoming the marker of choice for applications in other species. Technology that was developed for work in model species may provide useful tools for SNP discovery and genotyping in non-model organisms. However, SNP discovery can be expensive, labour intensive, and introduce ascertainment bias. In addition, the most efficient approaches to SNP discovery will depend on the research questions that the markers are to resolve as well as the focal species. We discuss advantages and disadvantages of several past and recent technologies for SNP discovery and genotyping and summarize a variety of SNP discovery and genotyping studies in ecology and evolution.

Temporal patterns of genetic variation in a salmon population undergoing rapid change in migration timing

Though genetic diversity is necessary for population persistence in rapidly changing environments, little is known about how climate-warming influences patterns of intra-population genetic variation. For a pink salmon population experiencing increasing temperatures, we used temporal genetic data (microsatellite = 1993, 2001, 2009; allozyme = 1979, 1981, 1983) to quantify the genetic effective population size (Ne ) and genetic divergence due to differences in migration timing and to estimate whether these quantities have changed over time. We predicted that temporal trends toward earlier migration timing and a corresponding loss of phenotypic variation would decrease genetic divergence based on migration timing and Ne . We observed significant genetic divergence based on migration timing and genetic heterogeneity between early- and late-migrating fish. There was also some evidence for divergent selection between early- and late-migrating fish at circadian rhythm genes, but results varied over time. Estimates of Ne from multiple methods were large (>1200) and Ne /Nc generally exceeded 0.2. Despite shifts in migration timing and loss of phenotypic variation, there was no evidence for changes in within-population genetic divergence or Ne over the course of this study. These results suggest that in instances of population stability, genetic diversity may be resistant to climate-induced changes in migration timing.

Variation in candidate genes CLOCK and ADCYAP1 does not consistently predict differences in migratory behavior in the songbird genus Junco

Recent studies exploring the molecular genetic basis for migratory variation in animals have identified polymorphisms in two genes ( CLOCK and ADCYAP1) that are linked to circadian rhythms and correlate with migratory propensity and phenology among individuals and populations. Results from these initial studies are mixed, however, and additional data are needed to assess the generality and diversity of the molecular mechanisms that regulate the biology of migration. We sequenced CLOCK and ADCYAP1 in 15 populations across the two species of the avian genus Junco, a North American lineage in which multiple recently diverged subspecies and populations range from sedentary to long-distance migrants. We found no consistent associations between allele length and migratory status across the genus for either CLOCK or ADCYAP1. However, within two subspecies groups, populations that migrate longer distances have longer CLOCK alleles on average. Additionally, there was a positive relationship between ADCYAP1 allele length and migratory restlessness (zugunruhe) among individuals within one of two captive populations studied—a result similar to those reported previously within captive blackcaps ( Sylvia atricapilla). We conclude that, while both ADCYAP1 and CLOCK may correlate with migratory propensity within or among certain populations or species, previously identified relationships between migratory behavior and sequence variants cannot be easily generalized across taxa.

Cost-effective genome-wide estimation of allele frequencies from pooled DNA in Atlantic salmon (Salmo salar L.)

Background

New sequencing technologies have tremendously increased the number of known molecular markers (single nucleotide polymorphisms; SNPs) in a variety of species. Concurrently, improvements to genotyping technology have now made it possible to efficiently genotype large numbers of genome-wide distributed SNPs enabling genome wide association studies (GWAS). However, genotyping significant numbers of individuals with large number of SNPs remains prohibitively expensive for many research groups. A possible solution to this problem is to determine allele frequencies from pooled DNA samples, such ‘allelotyping’ has been presented as a cost-effective alternative to individual genotyping and has become popular in human GWAS. In this article we have tested the effectiveness of DNA pooling to obtain accurate allele frequency estimates for Atlantic salmon (Salmo salar L.) populations using an Illumina SNP-chip.

Results

In total, 56 Atlantic salmon DNA pools from 14 populations were analyzed on an Atlantic salmon SNP-chip containing probes for 5568 SNP markers, 3928 of which were bi-allelic. We developed an efficient quality control filter which enables exclusion of loci showing high error rate and minor allele frequency (MAF) close to zero. After applying multiple quality control filters we obtained allele frequency estimates for 3631 bi-allelic loci. We observed high concordance (r > 0.99) between allele frequency estimates derived from individual genotyping and DNA pools. Our results also indicate that even relatively small DNA pools (35 individuals) can provide accurate allele frequency estimates for a given sample.

Conclusions

Despite of higher level of variation associated with array replicates compared to pool construction, we suggest that both sources of variation should be taken into account. This study demonstrates that DNA pooling allows fast and high-throughput determination of allele frequencies in Atlantic salmon enabling cost-efficient identification of informative markers for discrimination of populations at various geographical scales, as well as identification of loci controlling ecologically and economically important traits.

Genetic change for earlier migration timing in a pink salmon population

To predict how climate change will influence populations, it is necessary to understand the mechanisms, particularly microevolution and phenotypic plasticity, that allow populations to persist in novel environmental conditions. Although evidence for climate-induced phenotypic change in populations is widespread, evidence documenting that these phenotypic changes are due to microevolution is exceedingly rare. In this study, we use 32 years of genetic data (17 complete generations) to determine whether there has been a genetic change towards earlier migration timing in a population of pink salmon that shows phenotypic change; average migration time occurs nearly two weeks earlier than it did 40 years ago. Experimental genetic data support the hypothesis that there has been directional selection for earlier migration timing, resulting in a substantial decrease in the late-migrating phenotype (from more than 30% to less than 10% of the total abundance). From 1983 to 2011, there was a significant decrease--over threefold--in the frequency of a genetic marker for late-migration timing, but there were minimal changes in allele frequencies at other neutral loci. These results demonstrate that there has been rapid microevolution for earlier migration timing in this population. Circadian rhythm genes, however, did not show any evidence for selective changes from 1993 to 2009.

Detection of outlier loci and their utility for fisheries management

Genetics-based approaches have informed fisheries management for decades, yet remain challenging to implement within systems involving recently diverged stocks or where gene flow persists. In such cases, genetic markers exhibiting locus-specific (‘outlier’) effects associated with divergent selection may provide promising alternatives to loci that reflect genome-wide (‘neutral’) effects for guiding fisheries management. Okanagan Lake kokanee (Oncorhynchus nerka), a fishery of conservation concern, exhibits two sympatric ecotypes adapted to different reproductive environments; however, previous research demonstrated the limited utility of neutral microsatellites for assigning individuals. Here, we investigated the efficacy of an outlier-based approach to fisheries management by screening >11 000 expressed sequence tags for linked microsatellites and conducting genomic scans for kokanee sampled across seven spawning sites. We identified eight outliers among 52 polymorphic loci that detected ecotype-level divergence, whereas there was no evidence of divergence at neutral loci. Outlier loci exhibited the highest self-assignment accuracy to ecotype (92.1%), substantially outperforming 44 neutral loci (71.8%). Results were robust among-sampling years, with assignment and mixed composition estimates for individuals sampled in 2010 mirroring baseline results. Overall, outlier loci constitute promising alternatives for informing fisheries management involving recently diverged stocks, with potential applications for designating management units across a broad range of taxa.

The genetics of migration on the move

Across a range of organisms, related species or even populations of the same species exhibit strikingly different scales and patterns of movement. A significant proportion of the phenotypic variance in migratory traits is genetic, but the genes involved in shaping these phenotypes are still unknown. Although recent achievements in genomics will evolve migratory genetics research from a phenotypic to a molecular approach, fully sequenced and annotated genomes of migratory species are still lacking. Consequently, many of the genes involved in migration are unavailable as candidates. Migration is central to the life-history adaptations of many animals. Here, we review current understanding of the genetic architecture of migratory traits and discuss the significant implications this will have for other areas of biology, including population responses to climate change, speciation and conservation management.

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.

Detecting population structure in a high gene-flow species, Atlantic herring (Clupea harengus): direct, simultaneous evaluation of neutral vs putatively selected loci

In many marine fish species, genetic population structure is typically weak because populations are large, evolutionarily young and have a high potential for gene flow. We tested whether genetic markers influenced by natural selection are more efficient than the presumed neutral genetic markers to detect population structure in Atlantic herring (Clupea harengus), a migratory pelagic species with large effective population sizes. We compared the spatial and temporal patterns of divergence and statistical power of three traditional genetic marker types, microsatellites, allozymes and mitochondrial DNA, with one microsatellite locus, Cpa112, previously shown to be influenced by divergent selection associated with salinity, and one locus located in the major histocompatibility complex class IIA (MHC-IIA) gene, using the same individuals across analyses. Samples were collected in 2002 and 2003 at two locations in the North Sea, one location in the Skagerrak and one location in the low-saline Baltic Sea. Levels of divergence for putatively neutral markers were generally low, with the exception of single outlier locus/sample combinations; microsatellites were the most statistically powerful markers under neutral expectations. We found no evidence of selection acting on the MHC locus. Cpa112, however, was highly divergent in the Baltic samples. Simulations addressing the statistical power for detecting population divergence showed that when using Cpa112 alone, compared with using eight presumed neutral microsatellite loci, sample sizes could be reduced by up to a tenth while still retaining high statistical power. Our results show that the loci influenced by selection can serve as powerful markers for detecting population structure in high gene-flow marine fish species.

Extent and scale of local adaptation in salmonid fishes: review and meta-analysis

What is the extent and scale of local adaptation (LA)? How quickly does LA arise? And what is its underlying molecular basis? Our review and meta-analysis on salmonid fishes estimates the frequency of LA to be ∼55-70%, with local populations having a 1.2 times average fitness advantage relative to foreign populations or to their performance in new environments. Salmonid LA is evident at a variety of spatial scales (for example, few km to>1000 km) and can manifest itself quickly (6-30 generations). As the geographic scale between populations increases, LA is generally more frequent and stronger. Yet the extent of LA in salmonids does not appear to differ from that in other assessed taxa. Moreover, the frequency with which foreign salmonid populations outperform local populations (∼23-35%) suggests that drift, gene flow and plasticity often limit or mediate LA. The relatively few studies based on candidate gene and genomewide analyses have identified footprints of selection at both small and large geographical scales, likely reflecting the specific functional properties of loci and the associated selection regimes (for example, local niche partitioning, pathogens, parasites, photoperiodicity and seasonal timing). The molecular basis of LA in salmonids is still largely unknown, but differential expression at the same few genes is implicated in the convergent evolution of certain phenotypes. Collectively, future research will benefit from an integration of classical and molecular approaches to understand: (i) species differences and how they originate, (ii) variation in adaptation across scales, life stages, population sizes and environmental gradients, and (iii) evolutionary responses to human activities.

Correlated contemporary evolution of life history traits in New Zealand Chinook salmon, Oncorhynchus tshawytscha

Size at age and age at maturity are important life history traits, affecting individual fitness and population demography. In salmon and other organisms, size and growth rate are commonly considered cues for maturation and thus age at maturity may or may not evolve independently of these features. Recent concerns surrounding the potential phenotypic and demographic responses of populations facing anthropogenic disturbances, such as climate change and harvest, place a premium on understanding the evolutionary genetic basis for evolution in size at age and age at maturity. In this study, we present the findings from a set of common-garden rearing experiments that empirically assess the heritable basis of phenotypic divergence in size at age and age at maturity in Chinook salmon (Oncorhynchus tshawytscha) populations introduced to New Zealand. We found consistent evidence of heritable differences among populations in both size at age and age at maturity, often corresponding to patterns observed in the wild. Populations diverged in size and growth profiles, even when accounting for eventual age at maturation. By contrast, most, but not all, cases of divergence in age at maturity were driven by the differences in size or growth rate rather than differences in the threshold relationship linking growth rate and probability of maturation. These findings help us understand how life histories may evolve through trait interactions in populations exposed to natural and anthropogenic disturbances, and how we might best detect such evolution.

Clock genes and their genomic distributions in three species of salmonid fishes: Associations with genes regulating sexual maturation and cell cycling

Background

Clock family genes encode transcription factors that regulate clock-controlled genes and thus regulate many physiological mechanisms/processes in a circadian fashion. Clock1 duplicates and copies of Clock3 and NPAS2-like genes were partially characterized (genomic sequencing) and mapped using family-based indels/SNPs in rainbow trout (RT)(Oncorhynchus mykiss), Arctic charr (AC)(Salvelinus alpinus), and Atlantic salmon (AS)(Salmo salar) mapping panels.

Results

Clock1 duplicates mapped to linkage groups RT-8/-24, AC-16/-13 and AS-2/-18. Clock3/NPAS2-like genes mapped to RT-9/-20, AC-20/-43, and AS-5. Most of these linkage group regions containing the Clock gene duplicates were derived from the most recent 4R whole genome duplication event specific to the salmonids. These linkage groups contain quantitative trait loci (QTL) for life history and growth traits (i.e., reproduction and cell cycling). Comparative synteny analyses with other model teleost species reveal a high degree of conservation for genes in these chromosomal regions suggesting that functionally related or co-regulated genes are clustered in syntenic blocks. For example, anti-müllerian hormone (amh), regulating sexual maturation, and ornithine decarboxylase antizymes (oaz1 and oaz2), regulating cell cycling, are contained within these syntenic blocks.

Conclusions

Synteny analyses indicate that regions homologous to major life-history QTL regions in salmonids contain many candidate genes that are likely to influence reproduction and cell cycling. The order of these genes is highly conserved across the vertebrate species examined, and as such, these genes may make up a functional cluster of genes that are likely co-regulated. CLOCK, as a transcription factor, is found within this block and therefore has the potential to cis-regulate the processes influenced by these genes. Additionally, clock-controlled genes (CCGs) are located in other life-history QTL regions within salmonids suggesting that at least in part, trans-regulation of these QTL regions may also occur via Clock expression.

Clock polymorphism in Pacific salmon: evidence for variable selection along a latitudinal gradient

Seasonal timing of life-history events is often under strong natural selection. The Clock gene is a central component of an endogenous circadian clock that senses changes in photoperiod (day length) and mediates seasonal behaviours. Among Pacific salmonids (Oncorhynchus spp.), seasonal timing of migration and breeding is influenced by photoperiod. To expand a study of 42 North American Chinook salmon (Oncorhynchus tshawytscha) populations, we tested whether duplicated Clock genes contribute to population differences in reproductive timing. Specifically, we examined geographical variation along a similar latitudinal gradient in the polyglutamine domain (PolyQ) of OtsClock1a and OtsClock1b among 53 populations of three species: chum (Oncorhynchus keta), coho (Oncorhynchus kisutch) and pink salmon (Oncorhynchus gorbuscha). We found evidence for variable selection on OtsClock1b that corresponds to latitudinal variation in reproductive timing among these species. We evaluated the contribution of day length and a freshwater migration index to OtsClock1b PolyQ domain variation using regression trees and found that day length at spawning explains much of the variation in OtsClock1b allele frequency among chum and Chinook, but not coho and pink salmon populations. Our findings suggest that OtsClock1b mediates seasonal adaptation and influences geographical variation in reproductive timing in some of these highly migratory species.

Clock genes localize to quantitative trait loci for stage-specific growth in juvenile coho salmon, Oncorhynchus kisutch

In most organisms, an internal circadian clock coordinates the expression of biological rhythms and enables individuals to anticipate and respond to the seasonally changing environment. There is remarkable conservation of function in the molecular machinery underlying this circadian clock across taxa with 4 canonical proteins interacting to form an autoregulatory feedback loop: CLOCK, CRYPTOCHROME, PERIOD, and BMAL. We mapped duplicated copies of Clock and Cryptochrome in coho salmon (Oncorhynchus kisutch) to determine if these genes localize to quantitative trait loci (QTL) for hatch timing, weight, length, and growth rate measured throughout the juvenile life-history stage. We found that Cryptochrome2b mapped to a QTL region for growth (measured at 304 days post-hatching) on linkage group OKI06. The percentage of variation (PEV) explained by this QTL was 15.2%. Cryptochrome2b was also associated with a marginally nonsignificant QTL for length (measured at 395 days post-hatching). OtsClock1b mapped to a QTL region for growth rate (PEV 10.1%) and length (PEV 10.5%) on linkage group OKI24 (measured at 479 days posthatching). Neither gene localized to QTL for hatch timing or weight. Our findings indicate that the growth rate and length QTL associated with OtsClock1b and Cryptochrome2b are development stage-specific and may result from temporally differentiated gene expression patterns.

Current knowledge on the photoneuroendocrine regulation of reproduction in temperate fish species

Seasonality is an important adaptive trait in temperate fish species as it entrains or regulates most physiological events such as reproductive cycle, growth profile, locomotor activity and key life-stage transitions. Photoperiod is undoubtedly one of the most predictable environmental signals that can be used by most living organisms including fishes in temperate areas. This said, however, understanding of how such a simple signal can dictate the time of gonadal recruitment and spawning, for example, is a complex task. Over the past few decades, many scientists attempted to unravel the roots of photoperiodic signalling in teleosts by investigating the role of melatonin in reproduction, but without great success. In fact, the hormone melatonin is recognized as the biological time-keeping hormone in fishes mainly due to the fact that it reflects the seasonal variation in daylength across the whole animal kingdom rather than the existence of direct evidences of its role in the entrainment of reproduction in fishes. Recently, however, some new studies clearly suggested that melatonin interacts with the reproductive cascade at a number of key steps such as through the dopaminergic system in the brain or the synchronization of the final oocyte maturation in the gonad. Interestingly, in the past few years, additional pathways have become apparent in the search for a fish photoneuroendocrine system including the clock-gene network and kisspeptin signalling and although research on these topics are still in their infancy, it is moving at great pace. This review thus aims to bring together the current knowledge on the photic control of reproduction mainly focusing on seasonal temperate fish species and shape the current working hypotheses supported by recent findings obtained in teleosts or based on knowledge gathered in mammalian and avian species. Four of the main potential regulatory systems (light perception, melatonin, clock genes and kisspeptin) in fish reproduction are reviewed.

Contemporary isolation-by-distance, but not isolation-by-time, among demes of European Grayling (Thymallus thymallus, Linnaeus) with recent common ancestors

The development of isolation by distance (IBD) and isolation by time (IBT) was contrasted among demes of European grayling (Thymallus thymallus) that have diverged within the last 25 generations following colonization of a lake (Lesjaskogsvatnet). We find low but significant levels of genetic differentiation among spawning tributaries and a pattern of IBD among them. We do not, however, find evidence for IBT despite an up to four-week difference in spawning date between "warm/early" and "cold/late" spawning demes and differences in the incubation temperatures experienced by offspring. It appears that IBD has developed more rapidly than IBT in this system and that adaptive divergence has been initiated in the absence of IBT. Although analysis of selected loci could reveal reduced recombination in parts of the genome associated with temporal divergence, our analysis of neutral genetic data suggests that IBD is a more important isolating mechanism in the early stages of adaptive divergence in European grayling.

A latitudinal cline in the Chinook salmon (Oncorhynchus tshawytscha) Clock gene: evidence for selection on PolyQ length variants

A critical seasonal event for anadromous Chinook salmon (Oncorhynchus tshawytscha) is the time at which adults migrate from the ocean to breed in freshwater. We investigated whether allelic variation at the circadian rhythm genes, OtsClock1a and OtsClock1b, underlies genetic control of migration timing among 42 populations in North America. We identified eight length variants of the functionally important polyglutamine repeat motif (PolyQ) of OtsClock1b while OtsClock1a PolyQ was highly conserved. We found evidence of a latitudinal cline in average allele length and frequency of the two most common OtsClock1b alleles. The shorter 335 bp allele increases in frequency with decreasing latitude while the longer 359 bp allele increases in frequency at higher latitudes. Comparison to 13 microsatellite loci showed that 335 and 359 bp deviate significantly from neutral expectations. Furthermore, a hierarchical gene diversity analysis based on OtsClock1b PolyQ variation revealed that run timing explains 40.9 per cent of the overall genetic variance among populations. By contrast, an analysis based on 13 microsatellite loci showed that run timing explains only 13.2 per cent of the overall genetic variance. Our findings suggest that length polymorphisms in OtsClock1b PolyQ may be maintained by selection and reflect an adaptation to ecological factors correlated with latitude, such as the seasonally changing day length.