Conditionally Deleterious Mutation Load Accumulates in Genomic Islands of Local Adaptation but Can Be Purged with Sufficient Genotypic Redundancy.

AbstractLocal adaptation frequently evolves in patches or environments that are connected via migration. In these cases, genomic regions that are linked to a locally adapted locus experience reduced effective migration rates. Via individual-based simulations of a two-patch system, we show that this reduced effective migration results in the accumulation of conditionally deleterious mutations, but not universally deleterious mutations, adjacent to adaptive loci. When there is redundancy in the genetic basis of local adaptation (i.e., genotypic redundancy), turnover of locally adapted polymorphisms allows conditionally deleterious mutation load to be purged. The amount of mutational load that accumulates adjacent to locally adapted loci is dependent on redundancy, recombination rate, migration rate, population size, strength of selection, and the phenotypic effect size of adaptive alleles. Our results highlight the need to be cautious when interpreting patterns of local adaptation at the level of phenotype or fitness, as the genetic basis of local adaptation can be transient, and evolution may confer a degree of maladaptation to nonlocal environments.

Pelvic spine reduction affects diet but not gill raker morphology in two polymorphic brook stickleback (Culaea inconstans) populations.

Pelvic spine polymorphism occurs in several species in the stickleback family (Gasterosteidae). Given the similar phenotypic polymorphisms in multiple stickleback species, we sought to determine the extent of parallelism in the ecological correlates of pelvic spine reduction. Based on a metabarcoding analysis of brook stickleback gut contents in two polymorphic populations, we found that significant diet differences were associated with pelvic spine reduction, but we found no clear or consistent trend supporting a tendency for benthic feeding in pelvic-reduced brook sticklebacks. These results contrast with those found in threespine sticklebacks where pelvic spine reduction is often associated with a benthic diet. Hence, we found non-parallel consequences of spine polymorphism across species. Furthermore, a difference in gill raker morphology has been frequently observed between ecomorphs with different diets in many fish species. However, we found no evidence of any difference in gill raker morphology associated with pelvic spine polymorphism in brook sticklebacks.

Identification of a candidate sex determination gene in Culaea inconstans suggests convergent recruitment of an Amh duplicate in two lineages of stickleback.

Sex chromosomes vary greatly in their age and levels of differentiation across the tree of life. This variation is largely due to the rates of sex chromosome turnover in different lineages; however, we still lack an explanation for why sex chromosomes are so conserved in some lineages (e.g. mammals, birds) but so labile in others (e.g. teleosts, amphibians). To identify general mechanisms driving transitions in sex determination systems or forces which favour their conservation, we first require empirical data on sex chromosome systems from multiple lineages. Stickleback fishes are a valuable model lineage for the study of sex chromosome evolution due to variation in sex chromosome systems between closely-related species. Here, we identify the sex chromosome and a strong candidate for the master sex determination gene in the brook stickleback, Culaea inconstans. Using whole-genome sequencing of wild-caught samples and a lab cross, we identify AmhY, a male specific duplication of the gene Amh, as the candidate master sex determination gene. AmhY resides on Chromosome 20 in C. inconstans and is likely a recent duplication, as both AmhY and the sex-linked region of Chromosome 20 show little sequence divergence. Importantly, this duplicate AmhY represents the second independent duplication and recruitment of Amh as the sex determination gene in stickleback and the eighth example known across teleosts. We discuss this convergence in the context of sex chromosome turnovers and the role that the Amh/AmhrII pathway, which is crucial for sex determination, may play in the evolution of sex chromosomes in teleosts.

Unpacking Conditional Neutrality: Genomic Signatures of Selection on Conditionally Beneficial and Conditionally Deleterious Mutations.

