Allele Frequency Clines Research Articles

Circadian clock- and temperature-associated genes contribute to overall genomic differentiation along elevation in lichenized fungi.

Circadian regulation is linked to local environmental adaptation, and many species with broad climatic niches display variation in circadian genes. Here, we hypothesize that lichenizing fungi occupying different climate zones tune their metabolism to local environmental conditions with the help of their circadian systems. We study two species of the genus Umbilicaria occupying similar climatic niches (Mediterranean and the cold temperate) in different continents. Using homology to Neurospora crassa genes, we identify gene sets associated with circadian rhythms (11 core, 39 peripheral genes) as well as temperature response (37 genes). Nucleotide diversity of these genes is significantly correlated with mean annual temperature, minimum temperature of the coldest month and mean temperature of the coldest quarter. Furthermore, we identify altitudinal clines in allele frequencies in several non-synonymous substitutions in core clock components, for example, white collar-like, frh-like and various ccg-like genes. A dN/dS approach revealed a few significant peripheral clock- and temperature-associated genes (e.g. ras-1-like, gna-1-like) that may play a role in fine-tuning the circadian clock and temperature-response machinery. An analysis of allele frequency changes demonstrated the strongest evidence for differentiation above the genomic background in the clock-associated genes in U. pustulata. These results highlight the likely relevance of the circadian clock in environmental adaptation, particularly frost tolerance, of lichens. Whether or not the fungal clock modulates the symbiotic interaction within the lichen consortium remains to be investigated. We corroborate the finding of genetic variation in clock components along altitude-not only latitude-as has been reported in other species.

Archaic Introgression Shaped Human Circadian Traits.

When the ancestors of modern Eurasians migrated out of Africa and interbred with Eurasian archaic hominins, namely, Neanderthals and Denisovans, DNA of archaic ancestry integrated into the genomes of anatomically modern humans. This process potentially accelerated adaptation to Eurasian environmental factors, including reduced ultraviolet radiation and increased variation in seasonal dynamics. However, whether these groups differed substantially in circadian biology and whether archaic introgression adaptively contributed to human chronotypes remain unknown. Here, we traced the evolution of chronotype based on genomes from archaic hominins and present-day humans. First, we inferred differences in circadian gene sequences, splicing, and regulation between archaic hominins and modern humans. We identified 28 circadian genes containing variants with potential to alter splicing in archaics (e.g., CLOCK, PER2, RORB, and RORC) and 16 circadian genes likely divergently regulated between present-day humans and archaic hominins, including RORA. These differences suggest the potential for introgression to modify circadian gene expression. Testing this hypothesis, we found that introgressed variants are enriched among expression quantitative trait loci for circadian genes. Supporting the functional relevance of these regulatory effects, we found that many introgressed alleles have associations with chronotype. Strikingly, the strongest introgressed effects on chronotype increase morningness, consistent with adaptations to high latitude in other species. Finally, we identified several circadian loci with evidence of adaptive introgression or latitudinal clines in allele frequency. These findings identify differences in circadian gene regulation between modern humans and archaic hominins and support the contribution of introgression via coordinated effects on variation in human chronotype.

Landscape genetics of the protected Spanish Moon Moth in core, buffer, and peripheral areas of the Ordesa y Monte Perdido National Park (Central Pyrenees, Spain)

Conservation managers need to know the degree of connectivity showed by the populations to be preserved, especially when protected areas and/or species are involved. One of the conservation projects carried out by the Ordesa y Monte Perdido National Park (Central Pyrenees, Spain) is the monitoring of the protected Spanish Moon Moth, Graellsia isabellae (Lepidoptera, Saturniidae), in several sites within the actual park, buffer zone, and peripheral area. Here we studied the genetic diversity, geographical structure, and connectivity of this iconic insect in those areas with the aim of producing evidence-based recommendations that might help the National Park staff in their decision-making. For this, we non-lethally sampled 402 adult moths from 17 sites and worked at two geographic scales: Western/Central Pyrenees and the area monitored by the staff of the National Park. The multilocus genotypes obtained for nine nuclear microsatellite markers allowed us to quantify genetic variation, investigate population structure, and calculate recent migration rates. Our results revealed a large-scale (ca. 125 km) west–east cline in allele frequencies that causes low overall genetic differentiation (FST = 0.038) and similar levels of diversity among sites. Habitat connectivity revealed as an important element determining dispersal for G. isabellae, given the patchy distribution of the host plant (Pinus sylvestris) in the study area. Gene flow within and outside the National Park was proved, with a particular site of the buffer zone (Bujaruelo) acting as a source of migrants to other localities within and outside the National Park. This finding underlines the importance of considering buffer zones to preserve genetic diversity within protected areas, and that safeguarding the connectedness of pine patches is key to the conservation of this iconic moth.

