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

Gap-free genome assembly of anadromous Coilia nasus

The Chinese tapertail anchovy, Coilia nasus, is a socioeconomically important anadromous fish that migrates from near ocean waters to freshwater to spawn every spring. The analysis of genomic architecture and information of C. nasus were hindered by the previously released versions of reference genomes with gaps. Here, we report the assembly of a chromosome-level gap-free genome of C. nasus by incorporating high-coverage and accurate long-read sequence data with multiple assembly strategies. All 24 chromosomes were assembled without gaps, representing the highest completeness and assembly quality. We assembled the genome with a size of 851.67 Mb and used BUSCO to estimate the completeness of the assembly as 92.5%. Using a combination of de novo prediction, protein homology and RNA-seq annotation, 21,900 genes were functionally annotated, representing 99.68% of the total predicted protein-coding genes. The availability of gap-free reference genomes for C. nasus will provide the opportunity for understanding genome structure and function, and will also lay a solid foundation for further management and conservation of this important species.

Microsatellites as Molecular Markers with Applications in Exploitation and Conservation of Aquatic Animal Populations

A large number of species and taxa has been studied for genetic polymorphism. Microsatellites have been known as hypervariable neutral molecular markers with the highest resolution power in comparison with any other markers. However, the discovery of a new type of molecular marker—single nucleotide polymorphism (SNP) has put the existing applications of microsatellites to the test. To ensure good resolution power in studies of populations and individuals, a number of microsatellite loci from 14 to 20 was often used, which corresponds to about 200 independent alleles. Recently, these numbers have tended to be increased by the application of genomic sequencing of expressed sequence tags (ESTs), and the choice of the most informative loci for genotyping depends on the aims of research. Examples of successful applications of microsatellite molecular markers in aquaculture, fisheries, and conservation genetics in comparison with SNPs have been summarized in this review. Microsatellites can be considered superior markers in such topics as kinship and parentage analysis in cultured and natural populations, the assessment of gynogenesis, androgenesis and ploidization. Microsatellites can be coupled with SNPs for mapping QTL. Microsatellites will continue to be used in research on genetic diversity in cultured stocks, and also in natural populations as an economically advantageous genotyping technique.

Sex-differences in Phenology: A Tinbergian Perspective

Shifts in the timing of cyclic seasonal life-history events are among the most commonly reported responses to climate change, with differences in response rates among interacting species leading to phenological mismatches. Within a species, however, males and females can also exhibit differential sensitivity to environmental cues and may therefore differ in their responsiveness to climate change, potentially leading to phenological mismatches between the sexes. This occurs because males differ from females in when and how energy is allocated to reproduction, resulting in marked sex-differences in life-history timing across the annual cycle. In this review, we take a Tinbergian perspective and examine sex differences in timing of vertebrates from adaptive, ontogenetic, mechanistic, and phylogenetic viewpoints with the goal of informing and motivating more integrative research on sexually dimorphic phenologies. We argue that sexual and natural selection lead to sex-differences in life-history-timing and that understanding the ecological and evolutionary drivers of these differences is critical for connecting climate-driven phenological shifts to population resilience. Ontogeny may influence how and when sex differences in life-history timing arise because the early-life environment can profoundly affect developmental trajectory, rates of reproductive maturation, and seasonal timing. The molecular mechanisms underlying these organismal traits are relevant to identifying the diversity and genetic basis of population- and species-level responses to climate change, and promisingly, the molecular basis of phenology is becoming increasingly well-understood. However, because most studies focus on a single sex, the causes of sex-differences in phenology critical to population resilience often remain unclear. New sequencing tools and analyses informed by phylogeny may help generate hypotheses about mechanism as well as insight into the general "evolvability" of sex differences across phylogenetic scales, especially as trait and genome resources grow. We recommend that greater attention be placed on determining sex-differences in timing mechanisms and monitoring climate change responses in both sexes, and we discuss how new tools may provide key insights into sex-differences in phenology from all four Tinbergian domains.

Non-Lethal Sampling Supports Integrative Movement Research in Freshwater Fish

