High Gene Flow Species Research Articles

Restriction Site-Associated DNA Sequencing Reveals Local Adaptation Despite High Levels of Gene Flow in Sardinella lemuru (Bleeker, 1853) Along the Northern Coast of Mindanao, Philippines

Stock identification and delineation are important in the management and conservation of marine resources. These were highlighted as priority research areas for Bali sardinella (Sardinella lemuru) which is among the most commercially important fishery resources in the Philippines. Previous studies have already assessed the stocks of S. lemuru between Northern Mindanao Region (NMR) and Northern Zamboanga Peninsula (NZP), yielding conflicting results. Phenotypic variation suggests distinct stocks between the two regions, while mitochondrial DNA did not detect evidence of genetic differentiation for this high gene flow species. This paper tested the hypothesis of regional structuring using genome-wide single nucleotide polymorphisms (SNPs) acquired through restriction site-associated DNA sequencing (RADseq). We examined patterns of population genomic structure using a full panel of 3,573 loci, which was then partitioned into a neutral panel of 3,348 loci and an outlier panel of 31 loci. Similar inferences were obtained from the full and neutral panels, which were contrary to the inferences from the outlier panel. While the full and neutral panels suggested a panmictic population (global FST ∼ 0, p > 0.05), the outlier panel revealed genetic differentiation between the two regions (global FST = 0.161, p = 0.001; FCT = 0.263, p < 0.05). This indicated that while gene flow is apparent, selective forces due to environmental heterogeneity between the two regions play a role in maintaining adaptive variation. Annotation of the outlier loci returned five genes that were mostly involved in organismal development. Meanwhile, three unannotated loci had allele frequencies that correlated with sea surface temperature. Overall, our results provided support for local adaptation despite high levels of gene flow in S. lemuru. Management therefore should not only focus on demographic parameters (e.g., stock size and catch volume), but also consider the preservation of adaptive variation.

Spatial gradients of introgressed ancestry reveal cryptic connectivity patterns in a high gene flow marine fish.

Assessing genetic connectivity among populations in high gene flow species is sometimes insufficient to evaluate demographic connectivity. Genetic differentiation quickly becomes zero as soon as a few dozen migrants are exchanged per generation. This provides little information to determine whether migration can ensure demographic coupling. The resulting difficulties in delineating conservation units for the management of commercially exploited marine fish species are well illustrated in the case of the European sea bass (Dicentrarchus labrax). Previous attempts to assess connectivity patterns in the northeast Atlantic have been hampered by a lack of spatial genetic structure. In contrast, mark-recapture data suggested low migration rates between regional spawning areas. Here, we show how a spatial gradient of introgressed Mediterranean ancestry across the northeast Atlantic reflects cryptic patterns of genetic and demographic connectivity. Using a 1K SNP chip data set in 827 individuals sampled from Portugal to the North Sea, we found null overall genetic differentiation across the northeast Atlantic. We however detected a subtle latitudinal admixture gradient originating at the edge of the contact zone with the Mediterranean sea bass lineage. Two significant breaks in the ancestry gradient at the tip of Galicia and northern Brittany indicated barriers to effective dispersal between demographically distinct units. Moreover, a northward expansion signal in Irish and North Seas was revealed by the surfing of rare Mediterranean alleles at the edge of the species range. Our results show that introgressed ancestry gradients offer a powerful alternative to assess genetic and demographic connectivity when the neutral migration-drift balance is not informative.

Eco-evolutionary community turnover following environmental change.