The paradox of adaptive trait clines with nonclinal patterns in the underlying genes

Multivariate climate change presents an urgent need to understand how species adapt to complex environments. Population genetic theory predicts that loci under selection will form monotonic allele frequency clines with their selective environment, which has led to the wide use of genotype-environment associations (GEAs). This study used a set of simulations to elucidate the conditions under which allele frequency clines are more or less likely to evolve as multiple quantitative traits adapt to multivariate environments. Phenotypic clines evolved with nonmonotonic (i.e., nonclinal) patterns in allele frequencies under conditions that promoted unique combinations of mutations to achieve the multivariate optimum in different parts of the landscape. Such conditions resulted from interactions among landscape, demography, pleiotropy, and genetic architecture. GEA methods failed to accurately infer the genetic basis of adaptation under a range of scenarios due to first principles (clinal patterns did not evolve) or statistical issues (clinal patterns evolved but were not detected due to overcorrection for structure). Despite the limitations of GEAs, this study shows that a back-transformation of multivariate ordination can accurately predict individual multivariate traits from genotype and environmental data regardless of whether inference from GEAs was accurate. In addition, frameworks are introduced that can be used by empiricists to quantify the importance of clinal alleles in adaptation. This research highlights that multivariate trait prediction from genotype and environmental data can lead to accurate inference regardless of whether the underlying loci display clinal or nonmonotonic patterns.

Global biogeography of warning coloration in the butterfly Danaus chrysippus.

Warning coloration provides a textbook example of natural selection, but the frequent observation of polymorphism in aposematic species presents an evolutionary puzzle. We investigated biogeography and polymorphism of warning patterns in the widespread butterfly Danaus chrysippus using records from citizen science (n = 5467), museums (n = 8864) and fieldwork (n = 2586). We find that polymorphism in three traits controlled by known mendelian loci is extensive. Broad allele frequency clines, hundreds of kilometres wide, suggest a balance between long-range dispersal and predation of unfamiliar morphs. Mismatched clines for the white hindwing and forewing tip in East Africa are consistent with a previous finding that the black wingtip allele has spread recently in the region through hitchhiking with a heritable endosymbiont. Light/dark background coloration shows more extensive polymorphism. The darker genotype is more common in cooler regions, possibly reflecting a trade-off between thermoregulation and predator warning. Overall, our findings show how studying local adaptation at the global scale provides a more complete picture of the evolutionary forces involved.

The spatio-temporal dynamics of interacting genetic incompatibilities. Part I: the case of stacked underdominant clines.

We explore the interaction between two genetic incompatibilities (underdominant loci in diploid organisms) in a population occupying a one-dimensional space. We derive a system of partial differential equations describing the dynamics of allele frequencies and linkage disequilibrium between the two loci, and use a quasi-linkage equilibrium approximation in order to reduce the number of variables. We investigate the solutions of this system and demonstrate the existence of a solution in which the two clines in allele frequency remain stacked together. In the case of asymmetric incompatibilities (i.e. when one homozygote is favored over the other at each locus), these stacked clines propagate in the form of a traveling wave. We obtain an approximation for the speed of this wave which, in particular, is decreased by recombination between the two loci but is always larger than the speed of "one cline alone".

A Population Genomic Assessment of Three Decades of Evolution in a Natural Drosophila Population.

