Rapid genetic adaptation to a novel ecosystem despite a large founder event

Introduced and invasive species make excellent natural experiments for investigating rapid evolution. Here, we describe the effects of genetic drift and rapid genetic adaptation in pink salmon (Oncorhynchus gorbuscha) that were accidentally introduced to the Great Lakes via a single introduction event 31 generations ago. Using whole-genome resequencing for 134 fish spanning five sample groups across the native and introduced range, we estimate that the source population's effective population size was 146,886 at the time of introduction, whereas the founding population's effective population size was just 72-a 2040-fold decrease. As expected with a severe founder event, we show reductions in genome-wide measures of genetic diversity, specifically a 37.7% reduction in the number of SNPs and an 8.2% reduction in observed heterozygosity. Despite this decline in genetic diversity, we provide evidence for putative selection at 47 loci across multiple chromosomes in the introduced populations, including missense variants in genes associated with circadian rhythm, immunological response and maturation, which match expected or known phenotypic changes in the Great Lakes. For one of these genes, we use a species-specific agent-based model to rule out genetic drift and conclude our results support a strong response to selection occurring in a period gene (per2) that plays a predominant role in determining an organism's daily clock, matching large day length differences experienced by introduced salmon during important phenological periods. Together, these results inform how populations might evolve rapidly to new environments, even with a small pool of standing genetic variation.

Large genetic divergence underpins cryptic local adaptation across ecological and evolutionary gradients

Environmentally covarying local adaptation is a form of cryptic local adaptation in which the covariance of the genetic and environmental effects on a phenotype obscures the divergence between locally adapted genotypes. Here, we systematically document the magnitude and drivers of the genetic effect (V_G) for two forms of environmentally covarying local adaptation: counter- and cogradient variation. Using a hierarchical Bayesian meta-analysis, we calculated the overall effect size of V_G as 1.05 and 2.13 for populations exhibiting countergradient or cogradient variation, respectively. These results indicate that the genetic contribution to phenotypic variation represents a 1.05 to 2.13 s.d. change in trait value between the most disparate populations depending on if populations are expressing counter- or cogradient variation. We also found that while there was substantial variance among abiotic and biotic covariates, the covariates with the largest mean effects were temperature (2.41) and gamete size (2.81). Our results demonstrate the pervasiveness and large genetic effects underlying environmentally covarying local adaptation in wild populations and highlight the importance of accounting for these effects in future studies.

Incipient resistance to an effective pesticide results from genetic adaptation and the canalization of gene expression

The resistance of pest species to chemical controls has vast ecological, economic, and societal costs. In most cases, resistance is only detected after spreading throughout an entire population. Detecting resistance in its incipient stages, by comparison, provides time to implement preventative strategies. Incipient resistance can be detected by coupling standard toxicology assays with large-scale gene expression experiments. We apply this approach to a system where an invasive parasite, sea lamprey (Petromyzon marinus), has been treated with the highly effective pesticide 3-trifluoromethyl-4-nitrophenol (TFM) for 60 years. Toxicological experiments revealed that lamprey from treated populations did not have higher survival to TFM exposure than lamprey from untreated populations, demonstrating that full-fledged resistance has not yet evolved. In contrast, we find hundreds of genes differentially expressed in response to TFM in the population with the longest history of exposure, many of which relate to TFM's primary mode of action, the uncoupling of oxidative phosphorylation, and subsequent depletion of ATP. Three genes critical to oxidative phosphorylation, ATP5PB, PLCB1, and NDUFA9, were nearly fixed for alternative alleles in comparisons of SNPs between treated and untreated populations (FST  > 5 SD from the mean). ATP5PB encodes subunit b of ATP synthase and an additional subunit, ATP5F1B, was canalized for high expression in treated populations, but remained plastic in response to TFM treatment in individuals from the untreated population. These combined genomic and transcriptomic results demonstrate that an adaptive, genetic response to TFM is likely driving incipient resistance in a damaging pest species.

Thermal adaptation and phenotypic plasticity in a warming world: Insights from common garden experiments on Alaskan sockeye salmon

An important unresolved question is how populations of coldwater-dependent fishes will respond to rapidly warming water temperatures. For example, the culturally and economically important group, Pacific salmon (Oncorhynchus spp.), experience site-specific thermal regimes during early development that could be disrupted by warming. To test for thermal local adaptation and heritable phenotypic plasticity in Pacific salmon embryos, we measured the developmental rate, survival, and body size at hatching in two populations of sockeye salmon (Oncorhynchus nerka) that overlap in timing of spawning but incubate in contrasting natural thermal regimes. Using a split half-sibling design, we exposed embryos of 10 families from each of two populations to variable and constant thermal regimes. These represented both experienced temperatures by each population, and predicted temperatures under plausible future conditions based on a warming scenario from the downscaled global climate model (MIROC A1B scenario). We did not find evidence of thermal local adaptation during the embryonic stage for developmental rate or survival. Within treatments, populations hatched within 1 day of each other, on average, and among treatments, did not differ in survival in response to temperature. We did detect plasticity to temperature; embryos developed 2.5 times longer (189 days) in the coolest regime compared to the warmest regime (74 days). We also detected variation in developmental rates among families within and among temperature regimes, indicating heritable plasticity. Families exhibited a strong positive relationship between thermal variability and phenotypic variability in developmental rate but body length and mass at hatching were largely insensitive to temperature. Overall, our results indicated a lack of thermal local adaptation, but a presence of plasticity in populations experiencing contrasting conditions, as well as family-specific heritable plasticity that could facilitate adaptive change.