Late glacial history of the cold-adapted freshwater fish Cottus gobio,revealed by microsatellites

The distribution of genetic diversity at 10 highly polymorphic microsatellite loci within the European freshwater fish, Cottus gobio, L. was examined. The sampling range comprised a large geographical scale including lineages known to be highly divergent at both mitochondrial DNA (mtDNA) and allozymes. An analysis of genetic variability within populations showed that expected heterozygosity and allelic richness could be explained largely by current effective population sizes. Evidence was found, however, that historical processes predating the last major glaciation affected allelic richness. In addition to confirming the large-scale patterns from earlier studies, the microsatellite data revealed new insights into recent processes by analysing genetic structure within ancient lineages defined by mtDNA data. Stepwise mutation model (SMM) and nonSMM-based methods demonstrated a clear genetic structuring within the Northwestern European lineage comprising populations from Britain and Belgium, and within the Central European lineage populations from the rivers Danube, Elbe and Main. Supported by an analysis of genetic variability within populations these results showed that the bullhead populations most probably persisted throughout the last major glaciation within the British Isles and within the drainages of the rivers Elbe and Main. Such observations provide the first genetic evidence for a glacial refugium in such close proximity to the European glacial margins.     

Reduced genetic variation and strong genetic population structure in the freshwater killifish Valencia letourneuxi (Valenciidae) based on nuclear and mitochondrial markers

The genetic variation of the critically endangered Corfu killifish (Valencia letourneuxi), an endemic freshwater fish species of the western Balkans, was assessed for nine populations sampled in eight water systems in western continental Greece, the Peloponnese and the Ionian Island of Corfu, using mitochondrial and microsatellite markers. The analyses were based on data from three mtDNA regions (D‐loop, COI and 16S rRNA sequences) and 14 microsatellite loci. Samples from the congeneric species Valencia hispanica and the phylogenetically closely related species Aphanius fasciatus were also used in the study as outgroups. Both the mitochondrial and the microsatellite analyses revealed three distinct population groupings associated with the geographical distribution of the populations: one southern group occupying rivers draining to the Patraikos Gulf, the second one including the populations flowing into the Amvrakikos Gulf and the third, more northern group, including the other populations from rivers in Corfu Island and Epirus flowing into the Ionian Sea. Within these groupings there is limited genetic differentiation between populations; in addition, there is reduced intrapopulation genetic variation, evidenced by low heterozygosity values, number of alleles and haplotype diversity. In terms of taxonomic implications and appropriate management actions for conservation, our data suggest that the major population groups should be regarded at least as three distinct conservation units (CUs), with translocation and restocking actions to take place only within the geographical range of the CU concerned. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 334–349.     

Molecular Markers Reveal Limited Population Genetic Structure in a North American Corvid,Clark’s Nutcracker (Nucifraga columbiana)

The genetic impact of barriers and Pleistocene glaciations on high latitude resident species has not been widely investigated. The Clark’s nutcracker is an endemic North American corvid closely associated with Pinus-dominated forests. The nutcracker’s encompasses known barriers to dispersal for other species, and glaciated and unglaciated areas. Clark’s nutcrackers also irruptively disperse long distances in search of pine seed crops, creating the potential for gene flow among populations. Using the highly variable mitochondrial DNA control region, seven microsatellite loci, and species distribution modeling, we examined the effects of glaciations and dispersal barriers on population genetic patterns and population structure of nutcrackers. We sequenced 900 bp of mitochondrial control region for 169 individuals from 15 populations and analysed seven polymorphic microsatellite loci for 13 populations across the Clark’s nutcracker range. We used species distribution modeling and a range of phylogeographic analyses to examine evolutionary history. Clark’s nutcracker populations are not highly differentiated throughout their range, suggesting high levels of gene flow among populations, though we did find some evidence of isolation by distance and peripheral isolation. Our analyses suggested expansion from a single refugium after the last glacial maximum, but patterns of genetic diversity and paleodistribution modeling of suitable habitat were inconclusive as to the location of this refugium. Potential barriers to dispersal (e.g. mountain ranges) do not appear to restrict gene flow in Clark’s nutcracker, and postglacial expansion likely occurred quickly from a single refugium located south of the ice sheets.     

Unique mitochondrial DNA lineages in Irish stickleback populations: cryptic refugium or rapid recolonization?