Population genetics seeks to illuminate the forces shaping genetic variation, often based on a single snapshot of genomic variation. However, utilizing multiple sampling times to study changes in allele frequencies can help clarify the relative roles of neutral and non-neutral forces on short time scales. This study compares whole-genome sequence variation of recently collected natural population samples of Drosophila melanogaster against a collection made approximately 35 years prior from the same locality—encompassing roughly 500 generations of evolution. The allele frequency changes between these time points would suggest a relatively small local effective population size on the order of 10,000, significantly smaller than the global effective population size of the species. Some loci display stronger allele frequency changes than would be expected anywhere in the genome under neutrality—most notably the tandem paralogs Cyp6a17 and Cyp6a23, which are impacted by structural variation associated with resistance to pyrethroid insecticides. We find a genome-wide excess of outliers for high genetic differentiation between old and new samples, but a larger number of adaptation targets may have affected SNP-level differentiation versus window differentiation. We also find evidence for strengthening latitudinal allele frequency clines: northern-associated alleles have increased in frequency by an average of nearly 2.5% at SNPs previously identified as clinal outliers, but no such pattern is observed at random SNPs. This project underscores the scientific potential of using multiple sampling time points to investigate how evolution operates in natural populations, by quantifying how genetic variation has changed over ecologically relevant timescales.

A continent-wide high genetic load in African buffalo revealed by clines in the frequency of deleterious alleles, genetic hitchhiking and linkage disequilibrium.

A high genetic load can negatively affect population viability and increase susceptibility to diseases and other environmental stressors. Prior microsatellite studies of two African buffalo (Syncerus caffer) populations in South Africa indicated substantial genome-wide genetic load due to high-frequency occurrence of deleterious alleles. The occurrence of these alleles, which negatively affect male body condition and bovine tuberculosis resistance, throughout most of the buffalo’s range were evaluated in this study. Using available microsatellite data (2–17 microsatellite loci) for 1676 animals from 34 localities (from 25°S to 5°N), we uncovered continent-wide frequency clines of microsatellite alleles associated with the aforementioned male traits. Frequencies decreased over a south-to-north latitude range (average per-locus Pearson r = -0.22). The frequency clines coincided with a multilocus-heterozygosity cline (adjusted R2 = 0.84), showing up to a 16% decrease in southern Africa compared to East Africa. Furthermore, continent-wide linkage disequilibrium (LD) at five linked locus pairs was detected, characterized by a high fraction of positive interlocus associations (0.66, 95% CI: 0.53, 0.77) between male-deleterious-trait-associated alleles. Our findings suggest continent-wide and genome-wide selection of male-deleterious alleles driven by an earlier observed sex-chromosomal meiotic drive system, resulting in frequency clines, reduced heterozygosity due to hitchhiking effects and extensive LD due to male-deleterious alleles co-occurring in haplotypes. The selection pressures involved must be high to prevent destruction of allele-frequency clines and haplotypes by LD decay. Since most buffalo populations are stable, these results indicate that natural mammal populations, depending on their genetic background, can withstand a high genetic load.

Clinal genomic analysis reveals strong reproductive isolation across a steep habitat transition in stickleback fish

How ecological divergence causes strong reproductive isolation between populations in close geographic contact remains poorly understood at the genomic level. We here study this question in a stickleback fish population pair adapted to contiguous, ecologically different lake and stream habitats. Clinal whole-genome sequence data reveal numerous genome regions (nearly) fixed for alternative alleles over a distance of just a few hundred meters. This strong polygenic adaptive divergence must constitute a genome-wide barrier to gene flow because a steep cline in allele frequencies is observed across the entire genome, and because the cline center closely matches the habitat transition. Simulations confirm that such strong divergence can be maintained by polygenic selection despite high dispersal and small per-locus selection coefficients. Finally, comparing samples from near the habitat transition before and after an unusual ecological perturbation demonstrates the fragility of the balance between gene flow and selection. Overall, our study highlights the efficacy of divergent selection in maintaining reproductive isolation without physical isolation, and the analytical power of studying speciation at a fine eco-geographic and genomic scale.