It is common to look for signatures of local adaptation in genomes by identifying loci with extreme levels of allele frequency divergence among populations. This approach to finding genes associated with local adaptation often assumes antagonistic pleiotropy, wherein alternative alleles are strongly favored in alternative environments. Conditional neutrality has been proposed as an alternative to antagonistic pleiotropy, but conditionally neutral polymorphisms are transient, and it is unclear how much outlier signal would be maintained under different forms of conditional neutrality. Here, we use individual-based simulations and a simple analytical heuristic to show that a pattern that mimics local adaptation at the phenotypic level, where each genotype has the highest fitness in its home environment, can be produced by the accumulation of mutations that are neutral in their home environment and deleterious in nonlocal environments. Because conditionally deleterious mutations likely arise at a rate many times higher than conditionally beneficial mutations, they can have a significant cumulative effect on fitness even when individual effect sizes are small. We show that conditionally deleterious mutations driving nonlocal maladaptation may be undetectable by even the most powerful genome scans, as differences in allele frequency between populations are typically small. We also explore the evolutionary effects of conditionally beneficial mutations and find that they can maintain significant signals of local adaptation, and they would be more readily detectable than conditionally deleterious mutations using conventional genome scan approaches. We discuss implications for interpreting outcomes of transplant experiments and genome scans that are used to study the genetic basis of local adaptation.

Evolution of movement rate increases the effectiveness of marine reserves for the conservation of pelagic fishes.

Current debates about the efficacy of no-take marine reserves (MR) in protecting large pelagic fish such as tuna and sharks have usually not considered the evolutionary dimension of this issue, which emerges because the propensity to swim away from a given place, like any other biological trait, will probably vary in a heritable fashion among individuals. Here, based on spatially explicit simulations, we investigated whether selection to remain in MRs to avoid higher fishing mortality can lead to the evolution of more philopatric fish. Our simulations, which covered a range of life histories among tuna species (skipjack tuna vs. Atlantic bluefin tuna) and shark species (great white sharks vs. spiny dogfish), suggested that MRs were most effective at maintaining viable population sizes when movement distances were lowest. Decreased movement rate evolved following the establishment of marine reserves, and this evolution occurred more rapidly with higher fishing pressure. Evolutionary reductions in movement rate led to increases in within-reserve population sizes over the course of the 50 years following MR establishment, although this varied among life histories, with skipjack responding fastest and great white sharks slowest. Our results suggest the evolution of decreased movement can augment the efficacy of marine reserves, especially for species, such as skipjack tuna, with relatively short generation times. Even when movement rates did not evolve substantially over 50 years (e.g., given long generation times or little heritable variation), marine reserves were an effective tool for the conservation of fish populations when mean movement rates were low or MRs were large.

Interaction of ecological and angler processes: experimental stocking in an open access, spatially structured fishery.

Effective management of socioecological systems requires an understanding of the complex interactions between people and the environment. In recreational fisheries, which are prime examples of socioecological systems, anglers are analogous to mobile predators in natural predator-prey systems, and individual fisheries in lakes across a region are analogous to a spatially structured landscape of prey patches. Hence, effective management of recreational fisheries across large spatial scales requires an understanding of the dynamic interactions among ecological density dependent processes, landscape-level characteristics, and angler behaviors. We focused on the stocked component of the open access rainbow trout (Oncorhynchus mykiss) fishery in British Columbia (BC), and we used an experimental approach wherein we manipulated stocking densities in a subset of 34 lakes in which we monitored angler effort, fish abundance, and fish size for up to seven consecutive years. We used an empirically derived relationship between fish abundance and fish size across rainbow trout populations in BC to provide a measure of catch-based fishing quality that accounts for the size-abundance trade off in this system. We replicated our experimental manipulation in two regions known to have different angler populations and broad-scale access costs. We hypothesized that angler effort would respond to variation in stocking density, resulting in spatial heterogeneity in angler effort but homogeneity in catch-based fishing quality within regions. We found that there is an intermediate stocking density for a given lake or region at which angler effort is maximized (i.e., an optimal stocking density), and that this stocking density depends on latent effort and lake accessibility. Furthermore, we found no clear effect of stocking density on our measure of catch-based fishing quality, suggesting that angler effort homogenizes catch-related attributes leading to an eroded relationship between stocking density and catch-based fishing quality at the timescale of annual surveys. We conclude that declines in fishing quality resulting from understocking (due to declines in catch rate with low fish abundance) and overstocking (due to suppressed growth and limited recruitment at high density) give an optimal stocking rate that depends on accessibility and latent effort.

Supply-demand equilibria and the size-number trade-off in spatially structured recreational fisheries.