Repeated recolonization of freshwater environments following Pleistocene glaciations has played a major role in the evolution and adaptation of anadromous taxa. Located at the western fringe of Europe, Ireland and Britain were likely recolonized rapidly by anadromous fishes from the North Atlantic following the last glacial maximum (LGM). While the presence of unique mitochondrial haplotypes in Ireland suggests that a cryptic northern refugium may have played a role in recolonization, no explicit test of this hypothesis has been conducted. The three‐spined stickleback is native and ubiquitous to aquatic ecosystems throughout Ireland, making it an excellent model species with which to examine the biogeographical history of anadromous fishes in the region. We used mitochondrial and microsatellite markers to examine the presence of divergent evolutionary lineages and to assess broad‐scale patterns of geographical clustering among postglacially isolated populations. Our results confirm that Ireland is a region of secondary contact for divergent mitochondrial lineages and that endemic haplotypes occur in populations in Central and Southern Ireland. To test whether a putative Irish lineage arose from a cryptic Irish refugium, we used approximate Bayesian computation (ABC). However, we found no support for this hypothesis. Instead, the Irish lineage likely diverged from the European lineage as a result of postglacial isolation of freshwater populations by rising sea levels. These findings emphasize the need to rigorously test biogeographical hypothesis and contribute further evidence that postglacial processes may have shaped genetic diversity in temperate fauna.     

Phylogeography of Potamon ibericum (Brachyura: Potamidae) identifies Quaternary glacial refugia within the Caucasus biodiversity hot spot

Refugia are critical for the maintenance of biodiversity during the periods of Quaternary climatic oscillations. The long‐term persistence of refugial populations in a large continuous refugium has resulted in a homogenous pattern of genetic structure among populations, while highly structured evolutionary lineages characterize the restriction of refugial populations to smaller subrefugia. These mechanisms have resulted in the identification of hot spots of biodiversity within putative glacial refugia. We studied phylogeography of Potamon ibericum (Brachyura: Potamidae) in the drainages of the western Caucasus biodiversity hot spot (i.e., Colchis and the Caucasus) to infer spatial genetic structure and potential refugia for a freshwater crab in this region. These areas have traditionally considered as a refugium due to the presence of Tertiary relict species. We integrated population genetic data and historical demographic analysis from cytochrome oxidase subunit I sequences and paleoclimatic data from species distribution modeling (SDM). The results revealed the lack of phylogeographic structure and provided evidence for demographic expansion. The SDM presented a rather homogenous and large refugium that extended from northeast Turkey to Colchis during the last glacial period. In contrast to these findings, previous phylogeographic study on P. ibericum of the eastern Caucasus biodiversity hot spot (i.e., Hyrcania) identified multiple independent refugia. By combining these results, we explain the significance of this important western Palearctic hot spot of biological diversity in shaping the geographic distribution of intraspecific genetic diversity in a freshwater taxon.     

Mitochondrial DNA variation in European populations of Silurus glanis

Mitochondrial DNA diversity of 13 wild Silurus glanis populations (covering the entire range of the species) and eight hatchery populations was investigated. PCR-RFLP analysis of four regions of mitochondrial DNA (cytochrome b, control region, ND-5/6) was used. Nineteen haplotypes were found. Thirteen of them were private. The proportion of total genetic diversity attributable to population differentiation was almost 80%. Despite the existence of significant differentiation between populations for mtDNA variation, no consistent pattern of geographic structuring was revealed and nucleotide divergence among S. glanis populations was low. These phenomena are discussed with regard to the impact of glaciation events. The domesticated stocks show less genetic diversity than natural ones, possibly due to their mode of management. Analysis of three European catfish species S. glanis, S. aristotelis and Silurus triostegus (sampled in the Euphrates river) revealed several endonucleases which produced restriction phenotypes diagnostic for the three species.     

Unique genetic variation at a species' rear edge is under threat from global climate change

Global climate change is having a significant effect on the distributions of a wide variety of species, causing both range shifts and population extinctions. To date, however, no consensus has emerged on how these processes will affect the range-wide genetic diversity of impacted species. It has been suggested that species that recolonized from low-latitude refugia might harbour high levels of genetic variation in rear-edge populations, and that loss of these populations could cause a disproportionately large reduction in overall genetic diversity in such taxa. In the present study, we have examined the distribution of genetic diversity across the range of the seaweed Chondrus crispus, a species that has exhibited a northward shift in its southern limit in Europe over the last 40 years. Analysis of 19 populations from both sides of the North Atlantic using mitochondrial single nucleotide polymorphisms (SNPs), sequence data from two single-copy nuclear regions and allelic variation at eight microsatellite loci revealed unique genetic variation for all marker classes in the rear-edge populations in Iberia, but not in the rear-edge populations in North America. Palaeodistribution modelling and statistical testing of alternative phylogeographic scenarios indicate that the unique genetic diversity in Iberian populations is a result not only of persistence in the region during the last glacial maximum, but also because this refugium did not contribute substantially to the recolonization of Europe after the retreat of the ice. Consequently, loss of these rear-edge populations as a result of ongoing climate change will have a major effect on the overall genetic diversity of the species, particularly in Europe, and this could compromise the adaptive potential of the species as a whole in the face of future global warming.     