Pervasive Genomic Signatures of Local Adaptation to Altitude Across Highland Specialist Andean Hummingbird Populations.

Populations along steep environmental gradients are subject to differentiating selection that can result in local adaptation, despite countervailing gene flow, and genetic drift. In montane systems, where species are often restricted to narrow ranges of elevation, it is unclear whether the selection is strong enough to influence functional differentiation of subpopulations differing by a few hundred meters in elevation. We used targeted capture of 12 501 exons from across the genome, including 271 genes previously implicated in altitude adaptation, to test for adaptation to local elevations for 2 highland hummingbird species, Coeligena violifer (n = 62) and Colibri coruscans (n = 101). For each species, we described population genetic structure across the complex geography of the Peruvian Andes and, while accounting for this structure, we tested whether elevational allele frequency clines in single nucleotide polymorphisms (SNPs) showed evidence for local adaptation to elevation. Although the 2 species exhibited contrasting population genetic structures, we found signatures of clinal genetic variation with shifts in elevation in both. The genes with SNP-elevation associations included candidate genes previously discovered for high-elevation adaptation as well as others not previously identified, with cellular functions related to hypoxia response, energy metabolism, and immune function, among others. Despite the homogenizing effects of gene flow and genetic drift, natural selection on parts of the genome evidently optimizes elevation-specific cellular function even within elevation range-restricted montane populations. Consequently, our results suggest local adaptation occurring in narrow elevation bands in tropical mountains, such as the Andes, may effectively make them "taller" biogeographic barriers.

Molecular signatures of local adaptation to light in Norway spruce

Main conclusionTranscriptomic and exome capture analysis reveal an adaptive cline for shade tolerance in Norway spruce. Genes involved in the lignin pathway and immunity seem to play a potential role in contributing towards local adaptation to light.The study of natural variation is an efficient method to elucidate how plants adapt to local climatic conditions, a key process for the evolution of a species. Norway spruce is a shade-tolerant conifer in which the requirement of far-red light for growth increases latitudinally northwards. The objective of the study is to characterize the genetic control of local adaptation to light enriched in far-red in Norway spruce, motivated by a latitudinal gradient for the Red:Far-red (R:FR) ratio to which Norway spruce has been proven to be genetically adapted. We have established the genomic signatures of local adaptation by conducting transcriptomic (total RNA-sequencing) and genomic analyses (exome capture), for the identification of genes differentially regulated along the cline. RNA-sequencing revealed 274 differentially expressed genes in response to SHADE (low R:FR light), between the southern and northern natural populations in Sweden. Exome capture included analysis of a uniquely large data set (1654 trees) that revealed missense variations in coding regions of nine differentially expressed candidate genes, which followed a latitudinal cline in allele and genotype frequencies. These genes included five transcription factors involved in vital processes like bud-set/bud-flush, lignin pathway, and cold acclimation and other genes that take part in cell-wall remodeling, secondary cell-wall thickening, response to starvation, and immunity. Based on these results, we suggest that the northern populations might not only be able to adjust their growing season in response to low R:FR light, but they may also be better adapted towards disease resistance by up-regulation of the lignin pathway that is linked to immunity. This forms a concrete basis for local adaptation to light quality in Norway spruce, one of the most economically important conifer tree species in Sweden.

Hybrid zone genomics supports candidate species in Iberian Alytes obstetricans

Abstract While estimates of genetic divergence are increasingly used in molecular taxonomy, hybrid zone analyses can provide decisive evidence for evaluating candidate species. Applying a population genomic approach (RAD-sequencing) to a fine-scale transect sampling, we analyzed the transition between two Iberian subspecies of the common midwife toad (Alytes obstetricans almogavarii and A. o. pertinax) in Catalonia (northeastern Spain), which putatively diverged since the Plio-Pleistocene. Their hybrid zone was remarkably narrow, with extensive admixture restricted to a single locality (close to Tarragona), and congruent allele frequency clines for the mitochondrial (13 km wide) and the average nuclear genomes (16 km wide). We also fitted clines independently for 89 taxon-diagnostic SNPs: most of them behave like the nuclear background, but a subset (13%) is completely impermeable to gene flow and might be linked to barrier loci involved in hybrid incompatibilities. Assuming that midwife toads are able to disperse in the area of contact, we conclude that these taxa experience partial reproductive isolation and represent incipient species, i.e. Alytes almogavarii and Alytes obstetricans. Interestingly, their evolutionary age and mitochondrial divergence fall below the thresholds proposed in molecular systematics studies, emphasizing the difficulty of predicting the outcome of secondary contacts between young lineages entering the grey zone of speciation.

