Inverse relationship between F and microsatellite polymorphism in the marine fish, walleye pollock (Theragra chalcogramma): implications for resolving weak population structure

Patterns of persistence: Genetic and behavioral population complexity of winter flounder amid population declines

Winter flounder Pseudopleuronectes americanus (Walbaum 1792) are a coastal flatfish species of economic and cultural importance that have dwindled to <15, % of their historic abundance in the southern New England/Mid-Atlantic region of the United States, with evidence indicating near-extirpation of certain local populations. This species exhibits intricate behaviors in spawning and migration that contribute to population complexity and resilience. These behaviors encompass full or partial philopatry to natal estuaries, the generation of multiple pulses of larval delivery, and partial migration. The patterns of genetic diversity within and among estuaries and cohorts presented here carry important implications in understanding the susceptibility to demographic shocks, even if the full extent of genetic diversity within and among winter flounder stocks on the US East Coast remains unresolved. Our findings reveal connectivity between estuaries in Long Island, New York, suggesting the potential for genetic rescue of depleted subpopulations. Family reconstruction and relatedness analysis indicate that split cohorts and migration contingents are not the result of genetically distinct lineages. We found no evidence for genetic structure separating these groups, and in some instances, we were able to detect closely related individuals that belonged to different migratory contingents or cohorts. Characterizing the spatial and behavioral organization of this species at the population level is crucial for comprehending its potential for recovery, not only in terms of biomass but also in reinstating the complex population structure that supports resilience. The search for generality in winter flounder spawning and migration behavior remains elusive, but perhaps the lack of generalities within this species is what has allowed it to persist in the face of decades of environmental and anthropogenic stressors.

EST-Microsatellite Types and Structural Scenarios in European Hake Fisheries

A fishery’s structure and connectivity are priors to its effective management. A successful description of such processes depends on both the sampling design and the choice of adequate genetic markers. EST markers are perfusing the studies of marine metapopulations and are believed to provide access to functional polymorphisms. However, the assumed adaptive role of outlier EST loci might not be generalizable. EST-microsatellites represent the upper polymorphic boundary in these regions because of their high mutation rate. We have subclassified the polymorphisms of EST-microsatellites to assess their structural contribution in the European hake, a paradigmatic and highly mobile marine species (HMMS). Because of the counterbalanced forces between directional markers (15%) and balanced markers (23%), the whole marker set offers the same structural situation as the one observed with neutral markers (62%), i.e., k = 2 gene pools. In contrast to outlier EST- microsatellites, neutral EST subsets allow one to measure crucial population phenomena for fisheries’ management. The high inter-population divergence of outlier EST-microsatellites is compatible with drifted post-selection genomic regions rather than with ongoing local selective pressures. The structural scenario in hake is explainable by a limited gene flow across the Almería-Oran Front (AOF) and by the within-basin IBD pattern of connectivity plus drift-related demographic events. This study highlights how polymorphic properties of EST-microsatellite types can be useful to address mutually excluding research tasks in fisheries, i.e., to address its evolutionary history (directional markers or FAPS: Fossil Adaptive Polymorphic Systems); to delineate management units (neutral markers or NAPS: Non Adaptive Polymorphic Systems); or to ensure sustainability (balanced markers or APS: Adaptive Polymorphic Systems).

High gene flow in polar cod (Boreogadus saida) from West-Svalbard and the Eurasian Basin

The current and projected environmental change of the Arctic Ocean contrasts sharply with the limited knowledge of its genetic biodiversity. Polar cod Boreogadus saida (Lepechin, 1774) is an abundant circumpolar marine fish and ecological key species. The central role of polar cod in the Arctic marine food web warrants a better understanding of its population structure and connectivity. In this study, the genetic population structure of 171 juveniles, collected from several fjords off West-Svalbard (Billefjorden, Hornsund and Kongsfjorden), the northern Sophia Basin and the Eurasian Basin of the Arctic Ocean, was analysed using nine DNA microsatellite loci. Genetic analyses indicated moderate to high genetic diversity, but absence of spatial population structure and isolation-by-distance, suggesting ongoing gene flow between the studied sampling regions. High levels of connectivity may be key for polar cod to maintain populations across wide spatial scales. The adaptive capacity of the species will be increasingly important to face challenges such as habitat fragmentation, ocean warming and changes in prey composition. In view of a limited understanding of the population dynamics and evolution of polar cod, a valuable next step to predict future developments should be an integrated biological evaluation, including population genomics, a life-history approach, and habitat and biophysical dispersal modelling.

Genotyping-by-Sequencing of Gossypium hirsutum Races and Cultivars Uncovers Novel Patterns of Genetic Relationships and Domestication Footprints

Determining the genetic rearrangement and domestication footprints in Gossypium hirsutum cultivars and primitive race genotypes are essential for effective gene conservation efforts and the development of advanced breeding molecular markers for marker-assisted breeding. In this study, 94 accessions representing the 7 primitive races of G hirsutum, along with 9 G hirsutum and 12 Gossypium barbadense cultivated accessions were evaluated. The genotyping-by-sequencing (GBS) approach was employed and 146 558 single nucleotide polymorphisms (SNP) were generated. Distinct SNP signatures were identified through the combination of selection scans and association analyses. Phylogenetic analyses were also conducted, and we concluded that the Latifolium, Richmondi, and Marie-Galante race accessions were more genetically related to the G hirsutum cultivars and tend to cluster together. Fifty-four outlier SNP loci were identified by selection-scan analysis, and 3 SNPs were located in genes related to the processes of plant responding to stress conditions and confirmed through further genome-wide signals of marker-phenotype association analysis, which indicate a clear selection signature for such trait. These results identified useful candidate gene locus for cotton breeding programs.

Inferring genetic connectivity in real populations, exemplified by coastal and oceanic Atlantic cod

Estimates of migration are important for understanding the dynamics of natural populations. A statistic known as F_ST is often used to measure levels of genetic differentiation among natural populations. Equations that translate F_ST into estimates of migration are based on “ideal” populations, which are subject to many simplifying assumptions compared with real populations. Therefore, theoretical estimates of migration might not be realistic. We modeled populations of Atlantic cod in the North Sea and the adjacent Skagerrak region to compare how migration is related to the complexities of real populations, and how actual migration compares with predictions based on theory. Results are intended to help apply population genetic theory to practical situations.