Genetic diversity and population structure of the northern snakehead (<Emphasis Type="Italic">Channa argus</Emphasis> Channidae: Teleostei) in central China: implications for conservation and management

Major threats to freshwater fish diversity now include loss of native genetic diversity as a consequence of translocations of fishes between sites and from hatcheries to sites, and small effective population sizes resulting from overfishing and/or habitat loss. Ten polymorphic microsatellite markers were employed to evaluate genetic diversity, population genetic structure and gene flow amongst nine populations of the ecologically and economically important fish, the northern snakehead (Channa argus), in three river systems in central China. Multiple analyses revealed evidence of high genetic diversity and pronounced subdivision based on both regional separation and on river systems. A lack of evidence of genetic bottleneck over recent generations was consistent with the long-term stability of population size and contemporary distribution. The effective population sizes for most C. argus populations were small, suggesting the need for future conservation efforts focusing on these populations. Different lines of evidence point to the local enhancement of stocks by both aquaculture-reared fish and the transfer of wild fish. This study illustrates how human activities may affect genetic diversity and population genetic structure of C. argus populations, and highlights the need for new management regimes to protect native freshwater fish genetic diversity.     

Genetic admixture of burbot (Teleostei: Lota lota) in Lake Constance from two European glacial refugia

The burbot, Lota lota, is the only freshwater species of the codfish family and has a Holarctic distribution. Pleistocene glaciations caused significant geographical differentiation in the past, but its life history characterized by winter spawning migrations over large distances is likely to homogenize populations by contemporary gene flow. We investigated the population genetic structure of 541 burbots from Lake Constance and adjacent Rhine and Danube tributaries in Europe using the entire mitochondrial DNA (mtDNA) control region and 11 microsatellites. Microsatellites revealed considerable population divergence (F(ST) = 0.26) and evidenced recent bottlenecks in two Central European rivers. In accordance to previous evidence two main phylogeographic lineages (Atlantic and Danubian) were found co-occurring at similar frequencies in Lake Constance, where they currently undergo random mating as indicated by microsatellites. The Danubian lineage contributed only a small proportion to the lake's mtDNA diversity, and probably expanded within the lake shortly after its formation approximately 10,000-15,000 BP. The larger Atlantic haplotype diversity suggested a population expansion older than the lake itself. Levels of admixture at microsatellite loci were less obvious due to their high variability, and coalescence methods were used to estimate past admixture proportions. Our results reinforce a model of a two-step colonization of Europe by burbot from an ancestral Danubian refuge, and confirm the persistence of a secondary Atlantic refuge, as proposed to exist for other freshwater fish. We conclude that the present-day burbot population in Lake Constance bears the genetic signature of both contemporary gene flow and historical separation events.     

Mitochondrial and nuclear DNA phylogeography of European grayling (Thymallus thymallus): evidence for secondary contact zones in central Europe

Mitochondrial and microsatellite DNA markers were applied to infer the phylogeography, intraspecific diversity and dynamics of the distributional history of European grayling (Thymallus thymallus) with focus on its central and northern European distribution range. Phylogenetic and nested clade analyses revealed at least four major mtDNA lineages, which evolved in geographical isolation during the Pleistocene. These lineages should be recognized as the basic evolutionary significant units (ESUs) for grayling in central and northern Europe. In addition, and in contrast to previous work on grayling, the results of Bayesian analysis of individual admixture coefficients, two-dimensional scaling analysis and spatial analysis of molecular variance provided evidence for a high level of admixture among major lineages in contact zones between drainages (e.g. the low mountain range of Germany), most likely resulting from glacial perturbations and ancient river connections between drainages during the Pleistocene glaciations. Even within river systems, a high level of differentiation among populations was revealed as indicated by the microsatellite data. Grayling sampled from 29 sites displayed high levels of differentiation (overall F(ST) = 0.367), a high number of private alleles and high bootstrap support for the genetic distance-based population clusters across 12 loci. We specifically discuss our results in context of phylogeograpic studies on other European freshwater fish species with habitat preferences similar to those of grayling. Our study shows that both large-scale phylogeographical and detailed genetic analyses on a fine scale are mandatory for developing appropriate conservation guidelines of endangered species.     