Effects of landscapes and range expansion on population structure and local adaptation.

Understanding the origin and distribution of genetic diversity across landscapes is critical for predicting the future of organisms in changing climates. This study investigated how adaptive and demographic forces have shaped diversity and population structure in Pinus densata, a keystone species on Qinghai-Tibetan Plateau (QTP). We examined the distribution of genomic diversity across the range of P.densata using exome capture sequencing. We applied spatially explicit tests to dissect the impacts of allele surfing, geographic isolation and environmental gradients on population differentiation and forecasted how this genetic legacy may limit the persistence of P.densata in future climates. We found that allele surfing from range expansion could explain the distribution of 39% of the c. 48000 genotyped single nucleotide polymorphisms (SNPs). Uncorrected, these allele frequency clines severely confounded inferences of selection. After controlling for demographic processes, isolation-by-environment explained 9.2-19.5% of the genetic structure, with c. 4.0% of loci being affected by selection. Allele surfing and genotype-environment associations resulted in genomic mismatch under projected climate scenarios. We illustrate that significant local adaptation, when coupled with reduced diversity as a result of demographic history, constrains potential evolutionary response to climate change. The strong signal of genomic vulnerability in P.densata may be representative for other QTP endemics.

Strong divergent selection at multiple loci in two closely related species of ragworts adapted to high and low elevations on Mount Etna.

Recently diverged species present particularly informative systems for studying speciation and maintenance of genetic divergence in the face of gene flow. We investigated speciation in two closely related Senecio species, S.aethnensis and S.chrysanthemifolius, which grow at high and low elevations, respectively, on Mount Etna, Sicily and form a hybrid zone at intermediate elevations. We used a newly generated genome-wide single nucleotide polymorphism (SNP) dataset from 192 individuals collected over 18 localities along an elevational gradient to reconstruct the likely history of speciation, identify highly differentiated SNPs, and estimate the strength of divergent selection. We found that speciation in this system involved heterogeneous and bidirectional gene flow along the genome, and species experienced marked population size changes in the past. Furthermore, we identified highly-differentiated SNPs between the species, some of which are located in genes potentially involved in ecological differences between species (such as photosynthesis and UV response). We analysed the shape of these SNPs' allele frequency clines along the elevational gradient. These clines show significantly variable coincidence and concordance, indicative of the presence of multifarious selective forces. Selection against hybrids is estimated to be very strong (0.16-0.78) and one of the highest reported in literature. The combination of strong cumulative selection across the genome and previously identified intrinsic incompatibilities probably work together to maintain the genetic and phenotypic differentiation between these species - pointing to the importance of considering both intrinsic and extrinsic factors when studying divergence and speciation.

Clinal variation in growth cessation and FTL2 expression in Siberian spruce

