Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Population Genetic Patterns of Threatened European Mudminnow (Umbra krameri Walbaum, 1792) in a Fragmented Landscape: Implications for Conservation Management

  • Péter Takács ,

    takacs.peter@okologia.mta.hu

    Affiliation Balaton Limnological Institute, Centre for Ecological Research, MTA, Tihany, Hungary

  • Tibor Erős,

    Affiliation Balaton Limnological Institute, Centre for Ecological Research, MTA, Tihany, Hungary

  • András Specziár,

    Affiliation Balaton Limnological Institute, Centre for Ecological Research, MTA, Tihany, Hungary

  • Péter Sály,

    Affiliation Balaton Limnological Institute, Centre for Ecological Research, MTA, Tihany, Hungary

  • Zoltán Vitál,

    Affiliation Balaton Limnological Institute, Centre for Ecological Research, MTA, Tihany, Hungary

  • Árpád Ferincz,

    Affiliation Department of Aquaculture, Szent István University, Gödöllő, Hungary

  • Tamás Molnár,

    Affiliation Department of Nature Conservation, Kaposvár University, Kaposvár, Hungary

  • Zoltán Szabolcsi,

    Affiliation Institute of Forensic Medicine, Network of Forensic Science Institutes, Budapest, Hungary

  • Péter Bíró,

    Affiliation Balaton Limnological Institute, Centre for Ecological Research, MTA, Tihany, Hungary

  • Eszter Csoma

    Affiliation Department of Medical Microbiology, University of Debrecen, Debrecen, Hungary

Abstract

The European mudminnow (Umbra krameri) is a Middle Danubian endemic fish species, which is characterised by isolated populations living mainly in artificial habitats in the centre of its range, in the Carpathian Basin. For their long term preservation, reliable information is needed about the structure of stocks and the level of isolation. The recent distribution pattern, and the population genetic structure within and among regions were investigated to designate the Evolutionary Significant, Conservation and Management Units (ESUs, CUs, MUs) and to explore the conservation biological value of the shrinking populations. In total, eight microsatellite loci were studied in 404 specimens originating from eight regions. The results revealed a pronounced population structure, where strictly limited gene flow was detected among regions, as well as various strengths of connections within regions. Following the results of hierarchical structure analyses, two ESUs were supposed in the Carpathian Basin, corresponding to the Danube and Tisza catchments. Our results recommend designating the borders of CUs in an 80–90km range and 16 clusters should be set up as MUs for the 33 investigated populations. How these genetic findings can be used to better allocate conservation resources for the long term maintenance of the metapopulation structure of this threathened endemic fish is discussed.

Introduction

Although large floodplain rivers and their associated habitats (e.g., dead arms, backwaters, wetlands) represent a small fraction of the total land area, they nonetheless harbor an outstanding number of endangered and threatened species worldwide [1, 2]. However, large rivers in the developed regions of the world are highly influenced by anthropogenic impacts [3], [4] such as river regulation, construction of flood-control levees, and reservoirs, which disrupt the lateral and longitudinal connections of aquatic systems [5, 6]. Moreover, growing human populations, increased agricultural land use, and the development of infrastructure decreases and fragments the area of natural floodplain ecosystems [7].

Habitat degradation results in the erosion of genetic diversity [8, 9], in particular for those aquatic species which are characterised by special environmental needs, narrow ranges of distribution, and low dispersal ability [10, 11]. Therefore for certain conservation biological projects, such as habitat rehabilitation and reconstruction, the knowledge of population genetic structure is necessary if, for example, resettlement of these species is planned [12, 13]. The task of population genetic investigations in such cases is to assess the rate of gene flow among populations, to identify shrunken populations exposed to the adverse effects of inbreeding, and particularly to determine evolutionary significant units (ESU) within biogeographical patterns [1417]. This information can be used directly to elaborate species conservation and management plans [18, 19].

In Europe, up to 90% of the floodplain area is altered from its natural state [20]. Although less well known, the most severe degradation of floodplain systems has occurred in the Carpathian Basin. Until the middle of the 19th century, more than 21,000 km2 of land was flooded at least periodically in the basin [21]. However, of the ~15,000 km2 floodplain area of the Tisza River (the largest tributary of Danube River, with a drainage area of 157,000 km2), for example, approximately 4,400 km2 was permanently inundated [22] (see Fig 1a). River regulation began in the middle of the 19th century and more than 4,200 kms of levees have since been built, and hundreds of larger meanders were cut along the river. Therefore, the course of the Tisza River was shortened from the original 1419 to 966 km (i.e., by 32%) [20]. The length of the Berettyó River, a tributary of the Tisza River, was also shortened, to half of the original length (from 364 to 174 km) [23]. Concurrent with river regulation, draining of the marshlands was also carried out. The length of draining canals is currently greater than 40,000 km in the central area of the Carpathian Basin [24]. As the result of these large-scale and numerous regulations, the wetland area in the floodplain of the Tisza River was reduced to 539 km2. In the entire Carpathian Basin, the area of wetlands and backwaters directly connected to rivers was reduced to 700 km2 [21]. This broad alteration of the environment has consequently had a serious effect on the fauna. Several vertebrate species preferring these aquatic habitats became rare or extirpated from the basin [25], including some endemic fish species, such as the European mudminnow (Umbra krameri Walbaum, 1792), a Middle Danubian endemic fish species [26], which has suffered a serious population decline. This small (up to 100 mm in total length) stagnophilous species [27] was widely distributed in the marshy habitats of the Danube drainage basin from Vienna to the delta, and in the lower reaches of Dniestr drainage basin [28], prior to river regulation (Fig 1b). The centre of its distribution area is the inner Carpathian Basin [29, 30], where this species was so abundant, that it was used to feed domestic animals (e.g., swine) [31]. At present, only a few sporadic populations have survived the habitat loss and fragmentation [32, 33]. The unexpected invasion of a highly competitive, non-native fish species, the Amur sleeper (Perccottus glenii Dybowsky 1877), has further worsened the situation, especially in the Tiszanian area [34, 35]. Unfortunately, this voracious competitor has recently also appeared in the middle and western parts of the Carpathian Basin [3638]. Due to these impacts, mudminnow is known to have been extirpated from many locations [32], especially from the backwaters connected to the Hungarian Upper and Middle Tisza sections. In just the last decades, more than 30% population decline is estimated to have occurred [28], and European mudminnow has therefore been listed as vulnerable (VU) on the IUCN Red List [39] since 1996. Although the mudminnow is a strictly protected species in Hungary, only few studies have carried out a comprehensive analysis of the distribution of the mudminnow, and are in need of a more recent update [32]. In order to promote the conservation of the remaining stock, captive breeding and rearing techniques have been developed [40, 41], along with some pilot habitat revitalization and population translocation (transferring residual stocks from threatened habitats to newly established nature-like habitats) experiments [42, 43]. Nevertheless, we still do not have a comprehensive overview of the genetic structure of the remaining populations distributed across shrinking and strongly fragmented habitat patches in the Carpathian Basin. According to the findings of field surveys, many of the known populations now consist only of very few individuals, which indicates a high probability of inbreeding and the extreme risk of their extirpation [32, 42, 44]. We do not know, however, how the total genetic diversity of the mudminnow is distributed among populations. Regarding the rapid decrease of mudminnow populations and the limited amount of resources available for their conservation, it is especially important to explore population genetic patterns across the distribution area (i.e., toward the identification of ESUs), and to identify those populations containing most of the genetic diversity in order to set priorities for an effective action plan. For instance, imperative species conservation actions strongly demand information on which populations may be used for supplementing (i.e., serving as a pool of parent fish in artificial breeding programs) endangered populations and stocking reconstructed habitats in different geographic areas and on which populations have the genetic uniquity requiring them to be preserved without mixing them with other populations (i.e., the determination of conservation and management units, MUs and CUs) [17, 45].

