Genomic stability through time despite decades of exploitation in cod on both sides of the Atlantic.

The mode and extent of rapid evolution and genomic change in response to human harvesting are key conservation issues. Although experiments and models have shown a high potential for both genetic and phenotypic change in response to fishing, empirical examples of genetic responses in wild populations are rare. Here, we compare whole-genome sequence data of Atlantic cod (Gadus morhua) that were collected before (early 20th century) and after (early 21st century) periods of intensive exploitation and rapid decline in the age of maturation from two geographically distinct populations in Newfoundland, Canada, and the northeast Arctic, Norway. Our temporal, genome-wide analyses of 346,290 loci show no substantial loss of genetic diversity and high effective population sizes. Moreover, we do not find distinct signals of strong selective sweeps anywhere in the genome, although we cannot rule out the possibility of highly polygenic evolution. Our observations suggest that phenotypic change in these populations is not constrained by irreversible loss of genomic variation and thus imply that former traits could be reestablished with demographic recovery.

Whole-genome sequencing reveals fine-scale environment-associated divergence near the range limits of a temperate reef fish.

Environmental variation is increasingly recognized as an important driver of diversity in marine species despite the lack of physical barriers to dispersal and the presence of pelagic stages in many taxa. A robust understanding of the genomic and ecological processes involved in structuring populations is lacking for most marine species, often hindering management and conservation action. Cunner (Tautogolabrus adspersus) is a temperate reef fish with both pelagic early life-history stages and strong site-associated homing as adults; the species is also of interest for use as a cleaner fish in salmonid aquaculture in Atlantic Canada. We aimed to characterize genomic and geographic differentiation of cunner in the Northwest Atlantic. To achieve this, a chromosome-level genome assembly for cunner was produced and used to characterize spatial population structure throughout Atlantic Canada using whole-genome sequencing. The genome assembly spanned 0.72 Gbp and 24 chromosomes; whole-genome sequencing of 803 individuals from 20 locations from Newfoundland to New Jersey identified approximately 11 million genetic variants. Principal component analysis revealed four regional Atlantic Canadian groups. Pairwise FST and selection scans revealed signals of differentiation and selection at discrete genomic regions, including adjacent peaks on chromosome 10 across multiple pairwise comparisons (i.e. FST 0.5-0.75). Redundancy analysis suggested association of environmental variables related to benthic temperature and oxygen range with genomic structure. Results suggest regional scale diversity in this temperate reef fish and can directly inform the collection and translocation of cunner for aquaculture applications and the conservation of wild populations throughout the Northwest Atlantic.

Environmentally associated chromosomal structural variation influences fine-scale population structure of Atlantic Salmon (Salmo salar).

Chromosomal rearrangements (e.g., inversions, fusions, and translocations) have long been associated with environmental variation in wild populations. New genomic tools provide the opportunity to examine the role of these structural variants in shaping adaptive differences within and among wild populations of non-model organisms. In Atlantic Salmon (Salmo salar), variations in chromosomal rearrangements exist across the species natural range, yet the role and importance of these structural variants in maintaining adaptive differences among wild populations remains poorly understood. We genotyped Atlantic Salmon (n = 1429) from 26 populations within a highly genetically structured region of southern Newfoundland, Canada with a 220K SNP array. Multivariate analysis, across two independent years, consistently identified variation in a structural variant (translocation between chromosomes Ssa01 and Ssa23), previously associated with evidence of trans-Atlantic secondary contact, as the dominant factor influencing population structure in the region. Redundancy analysis suggested that variation in the Ssa01/Ssa23 chromosomal translocation is strongly correlated with temperature. Our analyses suggest environmentally mediated selection acting on standing genetic variation in genomic architecture introduced through secondary contact may underpin fine-scale local adaptation in Placentia Bay, Newfoundland, Canada, a large and deep embayment, highlighting the importance of chromosomal structural variation as a driver of contemporary adaptive divergence.

Genomics reveal population structure, evolutionary history, and signatures of selection in the northern bottlenose whale, Hyperoodon ampullatus.