Genetic data are commonly used to estimate connectivity between putative populations, but translating them to demographic dispersal rates is complicated. Theoretical equations that infer a migration rate based on the genetic estimator F_ST, such as Wright’s equation, F_ST ≈ 1/(4N_em + 1), make assumptions that do not apply to most real populations. How complexities inherent to real populations affect migration was exemplified by Atlantic cod in the North Sea and Skagerrak and was examined within an age-structured model that incorporated genetic markers. Migration was determined under various scenarios by varying the number of simulated migrants until the mean simulated level of genetic differentiation matched a fixed level of genetic differentiation equal to empirical estimates. Parameters that decreased the N_e/N_t ratio (where N_e is the effective and N_t is the total population size), such as high fishing mortality and high fishing gear selectivity, increased the number of migrants required to achieve empirical levels of genetic differentiation. Higher maturity-at-age and lower selectivity increased N_e/N_t and decreased migration when genetic differentiation was fixed. Changes in natural mortality, fishing gear selectivity, and maturity-at-age within expected limits had a moderate effect on migration when genetic differentiation was held constant. Changes in population size had the greatest effect on the number of migrants to achieve fixed levels of F_ST, particularly when genetic differentiation was low, F_ST ≈ 10⁻³. Highly variable migration patterns, compared with constant migration, resulted in higher variance in genetic differentiation and higher extreme values. Results are compared with and provide insight into the use of theoretical equations to estimate migration among real populations.

Widespread gene flow between oceans in a pelagic seabird species complex

Global-scale gene flow is an important concern in conservation biology as it has the potential to either increase or decrease genetic diversity in species and populations. Although many studies focus on the gene flow between different populations of a single species, the potential for gene flow and introgression between species is understudied, particularly in seabirds. The only well-studied example of a mixed-species, hybridizing population of petrels exists on Round Island, in the Indian Ocean. Previous research assumed that Round Island represents a point of secondary contact between Atlantic (Pterodroma arminjoniana) and Pacific species (Pterodroma neglecta and Pterodroma heraldica). This study uses microsatellite genotyping and tracking data to address the possibility of between-species hybridization occurring outside the Indian Ocean. Dispersal and gene flow spanning three oceans were demonstrated between the species in this complex. Analysis of migration rates estimated using bayesass revealed unidirectional movement of petrels from the Atlantic and Pacific into the Indian Ocean. Conversely, structure analysis revealed gene flow between species of the Atlantic and Pacific oceans, with potential three-way hybrids occurring outside the Indian Ocean. Additionally, geolocation tracking of Round Island petrels revealed two individuals travelling to the Atlantic and Pacific. These results suggest that interspecific hybrids in Pterodroma petrels are more common than was previously assumed. This study is the first of its kind to investigate gene flow between populations of closely related Procellariiform species on a global scale, demonstrating the need for consideration of widespread migration and hybridization in the conservation of threatened seabirds.

Genetic structure among greater white-fronted goose populations of the Pacific Flyway

An understanding of the genetic structure of populations in the wild is essential for long-term conservation and stewardship in the face of environmental change. Knowledge of the present-day distribution of genetic lineages (phylogeography) of a species is especially important for organisms that are exploited or utilize habitats that may be jeopardized by human intervention, including climate change. Here, we describe mitochondrial (mtDNA) and nuclear genetic (microsatellite) diversity among three populations of a migratory bird, the greater white-fronted goose (Anser albifrons), which breeds discontinuously in western and southwestern Alaska and winters in the Pacific Flyway of North America. Significant genetic structure was evident at both marker types. All three populations were differentiated for mtDNA, whereas microsatellite analysis only differentiated geese from the Cook Inlet Basin. In sexual reproducing species, nonrandom mate selection, when occurring in concert with fine-scale resource partitioning, can lead to phenotypic and genetic divergence as we observed in our study. If mate selection does not occur at the time of reproduction, which is not uncommon in long-lived organisms, then mechanisms influencing the true availability of potential mates may be obscured, and the degree of genetic and phenotypic diversity may appear incongruous with presumed patterns of gene flow. Previous investigations revealed population-specific behavioral, temporal, and spatial mechanisms that likely influence the amount of gene flow measured among greater white-fronted goose populations. The degree of observed genetic structuring aligns well with our current understanding of population differences pertaining to seasonal movements, social structure, pairing behavior, and resource partitioning.

Population genomics of an endemic Mediterranean fish: differentiation by fine scale dispersal and adaptation

The assessment of the genetic structuring of biodiversity is crucial for management and conservation. For species with large effective population sizes a low number of markers may fail to identify population structure. A solution of this shortcoming can be high-throughput sequencing that allows genotyping thousands of markers on a genome-wide approach while facilitating the detection of genetic structuring shaped by selection. We used Genotyping-by-Sequencing (GBS) on 176 individuals of the endemic East Atlantic peacock wrasse (Symphodus tinca), from 6 locations in the Adriatic and Ionian seas. We obtained a total of 4,155 polymorphic SNPs and we observed two strong barriers to gene flow. The first one differentiated Tremiti Islands, in the northwest, from all the other locations while the second one separated east and south-west localities. Outlier SNPs potentially under positive selection and neutral SNPs both showed similar patterns of structuring, although finer scale differentiation was unveiled with outlier loci. Our results reflect the complexity of population genetic structure and demonstrate that both habitat fragmentation and positive selection are on play. This complexity should be considered in biodiversity assessments of different taxa, including non-model yet ecologically relevant organisms.

Genetic diversity of Atlantic Bluefin tuna in the Mediterranean Sea: insights from genome-wide SNPs and microsatellites

Background

Elucidating the patterns of the Atlantic Bluefin tuna [ABFT, Thunnus thynnus (Linnaeus, 1758)] population structure constitutes a challenging task of great importance. Most of the unique challenges stem from its biology, as well as the attributes of the marine realm in which it disperses. Accurate information is urgently needed for stock assessment, and the identification of critical features to the persistence and adaptation of populations in order to formulate and adopt effective strategies for ABFT conservation and management. Conclusions of a great number of ABFT genetic studies on the Mediterranean Sea stock structure are rather controversial and not yet conclusive. In this study, ABFT genomic diversity was investigated in the Mediterranean Sea, which is the most important area for the species’ reproduction.