Population genetic analysis reveals a geographically limited transition zone between two genetically distinct Atlantic salmon lineages in Norway

Atlantic salmon is characterized by a high degree of population genetic structure throughout its native range. However, while populations inhabiting rivers in Norway and Russia make up a significant proportion of salmon in the Atlantic, thus far, genetic studies in this region have only encompassed low to modest numbers of populations. Here, we provide the first “in‐depth” investigation of population genetic structuring in the species in this region. Analysis of 18 microsatellites on >9,000 fish from 115 rivers revealed highly significant population genetic structure throughout, following a hierarchical pattern. The highest and clearest level of division separated populations north and south of the Lofoten region in northern Norway. In this region, only a few populations displayed intermediate genetic profiles, strongly indicating a geographically limited transition zone. This was further supported by a dedicated cline analysis. Population genetic structure was also characterized by a pattern of isolation by distance. A decline in overall genetic diversity was observed from the south to the north, and two of the microsatellites showed a clear decrease in number of alleles across the observed transition zone. Together, these analyses support results from previous studies, that salmon in Norway originate from two main genetic lineages, one from the Barents–White Sea refugium that recolonized northern Norwegian and adjacent Russian rivers, and one from the eastern Atlantic that recolonized the rest of Norway. Furthermore, our results indicate that local conditions in the limited geographic transition zone between the two observed lineages, characterized by open coastline with no obvious barriers to gene flow, are strong enough to maintain the genetic differentiation between them.     

Microsatellite data resolve phylogeographic patterns in European grayling, Thymallus thymallus, Salmonidae

The phylogeography of an endangered salmonid, European grayling (Thymallus thymallus), was studied based on analysis of 17 nuclear microsatellite DNA loci. In agreement with earlier mitochondrial DNA (mtDNA) studies, phylogenetic relationships of the populations suggested that northern Europe was colonized from two distinct Pleistocene refugia. Furthermore, microsatellites revealed highly supported grouping of mainland Swedish, Norwegian, Danish, German and Slovenian populations, suggesting that grayling from the northwestern and central Europe have descended from their southern conspecifics. The level of divergence between populations was substantial, even across short geographical distances. Although this was in part due to postglacial colonization patterns and contemporary barriers for gene flow, the high divergence estimates between hydrologically connected sampling locations implied efficient interpopulation reproductive isolation. Microsatellites revealed that the populations exhibited, on average, only 3.5 (+/-2.2) alleles per locus, indicating that T. thymallus has strikingly low levels of intrapopulation genetic diversity as compared with other freshwater fish species. Accordingly, as indicated by analysis of molecular variance (AMOVA), only 49.1-58.0% of the total grayling microsatellite diversity resided within populations. A latitudinal genetic diversity gradient, potentially resulting from glaciation-mediated founder events, was not evident. Alternatively, it is possible that grayling display limited dispersal behaviour/capability, leading to low long-term effective population sizes and, consequently, depauperate intrapopulation polymorphism. These findings have implications for conservation of T. thymallus. Importantly, they exemplify that microsatellites can be highly informative for intraspecific phylogeography studies dealing with substantial divergence scales.     

Migratory behaviour shapes spatial genetic structure of cyprinid fishes within the Lake Malawi catchment

