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.     

Population boundaries and genetic diversity in the endangered Mariana crow (Corvus kubaryi)

The Mariana crow (Corvus kubaryi) is an endangered species that is restricted to the islands of Guam and Rota in the Mariana archipelago. Predation by the introduced brown tree snake (Boiga irregularis) has decimated bird populations on Guam, and the crow population there is the last wild remnant of the endemic forest avifauna. The population on Guam is critically endangered and, despite intensive management, the population has continued to decline. Additional management options include intermixing the Guam and Rota populations, but such options are best evaluated within a population genetics framework. We used three types of molecular markers to assay genetic variation in the Mariana crow: mitochondrial DNA (mtDNA) sequences, minisatellites and microsatellites. The two populations could be differentiated by mtDNA sequencing and they differed in allele frequencies at nuclear markers. Thus, the populations could be designated as evolutionarily significant units. However, the Guam population is genetically more diverse than the Rota population, and its survival probability if managed separately is very low. All markers did indicate that the two populations are closely related and separated by a shallow genealogical division. Intermixing the populations is justified by two rationales. First, the apparent population differences may result from recent human activities. Second, a greater amount of genetic information may be preserved by joint management. The translocation of birds from Rota to Guam has begun, but strategies that will ensure maintenance of the variation in the Guam population warrant further exploration.     

Ancient,but not recent,population declines have had a genetic impact on alpine yellow-bellied toad populations,suggesting potential for complete recovery

Reduction in population size and local extinctions have been reported for the yellow-bellied toad, Bombina variegata, but the genetic impact of this is not yet known. In this study, we genotyped 200 individuals, using mtDNA cytochrome b and 11 nuclear microsatellites. We investigated fine-scale population structure and tested for genetic signatures of historical and recent population decline, using several statistical approaches, including likelihood methods and approximate Bayesian computation. Five major genetically divergent groups were found, largely corresponding to geography but with a clear exception of high genetic isolation in a highly touristic area. The effective sizes in the last few generations, as estimated from the random association among markers, never exceeded a few dozen of individuals. Our most important result is that several analyses converge in suggesting that genetic variation was shaped in all groups by a 7- to 45-fold demographic decline, which occurred between a few hundred and a few 1000 years ago. Remarkably, only weak evidence supports recent genetic impact related to human activities. We believe that the alpine B. variegata populations should be monitored and protected to stop their recent decline and to prevent local extinctions, with highest priority given to genetically isolated populations. Nonetheless, current genetic variation pattern, being mostly shaped in earlier times, suggests that complete recovery can be achieved. In general, our study is an example of how the potential for recovery should be inferred even under the co-occurrence of population decline, low genetic variation, and genetic bottleneck signals.     

Age, phylogeography and population structure of the microendemic banded spring snail, Mexipyrgus churinceanus

Recent theoretical and empirical studies of phylogeography and population structure indicate that many processes influence intraspecific evolutionary history. The present study represents the first examination of various forces influencing the spatial and temporal patterns of sequence variation in the freshwater Mexican banded spring snail, Mexipyrgus churinceanus. This snail occurs in one of the most critically endangered centres of freshwater endemism, the desert ecosystem of Cuatro Ciénegas. From cytochrome b mtDNA sequence variation, there is strong evidence of long-term isolation of three regions, suggesting that these regions represent evolutionarily distinct lineages. Molecular clock estimates of clade age indicate a time to most recent common ancestor of approximately 2.5 million years ago (Ma). The three regions differ considerably in the historical and demographic forces affecting population structure. The western populations have extremely low mtDNA diversity consistent with a severe bottleneck dating to 50,000 years before present (bp). The nearby Rio Mesquites drainage is characterized by fragmentation events, restricted gene flow with isolation by distance, and higher levels of mtDNA polymorphism. These patterns are consistent with the long-term stability of this drainage along with habitat heterogeneity and brooding contributing to population isolation and restricted gene flow. Southeastern populations show evidence of range expansion and a strong influence of genetic drift. Migration rates between drainages indicate very little gene flow between drainages except for asymmetric migration from the Rio Mesquites into both western and southeastern drainages.     

Genetic variation and phylogeography of free-living mouse species (genus Mus) in the Balkans and the Middle East

