Postglacial genetic differentiation of reproductive ecotypes of kokanee Oncorhynchus nerka in Okanagan Lake, British Columbia

Okanagan Lake, south-central interior of BC, contains two reproductive ecotypes of kokanee Oncorhynchus nerka ; individuals spawn in tributary streams ('stream-spawners') as well as on shoreline gravel areas ('beach-spawners'). We tested the hypothesis that these sympatric ecotypes comprise a single panmictic population by assaying variation in morphological traits and at allozyme, mitochondrial and minisatellite DNA loci in fish collected from three stream-spawning and two beach-spawning sites. No morphological traits consistently distinguished the reproductive ecotypes with the exception of the number of anal fin rays which was greater in stream-spawning kokanee. Four of 18 allozyme loci screened were polymorphic, but no significant allele frequency differences were detected among populations within ecotypes or between ecotypes. Similarly, allele frequencies at two minisatellite DNA loci were not significantly different among populations or between ecotypes. By contrast, significant differences in the frequencies of mitochondrial DNA restriction fragment length polymorphism (mtDNA RFLP) haplotypes were detected between stream- and beach-spawners, but not among populations within ecotypes. Further, two RFLPs that distinguished stream- and beach-spawning adults were found in juvenile kokanee sampled from the limnetic zone of Okanagan Lake. The two mtDNA RFLPs and a d-loop sequence variant appear to be unique to Okanagan Lake kokanee because we did not observe these haplotypes in sockeye salmon and kokanee sampled outside of Okanagan Lake. Our data suggest that: (i) there is restricted female-mediated gene flow between stream- and beach-spawning kokanee in Okanagan Lake, (ii) the forms have diverged within the lake basin since the retreat of the Wisconsinian glaciers (< ≊ 11 000 years ago), and (iii) distinct reproductive niches may promote divergence in north temperate freshwater fish faunas.     

Genetic evidence for ecological divergence in kokanee salmon

The evolution of locally adapted phenotypes among populations that experience divergent selective pressures is a central mechanism for generating and maintaining biodiversity. Recently, the advent of high‐throughput DNA sequencing technology has provided tools for investigating the genetic basis of this process in natural populations of nonmodel organisms. Kokanee, the freshwater form of sockeye salmon (Oncorhynchus nerka), occurs as two reproductive ecotypes, which differ in spawning habitat (tributaries vs. shorelines); however, outside of the spawning season the two ecotypes co‐occur in many lakes and lack diagnostic morphological characteristics. We used restriction site‐associated DNA (RAD) sequencing to identify 6145 SNPs and genotype kokanee from multiple spawning sites in Okanagan Lake (British Columbia, Canada). Outlier tests revealed 18 loci putatively under divergent selection between ecotypes, all of which exhibited temporally stable allele frequencies within ecotypes. Six outliers were annotated to sequences in the NCBI database, two of which matched genes associated with early development. There was no evidence for neutral genetic differentiation; however, outlier loci demonstrated significant structure with respect to ecotype and had high assignment accuracy in mixed composition simulations. The absence of neutral structure combined with a small number of highly divergent outlier loci is consistent with theoretical predictions for the early stages of ecological divergence. These outlier loci were then applied to a realistic fisheries scenario in which additional RAD sequencing was used to genotype kokanee collected by trawl in Okanagan Lake, providing preliminary evidence that this approach may be an effective tool for conservation and management.     

Two Reproductive Life History Types of Kokanee, Onchorynchus nerka, Exhibit Multivariate Morphometric and Protein Genetic Differentiation

