Rapid genetic adaptation to a novel environment despite a genome-wide reduction in genetic diversity

Genetic divergence and one-way gene flow influence contemporary evolution and ecology of a partially migratory fish

Recent work has revealed the importance of contemporary evolution in shaping ecological outcomes. In particular, rapid evolutionary divergence between populations has been shown to impact the ecology of populations, communities, and ecosystems. While studies have focused largely on the role of adaptive divergence in generating ecologically important variation among populations, much less is known about the role of gene flow in shaping ecological outcomes. After divergence, populations may continue to interact through gene flow, which may influence evolutionary and ecological processes. Here, we investigate the role of gene flow in shaping the contemporary evolution and ecology of recently diverged populations of anadromous steelhead and resident rainbow trout (Oncorhynchus mykiss). Results show that resident rainbow trout introduced above waterfalls have diverged evolutionarily from downstream anadromous steelhead, which were the source of introductions. However, the movement of fish from above to below the waterfalls has facilitated gene flow, which has reshaped genetic and phenotypic variation in the anadromous source population. In particular, gene flow has led to an increased frequency of residency, which in turn has altered population density, size structure, and sex ratio. This result establishes gene flow as a contemporary evolutionary process that can have important ecological outcomes. From a management perspective, anadromous steelhead are generally regarded as a higher conservation priority than resident rainbow trout, even when found within the same watershed. Our results show that anadromous and resident O. mykiss populations may be connected via gene flow, with important ecological consequences. Such eco-evolutionary processes should be considered when managing recently diverged populations connected by gene flow.

Testing for adaptive changes linked to range expansion following a single introduction of the fall webworm

Adaptive evolution following colonization can affect the impact of invasive species. The fall webworm (FWW) invaded China 40 years ago through a single introduction event involving a severe bottleneck and subsequently diverged into two genetic groups. The well-recorded invasion history of FWW, coupled with a clear pattern of genetic divergence, provides an opportunity to investigate whether there is any sign of adaptive evolution following the invasion. Based on genome-wide SNPs, we identified genetically separated western and eastern groups of FWW and correlated spatial variation in SNPs with geographical and climatic factors. Geographical factors explained a similar proportion of the genetic variation across all populations compared with climatic factors. However, when the two population groups were analysed separately, environmental factors explained more variation than geographical factors. SNP outliers in populations of the western group had relatively stronger response to precipitation than temperature-related variables. Functional annotation of SNP outliers identified genes associated with insect cuticle protein potentially related to desiccation adaptation in the western group and genes associated with lipase biosynthesis potentially related to temperature adaptation in the eastern group. Our study suggests that invasive species may maintain the evolutionary potential to adapt to heterogeneous environments despite a single invasion event. The molecular data suggest that quantitative trait comparisons across environments would be worthwhile.

Genomics reveals repeated landlocking of diadromous fish on an isolated island

Landlocking of diadromous fish in freshwater systems can have significant genomic consequences. For instance, the loss of the migratory life stage can dramatically reduce gene flow across populations, leading to increased genetic structuring and stronger effects of local adaptation. These genomic consequences have been well-studied in some mainland systems, but the evolutionary impacts of landlocking in island ecosystems are largely unknown. In this study, we used a genotyping-by-sequencing (GBS) approach to examine the evolutionary history of landlocking in common smelt (Retropinna retropinna) on Chatham Island, a small isolated oceanic island 800 kilometres east of mainland New Zealand. We examined the relationship between Chatham Island and mainland smelt and used coalescent analyses to test the number and timing of landlocking events on Chatham Island. Our genomic analysis, based on 21,135 SNPs across 169 individuals, revealed that the Chatham Island smelt was genomically distinct from the mainland New Zealand fish, consistent with a single ancestral colonisation event of Chatham Island in the Pleistocene. Significant genetic structure was also evident within the Chatham Island smelt, with a diadromous Chatham Island smelt group, along with three geographically structured landlocked groups. Coalescent demographic analysis supported three independent landlocking events, with this loss of diadromy significantly pre-dating human colonisation. Our results illustrate how landlocking of diadromous fish can occur repeatedly across a narrow spatial scale, and highlight a unique system to study the genomic basis of repeated adaptation.

Genetic approaches reveal a healthy population and an unexpectedly recent origin for an isolated desert spring fish

Foskett Spring in Oregon's desert harbors a historically threatened population of Western Speckled Dace (Rhinichthys klamathensis). Though recently delisted, the dace's recruitment depends upon regular removal of encroaching vegetation. Previous studies assumed that Foskett Dace separated from others in the Warner Valley about 10,000 years ago, thereby framing an enigma about the population's surprising ability to persist for so long in a tiny habitat easily overrun by plants. To investigate that persistence and the effectiveness of interventions to augment population size, we assessed genetic diversity among daces inhabiting Foskett Spring, a refuge at Dace Spring, and three nearby streams. Analysis revealed a robust effective population size (Ne) of nearly 5000 within Foskett Spring, though Ne in the Dace Spring refuge is just 10% of that value. Heterozygosity is slightly lower than expected based on random mating at all five sites, indicating mild inbreeding, but not at a level of concern. These results confirm the genetic health of Foskett Dace. Unexpectedly, genetic differentiation reveals closer similarity between Foskett Dace and a newly discovered population from Nevada's Coleman Creek than between Foskett Dace and dace elsewhere in Oregon. Demographic modeling inferred Coleman Creek as the ancestral source of Foskett Dace fewer than 1000 years ago, much more recently than previously suspected and possibly coincident with the arrival of large herbivores whose grazing may have maintained open water suitable for reproduction. These results solve the enigma of persistence by greatly shortening the duration over which Foskett Dace have inhabited their isolated spring.

A lack of genetic diversity and minimal adaptive evolutionary divergence in introduced Mysis shrimp after 50 years