This work presents a study of the distribution and pattern of variation throughout the ranges of three free-living mouse species of the genus Mus-M. macedonicus, M. spicilegus, and a M. cypriacus - based on sequencing of two segments of the mitochondrial DNA (mtDNA) control region. The study shows a similar level of variability in the three species and suggests their recent population expansion. The highest proportion of variation is found within populations indicating low genetic structuring. Phylogenetic analysis confirms the significant divergence of a mitochondrial lineage of M. macedonicus from Israel, recently described as a new subspecies, M. macedonicus spretoides. Conversely, no genetic hiatus is revealed between European and Asian populations of M. macedonicus macedonicus. Although phylogenetic relationships among M. spicilegus populations could not be unravelled precisely, the results suggest a recent westward expansion of the species. The mtDNA divergence between M. macedonicus and M. spicilegus is 7.3%, suggesting their split between c. 700,000 and 1 million years ago. These dates correspond with a coalescent estimate about 720,000 years ago. On the other hand, M. cypriacus appeared almost twice as divergent from the former species (4.5%) as from the latter (8.8%) suggesting a divergence of c. 430,000-610,000 years ago (coalescent approximately 490,000 years ago) and 830,000-1.2 million years ago (coalescent approximately 780,000 years ago), respectively. Approximate times of population expansion have also been estimated for all taxa and groups of populations. Existence of several glacial refuges and various colonization scenarios are discussed; since all estimated divergence times fall within interglacial periods it seems that climatic oscillations did not play a crucial role in the evolution of the three species.     

Phylogenetic and phylogeographic analysis of Iberian lynx populations

The Iberian lynx (Lynx pardinus), one of the world's most endangered cat species, is vulnerable due to habitat loss, increased fragmentation of populations, and precipitous demographic reductions. An understanding of Iberian lynx evolutionary history is necessary to develop rational management plans for the species. Our objectives were to assess Iberian lynx genetic diversity at three evolutionary timescales. First we analyzed mitochondrial DNA (mtDNA) sequence variation to position the Iberian lynx relative to other species of the genus LYNX: We then assessed the pattern of mtDNA variation of isolated populations across the Iberian Peninsula. Finally we estimated levels of gene flow between two of the most important remaining lynx populations (Do?ana National Park and the Sierra Morena Mountains) and characterized the extent of microsatellite locus variation in these populations. Phylogenetic analyses of 1613 bp of mtDNA sequence variation supports the hypothesis that the Iberian lynx, Eurasian lynx, and Canadian lynx diverged within a short time period around 1.53-1.68 million years ago, and that the Iberian lynx and Eurasian lynx are sister taxa. Relative to most other felid species, genetic variation in mtDNA genes and nuclear microsatellites were reduced in Iberian lynx, suggesting that they experienced a fairly severe demographic bottleneck. In addition, the effects of more recent reductions in gene flow and population size are being manifested in local patterns of molecular genetic variation. These data, combined with recent studies modeling the viability of Iberian lynx populations, should provide greater urgency for the development and implementation of rational in situ and ex situ conservation plans.     

Phylogeography and population history of the endangered golden sun moth (Synemon plana) revealed by allozymes and mitochondrial DNA analysis

A combination of allozyme and mitochondrial DNA markers were used to determine the contribution of recent and ancient causes of patterns of genetic variation within and among 46 populations of the endangered golden sun moth, Synemon plana. Allozyme analysis grouped the 46 populations into 5 major genetic clusters that corresponded closely with geographic location following a classic isolation-by-distance model. Phylogenetic analysis of 14 mtDNA haplotypes revealed two reciprocally monophyletic groups. One of these groups (containing 4 geographically distant populations) was clearly identified by allozyme analysis and represents a distinct evolutionary unit. The remaining 4 allozyme groups were not distinguishable by mtDNA analysis. The evidence suggests that the populations within these groups derived from a small founding population that underwent rapid demographic expansion in ancient times. This was followed by more recent population bottlenecks resulting from habitat fragmentation associated with the widespread introduction of agriculture into the region. The generally low levels of allozyme and nucleotide diversity within these populations support this hypothesis. This revised version was published online in July 2006 with corrections to the Cover Date.     

Contrasting patterns of mitochondrial and microsatellite genetic structure among Western European populations of tawny owls (Strix aluco)