Freshwater ecosystems and fishes are enormous resources for human uses and biodiversity worldwide. However, anthropogenic climate change and factors such as dams and environmental contaminants threaten these freshwater systems. One way that researchers can address conservation issues in freshwater fishes is via integrative non-lethal movement research. We review different methods for studying movement, such as with acoustic telemetry. Methods for connecting movement and physiology are then reviewed, by using non-lethal tissue biopsies to assay environmental contaminants, isotope composition, protein metabolism, and gene expression. Methods for connecting movement and genetics are reviewed as well, such as by using population genetics or quantitative genetics and genome-wide association studies. We present further considerations for collecting molecular data, the ethical foundations of non-lethal sampling, integrative approaches to research, and management decisions. Ultimately, we argue that non-lethal sampling is effective for conducting integrative, movement-oriented research in freshwater fishes. This research has the potential for addressing critical issues in freshwater systems in the future.

Genomic region associated with run timing has similar haplotypes and phenotypic effects across three lineages of Chinook salmon

Conserving life-history variation is a stated goal of many management programs, but the most effective means by which to accomplish this are often far from clear. Early- and late-migrating forms of Chinook salmon (Oncorhynchus tshawytscha) face unequal pressure from natural and anthropogenic forces that may alter the impacts of genetic variation underlying heritable migration timing. Genomic regions of chromosome 28 are known to be strongly associated with migration variation in adult Chinook salmon, but it remains unclear whether there is consistent association among diverse lineages and populations in large basins such as the Columbia River. With high-throughput genotyping (GT-seq) and phenotyping methods, we examined the association of genetic variation in 28 markers (spanning GREB1L to ROCK1 of chromosome 28) with individual adult migration timing characteristics gleaned from passive integrated transponder recordings of over 5000 Chinook salmon from the three major phylogeographic lineages that inhabit the Columbia River Basin. Despite the strong genetic differences among them in putatively neutral genomic regions, each of the three lineages exhibited very similar genetic variants in the chromosome 28 region that were significantly associated with adult migration timing phenotypes. This is particularly notable for the interior stream-type lineage, which exhibits an earlier and more constrained freshwater entry than the other lineages. In both interior stream-type and interior ocean-type lineages of Chinook salmon, heterozygotes of the most strongly associated linkage groups had largely intermediate migration timing relative to homozygotes, and results indicate codominance or possibly marginal partial dominance of the allele associated with early migration. Our results lend support to utilization of chromosome 28 variation in tracking and predicting run timing in these lineages of Chinook salmon in the Columbia River.

Population differences in Chinook salmon (Oncorhynchus tshawytscha) DNA methylation: Genetic drift and environmental factors

Local adaptation and phenotypic differences among populations have been reported in many species, though most studies focus on either neutral or adaptive genetic differentiation. With the discovery of DNA methylation, questions have arisen about its contribution to individual variation in and among natural populations. Previous studies have identified differences in methylation among populations of organisms, although most to date have been in plants and model animal species. Here we obtained eyed eggs from eight populations of Chinook salmon (Oncorhynchus tshawytscha) and assayed DNA methylation at 23 genes involved in development, immune function, stress response, and metabolism using a gene-targeted PCR-based assay for next-generation sequencing. Evidence for population differences in methylation was found at eight out of 23 gene loci after controlling for developmental timing in each individual. However, we found no correlation between freshwater environmental parameters and methylation variation among populations at those eight genes. A weak correlation was identified between pairwise DNA methylation dissimilarity among populations and pairwise F ST based on 15 microsatellite loci, indicating weak effects of genetic drift or geographic distance on methylation. The weak correlation was primarily driven by two genes, GTIIBS and Nkef. However, single-gene Mantel tests comparing methylation and pairwise F ST were not significant after Bonferroni correction. Thus, population differences in DNA methylation are more likely related to unmeasured oceanic environmental conditions, local adaptation, and/or genetic drift. DNA methylation is an additional mechanism that contributes to among population variation, with potential influences on organism phenotype, adaptive potential, and population resilience.