The successes of introduced populations in novel habitats often provide powerful examples of evolution and adaptation. In the 1950s, opossum shrimp (Mysis diluviana) individuals from Clearwater Lake in Minnesota, USA were transported and introduced to Twin Lakes in Colorado, USA by fisheries managers to supplement food sources for trout. Mysis were subsequently introduced from Twin Lakes into numerous lakes throughout Colorado. Because managers kept detailed records of the timing of the introductions, we had the opportunity to test for evolutionary divergence within a known time interval. Here, we used reduced representation genomic data to investigate patterns of genetic diversity, test for genetic divergence between populations, and for evidence of adaptive evolution within the introduced populations in Colorado. We found very low levels of genetic diversity across all populations, with evidence for some genetic divergence between the Minnesota source population and the introduced populations in Colorado. There was little differentiation among the Colorado populations, consistent with the known provenance of a single founding population, with the exception of the population from Gross Reservoir, Colorado. Demographic modeling suggests that at least one undocumented introduction from an unknown source population hybridized with the population in Gross Reservoir. Despite the overall low genetic diversity we observed, F ST outlier and environmental association analyses identified multiple loci exhibiting signatures of selection and adaptive variation related to elevation and lake depth. The success of introduced species is thought to be limited by genetic variation, but our results imply that populations with limited genetic variation can become established in a wide range of novel environments. From an applied perspective, the observed patterns of divergence between populations suggest that genetic analysis can be a useful forensic tool to determine likely sources of invasive species.

Diadromy Drives Elevated Rates of Trait Evolution and Ecomorphological Convergence in Clupeiformes (Herring, Shad, and Anchovies)

AbstractMigration can have a profound influence on rates and patterns of phenotypic evolution. Diadromy is the migration between marine and freshwater habitats for feeding and reproduction that can require individuals to travel tens to thousands of kilometers. The high energetic demands of diadromy are predicted to select for ecomorphological traits that maximize swimming and locomotor efficiency. Intraspecific studies have shown repeated instances of divergence among diadromous and nondiadromous populations in locomotor and foraging traits, which suggests that at a macroevolutionary scale diadromous lineages may experience convergent evolution onto one or multiple adaptive optima. We tested for differences in rates and patterns of phenotypic evolution among diadromous and nondiadromous lineages in Clupeiformes, a clade that has evolved diadromy more than 10 times. Our results show that diadromous clupeiforms show convergent evolution for some locomotor traits and faster rates of evolution, which we propose are adaptive responses to the locomotor demands of migration. We also find evidence that diadromous lineages show convergence into multiple regions of multivariate trait space and suggest that these respective trait spaces are associated with differences in migration and trophic ecology. However, not all locomotor traits and no trophic traits show evidence of convergence or elevated rates of evolution associated with diadromy. Our results show that long-distance migration influences the tempo and patterns of phenotypic evolution at macroevolutionary scales, but there is not a single diadromous syndrome.

Range-wide and temporal genomic analyses reveal the consequences of near-extinction in Swedish moose

Ungulate species have experienced severe declines over the past centuries through overharvesting and habitat loss. Even if many game species have recovered thanks to strict hunting regulation, the genome-wide impacts of overharvesting are still unclear. Here, we examine the temporal and geographical differences in genome-wide diversity in moose (Alces alces) over its whole range in Sweden by sequencing 87 modern and historical genomes. We found limited impact of the 1900s near-extinction event but local variation in inbreeding and load in modern populations, as well as suggestion of a risk of future reduction in genetic diversity and gene flow. Furthermore, we found candidate genes for local adaptation, and rapid temporal allele frequency shifts involving coding genes since the 1980s, possibly due to selective harvesting. Our results highlight that genomic changes potentially impacting fitness can occur over short time scales and underline the need to track both deleterious and selectively advantageous genomic variation.

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.

Population genomics analysis to identify ion and water transporter genes involved in the adaptation of Tibetan naked carps to brackish water

Understanding how evolutionary processes shape the genetic variations and influence the response of species to environmental alterations is critical for biodiversity conservation and molecular breeding. Gymnocypris przewalskii przewalskii is the only known cyprinid fish that dwells in the brackish water of Lake Qinghai on the Qinghai-Tibetan Plateau. To reveal the genetic basis of its adaptation to high salinity and alkalinity, whole-genome sequencing was performed in G. p. przewalskii and its freshwater relatives Gymnocypris eckloni and Gymnocypris przewalskii ganzihonensis. Compared with freshwater species, lower genetic diversity and higher linkage disequilibrium were observed in G. p. przewalskii. Selective sweep analysis identified 424 core-selective genes enriched in transport activities. Transfection analysis showed that genetic changes in the positively selected gene aquaporin 3 (AQP3) improved cell viability after salt treatment, suggesting its involvement in brackish water adaptation. Our analysis indicates that ion and water transporter genes experienced intensive selection, which might have contributed to the maintenance of high osmolality and ion content in G. p. przewalskii. The current study identified key molecules involved in the adaptation of fish to brackish water, providing valuable genomic resources for the molecular breeding of salt-tolerant fish.

How genomics can help biodiversity conservation

The availability of public genomic resources can greatly assist biodiversity assessment, conservation, and restoration efforts by providing evidence for scientifically informed management decisions. Here we survey the main approaches and applications in biodiversity and conservation genomics, considering practical factors, such as cost, time, prerequisite skills, and current shortcomings of applications. Most approaches perform best in combination with reference genomes from the target species or closely related species. We review case studies to illustrate how reference genomes can facilitate biodiversity research and conservation across the tree of life. We conclude that the time is ripe to view reference genomes as fundamental resources and to integrate their use as a best practice in conservation genomics.

Global invasion history and native decline of the common starling: insights through genetics

Few invasive birds are as globally successful as the Common or European Starling (Sturnus vulgaris). Native to the Palearctic, the starling has been intentionally introduced to North and South America, South Africa, Australia, and the Pacific Islands, enabling us to explore species traits that may contribute to its invasion success. Coupling the rich studies of life history and more recent explorations of genomic variation among invasions, we illustrate how eco-evolutionary dynamics shape the invasion success of this long-studied and widely distributed species. Especially informative is the comparison between Australian and North American invasions, because these populations colonized novel ranges concurrently and exhibit shared signals of selection despite distinct population histories. In this review, we describe population dynamics across the native and invasive ranges, identify putatively selected traits that may influence the starling’s spread, and suggest possible determinants of starling success world-wide. We also identify future opportunities to utilize this species as a model for avian invasion research, which will inform our understanding of species’ rapid evolution in response to environmental change.

Genetic variation of Cerastium alpinum L. from Babia Góra, a critically endangered species in Poland

Babia Góra massif is the only site of occurrence of the Cerastium alpinum L. in Poland, an arctic-alpine perennial plant with a wide distribution in North America, northwestern Asia, and Europe. To determine whether the isolated Polish populations are genetically distinct, we have performed an evaluation of C. alpinum from Babia Góra with the use of iPBS markers. A total number of 133 individuals of C. alpinum from seven populations representing four localizations of the species were analyzed, i.e., from Babia Góra (Poland), Alps (Switzerland), Nuolja massif (Sweden), and Kaffiøyra (Svalbard, Norway). Genetic analysis of all C. alpinum samples using eight PBS primers identified 262 bands, 79.4% of which were polymorphic. iPBS markers revealed low genetic diversity (average H_e = 0.085) and high population differentiation (F_ST = 0.617). AMOVA results confirmed that the majority of the genetic variation (62%) was recorded among populations. The grouping revealed by PCoA showed that C. alpinum from Svalbard is the most diverged population, C. alpinum from Switzerland and Sweden form a pair of similar populations, whereas C. alpinum from the Babia Góra form a heterogeneous group of four populations. Results of isolation by distance analysis suggested that the spatial distance is the most probable cause of the observed differentiation among populations. Although significant traces of a bottleneck effect were noted for all populations of C. alpinum from Babia Góra, the populations still maintain a low but significant level of genetic polymorphism. These results are of great importance for developing conservation strategies for this species in Poland.

Conservation Biology and Reproduction in a Time of Developmental Plasticity

The objective of this review is to ask whether, and how, principles in conservation biology may need to be revisited in light of new knowledge about the power of epigenetics to alter developmental pathways. Importantly, conservation breeding programmes, used widely by zoological parks and aquariums, may appear in some cases to reduce fitness by decreasing animals' abilities to cope when confronted with the 'wild side' of their natural habitats. Would less comfortable captive conditions lead to the selection of individuals that, despite being adapted to life in a captive environment, be better able to thrive if relocated to a more natural environment? While threatened populations may benefit from advanced reproductive technologies, these may actually induce undesirable epigenetic changes. Thus, there may be inherent risks to the health and welfare of offspring (as is suspected in humans). Advanced breeding technologies, especially those that aim to regenerate the rarest species using stem cell reprogramming and artificial gametes, may also lead to unwanted epigenetic modifications. Current knowledge is still incomplete, and therefore ethical decisions about novel breeding methods remain controversial and difficult to resolve.

Sympatric genetic divergence between early- and late-season weedy rice populations

Sympatric genetic divergence is the most appealing and controversial pattern in the theory of ecological speciation. Examples that support sympatric genetic divergence in plant species are extremely rare. Solid evidence of sympatric genetic divergence will provide deep insights for revealing the underlying mechanisms of ecological speciation. We analysed the total genomic DNA sequences of 120 weedy rice (WR; Oryza sativa f. spontanea) plants, representing three WR population pairs separately from three early- and late-season rice fields, in comparison with those of the co-occurring rice cultivars and other rice materials. We detected substantial genetic divergence within the pairs of the sympatric early- and late-season WR populations, although genetic divergence was unevenly distributed across the genomes. Restricted gene flow was determined between the sympatric WR populations, resulting in their distinct genetic structures. We also detected relatively low genetic diversity that was likely to be associated with stronger selection in early-season WR populations. Our findings provide strong evidence for sympatric genetic divergence between the WR populations in the same fields but in different seasons. We conclude that temporal isolation plays an important role in creating genetic divergence between sympatric populations/species in plants.

Genomic dynamics of brown trout populations released to a novel environment

Population translocations occur for a variety of reasons, from displacement due to climate change to human‐induced transfers. Such actions have adverse effects on genetic variation and understanding their microevolutionary consequences requires monitoring. Here, we return to an experimental release of brown trout (Salmo trutta) in order to monitor the genomic effects of population translocations. In 1979, fish from each of two genetically (F_ST = 0.16) and ecologically separate populations were simultaneously released, at one point in time, to a lake system previously void of brown trout. Here, whole‐genome sequencing of pooled DNA (Pool‐seq) is used to characterize diversity within and divergence between the introduced populations and fish inhabiting two lakes downstream of the release sites, sampled 30 years later (c. 5 generations). Present results suggest that while extensive hybridization has occurred, the two introduced populations are unequally represented in the lakes downstream of the release sites. One population, which is ecologically resident in its original habitat, mainly contributes to the lake closest to the release site. The other population, migratory in its natal habitat, is genetically more represented in the lake further downstream. Genomic regions putatively under directional selection in the new habitat are identified, where allele frequencies in both established populations are more similar to the introduced population stemming from a resident population than the migratory one. Results suggest that the microevolutionary consequences of population translocations, for example, hybridization and adaptation, can be rapid and that Pool‐seq can be used as an initial tool to monitor genome‐wide effects.

Population transcriptomics reveals the effect of gene flow on the evolution of range limits

One of the most important questions in evolutionary biology is how the spatial distribution of species is limited. Asymmetric gene flow from core populations is suggested to increase the number of poorly adapted immigrants in the populations at the range edge. Genetic load due to migration, i.e., migration load, should prevent adaptation to the local habitat, leading to decreases in distribution range via local extinction or the limiting range expansion. However, few experimental studies have examined the effects of immigration on fitness and natural selection within recipient populations. To investigate the influence of migration load on the evolution of distribution range, we performed field and laboratory observations as well as population transcriptomics for the common river snail, Semisulcospira reiniana. This species meets the conditions that migration from source populations can prevent local adaptation in a sink population because they inhabit the broader range of environments, including middle/upper reaches of a river and estuaries within a single river and they may be more vulnerable to being swept away by water currents due to lowered spontaneous (upward) locomotion activity. We found that river steepness was related to the lower distribution limit of S. reiniana, with a narrower distribution range in the steeper river. Population transcriptomic analysis showed that gene flow was heavily asymmetric from the upstream populations to downstream ones in the steep river, suggesting a greater migration load in the steep river. The number of genes putatively involved in adaptation to the local habitat was lower in the steep river than in the gentle river. Gene expression profiles suggested that individuals achieve better local adaptation in the gentle river. Laboratory experiments suggested that evolutionary differences in salinity tolerance among local populations were only found in the gentle river. Our results consistent with the hypothesis that migration load owing to asymmetric gene flow disturbs local adaptation and restricts the distribution range of river snails.

Genomic analysis unveils mechanisms of northward invasion and signatures of plateau adaptation in the Asian house rat

The Asian house rat (AHR), Rattus tanezumi, has recently invaded the northern half of China. The AHR is a highly adaptive rat species that has also successfully conquered the Qinghai-Tibet Plateau (QTP) and replaced the brown rat (BR), R. norvegicus, at the edge of the QTP. Here, we assembled a draft genome of the AHR and explored the mechanisms of its northward invasion and the genetic basis underlying plateau adaptation in this species. Population genomic analyses revealed that the northwardly invasive AHRs consisted of two independent and genetically distinct populations which might result from multiple independent primary invasion events. One invasive population exhibited reduced genetic diversity and distinct population structure compared with its source population, while the other displayed preserved genetic polymorphisms and little genetic differentiation from its source population. Genes involved in G-protein coupled receptors and carbohydrate metabolism may contribute to the local adaptation of northern AHRs. In particular, RTN4 was identified as a key gene for AHRs in the QTP that favours adaptation to high-altitude hypoxia. Coincidently, the physiological performance and transcriptome profiles of hypoxia-exposed rats both showed better hypoxia adaptation in AHRs than in BRs that failed to colonize the heart of the QTP, which may have facilitated the replacement of the BR population by the invading AHRs at the edge of the QTP. This study provides profound insights into the multiple origins of the northwardly invasive AHR and the great tolerance to hypoxia in this species.

A chromosome-anchored genome assembly for Lake Trout (Salvelinus namaycush)

Here, we present an annotated, chromosome‐anchored, genome assembly for Lake Trout (Salvelinus namaycush) – a highly diverse salmonid species of notable conservation concern and an excellent model for research on adaptation and speciation. We leveraged Pacific Biosciences long‐read sequencing, paired‐end Illumina sequencing, proximity ligation (Hi‐C) sequencing, and a previously published linkage map to produce a highly contiguous assembly composed of 7378 contigs (contig N50 = 1.8 Mb) assigned to 4120 scaffolds (scaffold N50 = 44.975 Mb). Long read sequencing data were generated using DNA from a female double haploid individual. 84.7% of the genome was assigned to 42 chromosome‐sized scaffolds and 93.2% of Benchmarking Universal Single Copy Orthologues were recovered, putting this assembly on par with the best currently available salmonid genomes. Estimates of genome size based on k‐mer frequency analysis were highly similar to the total size of the finished genome, suggesting that the entirety of the genome was recovered. A mitochondrial genome assembly was also produced. Self‐versus‐self synteny analysis allowed us to identify homeologs resulting from the salmonid specific autotetraploid event (Ss4R) as well as regions exhibiting delayed rediploidization. Alignment with three other salmonid genomes and the Northern Pike (Esox lucius) genome also allowed us to identify homologous chromosomes in related taxa. We also generated multiple resources useful for future genomic research on Lake Trout, including a repeat library and a sex‐averaged recombination map. A novel RNA sequencing data set for liver tissue was also generated in order to produce a publicly available set of annotations for 49,668 genes and pseudogenes. Potential applications of these resources to population genetics and the conservation of native populations are discussed.

Genetic load has potential in large populations but is realized in small inbred populations

Populations with higher genetic diversity and larger effective sizes have greater evolutionary capacity (i.e., adaptive potential) to respond to ecological stressors. We are interested in how the variation captured in protein-coding genes fluctuates relative to overall genomic diversity and whether smaller populations suffer greater costs due to their genetic load of deleterious mutations compared with larger populations. We analyzed individual whole-genome sequences (N = 74) from three different populations of Montezuma quail (Cyrtonyx montezumae), a small ground-dwelling bird that is sustainably harvested in some portions of its range but is of conservation concern elsewhere. Our historical demographic results indicate that Montezuma quail populations in the United States exhibit low levels of genomic diversity due in large part to long-term declines in effective population sizes over nearly a million years. The smaller and more isolated Texas population is significantly more inbred than the large Arizona and the intermediate-sized New Mexico populations we surveyed. The Texas gene pool has a significantly smaller proportion of strongly deleterious variants segregating in the population compared with the larger Arizona gene pool. Our results demonstrate that even in small populations, highly deleterious mutations are effectively purged and/or lost due to drift. However, we find that in small populations the realized genetic load is elevated because of inbreeding coupled with a higher frequency of slightly deleterious mutations that are manifested in homozygotes. Overall, our study illustrates how population genomics can be used to proactively assess both neutral and functional aspects of contemporary genetic diversity in a conservation framework while simultaneously considering deeper demographic histories.

Rapid genetic adaptation to recently colonized environments is driven by genes underlying life history traits

Background

Uncovering the mechanisms underlying rapid genetic adaptation can provide insight into adaptive evolution and shed light on conservation, invasive species control, and natural resource management. However, it can be difficult to experimentally explore rapid adaptation due to the challenges associated with propagating and maintaining species in captive environments for long periods of time. By contrast, many introduced species have experienced strong selection when colonizing environments that differ substantially from their native range and thus provide a “natural experiment” for studying rapid genetic adaptation. One such example occurred when sea lamprey (Petromyzon marinus), native to the northern Atlantic, naturally migrated into Lake Champlain and expanded their range into the Great Lakes via man-made shipping canals.

Results

Utilizing 368,886 genome-wide single nucleotide polymorphisms (SNPs), we calculated genome-wide levels of genetic diversity (i.e., heterozygosity and π) for sea lamprey collected from native (Connecticut River), native but recently colonized (Lake Champlain), and invasive (Lake Michigan) populations, assessed genetic differentiation between all populations, and identified candidate genes that responded to selection imposed by the novel environments. We observed a 14 and 24% reduction in genetic diversity in Lake Michigan and Lake Champlain populations, respectively, compared to individuals from the Connecticut River, suggesting that sea lamprey populations underwent a genetic bottleneck during colonization. Additionally, we identified 121 and 43 outlier genes in comparisons between Lake Michigan and Connecticut River and between Lake Champlain and Connecticut River, respectively. Six outlier genes that contained synonymous SNPs in their coding regions and two genes that contained nonsynonymous SNPs may underlie the rapid evolution of growth (i.e., GHR), reproduction (i.e., PGR, TTC25, STARD10), and bioenergetics (i.e., OXCT1, PYGL, DIN4, SLC25A15).

Conclusions

By identifying the genomic basis of rapid adaptation to novel environments, we demonstrate that populations of invasive species can be a useful study system for understanding adaptive evolution. Furthermore, the reduction in genome-wide levels of genetic diversity associated with colonization coupled with the identification of outlier genes underlying key life history traits known to have changed in invasive sea lamprey populations (e.g., growth, reproduction) illustrate the utility in applying genomic approaches for the successful management of introduced species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07553-x.

The inflated significance of neutral genetic diversity in conservation genetics

The current rate of species extinction is rapidly approaching unprecedented highs, and life on Earth presently faces a sixth mass extinction event driven by anthropogenic activity, climate change, and ecological collapse. The field of conservation genetics aims at preserving species by using their levels of genetic diversity, usually measured as neutral genome-wide diversity, as a barometer for evaluating population health and extinction risk. A fundamental assumption is that higher levels of genetic diversity lead to an increase in fitness and long-term survival of a species. Here, we argue against the perceived importance of neutral genetic diversity for the conservation of wild populations and species. We demonstrate that no simple general relationship exists between neutral genetic diversity and the risk of species extinction. Instead, a better understanding of the properties of functional genetic diversity, demographic history, and ecological relationships is necessary for developing and implementing effective conservation genetic strategies.

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.

Population genomics reveal rapid genetic differentiation in a recently invasive population of Rattus norvegicus

Background

Invasive species bring a serious effect on local biodiversity, ecosystems, and even human health and safety. Although the genetic signatures of historical range expansions have been explored in an array of species, the genetic consequences of contemporary range expansions have received little attention, especially in mammal species. In this study, we used whole-genome sequencing to explore the rapid genetic change and introduction history of a newly invasive brown rat (Rattus norvegicus) population which invaded Xinjiang Province, China in the late 1970s.

Results

Bayesian clustering analysis, principal components analysis, and phylogenetic analysis all showed clear genetic differentiation between newly introduced and native rat populations. Reduced genetic diversity and high linkage disequilibrium suggested a severe population bottleneck in this colonization event. Results of TreeMix analyses revealed that the introduced rats were derived from an adjacent population in geographic region (Northwest China). Demographic analysis indicated that a severe bottleneck occurred in XJ population after the split off from the source population, and the divergence of XJ population might have started before the invasion of XJ. Moreover, we detected 42 protein-coding genes with allele frequency shifts throughout the genome for XJ rats and they were mainly associated with lipid metabolism and immunity, which could be seen as a prelude to future selection analyses in the novel environment of XJ.

Conclusions

This study presents the first genomic evidence on genetic differentiation which developed rapidly, and deepens the understanding of invasion history and evolutionary processes of this newly introduced rat population. This would add to our understanding of how invasive species become established and aid strategies aimed at the management of this notorious pest that have spread around the world with humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12983-021-00387-z.

Investigating Diadromy in Fishes and Its Loss in an -Omics Era

Diadromy, the predictable movements of individuals between marine and freshwater environments, is biogeographically and phylogenetically widespread across fishes. Thus, despite the high energetic and potential fitness costs involved in moving between distinct environments, diadromy appears to be an effective life history strategy. Yet, the origin and molecular mechanisms that underpin this migratory behavior are not fully understood. In this review, we aim first to summarize what is known about diadromy in fishes; this includes the phylogenetic relationship among diadromous species, a description of the main hypotheses regarding its origin, and a discussion of the presence of non-migratory populations within diadromous species. Second, we discuss how recent research based on -omics approaches (chiefly genomics, transcriptomics, and epigenomics) is beginning to provide answers to questions on the genetic bases and origin(s) of diadromy. Finally, we suggest future directions for -omics research that can help tackle questions on the evolution of diadromy.

A century of intermittent eco-evolutionary feedbacks resulted in novel trait combinations in invasive Great Lakes alewives (Alosa pseudoharengus)

Species introductions provide opportunities to quantify rates and patterns of evolutionary change in response to novel environments. Alewives (Alosa pseudoharengus) are native to the East Coast of North America where they ascend coastal rivers to spawn in lakes and then return to the ocean. Some populations have become landlocked within the last 350 years and diverged phenotypically from their ancestral marine population. More recently, alewives were introduced to the Laurentian Great Lakes (~150 years ago), but these populations have not been compared to East Coast anadromous and landlocked populations. We quantified 95 years of evolution in foraging traits and overall body shape of Great Lakes alewives and compared patterns of phenotypic evolution of Great Lakes alewives to East Coast anadromous and landlocked populations. Our results suggest that gill raker spacing in Great Lakes alewives has evolved in a dynamic pattern that is consistent with responses to strong but intermittent eco-evolutionary feedbacks with zooplankton size. Following their initial colonization of Lakes Ontario and Michigan, dense alewife populations likely depleted large-bodied zooplankton, which drove a decrease in alewife gill raker spacing. However, the introduction of large, non-native zooplankton to the Great Lakes in later decades resulted in an increase in gill raker spacing, and present-day Great Lakes alewives have gill raker spacing patterns that are similar to the ancestral East Coast anadromous population. Conversely, contemporary Great Lakes alewife populations possess a gape width consistent with East Coast landlocked populations. Body shape showed remarkable parallel evolution with East Coast landlocked populations, likely due to a shared response to the loss of long-distance movement or migrations. Our results suggest the colonization of a new environment and cessation of migration can result in rapid parallel evolution in some traits, but contingency also plays a role, and a dynamic ecosystem can also yield novel trait combinations.

Whole genome resequencing reveals genomic regions associated with thermal adaptation in redband trout

Adaptation to local environments involves evolution of ecologically important traits and underlying physiological processes. Here, we used low coverage whole-genome resequencing (lcWGR) on individuals to identify genome regions involved in thermal adaptation in wild redband trout Oncorhynchus mykiss gairdneri, a subspecies of rainbow trout that inhabits ecosystems ranging from cold montane forests to high elevation deserts. This study includes allele frequency-based analyses for selective sweeps among populations, followed by multiple association tests for specific sets of phenotypes measured under thermal stress (acute and chronic survival/mortality; high or low cardiac performance groups). Depending on the groups in each set of analyses, sequencing reads covered 43%-75% of the genome at ≥15× and each analysis included millions of SNPs across the genome. In tests for selective sweeps among populations, a total of six chromosomal regions were significant. The further association tests for specific phenotypes revealed that the region on chromosome 4 was consistently the most significant and contains the cerk gene (ceramide kinase). This study provides insight into a potential genetic mechanism of local thermal adaptation and suggests cerk may be an important candidate gene. However, further validation of this cerk gene is necessary to determine if the association with cardiac performance results in a functional role to influence thermal performance when exposed to high water temperatures and hypoxic conditions.

Divergent and linked selection shape patterns of genomic differentiation between European and North American Atlantic salmon (Salmo salar)

As populations diverge many processes can shape genomic patterns of differentiation. Regions of high differentiation can arise due to divergent selection acting on selected loci, genetic hitchhiking of nearby loci, or through repeated selection against deleterious alleles (linked background selection); this divergence may then be further elevated in regions of reduced recombination. Atlantic salmon (Salmo salar) from Europe and North America diverged >600,000 years ago and despite some evidence of secondary contact, the majority of genetic data indicate substantial divergence between lineages. This deep divergence with potential gene flow provides an opportunity to investigate the role of different mechanisms that shape the genomic landscape during early speciation. Here, using 184,295 single nucleotide polymorphisms (SNPs) and 80 populations, we investigate the genomic landscape of differentiation across the Atlantic Ocean with a focus on highly differentiated regions and the processes shaping them. We found evidence of high (mean F_ST  = 0.26) and heterogeneous genomic differentiation between continents. Genomic regions associated with high trans-Atlantic differentiation ranged in size from single loci (SNPs) within important genes to large regions (1-3 Mbp) on four chromosomes (Ssa06, Ssa13, Ssa16 and Ssa19). These regions showed signatures consistent with selection, including high linkage disequilibrium, despite no significant reduction in recombination. Genes and functional enrichment of processes associated with differentiated regions may highlight continental differences in ocean navigation and parasite resistance. Our results provide insight into potential mechanisms underlying differences between continents, and evidence of near-fixed and potentially adaptive trans-Atlantic differences concurrent with a background of high genome-wide differentiation supports subspecies designation in Atlantic salmon.

Multiple decades of stocking has resulted in limited hatchery introgression in wild brook trout (Salvelinus fontinalis) populations of Nova Scotia

Many populations of freshwater fishes are threatened with losses, and increasingly, the release of hatchery individuals is one strategy being implemented to support wild populations. However, stocking of hatchery individuals may pose long-term threats to wild populations, particularly if genetic interactions occur between wild and hatchery individuals. One highly prized sport fish that has been heavily stocked throughout its range is the brook trout (Salvelinus fontinalis). In Nova Scotia, Canada, hatchery brook trout have been stocked since the early 1900s, and despite continued stocking efforts, populations have suffered declines in recent decades. Before this study, the genetic structure of brook trout populations in the province was unknown; however, given the potential negative consequences associated with hatchery stocking, it is possible that hatchery programs have adversely affected the genetic integrity of wild populations. To assess the influence of hatchery supplementation on wild populations, we genotyped wild brook trout from 12 river systems and hatchery brook trout from two major hatcheries using 100 microsatellite loci. Genetic analyses of wild trout revealed extensive population genetic structure among and within river systems and significant isolation-by-distance. Hatchery stocks were genetically distinct from wild populations, and most populations showed limited to no evidence of hatchery introgression (<5% hatchery ancestry). Only a single location had a substantial number of hatchery-derived trout and was located in the only river where a local strain is used for supplementation. The amount of hatchery stocking within a watershed did not influence the level of hatchery introgression. Neutral genetic structure of wild populations was influenced by geography with some influence of climate and stocking indices. Overall, our study suggests that long-term stocking has not significantly affected the genetic integrity of wild trout populations, highlighting the variable outcomes of stocking and the need to evaluate the consequences on a case-by-case basis.

Comparing raccoon major histocompatibility complex diversity in native and introduced ranges: Evidence for the importance of functional immune diversity for adaptation and survival in novel environments

The adaptive potential of invasive species is related to the genetic diversity of the invader, which is influenced by genetic drift and natural selection. Typically, the genetic diversity of invaders is studied with neutral genetic markers; however, the expectation of reduced diversity has not been consistently supported by empirical studies. Here, we describe and interpret genetic diversity at both neutral microsatellite loci and the immune-related MHC-DRB locus of native and invasive populations of raccoon to better understand of how drift and selection impact patterns of genetic diversity during the invasion process. We found that despite the loss of many MHC (major histocompatibility complex) alleles in comparison with native populations, functional MHC supertypes are preserved in the invasive region. In the native raccoon population, the number of supertypes within individuals was higher than expected under a neutral model. The high level of individual functional divergence may facilitate the adaptation to local conditions in the invasive range. In the invasive populations, we also detected increased population structure at microsatellites compared to the MHC locus, further suggesting that balancing selection is acting on adaptively important regions of the raccoon genome. Finally, we found that alleles known to exhibit resistance to rabies in the native range, Prlo-DRB*4, Prlo-DRB*16 and Prlo-DRB*102, were the most common alleles in the European populations, suggesting directional selection is acting on this locus. Our research shows empirical support for the importance of functional immune diversity for adaptation and survival in novel environments.

Impact of temperature shift on gill physiology during smoltification of Atlantic salmon smolts (Salmo salar L.)

Exposure to a temperature increase may disrupt smoltification and delay or stop the downstream migration of smolts. Thermal regimes are often different between a river and its tributaries, but the effects of a relative temperature shift are not well described. We used expression of smoltification genes coupled with gill Na⁺/K⁺-ATPase activity (NKA) and plasma cortisol and growth hormone (GH) levels to investigate the impact of a 5 °C difference between tributary and river on salmon juveniles. Responses to a temperature challenge were examined at four time points during the smoltification period, with juveniles reared under three regimes including control, early and late temperature increase. The temperature shifts reduced gill NKA, plasma GH and cortisol levels which indicate hypo-osmoregulation impairment and may reduce the survival of smolts. Out of the 22 genes examined, the expression of six genes was influenced by the temperature treatments, while changes in further eleven genes were influenced by the date of sampling. Genes usually known to be upregulated during smoltification were downregulated after the temperature increase, notably nkaα1b, nkcc1a and igf1r. Upregulation of some genes involved in the hormonal regulation and acid-base equilibrium in early June may indicate a switch towards desmoltification. This study gives further insights about the impact of temperature increase on the molecular processes underlying smoltification and possible responses to human-related water temperature increase. The data also suggest dual roles in the smoltification and desmoltification for GH and IGF1 and points to the implication of genes in the smoltification process, that have previously been unstudied (nbc) or with little data available (igf2).

Decimated little brown bats show potential for adaptive change

The degree to which species can rapidly adapt is key to survival in the face of climatic and other anthropogenic changes. For little brown bats (Myotis lucifugus), whose populations have experienced declines of over 90% because of the introduced fungal pathogen that causes white-nose syndrome (WNS), survival of the species may ultimately depend upon its capacity for adaptive change. Here, we present evidence of selectively driven change (adaptation), despite dramatic nonadaptive genomic shifts (genetic drift) associated with population declines. We compared the genetic makeups of wild survivors versus non-survivors of WNS, and found significant shifts in allele frequencies of genes associated with regulating arousal from hibernation (GABARB1), breakdown of fats (cGMP-PK1), and vocalizations (FOXP2). Changes at these genes are suggestive of evolutionary adaptation, given that WNS causes bats to arouse with unusual frequency from hibernation, contributing to premature depletion of fat reserves. However, whether these putatively adaptive shifts in allele frequencies translate into sufficient increases in survival for the species to rebound in the face of WNS is unknown.

Among-family variation in survival and gene expression uncovers adaptive genetic variation in a threatened fish

Variation in among-family transcriptional responses to different environmental conditions can help to identify adaptive genetic variation, even prior to a selective event. Coupling differential gene expression with formal survival analyses allows for the disentanglement of treatment effects, required for understanding how individuals plastically respond to environmental stressors, from the adaptive genetic variation responsible for differential survival. We combined these two approaches to investigate responses to an emerging conservation issue, thiamine (vitamin B₁ ) deficiency, in a threatened population of Atlantic salmon (Salmo salar). Thiamine is an essential vitamin that is increasingly limited in many ecosystems. In Lake Champlain, Atlantic salmon cannot acquire thiamine in sufficient quantities to support natural reproduction; fertilized eggs must be reared in hatcheries and treated with supplemental thiamine. We evaluated transcriptional responses (via RNA sequencing) to thiamine treatment across families and found 3,616 genes differentially expressed between control (no supplemental thiamine) and treatment individuals. Fewer genes changed expression equally across families (i.e., additively) than exhibited genotype × environment interactions in response to thiamine. Differentially expressed genes were related to known physiological effects of thiamine deficiency, including oxidative stress, cardiovascular irregularities and neurological abnormalities. We also identified 1,446 putatively adaptive genes that were strongly associated with among-family survival in the absence of thiamine treatment, many of which related to neurogenesis and visual perception. Our results highlight the utility of coupling RNA sequencing with formal survival analyses to identify candidate genes that underlie the among-family variation in survival required for an adaptive response to natural selection.

Increased Adaptive Variation Despite Reduced Overall Genetic Diversity in a Rapidly Adapting Invader

Invasive species often evolve rapidly following introduction despite genetic bottlenecks that may result from small numbers of founders; however, some invasions may not fit this “genetic paradox”. The invasive cane toad (Rhinella marina) displays high phenotypic variation across its introduced Australian range. Here, we used three genome-wide datasets to characterize their population structure and genetic diversity. We found that toads form three genetic clusters: 1) native range toads, 2) toads from the source population in Hawaii and long-established areas near introduction sites in Australia, and 3) toads from more recently established northern Australian sites. Although we find an overall reduction in genetic diversity following introduction, we do not see this reduction in loci putatively under selection, suggesting that genetic diversity may have been maintained at ecologically relevant traits, or that mutation rates were high enough to maintain adaptive potential. Nonetheless, toads encounter novel environmental challenges in Australia, and the transition between genetic clusters occurs at a point along the invasion transect where temperature rises and rainfall decreases. We identify environmentally associated loci known to be involved in resistance to heat and dehydration. This study highlights that natural selection occurs rapidly and plays a vital role in shaping the structure of invasive populations.

Assisted Evolution in Astrobiology-Convergence of Ecology and Evolutionary Biology within the Context of Planetary Colonization

In ecology and conservation biology, the concept of assisted evolution aims at the optimization of the resilience of organisms and populations to changing environmental conditions. What has hardly been considered so far is that this concept is also relevant for future astrobiological research, since in artificial extraterrestrial habitats (e.g., plants and insects in martian greenhouses) novel environmental conditions will also affect the survival and performance of organisms. The question therefore arises whether and how space-relevant organisms can be artificially adapted to the desired circumstances in advance. Based on several adaptation and acclimatization strategies in wild ecosystems of Earth, I discuss which methods can be considered for assisted evolution in the context of astrobiological research. This includes enhanced selective breeding, induction of epigenetic inheritance, and genetic engineering, as well as possible problems of these applications. This short overview article aims to stimulate an emerging discussion as to whether humans, which are already prominent drivers of Earth's evolution, should consider such interventions for future planetary colonization as well.

Discovery and validation of candidate smoltification gene expression biomarkers across multiple species and ecotypes of Pacific salmonids

Early marine survival of juvenile salmon is intimately associated with their physiological condition during smoltification and ocean entry. Smoltification (parr-smolt transformation) is a developmental process that allows salmon to acquire seawater tolerance in preparation for marine living. Traditionally, this developmental process has been monitored using gill Na⁺/K⁺-ATPase (NKA) activity or plasma hormones, but gill gene expression offers the possibility of another method. Here, we describe the discovery of candidate genes from gill tissue for staging smoltification using comparisons of microarray studies with particular focus on the commonalities between anadromous Rainbow trout and Sockeye salmon datasets, as well as a literature comparison encompassing more species. A subset of 37 candidate genes mainly from the microarray analyses was used for TaqMan quantitative PCR assay design and their expression patterns were validated using gill samples from four groups, representing three species and two ecotypes: Coho salmon, Sockeye salmon, stream-type Chinook salmon and ocean-type Chinook salmon. The best smoltification biomarkers, as measured by consistent changes across these four groups, were genes involved in ion regulation, oxygen transport and immunity. Smoltification gene expression patterns (using the top 10 biomarkers) were confirmed by significant correlations with NKA activity and were associated with changes in body brightness, caudal fin darkness and caudal peduncle length. We incorporate gene expression patterns of pre-smolt, smolt and de-smolt trials from acute seawater transfers from a companion study to develop a preliminary seawater tolerance classification model for ocean-type Chinook salmon. This work demonstrates the potential of gene expression biomarkers to stage smoltification and classify juveniles as pre-smolt, smolt or de-smolt.

Elevated Heterozygosity in Adults Relative to Juveniles Provides Evidence of Viability Selection on Eagles and Falcons