A recent study of mitochondrial phylogeography of tawny owls (Strix aluco) in western Europe suggested that this species survived the Pleistocene glaciations in three allopatric refugia located in Iberia, Italy, and the Balkans, and the latter was likely the predominant source of postglacial colonization of northern Europe. New data from seven microsatellite loci from 184 individual owls distributed among 14 populations were used to assess the genetic congruence between nuclear and mitochondrial DNA (mtDNA) markers. Microsatellites corroborated the major phylogeographical conclusions reached on the basis of the mtDNA sequences, but also showed important differences leading to novel inferences. Microsatellites corroborated the three major refugia and supported the Balkan origin of northern populations. When corrected for differences in effective population size, microsatellites and mtDNA yielded generally congruent overall estimates of population structure (N*ST=0.12 vs. RST=0.16); however, there was substantial heterogeneity in the RST among the seven nuclear loci that was not correlated with heterozygosity. Populations representing the Balkans postglacial expansion interact with populations from the other two refugia forming two clines near the Alps and the Pyrenees. In both cases, the apparent position of the contact zones differed substantially between markers due to the genetic composition of populations sampled in northern Italy and Madrid. Microsatellite data did not corroborate the lower genetic diversity of northern, recently populated regions as was found with mtDNA; this discrepancy was taken as evidence for a recent bottleneck recovery. Finally, this study suggests that congruence among genetic markers should be more likely in cases of range expansion into new areas than when populations interact across contact zones.     

Historical demography of freshwater mussels (Bivalvia: Unionidae): genetic evidence for population expansion and contraction during the late Pleistocene and Holocene

Genetic variation was examined in two endangered mussel species, Epioblasma brevidens and Epioblasma capsaeformis, and in a non‐listed species, Lampsilis fasciola, in the Clinch River, Tennessee, USA, by screening mitochondrial DNA (mtDNA) sequences and nuclear DNA microsatellites. Patterns of mtDNA polymorphism exhibited different trends in long‐term population sizes for each species during the late Pleistocene and Holocene (～20 000 ya to present); namely, E. brevidens has declined over time, E. capsaeformis has remained demographically stable, and L. fasciola has expanded. However, analyses using microsatellites did not exhibit similar trends, perhaps because homoplasy had eliminated long‐term population signatures for the loci examined. For both marker types, long‐term effective population size (N_e) was low in E. brevidens, intermediate in E. capsaeformis, and high in L. fasciola. Moderately diverged mtDNA lineages, perhaps indicative of secondary contact, were observed in E. brevidens and E. capsaeformis. Perhaps the most surprising result of this study was the high level of genetic variation observed at both mtDNA and microsatellite DNA markers for L. fasciola, variation seemingly contrary to the relatively small demes that currently reside in the Clinch River. However, the data are consistent with known demographic and life‐history traits of these three mussel species and their fish hosts, namely that they each use hosts with different dispersal capabilities, ranging from low, moderate, and high, respectively. The low divergence of mtDNA sequence variation reported in this and other recent mussel studies indicates that considerable extant population genetic variation probably originated during the late Pleistocene and Holocene. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 114 , 376–397.     

Genetic variation in original and colonizing Drosophila buzzatii populations analysed by microsatellite loci isolated with a new PCR screening method

A new polymerase chain reaction-based screening method for microsatellites is presented. Using this method, we isolated 12 microsatellite loci from Drosophila buzzatii, two of which were X-linked. We applied the other 10 microsatellite loci to the analysis of genetic variation in five natural populations of D. buzzatii. Two populations were from the species' original distribution in Argentina, whereas the other three were from Europe (two) and Australia that were colonized 200 and 65 years ago, respectively. Allelic variation was much larger in the original populations than in the colonizing ones and there was a tendency to decreased heterozygosity in the colonizing populations. We used three different statistical procedures for detecting population bottlenecks. All procedures suggested that the low variability in the populations in the Old World was not the result of the recent population decline, but was due to a founder effect followed by a population expansion. In fact, one procedure which detects population expansions and declines based on the genealogical history of microsatellite data suggested that an expansion had taken place in all the colonized populations.     

Genetic diversity and conservation in a small endangered horse population

