Local and systemic replicative fitness for viruses in specialist, generalist, and non-specialist interactions with salmonid hosts.

Host tissues represent diverse resources or barriers for pathogen replicative fitness. We tested whether viruses in specialist, generalist, and non-specialist interactions replicate differently in local entry tissue (fin), and systemic target tissue (kidney) using infectious hematopoietic necrosis virus (IHNV) and three salmonid fish hosts. Virus tissue replication was host specific, but one feature was shared by specialists and the generalist which was uncommon in the non-specialist interactions: high host entry and replication capacity in the local tissue after contact. Moreover, specialists showed increased replication in systemic target tissues early after host contact. By comparing ancestral and derived IHNV viruses, we also characterized replication tradeoffs associated with specialist and generalist evolution. Compared with the ancestral virus, a derived specialist gained early local replicative fitness in the new host but lost replicative fitness in the ancestral host. By contrast, a derived generalist showed small replication losses relative to the ancestral virus in the ancestral host but increased early replication in the local tissue of novel hosts. This study shows that the mechanisms of specialism and generalism are host specific and that local and systemic replication can contribute differently to overall within host replicative fitness for specialist and generalist viruses.

Correction: Demographic history shaped geographical patterns of deleterious mutation load in a broadly distributed Pacific Salmon.

[This corrects the article DOI: 10.1371/journal.pgen.1008348.].

Long-distance migration is a major factor driving local adaptation at continental scale in Coho salmon.

Inferring the genomic basis of local adaptation is a long-standing goal of evolutionary biology. Beyond its fundamental evolutionary implications, such knowledge can guide conservation decisions for populations of conservation and management concern. Here, we investigated the genomic basis of local adaptation in the Coho salmon (Oncorhynchus kisutch) across its entire North American range. We hypothesized that extensive spatial variation in environmental conditions and the species' homing behaviour may promote the establishment of local adaptation. We genotyped 7829 individuals representing 217 sampling locations at more than 100,000 high-quality RADseq loci to investigate how recombination might affect the detection of loci putatively under selection and took advantage of the precise description of the demographic history of the species from our previous work to draw accurate population genomic inferences about local adaptation. The results indicated that genetic differentiation scans and genetic-environment association analyses were both significantly affected by variation in recombination rate as low recombination regions displayed an increased number of outliers. By taking these confounding factors into consideration, we revealed that migration distance was the primary selective factor driving local adaptation and partial parallel divergence among distant populations. Moreover, we identified several candidate single nucleotide polymorphisms associated with long-distance migration and altitude including a gene known to be involved in adaptation to altitude in other species. The evolutionary implications of our findings are discussed along with conservation applications.

Demographic history shaped geographical patterns of deleterious mutation load in a broadly distributed Pacific Salmon.

A thorough reconstruction of historical processes is essential for a comprehensive understanding of the mechanisms shaping patterns of genetic diversity. Indeed, past and current conditions influencing effective population size have important evolutionary implications for the efficacy of selection, increased accumulation of deleterious mutations, and loss of adaptive potential. Here, we gather extensive genome-wide data that represent the extant diversity of the Coho salmon (Oncorhynchus kisutch) to address two objectives. We demonstrate that a single glacial refugium is the source of most of the present-day genetic diversity, with detectable inputs from a putative secondary micro-refugium. We found statistical support for a scenario whereby ancestral populations located south of the ice sheets expanded recently, swamping out most of the diversity from other putative micro-refugia. Demographic inferences revealed that genetic diversity was also affected by linked selection in large parts of the genome. Moreover, we demonstrate that the recent demographic history of this species generated regional differences in the load of deleterious mutations among populations, a finding that mirrors recent results from human populations and provides increased support for models of expansion load. We propose that insights from these historical inferences should be better integrated in conservation planning of wild organisms, which currently focuses largely on neutral genetic diversity and local adaptation, with the role of potentially maladaptive variation being generally ignored.

Heterogeneous genetic basis of age at maturity in salmonid fishes.

Understanding the genetic basis of repeated evolution of the same phenotype across taxa is a fundamental aim in evolutionary biology and has applications in conservation and management. However, the extent to which interspecific life-history trait polymorphisms share evolutionary pathways remains underexplored. Here, we address this gap by studying the genetic basis of a key life-history trait, age at maturity, in four species of Pacific salmonids (genus Oncorhynchus) that exhibit intra- and interspecific variation in this trait-Chinook Salmon, Coho Salmon, Sockeye Salmon, and Steelhead Trout. We tested for associations in all four species between age at maturity and two genome regions, six6 and vgll3, that are strongly associated with the same trait in Atlantic Salmon (Salmo salar). We also conducted a genome-wide association analysis in Steelhead to assess whether additional regions were associated with this trait. We found the genetic basis of age at maturity to be heterogeneous across salmonid species. Significant associations between six6 and age at maturity were observed in two of the four species, Sockeye and Steelhead, with the association in Steelhead being particularly strong in both sexes (p = 4.46 × 10-9 after adjusting for genomic inflation). However, no significant associations were detected between age at maturity and the vgll3 genome region in any of the species, despite its strong association with the same trait in Atlantic Salmon. We discuss possible explanations for the heterogeneous nature of the genetic architecture of this key life-history trait, as well as the implications of our findings for conservation and management.

