Skip to main content

Molecular genetic diversity and differentiation of Nile tilapia (Oreochromis niloticus, L. 1758) in East African natural and stocked populations

Abstract

Background

The need for enhancing the productivity of fisheries in Africa triggered the introduction of non-native fish, causing dramatic changes to local species. In East Africa, the extensive translocation of Nile tilapia (Oreochromis niloticus) is one of the major factors in this respect. Using 40 microsatellite loci with SSR-GBS techniques, we amplified a total of 664 individuals to investigate the genetic structure of O. niloticus from East Africa in comparison to Ethiopian and Burkina Faso populations.

Results

All three African regions were characterized by independent gene-pools, however, the Ethiopian population from Lake Tana was genetically more divergent (Fst = 2.1) than expected suggesting that it might be a different sub-species. In East Africa, the genetic structure was congruent with both geographical location and anthropogenic activities (Isolation By Distance for East Africa, R2 = 0.67 and Uganda, R2 = 0.24). O. niloticus from Lake Turkana (Kenya) was isolated, while in Uganda, despite populations being rather similar to each other, two main natural catchments were able to be defined. We show that these two groups contributed to the gene-pool of different non-native populations. Moreover, admixture and possible hybridization with other tilapiine species may have contributed to the genetic divergence found in some populations such as Lake Victoria. We detected other factors that might be affecting Nile tilapia genetic variation. For example, most of the populations have gone through a reduction in genetic diversity, which can be a consequence of bottleneck (G-W, < 0.5) caused by overfishing, genetic erosion due to fragmentation or founder effect resulting from stocking activities.

Conclusions

The anthropogenic activities particularly in the East African O. niloticus translocations, promoted artificial admixture among Nile Tilapia populations. Translocations may also have triggered hybridization with the native congenerics, which needs to be further studied. These events may contribute to outbreeding depression and hence compromising the sustainability of the species in the region.

Background

Nile tilapia, Oreochromis niloticus, is native to the Levant and African freshwater systems e.g., in the Western part of the continent (e.g., Senegal, Gambia, Niger, Benue, Chad) as well as to many of the East African rivers (e.g. R. Nile) and Rift Valley Lakes like, Albert, Turkana, George, Edward, Tanganyika, Kivu, etc. [1, 2]. Although O. niloticus is native to Africa, the cichlid is naturally absent in the world’s largest tropical freshwater body, Lake Victoria and the neighboring Lakes Kyoga and Nabugabo as well as many of the East African satellite lakes [2,3,4,5]. These lakes were naturally inhabited by two tilapiine species; O. variabilis (Nyasalapia) and O. esculentus (Ngege) [2,3,4]. For more than nine decades, O. niloticus has been intentionally dispersed worldwide, in particular for aquaculture and restocking programs [2, 6]. In East Africa, various fish introductions are reported, starting in the 1920s. For example, O. niloticus, and other tilapiines e.g. Athi River Tilapia (Tilapia spilurus nigra, Günther 1894) as well as Black bass (Micropterus salmoides), being initially translocated for enhancing fisheries productivity in water bodies naturally considered as unproductive like the southwestern Uganda high-altitude lakes [7, 8]. A case in point is Lake Bunyonyi which was stocked in the 1920s with individuals of O. niloticus from Lake Edward [7]. Similarly, in the 1950s, several tilapiine species were stocked into Lakes Victoria, Nabugabo and Kyoga to counteract the decline of native fish species (O. variabilis and O. esculentus) [3,4,5, 9]. The introduced species; O. niloticus, O. leucosticus (Blue-spotted tilapia), Coptodon zillii (Red-belly tilapia) and O. melanopleura, were all suspected to originate from Lake Albert [3, 4, 9, 10]. However, some introductions might have also originated from Lake Edward and Lake Turkana into Lake Victoria basin [2, 4, 8]. Following these introductions, the indigenous fish species in Lakes Victoria, Kyoga, and Nabugabo, significantly declined in the 1980s, coinciding with the dramatic increase in the stock size of the non-native O. niloticus [3, 4, 11]. The potential reasons for the declined native fish species (O. variabilis and O. esculentus) were suspected to a combination of factors including; competition, over fishing, as well as predation pressures from another introduced species, the Nile perch (Lates niloticus) [3, 5]. But one additional factor might have been hybridization between the native cichlids (O. variabilis and O. esculentus) and the introduced tilapiines, particularly O. niloticus [4, 9, 11]. Based on these events, the expanded distribution of O. niloticus in East Africa complicates the differentiation and identification of genetic units for management and conservation. For example, the population considered as non-native O. niloticus in Lakes Victoria and Kyoga might have genetically diverged via admixture and hybridization with the indigenous species [4, 9, 12]. The loss of indigenous O. mossambicus due to hybridization with the introduced O. niloticus has been reported in South Africa [13]. The situation in East Africa may have worsened with the recent boom of fish hatcheries and aquaculture production systems [14]. In this context, feral populations resulting from escapees might be an additional and serious threat to natural systems.

In nearly the last two decades, the East African countries have been developing measures for fisheries sustainable exploitation through the implementation of co-management strategies [15]. Nonetheless, conservation and management of the already admixed species might not be achieved if the genetic structure of the species in question is not well understood, as the stocks are difficult to define [16]. Therefore, with respect to the East African O. niloticus, as the species were potentially affected by various anthropogenic activities, a thorough characterization of the populations at the molecular level might be needed.

Based on the earliest studies, East African O. niloticus diversity has been studied using both traditional morphometric methods and molecular markers, which led to contradictory patterns in the species description. For example, using biometrics and counts, seven O. niloticus subspecies from African different regions or lakes were described [2]. Nevertheless, this classification was contradicted by subsequent studies using morphometric analyses accompanied by allozyme markers, which indicated that the O. niloticus strain from Lake Edward is closely related to that of the lower Nile (Egypt) [17]. Also, other investigations using restriction endonuclease mitochondrial DNA found that O. niloticus from Lake Tana is distinct, contrary to the earlier traditional morphometric and meristic analyses [18]. Other earlier molecular genetic studies employing allozymes and restriction fragment length polymorphism (RFLP) of mitochondrial DNA (mtDNA), as well as randomly amplified polymorphic DNA (RAPD) for investigating the demography of O. niloticus populations in East Africa, shed some additional light to these incongruences [1, 19, 20]. Some of these studies reported that O. niloticus populations from Albert Nile (the Egyptian stretch of River Nile) are distinct from the West African populations, also contrary to earlier morphometric studies [1]. Furthermore, these past investigations based on traditional markers indicated conflicting results amongst. For example, findings from a combination of allozymes and restriction fragment length polymorphism (RFLP) of mtDNA indicated that O. niloticus in from Lake Tana is clustered with Lake Edward and the Kenyan Lake Turkana system, which differs from the findings based on restriction endonuclease analysis of mtDNA [1, 18]. These results are inconsistent probably because of the different methodological approaches used that comprise different information content [21,22,23,24,25,26]. Additionally, the markers used so far have low resolving power to characterize variation within and between populations, and the genetic fingerprinting markers like RAPD cannot discern between homozygotes and heterozygotes [22]. The lack of methodologies with high discriminating power in the past studies, therefore, suggests that the genetic structure patterns of the East African O. niloticus are insufficiently documented.