Forest trees exhibit strong patterns of local adaptation in phenological traits along latitudinal gradients. Previous studies in spruce have shown that variation at genes from the photoperiodic pathway and the circadian clock are associated to these clines but it has been difficult to find solid evidence of selection for some of these genes. Here, we used growth cessation, gene expression, and single nucleotide polymorphism (SNP) data at two major candidate loci, FLOWERING LOCUS T/TERMINAL FLOWER1-Like2 (FTL2) and GIGANTEA (GI), as well as at background loci from a latitudinal gradient in Siberian spruce (Picea obovata) populations along the Ob River to test for clinal variation in growth cessation and at the two candidate genes. As in previous studies, there was a strong latitudinal cline in growth cessation that was accompanied by a significant cline in the expression of FTL2. Expression of FTL2 was significantly associated with allele frequencies at some of the GI’s SNPs. However, the cline in allele frequency at candidate genes was not as steep as in a Norway spruce cline and in a parallel Siberian spruce cline studied previously and nonsignificant when a correction for population structure was applied. A McDonald-Kreitman test did not detect decisive evidence of selection on GI (p value = 0.07) and could not be applied to FTL2 because of limited polymorphism. Nonetheless, polymorphisms contributed more to the increased neutrality index of PoGI than to that of control loci. Finally, comparing the results of two previously published studies to our new dataset led to the identification of strong candidate SNPs for local adaptation in FTL2 promoter and GI.

Differential introgression across newt hybrid zones: Evidence from replicated transects.

Genomic heterogeneity of divergence between hybridizing species may reflect heterogeneity of introgression, but also processes unrelated to hybridization. Heterogeneous introgression and its repeatability can be directly tested in natural hybrid zones by examining multiple transects. Here, we studied hybrid zones between the European newts Lissotriton montandoni and two lineages of Lissotriton vulgaris, with replicate transects within each zone. Over 1,000 nuclear genes located on a linkage map and mitochondrial DNA were investigated using geographical and genomic clines. Overall, the five transects were all similar, showing hallmarks of strong reproductive isolation: bimodal distribution of genotypes in central populations and narrow allele frequency clines. However, the extent of introgression differed between the zones, possibly as a consequence of their different ages, as suggested by the analysis of heterozygosity runs in diagnostic markers. In three transects genomic signatures of small-scale (~2km) zone movements were detected. We found limited overlap of cline outliers between transects, and only weak evidence of stronger differentiation of introgression between zones than between transects within zones. Introgression was heterogeneous across linkage groups, with patterns of heterogeneity similar between transects and zones. Predefined candidates for increased or reduced introgression exhibited only a subtle tendency in the expected direction, suggesting that interspecific differentiation is not a reliable indicator for the strength of introgression. These hierarchically sampled hybrid zones of apparently different ages show how introgression unfolds with time and offer an excellent opportunity to dissect the dynamics of hybridization and architecture of reproductive isolation at advanced stages of speciation.

RAD-Seq Reveals Patterns of Additive Polygenic Variation Caused by Spatially-Varying Selection in the American Eel (Anguilla rostrata).

The American Eel (Anguilla rostrata) has an exceptional life cycle characterized by panmictic reproduction at the species scale, random dispersal, and selection in a highly heterogeneous habitat extending from subtropical to subarctic latitudes. The genetic consequences of spatially-varying selection in this species have been investigated for decades, revealing subtle clines in allele frequency at a few loci that contrast with complete panmixia on the vast majority of the genome. Because reproduction homogenizes allele frequencies every generation, sampling size, and genomic coverage are critical to reach sufficient power to detect selected loci in this context. Here, we used a total of 710 individuals from 12 sites and 12,098 high-quality single nucleotide polymorphisms to re-evaluate the extent to which local selection affects the spatial distribution of genetic diversity in this species. We used environmental association methods to identify markers under spatially-varying selection, which indicated that selection affects ∼1.5% of the genome. We then evaluated the extent to which candidate markers collectively vary with environmental factors using additive polygenic scores. We found significant correlations between polygenic scores and latitude, longitude and temperature which are consistent with polygenic selection acting against maladapted genotypes in different habitats occupied by eels throughout their range of distribution. Gene functions associated with outlier markers were significantly enriched for the insulin signaling pathway, indicating that the trade-offs inherent to occupying such a large distribution range involve the regulation of metabolism. Overall, this study highlights the potential of the additive polygenic scores approach in detecting selective effects in a complex environment.

Two-locus clines on the real line with a step environment