Information on wildlife population structure, demographic history, and adaptations are fundamental to understanding species evolution and informing conservation strategies. To study this ecological context for a cetacean of conservation concern, we conducted the first genomic assessment of the northern bottlenose whale, Hyperoodon ampullatus, using whole-genome resequencing data (n = 37) from five regions across the North Atlantic Ocean. We found a range-wide pattern of isolation-by-distance with a genetic subdivision distinguishing three subgroups: the Scotian Shelf, western North Atlantic, and Jan Mayen regions. Signals of elevated levels of inbreeding in the Endangered Scotian Shelf population indicate this population may be more vulnerable than the other two subgroups. In addition to signatures of inbreeding, evidence of local adaptation in the Scotian Shelf was detected across the genome. We found a long-term decline in effective population size for the species, which poses risks to their genetic diversity and may be exacerbated by the isolating effects of population subdivision. Protecting important habitat and migratory corridors should be prioritized to rebuild population sizes that were diminished by commercial whaling, strengthen gene flow, and ensure animals can move across regions in response to environmental changes.

Environmental Change, If Unaccounted, Prevents Detection of Cryptic Evolution in a Wild Population.

AbstractDetecting contemporary evolution requires demonstrating that genetic change has occurred. Mixed effects models allow estimation of quantitative genetic parameters and are widely used to study evolution in wild populations. However, predictions of evolution based on these parameters frequently fail to match observations. Here, we applied three commonly used quantitative genetic approaches to predict the evolution of size at maturity in a wild population of Trinidadian guppies. Crucially, we tested our predictions against evolutionary change observed in common-garden experiments performed on samples from the same population. We show that standard quantitative genetic models underestimated or failed to detect the cryptic evolution of this trait as demonstrated by the common-garden experiments. The models failed because (1) size at maturity and fitness both decreased with increases in population density, (2) offspring experienced higher population densities than their parents, and (3) selection on size was strongest at high densities. When we accounted for environmental change, predictions better matched observations in the common-garden experiments, although substantial uncertainty remained. Our results demonstrate that predictions of evolution are unreliable if environmental change is not appropriately captured in models.

Resistance and resilience of genetic and phenotypic diversity to "black swan" flood events: A retrospective analysis with historical samples of guppies.

Rare extreme "black swan" disturbances can impact ecosystems in many ways, such as destroying habitats, depleting resources, and causing high mortality. In rivers, for instance, exceptional floods that occur infrequently (e.g., so-called "50-year floods") can strongly impact the abundance of fishes and other aquatic organisms. Beyond such ecological effects, these floods could also impact intraspecific diversity by elevating genetic drift or dispersal and by imposing strong selection, which could then influence the population's ability to recover from disturbance. And yet, natural systems might be resistant (show little change) or resilient (show rapid recovery) even to rare extreme events - perhaps as a result of selection due to past events. We considered these possibilities in two rivers where native guppies experienced two extreme floods - one in 2005 and another in 2016. For each river, we selected four sites and used archived "historical" samples to compare levels of genetic and phenotypic diversity before vs. after floods. Genetic diversity was represented by 33 neutral microsatellite markers, and phenotypic diversity was represented by body length and male melanic (black) colour. We found that genetic diversity and population structure was mostly "resistant" to even these extreme floods; whereas the larger impacts on phenotypic diversity were short-lived, suggesting additional "resilience". We discuss the determinants of these two outcomes for guppies facing floods, and then consider the general implications for the resistance and resilience of intraspecific variation to black swan disturbances.

Genomic basis of deep-water adaptation in Arctic Charr (Salvelinus alpinus) morphs.