Validation and association of candidate markers for adult migration timing and fitness in Chinook Salmon

Recent studies have begun to elucidate the genetic basis for phenotypic traits in salmonid species, but many questions remain before these candidate genes can be directly incorporated into conservation management. In Chinook Salmon (Oncorhynchus tshawytscha), a region of major effect for migration timing has been discovered that harbors two adjacent candidate genes (greb1L, rock1), but there has been limited work to examine the association between these genes and migratory phenotypes at the individual, compared to the population, level. To provide a more thorough test of individual phenotypic association within lineages of Chinook Salmon, 33 candidate markers were developed across a 220 Kb region on chromosome 28 previously associated with migration timing. Candidate and neutral markers were genotyped in individuals from representative collections that exhibit phenotypic variation in timing of arrival to spawning grounds from each of three lineages of Chinook Salmon. Association tests confirmed the majority of markers on chromosome 28 were significantly associated with arrival timing and the strongest association was consistently observed for markers within the rock1 gene and the intergenic region between greb1L and rock1. Candidate markers alone explained a wide range of phenotypic variation for Lower Columbia and Interior ocean-type lineages (29% and 78%, respectively), but less for the Interior stream-type lineage (5%). Individuals that were heterozygous at markers within or upstream of rock1 had phenotypes that suggested a pattern of dominant inheritance for early arrival across populations. Finally, previously published fitness estimates from the Interior stream-type lineage enabled tests of association with arrival timing and two candidate markers, which revealed that fish with homozygous mature genotypes had slightly higher fitness than fish with premature genotypes, while heterozygous fish were intermediate. Overall, these results provide additional information for individual-level genetic variation associated with arrival timing that may assist with conservation management of this species.

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.

Variation in DNA methylation is associated with migratory phenotypes of lake sturgeon Acipenser fulvescens in the St. Clair River, MI, USA

Lake sturgeon Acipenser fulvescens populations show a variety of movement patterns that are poorly understood. To compare two migratory phenotypes of A. fulvescens in the St. Clair River, MI, USA, multiple data types were analysed. Individual fish were classified into migratory phenotypes based on acoustic telemetry data collected 2012-2015. Acipenser fulvescens consistently showed movement from the St. Clair River upriver into Lake Huron or downriver into Lake St. Clair. The two migratory phenotypes were then compared for differences in morphometrics, genetics and epigenetics. Morphological differences based on linear measurements were not detected between phenotypes. Microsatellite data from 11 loci suggested one population with no genetic differentiation between migratory phenotypes. Our epigenetic results indicated that the migratory phenotypes are differentially methylated (P < 0.05), thus epigenetics may be associated with migratory differences in A. fulvescens. Only one restriction site was identified to be driving the differential methylation (P < 0.05). While little evidence at neutral loci occurred for genetic differentiation of A. fulvescens, DNA methylation may play a role in the observed movement pattern variation. When combined with microsatellite and morphometric analyses, our results suggested that different migratory patterns may reflect phenotypic plasticity, allowing A. fulvescens to acclimate to short-term environmental variability. Without an integrated approach, the role of epigenetics in the migratory phenotype of A. fulvescens may have been overlooked. Further characterization of migratory phenotypes could be important for management to conserve behavioural variation across the distribution of A. fulvescens and for design of stocking guidelines.

A bioenergetics evaluation of temperature-dependent selection for the spawning phenology by Snake River fall Chinook salmon