Viability selection yields adult populations that are more genetically variable than those of juveniles, producing a positive correlation between heterozygosity and survival. Viability selection could be the result of decreased heterozygosity across many loci in inbred individuals and a subsequent decrease in survivorship resulting from the expression of the deleterious alleles. Alternatively, locus-specific differences in genetic variability between adults and juveniles may be driven by forms of balancing selection, including heterozygote advantage, frequency-dependent selection, or selection across temporal and spatial scales. We use a pooled-sequencing approach to compare genome-wide and locus-specific genetic variability between 74 golden eagle (Aquila chrysaetos), 62 imperial eagle (Aquila heliaca), and 69 prairie falcon (Falco mexicanus) juveniles and adults. Although genome-wide genetic variability is comparable between juvenile and adult golden eagles and prairie falcons, imperial eagle adults are significantly more heterozygous than juveniles. This evidence of viability selection may stem from a relatively smaller imperial eagle effective population size and potentially greater genetic load. We additionally identify ~2000 single-nucleotide polymorphisms across the 3 species with extreme differences in heterozygosity between juveniles and adults. Many of these markers are associated with genes implicated in immune function or olfaction. These loci represent potential targets for studies of how heterozygote advantage, frequency-dependent selection, and selection over spatial and temporal scales influence survivorship in avian species. Overall, our genome-wide data extend previous studies that used allozyme or microsatellite markers and indicate that viability selection may be a more common evolutionary phenomenon than often appreciated.

Population Genomic Signatures of Genetic Structure and Environmental Selection in the Catadromous Roughskin Sculpin Trachidermus fasciatus

Understanding the patterns of genetic diversity and adaptation across species' range is crucial to assess its long-term persistence and determine appropriate conservation measures. The impacts of human activities on the genetic diversity and genetic adaptation to heterogeneous environments remain poorly understood in the marine realm. The roughskin sculpin (Trachidermus fasciatus) is a small catadromous fish, and has been listed as a second-class state protected aquatic animal since 1988 in China. To elucidate the underlying mechanism of population genetic structuring and genetic adaptations to local environments, RAD tags were sequenced for 202 individuals in nine populations across the range of T. fasciatus in China. The pairwise FST values over 9,271 filtered SNPs were significant except that between Dongying and Weifang. All the genetic clustering analysis revealed significant population structure with high support for eight distinct genetic clusters. Both the minor allele frequency spectra and Ne estimations suggested extremely small Ne in some populations (e.g., Qinhuangdao, Rongcheng, Wendeng, and Qingdao), which might result from recent population bottleneck. The strong genetic structure can be partly attributed to genetic drift and habitat fragmentation, likely due to the anthropogenic activities. Annotations of candidate adaptive loci suggested that genes involved in metabolism, development, and osmoregulation were critical for adaptation to spatially heterogenous environment of local populations. In the context of anthropogenic activities and environmental change, results of the present population genomic work provided important contributions to the understanding of genetic differentiation and adaptation to changing environments.

"Chromosomes and genes, spawned these fateful scenes": Rapid adaptation in an introduced fish

Humans by their very nature alter the distribution of species. Be it introduction of exotic species, habitat alterations or construction of barriers, anthropogenic changes provide novel experimental systems for the molecular ecologist to study evolutionary change. These events often provide a contradiction. Effective population sizes are generally low, and introduced populations are typically characterized by reduced diversity consistent with theoretical predictions of population bottlenecks and founder effects. However, despite reduced diversity, rapid change sometimes occurs. Identification of genomic regions associated with these rapid adaptive responses to novel selection pressures provides a window into genomic regions important in adaptive diversity, both in the novel and native ranges. These studies also provide an important means to estimate the pace of evolutionary change. In this issue, Willoughby et al. () compared the heterozygosity of steelhead (the anadromous form of rainbow trout Oncorhynchus mykiss) introduced into Lake Michigan in the late 1880s to the putative source population from the ancestral California range. After 25 generations of isolation in Lake Michigan, Willoughby et al. () found consistent genomewide reductions in genetic diversity as estimated by a measure of pooled heterozygosity. Despite this overall reduction in heterozygosity, three chromosomal regions showed signals of rapid adaptation and contained genes associated with osmoregulatory and wound-healing functions.

Long-term demographic and genetic effects of releasing captive-born individuals into the wild

Because of continued habitat destruction and species extirpations, the need to use captive breeding for conservation purposes has been increasing steadily. However, the long-term demographic and genetic effects associated with releasing captive-born individuals with varied life histories into the wild remain largely unknown. To address this question, we developed forward-time, agent-based models for 4 species with long-running captive-breeding and release programs: coho salmon (Oncorhynchus kisutch), golden lion tamarin (Leontopithecus rosalia), western toad (Anaxyrus boreas), and Whooping Crane (Grus americana). We measured the effects of supplementation by comparing population size and neutral genetic diversity in supplemented populations to the same characteristics in unaltered populations 100 years after supplementation ended. Releasing even slightly less fit captive-born individuals to supplement wild populations typically resulted in reductions in population sizes and genetic diversity over the long term when the fitness reductions were heritable (i.e., due to genetic adaptation to captivity) and populations continued to be regulated by density-dependent mechanisms over time. Negative effects for species with longer life spans and lower rates of population replacement were smaller than for species with shorter life spans and higher rates of population replacement. Programs that released captive-born individuals over fewer years or that avoided breeding individuals with captive ancestry had smaller reductions in population size and genetic diversity over the long term. Relying on selection in the wild to remove individuals with reduced fitness mitigated some negative demographic effects, but at a substantial cost to neutral genetic diversity. Our results suggest that conservation-focused captive-breeding programs should take measures to prevent even small amounts of genetic adaptation to captivity, quantitatively determine the minimum number of captive-born individuals to release each year, and fully account for the interactions among genetic adaptation to captivity, population regulation, and life-history variation.

Environmental extremes drive population structure at the northern range limit of Atlantic salmon in North America

Conservation of exploited species requires an understanding of both genetic diversity and the dominant structuring forces, particularly near range limits, where climatic variation can drive rapid expansions or contractions of geographic range. Here, we examine population structure and landscape associations in Atlantic salmon (Salmo salar) across a heterogeneous landscape near the northern range limit in Labrador, Canada. Analysis of two amplicon-based data sets containing 101 microsatellites and 376 single nucleotide polymorphisms (SNPs) from 35 locations revealed clear differentiation between populations spawning in rivers flowing into a large marine embayment (Lake Melville) compared to coastal populations. The mechanisms influencing the differentiation of embayment populations were investigated using both multivariate and machine-learning landscape genetic approaches. We identified temperature as the strongest correlate with genetic structure, particularly warm temperature extremes and wider annual temperature ranges. The genomic basis of this divergence was further explored using a subset of locations (n = 17) and a 220K SNP array. SNPs associated with spatial structuring and temperature mapped to a diverse set of genes and molecular pathways, including regulation of gene expression, immune response, and cell development and differentiation. The results spanning molecular marker types and both novel and established methods clearly show climate-associated, fine-scale population structure across an environmental gradient in Atlantic salmon near its range limit in North America, highlighting valuable approaches for predicting population responses to climate change and managing species sustainability.