The Old Kladruber horses arose in the 17th century as a breed used for ceremonial purposes. Currently, grey and black coat colour varieties exist as two sub-populations with different recent breeding history. As the population underwent historical bottlenecks and intensive inbreeding, loss of genetic variation is considered as the major threat. Therefore, genetic diversity in neutral and non-neutral molecular markers was examined in the current nucleus population. Fifty microsatellites, 13 single nucleotide polymorphisms (SNPs) in immunity-related genes, three mutations in coat colour genes and one major histocompatibility (MHC-DRA) gene were studied for assessing genetic diversity after 15 years of conservation. The results were compared to values obtained in a similar study 13 years ago. The extent of genetic diversity of the current population was comparable to other breeds, despite its small size and isolation. The comparison between 1997 and 2010 did not show differences in the extent of genetic diversity and no loss of allele richness and/or heterozygosity was observed. Genetic differences identified between the black and grey sub-populations observed 13 years ago persisted. Deviations from the Hardy–Weinberg equilibrium found in 19 microsatellite loci and in five SNP loci are probably due to selective breeding. No differences between neutral and immunity-related markers were found. No changes in the frequencies of markers associated with two diseases, melanoma and insect bite hypersensitivity, were observed, due probably to the short interval of time between comparisons. It, thus, seems that, despite its small size, previous bottlenecks and inbreeding, the molecular variation of Old Kladruber horses is comparable to other horse breeds and that the current breeding policy does not compromise genetic variation of this endangered population.     

The evolutionary dynamics of the lion Panthera leo revealed by host and viral population genomics

The lion Panthera leo is one of the world's most charismatic carnivores and is one of Africa's key predators. Here, we used a large dataset from 357 lions comprehending 1.13 megabases of sequence data and genotypes from 22 microsatellite loci to characterize its recent evolutionary history. Patterns of molecular genetic variation in multiple maternal (mtDNA), paternal (Y-chromosome), and biparental nuclear (nDNA) genetic markers were compared with patterns of sequence and subtype variation of the lion feline immunodeficiency virus (FIV(Ple)), a lentivirus analogous to human immunodeficiency virus (HIV). In spite of the ability of lions to disperse long distances, patterns of lion genetic diversity suggest substantial population subdivision (mtDNA Phi(ST) = 0.92; nDNA F(ST) = 0.18), and reduced gene flow, which, along with large differences in sero-prevalence of six distinct FIV(Ple) subtypes among lion populations, refute the hypothesis that African lions consist of a single panmictic population. Our results suggest that extant lion populations derive from several Pleistocene refugia in East and Southern Africa ( approximately 324,000-169,000 years ago), which expanded during the Late Pleistocene ( approximately 100,000 years ago) into Central and North Africa and into Asia. During the Pleistocene/Holocene transition ( approximately 14,000-7,000 years), another expansion occurred from southern refugia northwards towards East Africa, causing population interbreeding. In particular, lion and FIV(Ple) variation affirms that the large, well-studied lion population occupying the greater Serengeti Ecosystem is derived from three distinct populations that admixed recently.     

Limited genetic structure and evidence for dispersal among populations of the endangered Florida grasshopper sparrow, Ammodramus savannarum floridanus

The Florida grasshopper sparrow, Ammodramus savannarum floridanus, is a non-migratory, endangered subspecies endemic to the prairie region of south-central Florida. It has experienced significant population declines and is currently restricted to five locations. We found substantial levels of variation in microsatellites and mtDNA control region sequences, estimates of inbreeding genetic effective population sizes that were much larger than the estimated census size, and no evidence of inbreeding within five sampled populations (n = 105). We also found a lack of genetic structure among populations (F _ST = 0.0123 for microsatellites and θ = 0.008 for mtDNA), and evidence for dispersal between populations, with 7.6% of all individuals identified as immigrants to their population of capture. We suggest that the subspecies be managed as a single management unit on a regional scale rather than as multiple management units on a local subpopulation scale. There is still a limited opportunity to preserve much of the present genetic variation in this subspecies, if immediate measures are taken to reverse the current population decline before this variation is reduced by genetic drift.     

Mitochondrial DNA and human evolution

For the past seven years or so, much discussion and controversy in the field of human evolution has revolved around the application and interpretation of studies of human mitochondrial DNA variation, particularly the hypothesis that all mtDNA types in contemporary populations can be traced back to a single African ancestor who lived about 200,000 years ago. In this review I describe the evidence that led to this hypothesis, subsequent work, and where things stand now, particularly with respect to recent criticisms concerning the adequacy of phylogenetic analyses of the mtDNA data. I also describe a new method of analyzing mtDNA data that suggests that all human populations underwent a dramatic expansion some 40,000 years ago, possibly in association with revolutionary advances in human behavior, as well as an important implication of population expansions for mtDNA disease studies.     