Results

Analyzing genome-wide SNPs and microsatellites from ABFT samples collected throughout the Mediterranean Sea did not provide strong evidence of genetic structure, pointing towards the existence of a single panmictic unit. An alternative view would recognize a failure to reject the null hypothesis of a panmictic unit as an effect of the study’s sampling design, the type of markers used, and the effectiveness/suitability of analysis methods in respect to the species biological characteristics or any combination of the above.

Conclusions

Unravelling the drivers of ABFT population diversity would require the consideration of important aspects of the species spawning behavior for the determination of the appropriate sampling design. Novel approaches and methods of analysis that will bring together experts in genetics/-omics, ecology and oceanography are deemed necessary. Analyzing ABFT genetic data under the discipline of seascape genetics could provide the analysis framework under which major abiotic and biotic forces controlling ABFT recruitment could be identified, elucidating the complicated population dynamics of the species, while multiple and continuous fisheries monitoring should in all cases be considered as a prerequisite in order to achieve efficient and long-term ABFT conservation.

Level of genetic differentiation affects relative performances of expressed sequence tag and genomic SSRs

Microsatellites, also called simple sequence repeats (SSRs), are markers of choice to estimate relevant parameters for conservation genetics, such as migration rates, effective population size and kinship. Cross-amplification of SSRs is the simplest way to obtain sets of markers, and highly conserved SSRs have recently been developed from expressed sequence tags (EST) to improve SSR cross-species utility. As EST-SSRs are located in coding regions, the higher stability of their flanking regions reduces the frequency of null alleles and improves cross-species amplification. However, EST-SSRs have generally less allelic variability than genomic SSRs, potentially leading to differences in estimates of population genetic parameters such as genetic differentiation. To assess the potential of EST-SSRs in studies of within-species genetic diversity, we compared the relative performance of EST- and genomic SSRs following a multispecies approach on passerine birds. We tested whether patterns and levels of genetic diversity within and between populations assessed from EST- and from genomic SSRs are congruent, and we investigated how the relative efficiency of EST- and genomic SSRs is influenced by levels of differentiation. EST- and genomic SSRs ensured comparable inferences of population genetic structure in cases of strong genetic differentiation, and genomic SSRs performed slightly better than EST-SSRs when differentiation is moderate. However and interestingly, EST-SSRs had a higher power to detect weak genetic structure compared to genomic SSRs. Our study attests that EST-SSRs may be valuable molecular markers for conservation genetic studies in taxa such as birds, where the development of genomic SSRs is impeded by their low frequency.

Genotyping-by-sequencing of pear (Pyrus spp.) accessions unravels novel patterns of genetic diversity and selection footprints

Understanding of genetic diversity and marker-trait relationships in pears (Pyrus spp.) forms an important part of gene conservation and cultivar breeding. Accessions of Asian and European pear species, and interspecific hybrids were planted in a common garden experiment. Genotyping-by-sequencing (GBS) was used to genotype 214 accessions, which were also phenotyped for fruit quality traits. A combination of selection scans and association analyses were used to identify signatures of selection. Patterns of genetic diversity, population structure and introgression were also investigated. About 15 000 high-quality SNP markers were identified from the GBS data, of which 25% and 11% harboured private alleles for European and Asian species, respectively. Bayesian clustering analysis suggested negligible gene flow, resulting in highly significant population differentiation (F_st=0.45) between Asian and European pears. Interspecific hybrids displayed an average of 55% and 45% introgression from their Asian and European ancestors, respectively. Phenotypic (firmness, acidity, shape and so on) variation between accessions was significantly associated with genetic differentiation. Allele frequencies at large-effect SNP loci were significantly different between genetic groups, suggesting footprints of directional selection. Selection scan analyses identified over 20 outlier SNP loci with substantial statistical support, likely to be subject to directional selection or closely linked to loci under selection.

Characterization, development and multiplexing of microsatellite markers in three commercially exploited reef fish and their application for stock identification

Thirty-four microsatellite loci were isolated from three reef fish species; golden snapper Lutjanus johnii, blackspotted croaker Protonibea diacanthus and grass emperor Lethrinus laticaudis using a next generation sequencing approach. Both IonTorrent single reads and Illumina MiSeq paired-end reads were used, with the latter demonstrating a higher quality of reads than the IonTorrent. From the 1–1.5 million raw reads per species, we successfully obtained 10–13 polymorphic loci for each species, which satisfied stringent design criteria. We developed multiplex panels for the amplification of the golden snapper and the blackspotted croaker loci, as well as post-amplification pooling panels for the grass emperor loci. The microsatellites characterized in this work were tested across three locations of northern Australia. The microsatellites we developed can detect population differentiation across northern Australia and may be used for genetic structure studies and stock identification.

Secondary contact and asymmetrical gene flow in a cosmopolitan marine fish across the Benguela upwelling zone

The combination of oceanographic barriers and habitat heterogeneity are known to reduce connectivity and leave specific genetic signatures in the demographic history of marine species. However, barriers to gene flow in the marine environment are almost never impermeable which inevitably allows secondary contact to occur. In this study, eight sampling sites (five along the South African coastline, one each in Angola, Senegal and Portugal) were chosen to examine the population genetic structure and phylogeographic history of the cosmopolitan bluefish (Pomatomus saltatrix), distributed across a large South-east Atlantic upwelling zone. Molecular analyses were applied to mtDNA cytochrome b, intron AM2B1 and 15 microsatellite loci. We detected uncharacteristically high genetic differentiation (F_ST 0.15-0.20; P<0.001) between the fish sampled from South Africa and the other sites, strongly influenced by five outlier microsatellite loci located in conserved intergenic regions. In addition, differentiation among the remaining East Atlantic sites was detected, although mtDNA indicated past isolation with subsequent secondary contact between these East Atlantic populations. We further identified secondary contact, with unidirectional gene flow from South Africa to Angola. The directional contact is likely explained by a combination of the northward flowing offshore current and endogenous incompatibilities restricting integration of certain regions of the genome and limiting gene flow to the south. The results confirm that the dynamic system associated with the Benguela current upwelling zone influences species distributions and population processes in the South-east Atlantic.