• Genetic differences among freshwater fish populations are dependent on life‐history characteristics of the species, including the range of adult dispersal and the extent of homing to natal breeding grounds. However, the effects of variation in such characteristics on population genetic connectivity are rarely studied comparatively among closely related species.
• We studied population genetic structure within three congeneric cyprinid species from the Lake Malawi catchment that differ substantially in life‐history traits and conservation status, using a combination of microsatellite and mitochondrial DNA markers. Mpasa (Opsaridium microlepis) is a large (70 cm total length) migratory species that spawns in rivers, but as an adult is exclusively known from the main lake body. Sanjika (Opsaridium microcephalum), is a medium size (30 cm total length) species that exists in lake breeding, river‐lake migratory and apparently landlocked populations. Dwarf sanjika (Opsaridium tweddleorum) is a small non‐migratory species (15 cm total length) that persists in small tributaries surrounding the main lake and adjoining rivers.
• The results revealed striking differences among the three species in spatial genetic structuring. The river‐lake migratory mpasa showed only weak yet significant population genetic structure within the main Lake Malawi catchment, suggesting that there is no strong natal homing. The habitat‐generalist sanjika showed only weak spatial genetic differentiation at microsatellite loci within the Lake Malawi catchment, but moderate structure in mitochondrial DNA, potentially reflecting male‐biased dispersal. The river‐restricted dwarf sanjika showed strong genetic structure in both microsatellite and mitochondrial DNA, suggesting strictly limited dispersal at both adult and juvenile stages.
• We conclude that contrasting migration life histories have resulted in dramatically different patterns of population genetic structure among these congeneric species. The observed patterns demonstrate how divergent life‐history evolution may strongly influence broader patterns of population genetic connectivity in freshwater fish, with consequences for management and conservation. Specifically the results suggesting gene flow among Lake Malawi populations of mpasa, an IUCN red‐listed ‘Endangered’ species endemic to the lake catchment, imply that conservation initiatives operating at both local and catchment scales are needed to reverse local population decline.

     

Colonization, extinction, and phylogeographic patterning in a freshwater crustacean

Phylogeographic analyses have revealed the importance of Pleistocene vicariance events in shaping the distribution of genetic diversity in freshwater fishes. However, few studies have examined the patterning of variation in freshwater organisms with differing dispersal syndromes and life histories. The present investigation addresses this gap, examining the phylogeography of Sida crystallina, a species whose production of diapausing eggs capable of passive dispersal was thought to constrain its regional genetic differentiation. By contrast, the present analysis has revealed deep allozyme and cytochrome oxidase I mitochondrial DNA divergence between populations from North America and Europe. Moreover, North American populations are separated into four allopatric assemblages, whose distribution suggests their derivation from different Pleistocene refugia. These lineages show higher haplotype diversity and deeper sequence divergence than those of any fish from temperate North America. Its distinctive life history traits have evidently sheltered lineages of Sida from extinction, contributing to a remarkably comprehensive and high resolution phylogeographic record.     

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.     

Contemporary factors influencing genetic diversity in the Alaska humpback whitefish Coregonus clupeaformis complex

Thirteen microsatellite loci were used to address three hypotheses regarding genetic diversity in the humpback whitefish Coregonus clupeaformis complex in Alaska. The test results provided further insight into the factors influencing C. clupeaformis complex population structure and level of genetic variation. First, the microsatellite data did not provide evidence of two spatially distinct Beringian and Eurasian refugial groups as revealed in previous phylogeographic analyses of mitochondrial DNA variation. Rather, the population structure inferred from the microsatellite variation appears to reveal the influence of factors on a more recent time scale, including gene flow among the refugial groups and isolation of some anadromous and freshwater resident populations. Second, anadromous C. clupeaformis complex collections exhibited higher intra‐population genetic diversity than freshwater resident collections. This outcome is consistent with previous meta analyses suggesting that freshwater resident populations probably have smaller historical effective population sizes and less conspecific gene flow because the habitat tends to be smaller and supports fewer and smaller populations. Finally, the analysis of contemporary immigration rates was consistent with, but did not provide statistical support for, the hypothesis that gene flow among anadromous C. clupeaformis complex populations along coastal Alaska is influenced by the Alaska Coastal Current. Further studies are needed to evaluate gene flow among coastal Alaska C. clupeaformis complex populations.     

Genetic variation in a freshwater prawn species,Palaemon paucidens,in South Korea

The biogeography of freshwater biota in the Korean Peninsula has been affected by recent geological processes and anthropogenic activity. The freshwater prawn, Palaemon paucidens, can serve as a non-fish model organism suitable for assessing these factors, as it is found in all river systems in the Korean Peninsula and may have been introduced by humans in some regions. In this study, we investigated the geographical distribution of genetic variation and the genetic structure of P. paucidens populations using mitochondrial DNA sequences and genotypes identified from four microsatellite loci. Our results showed that populations from westward-flowing river systems that drain into the Yellow Sea have more genetic diversity than those from southward-flowing river systems, and that the highest genetic variance revealed by analysis of molecular variance (AMOVA) using both genetic markers was observed in river systems grouped as HAN + GEUM, NAKDONG + JEJU, YOUNGSAN, and SEOMJIN. These results suggest that HAN and GEUM, in which freshwater prawn populations have higher levels of genetic diversity, were the most recently isolated river systems from Asian continental systems. Therefore, populations from HAN and GEUM experienced less severe bottlenecks than those from YOUNGSAN and SEOMJIN.     