High water temperatures can increase the energetic cost for salmon to migrate and spawn, which can be important for Snake River fall-run Chinook salmon because they migrate great distances (>500 km) at a time when river temperatures (18-24°C) can be above their optimum temperatures (16.5°C). Average river temperatures and random combinations of migration and spawning dates were used to simulate fish travel times and determine the energetic consequences of different thermal experiences during migration. An energy threshold criterion (4 kJ/g) was also imposed on survival and spawning success, which was used to determine how prevailing temperatures might select against certain migration dates and thermal experiences, and in turn, explain the selection for the current spawning phenology of the population. Scenarios of tributary use for thermal refugia under increasing water temperatures (1, 2, and 3°C) were also run to determine which combinations of migration dates, travel rates, and resulting thermal experiences might be most affected by energy exhaustion. As expected, when compared to observations, the model under existing conditions and energy use could explain the onset, but not the end of the observed spawning migration. Simulations of early migrants had greater energy loss than late migrants regardless of the river temperature scenario, but higher temperatures disproportionately selected against a larger fraction of early-migrating fish, although using cold-water tributaries during migration provided a buffer against higher energy use at higher temperatures. The fraction of simulated fish that exceeded the threshold for migration success increased from 58% to 72% as average seasonal river temperatures over baseline temperatures increased. The model supports the conclusion that increases in average seasonal river temperatures as little as 1°C could impose greater thermal constraints on the fish, select against early migrants, and in turn, truncate the onset of the current spawning migration.

Genomewide association analyses of fitness traits in captive-reared Chinook salmon: Applications in evaluating conservation strategies

A novel application of genomewide association analyses is to use trait-associated loci to monitor the effects of conservation strategies on potentially adaptive genetic variation. Comparisons of fitness between captive- and wild-origin individuals, for example, do not reveal how captive rearing affects genetic variation underlying fitness traits or which traits are most susceptible to domestication selection. Here, we used data collected across four generations to identify loci associated with six traits in adult Chinook salmon (Oncorhynchus tshawytscha) and then determined how two alternative management approaches for captive rearing affected variation at these loci. Loci associated with date of return to freshwater spawning grounds (return timing), length and weight at return, age at maturity, spawn timing, and daily growth coefficient were identified using 9108 restriction site-associated markers and random forest, an approach suitable for polygenic traits. Mapping of trait-associated loci, gene annotations, and integration of results across multiple studies revealed candidate regions involved in several fitness-related traits. Genotypes at trait-associated loci were then compared between two hatchery populations that were derived from the same source but are now managed as separate lines, one integrated with and one segregated from the wild population. While no broad-scale change was detected across four generations, there were numerous regions where trait-associated loci overlapped with signatures of adaptive divergence previously identified in the two lines. Many regions, primarily with loci linked to return and spawn timing, were either unique to or more divergent in the segregated line, suggesting that these traits may be responding to domestication selection. This study is one of the first to utilize genomic approaches to demonstrate the effectiveness of a conservation strategy, managed gene flow, on trait-associated-and potentially adaptive-loci. The results will promote the development of trait-specific tools to better monitor genetic change in captive and wild populations.

Standing genetic diversity and selection at functional gene loci are associated with differential invasion success in two non-native fish species

Invasive species are expected to experience a unique combination of high genetic drift due to demographic factors while also experiencing strong selective pressures. The paradigm that reduced genetic diversity should limit the evolutionary potential of invasive species, and thus, their potential for range expansion has received little empirical support, possibly due to the choice of genetic markers. Our goal was to test for effects of genetic drift and selection at functional genetic markers as they relate to the invasion success of two paired invasive goby species, one widespread (successful) and one with limited range expansion (less successful). We genotyped fish using two marker types: single nucleotide polymorphisms (SNPs) in known-function, protein-coding genes and microsatellites to contrast the effects of neutral genetic processes. We identified reduced allelic variation in the invaded range for the less successful tubenose goby. SNPs putatively under selection were responsible for the observed differences in population structure between marker types for round goby (successful) but not tubenose goby (less successful). A higher proportion of functional loci experienced divergent selection for round goby, suggesting increased evolutionary potential in invaded ranges may be associated with round goby's greater invasion success. Genes involved in thermal tolerance were divergent for round goby populations but not tubenose goby, consistent with the hypothesis that invasion success for fish in temperate regions is influenced by capacity for thermal tolerance. Our results highlight the need to incorporate functional genetic markers in studies to better assess evolutionary potential for the improved conservation and management of species.

Genetic basis of adult migration timing in anadromous steelhead discovered through multivariate association testing