Understanding the Spatial Scale of Genetic Connectivity at Sea: Unique Insights from a Land Fish and a Meta-Analysis

Quantifying the spatial scale of population connectivity is important for understanding the evolutionary potential of ecologically divergent populations and for designing conservation strategies to preserve those populations. For marine organisms like fish, the spatial scale of connectivity is generally set by a pelagic larval phase. This has complicated past estimates of connectivity because detailed information on larval movements are difficult to obtain. Genetic approaches provide a tractable alternative and have the added benefit of estimating directly the reproductive isolation of populations. In this study, we leveraged empirical estimates of genetic differentiation among populations with simulations and a meta-analysis to provide a general estimate of the spatial scale of genetic connectivity in marine environments. We used neutral genetic markers to first quantify the genetic differentiation of ecologically-isolated adult populations of a land dwelling fish, the Pacific leaping blenny (Alticus arnoldorum), where marine larval dispersal is the only probable means of connectivity among populations. We then compared these estimates to simulations of a range of marine dispersal scenarios and to collated FST and distance data from the literature for marine fish across diverse spatial scales. We found genetic connectivity at sea was extensive among marine populations and in the case of A. arnoldorum, apparently little affected by the presence of ecological barriers. We estimated that ~5000 km (with broad confidence intervals ranging from 810-11,692 km) was the spatial scale at which evolutionarily meaningful barriers to gene flow start to occur at sea, although substantially shorter distances are also possible for some taxa. In general, however, such a large estimate of connectivity has important implications for the evolutionary and conservation potential of many marine fish communities.

Phylogenetics and population genetics of Plotosus canius (Siluriformes: Plotosidae) from Malaysian coastal waters

Plotosus canius (Hamilton, 1822) is a significant marine species in Malaysia from nutritional and commercial perspectives. Despite numerous fundamental research on biological characteristics of P. canius, there are various concerns on the level of population differentiation, genomic structure, and the level of genetic variability among their populations due to deficiency of genetic-based studies. Deficiency on basic contexts such as stock identification, phylogenetic relationship and population genetic structure would negatively impact their sustainable conservation. Hence, this study was conducted to characterize the genetic structure of P. canius for the first time through the application of mitochondrial Cytochrome Oxidase I (COI) gene, cross amplification of Tandanus tandanus microsatellites, and a total of 117 collected specimens across five selected populations of Malaysia. The experimental results of the mitochondrial analysis revealed that the haplotype diversity and nucleotide diversity varied from 0.395–0.771 and 0.033–0.65 respectively. Moreover, the statistical analysis of microsatellites addressed a considerable heterozygote insufficiency in all populations, with average observed heterozygosity (H_o) value of 0.2168, which was lower than the standard heterozygosity in marine populations (H_o = 0.79). This alongside the high Fis values estimation, high pairwise differentiation among populations and low within population variations are supposed to be associated with small sample size, and inbreeding system. Besides, the significant finding of this study was the sharing of common haplotype KR086940, which reflects a historical genetic connectivity between Peninsular Malaysia and Borneo populations due to the geological history of Southeast Asia during Pleistocene era. Demographic analyses showed that all populations were in an equilibrium state with no significant evidence of population expansion. To put it briefly, the current study has managed to provide an initial genomic database toward understanding of the genetic characterization, phylogenetic, molecular diversification and population structure in P. canius, and should be necessary highlighted for appropriate management and conservation of species. Further studies must be carried out involving more geographical and sampling sites, larger population size per site, and utilization of species specific microsatellites loci.

Contemporary genetic structure and postglacial demographic history of the black scorpionfish, Scorpaena porcus, in the Mediterranean and the Black Seas

Understanding the distribution of genetic diversity in the light of past demographic events linked with climatic shifts will help to forecast evolutionary trajectories of ecosystems within the current context of climate change. In this study, mitochondrial sequences and microsatellite loci were analysed using traditional population genetic approaches together with Bayesian dating and the more recent approximate Bayesian computation scenario testing. The genetic structure and demographic history of a commercial fish, the black scorpionfish, Scorpaena porcus, was investigated throughout the Mediterranean and Black Seas. The results suggest that the species recently underwent population expansions, in both seas, likely concomitant with the warming period following the Last Glacial Maximum, 20 000 years ago. A weak contemporaneous genetic differentiation was identified between the Black Sea and the Mediterranean Sea. However, the genetic diversity was similar for populations of the two seas, suggesting a high number of colonizers entered the Black Sea during the interglacial period and/or the presence of a refugial population in the Black Sea during the glacial period. Finally, within seas, an east/west genetic differentiation in the Adriatic seems to prevail, whereas the Black Sea does not show any structured spatial genetic pattern of its population. Overall, these results suggest that the Black Sea is not that isolated from the Mediterranean, and both seas revealed similar evolutionary patterns related to climate change and changes in sea level.

Genetic structure of the crown-of-thorns seastar in the Pacific Ocean, with focus on Guam

Population outbreaks of the corallivorous crown-of-thorns seastar (COTS), Acanthaster 'planci' L., are among the most important biological disturbances of tropical coral reefs. Over the past 50 years, several devastating outbreaks have been documented around Guam, an island in the western Pacific Ocean. Previous analyses have shown that in the Pacific Ocean, COTS larval dispersal may be geographically restricted to certain regions. Here, we assess the genetic structure of Pacific COTS populations and compared samples from around Guam with a number of distant localities in the Pacific Ocean, and focused on determining the degree of genetic structure among populations previously considered to be isolated. Using microsatellites, we document substantial genetic structure between 14 localities from different geographical regions in the Pacific Ocean. Populations from the 14 locations sampled were found to be structured in three significantly differentiated groups: (1) all locations immediately around Guam, as well as Kingman Reef and Swains Island; (2) Japan, Philippines, GBR and Vanuatu; and (3) Johnston Atoll, which was significantly different from all other localities. The lack of genetic differentiation between Guam and extremely distant populations from Kingman Reef and Swains Island suggests potential long-distance dispersal of COTS in the Pacific.

