Spatial and temporal variation in population genetic structure of wild Nile tilapia (Oreochromis niloticus) across Africa
© Bezault et al; licensee BioMed Central Ltd. 2011
Received: 15 July 2011
Accepted: 9 December 2011
Published: 9 December 2011
Reconstructing the evolutionary history of a species is challenging. It often depends not only on the past biogeographic and climatic events but also the contemporary and ecological factors, such as current connectivity and habitat heterogeneity. In fact, these factors might interact with each other and shape the current species distribution. However, to what extent the current population genetic structure reflects the past and the contemporary factors is largely unknown. Here we investigated spatio-temporal genetic structures of Nile tilapia (Oreochromis niloticus) populations, across their natural distribution in Africa. While its large biogeographic distribution can cause genetic differentiation at the paleo-biogeographic scales, its restricted dispersal capacity might induce a strong genetic structure at micro-geographic scales.
Using nine microsatellite loci and 350 samples from ten natural populations, we found the highest genetic differentiation among the three ichthyofaunal provinces and regions (Ethiopian, Nilotic and Sudano-Sahelian) (RST = 0.38 - 0.69). This result suggests the predominant effect of paleo-geographic events at macro-geographic scale. In addition, intermediate divergences were found between rivers and lakes within the regions, presumably reflecting relatively recent interruptions of gene flow between hydrographic basins (R ST = 0.24 - 0.32). The lowest differentiations were observed among connected populations within a basin (R ST = 0.015 in the Volta basin). Comparison of temporal sample series revealed subtle changes in the gene pools in a few generations (F = 0 - 0.053). The estimated effective population sizes were 23 - 143 and the estimated migration rate was moderate (m ~ 0.094 - 0.097) in the Volta populations.
This study revealed clear hierarchical patterns of the population genetic structuring of O. niloticus in Africa. The effects of paleo-geographic and climatic events were predominant at macro-geographic scale, and the significant effect of geographic connectivity was detected at micro-geographic scale. The estimated effective population size, the moderate level of dispersal and the rapid temporal change in genetic composition might reflect a potential effect of life history strategy on population dynamics. This hypothesis deserves further investigation. The dynamic pattern revealed at micro-geographic and temporal scales appears important from a genetic resource management as well as from a biodiversity conservation point of view.
The population genetic structure of living organisms is largely shaped by both historical and contemporary gene flow in the species range . Furthermore, the factors structuring the genetic diversity are also able to play a key role in the process of diversification, adaptation and speciation . Identifying the interactions between those factors is crucial to understand the evolutionary history of a species . On the one hand, the extrinsic factors, such as climatic and geological events, shape the large-scale population differentiation pattern [4–6]. On the other hand, both extrinsic (e.g. habitat heterogeneity), and intrinsic factors (e.g. dispersal capability, mating system and habitat preference) have an impact on gene pool composition at intra-population level or over short time periods [4, 7–9]. Species with wide distribution and high dispersal capability are supposed to exhibit genetic differentiation at the paleo-biogeographic scales, with limited micro-geographic structure. On the contrary, species with limited distribution and restricted dispersal capability and/or specific mating behaviour are supposed to exhibit strong genetic structure both at micro-geographic and temporal scales. However, species with a wide distribution area and specific life history traits (e.g. restricted dispersal, limited population size) are expected to exhibit more complex genetic diversity pattern.
Cichlid fishes are well known examples of complex population genetic structure among African ichthyofauna. Lake-wide studies have revealed the impact of paleo-historical lake level fluctuations on genetic diversity (e.g.). Small scale investigations have revealed the importance of both habitat heterogeneity, ecology and dispersal capability on population differentiation (e.g.). To date, however, very few studies have investigated the factors responsible for the population genetic structures both at the micro- and the macro-geographic scales at the same time. They often represent only a relatively restricted paleo-biogeographic scale .
The Nile tilapia, Oreochromis niloticus (Linnaeus, 1758), is an interesting model-species to study the interactions between intrinsic and extrinsic factors on the structure of the natural populations from a local and temporal to a broader biogeographic scale. This economically important fish has one of the largest natural distributions among African fresh-water fishes, covering the entire Nilo-Sudanian province (from Senegal to Nile basins), the Ethiopian Rift Valley province, the Kivu province, north Tanganyika province (Ruzizi) and the Northern part of the East African Rift Valley. This species shows an exceptional capacity of adaptation, which allowed its colonisation of a wide range of habitats from small forest rivers to large drainage and lakes, as well as alkaline pools with hot springs [13, 14]. The description of seven sub-species based on eco-morphology  largely reflects their adaptive divergences.
Due to its great interest for aquaculture and fisheries, the Nile tilapia and a few other tilapia species have been introduced outside their natural distributions . Introduced tilapias have become invasive especially in areas originally not containing any tilapiine cichlid, within as well as outside Africa. In area inhabited by congeneric tilapias, on the other hand, they have often led to hybridisation with the local allopatric species . However, between sympatric congeneric species of tilapias, signature of hybridisation has only been reported at evolutionary time scale, not at current ecological time scale. This is especially the case of the ancient mitochondrial introgression from O. aureus into O. niloticus populations restricted to West Africa, whereas the 2 species are also sympatric in the Nile River. This introgression has certainly happened during the drastic water level fluctuations of the Pleistocene . But neither recent natural hybridisation event nor successful translocation of any allopatric tilapia (i.e. Oreochromis spp.) within the natural distribution of O. niloticus has been reported so far .
As most of the cichlid fish, the Nile tilapia exhibits interesting and complex life history traits, and especially well-developed social behaviour [13, 14, 17–19]. During reproduction, males show strong territoriality and females provide elaborated parental care (i.e. maternal mouth-brooding and guarding) . Because of this reproduction behaviour and substrate affinity, the Nile tilapia is considered a rather sedentary species . In addition, the lekking behaviour of males to attract females for reproduction suggests the existence of a certain level of sexual selection [17, 20]. These life history traits are expected to strongly affect the population dynamics of this species, via limited dispersal or reduced effective population size.
From the paleo-geographic point of view, Africa has experienced severe hydrogeographic modifications since the Pleistocene. The East African Rift valley has been subject to many tectonic disruptions of the water basins with inversion of the course of some rivers in the Nile basin (sensu lato), whereas the Sudano-Sahelian region experienced dramatic climatic fluctuations with alternating humid and dry phases [21, 22]. By their drastic modifications of the extension and connectivity of the different water-basins, all these paleo-geographic and climatic events have undoubtedly affected (1) the distribution of fish species in these ichthyofaunal provinces [23, 24] and (2) their population genetic structure.
Previous studies have suggested an influence of paleo-geographic events on the historical distribution of Nile tilapia, based either on morphological traits  or on moderately polymorphic molecular markers, i.e. allozymes and mtDNA, [16, 25–28]. At the opposite, micro-geographic population structure has been reported in lacustrine populations of another tilapiine species, Sarotherodon melanotheron, with similar life history strategy, implying non-random mating . However, to date, the relative importance of these paleo-geographic events and the factors influencing the micro-geographic structure of the current populations is not well understood.
In this study, we investigated the spatio-temporal genetic structure of ten natural populations of O. niloticus, which cover the main part of the species' natural distribution in Africa. By investigating both spatial and temporal genetic diversity based on nine microsatellite loci, we aimed (1) to understand the current population genetic structure of O. niloticus in its natural habitat range in Africa and (2) to evaluate to what extend the population genetic structure was shaped by the paleo-geographic events and the current geographic connectivity at the different spatio-temporal scales. We also discuss potential roles of life history of the species in the species-range population genetic differentiation.
Spatial and temporal sampling
Natural populations of O. niloticus were sampled from ten geographic sites across Africa to maximise the representation of the biodiversity of the species, naturally spread over the vast Nilo-Sudanian and the Ethiopian ichthyofaunal provinces [23, 30]. In order to complement the analyses of population genetic differentiation from macro-geographic to micro-geographic and temporal scales, further sampling was focused along the Volta basin.
Population samples information
Lake or River
O. n. cancellatus
O. n. cancellatus
O. n. filoa
08° 51' N
O. n. vulcani
El Molo Bay
O. n. niloticus
O. n. niloticus
O. n. niloticus
O. n. niloticus
O. n. niloticus
O. n. niloticus
Micro-geographic and temporal sampling on the Volta basin was conducted at three sampling localities along the hydrographic basin, of which two were sampled successively multiple times between 2001 and 2003 (Table 1 and Figure 1). They represent 206 individual samples, which were used to estimate temporal variation of population genetic structure as well as effective size of these populations.
Genetic diversity per population.
A r ± SD
He ± SD
Hobs ± SD
2.15 ± 1.16
0.32 ± 0.28
0.29 ± 0.26
1.68 ± 0.60
0.17 ± 0.21
0.17 ± 0.24
2.31 ± 0.78
0.36 ± 0.23
0.34 ± 0.25
4.75 ± 1.40
0.72 ± 0.10
0.66 ± 0.18
4.94 ± 1.74
0.74 ± 0.11
0.69 ± 0.13
6.30 ± 3.41
0.67 ± 0.27
0.63 ± 0.29
3.00 ± 1.41
0.42 ± 0.24
0.39 ± 0.28
5.30 ± 1.77
0.65 ± 0.18
0.64 ± 0.19
5.79 ± 1.94
0.69 ± 0.19
0.58 ± 0.15
5.12 ± 2.08
0.64 ± 0.24
0.54 ± 0.21
4.49 ± 1.70
0.69 ± 0.17
0.63 ± 0.26
5.40 ± 1.79
0.66 ± 0.19
0.55 ± 0.15
5.90 ± 2.17
0.71 ± 0.17
0.65 ± 0.19
5.37 ± 2.51
0.64 ± 0.19
0.55 ± 0.24
5.35 ± 2.52
0.67 ± 0.16
0.52 ± 0.37
4.19 ± 2.10
0.55 ± 0.20
0.46 ± 0.27
4.84 ± 2.73
0.60 ± 0.25
0.53 ± 0.26
The populations of the Awash basin showed the lower gene diversity based on both allelic richness (A r = 1.68 - 2.31) and heterozygosity (Ho = 0.17 - 0.34). All other Nilotic and Sahelo-Sudanian populations presented an overall higher level of genetic diversity (A r = 4.75 - 6.30; Ho = 0.55 - 0.69), except for Niger revealing an intermediate level of diversity (A r = 3.00; Ho = 0.42).
Matrix of macro-geographic population differentiation based on R ST and F ST estimators
Macro-geographic populations structure & differentiation
The significance (P < 0.001) of the test suggested by Hardy et al.  revealed that SMM-based estimator (R ST , see Methods) is the most pertinent to describe genetic differentiation at macro-geographic scale The highest population differentiations were found between regions (Table 3), especially between Ethiopian and Sudano-Sahelian populations ( = 0.695 ± 0.08), while the Nilotic populations exhibit a high degree of genetic differentiation from either Ethiopian or Sudano-Sahelian populations ( = 0.55 ± 0.11 and = 0.40 ± 0.10, respectively; Figure 1). The hierarchical AMOVA conducted with populations clustered by biogeographic sub-regions revealed that 49.4% of the genetic variance was actually distributed between Ethiopian province and, Nilotic and Sudano-Sahelian regions, whereas 6.7% was distributed among population within regions, both with very high significance (P < 0.001). The alternative AMOVA model, with populations clustered according to taxonomy, showed comparatively less fitted variance partitioning, with lower percentage of variance explained by the clustering factor (47.3%). Although this result was highly significant, a larger part of the genetic variance was attributed to the between-populations within group component (9.1%). When samples were clustered by hydrographic basins, 49.8% of the genetic variability was partitioned among groups, and 3.5% was detected among populations within basins (R SC = 0.069, P < 0.001), revealing intra-basin heterogeneity.
Hierarchical partitioning of the genetic variance.
Variance components (%)
Among pops within groups
All populations in one group
All populations grouped by region 0
All populations grouped by basin 1
All populations grouped by taxonomy 2
within Sudano-Sahelian region
within Sudano-Sahel. grouped by basin3
within Ethiopian region (Awash)
within Nilotic region (Turkana & Nile)
within Volta basin
Micro-geographic & temporal scale
within Volta basin
within Volta grouped by site 4
within Kpandu (all temporal samples)
within Nyinuto (all temporal samples)
Overall, these results demonstrate that ichthyofaunal provinces definition and hydrographical basins boundaries are the two major structuring factors of O. niloticus populations at large geographic scales.
Micro-geographic and temporal comparisons
The analysis of micro-geographic and temporal variations has been conducted along the Volta basin. Based on factorial correspondence analysis (Figure 3e), three groups visibly emerged in the first factorial plan (F1 × F2), clustering the individuals according to their different geographic origins, Kou, Kpandu and Nyinuto. Temporal samples from Kpandu and Nyinuto were very close to each other within each location. Similarly, the STRUCTURE analysis reveals an optimal number of four clusters, representing preferentially the three different geographic locations, without strict discrimination (Figure 2b). The fourth genetic cluster was relatively equally distributed across all sites. No clear differentiation among temporal samples was observed at K = 4 or higher number of cluster (see Additional file 3).
Matrix of micro-geographic and temporal population differentiation based on F ST estimators.
The highest genetic differentiation was found between geographic sites in the Volta basin (FST = 0.05 - 0.13, = 0.09 ± 0.02; Figure 1), representing 6.5% of the genetic variance of the system (Table 4). The genetic differentiations within temporal series in Kpandu (F = 0 - 0.04, = 0.02 ± 0.02) and Nyinuto (F = 0 - 0.05, = 0.03 ± 0.03) were lower than observed between sample sites, but yet significant (Table 5). Partition of the genetic variance revealed that 2.1% of the genetic variance was due to genetic differentiation in time (across both Kpandu and Nyinuto; Table 4). To examine whether the genetic differentiation between temporal series was due to the effect of one sample site or both, we performed separate AMOVA on Kpandu and Nyinuto samples, respectively. Interestingly, we found that genetic variance due to among temporal series variability was significant for both sites: 1.3% for Kpandu (FST = 0.013, P < 0.05) and 7.4% for Nyinuto (FST = 0.074, P < 0.01).
This demonstrates that subtle temporal changes of gene pool composition need to be considered in population structuring of O. niloticus at micro-geographic scale.
Effective population size and immigration rate
Effective population size and dispersal rate.
Phylogeography of the Nile tilapia: past history and current genetic structure at macro-geographic scale
From the Miocene (i.e. 20-5 My BP) to the Pleistocene (i.e. 2 My to 10,000 yrs BP), numerous paleo-geographic (e.g. causing disruption of river drainages) and paleo-climatic events (e.g. humid and dry phase cycles), by changing the connectivity between the different hydrographical systems, have deeply affected African ichthyofauna [22–24, 40]. Biogeographic data and fossil records suggest a Nilotic origin of the Nile tilapia (Upper Pliocene ~ 5-2 My BP ) - compatible with allozyme data . The present study identified three genetic clusters of Nile tilapia populations corresponding to the major biogeographic subdivisions: (1) the Ethiopian Rift Valley ichthyofaunal province (i.e. the Awash system), (2) the Nilotic region (sensu lato, comprising the presently separated but previously connected Nile River and Turkana basin) and (3) the Sudano-Sahelian region (i.e. West-African basins) (Figure 1), the two later belonging to the Nilo-Sudanian ichthyofaunal province . The comparison of the AMOVA's results revealed a first differentiation of population between ichthyofaunal regions (49%), as also depicted by STRUCTURE analysis. Therein genetic variance was further partitioned between hydrographic basins (1-3% of genetic variance between basins independently of ichthyological regions), and further remaining significant molecular variance among populations within basins (3.5%). This pattern is consistent with those observed in the few other fish species widespread across the Nilo-Sudanian area and studied by genetic markers, Tilapia zillii, Oreochromis aureus and Clarias gariepinus. In particular, these studies also emphasise the deep split between the Nile and the Sudano-Sahelian populations. The westward spread of this species from the Nile to Sudano-Sahelian rivers has probably occurred during the Pleistocene [13, 14], presumably facing the raise of mega-paleo-lakes and drainage connectivity in the Saharan region [21, 24, 40].
The Ethiopian populations exhibited the lowest levels of polymorphism within each population and the largest genetic divergence from the other ones (Tables 2). The isolation of the Awash from the Nile system occurred by tectonic disruption of this river basin about 12,000-8,000 years ago , which can explain the substantial genetic divergence between these two basins ( = 0.55 ± 0.011 - Figure 1). The low genetic diversity in all three Awash populations (Table 2) and rather high differentiation ( = 0.32 ± 0.22 - Figure 1) among them suggest strong effects of genetic drift within this "islands system", which exhibits an extensive level of fragmentation between the Awash River and several relatively isolated highland lakes. Following Trewavas , two different morphological sub-species coexist in the Awash basin: O. n. cancellatus in the river and cold-water lakes (e.g. Lakes Hora and Koka), and O. n. filoa endemic to the hot spring areas (e.g. Lake Metahara). The lower genetic distance between Hora and Metahara (RST = 0.07) compared to Koka (RST = 0.43 and 0.47 respectively) does not reflect their taxonomic status (i.e. paraphyly of O. n. cancellatus). In addition, the description of O. n. filoa is mostly based on the reduced number of some meristic characters (i.e. number of vertebra, scales), which might be subject to developmental plasticity . These results, as well as the similar reproduction pattern and the absence of post-zygotic barrier [42, 43] support the hypothesis that O. n. filoa and O. n. cancellatus are two ecotypes, rather than distinct subspecies. However, further specific experiments, especially according to environment-induced developmental plasticity and mate choice preference, should be conducted to confirm this hypothesis.
In the Nilotic system sensu lato, disconnection occurred between the Nile river and the currently endorheic Lake Turkana basin about 8,000 years ago [21, 44]. Individuals from these two geographically distant populations have been previously identified as two distinct subspecies, O. n. niloticus and O. n. vulcani respectively . However, compared to the large geographic distance and separation time existing between these two populations, the genetic differentiation was expected to be higher than observed (RST = 0.26 - Figure 1) - with respect to other values observed in the present study. Connections between Lake Turkana and the Nile River have probably existed at least periodically until a few thousand years ago , as it seems to be reflected by nuclear and mtDNA data [16, 25–27].
The Sudano-Sahelian basins were isolated from the Nile River by multiple paleo-geographic events between 7,500 and 6,000 years ago . Although morphological analyses have clustered populations from the Nile River with those of the Sudano-Sahelian basins , previous allozyme data [25, 27] and this study ( = 0.38 ± 0.07) revealed a clear disconnection between these two biogeographic entities. This independent history is reinforced by the differential pattern of ancient introgression of mtDNA from O. aureus to O. niloticus involving the whole West African area, but not reported in any Nilotic population .
The Chad-Chari system appears to have played a key role between the Nilotic and Sudanian biogeographic entities, providing an unidirectional connectivity for the ichthyofauna from the Nile to the West African drainages, via the Niger, during the periods of growth of the mega-paleo-Lake Chad [21, 24, 40, 45]. The West African populations are the most geographically widespread, but also the least morphologically diversified . They also exhibited the lowest level of genetic differentiation between basins ( = 0.24 ± 0.12 - Figure 1). Such a pattern might be explained by the persistence of episodic connections among Senegal, Niger and Volta headwater between 6,000 and 4,500 years ago . Despite these recent gene flow events, all the Sudano-Sahelian populations could be strictly distinguished using microsatellites, even for populations at an intra-basin level, such as for the Volta River, contrary to the previous results obtained using allozyme or mitochondrial data .
As expected, no signature of anthropogenic impact on population differentiation was identified in this study. Paleo-geographic and climatic events appear to be important factors for the current pattern of the genetic structure of Nile tilapia populations within its natural distribution across Africa.
Microgeography and temporal variation
Within the Volta basin, most of the genetic variability was distributed within populations (93.7%), from the upstream Kou River, to the central Lake Volta (Kpandu) and downstream Volta River (Nyinuto). This suggests their recent differentiation and/or non-interrupted gene flow. The genetic differentiation within temporal series in Kpandu (= 0.02 ± 0.02) and Nyinuto (= 0.03 ± 0.03) was lower than between sample sites ( = 0.09 ± 0.02; Figure 1), but still significant and representing more than 2% of the genetic variance of the system. These results suggest a geographic differentiation between populations, and a quick turnover of the gene pools within population within only a few generations.
Our estimates of the effective population size from three different approaches consistently showed larger N e in Kpandu than in Nyinuto (Table 6). The estimated migration rate (m), obtained jointly with Neopen, provided similar values for both populations. Although the size of each sample for these estimates was limited (n = 28 - 30), the congruence of the different N e estimates together with relatively small confidence intervals indicates a good coherence and reliability of our estimations [39, 46]. The population of Kpandu, located at the shore of Lake Volta (one of the largest artificial lakes in the world), is believed to have a large census size and/or well connected to neighbouring populations. In contrast, the population of Nyinuto, located at the downstream limit of Nile tilapia distribution along the Volta River, is likely to have a more limited size and to be easily subject to fluctuations (e.g. variation of water level and salinity - E. Bezault pers. obs.), and/or more isolated, than a central basin population. Despite the potential limit associated to our sampling strategy, our N e estimates consistently showed relatively small effective population size of the natural Nile tilapia populations relative to expected census population size (i.e. fisheries records expressed in millions of tonnes per year for the Volta basin; E. K. Abban, pers. comm.) and a non-negligible dispersal rate.
The hypothesis of the effect of specific social and reproductive behaviour of the Nile tilapia emerges from the positive FIS observed in both populations and the significant relatedness detected in some temporal samples (i.e. 3 out of the 4 temporal samples with significant relatedness showing either a higher mean or variance relatedness distribution, revealing respectively overall related individuals or subgroups of related individuals). A shoaling behaviour involving relatives could for instance explain such a pattern. This behaviour has been observed in other populations of O. niloticus[47, 48] as well as in another tilapia species, Sarotherodon melanotheron[29, 47]. These authors suggested this pattern to be due to direct inbreeding, implying preferential reproduction between individuals related among themselves, as observed in other cichlid species . Furthermore, it has been claimed that this behaviour appears preferentially in open-water habitats, and especially in lagoon systems . In the present study, the pattern of inbreeding observed in both lake and river habitats cannot be totally explained by shoaling behaviour among relatives, due also to the absence of sub-structure related to seine net catches (i.e. potential shoal) within temporal samples. A more direct effect of a non-panmictic reproductive system should be considered. Polygynandric mating system, generally exhibited, by maternal mouthbrooding cichlids, as O. niloticus, is generally associated with sexual selection and pronounced between-individual variations in reproduction success [17, 20]. Such variance of fitness has been reported for male Nile tilapia but only in captivity . The rapid turnover in gene pool and the restricted effective population sizes, observed in our study, could be explained by the differential contribution of breeders to within and between cohorts over generations.
Furthermore the non-negligible dispersal rate detected in the present study opens up new research perspectives on cichlids. While sophisticated reproduction strategy and extended parental care are expected to increase territoriality and consequently prevent adult migration, our results suggest that the exceptional adaptive plasticity and riverine affinity of the tilapia could also be seen in its dispersal leading to larger migration rate than more specialised (stenotypic) cichlid species from East African Great Lakes [11, 20, 51]. This hypothesis is corroborated by the population genetic analysis of the riverine haplochromines generalist, Pseudocrenilabrus multicolour victoriae, presenting relatively large distribution area with geographic structure being the main driver of population structure and gene-flow . It has also been demonstrated in this species important change of gene pool composition at short temporal scale (i.e. 2 years), which has been associated with drastic fluctuation of hydrographic regime . Larger estimates of effective population size (including higher confidence interval reaching infinity) have been obtained for P. m. victoriae in this study compared to effective population size estimated here for O. niloticus. Unfortunately no census population size estimate for any of these species is available to effectively compare the N e /N ratio, however, in the case of P. m. victoriae it cannot be ruled out that the larger effective population size observed is influenced by flooding regime.
The distribution and pattern of genetic diversity of a species over its natural range is the result of multiple and interlinked factors, which acted through space and time. For the Nile tilapia, O. niloticus, the genetic pattern observed from large geographic to temporal scales first reflects the hierarchical impact of paleo-geographic and climatic events as well as the impact of hydrographic connectivity on population genetic structure. Furthermore, the hierarchical pattern of genetic differentiation seems to reflect only bio-geographic events, without revealing any signature of anthropogenic impact. Compared to the most studied, highly radiating cichlid lineages, the morphologically and ecologically generalist species, as especially the tilapiines, [13, 19], might have distinct population dynamic characteristics. The relatively small effective population size associated with a moderate dispersal and the rapid change in gene pool composition in only a few generations tend to reinforce our current knowledge about the biology of this species, and they might reflect a potential impact of natural history strategy on population dynamics. Further investigations are needed to confirm this hypothesis. In addition, more research emphasis should be put on micro-geographic population structures to evaluate a more precise relationship between census and effective population size, the extent of dispersal among the natural population of Nile tilapia, as well as the balance between selection and gene-flow across different habitats. Furthermore, the possibility that Nile tilapia populations may be organised in meta-populations within larger and more connected systems (like within the Lake Volta) should be further investigated, from a genetic resource as well as from a conservation point of view.
Sampling and analysis strategy
We collected 350 samples of O. niloticus from ten geographic sites across Africa (Table 1 and Figure 1). Our samplings aimed 1) to cover the macro-geographic area in the natural range of this species as extensive as possible and include different ichthyofaunal provinces and regions and therein hybrographic basins , and 2) to include successive temporal sampling through the short period of time in local populations, which allow us to survey the temporal genetic stability at the micro-geographic level. In both cases, a special attention was paid to avoid a potential influence of anthropogenic fish translocations and introductions, which could cause introgression and hybridisation. For that purpose, sampling sites without any record of Oreochromis spp. introduction were selected. In addition, all sampled individuals (considering only fully mature specimens, > 120 mm) were visually checked according to specific morphological characters  and confirmed the absence of sign of hybridisation for specimens analysed [16, 42, 43]. Sampling of Volta and Awash basins were done between 2001 and 2004, whereas samples from Lake Turkana, Lake Manzala, Niger and Senegal were a subset of populations already investigated for enzyme polymorphism and mtDNA [16, 25, 26], which have been considered pure O. niloticus at the nuclear genome level .
The sampling locations covered six hydrographic basins in Africa, belonging to the Sudano-Sahelian (Niger, Senegal and Volta) and Nilotic (Turkana, Nile) regions of the Nilo-Sudanian province and Ethiopian Rift valley (Awash basin) ichthyofaunal province [23, 30, 31] and represented four sub-species (O. n. niloticus, O. n. vulcani, O. n. cancellatus and O. n. filoa) according to their eco-morphological characters .
Analysis of micro-geographic and temporal variation was focussed on the Volta basin, where temporal sample series were collected repeatedly over a period of 18 months in two of the three geographic populations: Kpandu and Nyinuto (Table 1). These temporal samples were collected precisely at the same geographic location (verified by GPS coordinates) with the same standardised fishing protocol using seine nets. To minimise the potential impact of shoals' capture (i.e. biased sampling among relatives), each temporal sample was composed of individuals taken from multiple seine catches.
For our analyses, we considered two different datasets. (1) The macro-geographic analyses were conducted including all samples from the ten geographic populations considered (n = 350; temporal sample series from the two sites of the Volta basin pooled respectively into a single population-sample per site, i.e. Kpandu and Nyinuto). (2) The micro-geographic and temporal analyses were conducted at the scale of the entire Volta basin, including the samples from the 3 sites along the basin (n = 206): one sample from the Kou River (upstream Volta basin), four temporal samples from Kpandu (on the Lake Volta) and three temporal samples from Nyinuto (downstream on the Volta River) (Table 1 and Figure 1).
Nine unlinked microsatellites were selected from the O. niloticus genomic DNA library . Genomic DNA was extracted from fin clips or muscle fragments stored in 95% ethanol using a phenol-chloroform procedure . DNA samples were amplified by PCR using the indirect fluorescent tagging procedure described by Schuelke . The amplification was performed using standard PCR-mix concentration and a touch-down PCR procedure, allowing standardised conditions over the loci analysed (see details in Additional file 1), using a 384-well-plate technology (Genotyping Platform, Genopole Montpellier Languedoc-Roussillon, CIRAD). PCR products were detected on an automated Li-Cor gel sequencer (IR2, Lincoln, Neb.), and genotyped were accessed by eyes.
Genetic diversity Analysis
The total number of alleles (A), observed (Ho) and non-biased expected (Hnb) heterozygosities were calculated for each locus separately and across all loci using GENETIX 4.03  and FSTAT 2.9.3 . Allelic richness (Ar) was estimated for each population using rarefaction algorithm implemented in HP-rare , assuming an even number of genes per population, n = 20, based on the smallest population sample size in the study. For each sample, the fixation index, FIS, was estimated across all loci, and statistical significance was assessed with permutation tests (n = 1000) using GENETIX. Corrections for multiple tests were performed using the sequential Bonferroni procedure .
Heterozygote deficits, as indicated by FIS values significantly larger than 0, might reflect (1) sampling issues (e.g. Wahlund effect) or (2) impact of specific life history traits including a low dispersal or any social behaviour leading to a particular mating system. However, technical problems, such as null alleles, can also affect FIS. Thus, we first ran MICRO-CHECKER to estimate if null alleles might affect the fixation index. Potential null alleles were identified at three loci (UNH142, UNH162, and UNH211), however, excluding these three loci had very minor effects on FIS (data not shown).
We used analysis of molecular variance, AMOVA , implemented in ARLEQUIN 2.0  to test the potential population substructure due to different seine catches within the different temporal and local samples. We tested the possibility of non-random mating and inbreeding within population and cohorts as a potential cause of positive FIS values by calculating the degree of relatedness within samples using the pairwise relatedness index rxy of Queller & Goodnight  implemented in IDENTIX. The observed arithmetic mean of rxy (Mrxy) in a sample was compared to that expected for a panmictic population. The observed variance of rxy (Vrxy) was also compared with the expectation from a panmictic population to test for the presence of multiple family groups within samples (which would increase the variance of rxy) [64, 65]. For the sample pooled by station (Kp and Ny), the null distribution was generated by permuting genotypes within the population (1,000 permutations), following Castric et al.. For the temporal samples, the null distribution was calculated based on 1,000 randomised sub-samples of the same size as the tested sample, taken from the total samples from the given geographic site (i.e. Kp or Ny) .
Population clustering and genetic differentiation
Factorial correspondence analyses were performed using R  and the package ADE4 v1.4-9  to investigate the relationship between populations from multi-locus genotypes without defining groups a priori. Successive correspondence analyses were performed to group populations from the largest to the finest geographic scales investigated.
Similarly, genetic population clustering was assessed using STRUCTURE v.2.3.1[32, 33] and conducted at i) macro-geographic scale and ii) micro-geographic and temporal scales. Simulations were conducted assuming an admixture model, with 10,000 burn-in steps and 100,000 MCMC steps. Ten independent runs were performed for each number of clusters assumed (K), with K varying between 2 and z+1, with z being the number of samples including in the analysis (i.e. respectively 11 for macro-geographic and 9 for micro-geographic and temporal analysis). The optimal number of clusters was assessed based on correction proposed by Evanno et al. .
Genetic parameters, providing a measure of the genetic differentiation existing between two populations, can be calculated using different estimators (for a discussion see ). The F-statistics , based on the infinite alleles model (IAM), relies solely on allele identity information. As it is an easy target for homoplasy, leading to the underestimation of genetic differentiation among highly structured populations (where μ > m), this parameter generally is used to compare populations at small geographic scale. In contrast, the analogous R-statistics  based on the stepwise mutation model (SMM), takes into account allele size differences. As it typically suffers from high sampling variances, it is used for comparisons of large geographic scales. Here, as we were faced with distinct scales, we applied the test suggested by Hardy et al. to choose the most suitable estimators. This test indicates whether or not allele sizes provide information on population differentiation, helping to choose between F- statistics and R- statistics in the most objective way. We applied this test implemented in SPAGeDi  for both datasets, and used the most pertinent statistics. Matrices of pairwise differentiation between populations (R-statistics or F- statistics) were computed using ARLEQUIN. The averaged genetic differentiation calculated over pairwise value between or within group of populations were noted or .
Further, AMOVA  was also used to investigate how genetic variability was distributed in space and time. Considering macro-geographic scale, successive hierarchical levels were defined based alternatively on i) biogeographic regions and provinces (3 groups), ii) hydrographic basins (6 groups) and iii) sub-species taxonomy (4 groups - see Table 4 for details). Considering micro-geographic scales, hierarchical AMOVA was computed considering temporal samples within localities at i) the Volta basin and ii) the locality levels.
Estimating effective population size and migration rate
We estimated effective population size (N e ) and migration rate (m) based on the spatio-temporal samples from the Volta basin. In order to reduce the sampling biases due to the small sample size, we excluded a temporal sample with the smallest sample size (KpM2, n = 13) from the estimation, and then retained only larger samples (n ≥ 28). We used three different approaches to estimate N e (see discussion in [39, 46, 73]). Two of them are based on temporal variation of allele frequencies within populations. Both use likelihood methods and are implemented in MLNe 2.3 , but the first method (N e -CLOSED) assumes no immigration , assumption likely violated in the current context. The second method (N e -open), relaxes this assumption and allows the joint estimation of N e and m. For these approaches, the three temporal samples were considered to be from successive generations, assuming a generation time of about 6 months, as generally accepted for the Nile tilapia . For the estimation of Neopen of a given focal population, we pooled the samples from the other populations from the basin (i.e. Ko + Ny for Kp, and Ko + Kp for Ny) to mimic the allele frequencies in the source population of all potential immigrants. The third approach (LD-N e ), is based on the linkage disequilibrium method  and implemented in NE ESTIMATOR. In this method, no migration is assumed. LD-N e was calculated for each temporal sample and then averaged over them. A point estimate and the 95% confidence intervals were obtained for each population by each method.
We would like to thank Eddie Abban and the team of the Water Research Institute (Ghana), Abebe Getahun (University of Addis Ababa), Pascal Bonnet & Géraud Laval (CIRAD/ILRI), the Office des Parcs Nationaux du Sénégal, Gabriele Hörstgen-Schwark (University of Göttingen), and Brendan McAndrew (University of Stirling) for their cooperation in obtaining samples. We also thank Stéphane Mauger, Gaëlle Chupot, René Guyomard (INRA), and Claire Billot (CIRAD) for their laboratory assistance; Agathe Bezault for the illustrations; Wendy Brand-Williams (INRA) for English improvement, and especially, Nicolas Poulet (ONEMA) for his help in statistical analyses and discussions about demography, Bernard Chevassus (INRA) and Mark Jobling (University of Leicester) for stimulating discussions, and Dylan Fraser (Dalhousie University), Irene Keller (EAWAG), Marta Barluenga (Museum of Natural History, Madrid), the editor Michael Hansen and the anonymous reviewers for their constructive comments. This study was supported by grants from the CIRAD, INRA & MENRT (France) to EB.
- Holsinger KE, Weir BS: Genetics in geographically structured populations: defining, estimating and interpreting FST. Nature Review Genetics. 2009, 10 (9): 639-650. 10.1038/nrg2611.View ArticleGoogle Scholar
- Schluter D: The Ecology of Adaptive Radiation. Oxford: Oxford University Press. 2000Google Scholar
- Weiss S, Persat H, Eppe R, Schlotterer C, Uiblein F: Complex patterns of colonization and refugia revealed for European grayling Thymallus thymallus, based on complete sequencing of the mitochondrial DNA control region. Mol Ecol. 2002, 11 (8): 1393-1407. 10.1046/j.1365-294X.2002.01544.x.View ArticlePubMedGoogle Scholar
- Tibbets CA, Dowling TE: Effects of intrinsic and extrinsic factors on population fragmentation in three species of North American minnows (Teleostei: Cyprinidae). Evolution. 1996, 50 (3): 1280-1292. 10.2307/2410668.View ArticleGoogle Scholar
- McClenaghan LR, Smith MH, Smith MW: Biochemical genetics of mosquitofish. IV. Changes of alleles frequencies through time and space. Evolution. 1985, 39: 451-160. 10.2307/2408377.View ArticleGoogle Scholar
- Bernatchez L, Wilson CC: Comparative phylogeography of nearctic and palearctic fishes. Mol Ecol. 1998, 7 (4): 431-452. 10.1046/j.1365-294x.1998.00319.x.View ArticleGoogle Scholar
- Laroche J, Durand JD, Bouvet Y, Guinand B, Brohon B: Genetic structure and differentiation among populations of two cyprinids, Leuciscus cephalus and Rutilus rutilus, in a large European river. Can J Fish Aquat Sci. 1999, 56 (9): 1659-1667.View ArticleGoogle Scholar
- Heithaus MR, Laushman RH: Genetic variation and conservation of stream fishes: influence of ecology, life history, and water quality. Can J Fish Aquat Sci. 1997, 54 (8): 1822-1836.View ArticleGoogle Scholar
- Salgueiro P, Carvalho G, Collares-Pereira MJ, Coelho MM: Microsatellite analysis of genetic population structure of the endangered cyprinid Anaecypris hispanica in Portugal: implications for conservation. Biol Conserv. 2003, 109 (1): 47-56. 10.1016/S0006-3207(02)00132-5.View ArticleGoogle Scholar
- Egger B, Koblmuller S, Sturmbauer C, Sefc KM: Nuclear and mitochondrial data reveal different evolutionary processes in the Lake Tanganyika cichlid genus Tropheus. BMC Evol Biol. 2007, 7: 1-14. 10.1186/1471-2148-7-1.View ArticleGoogle Scholar
- Wagner CE, McCune AR: Contrasting patterns of spatial genetic structure in sympatric rock-dwelling cichlid fishes. Evolution. 2009, 63 (5): 1312-1326. 10.1111/j.1558-5646.2009.00612.x.View ArticlePubMedGoogle Scholar
- Crispo E, Chapman LJ: Population genetic structure across dissolved oxygen regimes in an African cichlid fish. Mol Ecol. 2008, 17 (9): 2134-2148. 10.1111/j.1365-294X.2008.03729.x.View ArticlePubMedGoogle Scholar
- Trewavas E: Tilapiine fishes of the genera Sarotherodon, Oreochromis and Danakilia. 1983, London, UK: British Museum Natural HistoryGoogle Scholar
- Philippart J-C, Ruwet JC: Ecology and distribution of tilapias. The biology and culture of tilapias. Edited by: Pullin RSV, Lowe-McConnell RH. 1982, Manila, Philippines: ICLARM, 7: 15-60.Google Scholar
- Eknath AE, Hulata G: Use and exchange of genetic resources of Nile tilapia (Oreochromis niloticus). Rev Aquac. 2009, 1 (3-4): 197-213. 10.1111/j.1753-5131.2009.01017.x.View ArticleGoogle Scholar
- Rognon X, Guyomard R: Large extent of mitochondrial DNA transfer from Oreochromis aureus to O. niloticus in West Africa. Mol Ecol. 2003, 12 (2): 435-445. 10.1046/j.1365-294X.2003.01739.x.View ArticlePubMedGoogle Scholar
- Keenleyside MHA: Parental care. Cichlid Fishes: Behaviour, ecology and evolution. Edited by: Keenleyside M. 1991, London, UK: Chapman & Hall, 191-208.Google Scholar
- Lowe-McConnell RH: Ecology of cichlids in South American and African waters, excluding the African Great Lakes. Cichlid Fishes: Behaviour, ecology and evolution. Edited by: Keenleyside M. 1991, London, UK: Chapman & Hall, 60-83.Google Scholar
- Klett V, Meyer A: What, if Anything, is a Tilapia? Mitochondrial ND2 Phylogeny of Tilapiines and the Evolution of Parental Care Systems in the African Cichlid Fishes. Mol Biol Evol. 2002, 19 (6): 865-883.View ArticlePubMedGoogle Scholar
- Kocher TD: Adaptative evolution and explosive speciation: the cichlid fish model. Nat Rev Genet. 2004, 5 (4): 288-298. 10.1038/nrg1316.View ArticlePubMedGoogle Scholar
- Beadle LG: The inland waters of tropical Africa. An introduction to tropical limnology. 1981, New York, USA.: Langman Inc, 2Google Scholar
- Grove AT: The physical evolution of the river basins. The Niger and its neighbours Environmental history and hydrobiology, human use and health hazard of the major West African rivers. Edited by: Grove AT. 1985, Netherlands: A.A. Balkena Publishers, 21-60.Google Scholar
- Roberts TR: Geographical distribution of African freshwater fishes. Zool J Linn Soc. 1975, 57: 249-319. 10.1111/j.1096-3642.1975.tb01893.x.View ArticleGoogle Scholar
- Lévêque C: Diversity dynamics and conservation: the freshwater fish of tropical Africa. 1997, Cambridge, UK: Cambridge University PressGoogle Scholar
- Rognon X, Andriamanga M, McAndrew B, Guyomard R: Allozyme variation in natural and cultured populations in two tilapia species: Oreochromis niloticus and Tilapia zillii. Heredity. 1996, 76: 640-650. 10.1038/hdy.1996.91.View ArticleGoogle Scholar
- Rognon X, Guyomard R: Mitochondrial DNA differentiation among East and West African Nile tilapia populations. J Fish Biol. 1997, 51 (1): 204-207. 10.1111/j.1095-8649.1997.tb02525.x.View ArticlePubMedGoogle Scholar
- Agnèse JF, Adépo-Gourène B, Abban EK, Fermon Y: Genetic differentiation among natural populations of the Nile tilapia Oreochromis niloticus (Teleostei, cichlidae). Heredity. 1997, 79: 88-96. 10.1038/hdy.1997.126.View ArticlePubMedGoogle Scholar
- Seyoum S, Kornfield I: Taxonomic notes on the Oreochromis niloticus subspecies-complex (Pisces, Cichlidae), with a description of a new subspecies. Can J Zool. 1992, 70 (11): 2161-2165. 10.1139/z92-291.View ArticleGoogle Scholar
- Pouyaud L, Desmarais E, Chenuil A, Agnèse JF, Bonhomme F: Kin cohesiveness and possible inbreeding in the mouthbrooding tilapia Sarotherodon melanotheron (Pisces, Cichlidae). Mol Ecol. 1999, 8: 803-812. 10.1046/j.1365-294X.1999.00632.x.View ArticleGoogle Scholar
- Paugy D, Zaiss R, Troubat JJ: Faunafri - World Wide Web electronic publication. 2008, [http://www.poissons-afrique.ird.fr/faunafri/]Google Scholar
- Paugy D: The Ethiopian subregion fish fauna: an original patchwork with various origins. Hydrobiologia. 2010, 649 (1): 301-315. 10.1007/s10750-010-0273-0.View ArticleGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155 (2): 945-959.PubMed CentralPubMedGoogle Scholar
- Falush D, Stephens M, Pritchard JK: Inference of population structure using multilocus genotype data: Linked loci and correlated allele frequencies. Genetics. 2003, 164 (4): 1567-1587.PubMed CentralPubMedGoogle Scholar
- Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software structure: a simulation study. Mol Ecol. 2005, 14 (8): 2611-2620. 10.1111/j.1365-294X.2005.02553.x.View ArticlePubMedGoogle Scholar
- Hardy OJ, Charbonnel N, Freville H, Heuertz M: Microsatellite allele sizes: A simple test to assess their significance on genetic differentiation. Genetics. 2003, 163 (4): 1467-1482.PubMed CentralPubMedGoogle Scholar
- Wang JL, Whitlock MC: Estimating effective population size and migration rates from genetic samples over space and time. Genetics. 2003, 163 (1): 429-446.PubMed CentralPubMedGoogle Scholar
- Hill WG: Estimation of Effective Population-Size from Data on Linkage Disequilibrium. Genet Res. 1981, 38 (3): 209-216. 10.1017/S0016672300020553.View ArticleGoogle Scholar
- Wang JL: A pseudo-likelihood method for estimating effective population size from temporally spaced samples. Genet Res. 2001, 78 (3): 243-257. 10.1017/S0016672301005286.View ArticlePubMedGoogle Scholar
- Fraser DJ, Hansen MM, Ostergaard S, Tessier N, Legault M, Bernatchez L: Comparative estimation of effective population sizes and temporal gene flow in two contrasting population systems. Mol Ecol. 2007, 16 (18): 3866-3889. 10.1111/j.1365-294X.2007.03453.x.View ArticlePubMedGoogle Scholar
- Drake NA, Blench RM, Armitage SJ, Bristow CS, White KH: Ancient watercourses and biogeography of the Sahara explain the peopling of the desert. Proc Natl Acad Sci USA. 2011, 108 (2): 458-462. 10.1073/pnas.1012231108.PubMed CentralView ArticlePubMedGoogle Scholar
- Rognon X, Teugels GG, Guyomard R, Galbusera P, Andriamanga M, Volckaert F, Agnese JF: Morphometric and allozyme variation in the African catfishes Clarias gariepinus and C. anguillaris. J Fish Biol. 1998, 53 (1): 192-207.Google Scholar
- Bezault E: Etude du système de déterminisme du sexe au sein de populations naturelles de tilapia du Nil, Oreochromis niloticus : Importance des composantes génétiques et environnementales. PhD thesis. 2005, Orsay, France: Université Paris XIGoogle Scholar
- Bezault E, Clota F, Derivaz M, Chevassus B, Baroiller JF: Sex determination and temperature-induced sex differentiation in three natural populations of Nile tilapia (Oreochromis niloticus) adapted to extreme temperature conditions. Aquaculture. 2007, 272: S3-S16.View ArticleGoogle Scholar
- Grove AT, Goudie AS: Late quaternary lake levels in the rift valley of southern Ethiopia and elsewhere in tropical Africa. Nature. 1971, 234: 403-405. 10.1038/234403a0.View ArticleGoogle Scholar
- Greenwood PH: Fish fauna of the Nile. The Nile, biology of the ancient river. Edited by: Rzóska J. 1976, The Hague, Netherlands, 127-141.View ArticleGoogle Scholar
- Palstra FP, Ruzzante DE: Genetic estimates of contemporary effective population size: what can they tell us about the importance of genetic stochasticity for wild population persistence?. Mol Ecol. 2008, 17: 3428-3447. 10.1111/j.1365-294X.2008.03842.x.View ArticlePubMedGoogle Scholar
- Agnèse JF, Adépo-Gourène B, Nyingi D: Functional microsatellite and possible selective sweep in natural populations of the black-chinned tilapia Sarotherodon melanotheron (Teleostei, Cichlidae). Marine Genomics. 2009, 1 (3-4): 103-107.View ArticleGoogle Scholar
- Nyingi D, De Vos L, Aman R, Agnèse J-F: Genetic characterization of an unknown and endangered native population of the Nile tilapia Oreochromis niloticus (Linnaeus, 1758) (Cichlidae; Teleostei) in the Loboi Swamp (Kenya). Aquaculture. 2009, 297 (1-4): 57-63. 10.1016/j.aquaculture.2009.09.017.View ArticleGoogle Scholar
- Thunken T, Bakker TCM, Baldauf SA, Kullmann H: Active inbreeding in a cichlid fish and its adaptive significance. Curr Biol. 2007, 17 (3): 225-229. 10.1016/j.cub.2006.11.053.View ArticlePubMedGoogle Scholar
- Fessehaye Y, El-Bialy Z, Rezk MA, Crooijmans R, Bovenhuis H, Komen H: Mating systems and male reproductive success in Nile tilapia (Oreochromis niloticus) in breeding hapas: A microsatellite analysis. Aquaculture. 2006, 256 (1-4): 148-158. 10.1016/j.aquaculture.2006.02.024.View ArticleGoogle Scholar
- Danley PD, Markert JA, Arnegard ME, Kocher TD: Divergence with gene flow in the rock-dwelling cichlids of Lake Malawi. Evolution. 2000, 54 (5): 1725-1737.View ArticlePubMedGoogle Scholar
- Crispo E, Chapman LJ: Temporal Variation in Population Genetic Structure of a Riverine African Cichlid Fish. J Hered. 2010, 101 (1): 97-106. 10.1093/jhered/esp078.View ArticlePubMedGoogle Scholar
- Lee WJ, Kocher TD: Microsatellites DNA markers for genetics mapping in Oreochromis niloticus. J Fish Biol. 1996, 49 (1): 169-171.Google Scholar
- Estoup A, Martin O: Marqueurs microsatellites: isolement à l'aide de sondes non-radioactives, caractérisation et mise au point. 1996, [http://www.agroparistech.fr/svs/genere/microsat/microsat.htm]Google Scholar
- Schuelke M: An economic method for the fluorescent labeling of PCR fragments. Nature Biotechnology. 2000, 18 (2): 233-234. 10.1038/72708.View ArticlePubMedGoogle Scholar
- Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F: GENETIX 4.02, logiciel sous Windows pour la génétique des populations. 2001, Montpellier, France: Laboratoire Génome, Populations, Interactions, CNRS UMR 5000, Université de Montpellier IIGoogle Scholar
- Goudet J: FSTAT, a program to estimate and test gene diversities and fixation indices (version 2.9.3). Lausanne, Switzerland: University of Lausanne. 2001Google Scholar
- Kalinowski ST: HP-RARE 1.0: a computer program for performing rarefaction on measures of allelic richness. Mol Ecol Notes. 2005, 5 (1): 187-189. 10.1111/j.1471-8286.2004.00845.x.View ArticleGoogle Scholar
- Rice WR: Anylizing tables of statistical tests. Evolution. 1989, 43 (1): 223-225. 10.2307/2409177.View ArticleGoogle Scholar
- van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P: MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004, 4 (3): 535-538. 10.1111/j.1471-8286.2004.00684.x.View ArticleGoogle Scholar
- Excoffier L, Smouse PE, Quattro JM: Analysis of molecular variance inferred from metric distances among DNA haplotypes - application to human mitochondrial-DNA restriction data. Genetics. 1992, 131 (2): 479-491.PubMed CentralPubMedGoogle Scholar
- Schneider S, Roessli D, Excoffier L: Arlequin Ver 2.000: A software for population genetics data analysis. 2000, Geneva, SwitzerlandGoogle Scholar
- Queller DC, Goodnight KF: Estimating relatedness using genetic markers. Evolution. 1989, 43: 258-275. 10.2307/2409206.View ArticleGoogle Scholar
- Belkhir K, Castric V, Bonhomme F: IDENTIX, a software to test for relatedness in a population using permutation methods. Mol Ecol Notes. 2002, 2 (4): 611-614. 10.1046/j.1471-8286.2002.00273.x.View ArticleGoogle Scholar
- Castric V, Bernatchez L, Belkhir K, Bonhomme F: Heterozygote deficiencies in small lacustrine populations of brook charr Salvelinus Fontinalis Mitchill (Pisces, Salmonidae): a test of alternative hypotheses. Heredity. 2002, 89: 27-35.View ArticlePubMedGoogle Scholar
- Behrmann-Godel J, Gerlach G, Eckmann R: Kin and population recognition in sympatric Lake Constance perch (Perca fluviatilis L.): can assortative shoaling drive population divergence?. Behav Ecol Sociobiol. 2006, 59 (4): 461-468. 10.1007/s00265-005-0070-3.View ArticleGoogle Scholar
- R Development Core Team: R: A language and environment for statistical computing, reference index version 2.2.1. 2005, Vienna, Austria: R Foundation for Statistical ComputingGoogle Scholar
- Chessel D, Dufour AB, Thioulouse J: The ade4 package - I: One-table methods. R News. 2004, 4: 5-10.Google Scholar
- Balloux F, Lugon-Moulin N: The estimation of population differentiation with microsatellite markers. Mol Ecol. 2002, 11 (2): 155-165. 10.1046/j.0962-1083.2001.01436.x.View ArticlePubMedGoogle Scholar
- Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.View ArticleGoogle Scholar
- Slatkin M: A measure of population subdivision based on microsatellite allele frequencies. Genetics. 1995, 139 (1): 457-462.PubMed CentralPubMedGoogle Scholar
- Hardy OJ, Vekemans X: SPAGeDi: a versatile computer program to analyse spatial genetic structure at the individual or population levels. Mol Ecol Notes. 2002, 2 (4): 618-620. 10.1046/j.1471-8286.2002.00305.x.View ArticleGoogle Scholar
- Wang JL: Estimation of effective population sizes from data on genetic markers. Philosophical Transactions of the Royal Society B-Biological Sciences. 2005, 360 (1459): 1395-1409. 10.1098/rstb.2005.1682.PubMed CentralView ArticleGoogle Scholar
- Peel D, Ovenden JR, Peel SL: Neestimator: Software for Estimating Effective Population Size (version 1.3). 2004, [http://www.dpi.qld.gov.au/cps/rde/dpi/hs.xsl/28_6908_ENA_HTML.htm]Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.