Next Article in Journal
The Effect of Grapevine Variety and Wine Region on the Primer Parameters of Wine Based on 1H NMR-Spectroscopy and Machine Learning Methods
Next Article in Special Issue
Genetic Variability within the Murgese Horse Breed Inferred from Genealogical Data and Morphometric Measurements
Previous Article in Journal
Taxonomic Diversity of Decapod and Stomatopod Crustaceans Associated with Pocilloporid Corals in the Central Mexican Pacific
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Defining Management Units for Wild Nile Tilapia Oreochromis niloticus from Nine River Basins in Ghana

by
Gifty Anane-Taabeah Attu
1,2,
Emmanuel A. Frimpong
2,* and
Eric M. Hallerman
2
1
Department of Fisheries and Watershed Management, Kwame Nkrumah University of Science and Technology, Private Mail Bag, University Post Office, Kumasi, Ghana
2
Department of Fish and Wildlife Conservation, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA
*
Author to whom correspondence should be addressed.
Diversity 2022, 14(2), 73; https://doi.org/10.3390/d14020073
Submission received: 19 November 2021 / Revised: 13 January 2022 / Accepted: 14 January 2022 / Published: 21 January 2022

Abstract

:
Despite the global importance of the Nile tilapia Oreochromis niloticus, especially to aquaculture, knowledge of genetic variability within native populations is still limited. While several studies have assessed genetic differentiation across the major drainage basins of Africa, relatively little effort has focused on characterizing genetic differentiation at finer scales. We assessed genetic variation in O. niloticus within and among nine drainage basins in Ghana using nuclear microsatellite DNA markers as the basis for identifying potential units of conservation among wild populations. We screened 312 wild individuals using eight nuclear microsatellite DNA markers. We found moderate genetic diversity within and differentiation among all wild populations studied, with strong signals of recent demographic bottlenecks in several populations. Genetic structure among 11 populations suggested the presence of up to ten management units (MUs). In particular, the Black Volta and the Tano–Asuhyea populations, which were the most genetically distinct and geographically isolated and may be most at risk of loss of genetic diversity over time, may well represent evolutionary significant units. Therefore, at the minimum, the Black Volta and Tano–Asuhyea populations should be prioritized for conservation actions to sustain them over the long-term.

1. Introduction

The Nile tilapia (Oreochromis niloticus) is the third most cultured species worldwide, with a 2018 production of 4.5 million metric tons [1]. Given the importance of O. niloticus to aquaculture and commercial fisheries, particularly in Africa, it is surprising that our knowledge of genetic diversity within and among populations is still rather limited. Most research has focused on characterizing aquaculture performance differences among geographically defined strains [2,3,4], with comparatively little attention to population genetic differentiation in the wild.
Several early studies have attempted to cover multiple major drainage basins in Africa where O. niloticus occurs. Allozyme markers showed genetic differentiation among western and eastern African populations [5,6], suggesting the effects of paleogeographic events. Allozyme and mitochondrial restriction fragment length polymorphisms showed genetic differentiation of western African, Ethiopian Rift Valley, and Nile drainage populations [7]. In addition to the effects of paleogeographic events, microsatellite DNA variation indicated more recent, isolation-driven differentiation between and within drainage basins [8]. Our understanding of differentiation among O. niloticus populations can be advanced by conducting drainage-scale, local and regional studies, which incorporate multiple, highly variable molecular markers [9,10,11]. This is especially useful for identifying populations that may be demographically isolated from other populations (as management units, or MUs; [12]) or populations that are genetically distinct, have unique adaptations, and require conservation management (evolutionarily significant units, or ESUs; [13]).
Frimpong et al. [14] showed the importance of fine-scale population genetic studies within individual drainage basins as opposed to non-systematic sampling among multiple large rivers. The authors screened five microsatellite loci among collections of O. niloticus from three major tributaries of the Volta River system within Ghana (the Oti, White Volta, and Afram rivers). Their analysis revealed that the three wild populations showed moderate levels of genetic differentiation (Afram and White Volta, FST = 0.049; Afram and Oti, FST = 0.027; and White Volta and Oti, FST = 0.054), suggesting that widening the sampling effort to include multiple sites within different rivers would reveal patterns of significant genetic variation within and among populations. In Ghana, the largest populations of O. niloticus occur within the Volta system, which drains about 70% of the country. Not surprisingly, many genetics studies have focused on O. niloticus populations in the Volta system in West Africa [4,15,16]. However, O. niloticus also occurs in almost all Ghanaian rivers, including the Pra, Densu, Ankobra, and Bia systems [17,18], and should be studied holistically.
Given the considerable historical hydrological connectivity apparent within the Volta River basin [19] and the moderate genetic differentiation observed among O. niloticus populations from the Oti, White Volta, and Afram rivers within the Volta River system [14], we expect considerable genetic differentiation among populations among different river basins in Ghana. Genetic differentiation may be the result of both non-selective (isolation, random genetic drift) and adaptive (natural selection, coadaptation) population genetic processes. Many rivers in West Africa are stressed from increased anthropogenic activities, such as alluvial mining, overfishing, water abstraction, and deforestation within the catchment and riparian zones, which pose significant threats to aquatic life. Without timely intervention, uncharacterized populations of O. niloticus could become extirpated, posing the loss of valuable genetic resources. Understanding the genetic variability within and among populations is vital for promoting long-term persistence of O. niloticus, especially because of growing threats to their persistence in the wild. Rognon and Guyomard [20] reported widespread genetic introgression by O. aureus into several West African O. niloticus populations. Similar results were reported for wild tilapia samples collected in Ghana, with evidence showing that pure populations of O. niloticus still persist in Ghanaian rivers [21]. Analysis of mitochondrial DNA sequence data from wild and farmed Oreochromis spp. collected in Ghana revealed the presence of non-native tilapia on some aquaculture facilities, including individuals showing introgression from non-native Oreochromis mossambicus, from which some stocks have escaped into the wild and interbred with native populations [22]. Given these threats, it is imperative that we assess population genetic data for wild populations of O. niloticus in a timely manner to inform conservation planning and management. This study was conducted to assess population genetic variation in O. niloticus within and among drainage basins in Ghana using nuclear microsatellite DNA markers to identify potential units of conservation among wild populations.

2. Materials and Methods

2.1. Fish Sampling and DNA Extraction

Fish morphologically identified as Nile tilapia Oreochromis niloticus were sampled from eight major rivers and one coastal lagoon in Ghana: Afram River, Oti River, White Volta River (two sites), Black Volta River (three sites), Pra River, Tano River (two sites), Ankobra River, Lower Volta River, and Juen Lagoon (Figure 1). A total of 312 individual fish collected from December 2014 through July 2017 (Table 1) were analyzed. Due to small sample sizes in the two Black Volta sites (Talewona and Lawra), all Black Volta individuals were combined and analyzed as one population represented by only one site location, Kantu, as seen in Figure 1. Total genomic DNA from fin-clips previously stored in paper envelopes was extracted at Virginia Tech, USA, following the methods described by Anane-Taabeah et al. [22].