In the present study, we utilize nuclear microsatellite markers, simple sequence repeats (SSRs), to typify the O. niloticus in East Africa using next-generation sequencing. SSR loci have been proven robust when investigating the genetic structure of O. niloticus, particularly, using SSR genotyping by sequencing (SSR-GBS) [27]. SSR-GBS approaches are useful because they reduce size homoplasy, which is one of the constraints of traditional SSR fragment length analysis [28, 29]. However, SSR-GBS is not without drawbacks [30]. For example, the presence of stutter complicates allele calling for di-nucleotides, null alleles due to mutation on primer binding sites, and it does not recover genomic information hence overestimating events that had a small impact on the gene-pool. Although generally, the use of SSR fragment length analysis can yield information for delineating populations, the recent studies in East Africa that have used this approach on O. niloticus were limited to few water bodies in Kenya, with the broader scope of the African Great Lakes missing [31,32,33]. It is important to conduct a comparative study of various water bodies where O. niloticus is present (native and non-native with possible admixture). Such research would provide information on the genetic structure and diversity which would establish a firm base for management and conservation of these cichlids [34].

Here, we explicitly investigate the genetic structure of O. niloticus, in East Africa including some populations from Ethiopia and West Africa (Burkina Faso), representing the Sub-Saharan African Great Lakes. We compare natural/native with introduced/non-native O. niloticus populations, including other populations from aquaculture systems. With this approach, we investigate the impact of anthropogenic activities, particularly the translocations, on the O. niloticus’ gene pool. This is especially important to evaluate the genetic integrity of native stocks. We hypothesized that anthropogenic activities have affected the genetic divergence of O. niloticus populations, particularly in environments where the species was introduced. We also predict that the geographical context exhibited by aquatic interconnectivity may influence the genetic homogeneity of cichlid in such environments. We test these hypotheses by answering the following research questions: 1) Does the genetic structure of the East African O. niloticus populations differ from those outside the region? 2) To what extent does the genetic structure of the East African O. niloticus populations reflect the geography and anthropogenic activities associated with the pathways of the translocation?

Results

Variability of SSR loci

In total, 13,530,228 paired reads were produced for genotyping, from which 9,579,578 passed the quality control steps, which were later used for allele calling. Genetic variation results for the 40 SSR loci are presented in the Additional file 1: Table S2. The number of alleles per locus had a mean value of 33.8 ± 20.5, ranging from seven to 84, with a total of 1352 alleles generated across all loci. Overall, 80% of the loci exhibited expected heterozygosity (He) values greater than 0.5. Polymorphic Information Content (PIC) was generally congruent with He, with 78% of loci indicating values of greater than 0.5 (Additional file 1: Table S2).

Genetic structure

The UPGMA dendrogram showed that all East African populations were more similar to each other than to the other regions (Fig. 2). In this case, the three Ethiopian populations (Hashenge, Ziway, and Chamo) formed the most distant group followed by Burkina Faso and the other Ethiopian water body, Lake Tana. Among the East African natives, the largest separation was between the Kenyan, Lake Turkana, and the Ugandan water bodies. In Uganda, with exception of Lake Victoria, the non-native lakes and fish farms grouped with a native population: the southern Ugandan high-altitude Lakes (Kayumbu and Mulehe) with a group comprised by Lakes George, Edward and Kazinga Chanel; Lake Kyoga populations- and Sindi Farm with River Nile; and Bagena and Rwitabingi farms with Albert. Four subpopulations of Lake Victoria (Gaba, Masese, Kakyanga, Kamuwunga) formed a sister group to the River Nile one. The Lake Victoria subpopulation Sango Bay showed the highest degree of divergence in Uganda.

Neighbor network results showed a similar pattern to the UPGMA dendrogram both at regional and local levels (Fig. 3). In this case, however, Burkina Faso was observed to be closer to the Ugandan populations. In general, network results reflected two Ugandan catchment groups: the George, Kazinga Channel, and Edward group together with the non-native Ugandan highland lakes, and on the other end, Albert and River Nile systems together with the non-native Lake Kyoga and all fish farms. Interestingly Lake Victoria exhibited an intermediate position between both groups with the subpopulation from Sango Bay showing a long branch, suggesting high genetic differentiation. Overall, most of the non-native populations (including farms) showed longer branches than the natives (Fig. 3).

Genetic distance between individuals which was visualized through principal coordinates analysis (PCoA), analysis showed a separation of population groups based on geographic regions (Fig. 4a). Samples formed four groups when analyzed at the regional/country level (Fig. 4a): two groups with individuals from Ethiopia, one with individuals from East Africa, and another intermediate group with samples from both regions. The composition of these groups was clearer when the distance between the native individuals was plotted (Fig. 4b). At this level, Lake Turkana clustered with Burkina Faso, and a division between the three Ethiopian Lakes (Hashenge, Chamo, and Ziway) and Lake Tana was clearly observed. Amongst East African populations, the separation between Lake Turkana and the remaining native populations was evident (Fig. 4b). Individuals found in the Ugandan native populations were divided into two main groups (Fig. 5a). One group was composed of Lake Albert and River Nile individuals while the other by Lake Edward, Kazinga Channel, and Lake George. This division was less evident when individuals from non-native and fish farm populations were included in the analysis (Fig. 5b). Here, some individuals from Sango Bay formed a separate group from the remaining Ugandan individuals. A further group composed of Lake Hashenge individuals was found when only Ethiopian individuals were plotted (Fig. 5c). Substructure within the same lake was only evident for Lakes Victoria and Kyoga (Fig. 6).

The Bayesian analysis with STRUCTURE was portrayed based on the optimal K values. For all populations, the best K was 10, all native populations, K = 7, East African native populations, K = 2, Ugandan native populations, K = 2, and all Ugandan populations including farms, K = 4 (Additional file 1: Figure S2). O. niloticus populations from each African region were assigned to different groups (Fig. 7a). Within each region, the same assignments were observed with Lakes Tana and Turkana isolated from the rest of Ethiopians and East African populations, respectively (Fig. 7a). Among the Ugandan native populations, clustering was also congruent with the two water systems, as indicated earlier by both network and PCoA analyses, see Fig. 7b and c. However, there were cases where the non-native populations showed independent clusters from the native. For example, in all analyses, Lake Victoria clusters differed from other populations even when only Ugandan O. niloticus were included in the analysis (Fig. 7c). Apparently, admixture was more evident amongst the East African populations but mostly detected when only non-native populations were considered (Fig. 7c).

Gene flow between population

Results from recent migration rates estimated with BayesAss indicated that Lakes Kyoga and George were the main sources of migration (Fig. 8), with values for other populations generally falling below (< 2%). Noticeable gene flow was from Lakes Kyoga to Victoria and George to Edward (27%), Kyoga to Albert (25%), Kyoga to Bagena farm (23%), Kyoga to Sindi farm, River Nile and Rwitabingi farm (22%), George to Kazinga Channel (21%) and finally George to Mulehe (20.4%) (Fig. 8). Migration rates estimated through Genalex were congruent with BayesAss, but with the difference that the O. niloticus population from Lake Victoria was also a source of migrants (Additional file 1: Table S3).

Genetic differentiation, diversity, and isolation by distance

Genetic differentiation of O. niloticus was consistent with the STRUCTURE results. For instance, the Fst values clearly demonstrated that the East African O. niloticus populations are genetically distant from the Ethiopian and West African populations (Fig. 9a). Despite O. niloticus populations from River Nile and Lake Kyoga showing relatively high Fst values, results from the East African populations generally showed low genetic differentiation. Also, the East African O. niloticus populations were genetically more diverse when compared to either Ethiopian or Burkina Faso (Fig. 9b-d). Based on all statistics, the non-native Lake Victoria and native Lake Turkana O. niloticus populations were the most genetically diverse. On the other hand, Lake Kyoga and River Nile O. niloticus populations were consistently the least diverse even when investigated at the subpopulation level (additional file 1: Figure S4).