Migration traits are presumed to be complex and to involve interaction among multiple genes. We used both univariate analyses and a multivariate random forest (RF) machine learning algorithm to conduct association mapping of 15 239 single nucleotide polymorphisms (SNPs) for adult migration-timing phenotype in steelhead (Oncorhynchus mykiss). Our study focused on a model natural population of steelhead that exhibits two distinct migration-timing life histories with high levels of admixture in nature. Neutral divergence was limited between fish exhibiting summer- and winter-run migration owing to high levels of interbreeding, but a univariate mixed linear model found three SNPs from a major effect gene to be significantly associated with migration timing (p < 0.000005) that explained 46% of trait variation. Alignment to the annotated Salmo salar genome provided evidence that all three SNPs localize within a 46 kb region overlapping GREB1-like (an oestrogen target gene) on chromosome Ssa03. Additionally, multivariate analyses with RF identified that these three SNPs plus 15 additional SNPs explained up to 60% of trait variation. These candidate SNPs may provide the ability to predict adult migration timing of steelhead to facilitate conservation management of this species, and this study demonstrates the benefit of multivariate analyses for association studies.

Climate variables explain neutral and adaptive variation within salmonid metapopulations: the importance of replication in landscape genetics

Understanding how environmental variation influences population genetic structure is important for conservation management because it can reveal how human stressors influence population connectivity, genetic diversity and persistence. We used riverscape genetics modelling to assess whether climatic and habitat variables were related to neutral and adaptive patterns of genetic differentiation (population-specific and pairwise FST ) within five metapopulations (79 populations, 4583 individuals) of steelhead trout (Oncorhynchus mykiss) in the Columbia River Basin, USA. Using 151 putatively neutral and 29 candidate adaptive SNP loci, we found that climate-related variables (winter precipitation, summer maximum temperature, winter highest 5% flow events and summer mean flow) best explained neutral and adaptive patterns of genetic differentiation within metapopulations, suggesting that climatic variation likely influences both demography (neutral variation) and local adaptation (adaptive variation). However, we did not observe consistent relationships between climate variables and FST across all metapopulations, underscoring the need for replication when extrapolating results from one scale to another (e.g. basin-wide to the metapopulation scale). Sensitivity analysis (leave-one-population-out) revealed consistent relationships between climate variables and FST within three metapopulations; however, these patterns were not consistent in two metapopulations likely due to small sample sizes (N = 10). These results provide correlative evidence that climatic variation has shaped the genetic structure of steelhead populations and highlight the need for replication and sensitivity analyses in land and riverscape genetics.

Integration of Random Forest with population-based outlier analyses provides insight on the genomic basis and evolution of run timing in Chinook salmon (Oncorhynchus tshawytscha)

Anadromous Chinook salmon populations vary in the period of river entry at the initiation of adult freshwater migration, facilitating optimal arrival at natal spawning. Run timing is a polygenic trait that shows evidence of rapid parallel evolution in some lineages, signifying a key role for this phenotype in the ecological divergence between populations. Studying the genetic basis of local adaptation in quantitative traits is often impractical in wild populations. Therefore, we used a novel approach, Random Forest, to detect markers linked to run timing across 14 populations from contrasting environments in the Columbia River and Puget Sound, USA. The approach permits detection of loci of small effect on the phenotype. Divergence between populations at these loci was then examined using both principle component analysis and FST outlier analyses, to determine whether shared genetic changes resulted in similar phenotypes across different lineages. Sequencing of 9107 RAD markers in 414 individuals identified 33 predictor loci explaining 79.2% of trait variance. Discriminant analysis of principal components of the predictors revealed both shared and unique evolutionary pathways in the trait across different lineages, characterized by minor allele frequency changes. However, genome mapping of predictor loci also identified positional overlap with two genomic outlier regions, consistent with selection on loci of large effect. Therefore, the results suggest selective sweeps on few loci and minor changes in loci that were detected by this study. Use of a polygenic framework has provided initial insight into how divergence in a trait has occurred in the wild.

Geo-Referenced, Abundance Calibrated Ocean Distribution of Chinook Salmon (Oncorhynchus tshawytscha) Stocks across the West Coast of North America

Understanding seasonal migration and localized persistence of populations is critical for effective species harvest and conservation management. Pacific salmon (genus Oncorhynchus) forecasting models predict stock composition, abundance, and distribution during annual assessments of proposed fisheries impacts. Most models, however, fail to account for the influence of biophysical factors on year-to-year fluctuations in migratory distributions and stock-specific survival. In this study, the ocean distribution and relative abundance of Chinook salmon (O. tshawytscha) stocks encountered in the California Current large marine ecosystem, U.S.A were inferred using catch-per-unit effort (CPUE) fisheries and genetic stock identification data. In contrast to stock distributions estimated through coded-wire-tag recoveries (typically limited to hatchery salmon), stock-specific CPUE provides information for both wild and hatchery fish. Furthermore, in contrast to stock composition results, the stock-specific CPUE metric is independent of other stocks and is easily interpreted over multiple temporal or spatial scales. Tests for correlations between stock-specific CPUE and stock composition estimates revealed these measures diverged once proportional contributions of locally rare stocks were excluded from data sets. A novel aspect of this study was collection of data both in areas closed to commercial fisheries and during normal, open commercial fisheries. Because fishing fleet efficiency influences catch rates, we tested whether CPUE differed between closed area (non-retention) and open area (retention) data sets. A weak effect was indicated for some, but not all, analyzed cases. Novel visualizations produced from stock-specific CPUE-based ocean abundance facilitates consideration of how highly refined, spatial and genetic information could be incorporated in ocean fisheries management systems and for investigations of biogeographic factors that influence migratory distributions of fish.

Genes predict long distance migration and large body size in a migratory fish, Pacific lamprey

Elucidation of genetic mechanisms underpinning migratory behavior could help predict how changes in genetic diversity may affect future spatiotemporal distribution of a migratory species. This ability would benefit conservation of one such declining species, anadromous Pacific lamprey (Entosphenus tridentatus). Nonphilopatric migration of adult Pacific lamprey has homogenized population-level neutral variation but has maintained adaptive variation that differentiates groups based on geography, run-timing and adult body form. To investigate causes for this adaptive divergence, we examined 647 adult lamprey sampled at a fixed location on the Columbia River and radiotracked during their subsequent upstream migration. We tested whether genetic variation [94 neutral and adaptive single nucleotide polymorphisms (SNPs) previously identified from a genomewide association study] was associated with phenotypes of migration distance, migration timing, or morphology. Three adaptive markers were strongly associated with morphology, and one marker also correlated with upstream migration distance and timing. Genes physically linked with these markers plausibly influence differences in body size, which is also consistently associated with migration distance in Pacific lamprey. Pacific lamprey conservation implications include the potential to predict an individual's upstream destination based on its genotype. More broadly, the results suggest a genetic basis for intrapopulation variation in migration distance in migratory species.

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.

Clock gene evolution: seasonal timing, phylogenetic signal, or functional constraint?

Genetic determinants of seasonal reproduction are not fully understood but may be important predictors of organism responses to climate change. We used a comparative approach to study the evolution of seasonal timing within a fish community in a natural common garden setting. We tested the hypothesis that allelic length variation in the PolyQ domain of a circadian rhythm gene, Clock1a, corresponded to interspecific differences in seasonal reproductive timing across 5 native and 1 introduced cyprinid fishes (n = 425 individuals) that co-occur in the Rio Grande, NM, USA. Most common allele lengths were longer in native species that initiated reproduction earlier (Spearman's r = -0.70, P = 0.23). Clock1a allele length exhibited strong phylogenetic signal and earlier spawners were evolutionarily derived. Aside from length variation in Clock1a, all other amino acids were identical across native species, suggesting functional constraint over evolutionary time. Interestingly, the endangered Rio Grande silvery minnow (Hybognathus amarus) exhibited less allelic variation in Clock1a and observed heterozygosity was 2- to 6-fold lower than the 5 other (nonimperiled) species. Reduced genetic variation in this functionally important gene may impede this species' capacity to respond to ongoing environmental change.