Genomic and phenotypic effects of inbreeding across two different hatchery management regimes in Chinook salmon.

Genomic approaches permit direct estimation of inbreeding and its effect on fitness. We used genomic-based estimates of inbreeding to investigate their relationship with eight adult traits in a captive-reared Pacific salmonid that is released into the wild. Estimates were also used to determine whether alternative broodstock management approaches reduced risks of inbreeding. Specifically, 1,100 unlinked restriction-site associated (RAD) loci were used to compare pairwise relatedness, derived from a relationship matrix, and individual inbreeding, estimated by comparing observed and expected homozygosity, across four generations in two hatchery lines of Chinook salmon that were derived from the same source. The lines are managed as "integrated" with the founding wild stock, with ongoing gene flow, and as "segregated" with no gene flow. While relatedness and inbreeding increased in the first generation of both lines, possibly due to population subdivision caused by hatchery initiation, the integrated line had significantly lower levels in some subsequent generations (relatedness: F2 -F4 ; inbreeding F2 ). Generally, inbreeding was similar between the lines despite large differences in effective numbers of breeders. Inbreeding did not affect fecundity, reproductive effort, return timing, fork length, weight, condition factor, and daily growth coefficient. However, it delayed spawn timing by 1.75 days per one standard deviation increase in F (~0.16). The results indicate that integrated management may reduce inbreeding but also suggest that it is relatively low in a small, segregated hatchery population that maximized number of breeders. Our findings demonstrate the utility of genomics to monitor inbreeding under alternative management strategies in captive breeding programs.

Functional Annotation of All Salmonid Genomes (FAASG): an international initiative supporting future salmonid research, conservation and aquaculture.

We describe an emerging initiative - the 'Functional Annotation of All Salmonid Genomes' (FAASG), which will leverage the extensive trait diversity that has evolved since a whole genome duplication event in the salmonid ancestor, to develop an integrative understanding of the functional genomic basis of phenotypic variation. The outcomes of FAASG will have diverse applications, ranging from improved understanding of genome evolution, to improving the efficiency and sustainability of aquaculture production, supporting the future of fundamental and applied research in an iconic fish lineage of major societal importance.

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.

Genetic variation underlying resistance to infectious hematopoietic necrosis virus in a steelhead trout (Oncorhynchus mykiss) population.

Understanding the mechanisms of host resistance to pathogens will allow insights into the response of wild populations to the emergence of new pathogens. Infectious hematopoietic necrosis virus (IHNV) is endemic to the Pacific Northwest and infectious to Pacific salmon and trout (Oncorhynchus spp.). Emergence of the M genogroup of IHNV in steelhead trout O. mykiss in the coastal streams of Washington State, between 2007 and 2011, was geographically heterogeneous. Differences in host resistance due to genetic change were hypothesized to be a factor influencing the IHNV emergence patterns. For example, juvenile steelhead trout losses at the Quinault National Fish Hatchery (QNFH) were much lower than those at a nearby facility that cultures a stock originally derived from the same source population. Using a classical quantitative genetic approach, we determined the potential for the QNFH steelhead trout population to respond to selection caused by the pathogen, by estimating the heritability for 2 traits indicative of IHNV resistance, mortality (h² = 0.377 (0.226 - 0.550)) and days to death (h² = 0.093 (0.018 - 0.203)). These results confirm that there is a genetic basis for resistance and that this population has the potential to adapt to IHNV. Additionally, genetic correlation between days to death and fish length suggests a correlated response in these traits to selection. Reduction of genetic variation, as well as the presence or absence of resistant alleles, could affect the ability of populations to adapt to the pathogen. Identification of the genetic basis for IHNV resistance could allow the assessment of the susceptibility of other steelhead populations.

Genetic evidence for landscape effects on dispersal in the army ant Eciton burchellii.