Population genetics of the endangered Wattled Curassow (Crax globulosa,Cracidae, Aves) of the Colombian–Peruvian Amazon using DNA microsatellites and ND2 mitochondrial sequences

The Wattled Curassow (Crax globulosa, Cracidae, Aves) is a large bird living in the Western Amazon basin and a critically endangered species in the Colombian and in the Peruvian Amazon. We carried out the first population genetics analysis of this species employing six nuclear microsatellite markers and sequences of the mtND2 gene. The main results are as follows. (1) The levels of gene diversity were high for the overall population as well as for each of the three islands for both microsatellites and mtDNA. (2) A small amount of genetic differentiation among populations was found with both types of markers (F_ST = 0.027 for microsatellites and N_ST = 0.17 for mitochondrial sequences). (3) Using microsatellites, the Geneclass 2.0 software detected a low correct assignment of individuals to their respective populations. The Structure software only detected one gene pool for the entire area studied. These results are relevant for conservation efforts of this critically endangered species.     

Relative information content of polymorphic microsatellites and mitochondrial DNA for inferring dispersal and population genetic structure in the olive sea snake, Aipysurus laevis

Polymorphic microsatellites are widely considered more powerful for resolving population structure than mitochondrial DNA (mtDNA) markers, particularly for recently diverged lineages or geographically proximate populations. Weaker population subdivision for biparentally inherited nuclear markers than maternally inherited mtDNA may signal male-biased dispersal but can also be attributed to marker-specific evolutionary characteristics and sampling properties. We discriminated between these competing explanations with a population genetic study on olive sea snakes, Aipysurus laevis. A previous mtDNA study revealed strong regional population structure for A. laevis around northern Australia, where Pleistocene sea-level fluctuations have influenced the genetic signatures of shallow-water marine species. Divergences among phylogroups dated to the Late Pleistocene, suggesting recent range expansions by previously isolated matrilines. Fine-scale population structure within regions was, however, poorly resolved for mtDNA. In order to improve estimates of fine-scale genetic divergence and to compare population structure between nuclear and mtDNA, 354 olive sea snakes (previously sequenced for mtDNA) were genotyped for five microsatellite loci. F statistics and Bayesian multilocus genotype clustering analyses found similar regional population structure as mtDNA and, after standardizing microsatellite F statistics for high heterozygosities, regional divergence estimates were quantitatively congruent between marker classes. Over small spatial scales, however, microsatellites recovered almost no genetic structure and standardized F statistics were orders of magnitude smaller than for mtDNA. Three tests for male-biased dispersal were not significant, suggesting that recent demographic expansions to the typically large population sizes of A. laevis have prevented microsatellites from reaching mutation-drift equilibrium and local populations may still be diverging.     

Genetic structure and signature of population decrease in the critically endangered freshwater cyprinid <Emphasis Type="Italic">Chondrostoma lusitanicum</Emphasis>

The endemic and critically endangered cyprinid Chondrostoma lusitanicum has a very restricted distribution range. In order to estimate genetic diversity, characterize population structure and infer the demographic history, we examined six microsatellite loci and cytochrome b (mtDNA) sequences from samples taken throughout C. lusitanicum’s geographical range. Estimates of genetic diversity were low in all samples (average He < 0.35). The microsatellite data pointed to a major difference between northern (Samarra and Tejo drainages) and southern (Sado and Sines drainages) samples. This separation was not so clear with mtDNA, since one sample from the Tejo drainage grouped with the southern samples. This could be related with ancestral polymorphism or with admixture events between northern and southern sites during the late Pleistocene. Nevertheless, both markers indicate high levels of population differentiation in the north (for microsatellites F _ST >  0.23; and for mtDNA Φ_ST > 0.74) and lower levels in the south (F _ST < 0.05; Φ_ST < 0.40). With microsatellites we detected strong signals of a recent population decrease in effective size, by more than one order of magnitude, starting in the last centuries. This is consistent with field observations reporting a severe anthropogenic-driven population decline in the last decades. On the contrary mtDNA suggested a much older expansion. Overall, these results suggest that the distribution of genetic diversity in C. lusitanicum is the result of both ancient events related with drainage system formation, and recent human activities. The potential effect of population substructure generating genetic patterns similar to a population decrease is discussed, as well as the implications of these results for the conservation of C. lusitanicum. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