The post-glacial colonization of Gander Lake in Newfoundland, Canada, by Arctic Charr (Salvelinus alpinus) provides the opportunity to study the genomic basis of adaptation to extreme deep-water environments. Colonization of deep-water (>50 m) habitats often requires extensive adaptation to cope with novel environmental challenges from high hydrostatic pressure, low temperature, and low light, but the genomic mechanisms underlying evolution in these environments are rarely known. Here, we compare genomic divergence between a deep-water morph adapted to depths of up to 288 m and a larger, piscivorous pelagic morph occupying shallower depths. Using both a SNP array and resequencing of whole nuclear and mitochondrial genomes, we find clear genetic divergence (FST  = 0.11-0.15) between deep and shallow water morphs, despite an absence of morph divergence across the mitochondrial genome. Outlier analyses identified many diverged genomic regions containing genes enriched for processes such as gene expression and DNA repair, cardiac function, and membrane transport. Detection of putative copy number variants (CNVs) uncovered 385 genes with CNVs distinct to piscivorous morphs, and 275 genes with CNVs distinct to deep-water morphs, enriched for processes associated with synapse assembly. Demographic analyses identified evidence for recent and local morph divergence, and ongoing reductions in diversity consistent with postglacial colonization. Together, these results show that Arctic Charr morph divergence has occurred through genome-wide differentiation and elevated divergence of genes underlying multiple cellular and physiological processes, providing insight into the genomic basis of adaptation in a deep-water habitat following postglacial recolonization.

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 FST  = 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.

Spatio-temporal dynamics of density-dependent dispersal during a population colonisation.

Predicting population colonisations requires understanding how spatio-temporal changes in density affect dispersal. Density can inform on fitness prospects, acting as a cue for either habitat quality, or competition over resources. However, when escaping competition, high local density should only increase emigration if lower-density patches are available elsewhere. Few empirical studies on dispersal have considered the effects of density at the local and landscape scale simultaneously. To explore this, we analyze 5 years of individual-based data from an experimental introduction of wild guppies Poecilia reticulata. Natal dispersal showed a decrease in local density dependence as density at the landscape level increased. Landscape density did not affect dispersal among adults, but local density-dependent dispersal switched from negative (conspecific attraction) to positive (conspecific avoidance), as the colonisation progressed. This study demonstrates that densities at various scales interact to determine dispersal, and suggests that dispersal trade-offs differ across life stages.

Eco-Evolutionary Feedbacks Predict the Time Course of Rapid Life-History Evolution.

Organisms can change their environment and in doing so change the selection they experience and how they evolve. Population density is one potential mediator of such interactions because high population densities can impact the ecosystem and reduce resource availability. At present, such interactions are best known from theory and laboratory experiments. Here we quantify the importance of such interactions in nature by transplanting guppies from a stream where they co-occur with predators into tributaries that previously lacked both guppies and predators. If guppies evolve solely because of the immediate reduction in mortality rate, the strength of selection and rate of evolution should be greatest at the outset and then decline as the population adapts to its new environment. If indirect effects caused by the increase in guppy population density in the absence of predation prevail, then there should be a lag in guppy evolution because time is required for them to modify their environment. The duration of this lag is predicted to be associated with the environmental modification caused by guppies. We observed a lag in life-history evolution associated with increases in population density and altered ecology. How guppies evolved matched predictions derived from evolutionary theory that incorporates such density effects.

Chromosome polymorphisms track trans-Atlantic divergence and secondary contact in Atlantic salmon.

Pleistocene glaciations drove repeated range contractions and expansions shaping contemporary intraspecific diversity. Atlantic salmon (Salmo salar) in the western and eastern Atlantic diverged >600,000 years before present, with the two lineages isolated in different southern refugia during glacial maxima, driving trans-Atlantic genomic and karyotypic divergence. Here, we investigate the genomic consequences of glacial isolation and trans-Atlantic secondary contact using 108,870 single nucleotide polymorphisms genotyped in 80 North American and European populations. Throughout North America, we identified extensive interindividual variation and discrete linkage blocks within and between chromosomes with known trans-Atlantic differences in rearrangements: Ssa01/Ssa23 translocation and Ssa08/Ssa29 fusion. Spatial genetic analyses suggest independence of rearrangements, with Ssa01/Ssa23 showing high European introgression (>50%) in northern populations indicative of post-glacial trans-Atlantic secondary contact, contrasting with low European ancestry genome-wide (3%). Ssa08/Ssa29 showed greater intrapopulation diversity, suggesting a derived chromosome fusion polymorphism that evolved within North America. Evidence of potential selection on both genomic regions suggests that the adaptive role of rearrangements warrants further investigation in Atlantic salmon. Our study highlights how Pleistocene glaciations can influence large-scale intraspecific variation in genomic architecture of northern species.