Inhibited dispersal, leading to reduced gene flow, threatens populations with inbreeding depression and local extinction. Fragmentation may be especially detrimental to social insects because inhibited gene flow has important consequences for cooperation and competition within and among colonies. Army ants have winged males and permanently wingless queens; these traits imply male-biased dispersal. However, army ant colonies are obligately nomadic and have the potential to traverse landscapes. Eciton burchellii, the most regularly nomadic army ant, is a forest interior species: colony raiding activities are limited in the absence of forest cover. To examine whether nomadism and landscape (forest clearing and elevation) affect population genetic structure in a montane E. burchellii population, we reconstructed queen and male genotypes from 25 colonies at seven polymorphic microsatellite loci. Pairwise genetic distances among individuals were compared to pairwise geographical and resistance distances using regressions with permutations, partial Mantel tests and random forests analyses. Although there was no significant spatial genetic structure in queens or males in montane forest, dispersal may be male-biased. We found significant isolation by landscape resistance for queens based on land cover (forest clearing), but not on elevation. Summed colony emigrations over the lifetime of the queen may contribute to gene flow in this species and forest clearing impedes these movements and subsequent gene dispersal. Further forest cover removal may increasingly inhibit Eciton burchellii colony dispersal. We recommend maintaining habitat connectivity in tropical forests to promote population persistence for this keystone species.

Genome-wide association analysis of the resistance to infectious hematopoietic necrosis virus in two rainbow trout aquaculture lines confirms oligogenic architecture with several moderate effect quantitative trait loci.

Infectious hematopoietic necrosis (IHN) is a disease of salmonid fish that is caused by the IHN virus (IHNV), which can cause substantial mortality and economic losses in rainbow trout aquaculture and fisheries enhancement hatchery programs. In a previous study on a commercial rainbow trout breeding line that has undergone selection, we found that genetic resistance to IHNV is controlled by the oligogenic inheritance of several moderate and many small effect quantitative trait loci (QTL). Here we used genome wide association analyses in two different commercial aquaculture lines that were naïve to previous exposure to IHNV to determine whether QTL were shared across lines, and to investigate whether there were major effect loci that were still segregating in the naïve lines. A total of 1,859 and 1,768 offspring from two commercial aquaculture strains were phenotyped for resistance to IHNV and genotyped with the rainbow trout Axiom 57K SNP array. Moderate heritability values (0.15-0.25) were estimated. Two statistical methods were used for genome wide association analyses in the two populations. No major QTL were detected despite the naïve status of the two lines. Further, our analyses confirmed an oligogenic architecture for genetic resistance to IHNV in rainbow trout. Overall, 17 QTL with notable effect (≥1.9% of the additive genetic variance) were detected in at least one of the two rainbow trout lines with at least one of the two statistical methods. Five of those QTL were mapped to overlapping or adjacent chromosomal regions in both lines, suggesting that some loci may be shared across commercial lines. Although some of the loci detected in this GWAS merit further investigation to better understand the biological basis of IHNV disease resistance across populations, the overall genetic architecture of IHNV resistance in the two rainbow trout lines suggests that genomic selection may be a more effective strategy for genetic improvement in this trait.

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.

Assortative mating for reproductive timing affects population recruitment and resilience in a quantitative genetic model.

Quantitative models that simulate the inheritance and evolution of fitness-linked traits offer a method for predicting how environmental or anthropogenic perturbations can affect the dynamics of wild populations. Random mating between individuals within populations is a key assumption of many such models used in conservation and management to predict the impacts of proposed management or conservation actions. However, recent evidence suggests that non-random mating may be underestimated in wild populations and play an important role in diversity-stability relationships. Here we introduce a novel individual-based quantitative genetic model that incorporates assortative mating for reproductive timing, a defining attribute of many aggregate breeding species. We demonstrate the utility of this framework by simulating a generalized salmonid lifecycle, varying input parameters, and comparing model outputs to theoretical expectations for several eco-evolutionary, population dynamic scenarios. Simulations with assortative mating systems resulted in more resilient and productive populations than those that were randomly mating. In accordance with established ecological and evolutionary theory, we also found that decreasing the magnitude of trait correlations, environmental variability, and strength of selection each had a positive effect on population growth. Our model is constructed in a modular framework so that future components can be easily added to address pressing issues such as the effects of supportive breeding, variable age structure, differential selection by sex or age, and fishery interactions on population growth and resilience. With code published in a public Github repository, model outputs may easily be tailored to specific study systems by parameterizing with empirically generated values from long-term ecological monitoring programs.

Population-size history inferences from the coho salmon (Oncorhynchus kisutch) genome.