The Structure of Genetic Diversity in Eelgrass (Zostera marina L.) along the North Pacific and Bering Sea Coasts of Alaska

Eelgrass (Zostera marina) populations occupying coastal waters of Alaska are separated by a peninsula and island archipelago into two Large Marine Ecosystems (LMEs). From populations in both LMEs, we characterize genetic diversity, population structure, and polarity in gene flow using nuclear microsatellite fragment and chloroplast and nuclear sequence data. An inverse relationship between genetic diversity and latitude was observed (heterozygosity: R2 = 0.738, P < 0.001; allelic richness: R2 = 0.327, P = 0.047), as was significant genetic partitioning across most sampling sites (θ = 0.302, P < 0.0001). Variance in allele frequency was significantly partitioned by region only in cases when a population geographically in the Gulf of Alaska LME (Kinzarof Lagoon) was instead included with populations in the Eastern Bering Sea LME (θp = 0.128-0.172; P < 0.003), suggesting gene flow between the two LMEs in this region. Gene flow among locales was rarely symmetrical, with notable exceptions generally following net coastal ocean current direction. Genetic data failed to support recent proposals that multiple Zostera species (i.e. Z. japonica and Z. angustifolia) are codistributed with Z. marina in Alaska. Comparative analyses also failed to support the hypothesis that eelgrass populations in the North Atlantic derived from eelgrass retained in northeastern Pacific Last Glacial Maximum refugia. These data suggest northeastern Pacific populations are derived from populations expanding northward from temperate populations following climate amelioration at the terminus of the last Pleistocene glaciation.

Multilocus evidence for globally distributed cryptic species and distinct populations across ocean gyres in a mesopelagic copepod

Zooplanktonic taxa have a greater number of distinct populations and species than might be predicted based on their large population sizes and open-ocean habitat, which lacks obvious physical barriers to dispersal and gene flow. To gain insight into the evolutionary mechanisms driving genetic diversification in zooplankton, we developed eight microsatellite markers to examine the population structure of an abundant, globally distributed mesopelagic copepod, Haloptilus longicornis, at 18 sample sites across the Atlantic and Pacific Oceans (n = 761). When comparing our microsatellite results with those of a prior study that used a mtDNA marker (mtCOII, n = 1059, 43 sample sites), we unexpectedly found evidence for the presence of a cryptic species pair. These species were globally distributed and apparently sympatric, and were separated by relatively weak genetic divergence (reciprocally monophyletic mtCOII lineages 1.6% divergent; microsatellite FST ranging from 0.28 to 0.88 across loci, P < 0.00001). Using both mtDNA and microsatellite data for the most common of the two species (n = 669 for microsatellites, n = 572 for mtDNA), we also found evidence for allopatric barriers to gene flow within species, with distinct populations separated by continental landmasses and equatorial waters in both the Atlantic and Pacific Ocean basins. Our study shows that oceanic barriers to gene flow can act as a mechanism promoting allopatric diversification in holoplanktonic taxa, despite the high potential dispersal abilities and pelagic habitat for these species.

Multilocus Bayesian Estimates of Intra-Oceanic Genetic Differentiation, Connectivity, and Admixture in Atlantic Swordfish (Xiphias gladius L.)

Previous genetic studies of Atlantic swordfish (Xiphias gladius L.) revealed significant differentiation among Mediterranean, North Atlantic and South Atlantic populations using both mitochondrial and nuclear DNA data. However, limitations in geographic sampling coverage, and the use of single loci, precluded an accurate placement of boundaries and of estimates of admixture. In this study, we present multilocus analyses of 26 single nucleotide polymorphisms (SNPs) within 10 nuclear genes to estimate population differentiation and admixture based on the characterization of 774 individuals representing North Atlantic, South Atlantic, and Mediterranean swordfish populations. Pairwise F_ST values, AMOVA, PCoA, and Bayesian individual assignments support the differentiation of swordfish inhabiting these three basins, but not the current placement of the boundaries that separate them. Specifically, the range of the South Atlantic population extends beyond 5°N management boundary to 20°N-25°N from 45°W. Likewise the Mediterranean population extends beyond the current management boundary at the Strait of Gibraltar to approximately 10°W. Further, admixture zones, characterized by asymmetric contributions of adjacent populations within samples, are confined to the Northeast Atlantic. While South Atlantic and Mediterranean migrants were identified within these Northeast Atlantic admixture zones no North Atlantic migrants were identified respectively in these two neighboring basins. Owing to both, the characterization of larger number of loci and a more ample spatial sampling coverage, it was possible to provide a finer resolution of the boundaries separating Atlantic swordfish populations than previous studies. Finally, the patterns of population structure and admixture are discussed in the light of the reproductive biology, the known patterns of dispersal, and oceanographic features that may act as barriers to gene flow to Atlantic swordfish.

High connectivity in Rastrelliger kanagurta: influence of historical signatures and migratory behaviour inferred from mtDNA cytochrome b

Phylogeographic patterns and population structure of the pelagic Indian mackerel, Rastrelliger kanagurta were examined in 23 populations collected from the Indonesian-Malaysian Archipelago (IMA) and the West Indian Ocean (WIO). Despite the vast expanse of the IMA and neighbouring seas, no evidence for geographical structure was evident. An indication that R. kanagurta populations across this region are essentially panmictic. This study also revealed that historical isolation was insufficient for R. kanagurta to attain migration drift equilibrium. Two distinct subpopulations were detected between the WIO and the IMA (and adjacent populations); interpopulation genetic variation was high. A plausible explanation for the genetic differentiation observed between the IMA and WIO regions suggest historical isolation as a result of fluctuations in sea levels during the late Pleistocene. This occurrence resulted in the evolution of a phylogeographic break for this species to the north of the Andaman Sea.

Genotyping-by-sequencing for Populus population genomics: an assessment of genome sampling patterns and filtering approaches