Co‐occurring species often differ in intraspecific genetic diversity, which in turn can affect adaptation in response to environmental change. Specifically, the simultaneous evolutionary responses of co‐occurring species to temporal environmental change may influence community dynamics. Local adaptation along environmental gradients combined with gene flow can enhance genetic diversity of traits within populations. Quantitative genetic theory shows that having greater gene flow results in (a) lower equilibrium population size due to maladaptive immigrant genotypes (migration load), but (b) faster adaptation to changing environments. Here, I build off this theory to study community dynamics of locally adapted species in response to temporal environmental changes akin to warming temperatures. Although an abrupt environmental change leaves all species initially maladapted, high gene flow species subsequently adapt faster due to greater genetic diversity. As a result, species can transiently reverse their relative abundances, but sometimes only after long lag periods. If constant temporal environmental change is applied, the community exhibits a shift toward stable dominance by species with intermediate gene flow. Notably, fast‐adapting high gene flow species can increase in absolute abundance under environmental change (although often only for a transient period) because the change suppresses superior competitors with lower gene flow. This eco‐evolutionary competitive release stabilizes ecosystem function. The eco‐evolutionary community turnover studied here parallels the purely ecological successional dynamics following disturbances. My results demonstrate how interspecific variation in life history can have far‐reaching impacts on eco‐evolutionary community response to environmental change.

Genome architecture enables local adaptation of Atlantic cod despite high connectivity.

Adaptation to local conditions is a fundamental process in evolution; however, mechanisms maintaining local adaptation despite high gene flow are still poorly understood. Marine ecosystems provide a wide array of diverse habitats that frequently promote ecological adaptation even in species characterized by strong levels of gene flow. As one example, populations of the marine fish Atlantic cod (Gadus morhua) are highly connected due to immense dispersal capabilities but nevertheless show local adaptation in several key traits. By combining population genomic analyses based on 12K single nucleotide polymorphisms with larval dispersal patterns inferred using a biophysical ocean model, we show that Atlantic cod individuals residing in sheltered estuarine habitats of Scandinavian fjords mainly belong to offshore oceanic populations with considerable connectivity between these diverse ecosystems. Nevertheless, we also find evidence for discrete fjord populations that are genetically differentiated from offshore populations, indicative of local adaptation, the degree of which appears to be influenced by connectivity. Analyses of the genomic architecture reveal a significant overrepresentation of a large ~5Mb chromosomal rearrangement in fjord cod, previously proposed to comprise genes critical for the survival at low salinities. This suggests that despite considerable connectivity with offshore populations, local adaptation to fjord environments may be enabled by suppression of recombination in the rearranged region. Our study provides new insights into the potential of local adaptation in high gene flow species within fine geographical scales and highlights the importance of genome architecture in analyses of ecological adaptation.

Temporal genetic change in North American Pacific oyster populations suggests caution in seascape genetics analyses of high gene-flow species

Temporal genetic change in North American Pacific oyster populations suggests caution in seascape genetics analyses of high gene-flow species

A Parallel Population Genomic and Hydrodynamic Approach to Fishery Management of Highly-Dispersive Marine Invertebrates: The Case of the Fijian Black-Lip Pearl Oyster Pinctada margaritifera.

Fishery management and conservation of marine species increasingly relies on genetic data to delineate biologically relevant stock boundaries. Unfortunately for high gene flow species which may display low, but statistically significant population structure, there is no clear consensus on the level of differentiation required to resolve distinct stocks. The use of fine-scale neutral and adaptive variation, considered together with environmental data can offer additional insights to this problem. Genome-wide genetic data (4,123 SNPs), together with an independent hydrodynamic particle dispersal model were used to inform farm and fishery management in the Fijian black-lip pearl oyster Pinctada margaritifera, where comprehensive fishery management is lacking, and the sustainability of exploitation uncertain. Weak fine-scale patterns of population structure were detected, indicative of broad-scale panmixia among wild oysters, while a hatchery-sourced farmed population exhibited a higher degree of genetic divergence (Fst = 0.0850–0.102). This hatchery-produced population had also experienced a bottleneck (NeLD = 5.1; 95% C.I. = [5.1–5.3]); compared to infinite NeLD estimates for all wild oysters. Simulation of larval transport pathways confirmed the existence of broad-scale mixture by surface ocean currents, correlating well with fine-scale patterns of population structuring. Fst outlier tests failed to detect large numbers of loci supportive of selection, with 2–5 directional outlier SNPs identified (average Fst = 0.116). The lack of biologically significant population genetic structure, absence of evidence for local adaptation and larval dispersal simulation, all indicate the existence of a single genetic stock of P. margaritifera in the Fiji Islands. This approach using independent genomic and oceanographic tools has allowed fundamental insights into stock structure in this species, with transferability to other highly-dispersive marine taxa for their conservation and management.

Methodological assessment of 2b-RAD genotyping technique for population structure inferences in yellowfin tuna (Thunnus albacares)

Global population genetic structure of yellowfin tuna (Thunnus albacares) is still poorly understood despite its relevance for the tuna fishery industry. Low levels of genetic differentiation among oceans speak in favour of the existence of a single panmictic population worldwide of this highly migratory fish. However, recent studies indicated genetic structuring at a much smaller geographic scales than previously considered, pointing out that YFT population genetic structure has not been properly assessed so far. In this study, we demonstrated for the first time, the utility of 2b-RAD genotyping technique for investigating population genetic diversity and differentiation in high gene-flow species. Running de novo pipeline in Stacks, a total of 6772 high-quality genome-wide SNPs were identified across Atlantic, Indian and Pacific population samples representing all major distribution areas. Preliminary analyses showed shallow but significant population structure among oceans (FST=0.0273; P-value<0.01). Discriminant Analysis of Principal Components endorsed the presence of genetically discrete yellowfin tuna populations among three oceanic pools. Although such evidence needs to be corroborated by increasing sample size, these results showed the efficiency of this genotyping technique in assessing genetic divergence in a marine fish with high dispersal potential.

Population differentiation determined from putative neutral and divergent adaptive genetic markers in Eulachon (Thaleichthys pacificus, Osmeridae), an anadromous Pacific smelt.

Twelve eulachon (Thaleichthys pacificus, Osmeridae) populations ranging from Cook Inlet, Alaska and along the west coast of North America to the Columbia River were examined by restriction-site-associated DNA (RAD) sequencing to elucidate patterns of neutral and adaptive variation in this high geneflow species. A total of 4104 single-nucleotide polymorphisms (SNPs) were discovered across the genome, with 193 putatively adaptive SNPs as determined by F(ST) outlier tests. Estimates of population structure in eulachon with the putatively adaptive SNPs were similar, but provided greater resolution of stocks compared with a putatively neutral panel of 3911 SNPs or previous estimates with 14 microsatellites. A cline of increasing measures of genetic diversity from south to north was found in the adaptive panel, but not in the neutral markers (SNPs or microsatellites). This may indicate divergent selective pressures in differing freshwater and marine environments between regional eulachon populations and that these adaptive diversity patterns not seen with neutral markers could be a consideration when determining genetic boundaries for conservation purposes. Estimates of effective population size (N(e)) were similar with the neutral SNP panel and microsatellites and may be utilized to monitor population status for eulachon where census sizes are difficult to obtain. Greater differentiation with the panel of putatively adaptive SNPs provided higher individual assignment accuracy compared to the neutral panel or microsatellites for stock identification purposes. This study presents the first SNPs that have been developed for eulachon, and analyses with these markers highlighted the importance of integrating genome-wide neutral and adaptive genetic variation for the applications of conservation and management.

A 40-year-old divided highway does not prevent gene flow in the alpine newt Ichthyosaura alpestris

Roads are of major concern in conservation biology, as they are known to restrict animal movements through landscape fragmentation, and may therefore impact genetic patterns in native terrestrial organisms. We assessed the effect of two large-scale transportation infrastructures (LTIs), a 40-year-old highway and a 30 year-old high-speed railway, on the spatial genetic structure of the alpine newt Ichthyosaura alpestris, a highly nomadic amphibian. Genetic data were gathered following a targeted individual-based sampling scheme and analysed using both overlay and correlative methods. While simulations suggested that the highway may be old enough for a significant barrier effect to be detected, LTIs were never detected as barriers to gene flow: inferred genetic boundaries rather coincided with transition zones between major landscape entities. Furthermore, spatial principal component analysis, a method designed to reveal cryptic genetic spatial patterns in high gene flow species, counter-intuitively suggested that the highway may act as a potential dispersal corridor in low-quality habitats, thus challenging traditional hypotheses on road impacts in amphibians. Our study showed that considering local interactions between species, infrastructures and landscape-specific characteristics is essential for better understanding the potential impacts of roads on movement patterns in terrestrial organisms.

Coastal pollution limits pelagic larval dispersal

The ecological impact of large coastal human populations on marine ecosystems remains relatively unknown. Here, we examine the population structure of Patiria miniata, the bat star, and correlate genetic distances with a model based on flow rates and proximity to P. miniata populations for the four major stormwater runoff and wastewater effluent sources of the Southern California Bight. We show that overall genetic connectivity is high (F(ST)~0.005); however, multivariate analyses show that genetic structure is highly correlated with anthropogenic inputs. The best models included both stormwater and wastewater variables and explained between 26.55 and 93.69% of the observed structure. Additionally, regressions between allelic richness and distance to sources show that populations near anthropogenic pollution have reduced genetic diversity. Our results indicate that anthropogenic runoff and effluent are acting as barriers to larval dispersal, effectively isolating a high gene flow species that is virtually free of direct human impact.

A Method for Detecting Population Genetic Structure in Diverse, High Gene-Flow Species

Detecting small amounts of genetic subdivision across geographic space remains a persistent challenge. Often a failure to detect genetic structure is mistaken for evidence of panmixia, when more powerful statistical tests may uncover evidence for subtle geographic differentiation. Such slight subdivision can be demographically and evolutionarily important as well as being critical for management decisions. We introduce here a method, called spatial analysis of shared alleles (SAShA), that detects geographically restricted alleles by comparing the spatial arrangement of allelic co-occurrences with the expectation under panmixia. The approach is allele-based and spatially explicit, eliminating the loss of statistical power that can occur with user-defined populations and statistical averaging within populations. Using simulated data sets generated under a stepping-stone model of gene flow, we show that this method outperforms spatial autocorrelation (SA) and Phi(ST) under common real-world conditions: at relatively high migration rates when diversity is moderate or high, especially when sampling is poor. We then use this method to show clear differences in the genetic patterns of 2 nearshore Pacific mollusks, Tegula funebralis (= Chlorostoma funebralis) and Katharina tunicata, whose overall patterns of within-species differentiation are similar according to traditional population genetics analyses. SAShA meaningfully complements Phi(ST)/F(ST), SA, and other existing geographic genetic analyses and is especially appropriate for evaluating species with high gene flow and subtle genetic differentiation.

Ocean currents help explain population genetic structure.

Management and conservation can be greatly informed by considering explicitly how environmental factors influence population genetic structure. Using simulated larval dispersal estimates based on ocean current observations, we demonstrate how explicit consideration of frequency of exchange of larvae among sites via ocean advection can fundamentally change the interpretation of empirical population genetic structuring as compared with conventional spatial genetic analyses. Both frequency of larval exchange and empirical genetic difference were uncorrelated with Euclidean distance between sites. When transformed into relative oceanographic distances and integrated into a genetic isolation-by-distance framework, however, the frequency of larval exchange explained nearly 50 per cent of the variance in empirical genetic differences among sites over scales of tens of kilometres. Explanatory power was strongest when we considered effects of multiple generations of larval dispersal via intermediary locations on the long-term probability of exchange between sites. Our results uncover meaningful spatial patterning to population genetic structuring that corresponds with ocean circulation. This study advances our ability to interpret population structure from complex genetic data characteristic of high gene flow species, validates recent advances in oceanographic approaches for assessing larval dispersal and represents a novel approach to characterize population connectivity at small spatial scales germane to conservation and fisheries management.

Population genomics of marine fishes: identifying adaptive variation in space and time

Studies of adaptive evolution have experienced a recent revival in population genetics of natural populations and there is currently much focus on identifying genomic signatures of selection in space and time. Insights into local adaptation, adaptive response to global change and evolutionary consequences of selective harvesting can be generated through population genomics studies, allowing the separation of the effects invoked by neutral processes (drift-migration) from those due to selection. Such knowledge is important not only for improving our basic understanding of natural as well as human-induced evolutionary processes, but also for predicting future trajectories of biodiversity and for setting conservation priorities. Marine fishes possess a number of features rendering them well suited for providing general insights into adaptive genomic evolution in natural populations. These include well-described population structures, substantial and rapidly developing genomic resources and abundant archived samples enabling temporal studies. Furthermore, superior possibilities for conducting large-scale experiments under controlled conditions, due to the economic resources provided by the large and growing aquaculture industry, hold great promise for utilizing recent technological developments. Here, we review achievements in marine fish genomics to date and highlight potential avenues for future research, which will provide both general insights into evolution in high gene flow species, as well as specific knowledge which can lead to improved management of marine organisms.

Remnant populations of the regal fritillary (Speyeria idalia) in Pennsylvania: Local genetic structure in a high gene flow species

The Regal Fritillary butterfly, Speyeria idalia (Drury) (Lepidoptera: Nymphalidae), has been described as a high gene flow species. Supporting this assertion, previous studies in the Great Plains, where it is still relatively widespread, have found evidence of gene flow across hundreds of kilometers. Using mitochondrial and microsatellite loci, we examined the spatial genetic structure of a very isolated Pennsylvania population of these butterflies that occupies three separate meadows located within ten kilometers of each other. We found restricted gene flow and a distinct structure, with each meadow having a unique genetic signature. Our findings indicate that even a species that normally exhibits high gene flow may show fine-scale genetic subdivision in areas where populations have been largely extirpated.

Landscape scale genetic effects of habitat fragmentation on a high gene flow species: Speyeria idalia (Nymphalidae).

Detection of the genetic effects of recent habitat fragmentation in natural populations can be a difficult task, especially for high gene flow species. Previous analyses of mitochondrial DNA data from across the current range of Speyeria idalia indicated that the species exhibited high levels of gene flow among populations, with the exception of an isolated population in the eastern portion of its range. However, some populations are found on isolated habitat patches, which were recently separated from one another by large expanses of uninhabitable terrain, in the form of row crop agriculture. The goal of this study was to compare levels of genetic differentiation and diversity among populations found in relatively continuous habitat to populations in both recently and historically isolated habitat. Four microsatellite loci were used to genotype over 300 individuals from five populations in continuous habitat, five populations in recently fragmented habitat, and one historically isolated population. Results from the historically isolated population were concordant with previous analyses and suggest significant differentiation. Also, microsatellite data were consistent with the genetic effects of habitat fragmentation for the recently isolated populations, in the form of increased differentiation and decreased genetic diversity when compared to nonfragmented populations. These results suggest that given the appropriate control populations, microsatellite markers can be used to detect the effects of recent habitat fragmentation in natural populations, even at a large geographical scale in high gene flow species.

Separating the wheat from the chaff: patterns of genetic differentiation in high gene flow species

In many marine species, high levels of gene flow ensure that the genetic signal from population differentiation is weak. As a consequence, various errors associated with estimating population genetic parameters that might normally be safely Ignored assume a relatively greater importance. This fact has important implications for the use of genetic data to address two common questions in fishery conservation and management: (1) How many stocks of a given species are there? and (2) How much gene flow occurs among stocks? This article discusses strategies to maximize the signal:nolse ratio in genetic studies of marine species and suggests a quantitative method to correct for bias due to a common sampling problem. For many marine species, however, genetic methods alone cannot fully resolve these key management questions because the amount of migration necessary to eliminate most genetic evidence of stock structure (only a handful of Individuals per generation) will generally be Inconsequential as a force for rebuilding depleted populations on a time scale of interest to humans. These limitations emphasize the importance of understanding the biology and life history of the target speciesfirst, to guide design of the sampling program, and second, so that additional information can be used to supplement Indirect estimates of migration rates based on genetic data.

Source regions for recruitment of Calanus finmarchicus to Georges Bank: evidence from molecular population genetic analysis of mtDNA

Molecular population genetic analysis has provided evidence that the copepod, Calanus finmarchicus, of the Labrador Current, Gulf of St Lawrence, Scotian Shelf, Gulf of Maine, and Georges Bank constitute a single, interbreeding population. The DNA sequence of a 350 base pair portion of the mitochondrial large subunit (16S) ribosomal RNA (rRNA) gene was determined for a total of 72 individuals collected in 1992, and 110 individuals collected in 1993 from these regions. There was significant heterogeneity in haplotype frequencies among the samples collected in 1992, but this heterogeneity did not resolve into regional patterns. The only regional differences seen were between pooled samples of the western N. Atlantic and those of the Norwegian Sea. There were no significant differences in haplotype frequencies among the samples collected in 1993, and fewer haplotypes were observed in these samples. Intraspecific molecular variation was typical of other marine species: there were 29 haplotypes among the 182 individuals sequenced. The frequency distribution of the haplotypes was highly skewed: 128 individuals shared one haplotype and 19 individuals were unique. There were 24 variable sites among the 350 bases sequenced; estimated nucleotide diversity was 0.0042. The genetic character of C. finmarchicus populations in the western N. Atlantic was stable over time in that three of the haplotypes (including the most abundant) occurred in both 1992 and 1993. However, haplotype frequencies differed significantly between the two years. The lack of regional structure in the 1992 samples and the genetic homogeneity of samples collected in 1993 across the domain from the Labrador Current to the Gulf of St Lawrence to Georges Bank and the Gulf of Maine indicated that there is significant gene flow across this region. The persistent genetic pattern suggests that the Gulf of St Lawrence may be an important source region for recruitment of C. finmarchicus to Georges Bank. Determination of zooplankton dispersal patterns within high gene flow species will provide information that may not be determined by conventional oceanographic analyses.

A MULTISPECIES APPROACH TO THE ANALYSIS OF GENE FLOW IN MARINE SHORE FISHES.

Ten species of marine shore fishes with a wide range of life-history strategies were collected from four areas in southern California, U.S.A., and Baja California, Mexico, and examined for patterns of genetic differentiation. Multilocus D and FST values (based on 32-42 presumptive gene loci in each species) were both negatively correlated with estimated dispersal capability. These results were robust to variations in the number and type of loci used in the analysis and are compatible with the hypothesis that levels of genetic differentiation in these shore fishes are determined primarily by gene flow and genetic drift. There is no a priori reason to expect the observed correlation to result from natural selection or historical factors. The findings thus suggest that populations of these shore fishes are in at least a quasi-equilibrium with respect to migration, mutation, and genetic drift. Present data were also used to compare estimates of mNe obtained by three different methods. Estimates based on FST values calculated by the methods of Nei and Chesser (FST(N) ) and Weir and Cockerham (FST(W) ) were highly correlated, but FST(N) ≤ FST(W) for every species, leading to generally higher mNe estimates for Nei and Chesser's method. Estimates of mNe based on the frequency of private alleles (Slatkin, 1985a) were not as strongly correlated with dispersal capability as were FST and D values. A low incidence of private alleles in many species may be responsible for this relatively weak correlation and may limit the general usefulness of Slatkin's method. In spite of their sensitivity to natural selection, FST and D may be better indicators of relative gene flow levels for high gene flow species.