2.2. Polymerase Chain Reaction and Genotyping

All fish samples were genotyped using polymerase chain reaction (PCR) at ten polymorphic microsatellite DNA loci developed for O. niloticus, which have shown high allelic variation in many studies [10,11,14,22,23]. Technical information for the respective microsatellite markers is presented in Table S1.
The final volume of 11 µL PCR amplification reactions consisted of 1 µL of 50–100 ng of genomic DNA, 0.1 µL of 5 U/µL of Taq DNA polymerase (Promega, Madison, WI, USA), 2 µL of 5× PCR buffer, 1 µL of 25 mM MgCl2, 0.4 µL of 2.5 mM dNTP mix, 0.4 µL of 1× bovine serum albumin, and 0.4 µL each of 5 µM of forward and reverse primers and 5.3 µL of dH2O. The following thermal cycling conditions were used: 94 °C for 3 min; 35 cycles of 94 °C for 40 s, 54–58 °C depending on marker for 30 s, and 72 °C for 1 min; a final extension of 72 °C for 5 min; and a 4 °C hold. PCR products were visualized using agarose gel electrophoresis to confirm amplification and amplicon sizes prior to genotyping with an ABI3730 automated DNA sequencer at the Virginia Tech Biocomplexity Institute (Blacksburg, VA, USA) or Cornell University (Ithaca, NY, USA). Amplification products were visualized with GeneMarker version 2.6.4 (SoftGenetics, State College, PA, USA) and scored by eye.

2.3. Data Analysis