Coho salmon (Oncorhynchus kisutch) are a culturally and economically important species that return from multiyear ocean migrations to spawn in rivers that flow to the Northern Pacific Ocean. Southern stocks of coho salmon in Canada and the United States have significantly declined over the past quarter century, and unfortunately, conservation efforts have not reversed this trend. To assist in stock management and conservation efforts, we generated a chromosome-level genome assembly. We also resequenced the genomes of 83 coho salmon across the North American range to identify nucleotide variants and understand the demographic histories of these salmon by modeling effective population size from genome-wide data. From demographic history modeling, we observed reductions in effective population sizes between 3,750 and 8,000 years ago for several northern sampling sites, which may correspond to bottleneck events during recolonization after glacial retreat.

Temporal variation in selection on body length and date of return in a wild population of coho salmon, Oncorhynchus kisutch.

BACKGROUND: A number of studies have measured selection in nature to understand how populations adapt to their environment; however, the temporal dynamics of selection are rarely investigated. The aim of this study was to assess the temporal variation in selection by comparing the mode, direction and strength of selection on fitness-related traits between two cohorts of coho salmon (Oncorhynchus kisutch). Specifically, we quantified individual reproductive success and examined selection on date of return and body length in a wild population at Big Beef Creek, Washington (USA). RESULTS: Reproductive success and the mode, direction and strength of selection on date of return and body length differed between two cohorts sampled in 2006 and 2007. Adults of the first brood year had greater success over those of the second. In 2006, disruptive selection favored early and late returning individuals in 2-year-old males, and earlier returning 3-year-old males had higher fitness. No evidence of selection on date of return was detected in females. In 2007, selection on date of return was not observed in males of either age class, but stabilizing selection on date of return was observed in females. No selection on body length was detected in males of both age classes in 2006, and large size was associated with higher fitness in females. In 2007, selection favored larger size in 3-year-old males and intermediate size in females. Correlational selection between date of return and body length was observed only in 2-year-old males in 2006. CONCLUSIONS: We found evidence of selection on body length and date of return to the spawning ground, both of which are important fitness-related traits in salmonid species, but this selection varied over time. Fluctuation in the mode, direction and strength of selection between two cohorts was likely to be due to factors such as changes in precipitation, occurrence of catastrophic events (flooding), the proportion of younger- versus older-maturing males, sex ratio and densities of spawners.

Variation in within-host replication kinetics among virus genotypes provides evidence of specialist and generalist infection strategies across three salmonid host species.

Theory of the evolution of pathogen specialization suggests that a specialist pathogen gains high fitness in one host, but this comes with fitness loss in other hosts. By contrast, a generalist pathogen does not achieve high fitness in any host, but gains ecological fitness by exploiting different hosts, and has higher fitness than specialists in nonspecialized hosts. As a result, specialist pathogens are predicted to have greater variation in fitness across hosts, and generalists would have lower fitness variation across hosts. We test these hypotheses by measuring pathogen replicative fitness as within-host viral loads from the onset of infection to the beginning of virus clearance, using the rhabdovirus infectious hematopoietic necrosis virus (IHNV) in salmonid fish. Based on field prevalence and virulence studies, the IHNV subgroups UP, MD, and L are specialists, causing infection and mortality in sockeye salmon, steelhead, and Chinook salmon juveniles, respectively. The UC subgroup evolved naturally from a UP ancestor and is a generalist infecting all three host species but without causing severe disease. We show that the specialist subgroups had the highest peak and mean viral loads in the hosts in which they are specialized, and they had low viral loads in nonspecialized hosts, resulting in large variation in viral load across hosts. Viral kinetics show that the mechanisms of specialization involve the ability to both maximize early virus replication and avoid clearance at later times, with different mechanisms of specialization evident in different host-virus combinations. Additional nuances in the data included different fitness levels for nonspecialist interactions, reflecting different trade-offs for specialist viruses in other hosts. The generalist UC subgroup reached intermediate viral loads in all hosts and showed the smallest variation in fitness across hosts. The evolution of the UC generalist from an ancestral UP sockeye specialist was associated with fitness increases in steelhead and Chinook salmon, but only slight decreases in fitness in sockeye salmon, consistent with low- or no-cost generalism. Our results support major elements of the specialist-generalist theory, providing evidence of a specialist-generalist continuum in a vertebrate pathogen. These results also quantify within-host replicative fitness trade-offs resulting from the natural evolution of specialist and generalist virus lineages in multi-host ecosystems.

Virulence and Infectivity of UC, MD, and L Strains of Infectious Hematopoietic Necrosis Virus (IHNV) in Four Populations of Columbia River Basin Chinook Salmon.