Small-scale intraspecific patterns of adaptive immunogenetic polymorphisms and neutral variation in Lake Superior lake trout.

Many fishes express high levels of intraspecific variability, often linked to resource partitioning. Several studies show that a species' evolutionary trajectory of adaptive divergence can undergo reversals caused by changes in its environment. Such a reversal in neutral genetic and morphological variation among lake trout Salvelinus namaycush ecomorphs appears to be underway in Lake Superior. However, a water depth gradient in neutral genetic divergence was found to be associated with intraspecific diversity in the lake. To investigate patterns of adaptive immunogenetic variation among lake trout ecomorphs, we used Illumina high-throughput sequencing. The population's genetic structure of the major histocompatibility complex (MHC Class IIß exon 2) and 18 microsatellite loci were compared to disentangle neutral and selective processes at a small geographic scale. Both MHC and microsatellite variation were partitioned more by water depth stratum than by ecomorph. Several metrics showed strong clustering by water depth in MHC alleles, but not microsatellites. We report a 75% increase in the number of MHC alleles shared between the predominant shallow and deep water ecomorphs since a previous lake trout MHC study at the same locale (c. 1990s data). This result is consistent with the reverse speciation hypothesis, although adaptive MHC polymorphisms persist along an ecological gradient. Finally, results suggested that the lake trout have multiple copies of the MHC II locus consistent with a historic genomic duplication event. Our findings indicated that conservation approaches for this species could focus on managing various ecological habitats by depth, in addition to regulating the fisheries specific to ecomorphs.

Ancient chromosomal rearrangement associated with local adaptation of a postglacially colonized population of Atlantic Cod in the northwest Atlantic.

Intraspecific diversity is central to the management and conservation of exploited species, yet knowledge of how this diversity is distributed and maintained in the genome of many marine species is lacking. Recent advances in genomic analyses allow for genome-wide surveys of intraspecific diversity and offer new opportunities for exploring genomic patterns of divergence. Here, we analysed genome-wide polymorphisms to measure genetic differentiation between an offshore migratory and a nonmigratory population and to define conservation units of Atlantic Cod (Gadus morhua) in coastal Labrador. A total of 141 individuals, collected from offshore sites and from a coastal site within Gilbert Bay, Labrador, were genotyped using an ~11k single nucleotide polymorphism array. Analyses of population structure revealed strong genetic differentiation between migratory offshore cod and nonmigratory Gilbert Bay cod. Genetic differentiation was elevated for loci within a chromosomal rearrangement found on linkage group 1 (LG1) that coincides with a previously found double inversion associated with migratory and nonmigratory ecotype divergence of cod in the northeast Atlantic. This inverted region includes several genes potentially associated with adaptation to differences in salinity and temperature, as well as influencing migratory behaviour. Our work provides evidence that a chromosomal rearrangement on LG1 is associated with parallel patterns of divergence between migratory and nonmigratory ecotypes on both sides of the Atlantic Ocean.

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.

Genetic and phenotypic variation along an ecological gradient in lake trout Salvelinus namaycush.

BACKGROUND: Adaptive radiation involving a colonizing phenotype that rapidly evolves into at least one other ecological variant, or ecotype, has been observed in a variety of freshwater fishes in post-glacial environments. However, few studies consider how phenotypic traits vary with regard to neutral genetic partitioning along ecological gradients. Here, we present the first detailed investigation of lake trout Salvelinus namaycush that considers variation as a cline rather than discriminatory among ecotypes. Genetic and phenotypic traits organized along common ecological gradients of water depth and geographic distance provide important insights into diversification processes in a lake with high levels of human disturbance from over-fishing. RESULTS: Four putative lake trout ecotypes could not be distinguished using population genetic methods, despite morphological differences. Neutral genetic partitioning in lake trout was stronger along a gradient of water depth, than by locality or ecotype. Contemporary genetic migration patterns were consistent with isolation-by-depth. Historical gene flow patterns indicated colonization from shallow to deep water. Comparison of phenotypic (Pst) and neutral genetic variation (Fst) revealed that morphological traits related to swimming performance (e.g., buoyancy, pelvic fin length) departed more strongly from neutral expectations along a depth gradient than craniofacial feeding traits. Elevated phenotypic variance with increasing water depth in pelvic fin length indicated possible ongoing character release and diversification. Finally, differences in early growth rate and asymptotic fish length across depth strata may be associated with limiting factors attributable to cold deep-water environments. CONCLUSION: We provide evidence of reductions in gene flow and divergent natural selection associated with water depth in Lake Superior. Such information is relevant for documenting intraspecific biodiversity in the largest freshwater lake in the world for a species that recently lost considerable genetic diversity and is now in recovery. Unknown is whether observed patterns are a result of an early stage of incipient speciation, gene flow-selection equilibrium, or reverse speciation causing formerly divergent ecotypes to collapse into a single gene pool.

Don't bet against the natal homing abilities of marine fishes.

Whether marine fishes are capable of homing to their natal areas has long been something of an enigma. For some estuarine species or sharks (which have extended nondispersal juvenile stages or are born as relatively large, fully formed juveniles), the answer is clearly 'yes' (Thorrold et al. ; Feldheim et al. ), but for most marine fishes, the issue is much more mysterious. Many species have free-floating eggs, and most have pelagic, passively dispersing larvae. It is challenging to imagine how adult fish might navigate to a region of the ocean they experienced only as eggs or larvae, and easier to assume that such dispersal leads inexorably to high gene flow, and even panmixia. One way to resolve the conundrum would be to track fish from hatching to reproduction, but for marine fishes with tiny eggs and drifting larvae, this is notoriously difficult to do (Bradbury & Laurel ). In this issue of Molecular Ecology, Bonanomi et al. () use a creative approach to solve this challenge for Atlantic cod (Gadus morhua) populations that mingle in the vicinity of Greenland. They show that cod that disperse more than a 1000 km away from Iceland as eggs and larvae, then spend years growing on the far side of Greenland, while mixing with two local populations, return as adults to spawning areas near Iceland - and further, that this behaviour has remained stable over more than six decades. They manage this feat with a clever use of historical cod tracking data, modern genomic data and genetic analysis of decades-old DNA obtained from archived materials. Their results have important implications for our view of the biocomplexity of marine fish populations, and how we should manage them.

Selection analysis on the rapid evolution of a secondary sexual trait.

Evolutionary analyses of population translocations (experimental or accidental) have been important in demonstrating speed of evolution because they subject organisms to abrupt environmental changes that create an episode of selection. However, the strength of selection in such studies is rarely measured, limiting our understanding of the evolutionary process. This contrasts with long-term, mark-recapture studies of unmanipulated populations that measure selection directly, yet rarely reveal evolutionary change. Here, we present a study of experimental evolution of male colour in Trinidadian guppies where we tracked both evolutionary change and individual-based measures of selection. Guppies were translocated from a predator-rich to a low-predation environment within the same stream system. We used a combination of common garden experiments and monthly sampling of individuals to measure the phenotypic and genetic divergence of male coloration between ancestral and derived fish. Results show rapid evolutionary increases in orange coloration in both populations (1 year or three generations), replicating the results of previous studies. Unlike previous studies, we linked this evolution to an individual-based analysis of selection. By quantifying individual reproductive success and survival, we show, for the first time, that males with more orange and black pigment have higher reproductive success, but males with more black pigment also have higher risk of mortality. The net effect of selection is thus an advantage of orange but not black coloration, as reflected in the evolutionary response. This highlights the importance of considering all components of fitness when understanding the evolution of sexually selected traits in the wild.

Critical review of NGS analyses for de novo genotyping multigene families.