Prior to analysis, we screened the data following the methods described in Anane-Taabeah et al. [22]. We checked for genotyping errors in the data using the program MICROCHECKER [24]. We then screened populations for deviation from linkage equilibrium [25] and Hardy–Weinberg equilibrium (HWE) at each locus using the program ARLEQUIN, version 3.1 [26]. The significance of departures from HWE was assessed using Fisher’s exact test with a Markov chain of 1,000,000 steps and 100,000 dememorization steps [27]. The significance of pairwise linkage tests was determined using the likelihood ratio test with 10,000 permutations [28]. Finally, we used the sequential Bonferroni correction to account for possible Type 1 errors associated with multiple pairwise evaluations [29].
We quantified genetic variation at the ten microsatellite loci using number of alleles per locus, observed and expected heterozygosities, allele frequencies at each locus, and number of private alleles. Allelic richness was calculated using FSTAT [30]. We used the Garza–Williamson index (M-ratio, the ratio of allelic richness to allelic range, the latter being the difference between the largest allele and the smallest allele per locus, [31]) to assess random genetic drift and FIS inbreeding coefficients to assess inbreeding within populations. We also estimated effective population size (Ne) using the program NeEstimator version 2.1 [32] and tested for recent bottleneck using the program BOTTLENECK [33].
We quantified population differentiation and structure using multiple methods. We used the program ARLEQUIN to calculate the FST [34] metric of differentiation and perform analysis of molecular variance (AMOVA, [35]). The null hypothesis was that all populations belong to a single group (no population differentiation) evident by high within-population variation, and the alternate hypothesis was significant population differentiation evident by high among-population variation. We calculated locus-by-locus genic differentiation using the Fisher’s exact G test executed in Genepop on the web version 4.2 (available at: http://genepop.curtin.edu.au/genepop_op3.html; accessed on 12 March, 2018, [36]). We also calculated locus-by-locus FST estimates in ARLEQUIN. We further estimated GST, a standardized measure of differentiation developed by Hendrick [37] and Jost’s D [38], an actual measure of differentiation, i.e., metrics which are not limited by the level of heterozygosity as is FST. Both GST and D were calculated using the program POPTREE2 [39]. We then compared the locus-by-locus FST estimates to the Fisher’s exact G test, GST, and D estimates.
We assessed population structure among populations using the program STRUCTURE version 2.3.4 [40]. We used both the admixture and no-admixture models both using location as a prior and not doing so in order to assess underlying patterns in the structure of the populations. We executed the models with 100,000 burn-in and 500,000 MCMC replicates with the ancestry-correlated option using default settings. We set the number of clusters K sequentially from 1 to 11 populations to aid in individual assignment analyses. The most likely K value was selected using both the Evanno et al. [41] and Pritchard and Wen [42] methods implemented online in STRUCTURE HARVESTER [43]. The Clumpak program was used to visualize the Structure plots [44]. We also tested for the likelihood of isolation by distance (IBD), that is, whether genetic dissimilarity among populations increased with increasing geographic distances among sites. We performed a Mantel test to assess IBD using Genepop on the web [36] and R© statistical software version 3.5.1 [45]. The general equation for a Mantel test statistic is:
Z m = i = 1 n j = 1 n g i j   ×   d i j
where gij and dij are the genetic and geographic distances, respectively, between populations i and j, given n populations [46]. For genetic distances, we used pairwise FST values transformed as FST(1 − FST); and for the geographic distances, we used both the raw distances between sites (km) and the log (ln)-transformed distances.

3. Results

3.1. Genetic Diversity within Wild Populations

We found that all ten microsatellite loci screened were polymorphic. Table S2 provides the allele frequencies for each of the eight microsatellite loci across sites. After analysis with MICROCHECKER, we found evidence of segregation of null alleles at locus UNH925 for 7 of the 11 populations from the wild, and departure from HWE consistently across all sites. We also found evidence of segregation of null alleles consistent with departure from HWE at locus UNH130 in some populations. We therefore excluded data from loci UNH925 and UNH130 from downstream analysis.
After Bonferroni correction (α = 0.05/28 pairs of comparison; adjusted α = 0.0018), significant linkage disequilibrium was observed at one pair of loci for the Tano (Elubo) and four pairs of loci for the Black Volta populations. Significant departures from HWE (α = 0.05) were observed in the Afram population at UNH180; White Volta–Binaba at UNH123, UNH858, UNH898, and UNH991; White Volta (Kulugu) at UNH123, UNH178, and UNH898; Oti at UNH123 and UNH180; Ankobra at UNH178; Tano (Asuhyea) at UNH178, Juen at UNH991; Black Volta at all loci except UNH991; and Lower Volta at UNH123, UNH180, and UNH934. Because such departures were distributed among the respective populations, data for all loci were retained in the analysis.
The summary statistics quantifying genetic variation across the eight loci for each site are provided in Table 1. Mean observed and expected heterozygosities were moderate to high across sites and ranged between 0.60 and 0.82 and 0.58 and 0.79, respectively (Table 1). The Black Volta, Ankobra, and Tano–Asuhyea populations had the lowest values, while the Afram and Tano-Elubo populations showed the highest heterozygosities. Similarly, the mean number of alleles per locus was moderate to high across sites and ranged from 5.88 to 9.75; the Tano–Asuhyea population showed the lowest number of alleles per locus, while Black Volta, Afram, and Lower Volta populations showed the highest mean number of alleles across the eight loci (≥9.0). Private alleles were observed at all loci for all sites except White Volta–Binaba and Tano–Asuhyea; numbers of private alleles ranged from one (Ankobra) to eight (Black Volta). The highest numbers of private alleles were observed at loci UNH123 (9) and UNH180 (7; Table S2). Locus-by-locus values of genetic diversity metrics are presented in Table S3.
Effective population sizes estimated using the linkage disequilibrium method ranged from 21 to infinity (Table 2). Seven populations—Aframso, White Volta–Binaba, Oti, Tano-Elubo, Juen, Black Volta, and Lower Volta rivers—had Ne estimates ranging between 21.4 and 178.2. Lower confidence interval estimates were sufficiently small that inbreeding and random drift could be strong population genetic forces within these populations. Upper confidence interval estimates for five populations were in the tens or low hundreds, suggesting a limited scope for adaptation over the long term. Upper confidence interval estimates for six populations, indicated as infinity, might better be regarded as undefined. These estimates could be bounded by the addition of more marker loci or more individuals to the analysis.
The results of the bottleneck tests show that the stepwise mutation model (SMM) was the more conservative in characterizing recent demographic declines within the populations than the infinite allele model (IAM). Three out of the four tests, including the sign test, standard differences test, and the Wilcoxon two-tailed test, under the SMM showed several bottlenecked populations. The sign test suggests that the White Volta–Kulugu, Oti, Pra, Ankobra, Black Volta and Lower Volta populations were bottlenecked (Table 3). The standard differences test suggested that all populations were bottlenecked, while the Wilcoxon two-tailed test suggested that all but the White Volta–Binaba and Tano–Asuhyea populations were bottlenecked. Mean M-ratios ranged between 0.37 and 0.49 for all populations across the eight loci; these low ratios (<0.7) supported the inference that genetic bottlenecks had occurred recently. Overall FIS ranged from −0.02 to 0.16 (Table 1); positive values of FIS indicate an excess of homozygotes, suggesting the occurrence of inbreeding. Locus-by-locus FIS values are reported in Table S3.

3.2. Genetic Differentiation among Wild Populations

Pairwise FST estimates revealed low to moderate genetic differentiation among sites (Table S4), which was statistically significant (p = 0.000). In general, Tano–Asuhyea was differentiated from all but the Tano–Elubo and Juen populations. The lowest differentiation was observed between the Oti and Afram populations (FST = 0.02), while the highest differentiation was observed between those of Tano–Asuhyea and Oti (FST = 0.18). The AMOVA results also support moderate to high differentiation among populations, with over 8% of the variance explained by differences among populations (Table S5). Genic differentiation (Fisher’s exact G test) and the locus-by-locus FST estimates all showed highly significant statistical differences (p < 0.000) across all loci for all sites combined (Table S4). Similarly, both Hendrick’s GST’ and Jost’s D supported high genetic differentiation among populations (Table S6).
STRUCTURE analysis showed significant differentiation among the populations at the 11 sites. After applying multiple analytic approaches as described in Methods above, the best-supported model was no admixture using location as a prior. The admixture model showed the highest support for grouping populations into four clusters (K = 4, results not shown), with several apparent migrants and admixed individuals at several sites. Results from use of the no-admixture model using location information, while generally comparable to those of the admixture model, were better supported in terms of higher LnP(D) values. The best-supported model for genetic structure using the Evanno et al. [41] method was K = 5, suggesting five clusters within and among individuals from the 11 sites (Figure 2). Genetic differentiation is often hierarchically organized, and it is often unrealistic to expect that there is one ‘true’ K that is best for modeling a particular data set (J. Pritchard, quoted in [47]). Therefore, viewing and reporting plots for multiple K-values is appropriate [48] because different values of K can give insights into different levels of structure. The LnP(D) values suggested that the most probable number of clusters was K = 10 (Figure 2, Table S7). STRUCTURE plots for both K = 5 and K = 10 are presented for comparison (Figure 2).
At K = 5 (Figure 2, top plot), the first cluster (shown in magenta) contained individuals from the two White Volta sites (Q, the level of individual membership in the inferred cluster, ranged between 0.92 and 0.99) and the Oti River (Q = 0.11). The second cluster (shown in blue) contained individuals mainly from Afram River (Q = 0.98) and Oti River (Q = 0.88). The third cluster (shown in purple) was almost exclusive to the Black Volta River (Q = 0.98). The fourth cluster (shown in orange) contained individuals from Pra River (Q = 0.97), Ankobra River (Q = 0. 98), Juen (Q = 0.20), and Lower Volta (Q = 0.91). Nearly 10% of individuals from the Lower Volta shared genotypes with individuals in the White Volta–Oti and Afram–Oti cluster, while almost no individuals shared genotypes with individuals from the Black Volta (Q = 0.0001). The fifth cluster (shown in green) contained individuals from Tano–Asuhyea River (Q ranged between 0.98 and 1.00), Tano–Elubo (Q = 0.97), and Juen (Q = 0.76). There was no signal of admixture or introgression of genes from the Black Volta population into individuals in the fifth cluster (Q = 0.0000), even though very weak signals were detected in the reverse direction (Q ranged between 0.0004 and 0.0079).
At K = 10 (Figure 2, bottom plot), each population clustered distinctly, with the exception of the two White Volta sites, the Pra, the Lower Volta, and the Black Volta. The two White Volta sites clustered together, with Q values similar to those observed when K = 5. The Pra and the Lower Volta clustered together. Interestingly, individuals in the Black Volta population grouped into two clusters. The first Black Volta cluster comprised 30 individuals originally sampled from Kantu. The second cluster comprised nine individuals sampled from Lawra and Talewona, which also shared genotypes with some individuals from Juen Lagoon (Q ranged between 0.14 and 0.20).
The Mantel test showed a significant positive relationship between genetic distances and geographic distances and provided modest support for isolation by distance (IBD) (adjusted r2 = 0.08, and p = 0.023, Figure 3). However, the genetic differentiation across space showed a weak relationship between genetic distances and geographic distances, even after log-transforming the geographic distances (adjusted r2 = 0.10, and p = 0.012) (Figure S1).

4. Discussion

4.1. Genetic Diversity within Wild O. niloticus Populations: Evidence of Recent Bottlenecks

The genetic diversity observed among wild populations of Nile tilapia in this study was higher than that found in a similar study conducted in Ghana narrowly within the Volta River basin [15]. However, the finding that significant departures from Hardy–Weinberg equilibrium (HWE) occurred at several loci for all populations except Pra River and Ankobra, and significant linkage disequilibrium (LD) observed, especially in the Black Volta River, suggests the occurrence of recent genetic bottlenecks within the populations. Further, the results of bottleneck tests conducted suggest that even the Ankobra and Pra populations were bottlenecked, and the mean M-ratio estimates in this study were below the empirical 0.7 threshold value used to assess population bottlenecks within populations [31]. It is plausible that there were natural historical bottlenecks and that the resulting demographic effect on populations was exacerbated with recent anthropogenic impacts. Effective population size estimates also suggested declining population numbers at several sites, which is alarming for a highly fecund species such as O. niloticus. Inbreeding among relatives is also likely due to the low to moderate and positive values of the inbreeding index (FIS ranged between 0.04 and 0.16). Deviations from HWE may be due to multiple causes, including genetic drift, inbreeding, or mixing of differentiated populations [21]. The inference of segregation of null alleles at three loci (UNH178, UNH203, and UNH898) in the Black Volta population also could have contributed to the departures from HWE and linkage equilibrium observed. Nonetheless, the Black Volta population requires further assessment and monitoring in order to understand the genetic dynamics within the population. This is especially pertinent since many private alleles were observed in the Black Volta population, as well as private mitochondrial DNA (mtDNA) haplotypes; and they are likely the main ancestral population of O. niloticus in Ghana [22].
The recent genetic bottlenecks inferred within the study populations may be the consequence of anthropogenic impacts upon the river systems in question. Many tilapia populations have been subject to heavy exploitation. Further, within the last decade, drastic habitat alteration has occurred in almost all freshwater habitats in Ghana, especially due to illegal alluvial gold mining. Turbidity levels of many rivers in Ghana are significantly above recommended limits for aquatic life, and cichlids in particular do not thrive under such environmental conditions. Dankwa et al. [49] found that in heavily mined areas (with turbidity levels between 100 and 1110 NTU), cichlids were completely absent, as were their planktonic prey. Physical and environmental impacts can affect reproduction and long-term persistence of O. niloticus in Ghana. High turbidity could affect the selection of mates, which also could result in inbreeding. Management effort should therefore focus on addressing these ecological stressors to ensure the long-term persistence of O. niloticus.

4.2. Genetic Differentiation among Wild Populations: Barriers to Dispersal and a Unique Dispersal Mechanism

The metrics used to assess genetic differentiation and structure among populations showed moderate genetic differentiation as a function of geographic location. The pairwise FST estimates showed that the Afram and Oti populations were the least differentiated. This finding was supported by the clustering of Afram River and Oti River populations into the same group by the STRUCTURE analysis. Even though moderately differentiated, some level of gene flow was apparent between O. niloticus populations in the Afram and Oti rivers across a distance of about 254 km (158 mi) between the two populations.
The clustering of the Oti and the White Volta populations was not as surprising because the distance between the Oti site and the farthest of the two White Volta sites (Binaba) is relatively short, 172 km (106 mi). The distance between the Afram site and the Lower Volta site is about 195 km (121 mi). The finding that some level of genetic exchange among Nile tilapia occurred between the Lower Volta and all the rivers within the Volta basin (except the Black Volta) suggests that isolation-by-distance is not as important in mediating gene flow within the Volta basin as the presence of physical barriers such as the 57-year-old Akosombo Dam. The presence of the Akosombo Dam and other irrigation reservoirs also could explain why the Black Volta River populations are highly differentiated from all other populations within the Volta basin. The results of the Mantel test further support the interpretation that isolation by distance alone could not explain the genetic structuring observed across space for the populations. For instance, significant genetic structuring was observed between the Black Volta sites even though two of the sites (Kantu and Talewona) are only 13 km (about 8 mi) apart.
On the other hand, both isolation-by-distance and physical barriers could be operating together to limit gene flow between the Black Volta and all sites outside the Volta basin. However, we observed that the Black Volta was more differentiated from the Tano–Asuhyea (FST = 0.143) than from the Tano–Elubo (FST = 0.083) population, which suggests some genetic exchange between the latter two population pairs. This result corroborates findings that the Black Volta population shared mitochondrial COI haplotypes with Tano–Elubo population and samples from Cote d’Ivoire [21,22] and suggests O. niloticus movement, likely through human intervention, between Ghana and Cote d’Ivoire at the shared border along the Black Volta River.
The pairwise FST estimates also revealed that the Tano–Asuhyea population was very distinct from all populations outside the Tano basin and moderately differentiated from the Tano–Elubo and Juen Lagoon populations. The apparent genetic differentiation from other populations, relatively low allelic richness, and the lack of private alleles within the Tano–Asuhyea population suggest the loss of diversity within the population due to genetic drift. This interpretation is supported by the low mean M-ratio estimated for the Tano–Asuhyea population (0.37) and significant bottleneck test results (Table 3). Occurring in the headwaters of the Tano River, the Tano–Asuhyea population is relatively isolated and may be exchanging relatively few migrants with the Tano–Elubo population.
STRUCTURE analysis also showed clustering of individuals from the Tano basin (Tano–Asuhyea, Tano–Elubo, and Juen Lagoon) and individuals from the Pra, Ankobra, and the Lower Volta rivers. Given that the Pra, Ankobra, Tano, and the Lower Volta rivers all run parallel to one another and drain into the Atlantic Ocean, oceanic dispersal appears to be a plausible mechanism for the genetic exchange apparent among O. niloticus populations. It appears that these populations are migrating from one freshwater environment through brackish water into another freshwater environment. Oceanic dispersal is possible because tilapias, in general, are euryhaline [50,51], and the salinity tolerance in O. niloticus has been documented [52,53,54]. Genetic data have been useful in detecting oceanic dispersal for species presumed to be restricted to freshwater environments, such as river sharks [55]. Oceanic dispersal also may explain the colonization of O. niloticus in the Ankobra and Pra rivers (from the Tano River [26]) and further explain why Tilapia pra is restricted to the Ankobra, Tano, and Pra rivers [18].

4.3. Defining Conservation Units for Wild Oreochromis niloticus from Nine River Basins in Ghana

Taken together, the genetic diversity and differentiation observed within and among the wild populations of Nile tilapia provide support for developing effective management and conservation plans for preserving these populations. The STRUCTURE results suggest that of the eleven populations studied, there may be as many as ten populations requiring separate management, which should be designated as management units (MUs). Additionally, the genetic distinctiveness of the Tano–Asuhyea and Black Volta populations and their isolation from other populations suggests the possibility that these populations may represent evolutionary significant units (ESUs) within O. niloticus populations in Ghana. Further studies focusing upon elucidation of any adaptive differences among candidate ESUs should help establish this possibility.
Conserving the Black Volta population in particular is important because of earlier findings from microsatellite data, which showed a lack of evidence of admixed individuals in the population, and from mtDNA analysis, which revealed that escaped, non-native, farmed tilapias are restricted to the Lower Volta and Afram rivers within the Volta basin [24]. The Lower Volta River is downstream of aquaculture operations on the Volta Lake, while there are aquaculture facilities in close proximity to the Afram River. On the other hand, the Black Volta River is relatively isolated from aquaculture farms even though cage farming in irrigation reservoirs is increasing in the region. Given the evidence of a possible human-mediated movement of O. niloticus between Ghana and Cote d’Ivoire at the shared border near the Black Volta River, if no conservation measures are instituted, it is only a matter of time before pure strains of O. niloticus in the Black Volta River are lost. Additionally, the possibility of ocean dispersal as a mechanism for migration of O. niloticus to otherwise isolated basins, restricting the dispersal of non-native tilapia populations in Ghana and along the West African coastline would be very challenging if non-native fish continue to escape from farms. The results from this study also provide crucial baseline information necessary for continuous monitoring of rivers close to aquaculture operations to conserve native O. niloticus populations.

4.4. Conclusions and Implications

The objective of this study was to assess genetic variation in O. niloticus within and among drainage basins in Ghana using nuclear microsatellite DNA markers to aid in detecting potential units of conservation within wild populations. We found moderate to high genetic diversity within and differentiation among wild populations, and we also detected strong signals of recent demographic and genetic bottlenecks in several populations. We also found significant genetic structure among the populations, suggesting the presence of up to 10 management units (MUs) within the 11 wild populations studied. In particular, the Black Volta and the Tano–Asuhyea populations, which were the most genetically distinct and geographically isolated and may be most at risk of loss of genetic diversity over time, may potentially represent evolutionary significant units. Therefore, at the minimum, both the Black Volta and Tano–Asuhyea populations should be designated as management units (MUs) and prioritized for conservation to sustain them in the long term.
Previous studies using allozyme [5] or mitochondrial DNA [20,22] markers showed introgression of O. aureus and O. mossambicus into some but not all western African O. niloticus populations. The microsatellite markers that we used are not in themselves powerful for detecting hybridization among Oreochromis species; mitochondrial and single-nucleotide polymorphism markers would be needed for assessment of interspecific hybridization.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/d14020073/s1, Table S1: Technical details for amplification of ten microsatellite loci for Ghanaian tilapias (Oreochromis spp.) collected from eight rivers and one coastal lagoon from December 2014 to July 2017 in Ghana, Table S2: Allele frequencies (%) across 8 nuclear microsatellite DNA loci for wild O. niloticus populations sampled from 11 sites (9 river basins) in Ghana from 2014 to 2017. AF = Afram, WB = White Volta–Binaba, WK = White Volta–Kulugu, OT = Oti, AN = Ankobra, TA = Tano–Asuhyea, TE = Tano–Elubo, JU = Juen, BV = Black Volta, and LV = Lower Volta. Private alleles are indicated in Bold, Table S3: Locus-by-locus genetic diversity metrics for all populations. Values in bold are significantly different from each other (p < 0.05), Table S4: Pairwise FST values from nuclear microsatellite DNA sequences for wild O. niloticus populations sampled from 11 sites (nine river basins) in Ghana from December 2014 to July 2017. All FST values were significant, Table S5: AMOVA for 8 nuclear DNA microsatellites loci in wild tilapia populations collected from 11 sites in Ghana from December 2014 to July 2017, Table S6: Genetic divergence statistics for each locus of a given sample size (n). Genic differentiation (G) from Fisher’s exact test, Hendrick’s GST, and Jost’s D are compared with FST values, Table S7: STRUCTURE results with estimates for probable K clusters for O. niloticus from Ghana. K = 10 as the most likely number of clusters is indicated in bold, Figure S1: Relationship between genetic distances and logarithm of geographic distances among Oreochromis niloticus collected from 11 sites in Ghana from December 2014 to July 2017.

Author Contributions

Conceptualization, G.A.-T.A., E.A.F., and E.M.H., methodology, G.A.-T.A., E.A.F., and E.M.H.; application of software, G.A.-T.A.; formal analysis, G.A.-T.A.; validation, G.A.-T.A., E.A.F., and E.M.H.; investigation, G.A.-T.A.; resources, G.A.-T.A., E.A.F., and E.M.H.; data curation, G.A.-T.A.; writing—original draft preparation, G.A.-T.A.; writing—review and editing, G.A.-T.A., E.A.F., and E.M.H.; visualization, G.A.-T.A., E.A.F., and E.M.H.; supervision, E.A.F., and E.M.H.; project administration, G.A.-T.A., E.A.F., and E.M.H.; funding acquisition, G.A.-T.A., and E.A.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the AquaFish Innovation Lab, supported by the U.S. Agency for International Development (USAID) under award number CA/LWA EPP-A-00-06-0012-00, and the Rufford Foundation (UK), under grant number 22083-1. G.A.-T.A. was supported in part by the Global Change Center (Virginia Tech Graduate School); Aquaculture Association of Ghana; Raanan Fish Feed, Ghana; Consulting in Aquaculture Research and Training for Development (CART4D); and the Burd Sheldon McGinnes Graduate Fellowship through the Department of Fish and Wildlife Conservation, Virginia Tech.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Ethics Committee of Kwame Nkrumah University of Science and Technology in a written letter in September 2014 (Protocol code number is not applicable).

Data Availability Statement

The data from this study have not been made publicly available yet. Some of the data can be accessed with written request to the authors.

Acknowledgments

We thank all the fishers, local guides, and personnel from the Fisheries Commission, Ghana, who collaborated on this project. We are grateful to the following individuals for their immense support in the field work in Ghana: Akwasi Ampofo-Yeboah, Nathaniel Gyasi Adjei, Francis Adjei, Anthony Aliebe, Abigail Ebachie Tarchie, and Emmanuel Nyamekye. We thank Miluska Olivera-Hyde and Sheila Harris (formerly of Virginia Tech) and Clay Ferguson of Virginia Tech assisted with laboratory work. Jess Jones and Donald Orth reviewed earlier versions of the manuscript and provided useful feedback.

Conflicts of Interest

The authors declare no conflict of interest, except those which may arise through the funding relationships described above.

References

  1. Food and Agriculture Organization (FAO) of the United Nations. The State of World Fisheries and Aquaculture 2020: Sustainability in Action; FAO: Rome, Italy, 2020. [Google Scholar] [CrossRef]
  2. Eknath, A.E.; Tayamen, M.M.; Palada-de Vera, M.S.; Danting, J.C.; Reyes, R.A.; Dionisio, E.E.; Capili, J.B.; Bolivar, H.L.; Abella, T.A.; Circa, A.V.; et al. Genetic improvement of farmed tilapias: The growth performance of eight strains of Oreochromis niloticus tested in different farm environments. Aquaculture 1993, 111, 171–188. [Google Scholar] [CrossRef]
  3. Eknath, A.E.; Dey, M.M.; Rye, M.; Gjerde, B.; Abella, T.A.; Sevilleja, R.; Tayamen, M.M.; Reyes, R.A.; Bentsen, H.B. Selective breeding of Nile tilapia for Asia. In Proceedings of the 6th World Congress of Genetics Applied to Livestock Production, Armidale, Australia, 11–16 January 1998; ICLARM Contribution No. 1397. International Center for Living Aquatic Resources Management: Manila, Philippines, 1998. [Google Scholar]
  4. Falk, T.M.; Abban, E.K. Genetic diversity of the Nile tilapia Oreochromis niloticus (Teleostei, Cichlidae) from the Volta System in Ghana. In Biodiversity, Management and Utilization of West African Fishes, WorldFish Center Conference Proceedings; Abban, E.K., Casal, C.M.V., Dugan, P., Falk, T.M., Eds.; WorldFish Center: Penang, Malaysia, 2004; pp. 13–15. [Google Scholar]
  5. 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. [Google Scholar] [CrossRef] [Green Version]
  6. Vreven, E.J.; Adepo-Gourene, B.; Agnèse, J.-F.; Teugels, G.G. Morphometric and allozyme variation in natural populations and cultured strains of Nile tilapia Oreochromis niloticus (Telesotei, Cichlidae). Belg. J. Zool. 1998, 128, 23–34. [Google Scholar]
  7. Agnèse, J.-F.; Adepo-Gourene, B.; Abban, E.K.; Fermon, Y. Genetic differentiation among natural populations of the Nile tilapia Oreochromis niloticus (Teleostei, Cichlidae). Heredity 1997, 79, 89–96. [Google Scholar] [CrossRef]
  8. Bezault, E.; Balaresque, P.; Toguyeni, A.; Fermon, A.; Araki, H.; Baroiller, J.-F.; Rognon, X. Spatial and temporal variation in population genetic structure of wild Nile tilapia (Oreochromis niloticus) across Africa. BMC Genet. 2011, 12, 102. [Google Scholar] [CrossRef] [PubMed]
  9. Lee, W.J.; Kocher, T.D. Microsatellite DNA markers for genetic mapping in Oreochromis niloticus. J. Fish Biol. 1996, 49, 169–171. [Google Scholar] [CrossRef]
  10. Carleton, K.L.; Streelman, J.T.; Lee, B.-Y.; Garnhart, N.; Kidd, M.; Kocher, T.D. Rapid isolation of CA microsatellites from the tilapia genome. Anim. Genet. 2002, 33, 140–144. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Streelman, J.T.; Kocher, T.D. Microsatellite variation associated with prolactin expression and growth of salt-challenged tilapia. Physiol. Genomics 2002, 9, 1–4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Moritz, C. Defining ‘evolutionary significant units’ for conservation. Trends Ecol. Evol. 1994, 9, 373–375. [Google Scholar] [CrossRef]
  13. Waples, R.S. Pacific salmon, Oncorynchus spp., and the definition of “species” under the Endangered Species Act. Mar. Fish. Rev. 1991, 53, 11–22. [Google Scholar]
  14. Frimpong, E.A.; Amisah, S.; Anane-Taabeah, G.; Ampofo-Yeboah, A.; Hallerman, E. Identifying Local Strains of Oreochromis Niloticus that Are Adapted to Future Climate Conditions; Report Submitted to United States Agency for International Development (USAID); Aquaculture Innovation Lab, Oregon State University: Corvallis, OR, USA, 2016; p. 27. [Google Scholar]
  15. Mireku, K.K.; Kassam, D.; Changadeya, W.; Attipoe, F.Y.K.; Adinortey, C.A. Assessment of genetic variations of Nile tilapia (Oreochromis niloticus L.) in the Volta Lake of Ghana using microsatellite markers. Afr. J. Biotechnol. 2017, 16, 312–321. [Google Scholar] [CrossRef] [Green Version]
  16. Lind, C.E.; Agyakwah, S.K.; Attipoe, F.Y.; Nugent, C.; Crooijmans, R.P.M.A.; Toguyeni, A. Genetic diversity of Nile tilapia (Oreochromis niloticus) throughout West Africa. Sci. Rep. 2019, 9, 16767. [Google Scholar] [CrossRef] [Green Version]
  17. Ansah, Y.B. Characterization of Pond Effluents and Biological and Physicochemical Assessment of Receiving Waters in Central Ghana. Master’s Thesis, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA, 2010. [Google Scholar]
  18. Dunz, A.R.; Schliewen, U. Description of a new species of Tilapia Smith, 1840 (Teleosti: Cichlidae) from Ghana. Zootaxa 2010, 2548, 1–21. [Google Scholar] [CrossRef]
  19. McCartney, M.; Forkuor, G.; Sood, A.; Amisigo, B.; Hattermann, F.; Muthuwatta, L. The Water Resource Implications of Changing Climate in the Volta River Basin; IWMI Research Report; International Water Management Institute: Colombo, Sri Lanka, 2012; Volume 146, p. 40. [Google Scholar] [CrossRef] [Green Version]
  20. Rognon, X.; Guyomard, R. Large extent of mitochondrial DNA transfer from Oreochromis aureus to O. niloticus in West Africa. Mol. Ecol. 2003, 12, 435–445. [Google Scholar] [CrossRef]
  21. Anane-Taabeah, G. Characterization of the Molecular Genetic Variation in Wild and Farmed Nile Tilapia (Oreochromis niloticus) in Ghana for Conservation and Aquaculture Development. Ph.D. Thesis, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA, 2018. [Google Scholar]
  22. Anane-Taabeah, G.; Frimpong, E.A.; Hallerman, E. Aquaculture-mediated invasion of the Genetically Improved Farmed Tilapia (GIFT) into the Lower Volta Basin of Ghana. Diversity 2019, 11, 188. [Google Scholar] [CrossRef] [Green Version]
  23. Cnaani, A.; Hallerman, E.M.; Ron, M.; Weller, J.I.; Kashi, Y.; Gall, G.A.E.; Hulata, G. A chromosomal region with quantitative trait loci affecting cold tolerance and fish size in an F2 tilapia hybrid. Aquaculture 2003, 223, 117–128. [Google Scholar] [CrossRef]
  24. Van Oosterhout, C.; Hutchinson, W.F.; Willis, D.P.M.; Shipley, P. Microchecker: Software for identifying and correcting genotyping errors in microsatellite data. Mol. Ecol. Notes 2004, 4, 535–538. [Google Scholar] [CrossRef]
  25. Sokal, R.R.; Rohlf, F.J. Biometry: The Principles and Practice of Statistics in Biological Research, 3rd ed.; W.H. Freeman: New York, NY, USA, 1994. [Google Scholar]
  26. Excoffier, L.; Laval, G.; Schneider, S. ARLEQUIN, version 3.0: An integrated software package for population genetics data analysis. Evol. Bioinform. 2005, 1, 47–50. [Google Scholar] [CrossRef] [Green Version]
  27. Guo, S.; Thompson, E. Performing the exact test of Hardy–Weinberg proportion for multiple alleles. Biometrics 1992, 48, 361–372. [Google Scholar] [CrossRef] [PubMed]
  28. Slatkin, M.; Excoffier, L. Testing for linkage disequilibrium in genotypic data using the EM algorithm. Heredity 1996, 76, 377–383. [Google Scholar] [CrossRef] [Green Version]
  29. Rice, W.R. Analyzing tables of statistical tests. Evolution 1989, 43, 223–225. [Google Scholar] [CrossRef]
  30. Goudet, J. Fstat (ver. 2.9.4), a Program to Estimate and Test Population Genetics Parameters; Lausanne, Switzerland, 2003; Available online: http://www.unil.ch/izea/softwares/fstat (accessed on 15 November 2021).
  31. Garza, J.C.; Williamson, E.G. Detection of reduction in population size using data from microsatellite loci. Mol. Ecol. 2001, 10, 305–318. [Google Scholar] [CrossRef] [PubMed]
  32. Do, C.; Waples, R.S.; Peel, D.; Macbeth, G.M.; Tillet, B.J.; Ovenden, J.R. NeEstimator v2: Re-implementation of software for the estimation of contemporary effective population size (Ne) from genetic data. Mol. Ecol. Resour. 2014, 14, 209–214. [Google Scholar] [CrossRef]
  33. Cornuet, J.M.; Luikart, G. Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics 1997, 144, 2001–2014. [Google Scholar] [CrossRef] [PubMed]
  34. Wright, S. Evolution and the Genetics of Populations, Vol. 4: Variability within and among Natural Populations; University of Chicago Press: Chicago, IL, USA, 1978. [Google Scholar]
  35. Excoffier, L.; Smouse, P.E.; Quattro, J.M. Analysis of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data. Genetics 1992, 131, 479–491. [Google Scholar] [CrossRef]
  36. Raymond, M.; Rousset, F. GENEPOP (version 1.2): Population genetics software for exact tests and ecumenicism. J. Hered. 1995, 86, 248–249. [Google Scholar] [CrossRef]
  37. Hedrick, P.W. A standardized genetic differentiation measure. Evolution 2005, 59, 1633–1638. [Google Scholar] [CrossRef]
  38. Jost, L. GST and its relatives do not measure differentiation. Mol. Ecol. 2008, 17, 4015–4026. [Google Scholar] [CrossRef]
  39. Takezaki, N.; Nei, M.; Tamura, K. POPTREE2: Software for constructing population trees from allele frequency data and computing some other population statistics with Windows interface. Mol. Biol. Evol. 2010, 27, 747–752. [Google Scholar] [CrossRef] [Green Version]
  40. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar] [CrossRef]
  41. Evanno, G.; Regnaut, S.; Goudet, J. Detecting the number of clusters of individuals using the software structure: A simulation study. Mol. Ecol. 2005, 14, 2611–2620. [Google Scholar] [CrossRef] [Green Version]
  42. Pritchard, J.K.; Wen, W. Documentation for STRUCTURE Software: Version 2; University of Chicago Press: Chicago, IL, USA, 2003; Available online: http://web.stanford.edu/group/pritchardlab/software/structure2_1.html (accessed on 12 March 2018).
  43. Earl, D.A.; von Holdt, B.M. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 2012, 4, 359–361. [Google Scholar] [CrossRef]
  44. Kopelman, N.M.; Mayzel, J.; Jakobsson, M.; Rosenberg, N.A.; Mayrose, I. Clumpak: A program for 676 identifying clustering modes and packaging population structure inferences across K. Mol. Ecol. Resour. 2015, 15, 1179–1191. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing; Vienna, Austria, 2018; Available online: https://www.R-project.org/ (accessed on 12 March 2018).
  46. Mantel, N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27, 209–220. [Google Scholar]
  47. Benestan, L.M.; Ferchaud, A.L.; Hohenlohe, P.A.; Garner, B.A.; Naylor, G.J.; Baums, I.B.; Schwartz, M.K.; Kelley, J.L.; Luikart, G. Conservation genomics of natural and managed populations: Building a conceptual and practical framework. Mol. Ecol. 2016, 25, 2967–2977. [Google Scholar] [CrossRef]
  48. Gilbert, K.J.; Andrew, R.L.; Bock, D.G.; Franklin, M.T.; Kane, N.C.; Moore, J.-S.; Moyers, B.T.; Renaut, S.; Rennison, D.J.; Veen, T.; et al. Recommendations for utilizing and reporting population genetic analyses: The reproducibility of genetic clustering using the program STRUCTURE. Mol. Ecol. 2012, 21, 4925–4930. [Google Scholar] [CrossRef]
  49. Dankwa, H.R.; Biney, C.A.; De Graft-Johnson, A.A. Impact of mining operations on the ecology of River Offin in Ghana. West Afr. J. Appl. Ecol. 2005, 7, 19–30. [Google Scholar] [CrossRef] [Green Version]
  50. Philippart, J.C.; Ruwet, J.C. Ecology and distribution of tilapias. The Biology and Culture of Tilapias. In ICLARM (International Center for Living Aquatic Resources Management) Conference Proceedings; Pullin, R.S.V., Lowe-McConnell, R.H., Eds.; International Center for Living Aquatic Resources Management: Bellagio, Italy, 1982; Volume 7, pp. 15–59. [Google Scholar]
  51. Stickney, R. Tilapia tolerance in saline waters: A review. Prog. Fish-Cult. 1986, 48, 61–67. [Google Scholar] [CrossRef]
  52. Watanabe, W.O.; Kuo, C.-M.; Huang, M.-C. Salinity tolerance of Nile tilapia fry (Oreochromis niloticus), spawned and hatched at various salinities. Aquaculture 1985, 48, 159–176. [Google Scholar] [CrossRef]
  53. Likongwe, J.S.; Stecko, T.D.; Stauffer, J.R., Jr.; Carline, R.F. Combined effects of water temperature and salinity on growth and feed utilization of juvenile Nile tilapia Oreochromis niloticus (Linneaus). Aquaculture 1996, 146, 37–46. [Google Scholar] [CrossRef]
  54. Kamal, A.H.M.M.; Mair, G.C. Salinity tolerance in superior genotypes of tilapia, Oreochromis niloticus, Oreochromis mossambicus and their hybrids. Aquaculture 2005, 247, 189–201. [Google Scholar] [CrossRef]
  55. Li, C.; Corrigan, S.; Yang, L.; Straube, N.; Harris, M.; Hofreiter, M.; White, W.T.; Naylor, G.J.P. DNA capture reveals transoceanic gene flow in endangered river sharks. Proc. Natl. Acad. Sci. USA 2015, 112, 13302–13307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Sampling locations for wild Oreochromis niloticus collected in Ghana at 11 sites from eight rivers and one coastal lagoon from December 2014 to July 2017. Sampling sites are indicated with red triangles.