Continuing advances in nucleotide sequencing technology are inspiring a suite of genomic approaches in studies of natural populations. Researchers are faced with data management and analytical scales that are increasing by orders of magnitude. With such dramatic advances comes a need to understand biases and error rates, which can be propagated and magnified in large-scale data acquisition and processing. Here we assess genomic sampling biases and the effects of various population-level data filtering strategies in a genotyping-by-sequencing (GBS) protocol. We focus on data from two species of Populus, because this genus has a relatively small genome and is emerging as a target for population genomic studies. We estimate the proportions and patterns of genomic sampling by examining the Populus trichocarpa genome (Nisqually-1), and demonstrate a pronounced bias towards coding regions when using the methylation-sensitive ApeKI restriction enzyme in this species. Using population-level data from a closely related species (P. tremuloides), we also investigate various approaches for filtering GBS data to retain high-depth, informative SNPs that can be used for population genetic analyses. We find a data filter that includes the designation of ambiguous alleles resulted in metrics of population structure and Hardy-Weinberg equilibrium that were most consistent with previous studies of the same populations based on other genetic markers. Analyses of the filtered data (27,910 SNPs) also resulted in patterns of heterozygosity and population structure similar to a previous study using microsatellites. Our application demonstrates that technically and analytically simple approaches can readily be developed for population genomics of natural populations.

Testing the role of meander cutoff in promoting gene flow across a riverine barrier in ground skinks (Scincella lateralis)

Despite considerable attention, the long-term impact of rivers on species diversification remains uncertain. Meander loop cutoff (MLC) is one river phenomenon that may compromise a river's diversifying effects by passively transferring organisms from one side of the river to the other. However, the ability of MLC to promote gene flow across rivers has not been demonstrated empirically. Here, we test several predictions of MLC-mediated gene flow in populations of North American ground skinks (Scincella lateralis) separated by a well-established riverine barrier, the Mississippi River: 1) individuals collected from within meander cutoffs should be more closely related to individuals across the river than on the same side, 2) individuals within meander cutoffs should contain more immigrants than individuals away from meander cutoffs, 3) immigration rates estimated across the river should be highest in the direction of the cutoff event, and 4) the distribution of alleles native to one side of the river should be better predicted by the historical rather than current path of the river. To test these predictions we sampled 13 microsatellite loci and mitochondrial DNA from ground skinks collected near three ancient meander loops. These predictions were generally supported by genetic data, although support was stronger for mtDNA than for microsatellite data. Partial support for genetic divergence of samples within ancient meander loops also provides evidence for the MLC hypothesis. Although a role for MLC-mediated gene flow was supported here for ground skinks, the transient nature of river channels and morphologies may limit the long-term importance of MLC in stemming population divergence across major rivers.

Gene flow at major transitional areas in sea bass (Dicentrarchus labrax) and the possible emergence of a hybrid swarm

The population genetic structure of sea bass (Dicentrarchus labrax) along a transect from the Atlantic Ocean (AO) to the Eastern Mediterranean (EM) Sea differs from that of most other marine taxa in this area. Three populations (AO, Western Mediterranean [WM], EM) are recognized today, which were originally two allopatric populations. How two ancestral genetic units have evolved into three distinct units has not been addressed yet. Therefore, to investigate mechanisms that lead to the emergence of the central WM population, its current status, and its connectivity with the two parental populations, we applied 20 nuclear loci that were either gene associated or gene independent. Results confirmed the existence of three distinct gene pools, with higher differentiation at two transitional areas, the Almeria-Oran Front (AOF) and of the Siculo-Tunisian Strait (STS), than within any population. Significant linkage disequilibrium and heterozygote excess indicated that the STS is probably another tension zone, as already described for the AOF. Neutrality tests fail to reveal marker loci that could be driven by selection within or among metapopulations, except for locus DLA0068. Collectively, results support that the central WM population arose by trapping two tensions zones at distinct geographic locations of limited connectivity. Population assignment further revealed that WM individuals were more introgressed than individuals from the other two metapopulations. This suggests that this population might result from hybrid swarming, and was or is still seeded by genes received through the filter of each tension zone.

Genetic stock structure of walleye pollock (Theragra chalcogramma) inferred by PCR-RFLP analysis of the mitochondrial DNA and SNP analysis of nuclear DNA

Walleye pollock, Theragra chalocogramma, is one of the most important species in the North Pacific and Bering Sea ecosystems. However genetic population structuring of walleye pollock is uncertain. In the present study, genetic variation of walleye pollock collected in several spawning areas ranging from the Japan Sea to the Gulf of Alaska was investigated by DNA analysis. Three regions of the spacer control region, the ND5 and ND6 region (ND complex), and the ND1 and 16S rRNA region (rDNA complex) were amplified using the polymerase chain reaction (PCR). Restriction fragment length polymorphism (RFLP) was conducted on these PCR products and composite haplotypes were calculated. Furthermore, several nuclear DNA regions (actin, Calmodulin, S7 ribosomal protein, creatin kinase, and SypI gene) were investigated to study the stock structure of walleye pollock. It was considered that Calmodulin gene was one of good genetic marker, therefore we conducted SNP analysis for Calmodulin gene by SnaPshot kits. In RFLP analyses, there were no area-specific fragment patterns in the three regions, control region, ND complex and rDNA complex of mtDNA. However compositions of the fragment patterns for the three digested sets, control region/HinfI, rDNA complex/MspI and ND complex/MspI indicated that there are significant differences between around the Japan (Sado-Funka Bay-Wakkanai-Rausu) and the Bering Sea (Western Bering Sea-Nabarin-Atka I.-Bogoslof I). Furthermore, in the case of haplotype frequency, composition showed also significant genetic difference between two areas. Moreover, in Calmodulin analyses, haplotype compositions were changing from western area to eastern area gradually and the results of AMOVA analysis showed that there are interesting differences between western Pacific, western Bering Sea, and eastern Bering Sea. Judging from these results, it was considered that there are three populations of walleye pollock in the Northern Ocean. However, area-specific pattern was not found in some populations in the Northern Ocean. Therefore, we suggested that these populations were related by weak gene flow, and the walleye pollock was formed with meta-population around the Japan and the Bering Sea.