Recreational fishing effort varies across complex inland landscapes (e.g., lake-districts) and appears influenced by both angler preferences and qualities of the fishery resource, like fish size and abundance. However, fish size and abundance have an ecological trade-off within a population, thereby structuring equal-quality isopleths expressing this trade-off across the fishing landscape. Since expressed preferences of recreational anglers (i.e., site-selection of high-quality fishing opportunities among many lakes) can be analogous to optimal foraging strategies of natural predators, adopting such concepts can aid in understanding scale-dependence in fish-angler interactions and impacts of fishing across broad landscapes. Here, we assumed a fish supply-angler demand equilibria and adapted a novel bivariate measure of fishing quality based on fish size and catch rates to assess how recreational anglers influence fishing quality among a complex inland landscape. We then applied this metric to evaluate (1) angler preferences for caught and released fish compared to harvested fish, (2) the nonlinear size-numbers trade-off with uncertainty in both traits, and (3) the spatial-scale of the equilibria across 62 lakes and four independent management regions in British Columbia's (BC) rainbow trout Oncorhynchus mykiss fishery. We found anglers had low preference for caught and released fish (~10% of the value compared to harvested fish), which modified anglers' perception of fishing quality. Hence, fishing quality and angler effort was not influenced simply by total fish caught, but largely by harvested fish catch rates. Fishing quality varied from BC's northern regions (larger fish and more abundant) compared to southern regions (smaller fish and less abundant) directly associated with a 2.5 times increase in annual fishing effort in southern regions, suggesting that latent fishing pressure can structure the size-numbers trade-off in rainbow trout populations. The presence of two different equal-quality isopleths suggests at least two effective landscapes support co-occurring ideal free distributions of recreational fishing effort in BC's rainbow fishery. Anglers' expressed preferences among lakes interacted with density dependent growth and survival within lakes to structure a size-numbers trade-off influencing how anglers perceive fishing quality and, ultimately, distribute across complex inland landscapes.

Rainbow trout in seasonal environments: phenotypic trade-offs across a gradient in winter duration.

Survival through periods of resource scarcity depends on the balance between metabolic demands and energy storage. The opposing effects of predation and starvation mortality are predicted to result in trade-offs between traits that optimize fitness during periods of resource plenty (e.g., during the growing season) and those that optimize fitness during periods of resource scarcity (e.g., during the winter). We conducted a common environment experiment with two genetically distinct strains of rainbow trout to investigate trade-offs due to (1) the balance of growth and predation risk related to foraging rate during the growing season and (2) the allocation of energy to body size prior to the winter. Fry (age 0) from both strains were stocked into replicate natural lakes at low and high elevation that differed in winter duration (i.e., ice cover) by 59 days. Overwinter survival was lowest in the high-elevation lakes for both strains. Activity rate and growth rate were highest at high elevation, but growing season survival did not differ between strains or between environments. Hence, we did not observe a trade-off between growth and predation risk related to foraging rate. Growth rate also differed significantly between the strains across both environments, which suggests that growth rate is involved in local adaptation. There was not, however, a difference between strains or between environments in energy storage. Hence, we did not observe a trade-off between growth and storage. Our findings suggest that intrinsic metabolic rate, which affects a trade-off between growth rate and overwinter survival, may influence local adaptation in organisms that experience particularly harsh winter conditions (e.g., extended periods trapped beneath the ice in high-elevation lakes) in some parts of their range.

Identifying designatable units for intraspecific conservation prioritization: a hierarchical approach applied to the lake whitefish species complex (Coregonus spp.).

The concept of the designatable unit (DU) affords a practical approach to identifying diversity below the species level for conservation prioritization. However, its suitability for defining conservation units in ecologically diverse, geographically widespread and taxonomically challenging species complexes has not been broadly evaluated. The lake whitefish species complex (Coregonus spp.) is geographically widespread in the Northern Hemisphere, and it contains a great deal of variability in ecology and evolutionary legacy within and among populations, as well as a great deal of taxonomic ambiguity. Here, we employ a set of hierarchical criteria to identify DUs within the Canadian distribution of the lake whitefish species complex. We identified 36 DUs based on (i) reproductive isolation, (ii) phylogeographic groupings, (iii) local adaptation and (iv) biogeographic regions. The identification of DUs is required for clear discussion regarding the conservation prioritization of lake whitefish populations. We suggest conservation priorities among lake whitefish DUs based on biological consequences of extinction, risk of extinction and distinctiveness. Our results exemplify the need for extensive genetic and biogeographic analyses for any species with broad geographic distributions and the need for detailed evaluation of evolutionary history and adaptive ecological divergence when defining intraspecific conservation units.