Figure 1. Sampling locations for wild Oreochromis niloticus collected in Ghana at 11 sites from eight rivers and one coastal lagoon from December 2014 to July 2017. Sampling sites are indicated with red triangles.
Diversity 14 00073 g001
Figure 2. Analyses of wild samples of O. niloticus sampled from Ghana with STRUCTURE results showing the proportion of individuals from each wild O. niloticus population’s ancestry inferred from each of K = 5 (top plot) and K = 10 (bottom plot) clusters (MCMC = 500,000).
Figure 2. Analyses of wild samples of O. niloticus sampled from Ghana with STRUCTURE results showing the proportion of individuals from each wild O. niloticus population’s ancestry inferred from each of K = 5 (top plot) and K = 10 (bottom plot) clusters (MCMC = 500,000).
Diversity 14 00073 g002
Figure 3. Relationship between genetic distance and geographic distance among wild Oreochromis niloticus collected from 11 sites in Ghana from December 2014 and July 2017.
Figure 3. Relationship between genetic distance and geographic distance among wild Oreochromis niloticus collected from 11 sites in Ghana from December 2014 and July 2017.
Diversity 14 00073 g003
Table 1. Summary of genetic variation among eight microsatellite DNA loci examined in wild tilapia populations from 11 sites (nine river basins) collected in Ghana from December 2014 to July 2017. N = number of individuals genotyped per locus, Ho = observed heterozygosity, He = expected heterozygosity, A = mean number of observed alleles per locus, Ar = allelic richness, PrA = number of private alleles, Range = base-pair differences between the shortest and longest microsatellite alleles observed per locus, M-ratio = approximate ratio of A and Range, and FIS = inbreeding coefficient.
Table 1. Summary of genetic variation among eight microsatellite DNA loci examined in wild tilapia populations from 11 sites (nine river basins) collected in Ghana from December 2014 to July 2017. N = number of individuals genotyped per locus, Ho = observed heterozygosity, He = expected heterozygosity, A = mean number of observed alleles per locus, Ar = allelic richness, PrA = number of private alleles, Range = base-pair differences between the shortest and longest microsatellite alleles observed per locus, M-ratio = approximate ratio of A and Range, and FIS = inbreeding coefficient.
BasinSampling LocationYear SampledPopulation IDNHoHeAArPrARangeM-RatioFIS
Afram RiverAframso2014; 2015AF300.750.739.006.75215.550.49−0.03
White Volta RiverBinaba2014; 2015WB300.650.727.385.87016.750.410.10
White Volta RiverKulugu2014; 2015WK300.630.708.386.27317.630.440.09
Oti RiverSabare2014; 2015OT300.670.688.636.07419.880.410.02
Pra RiverTwifo-Praso2017PR300.670.658.386.10218.000.46−0.03
Ankobra RiverAnkobra2017AN300.600.586.634.90116.750.39−0.02
Tano RiverAsuhyea2017TA300.650.645.884.78015.630.37−0.02
Tano RiverElubo2017TE180.820.798.387.48518.500.45−0.04
Juen LagoonJehwi-Wharf2017JU120.690.747.637.63618.130.390.07
Black Volta RiverKantu **2017BV390.600.719.256.37820.750.450.16
Lower Volta RiverNotreku-Akuse2017LV330.690.729.756.77418.130.390.05
** The majority of Black Volta River samples were obtained from Kantu (N = 30). Samples were also collected from Lawra (N = 3) and Talewona (N = 6).
Table 2. Effective population size (Ne) estimates and 95% confidence intervals for wild O. niloticus population samples 1 collected in Ghana. ∞ = infinite. AF = Afram, WB = White Volta–Binaba, WK = White Volta–Kulugu, OT = Oti, AN = Ankobra, TA = Tano–Asuhyea, TE = Tano–Elubo, JU = Juen, BV = Black Volta, and LV = Lower Volta.
Table 2. Effective population size (Ne) estimates and 95% confidence intervals for wild O. niloticus population samples 1 collected in Ghana. ∞ = infinite. AF = Afram, WB = White Volta–Binaba, WK = White Volta–Kulugu, OT = Oti, AN = Ankobra, TA = Tano–Asuhyea, TE = Tano–Elubo, JU = Juen, BV = Black Volta, and LV = Lower Volta.
Population IDEstimated Ne
AF43.3 (25.1, 102.6)
WB21.4 (14.0, 36.2)
WK∞ (104.4, ∞)
OT46.4 (23.7, 182.2)
PR1184.4 (52.8, ∞)
AN∞ (58.1, ∞)
TA∞ (120.3, ∞)
TE28.1 (15.8, 74.1)
JU50.1 (12.4, ∞)
BV22.8 (14.9, 37.6)
LV178.2 (43.2, ∞)
1 AF = Afram, WB = White Volta–Binaba, WK = White Volta–Kulugu, OT = Oti, AN = Ankobra, TA = Tano–Asuhyea, TE = Tano–Elubo, JU = Juen, BV = Black Volta, and LV = Lower Volta.
Table 3. Summary statistics and p-values for bottleneck tests conducted for wild O. niloticus populations 1. Tests with significant p-values (α = 0.05) are indicated in bold. Sign = sign test, stdv = standard differences test, W_1t = Wilcoxon one-tailed test, W_2t = Wilcoxon two-tailed test. IAM = infinite allele model, SMM = stepwise mutation model.
Table 3. Summary statistics and p-values for bottleneck tests conducted for wild O. niloticus populations 1. Tests with significant p-values (α = 0.05) are indicated in bold. Sign = sign test, stdv = standard differences test, W_1t = Wilcoxon one-tailed test, W_2t = Wilcoxon two-tailed test. IAM = infinite allele model, SMM = stepwise mutation model.
Populationp_sign_IAMp_sign_SMMp_stdv_IAMp_stdv_SMMp_W_1t_IAMp_W_1t_SMMp_W_2t_IAMp_W_2t_SMM
AF0.3090.0580.1830.0000.1560.9900.3130.027
WB0.0970.1860.0350.0270.0200.7700.0390.547
WK0.5800.0010.2920.0000.1561.0000.3130.004
OT0.4290.0010.4890.0000.5781.0000.9450.004
PR0.4070.0010.0470.0000.7271.0000.6410.004
AN0.4550.0100.1830.0000.6800.9960.7420.012
TA0.2520.4260.0650.0070.0270.8090.0550.461
TE0.1060.0540.0590.0310.0100.9860.0200.039
JU0.5950.0520.4630.0010.5270.9861.0000.039
BV0.5740.0010.4690.0000.3201.0000.6410.004
LV0.1730.0100.3030.0000.6800.9940.7420.020
1 AF = Afram, WB = White Volta–Binaba, WK = White Volta–Kulugu, OT = Oti, AN = Ankobra, TA = Tano–Asuhyea, TE = Tano–Elubo, JU = Juen, BV = Black Volta, and LV = Lower Volta.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Anane-Taabeah Attu, G.; Frimpong, E.A.; Hallerman, E.M. Defining Management Units for Wild Nile Tilapia Oreochromis niloticus from Nine River Basins in Ghana. Diversity 2022, 14, 73. https://doi.org/10.3390/d14020073

AMA Style

Anane-Taabeah Attu G, Frimpong EA, Hallerman EM. Defining Management Units for Wild Nile Tilapia Oreochromis niloticus from Nine River Basins in Ghana. Diversity. 2022; 14(2):73. https://doi.org/10.3390/d14020073

Chicago/Turabian Style

Anane-Taabeah Attu, Gifty, Emmanuel A. Frimpong, and Eric M. Hallerman. 2022. "Defining Management Units for Wild Nile Tilapia Oreochromis niloticus from Nine River Basins in Ghana" Diversity 14, no. 2: 73. https://doi.org/10.3390/d14020073

APA Style

Anane-Taabeah Attu, G., Frimpong, E. A., & Hallerman, E. M. (2022). Defining Management Units for Wild Nile Tilapia Oreochromis niloticus from Nine River Basins in Ghana. Diversity, 14(2), 73. https://doi.org/10.3390/d14020073

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop