Genomics detects population structure within and between ocean basins in a circumpolar seabird: The white‐chinned petrel

The Southern Ocean represents a continuous stretch of circumpolar marine habitat, but the potential physical and ecological drivers of evolutionary genetic differentiation across this vast ecosystem remain unclear. We tested for genetic structure across the full circumpolar range of the white‐chinned petrel (Procellaria aequinoctialis) to unravel the potential drivers of population differentiation and test alternative population differentiation hypotheses. Following range‐wide comprehensive sampling, we applied genomic (genotyping‐by‐sequencing or GBS; 60,709 loci) and standard mitochondrial‐marker approaches (cytochrome b and first domain of control region) to quantify genetic diversity within and among island populations, test for isolation by distance, and quantify the number of genetic clusters using neutral and outlier (non‐neutral) loci. Our results supported the multi‐region hypothesis, with a range of analyses showing clear three‐region genetic population structure, split by ocean basin, within two evolutionary units. The most significant differentiation between these regions confirmed previous work distinguishing New Zealand and nominate subspecies. Although there was little evidence of structure within the island groups of the Indian or Atlantic oceans, a small set of highly‐discriminatory outlier loci could assign petrels to ocean basin and potentially to island group, though the latter needs further verification. Genomic data hold the key to revealing substantial regional genetic structure within wide‐ranging circumpolar species previously assumed to be panmictic.     

Disentangling the effects of geographic peripherality and habitat suitability on neutral and adaptive genetic variation in Swiss stone pine

It is generally accepted that the spatial distribution of neutral genetic diversity within a species’ native range mostly depends on effective population size, demographic history, and geographic position. However, it is unclear how genetic diversity at adaptive loci correlates with geographic peripherality or with habitat suitability within the ecological niche. Using exome‐wide genomic data and distribution maps of the Alpine range, we first tested whether geographic peripherality correlates with four measures of population genetic diversity at > 17,000 SNP loci in 24 Alpine populations (480 individuals) of Swiss stone pine (Pinus cembra) from Switzerland. To distinguish between neutral and adaptive SNP sets, we used four approaches (two gene diversity estimates, F_ST outlier test, and environmental association analysis) that search for signatures of selection. Second, we established ecological niche models for P. cembra in the study range and investigated how habitat suitability correlates with genetic diversity at neutral and adaptive loci. All estimates of neutral genetic diversity decreased with geographic peripherality, but were uncorrelated with habitat suitability. However, heterozygosity (H_e) at adaptive loci based on Tajima's D declined significantly with increasingly suitable conditions. No other diversity estimates at adaptive loci were correlated with habitat suitability. Our findings suggest that populations at the edge of a species' geographic distribution harbour limited neutral genetic diversity due to demographic properties. Moreover, we argue that populations from suitable habitats went through strong selection processes, are thus well adapted to local conditions, and therefore exhibit reduced genetic diversity at adaptive loci compared to populations at niche margins.     

Riverscape genomics of a threatened fish across a hydroclimatically heterogeneous river basin

Understanding how natural selection generates and maintains adaptive genetic diversity in heterogeneous environments is key to predicting the evolutionary response of populations to rapid environmental change. Detecting selection in complex spatial environments remains challenging, especially for threatened species where the effects of strong genetic drift may overwhelm signatures of selection. We carried out a basinwide riverscape genomic analysis in the threatened southern pygmy perch (Nannoperca australis), an ecological specialist with low dispersal potential. High‐resolution environmental data and 5162 high‐quality filtered SNPs were used to clarify spatial population structure and to assess footprints of selection associated with a steep hydroclimatic gradient and with human disturbance across the naturally and anthropogenically fragmented Murray–Darling Basin (Australia). Our approach included F_ST outlier tests to define neutral loci, and a combination of spatially explicit genotype–environment association analyses to identify candidate adaptive loci while controlling for the effects of landscape structure and shared population history. We found low levels of genetic diversity and strong neutral population structure consistent with expectations based on spatial stream hierarchy and life history. In contrast, variables related to precipitation and temperature appeared as the most important environmental surrogates for putatively adaptive genetic variation at both regional and local scales. Human disturbance also influenced the variation in candidate loci for adaptation, but only at a local scale. Our study contributes to understanding of adaptive evolution along naturally and anthropogenically fragmented ecosystems. It also offers a tangible example of the potential contributions of landscape genomics for informing in situ and ex situ conservation management of biodiversity.     