thumbnail
Fig 1.

(a) Recent river network of the Carpathian Basin, with 33 sampling sites across eight sampling regions. Sites indicated by different colours belong to different regions. Open circles indicate further known (but not analysed) stocks of mudminnow. Periodically and permanently flooded areas before the beginning of river regulation works (in the mid. 19th century) are indicated by light and dark grey patches respectively. For detailed information, refer to the text and Table 1. (b) Distribution area of European muddminow (red coloured area). The location of the Carpathian Basin in Europe is indicated by a dotted rectangle.

https://doi.org/10.1371/journal.pone.0138640.g001

Therefore, the aims of the current study were: (i) to explore the recent distribution of mudminnow populations in the Hungarian part of the Carpathian Basin, the area inhabited by the majority of the known stocks; (ii) to reveal spatial patterns in the population genetic structure and among sites’ migratory rates; (iii) to assess the conservation biological value of extremely small populations; and (iv) to designate ESUs, CUs and MUs.

Materials and Methods

Ethics Statement

This study was carried out following relevant national and international guidelines pertaining to the care and welfare of fish. Our surveys (e.g., tissue sampling during field surveys) are not qualified as animal tests by the operative Hungarian law (No: 40/2013.(II.14.). However, since the studied species is strictly protected in Hungary, any procedure to be applied in connection with it is controlled by the order of government No: 348/2006.(XII.23.), and is subject to authorisation from the National Inspectorate for Environment, Nature and Water. The collection and storage of mudminnow tissue samples was authorized by this bureau (permission numbers: 14/881/5/2011, 14/678-9/2012). During the sampling procedure, all efforts were made to minimize suffering. Field studies did not involve other fish species that were endangered according to the IUCN Red List of Threatened Species v. 2014.3 (www.iucnredlist.org).

Sample collection, DNA extraction and microsatellite amplification

Studied areas and sites were selected based on the last published distribution data of mudminnow [32]. Specimens were collected by electrofishing between 2011 and 2013, from each region of the Carpathian Basin where the presence of this species had previously been noted (Table 1). To minimize suffering, individuals collected for this study (N = 5–20 per site, depending on the stock size) were narcotized using clove oil. The active agent of clove oil is eugenol (4-allyl-2-ethoxyphenol), which is considered non-carcinogenic, non-mutagenic, and a "Generally Recognized As Safe" (GRAS) substance by the U.S. Food and Drug Administration [46, 47]. The aqueous emulsion of clove oil was deposited into the dangerous refuse containers of the Balaton Limnological Institute after use. For the genetic survey, an approximately 2mm2 anal fin clip was sampled from each specimen and stored in 96% ethanol at -20°C until DNA extraction. After fin tissue sampling, fish were released to the habitat they were captured from. DNA was isolated with DNeasy Blood and Tissue kit (Qiagen, Germany), using 10–20 mg of tissue, as per the manufacturer’s instructions. Quality and quantity of the extracted DNA were verified using a NanoDrop 2000c Spectrophotometer (Thermo Scientific, USA).

thumbnail
Table 1. Name, location, code, and geographical position of the sampling sites.

https://doi.org/10.1371/journal.pone.0138640.t001

Nine microsatellite markers previously developed by Winkler and Weiss [48] were investigated by multiplex PCR. 50 ng of DNA was amplified in a multiplex PCR reaction with three primer pairs in a PCR reaction using Type-it Microsatellite PCR Kit (Qiagen, USA), and fluorescently labelled forward primers according to the protocol of the kit (details and multiplex combinations are provided in Table 2). The final concentration of each primer was 0.2 μM, and the annealing temperature was 60°C in each PCR. Amplified products were detected using an ABI 3130 sequencer (Applied Biosystems, USA). Peak Scanner Software v1.0 was used to analyse the data. The functionality of all primer pairs and the reliability of multiplex PCR were tested in singleplex reactions using the same PCR conditions and reagents.

thumbnail
Table 2. Primer sequences, fluorescent dyes and combinations of primers in three (1., 2., 3.) multiplex PCR, number of alleles, locus-specific Fst values, number of populations out of the 33 in which deviations from Hardy-Weinberg Equilibrium (HWE) were detected.

(*UkrTet2 UkrTet3 sequences are identical, therefore locus UkrTet2 was omitted from the further analyses, for more details refer to text), and Mean ± SD estimated frequency of null alleles.

https://doi.org/10.1371/journal.pone.0138640.t002

Statistical analyses

Genetic variation.

Microsatellite markers (Table 2) and primers published by Winkler and Weiss [48] were checked in GenBank before PCR and microsatellite analysis. Based on the available sequences, it was suggested that Umbra krameri microsatellites UkrTet2 (GenBank: FJ228219.1) and UkrTet3 (GenBank: FJ228220.1) sequences are identical, the reverse complements of each other. The DNA strand direction was mistakenly reversed for these. Data analysis of 156 fish samples (the first 14 populations) further strengthened this hypothesis. Allele numbers (n = 18) and allele frequencies, allele specific Fst values (Fst = 0.156) and the proportion of null alleles (prop. = 0.040372339) were equal for the UkrTet2 and UkrTet3 loci. Consequently, the UkrTet2 marker was left out of further statistical analyses.

Fisher’s exact test of linkage disequilibrium and tests for deviations from Hardy-Weinberg equilibrium were conducted using GENEPOP 4.2.2 [49] for each locus using a Markov chain of 10,000 dememorization steps, 20 batches, and with 5,000 iterations per batch. MICROCHECKER 2.2.3. [50] was used to estimate the frequency of null alleles for each locus. Number of alleles, and locus specific Fst values, were calculated using GenAlEx 6.5 [51]. Similarly, this software was used to estimate the mean of allelic richness (Ar), Shannon's Information Index (I), observed and expected heterozygosity (Ho, He), fixation index (F) and list the private alleles for each population.

Population structure.

The Lynch & Ritland [52] estimator was used to calculate between-specimen pairwise relatedness. From this semimatrix, mean within-population pairwise values were calculated with 999 permutations, and 1,000 bootstraps in GenAlEx 6.5 [51]. Because in some cases null alleles were detected, their effect on pairwise Fst values were corrected using the ENA procedure of Chapuis and Estoup [53, 54]. Genetic distances among populations were calculated using the Cavalli-Sforza and Edwards [55] estimator after INA correction [53], and presented in PCoA ordination. GenAlEx 6.5 [51] was used to perform the analysis of molecular variance (AMOVA), with 999 permutations for population differentiation and hierarchical partitioning of genetic variation among and within regions and populations (F-statistics: FrtH0: individuals are shuffled among regions, FsrH0: individuals are shuffled within regions, FstH0: populations are shuffled among regions, Fis–H0: individuals are shuffled within populations, and Fit–H0: individuals are shuffled in the whole sample).

Genetic population structure was inferred using the hierarchical approach [56] of the STRUCTURE analysis [57] to estimate the most probable number of genetic groups (clusters, K) for all analysed individuals. Namely, the STRUCTURE analyses was first run including all samples, then samples were separated by river drainage basin, and finally by region. (A, B, C, etc.) Values of K were investigated from 1 to 20, with a burn-in period of 100,000 followed by 100,000 MCMC iterations and 10 runs for each K using an admixture model with correlated allele frequencies. Results of these Bayesian statistics were evaluated by STRUCTURE HARVESTER [58], implementing the (deltaK) Evanno method [59]. Results of the 10 repetitions were combined using the software CLUMPP 1.1.2. [60]. For the genetic assignment of the studied individuals, Bayesian cross validation tests [61] were carried out on drainage basin, regional, and population levels using GeneClass2 [62] software. Cross validation tests were carried out similarly on the clusters defined at various levels by the hierarchical STRUCTURE analyses.

Migration and spatial analyses.

Migration rate estimation was performed by MIGRATE-N 3.2.15 [63, 64], using the maximum likelihood method under the following parameterization: 20 short chains (500 trees used out of the sampled 10,000) and 5 long chains (5,000 trees used out of the sampled 100,000). Missing data were not included. Theta values were generated from the Fst-calculation. Mutation rates among loci were estimated from the data. Migration was estimated for the original regional groups designated by catchment basin, and based on the results of hierarchical STRUCTURE analyses at the three levels.

Spatial genetic structure was assessed at the population level using Mantel tests [65], through the comparison of pairwise Fst data and pairwise straight line geographic distances (GGD) with 999 randomisations for the whole dataset and at drainage basin levels. On the individual level, two kinds of spatial autocorrelation [66, 67] computation were made. Spatial autocorrelation computation is suitable to accurately identify the scale at which population genetic structure is detectable. In this case, “r” coefficients were calculated using multiple distance class analyses. This method plots “r” as a function of increasing distance class sizes. The last distance class for which “r” is significant is considered to be the limit of detectable isolation by distance (IBD). In the second case, autocorrelation coefficients, “R”, were plotted as a function of discrete distance classes, partitioned so as to achieve a similar number of pairwise comparisons for each class [68]. Positive and significant “R” values indicate IBD and the “x” intercept provides an estimate of the extent of IBD. For more details see Peakall et al. [69]. All of the spatial autocorrelation computations, including the GGD semimatrix calculation from geocoordinates, were made in GenAlEx v6.5 [51].

Results and Discussion

Results

Distribution patterns.

As a result of our broad field investigations, European mudminnow was observed at more than 40 sampling sites from the Middle Danubian catchment (Fig 1). A total of 404 mudminnow specimens were sampled for population genetic research, from 33 sites, originating mainly from artificial habitats (i.e., ditches, canals, and ponds), across eight regions (Fig 1, Table 1). In several locations, no more than 5 individuals were able to be collected at one time, despite the intensive sampling effort.

Data quality and Genetic variation.

Results of Fisher’s exact tests did not reveal evidence for large allele drop-out for any locus, or for linkage disequilibrium between any pairs of loci. All microsatellite loci showed significant deviation from the Hardy-Weinberg equilibrium state if all 404 specimens were analysed as one entire population (404 samples), and we calculated high Fis values (>0.08) at all loci among the 404 samples. These results suggest strong genetic substructure caused by geographical isolation and/or non-random mating in subgroups. Characteristics of each microsatellite locus used in our work are provided in Table 2. The total number of alleles per locus varied between 6 and 25 (mean: 17.6). Locus specific Fst values ranged between 0.148 and 0.482. Chi-Square Tests for Hardy-Weinberg Equilibrium for each loci at each population showed significant values in 17 cases (see Table 2 and S1 Table).

A total of 136 alleles were observed for the eight loci used in the analysis. The whole raw dataset used for further analyses is indicated in S2 Table in GENEPOP format. Average allelic richness, such as observed heterozygosity, showed high-level differences among populations. These ranged from 2.75 to 10 and from 0.44 to 0.85 respectively, within the studied populations. Shannon's Information Index showed a high level of variability as well, ranging between 0.8 and 1.97 (Table 1). Deviation from HWE was significant (p<0.05) in two studied populations (E8, G2).

As MICROCHECKER analyses proved the presence of false homozygotes (null alleles) in some cases (i.e., for UkrTet1 at sites B3, D2, E3, F1, for UkrTet3 at site E3, for UkrTet4 at site B2), therefore null allele correction was made prior to further analyses (see S3 Table).

13 of the detected 136 alleles were unique to a single location (Table 3). The frequency of the unique alleles ranged between 0.031 and 0.071. Eight of the unique alleles were observed in sites of the Tisza River drainage basin, and five from the Danubian part of the study area. No unique alleles were found for small stocks (where N = 5), in populations A3, C1, E5 and F2.

thumbnail
Table 3. Summary of private alleles (allele sizes and frequencies) by population.

Codes of populations correspond to those of Table 1 and Fig 1.

https://doi.org/10.1371/journal.pone.0138640.t003

Pairwise genetic distances ranged between 0.198 and 0.821 (mean: 0.621). Pairwise Fst values ranged between -0.01 and 0.476 (mean: 0.195), and only five of the 528 pairwise comparisons (A3-4, B1-2, C1-2, F1-2 and G4-5) did not show significant differences (Table 4).

thumbnail
Table 4. Pairwise Fst values generated by FREENA for each pair of sites using ENA method (below the diagonal).

Italicized numbers indicate significant differentiations (p<0.05). Cavalli-Sforza and Edwards genetic distances for each pair of populations are indicated above the diagonal. Codes of sites correspond to those of Table 1.

https://doi.org/10.1371/journal.pone.0138640.t004

Mean within-population pairwise correlation relatedness values, “r”, ranged between 0.063 and 0.651, and were significant in all cases (Fig 2). These values correlated negatively with the observed heterozygosity and with Shannon's Information Index (Spearman’s rho: -0.7681, p<0.01, and -0.956, p<0.001 respectively). Therefore, higher “r” values refer to higher separation and lower genetic diversity (these features may indicate inbreeding). In some cases (e.g., populations G2 = 0.651, C2 = 0.498, D2 = 0.467), high values were found, but pronounced regional patterns were not. In almost all sampling regions (except regions A and B), one or more populations formed by strongly related individuals were observed.

thumbnail
Fig 2. Correlogram representing the mean within-population pairwise similarities using the Lynch & Ritland (1999) estimator.

Upper (U) and lower (L) confidence limits bound the 95% CI of the null hypothesis of 'No Difference' across the populations using 999 random permutations, and for estimates of r by bootstrapping 1,000 pairwise comparisons. Whiskers represent the highest and lowest values within a dataset for each distance class. Regions are separated by vertical grey lines. Codes of sites correspond to those of Table 1.

https://doi.org/10.1371/journal.pone.0138640.g002

The PCoA plot based on pairwise population genetic distances (Fig 3) showed clear separation along the “x” axis of eastern (Tisza drainage basin) and western (Danube drainage basin) sites. Moreover, the populations in the Tisza drainage basin showed higher similarity. In the case of these populations, the separation along the second axis is also important. Pronounced separation of neighbouring populations (i.e., populations from the same region) within a region can only be seen in the case of population E1. This population belongs to the Danube drainage basin Middle Hungarian region, but shows stronger similarity to the Mura region or to the C1, C2, and D4 populations of the Tisza drainage basin.

thumbnail
Fig 3. PCoA plot of pairwise population genetic distances (using Cavalli-Sforza and Edwards estimator after FREENA-INA null allele correction).

Percentage variances represented by axes are shown in parentheses. Codes and colours of sites correspond to those of Table 1 and Fig 1.

https://doi.org/10.1371/journal.pone.0138640.g003

AMOVA analysis showed that among-population and among-regional differences represented 12.8% and 7.8% of the total variation respectively, while within-population differences accounted for 78.9% of the total genetic variance. Results of the regional F-statistics, except Fis, showed significant within- and among-regional differentiation and population level differences.

Bayesian hierarchical clustering separated Danubian and Tiszaian mudminnow stocks in the first step (Fig 4). The only exception is the E1 population, which was assigned clearly to the Tiszanian clade. The Tiszanian group was split into 4 clusters; the E1 population, Tápió region (C1-2), and sites from NE Bihar region (D1-3) were separated from the others. Stocks form the Upper Tisza region (A1-4), Borsodi-mezőség region (B1-3), and SW of Bihar plain (D4-5) were separated only in the third clustering step, and the extent of separation is not as unambiguous as in the previous steps. The Danubian branch split to two larger groups; the Middle Hungarian region together with the Hanság-Szigetköz region (E+F) formed a group separated from the Balaton and Mura (G+H) regions. The E+F cluster split further into four minor groups (E2+E8; E3+E4+E5; E6+E7; F1-F4). Populations in the G+H cluster formed five minor groups. Distinct populations were found in this cluster, except G3-G4-G5 stocks, which all originated from the Kis-Balaton marshland area and formed one population.

thumbnail
Fig 4. Estimated population structure as inferred by three rounds of hierarchical structure analyses.

Each individual is represented by a thin horizontal line, which is partitioned into K-coloured segments representing individual’s estimated membership fractions in K clusters. Black lines separate individuals from different sampling sites. The most probable K for a sample, given by the arrows, is based on the results of StructureHarvester using the Evanno method. Codes of sampling sites correspond to those of Table 1. Red codes represented on the three hierarchical levels correspond to the codes used in Fig 6.

https://doi.org/10.1371/journal.pone.0138640.g004

Bayesian cross-validation showed that 98.3%, 89.6%, and 71.8% of the individuals were grouped correctly on drainage basin, region and population levels respectively (Fig 5). Most of the reclassifications were made within the same region, mainly between neighbouring sites (e.g., C1-C2).

thumbnail
Fig 5. Plots of the Bayesian cross-validation test for microsatellite data.

Correctly classified individuals are placed on the diagonal. The square colour corresponds to the proportion of individuals of posterior group assignment based on posterior probabilities. Rows correspond to actual sites (a priori), while columns correspond to inferred sites (posteriori). White thick lines separate the Tisza and Danube drainage basins, red rectangles show the regional detachments of populations. Codes of sites correspond to Table 1.

https://doi.org/10.1371/journal.pone.0138640.g005

Result of the cross validation procedure showed negligible linkage between the Danubian and Tisza drainage basins. Only three specimens (i.e., 0.7% of the total sample) were misclassified between the two drainage basins. Within a drainage basin, the percentage of misclassified cases is higher; five specimens of the total 156 (i.e., 3.2%), and 13 specimens of the 248 (i.e., 5.2%) were classified to other regions on the Tisza and Danube drainage basins respectively. The highest number of misclassified cases was found between the Middle-Hungarian and Hanság-Szigetköz regions (E and F regions) where 4.9% and 10% of the stocks were misclassified.

Population assignment results from GeneClass2 strongly support the structure revealed on the three hierarchical levels identified by Bayesian STRUCTURE analyses. 97.8, 96.0 and 94.6% of the individuals were classified correctly on the three hierarchical levels respectively (see: Fig 4). The third level of the hierarchical STRUCTURE analyses was therefore found to have classified the individuals much more reliably than the original population based classification.

The between-region (A-H) migration ranged between 0.28 and 2.31 individuals per generation. Mutation-scaled immigration rate is less than one individual per generation in most cases, therefore negligible gene flow is assumed among the eight studied regions in the Carpathian Basin, especially in the case of regions C and H.

Migration computations between the two clusters (corresponding to the Danube and Tisza catchments, designated on the first level of hierarchical STRUCTURE analysis) showed that the Danube system received 8.19 (Θ = 2.4941 and M = 3.286) and released 7.24 (Θ = 2.5105 and M = 2.887) migrating individuals per generation to the Tisza drainage basin alone.

On the second level of STRUCTURE analysis, the immigration rate (M) between the four clusters (Figs 4 and 6) found in the Tisza catchment ranged between 0.75 and 1.52. The highest level of mutation-scaled immigration rate was found between the “a” and “c” clusters. In the case of the two Danubian clusters (see e and f clusters on Fig 4) the migration rates were less than one in both directions, resulting in a 1.36 and 1.48 individual per generation migration rate. These groups are therefore practically separated. Pairwise between-region and between-clusters immigration rates are presented in Tables 59.

thumbnail
Fig 6. Migration within and among different clusters designated by hierarchical STRUCTURE analyses.

Sites in the green areas belong to the same evolutionarily significant units (ESUs) identified by the 1st level of hierarchical STRUCTURE analysis. Sites framed by red and blue lines belong to the same conservation units (CUs) and management units (MUs) designated by the 2nd and 3rd level by STRUCTURE analyses respectively. Arrows directions show between cluster migrations, where the migration rate is >1 individual per generation. Arrow colouring corresponds to hieararchical levels. Open circles show sampling sites. Codes of clusters correspond to Fig 4.

https://doi.org/10.1371/journal.pone.0138640.g006

thumbnail
Table 5. Migration computation among the eight hydrological regions preliminary defined, + receiving population (for region codes see Table 1).

https://doi.org/10.1371/journal.pone.0138640.t005

thumbnail
Table 6. Theta values (Θ) and immigration rates (M) of the six clusters defined in the second round of hierarchical STRUCTURE analysis (For cluster codes see Fig 4, + receiving population).

https://doi.org/10.1371/journal.pone.0138640.t006

thumbnail
Table 7. Within “a” cluster theta values (Θ) and immigration rates (M). Clusters defined in the third round of hierarchical STRUCTURE analysis.

+ receiving population (For cluster codes see Fig 4).

https://doi.org/10.1371/journal.pone.0138640.t007

thumbnail
Table 8. Within “e” cluster theta values (Θ) and immigration rates (M). Clusters defined in the third round of hierarchical STRUCTURE analysis.

+ receiving population (For cluster codes see Fig 4).

https://doi.org/10.1371/journal.pone.0138640.t008

thumbnail
Table 9. Within “f” cluster theta values (Θ) and immigration rates (M). Clusters defined in the third round of hierarchical STRUCTURE analysis.

+ receiving population (For cluster codes see Fig 4).

https://doi.org/10.1371/journal.pone.0138640.t009

A Mantel test conducted on the entire dataset at the population level showed no significant correlation between pairwise Fst and GGD semimatrices (Rxy = 0.061, p>0.05). The inverse was found for the two drainage basins, where the genetic differences showed significant correlations both for the Tisza (Rxy = 0.4196, p<0.01) and the Danubian drainage basins (Rxy = 0.2489, p<0.01) as well. Therefore, population level IBD can be detected only within the Danubian and Tisza drainage basins (see Fig 7)

thumbnail
Fig 7. Results of Mantel tests conducted on the entire dataset (A), on the Tisza (B), and on the Danube drainage basin (C).

https://doi.org/10.1371/journal.pone.0138640.g007

Spatial autocorrelation computation carried out on the entire dataset at the individual level showed that IBD was not appreciable beyond an approximate 340–360 km distance. There was a significant positive correlation for the first two distance classes, which intercepted the “x” axes at 124 km (Fig 8A and 8B) and was therefore significant within this range. Autocorrelation calculation at the drainage basin level showed similar results in both cases. The IBD is still detectable to approximately 150–160 km, but it is significant only within an 80.1 and 86.7 km range for the Tisza and Danube drainage basins respectively (see Fig 9).

thumbnail
Fig 8. Correlograms showing autocorrelations calculated for the entire dataset.

Coefficient “r” as a function of increasing cumulated distances (A). Correlogram “B” showing “R” as a function of distance classes designated by a similar number of pairwise comparisons for each class. Intercept value for the “x” axis is 124 km. Red dotted lines show the upper and lower confidence limits bounding the 95% CI of the null hypothesis of 'No Difference' across the populations using 999 random permutations, and for estimates of correlation coefficients by bootstrapping 1,000 pairwise comparisons for each distance class. Whiskers represent the highest and lowest values within a dataset.

https://doi.org/10.1371/journal.pone.0138640.g008

thumbnail
Fig 9. Correlograms showing the autocorrelation coefficients as a function of increasing distances (“r”) and distance classes (“R”) designated by a similar number of pairwise comparisons for each class, for the Tisza (A, C) and Danube drainage basin (B, D) respectively.

Red dotted lines show the upper and lower confidence limits bounding the 95% CI of the null hypothesis of 'No Difference' across the populations using 999 random permutations, and for estimates of r by bootstrapping 1,000 pairwise comparisons for each distance class. Whiskers represent the highest and lowest values within a dataset. Intercept values for the “x” axes are 80.1 and 86.7kms for the Tisza and Danube drainage basins respectively.

https://doi.org/10.1371/journal.pone.0138640.g009

Discussion

Tockner et al. [4] noted that fragmentation, water stress, land use change, and the introduction of non-native species are the most important threats to freshwater biodiversity. All these factors currently threaten mudminnow stocks, and their impacts are reflected in present day distribution and population genetic patterns. In this study, population genetic structure and dynamics of genetic mixing among populations in a strongly altered and fragmented landscape were analysed, and revealed some basic principles of optimal species conservation strategies.

Distribution patterns

The field surveys of the current work have provided new information about the present day distribution of the European mudminnow at the centre of its range. In the “C” and “F” regions, stable and large assemblages were found, contrary to earlier communication which noted that the mudminnow was extremely rare and had been extirpated from most habitats in these regions [32]. At the same time, due to the expansion of Amur sleeper (Perccottus glenii) [70], European mudminnow has disappeared from many locations in the Hungarian Upper Tisza region (i.e., region “A” of this study), where it had previously been abundant [44]. At present, these shrinking populations seem to be the most threatened ones within the Carpathian Basin.

In almost all regions, the largest assemblages were found in artificial habitats. Of the 33 sampling sites, only seven were in a near-pristine state (Table 1). Consequently, similar to other threatened stagnophilous fish, it seems that some of the human-altered or man-made habitats (e.g., irrigation and drainage canals) with rich macrophyte coverage and permanent water supply might be the last chance to provide refuge for the shrinking European mudminnow populations [27, 71].

Genetic and spatial structure

This study presents the first data on the population genetic structure of threatened and endemic European mudminnow for the centre of its distribution range. In general, remarkable genetic differentiation was found among the 33 populations originating from eight regions. The spatial structure revealed seems to be quite fixed with negligible migration among the studied regions (Table 5). The genetic structure corresponded clearly to past and recent geographical migratory routes (i.e., it reflected natural distribution patterns). The only identified anthropogenic alteration within this pattern is a result of a documented stock transfer from Middle Hungary to the Hanság-Szigetköz region (i.e., from region E to region F) in April 2005 [42]. Results of Bayesian cross validation tests (Fig 5), confirmed by migration tests (Fig 6), suggest that 10% and 5% of the specimens occurring here are closely related to Middle Hungarian (E) and Hanság-Szigetköz (F) populations respectively (Figs 5 and 6; Tables 59.). An alternative explanation for the mixed genetic composition of the Hanság-Szigetköz (F) region population could be the dispersion along the Danube valley route, as has been proven for many other species [72, 73]. In any case, SRUCTURE analysis showed (Fig 4) that the aforementioned stocking activity had only a moderate effect on the genetic pattern of mudminnow assemblages inhabiting this region.

Some of the studied populations showed complete isolation and a high level of relatedness (Figs 2 and 5). Populations which are characterised by special environmental requirements and restricted expansion ability may hold unique pieces of genetic diversity, and their loss can therefore cause a significant reduction in genetic diversity in general [74]. Therefore, to conserve the genetic diversity of this species, special attention must be paid to the preservation of these separated stocks (e.g., B2, D1, and G2 on Fig 2).

Although European mudminnow has a restricted dispersal ability [28], and historical processes therefore mainly underlie the pronounced spatial pattern in the genetic structure, our analyses also captured some human-induced effects related to the anthropogenic modification of the Carpathian Basin floodplains. Historical hydrological data [75, 76] suggest that prior to river regulation there were no notable barriers among some sub-drainage basins in the floodplain of the Tisza River (Fig 1). Periodic floods inundated almost the entire floodplain, sometimes for more than 3 months at a time (see: [22]), and migration could therefore take place freely among A-B-D regions in the Tisza drainage basin. As a result of river regulation, which began in 1846 and ended in the 1930s [77], the connections between regions became much more restricted. The strong distinctiveness of populations inhabiting the Tápió region (i.e., the “C” region here), outside the main floodplain, also supports this hypothesis (Fig 3).

According to data from the literature (cited in: [28, 39]), mudminnow may spawn at the age of one year for the first time. The separation of some Tiszanian regions could therefore have begun more than 100 generations ago. The separation of the Danubian and Tiszanian systems, as well as the mudminnow stocks, is presumably much older. The within-region separations are not complete in many cases, however, because the continuous flood-protection and regulation activity (i.e., dredging, flow direction shifts within the drain-canal systems) could have caused the connections between stocks to be re-established from time to time [78]. E1 is the only population among the 33 studied which seems to be misclassified by both the STRUCTURE and the pairwise population genetic distance analyses (Figs 3 and 4). Although this site currently belongs to the Danubian drainage system, the E1 population seems to be more related to those of the Tisza catchment. This relationship can be explained by paleohydrological changes caused by the emergence of Gödöllő Hilly Region in the Quaternary period [79], when the river system of this area was partly disjointed and the flow direction of the upper sections of the streams were diverted to the NW (i.e., to the Danube). At the same time, the two catchment areas were separated only by marshy areas referred to as “intra-valley drainage basin divides” [80], which these assemblages can traverse. This assumption is supported by the fact that the closest populations (C1-C2) to the Tápió system showed the highest similarities and bidirectional gene flow was found between the two areas. However, it should be noted that both assemblages showed high levels of isolation and inbreeding, and the connection between them is therefore presumably strictly limited. Similarly, contact between the G and H regions is plausible, because the isolation of the Middle and South Danubian catchments started only at the end of Pleistocene (by the formulation of the Lake Balaton basin) and this separation is still not complete [81]. Our results correspond with the findings of Brauer et al. [11], who investigated a South Australian endemic fish species with very similar environmental requirements to those of European mudminnow, the pygmy perch (Nannoperca obscura Klunzinger, 1872), and found similarly high levels of isolation, and designed CUs in very similar ranges. These results therefore show that preferred habitat type and dispersion abilities may have major roles in the formulation of these features rather than the geographic range (i.e., Australia vs. Europe) or the taxonomic position (i.e., Perciformes vs. Esociformes) of the studied species.

Implications for conservation and management

Although pairwise Fst computation showed significant differentiation in most cases, none of the small stocks differed considerably from their neighbouring populations (Table 5; Fig 5). Moreover, these small stocks did not carry any private alleles (Table 3). It is therefore likely that these small stocks belonged to a larger metapopulation system, and could have been isolated only recently. Therefore, special consideration for their preservation is necessary only if there is no larger population nearby presenting the same genetic features.

Based on the current data, it is recommended that any reconstructed or newly established habitats should be stocked by individuals originating from the same CU, namely from populations within an 80–90 km range. The 16 clusters separated in the third step of the hierarchical STRUCTURE analysis can be accepted as MUs for conservation, instead of populations. In many cases, the MUs agree with some strongly isolated populations (e.g., populations C1, E1, G1, and G2). These separate entities need special attention to prevent substantial genetic diversity loss. The ex situ preservation of these populations could also be considered, preferably by artificially recruiting and translocating them to fish-free revitalized or newly established nature-like habitats (see e.g., [43]).

Conclusions

The investigations reported here reveal unexpectedly high genetic diversity of this endemic fish species, despite its declining and fragmented populations. The studied stocks show high levels of genetic variability, and this pattern seems to have been only slightly influenced by anthropogenic impacts (i.e., resettlements) so far. Continued habitat loss, along with the invasion of non-native competitors, strongly threaten this habitat specialist, endemic fish. Climate change increases the probability and severity of dry periods in the area, and represents a high risk for European mudminnow with its shrunken and fragmented population structure [82]. Consequently, there is a strong need for implementing comprehensive conservation management programs. Nevertheless, the extremely small and/or isolated populations (e.g., in region “C”) would require special attention, and thus should be prioritized for conservation.

Supporting Information

S1 Table. Summary of Chi-Square Tests for Hardy-Weinberg Equilibrium.

https://doi.org/10.1371/journal.pone.0138640.s001

(DOCX)

S2 Table. Microsatellite data of Carpathian Umbra stocks (404 individuals, 33 populations, 8 regions).

https://doi.org/10.1371/journal.pone.0138640.s002

(DOCX)

S3 Table. Results of null allele analyses using MICROCHECKER software.

https://doi.org/10.1371/journal.pone.0138640.s003

(DOCX)

Acknowledgments

We would like to thank our colleagues, Katalin Keresztessy, András Ambrus, László Antal, and Zoltán Sallai for their help in sample collections and to Stephanie CJ Palmer who corrected the grammar of the manuscript.

Author Contributions

Conceived and designed the experiments: PT TE PB. Performed the experiments: PT EC TE AS PS ZV ÁF. Analyzed the data: PT ZS TM EC. Contributed reagents/materials/analysis tools: PB EC. Wrote the paper: PT EC TE AS PS ZV ÁF TM ZS BP.

References

  1. 1. Williams JD, Dodd CK Jr (1980) Importance of wetlands to endangered and threatened species. In: Greeson PE, Clark JR., Clark JE (eds.) Wetland functions and values: the state of our understanding. American Water Resources Association, Minneapolis, pp. 565–575.
  2. 2. Buijse AD, Coops H, Staras M, Jans LH, Van Geest GJ, Grift RE, et al. (2002) Restoration strategies for river floodplains along large lowland rivers in Europe. Freshwater Biol 47: 889–907.
  3. 3. Meybeck M (2003) Global analysis of river systems: from Earth system controls to Anthropocene syndromes Phil. Trans. R. Soc. 358(1440): 1935–1955.
  4. 4. Tockner K, Pusch M, Borchardt D, Lorang MS (2010) Multiple stressors in coupled river–floodplain ecosystems. Freshwater Biol 55: 135–151.
  5. 5. Dynesius M, Nilsson C (1994) Fragmentation and flow regulation of river systems in the northern 3rd of the world. Science 266: 753–762. pmid:17730396
  6. 6. Kondolf MG, Boulton AJ, O'Daniel S, Poole GC, Rahel FJ, Stanley EH, et al. (2006) Process-Based Ecological River Restoration: Visualizing Three-Dimensional Connectivity and Dynamic Vectors to Recover Lost Linkages. Ecol Soc. 11(2): 5.
  7. 7. Nilsson C, Berggren K (2000) Alterations of riparian ecosystems caused by river regulation. BioScience 50: 783–792.
  8. 8. Allendorf FW, Luikart G (2009) Conservation and the genetics of populations. John Wiley & Sons.
  9. 9. Hudman SP, Gido KB (2013) Multi-scale effects of impoundments on genetic structure of creek chub (Semotilus atromaculatus) in the Kansas River basin. Freshwater Biol 58: 441–453.
  10. 10. Lelek A (1987) The freshwater fishes of Europe. Threatened fishes of Europe, Volume 9. AULA-Verlag, Wiesbaden, 343 pp.
  11. 11. Brauer CJ, Unmack PJ, Hammer MP, Adams M, Beheregaray LB (2013) Catchment-Scale Conservation Units Identified for the Threatened Yarra Pygmy Perch (Nannoperca obscura) in Highly Modified River Systems. PLoS ONE 8(12): e82953. pmid:24349405
  12. 12. Sarrazin F, Barbault R (1996) Reintroduction: challenges and lessons for basic ecology. Trends Ecol Evol 11: 474–478. pmid:21237930
  13. 13. Sanderson EW (2006) How many animals do we want to save? The many ways of setting population target levels for conservation. BioScience 56(11): 911–922.
  14. 14. Moritz C (1994) Defining ‘Evolutionarily Significant Units’ for conservation. Trends Ecol Evol 9: 373–375. pmid:21236896
  15. 15. Crandall KA, Bininda-Emonds ORP, Mace GM, Wayne RK (2000) Considering evolutionary processes in conservation biology. Trends Ecol Evol 15: 290–295. pmid:10856956
  16. 16. Fraser DF, Bernachez L (2001) Adaptive evolutionary conservation: towards a unified concept for defining conservation units. Mol Ecol 10: 2741–2752. pmid:11903888
  17. 17. Funk WC, McKay JK, Hohenlohe PA, Allendorf FW (2012) Harnessing genomics for delineating conservation units. Trends Ecol Evol 27(9): 489–496. pmid:22727017
  18. 18. Alves MJ, Coelho H, Collares-Pereira MJ, Coelho MM (2001) Mitochondrial DNA variation in the highly endangered cyprinid fish Anaecyprishispanica: importance for conservation. Heredity 87: 463–473. pmid:11737295
  19. 19. O’Brien SJ (1994) A role for molecular genetics in biological conservation. Proc P Natl Acad Sci USA 91: 5748–5755.
  20. 20. Tockner K, Uehlinger U, Robinson CT (2009) Rivers of Europe. Academic Press.
  21. 21. Nagy L (2013) Az elöntött öblözetek a mai Magyarország területén. In: Lajos (szerk.) Szlávik XXXI. MHT Vándorgyűlés Kiadványa. Gödöllő, Magyarország, 2013.07.03–2013.07.05. Budapest: Magyar Hidrológiai Társaság, Paper 17. ISBN: 963-8172-31-0 [in Hungarian]
  22. 22. Botár I, Károlyi Zs (1971) A Tisza szabályozása I-II. (Regulation works on the River Tisza I-II) Series of Hydraulic Engineering History in Hungary 3., Budapest pp. 65. [in Hungarian]
  23. 23. Dankó I (1979) A Berettyó-völgy középkori települési, közlekedési és árucsereviszonyai—A Hajdú-Bihar Megyei Levéltár évkönyve, 6: 5–18. [in Hungarian]
  24. 24. Martonné Erdős K (2004) Magyarország természeti földrajza I, Debreceni Egyetem Kossuth Egyetemi Kiadója, Debrecen, pp. 245. [in Hungarian]
  25. 25. Haraszty L (ed.) (2014) Natura 2000 fajok és élőhelyek Magyarországon. (Natura 2000 species and habitats in Hungary), Pro Vértes közalapítvány, Csákvár. pp. 955 [in Hungarian]
  26. 26. Gaudant J (2012) An attempt at the palaeontological history of the European mudminnows (Pisces, Teleostei, Umbridae). Neues Jahrb Geol P-A. 263/2, 2012, S. 93–109.
  27. 27. Pekárik L, Hajdú J, Košco J (2014) Identifying the key habitat characteristics of threatened European mudminnow (Umbra krameri Walbaum 1792) Fund Appl Limnol 184(2): 151–159.
  28. 28. Wanzenböck J (1995) Current knowledge on the European mud-minnow Umbra krameri Walbaum, 1792. (Proceedings of the first international workshop on Umbra krameri) Ann Naturhist Mus Wien 97B: 439–449
  29. 29. Kottelat M, Freyhof J (2007) Handbook of European freshwater fishes. Publications Kottelat, Cornol, Switzerland, 646 p.
  30. 30. Kuehne LM, Olden JD (2014) Ecology and Conservation of Mudminnow Species Worldwide. Fisheries 39/8: 341–351.
  31. 31. Herman O (1887) A magyar halászat könyve. Királyi Magyar Természettudományi Társulat, Budapest, pp.847 [in Hungarian]
  32. 32. Sallai Z (2005) A lápi póc (Umbra krameri) magyarországi elterjedése, élőhelyi körülményeinek és növekedési ütemének vizsgálata a kiskunsági Kolon-tóban. A Puszta 1(22): 113–172 [in Hungarian]
  33. 33. IUCN (2003). 2003 IUCN Red List of Threatened Species. IUCN, Gland, Switzerland. Available at: http://www.redlist.org.
  34. 34. Reshetnikov AN (2010) The current range of Amur sleeper Perccottus glenii Dybowski, 1877 (Odontobutidae, Pisces) in Eurasia. Russian Journal of Biological Invasions 1(2): 119–126.
  35. 35. Reshetnikov AN (2013) Spatio-temporal dynamics of the expansion of rotan Perccottus glenii from West-Ukrainian centre of distribution and consequences for European freshwater ecosystems. Aquat Invasions 8(2): 193–206.
  36. 36. Erős T, Takács P, Sály P, Specziár A, György ÁI, Bíró P (2008) Az amurgéb (Perccottus glenii Dybowski, 1877) megjelenése a Balaton vízgyűjtőjén. Halászat 101/2: 75–77. [in Hungarian with English abstract]
  37. 37. Takács P, Vitál Z (2012) New amur sleeper (Perccottus glenii Dybowski, 1877) locality beside the Hungarian lower Dabube section Halászat 105 (4): 6. [in Hungarian with English abstract]
  38. 38. Takács P, Czeglédi I, Ferincz Á (2015) First data of amur sleeper (Perccottus glenii Dybowski, 1877) from the Hungarian Drava system. Halászat 108: in press. [in Hungarian with English abstract]
  39. 39. Freyhof J (2013) Umbra krameri. The IUCN Red List of Threatened Species. Version 2014.3. <www.iucnredlist.org>. Downloaded on 03 December 2014.
  40. 40. Müller T, Balován B, Tatár S, Müllerné-Trenovszki M, Urbányi B, Demény F (2011) Propagation and rearing of European mudminnow (Umbra krameri) in the interest of natural stock maintenance. Pisces Hungarici 5: 15–20. [in Hungarian with English abstract]
  41. 41. Tatár S, Bajomi B, Tóth B, Demény F, Trenovszki M, Urbányi B, et al. (2015) Integrating reintroduction biology and restoration ecology: habitat reconstruction for marshland fishes. Oryx (in press)
  42. 42. Ambrus A, Sallai Z (2014) Distribution and conservation of the European mudminnow (Umbra krameri Walbaum, 1792) at the range of the Fertő-Hanság National Park (North West Hungary) Pisces Hungarici 8, 97–100. [in Hungarian with English abstract]
  43. 43. Bajomi B, Tatár S, Tóth B, Demény F, Müllerné-Trenovszki M, Urbányi B, et al. (2013) Captive-breeding, re-introduction and supplementation of the European mudminnow in Hungary. pp 15–20. In Soorae P. S. (ed.). Global Re-introduction Perspectives: 2013. Further case studies from around the globe. Gland, Switzerland: IUCN/SSC Re-introduction Specialist Group and Abu Dhabi, UAE: Environment Agency-Abu Dhabi. xiv + 282 pp.
  44. 44. Harka Á (1999) Adatok a lápi póc (Umbra krameri) újabb magyarországi lelőhelyeiről. Halászat 92(3): 119–120. [in Hungarian]
  45. 45. Palsbøll PJ, Berube M, Allendorf FW (2007) Identification of management units using population genetic data. Trends Ecol Evol, 22(1): 11–16. pmid:16982114
  46. 46. Nagababu E, Lakshmaiah N (1992) Inhibitory effect of eugenol on non-enzymatic lipid peroxidation in rat liver mitochondria. Biochemical Pharmacology 43:2393–2400. pmid:1319160
  47. 47. Anderson WG, McKinley RS, Colavecchia M (1997). The use of clove oil as an anesthetic for rainbow trout and its effects on swimming performance. North American Journal of Fisheries Management, 17(2), 301–307.
  48. 48. Winkler KA, Weiss S (2009) Nine new tetranucleotide microsatellite DNA markers for the European mudminnow Umbra krameri. Conserv Genet 10(4):1155–1157.
  49. 49. Rousset F (2008) Genepop'007: a complete reimplementation of the Genepop software for Windows and Linux. Mol Ecol Res 8: 103–106.
  50. 50. Van Oosterhout C, Hutchinson WF, Wills DP, Shipley P (2004). MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Molecular Ecology Notes, 4(3), 535–538.
  51. 51. Peakall R, Smouse PE (2012) GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics 28: 2537–2539. pmid:22820204
  52. 52. Lynch M, Ritland K (1999) Estimation of pairwise relatedness with molecular markers. Genetics 152:1753–1766. pmid:10430599
  53. 53. Chapuis MP, Estoup A (2007) Microsatellite null alleles and estimation of population differentiation. Mol Biol Evol 24: 621–631. pmid:17150975
  54. 54. Dempster AP, Laird NM, Rubin DB (1977) Maximum likelihood from incomplete data via the EM algorithm. J Roy. Stat. Soc. B 39:1–38.
  55. 55. Cavalli-Sforza LL, Edwards AWF (1967) Phylogenetic analysis. Models and estimation procedures. Am J Hum Genet 19(3 Pt 1): 233–257. pmid:6026583
  56. 56. Vähä JP, Erkinaro J, Niemela E, Primmer CR (2007) Life-history and habitat features influence the within-river genetic structure of Atlantic salmon. Mol Ecol 16: 2638–2654. pmid:17594436
  57. 57. Pritchard JK, Wen X, Falush D (2010) Documentation for STRUCTURE software: version 2.3.
  58. 58. Dent EA, von Holdt BM (2012) STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour 4: 359–361.
  59. 59. Evanno G, Regnaut S, Goudet J (2005) Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol 14: 2611–2620. pmid:15969739
  60. 60. Jakobsson M, Rosenberg NA (2007) CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics 23: 1801–1806. pmid:17485429
  61. 61. Rannala B, Mountain JL (1997) Detecting immigration by using multilocus genotypes. Proc. Natl. Acad. Sci. USA 94: 9197–9221. pmid:9256459
  62. 62. Piry S, Alapetite A, Cornuet JM, Paetkau D, Baudouin L, Estoup A (2004) GeneClass2: A Software for Genetic Assignment and First-Generation Migrant Detection. J Hered 95: 536–539. pmid:15475402
  63. 63. Beerli P, Felsenstein J (1999) Maximum-likelihood estimation of migration rates and effective population numbers in two populations using a coalescent approach. Genetics 152: 763–773. pmid:10353916
  64. 64. Beerli P, Felsenstein J (2001) Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach. Proc Natl Acad Sci USA 98: 4563–4568. pmid:11287657
  65. 65. Mantel N (1967) The detection of disease clustering and a generalized regression approach. Cancer Res 27: 209–220. pmid:6018555
  66. 66. Legendre P (1993) Spatial autocorrelation: trouble or new paradigm? Ecology 74: 1659–1673.
  67. 67. Smouse PE, Peakall R (1999) Spatial autocorrelation analysis of individual multiallele and multilocus genetic structure. Heredity 82: 561–573. pmid:10383677
  68. 68. Diniz-Filho JAF, De Campos Telles MP (2002) Spatial autocorrelation analysis and the identification of operational units for conservation in continuous populations. Conserv Biol 16: 924–935.
  69. 69. Peakall R, Ruibal M, Lindenmayer DB (2003) Spatial autocorrelation analysis offers new insights into gene flow in the Australian bush rat, Rattus fuscipes. Evolution. 57: 1182–1195. pmid:12836834
  70. 70. Harka Á, Sallai Z, Kosco J (2003) Az amurgéb (Perccottus glenii) terjedese a Tisza vizrendszereben, A Puszta, vol. 18, p. 49–56 [in Hungarian].
  71. 71. Hajdú J, Kováč V (2002) Fishfauna of several water bodies in the Large Danube Island area. Folia Faun Slovaca 7: 75–81.
  72. 72. Bij de Vaate A, Jazdzewski K, Ketelaars HA, Gollasch S, Van der Velde G (2002) Geographical patterns in range extension of Ponto-Caspian macroinvertebrate species in
  73. 73. Hewitt GM (1999) Post‐glacial re‐colonization of European biota. Biol J Linn Soc 68(1‐2): 87–112.
  74. 74. Dudgeon D, Arthington AH, Gessner MO, Kawabata ZI, Knowler DJ, et al. (2006) Freshwater biodiversity: importance, threats, status and conservation challenges. Biol Rev 81: 163–182. pmid:16336747
  75. 75. Gy Gábris, A Nádor (2007) Long-term fluvial archives in Hungary: response of the Danube and Tisza rivers to tectonic movements and climatic changes during the Quaternary: a review and new synthesis Quaternary Sci Rev 26/22-2: 2758–2782.
  76. 76. Timár G, Gábris Gy (2008) Estimation of water conductivity of the natural flood channels on the Tisza flood-plain, the Great Hungarian Plain. Geomorphology 98: 250–261.
  77. 77. Lászlóffy W (1982) A Tisza (vizi munkálatok és vízgazdálkodás a tiszai vízrendszerekben). Akadémiai Kiadó, Budapest, 610. [in Hungarian].
  78. 78. Zs Szabó (2010) GIS analysis of effects of inland water channels reconstruction on a plain sample area—Agricultural Informatics. 1(3):40–47 [in Hungarian with English abstract]
  79. 79. Leél-Őssy S (1953) A Rákosvidék geomorfológiája—Geomorphology of the Rákos region. Földrajzi Értesítő 2(1): 70–86. [in Hungarian]
  80. 80. Ruszkiczay-Rüdiger Zs, Fodor L, Horváth E, Telbisz T (2007) Folyóvízi, eolikus és neotektonikai hatások szerepe a Gödöllői-dombság felszínfejlődésében—DEM-alapú morfometriai vizsgálat. Földrajzi Közlemények 131/4: 319–342. [in Hungarian with English abstract]
  81. 81. Gábris Gy, Mari L (2007) A Zala folyó lefejezése—The Pleistocene beheading of the Zala River (West Hungary). Földrajzi értesítő 56(1–2): 39–50. [In Hungarian with an English summary]
  82. 82. Sallai Z, Müller T (2014) A lápi póc. In: Veszélyeztetett lápi halak megóvása (lápi póc, réticsík, széles kárász). Vármédia Print Kft, Gödöllő, pp. 11–84. ISBN 978-963-269-428-3. [In Hungarian]