The genotyping of highly polymorphic multigene families across many individuals used to be a particularly challenging task because of methodological limitations associated with traditional approaches. Next-generation sequencing (NGS) can overcome most of these limitations, and it is increasingly being applied in population genetic studies of multigene families. Here, we critically review NGS bioinformatic approaches that have been used to genotype the major histocompatibility complex (MHC) immune genes, and we discuss how the significant advances made in this field are applicable to population genetic studies of gene families. Increasingly, approaches are introduced that apply thresholds of sequencing depth and sequence similarity to separate alleles from methodological artefacts. We explain why these approaches are particularly sensitive to methodological biases by violating fundamental genotyping assumptions. An alternative strategy that utilizes ultra-deep sequencing (hundreds to thousands of sequences per amplicon) to reconstruct genotypes and applies statistical methods on the sequencing depth to separate alleles from artefacts appears to be more robust. Importantly, the 'degree of change' (DOC) method avoids using arbitrary cut-off thresholds by looking for statistical boundaries between the sequencing depth for alleles and artefacts, and hence, it is entirely repeatable across studies. Although the advances made in generating NGS data are still far ahead of our ability to perform reliable processing, analysis and interpretation, the community is developing statistically rigorous protocols that will allow us to address novel questions in evolution, ecology and genetics of multigene families. Future developments in third-generation single molecule sequencing may potentially help overcome problems that still persist in de novo multigene amplicon genotyping when using current second-generation sequencing approaches.

Human disturbance causes the formation of a hybrid swarm between two naturally sympatric fish species.

Most evidence for hybrid swarm formation stemming from anthropogenic habitat disturbance comes from the breakdown of reproductive isolation between incipient species, or introgression between allopatric species following secondary contact. Human impacts on hybridization between divergent species that naturally occur in sympatry have received considerably less attention. Theory predicts that reinforcement should act to preserve reproductive isolation under such circumstances, potentially making reproductive barriers resistant to human habitat alteration. Using 15 microsatellites, we examined hybridization between sympatric populations of alewife (Alosa pseudoharengus) and blueback herring (A. aestivalis) to test whether the frequency of hybridization and pattern of introgression have been impacted by the construction of a dam that isolated formerly anadromous populations of both species in a landlocked freshwater reservoir. The frequency of hybridization and pattern of introgression differed markedly between anadromous and landlocked populations. The rangewide frequency of hybridization among anadromous populations was generally 0-8%, whereas all landlocked individuals were hybrids. Although neutral introgression was observed among anadromous hybrids, directional introgression leading to increased prevalence of alewife genotypes was detected among landlocked hybrids. We demonstrate that habitat alteration can lead to hybrid swarm formation between divergent species that naturally occur sympatrically, and provide empirical evidence that reinforcement does not always sustain reproductive isolation under such circumstances.

Variable parallelism in the genomic basis of age at maturity across spatial scales in Atlantic Salmon.

Complex traits often exhibit complex underlying genetic architectures resulting from a combination of evolution from standing variation, hard and soft sweeps, and alleles of varying effect size. Increasingly, studies implicate both large-effect loci and polygenic patterns underpinning adaptation, but the extent that common genetic architectures are utilized during repeated adaptation is not well understood. Sea age or age at maturation represents a significant life history trait in Atlantic Salmon (Salmo salar), the genetic basis of which has been studied extensively in European Atlantic populations, with repeated identification of large-effect loci. However, the genetic basis of sea age within North American Atlantic Salmon populations remains unclear, as does the potential for a parallel trans-Atlantic genomic basis to sea age. Here, we used a large single-nucleotide polymorphism (SNP) array and low-coverage whole-genome resequencing to explore the genomic basis of sea age variation in North American Atlantic Salmon. We found significant associations at the gene and SNP level with a large-effect locus (vgll3) previously identified in European populations, indicating genetic parallelism, but found that this pattern varied based on both sex and geographic region. We also identified nonrepeated sets of highly predictive loci associated with sea age among populations and sexes within North America, indicating polygenicity and low rates of genomic parallelism. Despite low genome-wide parallelism, we uncovered a set of conserved molecular pathways associated with sea age that were consistently enriched among comparisons, including calcium signaling, MapK signaling, focal adhesion, and phosphatidylinositol signaling. Together, our results indicate parallelism of the molecular basis of sea age in North American Atlantic Salmon across large-effect genes and molecular pathways despite population-specific patterns of polygenicity. These findings reveal roles for both contingency and repeated adaptation at the molecular level in the evolution of life history variation.