The shape of allele-frequency clines maintained by migration–selection balance depends not only on the properties of migration and selection, but also on the dominance relations among alleles and on linkage to other loci under selection. We investigate a two-locus model in which two diallelic, recombining loci are subject to selection caused by an abrupt environmental change. The habitat is one-dimensional and unbounded, selection at each locus is modeled by step functions such that in one region one allele at each locus is advantageous and in the other deleterious. We admit an environmentally independent, intermediate degree of dominance at both loci, including complete dominance. First, we derive an explicit expression for the single-locus cline with dominance, thus generalizing classical results by Haldane (1948). We show that the slope of the cline in the center (at the step) or, equivalently, the width of the cline, is independent of the degree of dominance. Second, under the assumption of strong recombination relative to selection and migration, the first-order approximations of the allele-frequency clines at each of the loci and of the linkage disequilibrium are derived. This may be interpreted as the quasi-linkage-equilibrium approximation of the two-locus cline. Explicit asymptotic expressions for the clines are deduced as x→±∞. For equivalent loci, explicit expressions for the whole clines are derived. The influence of dominance and of linkage on the slope of the cline in the center and on a global measure of steepness are investigated. This global measure reflects the influence of dominance. Finally, the accuracy of the approximations and the dependence of the shape of the two-locus cline on the full range of recombination rates is explored by numerical integration of the underlying system of partial differential equations.

RAD sequencing reveals genomewide divergence between independent invasions of the European green crab (Carcinus maenas) in the Northwest Atlantic.

Genomic studies of invasive species can reveal both invasive pathways and functional differences underpinning patterns of colonization success. The European green crab (Carcinus maenas) was initially introduced to eastern North America nearly 200 years ago where it expanded northwards to eastern Nova Scotia. A subsequent invasion to Nova Scotia from a northern European source allowed further range expansion, providing a unique opportunity to study the invasion genomics of a species with multiple invasions. Here, we use restriction‐site‐associated DNA sequencing‐derived SNPs to explore fine‐scale genomewide differentiation between these two invasions. We identified 9137 loci from green crab sampled from 11 locations along eastern North America and compared spatial variation to mitochondrial COI sequence variation used previously to characterize these invasions. Overall spatial divergence among invasions was high (pairwise FST ~0.001 to 0.15) and spread across many loci, with a mean FST ~0.052 and 52% of loci examined characterized by FST values >0.05. The majority of the most divergent loci (i.e., outliers, ~1.2%) displayed latitudinal clines in allele frequency highlighting extensive genomic divergence among the invasions. Discriminant analysis of principal components (both neutral and outlier loci) clearly resolved the two invasions spatially and was highly correlated with mitochondrial divergence. Our results reveal extensive cryptic intraspecific genomic diversity associated with differing patterns of colonization success and demonstrates clear utility for genomic approaches to delineating the distribution and colonization success of aquatic invasive species.

Identifying patterns of dispersal, connectivity and selection in the sea scallop,Placopecten magellanicus,usingRADseq‐derivedSNPs

Understanding patterns of dispersal and connectivity among marine populations can directly inform fisheries conservation and management. Advances in high‐throughput sequencing offer new opportunities for estimating marine connectivity. We used restriction‐site‐associated DNA sequencing to examine dispersal and realized connectivity in the sea scallop Placopecten magellanicus, an economically important marine bivalve. Based on 245 individuals sampled rangewide at 12 locations from Newfoundland to the Mid‐Atlantic Bight, we identified and genotyped 7163 single nucleotide polymorphisms; 112 (1.6%) were identified as outliers potentially under directional selection. Bayesian clustering revealed a discontinuity between northern and southern samples, and latitudinal clines in allele frequencies were observed in 42.9% of the outlier loci and in 24.6% of neutral loci. Dispersal estimates derived using these clines and estimates of linkage disequilibrium imply limited dispersal; 373.1 ± 407.0 km (mean ± SD) for outlier loci and 641.0 ± 544.6 km (mean ± SD) for neutral loci. Our analysis suggests restricted dispersal compared to the species range (>2000 km) and that dispersal and effective connectivity differ. These observations support the hypothesis that limited effective dispersal structures scallop populations along eastern North America. These findings can help refine the appropriate scale of management and conservation in this commercially valuable species.