Repeat swimming performance and its implications for inferring the relative fitness of asexual hybrid dace (Pisces: Phoxinus) and their sexually reproducing parental species.

While theories explaining the evolution and maintenance of sex are abundant, empirical data on the costs and benefits of asexual relative to sexual reproduction are less common. Asexually reproducing vertebrates, while few, provide a rare opportunity to measure differences in fitness between asexual and sexual species. All known asexually reproducing vertebrates are of hybrid origin, and hybrid disadvantage (i.e., reduced hybrid fitness) is thought to facilitate long-term coexistence between asexual and sexual species. We used repeat swimming performance as a proxy for fitness to compare the fitness of asexual hybrid dace (Pisces: Phoxinus) and their sexually reproducing parental species, finescale dace (Phoxinus neogaeus) and northern redbelly dace (Phoxinus eos). We tested the prediction that, given the widespread coexistence of these hybrid and parental dace, the parental species should show equivalent and perhaps superior repeat performance relative to hybrids. A repeat constant acceleration test (U(max)) was conducted at both acclimation temperature (16 °C) and at an elevated temperature (25 °C) to simulate the combined influence of a repeat swim and acute temperature change that fish might experience in the wild. The asexual hybrids performed more poorly than at least one of the parental species. There was a negative effect of temperature on repeat swimming performance in all fish, and the repeat performance of hybrids was more severely affected by temperature than that of finescale dace. No difference in the effect of temperature on repeat performance was detected between hybrids and northern redbelly dace. These results suggest that hybrids suffer physiological costs relative to the parentals or at least that the hybrids do not gain advantage from hybrid vigor, which probably contributes to the coexistence of asexual and sexual species in this system.

Distribution of Phoxinus eos, Phoxinus neogaeus, and their asexually-reproducing hybrids (Pisces: Cyprinidae) in Algonquin Provincial Park, Ontario.

Hybrid Phoxinus are one of the few asexually reproducing vertebrates species. The distribution of hybrid Phoxinus among lakes in Algonquin Park, Ontario, was evaluated relative to the distribution of parental species and relative to physiochemical lake characteristics. No association between the distribution of hybrids and the distribution of parental species was found, suggesting that the hybrids can successfully coexist with either parental species. In addition, we found no association between hybrid distribution and the physiochemical characteristics of lakes, suggesting that the hybrids are generalists with respect to the ecological niches available in Algonquin Park. There was a difference between the physiochemical characteristics of lakes with and without the parental species P. neogaeus. The lakes containing P. neogaeus were lower elevation than the lakes containing the other parental species, P. eos. The difference in distribution between the parental species may therefore be due to different dispersal abilities, to later arrival following post-glacial dispersal, or to differences in ecology. These results suggest that asexual reproduction is a successful strategy for hybrid Phoxinus in Algonquin Park because these sperm-dependent asexual hybrids are able to survive and persist regardless of which parental species is present, and regardless of the physiochemical characteristics of their habitat.

Variation in the strength of male mate choice allows long-term coexistence of sperm-dependent asexuals and their sexual hosts.

In several asexual taxa, reproduction requires mating with related sexual species to stimulate egg development, even though genetic material is not incorporated from the sexuals (gynogenesis). In cases in which gynogens do not invest in male function, they can potentially have a twofold competitive advantage over sexuals because the asexuals avoid the cost of producing males. If unmitigated, however, the competitive success of the asexuals would ultimately lead to their own demise, following the extinction of the sexual species that stimulate egg development. We have studied a model of mate choice among sexual individuals and asexual gynogens, where males of the sexual species preferentially mate with sexual females over gynogenetic females, to determine if such mating preferences can stably maintain both gynogenetic and sexual individuals within a community. Our model shows that stable coexistence of gynogens and their sexual hosts can occur when there is variation among males in the degree of preference for mating with sexual females and when pickier males pay a higher cost of preference.