Phylogeography of two closely related Indo-Pacific butterflyfishes reveals divergent evolutionary histories and discordant results from mtDNA and microsatellites

Marine biogeographic barriers can have unpredictable consequences, even among closely related species. To resolve phylogeographic patterns for Indo-Pacific reef fauna, we conducted range-wide surveys of sister species, the scrawled butterflyfish (Chaetodon meyeri; N = 134) and the ornate butterflyfish (Chaetodon ornatissimus; N = 296), using mitochondrial DNA cytochrome b sequences and 10 microsatellite loci. The former is distributed primarily in the Indian Ocean but also extends to the Line Islands in the Central Pacific, whereas the latter is distributed primarily in the Central-West Pacific (including Hawaii and French Polynesia) but extends to the eastern margin of the Indian Ocean. Analyses of molecular variance and Bayesian STRUCTURE results revealed 1 range-wide group for C. meyeri and 3 groups for C. ornatissimus: 1) eastern Indian Ocean and western Pacific, 2) Central Pacific, and 3) Hawaii. Estimates of the last population expansion were much more recent for C. meyeri (61 500 to 95 000 years) versus C. ornatissimus (184 700 to 286 300 years). Despite similarities in ecology, morphology, life history, and a broadly overlapping distribution, these sister species have divergent patterns of dispersal and corresponding evolutionary history. The mtDNA and microsatellite markers did not provide concordant results within 1 of our study species (C. meyeri), or in 7 out of 12 other cases of marine fishes in the published literature. This discordance renews caution in relying on one or a few markers for reconstructing historical demography.

Common misconceptions in molecular ecology: echoes of the modern synthesis

The field of molecular ecology has burgeoned into a large discipline spurred on by technical innovations that facilitate the rapid acquisition of large amounts of genotypic data, by the continuing development of theory to interpret results, and by the availability of computer programs to analyse data sets. As the discipline grows, however, misconceptions have become enshrined in the literature and are perpetuated by routine citations to other articles in molecular ecology. These misconceptions hamper a better understanding of the processes that influence genetic variation in natural populations and sometimes lead to erroneous conclusions. Here, we consider eight misconceptions commonly appearing in the literature: (i) some molecular markers are inherently better than other markers; (ii) mtDNA produces higher F(ST) values than nDNA; (iii) estimated population coalescences are real; (iv) more data are always better; (v) one needs to do a Bayesian analysis; (vi) selective sweeps influence mtDNA data; (vii) equilibrium conditions are critical for estimating population parameters; and (viii) having better technology makes us smarter than our predecessors. This is clearly not an exhaustive list and many others can be added. It is, however, sufficient to illustrate why we all need to be more critical of our own understanding of molecular ecology and to be suspicious of self-evident truths.

Neither philopatric nor panmictic: microsatellite and mtDNA evidence suggests lack of natal homing but limits to dispersal in Pacific lamprey

Most species with lengthy migrations display some degree of natal homing; some (e.g. migratory birds and anadromous salmonids) show spectacular feats of homing. However, studies of the sea lamprey (Petromyzon marinus) indicate that this anadromous species locates spawning habitat based on pheromonal cues from larvae rather than through philopatry. Previous genetic studies in the anadromous Pacific lamprey (Entosphenus tridentatus) have both supported and rejected the hypothesis of natal homing. To resolve this, we used nine microsatellite loci to examine the population structure in 965 Pacific lamprey from 20 locations from central British Columbia to southern California and supplemented this analysis with mitochondrial DNA restriction fragment length polymorphism analysis on a subset of 530 lamprey. Microsatellite analysis revealed (i) relatively low but often statistically significant genetic differentiation among locations (97% pairwise F(ST) values were <0.04 but 73.7% were significant); and (ii) weak but significant isolation by distance (r(2) = 0.0565, P = 0.0450) but no geographic clustering of samples. The few moderate F(ST) values involved comparisons with sites that were geographically distant or far upstream. The mtDNA analysis--although providing less resolution among sites (only 4.7%F(ST) values were significant)--was broadly consistent with the microsatellite results: (i) the southernmost site and some sites tributary to the Salish Sea were genetically distinct; and (ii) southern sites showed higher haplotype and private haplotype richness. These results are inconsistent with philopatry, suggesting that anadromous lampreys are unusual among species with long migrations, but suggest that limited dispersal at sea precludes panmixia in this species.

Genetic structure and the North American postglacial expansion of the barnacle, Semibalanus balanoides

Population genetic characteristics are shaped by the life-history traits of organisms and the geologic history of their habitat. This study provides a neutral framework for understanding the population dynamics and opportunities for selection in Semibalanus balanoides, a species that figures prominently in ecological and evolutionary studies in the Atlantic intertidal. We used mitochondrial DNA (mtDNA) control region (N = 131) and microsatellite markers (∼40 individuals/site/locus) to survey populations of the broadly dispersing acorn barnacle from 8 sites spanning 800 km of North American coast and 1 site in Europe. Patterns of mtDNA sequence evolution were consistent with larger population sizes in Europe and population expansion at the conclusion of the last ice age, approximately 20 000 years ago, in North America. A significant portion of mitochondrial diversity was partitioned between the continents (ϕ(ST) = 0.281), but there was only weak structure observed from mtDNA within North America. Microsatellites showed significant structuring between the continents (F(ST) = 0.021) as well as within North America (F(ST) = 0.013). Isolation by distance in North America was largely driven by a split between populations south of Cape Cod and all others (P < 10(-4)). The glacial events responsible for generating allelic diversity at mtDNA and microsatellites may also be responsible for generating selectable variation at metabolic enzymes in S. balanoides.

Species delineation using Bayesian model-based assignment tests: a case study using Chinese toad-headed agamas (genus Phrynocephalus)

Background

Species are fundamental units in biology, yet much debate exists surrounding how we should delineate species in nature. Species discovery now requires the use of separate, corroborating datasets to quantify independently evolving lineages and test species criteria. However, the complexity of the speciation process has ushered in a need to infuse studies with new tools capable of aiding in species delineation. We suggest that model-based assignment tests are one such tool. This method circumvents constraints with traditional population genetic analyses and provides a novel means of describing cryptic and complex diversity in natural systems. Using toad-headed agamas of the Phrynocephalus vlangalii complex as a case study, we apply model-based assignment tests to microsatellite DNA data to test whether P. putjatia, a controversial species that closely resembles P. vlangalii morphologically, represents a valid species. Mitochondrial DNA and geographic data are also included to corroborate the assignment test results.