Infectious Hematopoietic Necrosis Virus (IHNV) infects juvenile salmonid fish in conservation hatcheries and aquaculture facilities, and in some cases, causes lethal disease. This study assesses intra-specific variation in the IHNV susceptibility of Chinook salmon (Oncorhynchus tshawytscha) in the Columbia River Basin (CRB), in the northwestern United States. The virulence and infectivity of IHNV strains from three divergent virus genogroups are measured in four Chinook salmon populations, including spring-run and fall-run fish from the lower or upper regions of the CRB. Following controlled laboratory exposures, our results show that the positive control L strain had significantly higher virulence, and the UC and MD strains that predominate in the CRB had equivalently low virulence, consistent with field observations. By several experimental measures, there was little variation in host susceptibility to infection or disease. However, a small number of exceptions suggested that the lower CRB spring-run Chinook salmon population may be less susceptible than other populations tested. The UC and MD viruses did not differ in infectivity, indicating that the observed asymmetric field prevalence in which IHNV detected in CRB Chinook salmon is 83% UC and 17% MD is not due to the UC virus being more infectious. Overall, we report little intra-species variation in CRB Chinook salmon susceptibility to UC or MD IHNV infection or disease, and suggest that other factors may instead influence the ecology of IHNV in the CRB.

Temperature Variation and Host Immunity Regulate Viral Persistence in a Salmonid Host.

Environmental variation has important effects on host-pathogen interactions, affecting large-scale ecological processes such as the severity and frequency of epidemics. However, less is known about how the environment interacts with host immunity to modulate virus fitness within hosts. Here, we studied the interaction between host immune responses and water temperature on the long-term persistence of a model vertebrate virus, infectious hematopoietic necrosis virus (IHNV) in steelhead trout (Oncorhynchus mykiss). We first used cell culture methods to factor out strong host immune responses, allowing us to test the effect of temperature on viral replication. We found that 15 ∘C water temperature accelerated IHNV replication compared to the colder 10 and 8 ∘C temperatures. We then conducted in vivo experiments to quantify the effect of 6, 10, and 15 ∘C water temperatures on IHNV persistence over 8 months. Fish held at 15 and 10 ∘C were found to have higher prevalence of neutralizing antibodies compared to fish held at 6 ∘C. We found that IHNV persisted for a shorter time at warmer temperatures and resulted in an overall lower fish mortality compared to colder temperatures. These results support the hypothesis that temperature and host immune responses interact to modulate virus persistence within hosts. When immune responses were minimized (i.e., in vitro) virus replication was higher at warmer temperatures. However, with a full potential for host immune responses (i.e., in vivo experiments) longer virus persistence and higher long-term virulence was favored in colder temperatures. We also found that the viral RNA that persisted at later time points (179 and 270 days post-exposure) was mostly localized in the kidney and spleen tissues. These tissues are composed of hematopoietic cells that are favored targets of the virus. By partitioning the effect of temperature on host and pathogen responses, our results help to better understand environmental drivers of host-pathogen interactions within hosts, providing insights into potential host-pathogen responses to climate change.

Simulating fishery-induced evolution in chinook salmon: the role of gear, location, and genetic correlation among traits.

Adaptation to human-modified ecosystems has been implicated in changing the life history of a number of wild animal populations, potentially contributing to their collapse. Fishing may be an important evolutionary force that can change the distribution of fitness-related traits; however, the magnitude and direction of the evolutionary response may be influenced by different management strategies. Most phenotypic traits subject to human-induced selection are simultaneously influenced by the environment and by genetic variation, and many traits are genetically correlated. Here, we evaluated the evolutionary outcomes of harvest activities on mean length and age at maturity in a fish population by coupling a multivariate quantitative genetic model with a Leslie life history matrix model. Lengths-at-ages were treated as genetically correlated characters parameterized from empirical data on chinook salmon (Oncorhynchus tshawytscha) populations. Using simulations, we explored the outcomes of 100 years of harvest using gill nets, which impose disruptive selection, or longlines, which impose minimum size selection, that targeted immature individuals in the high seas or maturing individuals in terminal spawning areas. Response in mean length and age depended on selection differentials imposed by harvest (which depended in turn on fishing location, gear type, and proportion of the population harvested) and on the genetic correlations between traits. Mean length was strongly influenced by the selection differential of the most abundant age class. Large differences in response were observed between the high-seas fishery, where the most abundant age was the youngest age vulnerable to harvest, compared to the terminal area fishery, where an older age class was most abundant. We observed a substantial difference in response between gill nets and longlines in the terminal fishery only. The evolution of mean age of mature individuals was less predictable, but generally increased as length decreased and decreased as length increased. The model presented here has potential for incorporating empirical data into fisheries forecasting and therefore provides a powerful means of integrating evolutionary considerations into harvest management.

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.