Multiple parasites are driving major histocompatibility complex polymorphism in the wild

Parasite mediated selection may result in arms races between host defence and parasite virulence. In particular, simultaneous infections from multiple parasite species should cause diversification (i.e. balancing selection) in resistance genes both at the population and the individual level. Here, we tested these ideas in highly polymorphic major histocompatibility complex (MHC) genes from three-spined sticklebacks (Gasterosteus aculeatus L.). In eight natural populations, parasite diversity (15 different species), and MHC class IIB diversity varied strongly between habitat types (lakes vs. rivers vs. estuaries) with lowest values in rivers. Partial correlation analysis revealed an influence of parasite diversity on MHC class IIB variation whereas general genetic diversity assessed at seven microsatellite loci was not significantly correlated with parasite diversity. Within individual fish, intermediate, rather than maximal allele numbers were associated with minimal parasite load, supporting theoretical models of self-reactive T-cell elimination. The optimal individual diversity matched those values female fish try to achieve in their offspring by mate choice. We thus present correlative evidence supporting the 'allele counting' strategy for optimizing the immunocompetence in stickleback offspring.     

Threatened freshwater pearl mussel <Emphasis Type="Italic">Margaritifera margaritifera</Emphasis> L. in NW Spain: low and very structured genetic variation in southern peripheral populations assessed using microsatellite markers

A genetic analysis of freshwater pearl mussel Margaritifera margaritifera populations from NW Spain, a peripheral area of its European distribution, was carried out using microsatellite markers. These populations were formerly reported as genetically differentiated on the basis of growth and longevity studies. Ten loci previously characterized in populations from central Europe were used to comparatively analyze the genetic variability at the southern edge of the species’ range. Iberian pearl mussel populations showed very low genetic variability and significant high genetic differentiation. Half of the total genetic diversity observed appeared to be distributed between populations, which suggested a highly structured adaptive potential in pearl mussel at the southern peripheral distribution of the species. Population distinctiveness was evidenced by assignment tests, which revealed a high accuracy of individual assignments to their population of origin. All data suggested low effective population size and major effects of genetic drift on population genetic structure. In order to avoid further loss of genetic variation in biologically distinctive populations from NW Spain, prioritization of genetic resources of this species is required for conservation and management.     

Historic speciation and recent colonization of Eurasian monkey gobies (Neogobius fluviatilis and N. pallasi) revealed by DNA sequences,microsatellites, and morphology

Aim Hidden diversity within an invasive ‘species’ can mask both invasion pathways and confound management goals. We assessed taxonomic status and population structure of the monkey goby Neogobius fluviatilis across Eurasia, comparing genetic variation across its native and invasive ranges. Location Native populations were analysed within the Black and Caspian Sea basins, including major river drainages (Dnieper, Dniester, Danube, Don and Volga rivers), along with introduced locations within the upper Danube and Vistula river systems. Methods DNA sequences and 10 nuclear microsatellite loci were analysed to test genetic diversity and divergence patterns of native and introduced populations; phylogenetic analysis of mtDNA cytochrome b and nuclear RAG‐1 sequences assessed taxonomic status of Black and Caspian Sea lineages. Multivariate analysis of morphology was used to corroborate phylogenetic patterns. Population genetic structure within each basin was evaluated with mtDNA and microsatellite data using F_ST analogues and Bayesian assignment tests. Results Phylogenetic analysis of mitochondrial and nuclear sequences discerned a pronounced genetic break between monkey gobies in the Black and Caspian Seas, indicating a long‐term species‐level separation dating to c. 3 million years. This pronounced separation further was confirmed from morphological and population genetic divergence. Bayesian inference showed congruent patterns of population structure within the Black Sea basin. Introduced populations in the Danube and Vistula River basins traced to north‐west Black Sea origins, a genetic expansion pattern matching that of other introduced Ponto‐Caspian gobiids. Main conclusions Both genetic and morphological data strongly supported two species of monkey gobies that were formerly identified as subspecies: N. fluviatilis in the Black Sea basin, Don and Volga Rivers, and the Kumo‐Manych Depression, and Neogobius pallasi in the Caspian Sea and Volga River delta. Genetic origins of introduced N. fluviatilis populations indicated a common invasion pathway shared with other introduced Ponto‐Caspian fishes and invertebrates.