Results

Assignment tests revealed two distinct nuclear DNA clusters with 95% (230/243) of the individuals being assigned to one of the clusters with > 90% probability. The nuclear genomes of the two clusters remained distinct in sympatry, particularly at three syntopic sites, suggesting the existence of reproductive isolation between the identified clusters. In addition, a mitochondrial ND2 gene tree revealed two deeply diverged clades, which were largely congruent with the two nuclear DNA clusters, with a few exceptions. Historical mitochondrial introgression events between the two groups might explain the disagreement between the mitochondrial and nuclear DNA data. The nuclear DNA clusters and mitochondrial clades corresponded nicely to the hypothesized distributions of P. vlangalii and P. putjatia.

Conclusions

These results demonstrate that assignment tests based on microsatellite DNA data can be powerful tools for distinguishing closely related species and support the validity of P. putjatia. Assignment tests have the potential to play a significant role in elucidating biodiversity in the era of DNA data. Nonetheless, important limitations do exist and multiple independent datasets should be used to corroborate results from assignment programs.

Influence of mutation rate on estimators of genetic differentiation--lessons from Arabidopsis thaliana

BACKGROUND

The analysis of molecular variation within and between populations is crucial to establish strategies for conservation as well as to detect the footprint of spatially heterogeneous selection. The traditional estimator of genetic differentiation (F(ST)) has been shown to be misleading if genetic diversity is high. Alternative estimators of F(ST) have been proposed, but their robustness to variation in mutation rate is not clearly established. We first investigated the effect of mutation and migration rate using computer simulations and examined their joint influence on Q(ST), a measure of genetic differentiation for quantitative traits. We further used experimental data in natural populations of Arabidopsis thaliana to characterize the effect of mutation rate on various estimates of population differentiation. Since natural species exhibit various degrees of self-fertilisation, we also investigated the effect of mating system on the different estimators.

RESULTS

If mutation rate is high and migration rate low, classical measures of genetic differentiation are misleading. Only Phi(ST), an estimator that takes the mutational distances between alleles into account, is independent of mutation rate, for all migration rates. However, the performance of Phi(ST) depends on the underlying mutation model and departures from this model cause its performance to degrade. We further show that Q(ST) has the same bias. We provide evidence that, in A. thaliana, microsatellite variation correlates with mutation rate. We thereby demonstrate that our results on estimators of genetic differentiation have important implications, even for species that are well established models in population genetics and molecular biology.

CONCLUSIONS

We find that alternative measures of differentiation like F'(ST) and D are not suitable for estimating effective migration rate and should not be used in studies of local adaptation. Genetic differentiation should instead be measured using an estimator that takes mutation rate into account, such as Phi(ST). Furthermore, in systems where migration between populations is low, such as A. thaliana, Q(ST) < F(ST) cannot be taken as evidence for homogenising selection as has been traditionally thought.

Mito-nuclear discordance in the degree of population differentiation in a marine goby

An increasing number of phylogeographic studies on marine species shows discordant patterns in the degree of population differentiation between nuclear and mitochondrial markers. To understand better which factors have the potential to cause these patterns of discordance in marine organisms, a population genetic study was realized on the sand goby Pomatoschistus minutus (Pallas 1770; Gobiidae, Teleostei). Sand gobies from eight European locations were genotyped at eight microsatellite markers. Microsatellites confirmed the global phylogeographical pattern of P. minutus observed with mitochondrial DNA (mtDNA) markers and nuclear allozyme markers. Three groups consistent with the mitochondrial lineages were defined (the Mediterranean, Iberian and North Atlantic groups) and indications of a recent founder event in the northern Baltic Sea were found. Nevertheless, differences in the degree of population differentiation between the nuclear and mitochondrial markers were large (global F(ST)-values for microsatellites=0.0121; for allozymes=0.00831; for mtDNA=0.4293). Selection, sex-biased dispersal, homoplasy and a high effective population size are generally accepted as explanations for this mitonuclear discrepancy in the degree of population differentiation. In this study, selection on mtDNA and microsatellites, male-biased dispersal and homoplasy on microsatellite markers are unlikely to be a main cause for this discrepancy. The most likely reason for the discordant pattern is a recent demographical expansion of the sand goby, resulting in high effective population sizes slowing down the differentiation of nuclear DNA.

Structural and functional connectivity of marine fishes within a semi-enclosed Newfoundland fjord

The interplay between structural connectivity (i.e. habitat continuity) and functional connectivity (i.e. dispersal probability) in marine fishes was examined in a coastal fjord (Holyrood Pond, Newfoundland, Canada) that is completely isolated from the North Atlantic Ocean for most of the year. Genetic differentiation was described in three species (rainbow smelt Osmerus mordax, white hake Urophycis tenuis and Atlantic cod Gadus morhua) with contrasting life histories using seven to 10 microsatellite loci and a protein-coding locus, PanI (G. morhua). Analysis of microsatellite differentiation indicated clear genetic differences between the fjord and coastal regions; however, the magnitude of difference was no more elevated than adjacent bays and was not enhanced by the fjord's isolation. Osmerus mordax was characterized by the highest structure overall with moderate differentiation between the fjord and St Mary's Bay (F(ST)c.0.047). In contrast, U. tenuis and G. morhua displayed weak differentiation (F(ST) < 0.01). Nonetheless, these populations did demonstrate high rates (< 75%) of Bayesian self-assignment. Furthermore, elevated differentiation was observed at the PanI locus in G. morhua between the fjord and other coastal locations. Interestingly, locus-specific genetic differentiation and expected heterozygosity were negatively associated in O. mordax, in contrast to the positive associations observed in U. tenuis and G. morhua. Gene flow in these species is apparently unencumbered by limited structural connectivity, yet the observed differentiation suggests that population structuring exists over small scales despite high dispersal potential.