Results from the Garza-Williamson index (G-W), generally indicated that nearly all of the studied populations went through a bottleneck, apart from the Ethiopian Lake Tana (Fig. 10a). In the analysis, only Lake Tana exhibited G-W values > 0.5 (0.56 ± 0.44). Regarding population genetic diversity, however, Lakes Victoria and Turkana showed the highest number of private alleles (Fig. 10b).

When we partitioned Lake Victoria to assess the genetic diversity patterns within the water body, generally one sub-population was distinguished from the others (Fig. 11). Sango Bay, in particular, was isolated based on Fst values, and consistently exhibited higher genetic diversity indices (Na, He and Ar) (Fig. 11).

Mantel tests for isolation by distance (IBD) across all samples showed a positive correlation between geographical and genetic distance (R2 = 0.30) (Fig. 12a). However, the strong correlation (R2 = 0.67) between the populations was only found when Burkina Faso was excluded from the analysis (Fig. 12b). The genetic differentiation between the East African and the Ethiopian populations appears to inflate this correlation. Similarly, a strong IBD was also found amongst East African populations (Fig. 12c), which was not the case when only Ugandan populations (excluding Turkana) were considered (Fig. 12d).

Discussion

Fisheries and fishery products are vital in the developing world but heavily threatened through various anthropogenic activities which may compromise the continuity of the resources [35]. One aspect of the anthropogenic threats is the change or alteration of the natural genetic structure of fish stocks through admixture [36, 37]. Understanding the admixture of stocks is only possible if the source populations can be differentiated using genetic markers. We show the importance of SSR-GBS for a deeper understanding of population dynamics, in particular, the East African O. niloticus, towards the alignment of management and conservation strategies. In this study, we investigated the phylogeographical patterns and we found large differences between lakes (e. g. Lake Tana) and also differences between natural water catchments that allow populations to be identified. Here, we discuss the current state of O. niloticus in reference to phylogeographical patterns and anthropogenic activities.

Phylogeography of east African O. niloticus

In all analyses, we found a clear differentiation among all three African regions included in this study (East Africa, Burkina Faso, and Ethiopia), indicating a low degree of connectivity amidst them and highlighting the high level of differentiation between regions. Lake Tana was completely distinct from the remaining populations. This applies not only to the Ethiopian populations but also to the East African ones. So, the genetic distance in Ethiopia is higher than between the East African and West Africa populations, indicating a divergence higher than we would expect within a species. These results are consistent with previous genetic reports [18], but not the findings of the subspecies treatment based on the traditional morphometric and meristics [2]. This high level of differentiation argues for a revision of the species delimitation for these populations.

Lake Tana lies in the Ethiopian mountains and is isolated from the Lakes in the Rift valley [38]. This might explain the high degree of differentiation of this lake because of the lack of connectivity and divergent ecological conditions. Contrary, Lake Hashenge which is also in the Ethiopian mountains is related to the Rift Valley lakes. Lake Hashenge is reported to have been stocked with O. niloticus following mass mortalities of the native species [39]. The native status of this lake is unclear since it could have been restocked with O. niloticus that originated from the Rift Valley Lakes. Besides that, we see a slight differentiation in PCoA between Lake Hashenge and the Rift Valley Lakes in Ethiopia, which may reflect an unsampled source of stocking or differentiation accumulated because of the high degree of isolation of the lake.

In East Africa, genetic structure reflected different catchments. The population from Lake Turkana was genetically distinct from the Ugandan populations which is expected given its high geographical isolation [40]. Our findings concur with the previous works that treated the Turkana population as a different subspecies (O. vulcani) [2]. The high diversity and number of private alleles found in Lake Turkana can be a consequence of this isolation. The East African arid, Lake Turkana, naturally is also characterized by a remarkable genetic diversity. One factor might be introgression perhaps from anthropogenic activities or influx of gene flow from River Omo (Ethiopia). However, this is not clear and a better sampling from the region needs to be included to evaluate the extent of the observed current genetic structure of the population.

In Uganda, despite the high degree of connectivity and proximity between the water bodies, O. niloticus populations were clearly structured. These reflected three main groups: 1) (Lakes George and Edward, as well as Kazinga Channel, 2) Lake Albert, River Nile, and Kyoga and 3) Lake Victoria system. The 2nd and 3rd groups are discussed in more detail under anthropogenic activities subsection. The 1st group, Lakes George and Edward are connected via the Kazinga channel which also explains the high natural migration rates between these populations. The different genetic structure between the western Rift Valley Lakes (Edward-George-Kazinga Channel and Albert) was conserved despite being connected through River Semliki that flows from Lake Edward and Albert [41]. The strong rapids and falls present in this river [41, 42], might constitute a strong barrier to gene flow, which maintains these systems apart. These findings are congruent with recent work on O. niloticus geometric morphometrics [43] but do not concur with past studies [2, 20]. This incongruity might be associated with different methodological approaches utilized between the earliest studies and the current one. For example, using morphometric and meristics methods, O. niloticus from the Edward-George system and Albert was treated as one subspecies; O. niloticus eduardianus [2]. However, inference from traditional morphometrics are weak due to the lack of informative characters [18]. Similarly, while we used SSR-GBS techniques, [20] employed random amplified polymorphic DNA (RAPD) markers, which due to their dominance genotypic nature, provide only part of the information content [22].

Anthropogenic activities-fish translocations

In East Africa, we know that O. niloticus was introduced into several water bodies through stocking activities. We were able to genetically track these translocation events to both non-native water bodies and fish farms. All genetic structure analyses and migration rates showed that the two Ugandan groups (the George-Edward complex and Lake Albert) contributed to the stocking of different water bodies. O. niloticus from the southwestern Ugandan high-altitude Lakes; Mulehe and Kayumbu, originated from the Western Rift Valley Lakes – Edward and George. For the 2nd group, Lake Kyoga and River Nile (Victoria Nile) are genetically similar to Lake Albert, suggesting that, the latter population might have contributed genes to the gene-pool of the former systems. Although Lake Kyoga is connected to Lake Albert via River Nile, their genetic similarity is unlikely related to the consequence of natural migration via water flow. The main reason here is the natural occurrence of Murchison Falls on the River Nile that acts as a barrier between the systems [3, 41]. For this matter, the genetic similarity between River Nile, Lakes Kyoga, and Albert populations may have resulted in stocking regimes using the latter as source [3].

Fish farms seem to have sourced fish seed from multiple populations, resulting in admixed stocks. Our results show that Lakes Albert, and Kyoga, as well as River Nile, contributed to the gene pool of the farmed populations (Figs. 3, 8 & 7c). Based on genetic distance, Lake Albert was the main contributor to Rwitabingi and Bagena farms while Kyoga to Sindi farm. However, we also observed a high amount of gene flow from Kyoga to Rwitabingi and all these farms appeared to be admixed with other populations including Lake Victoria. Apart from farms, evidence of admixture was probable in the East African natural populations, which seems to have been promoted by anthropogenic activities [3, 9]. This is supported by the fact that when non-native populations were unconsidered in the STRUCTURE and PCoA analyses, signals of admixture were minimal, and clear genetic structure assignments could be observed. In East African, admixture in O. niloticus populations may stem from three main processes: 1) translocation from multiple sources into the non-native water bodies, 2) back translocation from non-native to native populations, and 3) hybridization of O. niloticus with congeneric species promoted by translocations.