MITOCHONDRIAL DNA VARIATION IN THE FIELD VOLE (MICROTUS AGRESTIS): REGIONAL POPULATION STRUCTURE AND COLONIZATION HISTORY

Restriction fragment length polymorphism analysis of mitochondrial DNA (mtDNA) was used to examine the genetic structure among field voles (Microtus agrestis) from southern and central Sweden. A total of 57 haplotypes was identified in 158 voles from 60 localities. Overall mtDNA diversity was high, but both haplotype and nucleotide diversity exhibited pronounced geographic heterogeneity. Phylogenetic analyses revealed a shallow tree with seven primary mtDNA lineages separated by sequence divergences ranging from 0.6% to 1.0%. The geographic structure of mtDNA diversity and lineage distribution was complex but strongly structured and deviated significantly from an equilibrium situation. The extensive mtDNA diversity observed and the recent biogeographic history of the region suggests that the shallow mtDNA structure in the field vole cannot be explained solely by stochastic lineage sorting in situ or isolation by distance. Instead, the data suggest that the genetic imprints of historical demographic conditions and vicariant geographic events have been preserved and to a large extent determine the contemporary geographic distribution of mtDNA variation. A plausible historical scenario involves differentiation of mtDNA lineages in local populations in glacial refugia, a moving postglacial population structure, and bottlenecks and expansions of mtDNA lineages during the postglacial recolonization of Sweden. By combining the mtDNA data with an analysis of Y-chromosome variation, a specific population unit was identified in southwestern Sweden. This population, defined by a unique mtDNA lineage and fixation of a Y-chromosome variant, probably originated in a population bottleneck in southern Sweden about 12,000 to 13,000 calendar years ago.     

Homing behaviour facilitates subtle genetic differentiation among river populations of Alosa sapidissima: microsatellites and mtDNA

Significant but subtle differentiation was detected for both microsatellite DNA and mitochondrial DNA among four populations of American shad Alosa sapidissima . The data indicate that straying among rivers is sufficient to permit only marginal population differentiation in this species, but suggest that individual river populations should be managed as distinct stocks. Comparison of the Hudson and Columbia populations, the latter derived from the former over 100 years ago, revealed only a slight reduction in microsatellite DNA variation for the founded population but halving of mitochondrial DNA, consistent with the haploid maternal inheritance of the latter marker. The depleted and endangered James River (Virginia) population and two other Atlantic coast populations exhibited similar levels of microsatellite DNA variation, but mtDNA diversity in the James River was marginally lower than in other Atlantic populations, again consistent with the low effective population size of mtDNA.     

Exploring neutral and adaptive processes in expanding populations of gilthead sea bream, Sparus aurata L., in the North-East Atlantic

Recent studies in empirical population genetics have highlighted the importance of taking into account both neutral and adaptive genetic variation in characterizing microevolutionary dynamics. Here, we explore the genetic population structure and the footprints of selection in four populations of the warm-temperate coastal fish, the gilthead sea bream (Sparus aurata), whose recent northward expansion has been linked to climate change. Samples were collected at four Atlantic locations, including Spain, Portugal, France and the South of Ireland, and genetically assayed using a suite of species-specific markers, including 15 putatively neutral microsatellites and 23 expressed sequence tag-linked markers, as well as a portion of the mitochondrial DNA (mtDNA) control region. Two of the putatively neutral markers, Bld-10 and Ad-10, bore signatures of strong directional selection, particularly in the newly established Irish population, although the potential 'surfing effect' of rare alleles at the edge of the expansion front was also considered. Analyses after the removal of these loci suggest low but significant population structure likely affected by some degree of gene flow counteracting random genetic drift. No signal of historic divergence was detected at mtDNA. BLAST searches conducted with all 38 markers used failed to identify specific genomic regions associated to adaptive functions. However, the availability of genomic resources for this commercially valuable species is rapidly increasing, bringing us closer to the understanding of the interplay between selective and neutral evolutionary forces, shaping population divergence of an expanding species in a heterogeneous milieu.