Evidence for the genetic basis and epistatic interactions underlying ocean‐ and river‐maturing ecotypes of Pacific Lamprey (Entosphenus tridentatus) returning to the Klamath River,California

Surveys of genomic variation have improved our understanding of the relationship between fitness‐related phenotypes and their underlying genetic basis. In some cases, single large‐effect genes have been found to underlie important traits; however, complex traits are expected to be under polygenic control and elucidation of multiple gene interactions may be required to fully understand the genetic basis of the trait. In this study, we investigated the genetic basis of the ocean‐ and river‐maturing ecotypes in anadromous Pacific lamprey (Entosphenus tridentatus). In Pacific lamprey, the ocean‐maturing ecotype is distinguished by advanced maturity of females (e.g., large egg mass) at the onset of freshwater migration relative to immature females of the river‐maturing ecotype. We examined a total of 219 adult Pacific lamprey that were collected at‐entry to the Klamath River over a 12‐month period. Each individual was genotyped at 308 SNPs representing known neutral and adaptive loci and measured at morphological traits, including egg mass as an indicator of ocean‐ and river‐maturing ecotype for females. The two ecotypes did not exhibit genetic structure at 148 neutral loci, indicating that ecotypic diversity exists within a single population. In contrast, we identified the genetic basis of maturation ecotypes in Pacific lamprey as polygenic, involving two unlinked gene regions that have a complex epistatic relationship. Importantly, these gene regions appear to show stronger effects when considered in gene interaction models than if just considered additive, illustrating the importance of considering epistatic effects and gene networks when researching the genetic basis of complex traits in Pacific lamprey and other species.     

Atypical panmixia in a European dolphin species (Delphinus delphis): implications for the evolution of diversity across oceanic boundaries

Despite the scarcity of geographical barriers in the ocean environment, delphinid cetaceans often exhibit marked patterns of population structure on a regional scale. The European coastline is a prime example, with species exhibiting population structure across well‐defined environmental boundaries. Here we undertake a comprehensive population genetic study on the European common dolphin (Delphinus delphis, based on 492 samples and 15 loci) and establish that this species shows exceptional panmixia across most of the study range. We found differentiation only between the eastern and western Mediterranean, consistent with earlier studies, and here use approximate Bayesian computations to explore different scenarios to explain the observed pattern. Our results suggest that a recent population bottleneck likely contributed significantly to the differentiation of the Eastern Mediterranean population (in Greek waters). This interpretation is consistent with independent census data that suggest a sharp population decline in the recent past. The implication is that an unperturbed population may currently show panmixia across the full study range. This exception to the more typical pattern of population structure seen for other regional dolphin species (and for common dolphin populations elsewhere in the world) suggests particular ecological or life‐history traits distinct to this species in European waters.     

Adaptation in temporally variable environments: stickleback armor in periodically breaching bar‐built estuaries

The evolutionary consequences of temporal variation in selection remain hotly debated. We explored these consequences by studying threespine stickleback in a set of bar‐built estuaries along the central California coast. In most years, heavy rains induce water flow strong enough to break through isolating sand bars, connecting streams to the ocean. New sand bars typically re‐form within a few weeks or months, thereby re‐isolating populations within the estuaries. These breaching events cause severe and often extremely rapid changes in abiotic and biotic conditions, including shifts in predator abundance. We investigated whether this strong temporal environmental variation can maintain within‐population variation while eroding adaptive divergence among populations that would be caused by spatial variation in selection. We used neutral genetic markers to explore population structure and then analysed how stickleback armor traits, the associated genes Eda and Pitx1 and elemental composition (%P) varies within and among populations. Despite strong gene flow, we detected evidence for divergence in stickleback defensive traits and Eda genotypes associated with predation regime. However, this among‐population variation was lower than that observed among other stickleback populations exposed to divergent predator regimes. In addition, within‐population variation was very high as compared to populations from environmentally stable locations. Elemental composition was strongly associated with armor traits, Eda genotype and the presence of predators, thus suggesting that spatiotemporal variation in armor traits generates corresponding variation in elemental phenotypes. We conclude that gene flow, and especially temporal environmental variation, can maintain high levels of within‐population variation while reducing, but not eliminating, among‐population variation driven by spatial environmental variation.     

Long‐term panmixia in a cosmopolitan Indo‐Pacific coral reef fish and a nebulous genetic boundary with its broadly sympatric sister species

Phylogeographical studies have shown that some shallow‐water marine organisms, such as certain coral reef fishes, lack spatial population structure at oceanic scales, despite vast distances of pelagic habitat between reefs and other dispersal barriers. However, whether these dispersive widespread taxa constitute long‐term panmictic populations across their species ranges remains unknown. Conventional phylogeographical inferences frequently fail to distinguish between long‐term panmixia and metapopulations connected by gene flow. Moreover, marine organisms have notoriously large effective population sizes that confound population structure detection. Therefore, at what spatial scale marine populations experience independent evolutionary trajectories and ultimately species divergence is still unclear. Here, we present a phylogeographical study of a cosmopolitan Indo‐Pacific coral reef fish Naso hexacanthus and its sister species Naso caesius, using two mtDNA and two nDNA markers. The purpose of this study was two‐fold: first, to test for broad‐scale panmixia in N. hexacanthus by fitting the data to various phylogeographical models within a Bayesian statistical framework, and second, to explore patterns of genetic divergence between the two broadly sympatric species. We report that N. hexacanthus shows little population structure across the Indo‐Pacific and a range‐wide, long‐term panmictic population model best fit the data. Hence, this species presently comprises a single evolutionary unit across much of the tropical Indian and Pacific Oceans. Naso hexacanthus and N. caesius were not reciprocally monophyletic in the mtDNA markers but showed varying degrees of population level divergence in the two nuclear introns. Overall, patterns are consistent with secondary introgression following a period of isolation, which may be attributed to oceanographic conditions of the mid to late Pleistocene, when these two species appear to have diverged.     

Genetic depletion at adaptive but not neutral loci in an endangered bird species

Many endangered species suffer from the loss of genetic diversity, but some populations may be able to thrive even if genetically depleted. To investigate the underlying genetic processes of population bottlenecks, we apply an innovative approach for assessing genetic diversity in the last known population of the endangered Pale‐headed Brushfinch (Atlapetes pallidiceps) in Ecuador. First, we measure genetic diversity at eleven neutral microsatellite loci and adaptive SNP variation in five Toll‐like receptor (TLR) immune system genes. Bottleneck tests confirm genetic drift as the main force shaping genetic diversity in this species and indicate a 99 % reduction in population size dating back several hundred years. Second, we compare contemporary microsatellite diversity with historic museum samples of A. pallidiceps, finding no change in genetic diversity. Third, we compare genetic diversity in the Pale‐headed Brushfinch with two co‐occurring‐related brushfinch species (Atlapetes latinuchus, Buarremon torquatus), finding a reduction of up to 91% diversity in the immune system genes but not in microsatellites. High TLR diversity is linked to decreased survival probabilities in A. pallidiceps. Low TLR diversity is thus probably an adaptation to the specific selection regime within its currently very restricted distribution (approximately 200 ha), but could severely restrict the adaptive potential of the species in the long run. Our study illustrates the importance of investigating both neutral and adaptive markers to assess the effect of population bottlenecks and for recommending specific management plans in endangered species.     

Spatial population genetic structure of a bacterial parasite in close coevolution with its host

Knowledge of a species’ population genetic structure can provide insight into fundamental ecological and evolutionary processes including gene flow, genetic drift and adaptive evolution. Such inference is of particular importance for parasites, as an understanding of their population structure can illuminate epidemiological and coevolutionary dynamics. Here, we describe the population genetic structure of the bacterium Pasteuria ramosa, a parasite that infects planktonic crustaceans of the genus Daphnia. This system has become a model for investigations of host–parasite interactions and represents an example of coevolution via negative frequency‐dependent selection (aka “Red Queen” dynamics). To sample P. ramosa, we experimentally infected a panel of Daphnia hosts with natural spore banks from the sediments of 25 ponds throughout much of the species range in Europe and western Asia. Using 12 polymorphic variable number tandem repeat loci (VNTR loci), we identified substantial genetic diversity, both within and among localities, that was structured geographically among ponds. Genetic diversity was also structured among host genotypes within ponds, although this pattern varied by locality, with P. ramosa at some localities partitioned into distinct host‐specific lineages, and other localities where recombination had shuffled genetic variation among different infection phenotypes. Across the sample range, there was a pattern of isolation by distance, and principal components analysis coupled with Procrustes rotation identified congruence between patterns of genetic variation and geography. Our findings support the hypothesis that Pasteuria is an endemic parasite coevolving closely with its host. These results provide important context for previous studies of this model system and inform hypotheses for future research.     

Population history in Arabidopsis halleri using multilocus analysis

A. halleri is a psuedometallophyte with a patchy distribution in Europe and is often spread by human activity. To determine the population history and whether this history is consistent with potential human effects, we surveyed nucleotide variation using 24 loci from 12 individuals in a large A. halleri population. The means of total and silent nucleotide variation (θ_W) are within the range expected for the species. The population genetic neutrality tests Tajima’s D and Wall’s B had significant composite results rejecting panmixia, and Approximate Bayesian Computation analysis revealed that a subdivision model better explained the variation than the standard neutral model, refugia (or admixture), bottleneck or change of population size models. A categorical regression analysis further supports the subdivision model, and under the subdivision model, the neutrality tests are no longer significant. The best support was for two source populations, a situation consistent with the mixing of two populations possibly mediated by human activity. This scenario might limit the genetic diversity and adaptive potential of the population. The non‐neutral population variation described here should be considered in bioinformatic searches for adaptation.     

On the roles of landscape heterogeneity and environmental variation in determining population genomic structure in a dendritic system

Dispersal and natural selection are key evolutionary processes shaping the distribution of phenotypic and genetic diversity. For species inhabiting complex spatial environments however, it is unclear how the balance between gene flow and selection may be influenced by landscape heterogeneity and environmental variation. Here, we evaluated the effects of dendritic landscape structure and the selective forces of hydroclimatic variation on population genomic parameters for the Murray River rainbowfish, Melanotaenia fluviatilis across the Murray–Darling Basin, Australia. We genotyped 249 rainbowfish at 17,503 high‐quality SNP loci and integrated these with models of network connectivity and high‐resolution environmental data within a riverscape genomics framework. We tested competing models of gene flow before using multivariate genotype–environment association (GEA) analysis to test for signals of adaptive divergence associated with hydroclimatic variation. Patterns of neutral genetic variation were consistent with expectations based on the stream hierarchy model and M. fluviatilis’ moderate dispersal ability. Models incorporating dendritic network structure suggested that landscape heterogeneity is a more important factor determining connectivity and gene flow than waterway distance. Extending these results, we also introduce a novel approach to controlling for the unique effects of dendritic network structure in GEA analyses of populations of aquatic species. We identified 146 candidate loci potentially underlying a polygenic adaptive response to seasonal fluctuations in stream flow and variation in the relative timing of temperature and precipitation extremes. Our findings underscore an emerging predominant role for seasonal variation in hydroclimatic conditions driving local adaptation and are relevant for informing proactive conservation management.     

Global lack of flyway structure in a cosmopolitan bird revealed by a genome wide survey of single nucleotide polymorphisms

Knowledge about population structure and connectivity of waterfowl species, especially mallards (Anas platyrhynchos), is a priority because of recent outbreaks of avian influenza. Ringing studies that trace large‐scale movement patterns have to date been unable to detect clearly delineated mallard populations. We employed 363 single nucleotide polymorphism markers in combination with population genetics and phylogeographical approaches to conduct a population genomic test of panmixia in 801 mallards from 45 locations worldwide. Basic population genetic and phylogenetic methods suggest no or very little population structure on continental scales. Nor could individual‐based structuring algorithms discern geographical structuring. Model‐based coalescent analyses for testing models of population structure pointed to strong genetic connectivity among the world's mallard population. These diverse approaches all support the conclusion that there is a lack of clear population structure, suggesting that the world's mallards, perhaps with minor exceptions, form a single large, mainly interbreeding population.     

Flyway structure in the circumpolar greater white‐fronted goose

Dispersal and migratory behavior are influential factors in determining how genetic diversity is distributed across the landscape. In migratory species, genetic structure can be promoted via several mechanisms including fidelity to distinct migratory routes. Particularly within North America, waterfowl management units have been delineated according to distinct longitudinal migratory flyways supported by banding data and other direct evidence. The greater white‐fronted goose (Anser albifrons) is a migratory waterfowl species with a largely circumpolar distribution consisting of up to six subspecies roughly corresponding to phenotypic variation. We examined the rangewide population genetic structure of greater white‐fronted geese using mtDNA control region sequence data and microsatellite loci from 23 locales across North America and Eurasia. We found significant differentiation in mtDNA between sampling locales with flyway delineation explaining a significant portion of the observed genetic variation (~12%). This is concordant with band recovery data which shows little interflyway or intercontinental movements. However, microsatellite loci revealed little genetic structure suggesting a panmictic population across most of the Arctic. As with many high‐latitude species, Beringia appears to have played a role in the diversification of this species. A common Beringian origin of North America and Asian populations and a recent divergence could at least partly explain the general lack of structure at nuclear markers. Further, our results do not provide strong support for the various taxonomic proposals for this species except for supporting the distinctness of two isolated breeding populations within Cook Inlet, Alaska (A. a. elgasi) and Greenland (A. a. flavirostris), consistent with their subspecies status.     

Adaptive differentiation in floral traits in the presence of high gene flow in scarlet gilia (Ipomopsis aggregata)

Plant–pollinator interactions are thought to be major drivers of floral trait diversity. However, the relative importance of divergent pollinator‐mediated selection vs. neutral processes in floral character evolution has rarely been explored. We tested for adaptive floral trait evolution by comparing differentiation at neutral genetic loci to differentiation at quantitative floral traits in a putative Ipomopsis aggregata hybrid zone. Typical I. aggregata subsp. candida displays slender white tubular flowers that are typical of flowers pollinated by hawkmoths, and subsp. collina displays robust red tubular flowers typical of flowers pollinated by hummingbirds; yet, hybrid flower morphs are abundant across the East Slope of the Colorado Rockies. We estimated genetic differentiation (F_ST) for nuclear and chloroplast microsatellite loci and used a half‐sib design to calculate quantitative trait divergence (Q_ST) from collection sites across the morphological hybrid zone. We found little evidence for population structure and estimated mean F_ST to be 0.032. Q_ST values for several floral traits including corolla tube length and width, colour, and nectar volume were large and significantly greater than mean F_ST. We performed multivariate comparisons of neutral loci to genetic correlations within and between populations and found a strong signal for divergent selection, suggesting that specific combinations of floral display and reward traits may be the targets of selection. Our results show little support for historical subspecies categories, yet floral traits are more diverged than expected due to drift alone. Non‐neutral divergence for multivariate quantitative traits suggests that selection by pollinators is maintaining a correlation between display and reward traits.     

Regional environmental pressure influences population differentiation in turbot (Scophthalmus maximus)

Unravelling the factors shaping the genetic structure of mobile marine species is challenging due to the high potential for gene flow. However, genetic inference can be greatly enhanced by increasing the genomic, geographical or environmental resolution of population genetic studies. Here, we investigated the population structure of turbot (Scophthalmus maximus) by screening 17 random and gene‐linked markers in 999 individuals at 290 geographical locations throughout the northeast Atlantic Ocean. A seascape genetics approach with the inclusion of high‐resolution oceanographical data was used to quantify the association of genetic variation with spatial, temporal and environmental parameters. Neutral loci identified three subgroups: an Atlantic group, a Baltic Sea group and one on the Irish Shelf. The inclusion of loci putatively under selection suggested an additional break in the North Sea, subdividing southern from northern Atlantic individuals. Environmental and spatial seascape variables correlated marginally with neutral genetic variation, but explained significant proportions (respectively, 8.7% and 10.3%) of adaptive genetic variation. Environmental variables associated with outlier allele frequencies included salinity, temperature, bottom shear stress, dissolved oxygen concentration and depth of the pycnocline. Furthermore, levels of explained adaptive genetic variation differed markedly between basins (3% vs. 12% in the North and Baltic Sea, respectively). We suggest that stable environmental selection pressure contributes to relatively strong local adaptation in the Baltic Sea. Our seascape genetic approach using a large number of sampling locations and associated oceanographical data proved useful for the identification of population units as the basis of management decisions.     

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.     

Climate variables explain neutral and adaptive variation within salmonid metapopulations: the importance of replication in landscape genetics

Understanding how environmental variation influences population genetic structure is important for conservation management because it can reveal how human stressors influence population connectivity, genetic diversity and persistence. We used riverscape genetics modelling to assess whether climatic and habitat variables were related to neutral and adaptive patterns of genetic differentiation (population‐specific and pairwise F_ST) within five metapopulations (79 populations, 4583 individuals) of steelhead trout (Oncorhynchus mykiss) in the Columbia River Basin, USA. Using 151 putatively neutral and 29 candidate adaptive SNP loci, we found that climate‐related variables (winter precipitation, summer maximum temperature, winter highest 5% flow events and summer mean flow) best explained neutral and adaptive patterns of genetic differentiation within metapopulations, suggesting that climatic variation likely influences both demography (neutral variation) and local adaptation (adaptive variation). However, we did not observe consistent relationships between climate variables and F_ST across all metapopulations, underscoring the need for replication when extrapolating results from one scale to another (e.g. basin‐wide to the metapopulation scale). Sensitivity analysis (leave‐one‐population‐out) revealed consistent relationships between climate variables and F_ST within three metapopulations; however, these patterns were not consistent in two metapopulations likely due to small sample sizes (N = 10). These results provide correlative evidence that climatic variation has shaped the genetic structure of steelhead populations and highlight the need for replication and sensitivity analyses in land and riverscape genetics.     

Genetic demography at the leading edge of the distribution of a rabies virus vector

The common vampire bat, Desmodus rotundus, ranges from South America into northern Mexico in North America. This sanguivorous species of bat feeds primarily on medium to large‐sized mammals and is known to rely on livestock as primary prey. Each year, there are hotspot areas of D. rotundus‐specific rabies virus outbreaks that lead to the deaths of livestock and economic losses. Based on incidental captures in our study area, which is an area of high cattle mortality from D. rotundus transmitted rabies, it appears that D. rotundus are being caught regularly in areas and elevations where they previously were thought to be uncommon. Our goal was to investigate demographic processes and genetic diversity at the north eastern edge of the range of D. rotundus in Mexico. We generated control region sequences (441 bp) and 12‐locus microsatellite genotypes for 602 individuals of D. rotundus. These data were analyzed using network analyses, Bayesian clustering approaches, and standard population genetic statistical analyses. Our results demonstrate panmixia across our sampling area with low genetic diversity, low population differentiation, loss of intermediate frequency alleles at microsatellite loci, and very low mtDNA haplotype diversity with all haplotypes being very closely related. Our study also revealed strong signals of population expansion. These results follow predictions from the leading‐edge model of expanding populations and supports conclusions from another study that climate change may allow this species to find suitable habitat within the U.S. border.     

Single nucleotide polymorphism discovery in albacore and Atlantic bluefin tuna provides insights into worldwide population structure

The optimal management of the commercially important, but mostly over‐exploited, pelagic tunas, albacore (Thunnus alalunga Bonn., 1788) and Atlantic bluefin tuna (BFT; Thunnus thynnus L., 1758), requires a better understanding of population structure than has been provided by previous molecular methods. Despite numerous studies of both species, their population structures remain controversial. This study reports the development of single nucleotide polymorphisms (SNPs) in albacore and BFT and the application of these SNPs to survey genetic variability across the geographic ranges of these tunas. A total of 616 SNPs were discovered in 35 albacore tuna by comparing sequences of 54 nuclear DNA fragments. A panel of 53 SNPs yielded F_ST values ranging from 0.0 to 0.050 between samples after genotyping 460 albacore collected throughout the distribution of this species. No significant heterogeneity was detected within oceans, but between‐ocean comparisons (Atlantic, Pacific and Indian oceans along with Mediterranean Sea) were significant. Additionally, a 17‐SNP panel was developed in Atlantic BFT by cross‐species amplification in 107 fish. This limited number of SNPs discriminated between samples from the two major spawning areas of Atlantic BFT (F_ST = 0.116). The SNP markers developed in this study can be used to genotype large numbers of fish without the need for standardizing alleles among laboratories.     

Breaking RAD: an evaluation of the utility of restriction site‐associated DNA sequencing for genome scans of adaptation

Understanding how and why populations evolve is of fundamental importance to molecular ecology. Restriction site‐associated DNA sequencing (RADseq), a popular reduced representation method, has ushered in a new era of genome‐scale research for assessing population structure, hybridization, demographic history, phylogeography and migration. RADseq has also been widely used to conduct genome scans to detect loci involved in adaptive divergence among natural populations. Here, we examine the capacity of those RADseq‐based genome scan studies to detect loci involved in local adaptation. To understand what proportion of the genome is missed by RADseq studies, we developed a simple model using different numbers of RAD‐tags, genome sizes and extents of linkage disequilibrium (length of haplotype blocks). Under the best‐case modelling scenario, we found that RADseq using six‐ or eight‐base pair cutting restriction enzymes would fail to sample many regions of the genome, especially for species with short linkage disequilibrium. We then surveyed recent studies that have used RADseq for genome scans and found that the median density of markers across these studies was 4.08 RAD‐tag markers per megabase (one marker per 245 kb). The length of linkage disequilibrium for many species is one to three orders of magnitude less than density of the typical recent RADseq study. Thus, we conclude that genome scans based on RADseq data alone, while useful for studies of neutral genetic variation and genetic population structure, will likely miss many loci under selection in studies of local adaptation.