The first and third processes may explain partly the genetic variation found in the 3rd group; Lake Victoria (see above the three Ugandan groups). Although O. niloticus in Lake Victoria is generally isolated, based on the distance neighbor Network tree (Fig. 3), the population occupied an intermediate position between the above described; 1st and 2nd, Ugandan groups. Thus, it is clearly possible that multiple stockings might have contributed to the gene-pool indicated by the Lake Victoria population. For example, [2] suggests that introductions into Lake Victoria may have originated from Lake Edward, with other authors suggesting multiple sources [4, 5, 12, 44], which support our results. The highly diverse and differentiated gene-pool in Lake Victoria could have originated from the admixture of several lineages due to multiple sources.

On the other hand, possible hybridization of the introduced O. niloticus with the indigenous relative species (O. variabilis and O. esculentus) in Lake Victoria may explain some of the genetic variation patterns found in this lake. First, this lake together with Turkana showed values of private alleles up to four times higher than the remain populations. This genetic variation could have originated from introgression by species that have not been included in the analysis. Similarly, the probable hybridization may explain the high genetic diversity and divergent gene-pool detected in the system. Within Lake Victoria, the Sango Bay subpopulation appears to be an extreme case from this by showing the highest degree of genetic divergence. Remarkable genetic differentiation in Sango Bay was noticed only when compared with the remaining subpopulations within the lake, but also with the other East African populations. In this case, during the boom of the O. niloticus population in Lake Victoria [3,4,5, 45], a larger portion of the native species’ genetic materials may have been introduced into O. niloticus gene-pool. This is just a hypothesis since, in this study, we cannot directly test for hybridization because we did not include samples of O. niloticus congenerics. However, hybridization involving O. niloticus and other tilapiines has been reported to be relatively frequent and it needs to be considered [9, 33, 46, 47].

If admixture/hybridization shaped the gene-pool of Lake Victoria, it may have adaptive consequences and compromise the sustainability of O. niloticus. Although hybridization may lead to heterosis/hybrid vigor [48, 49], admixture is usually reported to have negative consequences [37, 50]. Introgression can contribute to outbreeding depression either by the introduction of maladaptive alleles or through the dilution of alleles important for local adaptation [51]. In more drastic scenarios, hybridization can result in genomic incompatibilities contributing to a fast reduction of population fitness [51]. Alternatively, the hybrids may potentially exhibit more fitness and subsequently extirpate the parental lines [46]. The observed genetic structure of O. niloticus populations in Lake Victoria was unexpected and has not been reported before, which calls for further investigations for taxonomic recognition.

Evidence for the second process of admixture was only found in Lake Albert. In the structure analysis, this population showed admixture with Lake Kyoga. We also found significant migrations from Lake Kyoga to Lake Albert. These results indicated that admixture with respect to translocations not only contributes to non-native populations but also to native ones. The sequence of gene flow from Lake Kyoga to Albert is not clear as none of the previous reports have indicated this. However, it is likely that aquaculture activities might be contributing to the observed gene flow between Lakes Kyoga and Albert.

Anthropogenic activities-consequences of overfishing

Some water bodies, especially Lake Kyoga and River Nile showed low genetic variability and evidence of bottleneck with respect to G-W estimations. Given the recent stocking of these water bodies, this pattern may be explained by the founder effects. Nevertheless other anthropogenic activities need to be considered as well. High loss of genetic diversity among populations, particularly, in fishes has been attributed to over-exploitation [52]. This might be the case for the L. Kyoga population. For example, although O. niloticus boosted the capture fisheries in the Lake Victoria basin (Lakes Victoria and Kyoga) following introductions, the species was subsequently overexploited between the 1970s and 80s [53, 54]. This was reflected in the dramatic decline of the stock sizes and increased fecundity, which are clear indicators of overfishing [54]. The low diversity in River Nile could be linked to low gene-flow connectivity with other water bodies due to hydro-electric power dams that have been constructed along the river (the upper Nile of the Ugandan side), which increases the effect of genetic drift. However, this needs to be assessed in further analyses, especially when additional samples are collected in sections of the lower Nile (below Murchison falls), where apparently there are no dams.

Implications for management and outlook

Overall, we found evidence that anthropogenic activities affected the gene-pool of the East African O. niloticus. The main consequence might have been admixture and potentially hybridization between different stocks and species respectively. In the long term, this may have negative effects on population fitness due to outbreeding depression and genetic swamping. Thus, management measures should inhibit any form of unauthorized spread of fish in the aquatic ecosystems. The Western or Albertine Rift Valley lakes (Edward-George) may be ideal broodstock sources for subsequent breeding programs and aquaculture, as these systems seem not to be admixed. To avoid an influx of feral populations, a proper environmental impact assessment should be prioritized before implementation. Genetic diversity might also have been affected by overfishing and the construction of hydropower dams, which should also be taken into consideration in future management options.

Conclusions

Our results were congruent with the hypothesis that anthropogenic activities affected the genetic structure of O. niloticus populations in East Africa. The genetic variation of some populations, especially from Lake Victoria, corresponded with possible hybridization of O. niloticus with native congeneric species, which may have been mediated by anthropogenic activities. This study also contributed to the knowledge of O. niloticus phylogeography in East Africa. In this case, we found several new genetic groups such as the populations from Lake Tana, Victoria and the two natural catchments in Uganda. Some of these may require further taxonomic exploration. Additionally, we show that gene-flow among the East African O. niloticus populations was not entirely from native to non-native environments, but also from non-native to native environments likely through aquaculture and restocking programs. Moreover, this study shows the importance of molecular markers, in particular, the use of SSR-GBS in cataloging populations. Further studies should include O. niloticus samples from other regions such as the lower Nile (below Murchison Falls), Lake Kivu (Rwanda), Tanganyika and Baringo as well as the congenerics for a more comprehensive picture.

Methods

Sampling/study areas

We collected O. niloticus specimens from three water body types: a) those where O. niloticus is native, b) where introduced, and c) from fish farms (Fig. 1), following our earlier sampling design [43]. Most samples were collected by local fishermen using gill nets set overnight. At Lake Turkana, a seine net was utilized. From Ethiopia and Burkina Faso, four and one native populations were sampled, respectively. Considering the large extent of Lake Victoria and multiple O. niloticus introductions into the world’s largest tropical freshwater body, we sampled five locations to assess possible genetic heterogeneity within the system (Fig. 1). Similarly, in other relatively large lakes like Lake Edward, Kyoga, and Albert, we sampled two locations each for subsequent subpopulation analyses (Table 1). A total of 664 samples were collected from 18 water bodies during several field excursions in 2016. From every single fish, a muscle tissue sample (approx. 30 mg) was extracted from the dorsal region, preserved in absolute ethanol contained in 2 ml Eppendorf tubes and later stored in a freezer until genotyping at the Institute for Integrative Nature Conservation Research-University of Natural Resources and Applied Life Sciences Vienna (BOKU), Austria. Sampling was conducted in collaboration with respective authorities per region and therefore no special permission was required. In all cases, the fish were already dead when obtained from the fishermen, therefore no special treatment for the animals was administered in the process. As contamination of the specimens was not likely during sampling with gill nets, great care and attention were provided for during seining on Lake Turkana. The non-native and farm populations were only sampled in Uganda. Here, we refer to the non-native populations like those found in the high-altitude satellite lakes of south-western Uganda (Lakes Mulehe and Kayumbu) as well as in lower altitude lakes (Lake Victoria and Kyoga) [43]. The three sampled fish farms include; Rwitabingi (located near River Nile and Lake Kyoga), Bagena and Sindi from South-western Uganda. The rest of the populations are regarded as native (Fig. 1; Table 1).

Fig. 1
figure1