Two reproductive types of kokanee are found in Okanagan Lake, British Columbia: one form that spawns in streams, and another that spawns approximately 2–4 weeks later along beaches of the lake. We examined the levels and patterns of genetic and morphometric variation among three populations (1 beach and 2 stream populations) to better understand life history differentiation. We assayed allozyme variation at 74 loci and identified 15 P_0.095 loci. Average F_ST was 0.041 among the three samples; the two stream-spawning populations grouped together in genetic distance analyses. We examined multivariate morphometric variation and fin size/shape variation using 35 truss characters. Populations did not sort morphologically by spawning type in principal component and relative warp analyses. Instead both analyses indicated that one stream-spawning population (Peachland Creek) was significantly more robust and had shallower caudal fins in comparison to the other samples. Second, clear multivariate shape differences between males and females were identified involving jaw size, mid-body dimensions, and caudal fin shape. A simple bivariate plot of tail `forkness' against fork length demonstrated that males had more forked tails and females had flattened tails. Level of differentiation of secondary sexual characteristics did not vary between the spawning types, although the beach-spawning population and one stream-spawning population did not show the strong, gender-discriminating variability in tail forkness. Although these two reproductive life history types of kokanee form discrete genetic groups, environmental differences are apparently insufficient to effect consistent differences in body shape, or fin size and shape between stream and beach-spawning morphs.     

Variability of kokanee and rainbow trout food habits,distribution, and population dynamics,in an ultraoligotrophic lake with no manipulative management

Crater Lake is a unique environment to evaluate the ecology of introduced kokanee and rainbow trout because of its otherwise pristine state, low productivity, absence of manipulative management, and lack of lotic systems for fish spawning. Between 1986 and 2004, kokanee displayed a great deal of variation in population demographics with a pattern that reoccurred in about 10 years. We believe that the reoccurring pattern resulted from density dependent growth, and associated changes in reproduction and abundance, driven by prey resource limitation that resulted from low lake productivity exacerbated by prey consumption when kokanee were abundant. Kokanee fed primarily on small-bodied prey from the mid-water column; whereas rainbow trout fed on large-bodied prey from the benthos and lake surface. Cladoceran zooplankton abundance may be regulated by kokanee. And kokanee growth and reproductive success may be influenced by the availability of Daphnia pulicaria, which was absent in zooplankton samples collected annually from 1990 to 1995, and after 1999. Distribution and diel migration of kokanee varied over the duration of the study and appeared to be most closely associated with prey availability, maximization of bioenergetic efficiency, and fish density. Rainbow trout were less abundant than were kokanee and exhibited less variation in population demographics, distribution, and food habits. There is some evidence that the population dynamics of rainbow trout were in-part related to the availability of kokanee as prey.     

MOLECULAR GENETIC EVIDENCE FOR PARALLEL LIFE-HISTORY EVOLUTION WITHIN A PACIFIC SALMON (SOCKEYE SALMON AND KOKANEE,ONCORHYNCHUS NERKA)

The Pacific salmon Oncorhynchus nerka typically occurs as a sea-run form (sockeye salmon) or may reside permanently in lakes (kokanee) thoughout its native North Pacific. We tested whether such geographically extensive ecotypic variation resulted from parallel evolutionary divergence thoughout the North Pacific or whether the two forms are monophyletic groups by examining allelic variation between sockeye salmon and kokanee at two minisatellite DNA repeat loci and in mitochondrial DNA (mtDNA) Bgl II restriction sites. Our examination of over 750 fish from 24 populations, ranging from Kamchatka to the Columbia River, identified two major genetic groups of North Pacific O. nerka: a “northwestern” group consisting of fish from Kamchatka, western Alaska, and northwestern British Columbia, and a “southern” group consisting of sockeye salmon and kokanee populations from the Fraser and Columbia River systems. Maximum-likelihood analysis accompanied by bootstrapping provided strong support for these two genetic groups of O. nerka; the populations did not cluster by migratory form, but genetic affinities were organized more strongly by geographic proximity. The two major genetic groups resolved in our study probably stem from historical isolation and dispersal of O. nerka from two major Wisconsinan glacial refugia in the North Pacific. There were significant minisatellite DNA allele frequency differences between sockeye salmon and kokanee populations from different parts of the same watershed, between populations spawning in different tributaries of the same lake, and also between sympatric populations spawning in the same stream at the same time. MtDNA Bgl II restriction site variation was significant between sockeye salmon and kokanee spawning in different parts of the same major watershed but not between forms spawning in closer degrees of reproductive sympatry. Patterns of genetic affinity and allele sharing suggested that kokanee have arisen from sea-run sockeye salmon several times independently in the North Pacific. We conclude that sockeye salmon and kokanee are para- and polyphyletic, respectively, and that the present geographic distribution of the ecotypes results from parallel evolutionary origins of kokanee from sockeye (divergences between them) thoughout the North Pacific.     

Results of the introduction of the resident form of nerka <Emphasis Type="Italic">Oncorhynchus nerka</Emphasis> in Tolmachevskoe Lake (Kamchatka)

Data on the present-day (2003–2007) state of the population of kokanee Oncorhynchus nerka introduced in Tolmachevskoe Lake more than 20 years ago are provided. Since the moment of introduction, the population of kokanee has underwent considerable changes: a decrease occurred in the specific rate of growth, average length and weight of fish in catches, and fecundity. The pattern of feeding changed from planktonic to benthic. Kokanee began to feed during spawning. Spawning shifted to later dates (from the beginning of September to the end of September and beginning of October).     

Microsatellite DNA Data Indicate Distinct Native Populations of Kokanee, Oncorhynchus nerka, Persist in the Lake Sammamish Basin, Washington

Large-scale introductions of resident and anadromous salmonids from exogenous sources and urbanization have led to major changes in, and concern for the fate of, indigenous fish populations of the Lake Sammamish/Lake Washington Basin. Specifically, introductions of kokanee (the resident form of Oncorhynchus nerka) from the Lake Whatcom Hatchery and sockeye (the anadromous form of O. nerka) from Baker Lake have caused uncertainty about the ancestry of the kokanee that currently spawn in the basin. We used nine microsatellite loci to investigate the inter-relationships of kokanee populations that spawn in streams in the Sammamish sub-basin, sockeye salmon populations that share spawning areas with the kokanee, Lake Whatcom Hatchery kokanee and Baker Lake sockeye, and an outgroup, Meadow Creek kokanee, from Lake Kootenay which drains into the upper Columbia River. We observed high levels of genetic variation (5–49 alleles per locus). Explicit tests of population sub-division revealed that collections from most spawning aggregations differed from each other. Observed allele frequency distributions strongly suggest that natural spawning kokanee in the basin are not descended from recent Lake Whatcom stock introductions. We found no compelling evidence to suggest that the kokanee sampled from spawning areas within the Lake Sammamish sub-basin have resulted from, or been altered substantially by, past introductions of non-native kokanee or sockeye.     

Combining the analyses of introgressive hybridisation and linkage mapping to investigate the genetic architecture of population divergence in the lake whitefish (Coregonus clupeaformis, Mitchill)

Adaptation and reproductive isolation, the engines of biological diversity, are still elusive when discussing the genetic bases of speciation. Namely, the number of genes and magnitude of selection acting positively or negatively on genomic traits implicated in speciation is contentious. Here, we describe the first steps of an ongoing research program aimed at understanding the genetic bases of population divergence and reproductive isolation in the lake whitefish (Coregonus clupeaformis). A preliminary linkage map originating from a hybrid cross between dwarf and normal ecotypes is presented, whereby some of the segregating AFLP markers were found to be conserved among natural populations. Maximum-likelihood was used to estimate hybrid indices from non-diagnostic markers at 998 AFLP loci. This allowed identification of the most likely candidate loci that have been under the influence of selection during the natural hybridisation of whitefish originating from different glacial races. As some of these loci could be identified on the linkage map, the possibility that selection of traits in natural populations may eventually be correlated to specific chromosomal regions was demonstrated. The future prospects and potential of these approaches to elucidate the genetic bases of adaptation and reproductive isolation among sympatric ecotypes of lake whitefish is discussed.     

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.     

Investigating the extent of parallelism in morphological and genomic divergence among lake trout ecotypes in Lake Superior

Understanding the emergence of species through the process of ecological speciation is a central question in evolutionary biology which also has implications for conservation and management. Lake trout (Salvelinus namaycush) is renowned for the occurrence of different ecotypes linked to resource and habitat use throughout North America. We aimed to unravel the fine genetic structure of the four lake trout ecotypes in Lake Superior. A total of 486 individuals from four sites were genotyped at 6822 filtered SNPs using RADseq technology. Our results revealed different extent of morphological and genetic differentiation within the different sites. Overall, genetic differentiation was weak but significant and was on average three times higher between sites (mean F_ST = 0.016) than between ecotypes within sites (mean F_ST = 0.005) indicating higher level of gene flow or a more recent shared ancestor between ecotypes within each site than between populations of the same ecotype. Evidence of divergent selection was also found between ecotypes and/or in association with morphological variation. Outlier loci found in genes related to lipid metabolism and visual acuity were of particular interest in this context of ecotypic divergence. However, we did not find clear indication of parallelism at the genomic level, despite the presence of phenotypic parallelism among some ecotypes from different sampling sites. Overall, the occurrence of different levels of both genomic and phenotypic differentiation between ecotypes within each site with several differentiated loci linked to relevant biological functions supports the presence of a continuum of divergence in lake trout.     

Genetic diversity of sockeye salmon <Emphasis Type="Italic">Oncorhynchus nerka</Emphasis> Walbaum of Kamchatka and the Commander Islands based on analysis of the variability of microsatellite DNA

The genetic diversity of the resident and migratory forms of sockeye salmon is investigated in 14 populations from various water bodies of Kamchatka and the Commander Islands by ten loci of microsatellite DNA. There are considerable differences in the frequencies of alleles among the populations of kokanee from Lake Kronotskoe, the residual form of sockeye salmon from Lake Kopylie, and other populations analyzed. Clustering of samples corresponds to their geographic position. No differences in the frequencies of alleles of the investigated loci are found between two forms of resident sockeye salmon from Kronotskoe Lake. In the sockeye salmon from the Commander Islands, a relatively low genetic diversity is found, as well as the greatest remoteness from the other Kamchatka group.     

Age structure and growth of resident kokanee <Emphasis Type="Italic">Oncorhynchus nerka</Emphasis> of natural and introduced populations in lakes of Kamchatka

The analysis of growth of kokanee Oncorhynchus nerka demonstrated that there are differences in several indices between the natural population of Kronotskoe Lake and artificial populations of Tolmachevskoe and Klyuchevoe lakes. The highest specific rate of growth was recorded for individuals in Kronotskoe Lake; the lowest, in Tolmachevskoe and Klyuchevoe lakes (in 2006). For the Tolmachev kokanee, a multiyear tendency for a decrease in the gain of scales in generations was found, which is determined by the persisting formation of the population and the depletion of food resources. A considerable effect on the process of formation of the Tolmachev population of kokanee is made by commercial activity and the related transformation of the lake into a reservoir.     

Parasite diversity,patterns of MHC II variation and olfactory based mate choice in diverging three-spined stickleback ecotypes

Ecological speciation has been the subject of intense research in evolutionary biology but the genetic basis of the actual mechanism driving reproductive isolation has rarely been identified. The extreme polymorphism of the major histocompatibility complex (MHC), probably maintained by parasite-mediated selection, has been proposed as a potential driver of population divergence. We performed an integrative field and experimental study using three-spined stickleback river and lake ecotypes. We characterized their parasite load and variation at MHC class II loci. Fish from lakes and rivers harbor contrasting parasite communities and populations possess different MHC allele pools that could be the result of a combined action of genetic drift and parasite-mediated selection. We show that individual MHC class II diversity varies among populations and is lower in river ecotypes. Our results suggest the action of homogenizing selection within habitat type and diverging selection between habitat types. Finally, reproductive isolation was suggested by experimental evidence: in a flow channel design females preferred assortatively the odor of their sympatric male. This demonstrates the role of olfactory cues in maintaining reproductive isolation between diverging fish ecotypes.     

Long- term effects of previous experience determine nutrient discrimination abilities in birds

Background

The evolution of reproductive traits, such as hybrid incompatibility (postzygotic isolation) and species recognition (prezygotic isolation), have shown their key role in speciation. Theoretical modeling has recently predicted that close linkage between genes controlling pre- and postzygotic reproductive isolation could accelerate the conditions for speciation. Postzygotic isolation could develop during the sympatric speciation process contributing to the divergence of populations. Using hybrid fitness as a measure of postzygotic reproductive isolation, we empirically studied population divergence in perch (Perca fluviatilis L.) from two genetically divergent populations within a lake.

Results

During spawning time of perch we artificially created parental offspring and F₁ hybrids of the two populations and studied fertilization rate and hatching success under laboratory conditions. The combined fitness measure (product of fertilization rate and hatching success) of F₁ hybrids was significantly reduced compared to offspring from within population crosses.

Conclusion

Our results suggest intrinsic genetic incompatibility between the two populations and indicate that population divergence between two populations of perch inhabiting the same lake may indeed be promoted by postzygotic isolation.     

Influences of Lake Trout (Salvelinus namaycush) and Mysis diluviana on kokanee (Oncorhynchus nerka) in Lake Pend Oreille,Idaho

Research on Lake Pend Oreille, Idaho, has focused on the influence of two potential limiting factors for kokanee Oncorhynchus nerka (Walbaum, 1792): competition for food with Mysis diluviana (Loven, 1862, hereafter Mysis) and predation by lake trout Salvelinus namaycush (Walbaum, 1792). Population fluctuations of Mysis and lake trout have resulted in substantial heterogeneity in food web conditions, apparently altering both bottom-up and top-down dynamics. Therefore, relative importance of predation and competition were evaluated as drivers of kokanee abundance, biomass, and production. A series of general linear models was used to evaluate relative influences of Mysis and lake trout on kokanee. Kokanee production was a density-dependent process and the collapse of Mysis corresponded to an increase in the modeled maximum annual production of kokanee from 224 tonnes to 408 tonnes. Lake trout also negatively influenced kokanee biomass. A Mysis-mediated, predator-induced kokanee biomass collapse occurred when lake trout and Mysis abundances were both high. Sustainable management of this fishery requires recognition that competition with Mysis will define the scope of kokanee production and therefore the scope of sustainable predation.

     

A combination of molecular markers identifies evolutionarily significant units in Cedrela odorata L. (Meliaceae) in Costa Rica

The necessity for conservation of the geneticcomponent of biodiversity is now widelyrecognised. A broad genetic base is required tomaintain evolutionary potential and thepopulation erosion occurring in much of theworld's forests threatens the genetic integrityof many tree species. Spanish Cedar (Cedrela odorata L.) has been under severepressure for generations and is now the focusof a study aimed at assessing the levels anddistribution of genetic diversity in remainingpopulations. Ten Costa Rican populations wereanalysed using chloroplast and AFLP markers.The overall level of diversity was as expectedfor an outcrossing, long-lived, woody species(H _T = 0.27). However, this concealeda deep divergence within the species, forchloroplast and AFLP (_CT = 0.83)markers. Populations were differentiated in twogroups that exhibited contrasting habitatpreferences and two ecotypes, wet and dry, wereidentified. Within the ecotypes, all but onepopulation were fixed for a single chloroplasthaplotype and within populations, total genomicdiversity levels were low (H _S= 0.03–0.13). Populations possessing the dryecotype maintained significantly more diversitythan those from wet regions. Within the wetecotype group, pairwise genetic distancebetween populations fitted an isolation bydistance model. The group was stronglysubdivided and showed isolation by distancearound the southern edge of the centralmountain ranges. The genetic divergence of thetwo ecotypes, observed at both organellar andnuclear loci, identifies evolutionarilysignificant units that, taken together withprevious studies of the species, provide arational basis on which to build a conservationpolicy for the species.     

Low genetic diversity and significant population structuring in the relict Amentotaxus argotaenia complex (Taxaceae) based on ISSR fingerprinting

Amentotaxus, a genus of the Taxaceae, represents an ancient lineage that has long existed in Eurasia. All Amentotaxus species experienced frequent population expansion and contraction over periodical glaciations in Tertiary and Quaternary. Among them, Amentotaxus argotaenia complex consists of three morphologically alike species, A. argotaenia, Amentotaxus yunnanensis, and Amentotaxus formosana. This complex is distributed in the subtropical region of mainland China and Taiwan where many Pleistocene refugia have been documented. In this study, genetic diversity and population structuring within and between species were investigated based on the inter-simple sequence repeats (ISSR) fingerprinting. Mean genetic diversity within populations was estimated in three ways: (1) the percentage of polymorphic loci out of all loci (P) (2) Neis unbiased expected heterozygosity (He), and (3) Shannons index of phenotypic diversity. For a total of 310 individuals of 15 populations sampled from the three species, low levels of ISSR genetic variation within populations were detected, with P=4.66–16.58%, He=0.0176–0.0645 and Hpop=0.0263–0.0939, agreeing with their seriously threatened status. AMOVA analyses revealed that the differences between species only accounted for 27.38% of the total variation, whereas differences among populations and within populations were 57.70 and 14.92%, respectively, indicating substantial isolation between the patch-like populations. A neighbor-joining tree identified a close affinity between A. yunnanensis and A. formosana. Genetic drift due to small population size, plus limited current gene flow, resulted in significant genetic structuring. Low levels of intrapopulational genetic variation and considerable interpopulational divergence were also attributable to demographic bottlenecks during and/or after the Pleistocene glaciations.     

Genetic diversity and structure of the endemic Caesalpinia hintonii complex (Leguminosae: Caesalpinioideae) in Mexico

The Caesalpinia hintonii complex is formed by five endemic species (C. hintonii, C. laxa, C. macvaughii, C. melanadenia and C. epifanioi) occurring in central Mexico. This species complex is under incipient genetic divergence as by-product of local adaptations in reproductive and morphological traits to different habitats. We estimate the genetic variation and structure of populations of this species complex to assess the extent of genetic differentiation among populations and related species along its geographic distribution. Estimations of genetic diversity and structure were done based on ten enzymes and 18 loci. Mean number of alleles per locus ranged from 1.5 to 1.9. Polymorphic loci ranged from 42.1 to 68.4. Observed (H_o: range 0.191–0.275) and expected (H_e: range 0.205–0.317) heterozygosities in this complex were higher compared with other endemic and legume species. Neis genetic diversity estimates showed that most genetic variation was found within (H_S = 0.325) rather than among populations (D_ST = 0.085). Populations of the species C. hintonii showed a considerable genetic differentiation (F_ST = 0.207). The results of genetic diversity and structure within and among populations are in accord with the great morphological differentiation described for this species complex.     

Associations between isozyme phenotypes and environment in the slender wild oat (Avena barbata) in Israel

Summary Collections from 31 populations of A. barbata from diverse habitats in Israel were assayed electrophoretically for seven enzyme systems. Phenotype frequencies were scored in nine enzyme zones, probably representing 27 loci, to determine isozyme variability within and among populations. Many different isozyme phenotypes were found in all of the populations; also the array of isozyme phenotypes found in each population differed distinctly from that found in each other population. Overlays of phenotypic frequencies on map locations showed that isozyme variability is distributed in mosaic patterns not related to geographical distance. Principal-component and multiple-regression analyses revealed that temperature and moisture-related variables are significantly correlated with particular isozyme phenotypes. Further, the mosaic patterns of isozyme variation were found to correspond closely to mosaic patterns of the habitat. This structuring of the genetic variability into multilocus combinations was attributed to the combined effects of directional and diversifying selection. Comparisons of patterns and extent of genetic variation in Israel and California led to the conclusion that the evolution of ecotypes, each adapted to a specific habitat and marked by a particular set of enzyme alleles, has proceeded further in Israel, where A. barbata is endemic, than in California, where it is a recent introduction.This study was supported in part by NSF Grant BMS-01113-A01. Seed collections were supported by a United States-Israel Binational Science Foundation Grant     

Countergradient variation and secondary sexual color: phenotypic convergence promotes genetic divergence in carotenoid use between sympatric anadromous and nonanadromous morphs of sockeye salmon (Oncorhynchus nerka)

Genetically distinct anadromous (sockeye) and nonanadromous (kokanee) morphs of the Pacific salmon, Oncorhynchus nerka, develop identical, brilliant red color at maturity during sympatric breeding in freshwater streams. The marine and lacustrine environments they occupy prior to maturity, however, appear to differ in the availability of dietary carotenoid pigments necessary to produce red coloration. We tested the hypothesis that kokanee, which occupy carotenoid-poor lakes, are more efficient at using the dietary pigments than are sockeye, which occupy the more productive North Pacific Ocean. In a 2-year controlled breeding study, flesh and skin color of mature and immature crosses fed a low-carotenoid diet were quantified with both a chromameter and by chemical extraction of carotenoid pigments. Results revealed striking countergradient variation in carotenoid use, with kokanee approximately three times more efficient at sequestering the pigments to the flesh musculature than similar age sockeye. This difference translated into virtually nonoverlapping differences between pure crosses in secondary sexual color at maturity, when the pigments are mobilized and transported to the skin. Kokanee crosses turned pinkish red over most of their body, whereas sockeye turned olive green. The olive green was similar to the breeding color of residuals in the wild, the progeny of anadromous sockeye that remain in fresh water and are believed to have given rise to kokanee on numerous independent occasions. Reciprocal hybrids were similar to each other and intermediate to the pure crosses, indicating additive genetic inheritance. Mate choice trials with sockeye males in the wild showed the ancestral morph strongly preferred red over green models. These results suggest a preference for red mates maintained in nonanadromous breeding populations drove the reevolution of the red phenotype in kokanee via more efficient use of dietary carotenoid pigments. This is a novel, yet hidden, mechanism by which sexual selection promotes the genetic differentiation of these sympatric populations.