Illustration of sample collection and sources in the African Great Lakes region; East Africa (Uganda and Kenya), Ethiopia and Burkina Faso; modified from our previous work [43]

Table 1 Details of the sampling sites and the total number of individuals collected per water body and location/site. The indigenous O. niloticus populations, are also herein referred to as natives and introduced, non-natives and farms are the pond culture systems
Fig. 2
figure2

Genetic structure of O. niloticus populations based on UPGMA dendrogram. Node values correspond to bootstrap values

Fig. 3
figure3

Genetic structure based on unrooted network tree illustrating population relationships based on genetic distance. a represents a network tree for all the populations and b for only the Ugandan populations. Dotted oval and rectangular shapes depict closely related genetic groups

Fig. 4
figure4

Genetic scatter plots of O. niloticus exhibited by Principal Coordinate Analysis (PCoA). a populations per region, b all indigenous populations. PCoA was constructed with respect to unbiased Nei’s genetic distance among individuals

Fig. 5
figure5

Genetic scatter plots of O. niloticus based-on Principal Coordinate Analysis (PCoA). a Ugandan native populations, b all Ugandan populations including non-natives, natives, and farms, and c all Ethiopian populations

Fig. 6
figure6

Genetic scatter plots of O. niloticus exhibited by PCoA within Lakes Victoria (a) and Kyoga (b) populations

Fig. 7
figure7

Bayesian clustering for genetic assignments of O. niloticus populations. a represents all populations, b all indigenous populations, and c all Ugandan populations including indigenous, non-indigenous and farms. Ks with a superscript symbol (¶) indicates the optimal K values based on STRUCTURE HARVESTER analyses

Fig. 8
figure8

Bayesian inference of recent migratory rates for the 13 East African O. niloticus populations. Oval light blue and rectangular light-blue shapes indicate natural and farm populations, respectively. The arrows contain percentage values showing the direction and magnitude of gene flow. Darker and thick arrows represent stronger gene flow, while thin, dotted arrows indicate weaker gene flow. Native and non-native populations are indicated by the letters, “N” and “NN”, respectively. This analysis is based on BayesAss program and for GenAlex program, see the Additional file 1: Table S3

Fig. 9
figure9

Genetic diversity and differentiation indices. a number of alleles, b allelic richness, c fixation index (Fst) and d expected heterozygosity

Fig. 10
figure10

Estimations of population bottleneck derived from Garza-Williamson Index (G-W) (a) and measure of genetic diversity based on private alleles (b)

Fig. 11
figure11

Genetic diversity of Lake Victoria within the population. a number of alleles, b allelic richness c expected heterozygosity and d private alleles

Fig. 12
figure12

Mantel tests for correlations between genetic distance (Fst) and Euclidean geographical distance (GGD in Km) for O. niloticus populations. a represents isolation by distance (IBD) between all populations, b all populations without Burkina Faso, c East African, and d only Ugandan populations

Genotyping

Genomic DNA extraction was conducted using magnetic beads based on the MagSi-DNA Vegetal kit (MagnaMedics, Geleen, Netherlands) and a magnetic separator, SL-MagSep96 (Steinbrenner, Germany) [27, 30]. We used microsatellite markers [27], to which we added 15 extra primers (Table 2, see also Additional file 1: Table S3). The SSR primers were designed and tested following our earlier work [27], using the same shotgun sequencing data present in the sequence read archive database (SRA) under the reference number SRX3398501. Screened primers were then grouped into three multiplexes and used to prepare amplicon SSR-GBS libraries using the same approach and specifications of [27]. The PCR products were then pooled and sent for paired-end 300 bp sequencing in Illumina MiSeq, at the Genomics Service Unit in Ludwig Maximillian Universität, München, Germany. The raw sequence data were deposited in the GenBank, sequence read archive database (SRA) under the project PRJNA550300 with the accession numbers, SRR9587388 to SRR9587270. Sequences generated by Illumina, were subsequently quality checked and controlled, which were later used for alleles calling as described in [27, 30] using the scripts from the SSR-GBS pipeline (https://github.com/mcurto/SSR-GBS-pipeline). The resulting codominant matrix and information for which sequences correspond to each allele can be found in the Additional file 2 (see the file named “Second_additional fileAllelesList & matrix_”).For subsequent analyses, all loci and samples with missing genotypes ≥50% were excluded, leaving a total number of 40 markers (Additional file 1: Tables S1, S3). Other studies have indicated that many SSR loci are not necessary in order to detect population structure [55, 56], so we did not see the need of developing additional markers to the 40 already in use.

Table 2 15 new primer pairs developed in the present study. The other 26 tested primers developed by [27] can be found in the additional file section, Additional file 1: Table S1

Genetic structure

Genetic structure was first assessed by calculating the genetic distance between individuals and then visualized through Principal Coordinate Analysis (PCoA), all conducted in GenAlex Version 6.5 [57]. Genetic similarity between populations was evaluated by plotting a Neighbor-Net tree based on Nei’s genetic distance [58], using the program, SplitsTree4 version, 4.14.8 [59]. We also constructed UPGMA dendrograms for making inferences on the hierarchical clustering using Nei’s genetic distance as implemented in Populations-1.2.32 [60]. Support values were estimated with 1000 bootstrap replicates based on loci resampling. Neighbor-Net tree and the UPGMA dendrogram were conducted with the inclusion of subpopulations, when applicable to evaluate possible substructure within the populations. Genetic structure was further investigated using the program, STRUCTURE Version 2.3.4 [61]. STRUCTURE clusters individuals into hypothetical populations through optimization of Hardy-Weinberg equilibrium [62]. STRUCTURE was run from K = 1–35 for 10,000 Markov chain Monte Carlo (MCMC) generations after a burn-in length of 10,000 generations [63], whereby each run was iterated 20 times. The program’s default settings for the admixture model and allele frequencies correlated were implemented. Detection of optimal K was done with STRUCTURE HARVESTER [64] using the delta K (ΔK) statistic, which is the second-order rate of change (InP(D)) across successive K values [63, 65]. In this context, STRUCTURE HARVESTER uses ΔK to identify the highest value and henceforth the best K. Results from multiple replicates were summarized using the online pipeline Clumpak program [66] available at http://clumpak.tau.ac.il/. Similar analyses were performed for Lake Victoria within populations.

Migration rates and number of migrants per generation (nm)

Recent migratory rates and the number of migrants per generation were determined as proxy estimates of gene flow among the O. niloticus populations. However, recent migratory rates were only estimated for the East African populations, since the corresponding water bodies are the most affected by anthropogenic activities such as fish translocations. Pairwise recent migration rates were estimated using BayesAss Version 3.0 [67]. Here, the program was run for 200, 000,000 iterations, discarding the first 100,000,000 generations and sampling every 1000th generation [68]. Only results with a 95% confidence interval of a fraction of migrants per population above 0.01 were considered significant. Recent migration rates were used because most of the fish translocations in the region, seemingly were recent. Additionally, we estimated the number of migrants (Nm) per generation between population pairs, to validate the recent migration rates using GenAlex program. Consequently, we present both, the percentage of migrants estimated in BayesAss and the number of migrants between population pairs against the fixation index (Fst) values.

Genetic diversity, differentiation, and isolation by distance (IBD)

Genetic diversity and differentiation indices between O. niloticus populations throughout East Africa and beyond were examined using the following indices: expected heterozygosity (He), observed heterozygosity (Ho), number of alleles (Na), allelic richness (Ar), fixation index (Fst), private alleles, and Garza-Williamson index (G-W). Na, Fst, G-W and He per population were analyzed using the program Arlequin Version 3.5 [69]. Ho, He, Na and PIC per locus were determined through Cervus version 3.0.7 [70]. Ar was analyzed using the rarefaction algorithm implemented in the Hp-rare program [71]. G-W was used to explore the possibility of bottlenecks amongst the populations. If G-W values are closer to zero, it implies that the populations went through a bottleneck, but when the values are close to one, the populations are in a stable phase [72]. To test whether the genetic diversity and differentiation of O. niloticus populations conform to isolation by distance (IBD), we plotted genetic distance (Fst) against the geographical distance (GGD in kilometers) and conducted correlation analyses using Mantel test (999 permutations) implemented in GenAlex Version 6.5 [57].

Availability of data and materials

Raw sequence data were submitted to the sequence read archive (SRA) database and can be accessed under the reference number PRJNA550300.

Abbreviations

Ar:

Allelic richness

G-W:

Garza-Williamson index

GenAlex:

Genetic Analysis in excel

He:

expected heterozygosity

Ho:

observed heterozygosity

IBD:

Isolation By Distance

MCMC:

Markov chain Monte Carlo

mtDNA:

mitochondrial DNA

Na:

Number of alleles

PCoA:

Principal Coordinates Analysis

PIC:

Polymorphic Information Content

RAPD:

Randomly Amplified Polymorphic DNA

RFLP:

Restriction FragmentLlength Polymorphism

SRA:

Sequence Read Archive

SSR-GBS:

Simple Sequence Repeat- Genotyping By Sequencing

SSRs:

Simple Sequence Repeats

UPGMA:

Unweighted Pair Group Method with Arithmetic Mean

References

  1. 1.

    Agnèse J-F, 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(1):88–96.

    PubMed  Article  Google Scholar 

  2. 2.

    Trewavas E. Tilapiine fishes of the genera Sarotherodon, Oreochromis and Danakilia. London: British Museum (Natural History); 1983.

  3. 3.

    Ogutu-Ohwayo R. The reduction in fish species diversity in lakes Victoria and Kyoga (East Africa) following human exploitation and the introduction of non-native fishes. J Fish Biol. 1990;37:207–8.

    Article  Google Scholar 

  4. 4.

    Balirwa JS. The evolution of the fishery of Oreochromis niloticus (Pisces: Cichlidae) in Lake Victoria. Hydrobiologia. 1992;232(1):85–9.

    Article  Google Scholar 

  5. 5.

    Kaufman L. Catastrophic change in species-rich freshwater ecosystems. BioScience. 1992;42(11):846–58.

    Article  Google Scholar 

  6. 6.

    Canonico GC, Arthington A, McCrary JK, Thieme ML. The effects of introduced tilapias on native biodiversity. Aquat Conserv Mar Freshwat Ecosyst. 2005;15(5):463–83.

    Article  Google Scholar 

  7. 7.

    Worthington EB. A report on the fisheries of Uganda investigated by the Cambridge expedition to the east African Lakes, 1930–31; 1932.

    Google Scholar 

  8. 8.

    Agnèse J-F, Adépo-Gourène B, Pouyaud L. Natural hybridization in tilapias. In: Agnèse Jean-François (ed.). Genetics and aquaculture in Africa. Paris: ORSTOM, 95-103. (Colloques et Séminaires). Génétique et Aquaculture en Afrique: Colloque, Abidjan (CIV); 1997. ISBN 2-7099-1390-9.

  9. 9.

    Welcomme R. Recent changes in the stocks of Tilapia in Lake Victoria. Nature. 1966;212(5057):52.

    Article  Google Scholar 

  10. 10.

    Mboya D, Manyala J, Ngugi C. Fish introductions and their impact on the biodiversity and the fisheries of Lake Victoria. Fisheries Department, Kisumu, Kenya, 11 pp. 2001.

  11. 11.

    Mwanja WW. Genetic biodiversity and evolution of two principal fisheries species groups, the Labeine and Tilapiine, of Lake Victoria region. East Africa: The Ohio State University; 2000.

    Google Scholar 

  12. 12.

    Pullin RS. Tilapia genetic resources for aquaculture. Proceedings of the workshop on Tilapia genetic resources for aquaculture, 23-24 March 1987. The WorldFish Center Working Papers 1988.

  13. 13.

    Zengeya TA, Booth AJ, Chimimba CT. Broad niche overlap between invasive Nile tilapia Oreochromis niloticus and indigenous congenerics in southern Africa: should we be concerned? Entropy. 2015;17(7):4959–73.

    Article  Google Scholar 

  14. 14.

    Rothuis A, Turenhout M, van Duijn A, Roem A, Rurangwa E, Katunzi E, et al. Aquaculture in East Africa: a regional approach. IMARES/LEI; 2014.

  15. 15.

    Nunan F, Hara M, Onyango P. Institutions and co-management in east African inland and Malawi fisheries: a critical perspective. World Dev. 2015;70:203–14.

    Article  Google Scholar 

  16. 16.

    Allendorf FW, Luikart G. Conservation and the genetics of populations. Mammalia. 2007;2007:189–97.

    Google Scholar 

  17. 17.

    Vreven EJ, Adépo-Gourène B, Agnèse J-F, Teugels G. Morphometric and allozyme variation in natural populations and cultured strains of the Nile tilapia Oreochromis niloticus (Teleostei, Cichlidae). Belgian J Zool. 1998;128(1):23–34.

    Google Scholar 

  18. 18.

    Seyoum S, Kornfield I. Identification of the subspecies of Oreochromis niloticus (Pisces: Cichlidae) using restriction endonuclease analysis of mitochondrial DNA. Aquaculture. 1992;102(1–2):29–42.

    CAS  Article  Google Scholar 

  19. 19.

    Mwanja WW, Booton G, Kaufman L, Chandler M, Fuerst P, Donaldson E, et al., editors. Population and stock characterization of Lake Victoria. Tilapine fishes based on RAPD markers. Aquaculture Biotechnology Symposium Prodeedings of the International Congress on the Biology of Fishes (ed, by EM Donaldson & DD Mackinlay); 1996.

  20. 20.

    Mwanja WW, Booton GC, Kaufman L, Fuerst PA. A profile of the introduced Oreochromis niloticus (Pisces: Teleostei) populations in Lake Victoria Region in relation to its putative origin of Lakes Edward and Albert (Uganda-E. Africa) based on random amplified polymorphic DNA analysis. African J Biotechnol. 2008;7(11):1769–73.

  21. 21.

    Miah G, Rafii MY, Ismail MR, Puteh AB, Rahim HA, Islam KN, et al. A review of microsatellite markers and their applications in rice breeding programs to improve blast disease resistance. Int J Mol Sci. 2013;14(11):22499–528.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  22. 22.

    Wirgin II, Waldman JR. What DNA can do for you. Fisheries. 1994;19(7):16–27.

    Article  Google Scholar 

  23. 23.

    Okumuş İ, Çiftci Y. Fish population genetics and molecular markers: II-molecular markers and their applications in fisheries and aquaculture. Turk J Fish Aquat Sci. 2003;3(1):51–79.

    Google Scholar 

  24. 24.

    Guichoux E, Lagache L, Wagner S, Chaumeil P, Léger P, Lepais O, et al. Current trends in microsatellite genotyping. Mol Ecol Resour. 2011;11(4):591–611.

    CAS  PubMed  Article  Google Scholar 

  25. 25.

    Hedrick PW. Conservation genetics: where are we now? Trends Ecol Evol. 2001;16(11):629–36.

    Article  Google Scholar 

  26. 26.

    Liu Z, Cordes J. DNA marker technologies and their applications in aquaculture genetics. Aquaculture. 2004;238(1):1–37.

    CAS  Article  Google Scholar 

  27. 27.

    Tibihika PD, Curto M, Dornstauder-Schrammel E, Winter S, Alemayehu E, Waidbacher H, et al. Application of microsatellite genotyping by sequencing (SSR-GBS) to measure genetic diversity of the east African Oreochromis niloticus. Conserv Genet. 2019:1–16.

  28. 28.

    Narum SR, Buerkle CA, Davey JW, Miller MR, Hohenlohe PA. Genotyping-by-sequencing in ecological and conservation genomics. Mol Ecol. 2013;22(11):2841–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  29. 29.

    Vartia S, Villanueva-Cañas JL, Finarelli J, Farrell ED, Collins PC, Hughes GM, et al. A novel method of microsatellite genotyping-by-sequencing using individual combinatorial barcoding. R Soc Open Sci. 2016;3(1):150565.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  30. 30.

    Curto M, Winter S, Seiter A, Schmid L, Scheicher K, Barthel LM, et al. Application of a SSR-GBS marker system on investigation of European Hedgehog species and their hybrid zone dynamics. Ecol Evol. 2019;9(5):2814–32.

    PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Ndiwa TC, Nyingi DW, Agnese J-F. An important natural genetic resource of Oreochromis niloticus (Linnaeus, 1758) threatened by aquaculture activities in Loboi drainage, Kenya. PloS One. 2014;9(9):1–15.

    Article  CAS  Google Scholar 

  32. 32.

    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.

    CAS  Article  Google Scholar 

  33. 33.

    Angienda PO, Lee HJ, Elmer KR, Abila R, Waindi EN, Meyer A. Genetic structure and gene flow in an endangered native tilapia fish (Oreochromis esculentus) compared to invasive Nile tilapia (Oreochromis niloticus) in Yala swamp. East Afr Conserv Genet. 2011;12(1):243–55.

    Article  Google Scholar 

  34. 34.

    Frankham R, Briscoe DA, Ballou JD. Introduction to conservation genetics: Cambridge university press; 2002.

    Google Scholar 

  35. 35.

    Welcomme R, Cowx IG, Coates D, Béné C, Funge-Smith S, Halls A, et al. Inland capture fisheries. Philosophical Transactions of the Royal Society B: Biological Sciences. 2010;365(1554):2881–96.

    Article  Google Scholar 

  36. 36.

    Todesco M, Pascual MA, Owens GL, Ostevik KL, Moyers BT, Hübner S, et al. Hybridization and extinction. Evol Appl. 2016;9(7):892–908.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  37. 37.

    Rhymer JM, Simberloff D. Extinction by hybridization and introgression. Annu Rev Ecol Syst. 1996;27(1):83–109.

    Article  Google Scholar 

  38. 38.

    Vijverberg J, Dejen E, Getahun A, Nagelkerke LA. The composition of fish communities of nine Ethiopian lakes along a north-south gradient: threats and possible solutions. Anim Biol. 2012;62(3):315–35.

    Article  Google Scholar 

  39. 39.

    Tsegay T. Characteristics and status of the Highland Lake Hashenge fishery, Tigray, northern Ethiopia. Int J Aquaculture Fishery Sci. 2017;3(4):71–6.

    Article  Google Scholar 

  40. 40.

    Odada EO, Olago DO, Bugenyi F, Kulindwa K, Karimumuryango J, West K, et al. Environmental assessment of the east African rift valley lakes. Aquat Sci. 2003;65(3):254–71.

    Article  Google Scholar 

  41. 41.

    Dumont HJ. The Nile: origin, environments, limnology and human use: Springer Science & Business Media; 2009.

  42. 42.

    Ogutu-Ohwayo R. The decline of the native fishes of lakes Victoria and Kyoga (East Africa) and the impact of introduced species, especially the Nile perch, Lates niloticus, and the Nile tilapia, Oreochromis niloticus. Environ Biol Fishes. 1990;27(2):81–96.

    Article  Google Scholar 

  43. 43.

    Tibihika PD, Waidbacher H, Masembe C, Curto M, Sabatino S, Alemayehu E, et al. Anthropogenic impacts on the contextual morphological diversification and adaptation of Nile tilapia (Oreochromis niloticus, L. 1758) in East Africa. Environ Biol Fish. 2018;101(3):363–81.

    Article  Google Scholar 

  44. 44.

    Fuerst P, Mwanja W, Kaufman L, editors. The genetic history of the introduced Nile tilapia of Lake Victoria (Uganda–E. Africa): the population structure of Oreochromis niloticus (Pisces: Cichlidae) revealed by DNA microsatellite markers. Tilapia Aquaculture in the 21st Century Proceedings from the Fifth International Symposium on Tilapia in Aquaculture; 2000.

  45. 45.

    Balirwa JS, Chapman CA, Chapman LJ, Cowx IG, Geheb K, Kaufman L, et al. Biodiversity and fishery sustainability in the Lake Victoria basin: an unexpected marriage? BioScience. 2003;53(8):703–15.

    Article  Google Scholar 

  46. 46.

    Deines A, Bbole I, Katongo C, Feder J, Lodge D. Hybridisation between native Oreochromis species and introduced Nile tilapia O. niloticus in the Kafue River, Zambia. Afr J Aquat Sci. 2014;39(1):23–34.

    Article  Google Scholar 

  47. 47.

    Shechonge A, Ngatunga BP, Tamatamah R, Bradbeer SJ, Harrington J, Ford AG, et al. Losing cichlid fish biodiversity: genetic and morphological homogenization of tilapia following colonization by introduced species. Conserv Genet. 2018;19(5):1199–209.

    PubMed  PubMed Central  Article  Google Scholar 

  48. 48.

    Benvenuto C, Tabone E, Vercken E, Sorbier N, Colombel E, Warot S, et al. Intraspecific variability in the parasitoid wasp Trichogramma chilonis: can we predict the outcome of hybridization? Evol Appl. 2012;5(5):498–510.

    PubMed  PubMed Central  Article  Google Scholar 

  49. 49.

    Dodson JJ, Bourret A, Barrette MF, Turgeon J, Daigle G, Legault M, et al. Intraspecific genetic admixture and the morphological diversification of an estuarine fish population complex. PLoS One. 2015;10(4):e0123172.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  50. 50.

    Allendorf FW, Leary RF, Spruell P, Wenburg JK. The problems with hybrids: setting conservation guidelines. Trends Ecol Evol. 2001;16(11):613–22.

    Article  Google Scholar 

  51. 51.

    Roberts DG, Gray CA, West RJ, Ayre DJ. Marine genetic swamping: hybrids replace an obligately estuarine fish. Mol Ecol. 2010;19(3):508–20.

    PubMed  Article  Google Scholar 

  52. 52.

    Allendorf FW, England PR, Luikart G, Ritchie PA, Ryman N. Genetic effects of harvest on wild animal populations. Trends Ecol Evol. 2008;23(6):327–37.

    PubMed  Article  Google Scholar 

  53. 53.

    Ogutu-Ohwayo R, Odongkara K, Okello W, Mbabazi D, Wandera S, Ndawula L, et al. Variations and changes in habitat, productivity, composition of aquatic biota and fisheries of the Kyoga lake system: lessons for management. Afr J Aquatic Sci. 2013;38(sup1):1–14.

    Article  Google Scholar 

  54. 54.

    Njiru M, Getabu A, Jembe T, Ngugi C, Owili M, Van der Knaap M. Management of the Nile tilapia (Oreochromis niloticus (L.)) fishery in the Kenyan portion of Lake Victoria, in light of changes in its life history and ecology. Lakes Reserv Res Manag. 2008;13(2):117–24.

    Article  Google Scholar 

  55. 55.

    Neophytou C, Torutaeva E, Winter S, Meimberg H, Hasenauer H, Curto M. Analysis of microsatellite loci in tree of heaven (Ailanthus altissima (Mill.) Swingle) using SSR-GBS. Tree Genet Genomes. 2018;14(6):82.

    Article  Google Scholar 

  56. 56.

    Arthofer W, Heussler C, Krapf P, Schlick-Steiner BC, Steiner FM. Identifying the minimum number of microsatellite loci needed to assess population genetic structure: a case study in fly culturing. Taylor Francis. 2018;12(1):13–22.

    Google Scholar 

  57. 57.

    Peakall R, Smouse PE. GENALEX 6: genetic analysis in excel. Population genetic software for teaching and research. Mol Ecol Resour. 2006;6(1):288–95.

    Article  Google Scholar 

  58. 58.

    Nei M. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics. 1978;89(3):583–90.

    CAS  PubMed  PubMed Central  Google Scholar 

  59. 59.

    Huson DH, Bryant D. Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2005;23(2):254–67.

    Article  CAS  Google Scholar 

  60. 60.

    Langella O. Populations 1.2.32: population genetic software (individuals or population distances, phylogenetic trees). http://www.bioinformatics.org/project/?group_id=84. 2011.

  61. 61.

    Hubisz MJ, Falush D, Stephens M, Pritchard JK. Inferring weak population structure with the assistance of sample group information. Mol Ecol Resour. 2009;9(5):1322–32.

    PubMed  PubMed Central  Article  Google Scholar 

  62. 62.

    Porras-Hurtado L, Ruiz Y, Santos C, Phillips C, Carracedo Á, Lareu MV. An overview of STRUCTURE: applications, parameter settings, and supporting software. Front Genet. 2013;4:1–3.

  63. 63.

    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–20.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  64. 64.

    Earl DA. vonHoldt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4(2):359–61.

    Article  Google Scholar 

  65. 65.

    Chaluvadi SR, Khanam S, Aly MA, Bennetzen JL. Genetic diversity and population structure of native and introduced date palm (Phoenix dactylifera) germplasm in the United Arab Emirates. Trop Plant Biol. 2014;7(1):30–41.

    Article  Google Scholar 

  66. 66.

    Kopelman NM, Mayzel J, Jakobsson M, Rosenberg NA, Mayrose I. Clumpak: a program for identifying clustering modes and packaging population structure inferences across K. Mol Ecol Resour. 2015;15(5):1179–91.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  67. 67.

    Wilson GA, Rannala B. Bayesian inference of recent migration rates using multilocus genotypes. Genetics. 2003;163(3):1177–91.

    PubMed  PubMed Central  Google Scholar 

  68. 68.

    Curto M, Puppo P, Kratschmer S, Meimberg H. Genetic diversity and differentiation patterns in Micromeria from the Canary Islands are congruent with multiple colonization dynamics and the establishment of species syngameons. BMC Evol Biol. 2017;17(198):1–16.

    Google Scholar 

  69. 69.

    Excoffier L, Lischer HE. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and windows. Mol Ecol Resour. 2010;10(3):564–7.

    Article  Google Scholar 

  70. 70.

    Kalinowski ST, Taper ML, Marshall TC. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol Ecol. 2007;16(5):1099–106.

    PubMed  Article  Google Scholar 

  71. 71.

    Kalinowski ST. Hp-rare 1.0: a computer program for performing rarefaction on measures of allelic richness. Mol Ecol Resour. 2005;5(1):187–9.

    CAS  Article  Google Scholar 

  72. 72.

    Excoffier L, Estoup A, Cornuet J-M. Bayesian analysis of an admixture model with mutations and arbitrarily linked markers. Genetics. 2005;169(3):1727–38.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

Download references

Acknowledgments

We are indebted to Kachwekano Zonal Agricultural Research and Development Institute, National Agricultural Research Organization, Uganda for their great support during field activities. We are grateful to the local governments, in particular, the respective field Fisheries Officers from East Africa and the Ethiopian Fishery and Aquatic Life Research Centre, for their assistance and advice during field excursions. Special thanks to Ms. Howot Teshome and Dr. Mano Komandan who transported our Ethiopian and Burkina Faso samples respectively to Vienna, Austria.

Funding

This work was funded by the Austrian Partnership Program in Higher Education and Research for Development (APPEAR) through the STRECAFISH project, a program of the Austrian Development Cooperation (ADC) and implemented by the Austrian Agency for International Cooperation in Education and Research (OeAD-GmbH).

Author information

Affiliations

Authors

Contributions

PDT, MC, and HM developed the research concept and performed experimental design, laboratory work, bioinformatics, and data analysis with contributions of EA, CM, and HW. PDT, EA, HW, PA, and CM conducted fieldwork. HM provided laboratory space/analytical tools. PDT and MC led the writing process with substantial contributions of the other authors. All co-authors read and approved the final version of the manuscript.

Corresponding author

Correspondence to Papius Dias Tibihika.

Ethics declarations

Ethics approval and consent to participate

Field excursions were conducted together with respective authorities per region and therefore no special permission was required. Throughout the sampling regimes, fish were always dead when obtained from the local fishermen, thus no special treatment for the animals was administered in the process.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Final List of 40 SSR primers utilized in PCR reactions. Table S2. Variation of the final set of 40 SSR loci. Loci Ti1–Ti35 were taken from Tibihika et al. [28] while loci Ti39–Ti61 were developed in the current study. Table S3. List of the linkage group and location/position for each of the primer pairs used. These primers were BLAST against the genome assembly GCF_001858045.2 available in GenBank. The field reference refers to the accession number of the matching Linkage group. Figure S1. Genetic structure based on PCoA within the East African Western Rift Valley lakes, Albert (a) and Edward (b) metapopulations. Figure S2. STRUCTURE HARVESTER analyses for depicting the optimal K values, which is derived from STRUCTURE results. a represents the best K for all populations (K = 10), b best K for all native populations (K = 7), c best K for East African native populations (K = 2), d best K for only Ugandan native populations (K = 2), e best K for all Ugandan populations including farms (K = 4) and f best K for Lake Victoria within populations (K = 2). Figure S3. Bayesian clustering for genetic assignments of O. niloticus populations. a represents East African native, b Only Uganda native populations and c Lake Victoria subpopulations. K values with a superscript symbol (¶) indicate the optimal K value based on STRUCTURE HARVESTER analyses. Table S3. Matrix for the number of migrants per generation between population pairs as evidenced by Fst values. Bold values (> 10) indicate the vital number of migrants explained by Fst values: Note the variations in boldness which explains the level of important migrations between population pairs. Figure S4. Within-population genetic diversity of Albert, Edward, Kyoga, and Victoria populations. Vic = Victoria, Na = mean number of alleles, No. = number, Pa = mean number of private alleles and He = mean expected heterozygosity.

Additional file 2.

Allele list.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Tibihika, P.D., Curto, M., Alemayehu, E. et al. Molecular genetic diversity and differentiation of Nile tilapia (Oreochromis niloticus, L. 1758) in East African natural and stocked populations. BMC Evol Biol 20, 16 (2020). https://doi.org/10.1186/s12862-020-1583-0

Download citation

Keywords

  • Cichlids
  • Fish translocations
  • Genetic structure
  • Gene flow
  • Bottleneck