Genetic Identification and Diversity of Stocks of The African Bonytongue, Heterotis Niloticus (Osteoglossiformes: Arapaiminae), in Nigeria, West Africa

doi:10.21203/rs.3.rs-944303/v1

Download PDF

Research Article

Genetic Identification and Diversity of Stocks of The African Bonytongue, Heterotis Niloticus (Osteoglossiformes: Arapaiminae), in Nigeria, West Africa

https://doi.org/10.21203/rs.3.rs-944303/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Inland fisheries are an important source of protein and income for people in Africa, and sustainable fisheries management requires knowledge of regional stocks. Presently, information regarding genetic diversity of African freshwater fishes is very limited. Here we present findings from a spatial analysis of genetic patterns in the African bonytongue, Heterotis niloticus, an important component of inland fisheries in West Africa. The species is common throughout Western and Central Africa and the Nile Basin; and has been introduced for aquaculture in locations across Africa. Nigeria has the largest fishery for African bonytongues, representing ~86% of the global total. Recent declines in yields at some Nigerian locations, however, suggest current levels of exploitation may be unsustainable. Habitat degradation also may be impacting some stocks. Despite its commercial and nutritional importance, the African bonytongue has been the subject of scant research to support management. We examined patterns of genetic diversity in natural populations of H. niloticus at four locations in Nigeria, including Kainji Lake, a reservoir on the Niger River in north-central Nigeria, and three southtern localities (Ethiope River, Igbokoda River, and Epe Lagoon), as well fish from the Ouémé River delta near Porto Novo, Benin. Eighty-five specimens were genotyped for nine microsatellite-loci. Genetic diversity estimates were highest at Kainji Lake, and substantially lower at southern localities. High levels of genetic differentiation were detected between samples from Kainji Lake and those from the southern localities. Low, yet significant F_ST values were observed among samples from southern Nigerian localities that were more differentiated from the sample from nearby coastal Benin. We thus recommend that African bonytongues from the five locations be considered distinct genetic stocks and managed accordingly.

Genetic

diversity

bonytongue

exploitation

Inland fisheries provide an important source of protein, employment and income for people in Africa [1-4]. Capture inland fisheries in Africa produced an estimated 3,004,437 metric tons in 2018, an amount only surpassed by that from Asia at 7,953,840 metric tons [5]. Inland fisheries in Africa, however, are facing growing threats that include overexploitation, habitat degradation, and introduction of exotic species that can impact native fish stocks [6,7]. Given that human populations and economies are growing rapidly in Africa, impacts from these threats are expected to increase [7,8]. The human population in sub-Saharan Africa is projected to double by 2050 [9], and it is estimated that one in every four people will be undernourished in this region [10]. Inland fisheries will continue to play a critical role in supporting food security and livelihoods in growing populations and economies of the continent. Research that documents patterns of genetic diversity of inland fisheries species is needed for identification of management and/or conservation units and to monitor impacts from fishing [11]. At present, genetic studies are very limited for most African freshwater fishes, even though such studies have been long recognized as a research priority for the sustainable exploitation of these valuable resources [12].

With ~206 M people, Nigeria is the most populous country in Africa, ranks 7th among nations globally, and is one of nine countries that collectively will contribute more than half the projected growth of the global population between now and 2050 [9]. Freshwater fish are important in the diet of Nigerians [2], and demand for fish has increased considerably with population growth. An estimated 420,078 metric tons of freshwater fish were caught in Nigeria in 2017 and 392,188 metric tons in 2018, ranking first and second, respectively, among harvests by African countries [13]. These amounts are nearly four-fold greater than the amount harvested in 1980, when Nigeria’s population was ~73.4 M, and 107,530 metric tons of freshwater fishes were captured. Nigeria is one of the seven countries that are driving most of the growth in global inland fisheries [2], and by 2050 could be one of the world's largest consumers of fish [14]. Abundance of many freshwater fishes in Nigeria, however, appears to be declining due to overfishing [15]. Of the 266 species documented in 2005, 47 (17%) were reported to be critically endangered, 15 (5%) endangered, 8 (3%) vulnerable, 23 (8%) near threatened, and the status of 129 (48%) was unknown [6].

The African bonytongue, Heterotis niloticus (Cuvier, 1829), is a freshwater fish widely distributed in tropical rivers and freshwater lakes of Western and Central Africa and the Nile Basin where it supports important commercial and subsistence fisheries [16]. The African bonytongue has been introduced for aquaculture at several locations across Africa [17] and is considered one of the top ten species for future use in global aquaculture [18]. The greatest harvest of African bonytongues occurs in Nigeria [13]. Between 1990 and 2019, ~362,534 metric tons of this fish were caught in Nigeria, corresponding to 6.6% of the national freshwater fish yield, and ~83.4% of the total yield for this species from the entire continent. Recent bonytongue yields represent a sharp increase in comparison to previous decades. For example, in 2010, 15,392 metric tons were caught in Nigeria, representing 85.2% of the total yield for this species from the continent, and 5.2% of the national freshwater fish yield; whereas in 2000, only 7,915 metric tons were caught in Nigeria. Current levels of exploitation of this species may be unsustainable in regions of Nigeria where catches appear to be declining [19].

The only previous research on the population genetics of the African bonytongue was conducted in Benin, the country on the western border of Nigeria [20]. That study detected high levels of genetic differentiation among fish sampled from three river basins: (1) the Niger River at Malanville, a city in northern Benin located upstream from where the river enters Nigeria; (2) a location on the lower Mono River in coastal south-west Benin; and (3) nine localities in the Ouémé-Sô floodplain system in south-central Benin. Within the Ouémé–So river-floodplain system, however, low levels of genetic differentiation were detected among the localities sampled, suggesting seasonal flooding facilitates gene flow. Low genetic variability was detected for African bonytongues at Malanville, in the Niger River, compared to fish from the Mono and Ouémé-Sô basins [20]. Malanville is a fishing community where some stocks may be overexploited [21]. Therefore, it is important to examine other localities in the Niger Basin to determine whether low genetic diversity is a characteristic of H. niloticus in that system.

Given the significance of the African bonytongue fishery in Nigeria, and its potential for genetic differentiation, it is critical to evaluate genetic patterns within exploited stocks and to identify appropriate management units. Given that fishing can erode genetic diversity of inland fish stocks, estimation of genetic diversity is important for assessment of population viability and evolutionary potential. Herein, we examined genetic diversity and differentiation among natural populations of H. niloticus from four locations in Nigeria where fishing activity is intense. These include Kainji Lake, in north central Nigeria, a reservoir on the Niger River formed by the Kainji Dam, and three localities in southwestern Nigeria close to the Atlantic coast: Ethiope River, Igbokoda River, and Epe Lagoon (Figure 1). We also analyzed samples from the lower Ouémé River near Porto Novo in coastal southern Benin. Kainji Lake is approximately 230 Km downstream from Malanville and approximately 700 Km upstream from the massive Niger River delta, which forms ~240 Km before the river empties into the Atlantic. Within the delta, the Niger River channel divides to form an intricate network of waterways surrounded by extensive floodplains that spread along the coast for ~320 Km and cover an area of ~36,000 Km². The Ethiope and Igbokoda rivers are located in the western portion of the delta. Epe Lagoon, which is close to the Igbokoda River, and Porto Novo are located west of the delta. Given the large distance separating Kainji Lake and the southern localities, we hypothesized there would be very little gene flow, and therefore relatively large genetic differentiation between samples from the two areas. Furthermore, Kainji Dam, which was completed in 1968, provides an effective barrier to upstream fish movement, although the spillway may allow for some degree of downstream passage. Relatively low levels of genetic differentiation are expected among southern coastal localities where high fluvial connectivity and annual flooding could facilitate fish dispersal.

Sampling and DNA extractions

Samples of wild H. niloticus specimens were obtained from four inland waters in Nigeria: Kainji Lake (including samples from an area just above the lake and a location below the dam), Igbokoda River, Epe Lagoon and Ethiope River (Figure 1). Additional specimens were obtained from the Ouémé River floodplain near Porto Novo, Benin. Dead fish were purchased from local commercial fishermen, and in some cases only pectoral fin clips were obtained from fish at the point of sales to consumers. Animal ethics approval was not required for this study because live fish samples were not collected. A pectoral fin clip from each specimen was stored in an Eppendorf tube with 70% ethanol. Genomic DNA was extracted from the pectoral fin clips of the tissue samples using the Qiagen DNeasy Blood & Tissue Kit (Valencia CA).

Microsatellites amplification and genotyping

Samples were genotyped for nine microsatellite loci (Table SM1) reported for the African bonytongue [20,22]. Amplification of the extracted DNA was accomplished through multiplexing using Qiagen’s Multiplex PCR kit (Multiplex PCR Kit, QIAGEN, Valencia, CA) carried out in 10 µl reaction volumes containing 50–1000 ng genomic DNA, 5 µl QIAGEN Multiplex PCR Master Mix (containing HotStarTaq® DNA Polymerase, Multiplex PCR buffer and dNTP Mix), 1 µl of the forward (fluorescent labelled) and reverse primer mix and 3 µl of nuclease-free water (Qiagen Sample and Assay Technologies, Germany). Amplifications were carried out using a BIO-RAD Thermal Cycler (T100TM) programmed according to standard multiplex PCR protocols as follows: initial activation step: 15 min at 95°C, 35 cycles of denaturation at 94°C for 30 s, annealing at 57°C for 90 s, extension at 72°C for 60 s, and a final extension step 60°C for 30 min. PCR products were run on 2% agarose gel electrophoresis to check for amplification of all loci prior to genotyping, and positive PCR products were then run on an ABI 3130xl automated sequencer using Genescan 400 HD ROX as size standard (Applied Biosystems). Allele calling and sizing for each sample was obtained using Genemarker software (version 1.75).

Genetic diversity and population differentiation analyses

To convert input files between software packages, we used PGDSpider v. 2.1.1.3 [23]. Tests of linkage disequilibrium (LD) and test of conformity to the expectations of Hardy-Weinberg Equilibrium proportions (HWP) were performed using GENEPOP on the Web 4.7 [24,25]. To control for the occurrence of false positives due to multiple comparisons, we used the Benjamini-Hochberg procedure [26] at a false discovery rate (FDR) of 0.05 for assessing the significance of LD and HWP p-values. All GENEPOP analyses were performed with a dememorization number of 10,000 and 1,000 batches with 10,000 iterations per batch. GenAIEx v. 6.5 [27,28] was used to calculate the following genetic diversity statistics: number of private alleles (N_P), number of private alleles average number of alleles (Na), average number of effective alleles (Ne), Shannon’s information index (I), mean observed heterozygosity (H_O), mean expected heterozygosity (He), unbiased expected heterozygosity (uHe), and inbreeding coefficient (F_IS). We used MICRO-CHECKER v. 2.2.3 [29] to check for the potential presence of null alleles and scoring errors.

Population genetic structure was analyzed using several different approaches. Pairwise population F_ST values, with and without correction for null alleles, and their corresponding 95% confidence intervals (CI) were estimated using the program FreeNA [30]. F_ST values in which the 95% CIs excluded zero were deemed statistically significant. Pairwise population F_ST values were also estimated using the program Arlequin v.3.5.2.2 [31], which does not correct for null alleles, and their corresponding p-values were estimated with 10,000 permutations. The Benjamini-Hochberg procedure at a false discovery rate (FDR) of 0.05 was used to determine whether these F_ST values were statistically significant. A discriminant analysis of principal component analysis (DAPC) with the R package adegenet v.2.1.1 [32,33] was conducted to visualize levels of genetic differentiation among the samples from the different locations, defining sampling locations as a priori groups. DAPC is a non-model-based method that maximizes the differences between groups while minimizing variation within groups [34]; and distances between genetic clusters using a priori groupings reflect underlying F_ST [35]. STRUCTURE 2.2.3 [36], which performs model-based clustering with a Bayesian approach, was also used to examine population subdivision. Two models were used: admixture with correlated allele frequencies, and no admixture with correlated allele frequencies. K values from 1 to 5 were tested in ten iterations, with 500,000 steps and a burn-in of 125,000 steps, and all other settings were set to default. GenAlEx was used to construct a genetic distance matrix, from which a principal coordinates analysis (PCoA) was performed to identify population clusters.

Finally, we analyzed the relationship between genetic differentiation and geographic distance by testing for isolation by distance (IBD) using ISOLDE in GENEPOP on the Web 4.7. Two kinds of analyses were performed: (1) between localities; and, (2) between individuals. We used Google Earth Pro to measure distances between sample locations. For the distances between locations within the southern rivers, we measured the shortest straight-line distance. For distances between the Kainji Lake and all other localities, we followed the contour of the Niger River to the closest point to the Ethiope River, and from there we used straight lines. Statistical significance based on the Spearman’s rank correlation coefficient was evaluated using Mantel tests [37].

Genetic diversity

Genotyping scores for all individuals are shown in Dataset S1. Complete genotypes for the nine microsatellite-loci were obtained for 83 specimens: 23 from Kainji Lake, 19 from Ethiope River, 15 from Igbokoda, 20 from Epe Lagoon, and six from Porto Novo (Table 1). No linkage disequilibrium was detected among loci. A total of 98 individual alleles were observed across these 83 specimens; and the number of alleles per locus ranged from 4 in locus Hn45 to 19 in Hn30 (Table SM2). Departures of Hardy-Weinberg equilibrium were detected for locus Hn47 in Kainji Lake, and locus Hn11 in Igbokoda. According to MICROCHECKER, only locus Hn11 in the Igbokoda sample shows signs of a null allele, which is suggested by a general excess of homozygotes for most allele size classes. Excess of homozygotes at this locus, however, may reflect inbreeding at this locality (see below).

Table 1

Genetic diversity estimates for *H. niloticus* from five localities in Nigeria and Benin
Locality	N	N_A	N_P	Na	Ne	I	Ho	He	uHe	F
Kainji	23	75	21	8.33	4.56	1.60	0.73	0.70	0.72	-0.03
Ethiope River	19	44	2	4.89	2.58	1.01	0.50	0.50	0.51	-0.01
Igbokoda	15	51	4	5.67	2.55	1.11	0.47	0.54	0.56	0.18
Epe Lagoon	20	46	3	5.11	2.65	0.97	0.44	0.47	0.48	0.06
Porto Novo	6	32	3	3.56	2.49	0.93	0.57	0.49	0.53	-0.16
N = number of individuals; N_A = number of alleles; N_P = number of private alleles; Na = average number of alleles per locus; Ne = average number of effective alleles per locus; I = Shannon’s information index; Ho = mean observed heterozygosity; He = mean expected heterozygosity; uHe = mean unbiased expected heterozygosity; F = inbreeding coefficient (fixation index).

The sample from Kainji Lake revealed the highest genetic diversity among the five populations examined (Table 1): average number of alleles (Na) = 8.33; average number of effective alleles (Ne) = 4.56; Shannon’s information index (I) = 1.6; observed heterozygosity (Ho) = 0.73; expected heterozygosity (He) = 0.70; and unbiased expected heterozygosity (uHe) = 0.72. The highest number of private alleles was also observed at Kainji Lake: 21 out of 75 alleles in this locality. This sample had a low inbreeding coefficient (F_IS = –0.03). Genetic diversity in the samples from the three southern Nigerian populations were substantially lower: Na ranged between 4.89 (Ethiope River) and 5.67 (Igbokoda); Ne between 2.55 (Igbokoda) and 2.65 (Epe Lagoon); I between 0.97 (Epe Lagoon) and 1.11 (Igbokoda); Ho between 0.44 (Epe Lagoon) and 0.50 (Ethiope River); He between 0.47 (Epe Lagoon) and 0.54 (Igbokoda); uHe between 0.48 (Epe Lagoon) and 0.56 (Igbokoda). In these localities, F_IS ranged between –0.01 (Ethiope River) and 0.18 (Igbokoda). Excluding locus Hn11 in Igbokoda, for which null alleles were suggested in this locality, F_IS was 0.12, still the highest value. Thus, it is possible that the excess of homozygotes in Hn11 observed in Igbokoda results from inbreeding rather than null alleles. Between two and four private alleles were observed in the southern Nigerian populations (out of 44, 51, and 46 alleles observed in Ethiope River, Igbokoda, and Epe Lagoon, respectively). The sample from the Ouémé River in Porto Novo, Benin, had the lowest allelic diversity (3.56, 2.49 and 0.93 for Na, Ne and I, respectively); although Ho (0.57) was higher than the southern Nigerian localities, and He (0.49) and uHe (0.53) were within the range observed in these localities. Three private alleles were observed at this locality (out of 32). The lowest inbreeding coefficient was estimated for the sample of this locality (F_IS = -0.16). Estimates of genetic diversity per locus per population is presented in Table SM2.

Population genetic structure

F _ST values obtained with the program FreeNA were practically identical with and without correction for null alleles, indicating that null alleles, if present, did not bias results. The highest F_ST values were obtained for the pairwise comparisons that included Kainji Lake (Table 2). Average F_ST for these comparisons correcting for null alleles was 0.18. Average F_ST for pairwise comparisons between southern Nigerian localities correcting for null alleles was 0.05. Average F_ST for the pairwise comparisons between southern Nigeria localities and the Ouémé River in Benin was 0.11. All pairwise 95% CI excluded zero; thus, they are considered to be statistically significant. Similarly, all pairwise comparisons were significant according to the FDR test (FDR ≤ 0.05), using the p-values obtained for the pairwise comparisons without correcting for null alleles.

Table 2

Pairwise F_ST estimations with (above the diagonal) and without correction for null alleles (95% confidence intervals in brackets; p-values below brackets; * = <0.0001)
Population	Kainji Lake	Ethiope River	Igbokoda	Epe Lagoon	Porto Novo
Kainji Lake		0.17 [0.11–0.24]	0.15 [0.11–0.19]	0.21 [0.15–0.29]	0.19 [0.09–0.30]
Ethiope River	0.18 [0.11–0.24] *		0.03 [0.01–0.04]	0.08 [0.03–0.14]	0.09 [0.04–0.14]
Igbokoda	0.17 [0.12–0.21] *	0.02 [0.004–0.03] *		0.05 [0.02–0.08]	0.12 [0.06–0.17]
Epe Lagoon	0.22 [0.15–0.29] *	0.08 [0.03–0.14] *	0.04 [0.01–0.07] *		0.13 [0.06–0.23]
Porto Novo	0.19 [0.09–0.30] *	0.09 [0.04–0.14] 0.0001	0.12 [0.06–0.17] 0.0001	0.13 [0.06–0.23] *

Distances between clusters in DAPC plots reflect F_ST results (Figure 2), with the greatest distances observed between fish from Kainji Lake and fish from the southern localities. When the sample from Kainji Lake was removed, the greatest distances observed were between fish from the Ouémé River locality in Benin and fish from the three southern Nigeria localities. Fish from Epe Lagoon appear to be more genetically differentiated from conspecifics obtained from the other two southern Nigerian localities, and some overlap is observed between individuals from Igbokoda and Ethiope River. Analyses in STRUCTURE and PCA (Figures 3 and 4) also show a clear separation between fish from Kainji Lake and fish from the other localities; however, they do not show a separation among fish samples from the different southern localities. The STRUCTURE plot in Figure 3 includes two additional specimens from the Kainji Lake sample for which partial genotypes were obtained (one with four and one with five microsatellite-loci), and they grouped with all other specimens from Kainji Lake. The sample from Kainji Lake contained 3–4 individuals that were collected at Faku, just below the dam. For analyses excluding the Kainji Lake sample, STRUCTURE and PCA plots do not show evidence of genetic differentiation among fish from the southern localities. IBD analysis using F_ST values between populations produced an R² = 0.65 with a p-value = 0.049 for the analysis including all localities (Figure 5a), whereas the analysis excluding Kainji Lake resulted in an R² = 0.22 with a p-value = 0.21 (Figure 5b). IBD analyses using genetic differentiation between individuals resulted in an R² = 0.09 with a p-value < 0.0001 for the analysis including all localities (Figure 5c), whereas the analysis excluding individuals from Kainji Lake resulted in an R² = 0.005 with a p-value = 0.21 (Figure 5d).

Even though genetic evidence has long been recognized as crucial information for sustainable management of fisheries [11,38], such information for African inland fisheries remains very limited [12]. The present study advances understanding of population genetics of the African bonytongue, an important fish for commerce and subsistence throughout a large region of Africa. We documented patterns of genetic diversity and structure of African bonytongue in Nigeria, where > 85% of the global capture of this fish occurs, and one of the top seven countries driving growth of inland fisheries globally. The knowledge generated by the present study is important for guiding management and conservation strategies for this species.

The greatest genetic differentiation was observed between fish from Kainji Lake and those collected in the southern localities (average F_ST = 0.18). Although high, this level of differentiation is smaller than that observed (average F_ST = 0.27) between fish from Malanville, also in the Niger River, and southern Benin localities in Mono River and the Ouemé- Sô river floodplain system [20]. Differentiation between fish from Kainji Lake and those collected in the southern localities was supported by results from all genetic differentiation analyses, as well as the high number of private alleles (21) observed in Kainji Lake. The large geographic distance between Kainji Lake and southern populations is probably an important factor restricting gene flow of the African bonytongue between the two regions. A pattern of IBD was detected only when analyses included the sample of Kainji Lake. In addition, movements of the African bonytongue were affected by the construction of the Kainji Dam, completed in 1968. The dam constitutes an effective barrier for upstream movement of fishes from the river below, although fish can move downstream from the lake through the spillway. Nonetheless, the dam appears to have impacted downstream populations after its construction. Bonytongues were rarely caught at Faku, just below the dam, and Awuru, ~25 Km downstream from the dam, in 1973 and 1974 [39]. Our sample from Kainji Lake included a few specimens from Faku, and our genetic analyses grouped them in the same cluster as fish obtained from upstream of the lake, suggesting gene flow occurs from the lake to reaches just below the dam. The African bonytongue also became less frequent in landings in the second half of 1968 and in 1969 at Jebba and Pategi, ~100 Km and ~220 Km below the dam, respectively [40]. Between 1970 and 1972, however, bonytongues comprised a large percentage of the fish biomass landed at these locations, as well as at Lokoja, ~400 Km downstream the dam. This proportion dropped substantially in these three below-dam localities between 1974 and 1975 [39], which may be due to local exploitation. It is important to investigate genetic differentiation of African bonytongues at these below-dam localities, where this species has historically been a valuable fishery component.

Comparisons between microsatellite results for bonytongues from Malanville in Benin [20] and Kainji Lake (this study) suggest restricted gene flow between these two localities on the Niger River separated by a river distance of ~230 Km. Microsatellite datasets for these samples were not analyzed together because of potential bias when combining datasets from different studies [41]. Marked differences in mean expected heterozygosity (He) for these samples, however, suggest highly restricted gene flow between these two areas of the Niger River. He was 0.57 in Malanville (n = 12), whereas He was 0.73 in the sample from Kainji Lake (n = 23). Mean He is expected to be fairly insensitive to potential bias due to these differences in sample size [42,43], and can be useful for comparisons between different studies [44].

Among southern locations, gene flow appears to be more restricted between the lower Ouémé River in Benin (Porto Novo) and the southern Nigeria rivers sampled (average F_ST = 0.11). This level of differentiation is higher than the one reported in southern Benin between fish from Mono River and localities in the Ouemé- Sô river floodplain system (average F_ST = 0.09). Average divergence among southern Nigeria localities (average F_ST = 0.05) is higher than the reported between localities within the Ouemé- Sô river floodplain system (average F_ST = 0.03). A pattern of IBD was not detected among southern localities in this study. An ocean connection and the presence of brackish waters in Lagos Lagoon, which is connected to the western portion of Epe Lagoon, may represent effective barriers to gene flow for H. niloticus between Porto Novo and Epe Lagoon. According to our results, gene flow between upstream areas of the Ouémé River Basin and rivers connecting to Epe Lagoon also appears restricted. Significant pairwise F_ST values among the southern Nigerian localities suggest restricted gene flow among them, despite potential surface water connections provided by a complex network of waterways and annual flooding. Igbokoda River and Epe Lagoon are connected through a water channel that joins the Lekki Lagoon, which is connected to the eastern portion of Epe Lagoon. Thus, it appears that although movement of fish between the two localities may occur, it is not enough to prevent subtle genetic differentiation. Igbokoda and Ethiope rivers are separated by > 150 Km. Flooding and water channels may facilitate movement of fish between both localities, although the lower portion of the Igbokoda River has brackish water that may act as a barrier to dispersal. STRUCTURE analyses did not detect subdivision among fish from southern locations, but it has been shown that this method can fail to detect substructure in cases where F_ST values are low [45] or even relatively high. For example, in the previous study from Benin, most STRUCTURE analyses did not detect genetic differentiation between fish from the Mono River and those from localities in the Ouemé- Sô River floodplain system, even though all F_ST pairwise values were significant (average F_ST was 0.09) and other methods also clearly showed this differentiation [20].

The observed patterns of genetic differentiation indicate that African bonytongue from each of the localities examined in this study correspond to differentiated populations (i.e., genetic stocks), and should thus be treated separately for conservation and management. The Kainji Lake population is highly differentiated from the southern populations. In the south, the sample from Porto Novo is highly differentiated from southern Nigerian populations. Porto Novo is part of the Ouemé- Sô river floodplain system in Benin, for which Hurtado et al. [20] report low levels of genetic differentiation among fish collected from this system, separated by up to ~75 km (average F_ST = 0.03). In that study, the fish from the Ouemé- Sô river floodplain system were highly differentiated from those collected in the Mono River and Malanville, Niger River. The southern Nigerian populations sampled in this study show subtle but statistically significant differentiation; we thus recommend that they be managed as local stocks. These genetically distinct populations could constitute valuable genetic resources for future use in aquaculture.

Multiple activities threaten the sustainability of the African bonytongue stocks we have identified in Nigeria. Overfishing compromises the sustainability of fishery resources in Kainji Lake [46,47], and illegal fishing activities, including the use of prohibited gear (e.g. small mesh size nets and destructive fishing gear), fish poisoning, and explosives, have exacerbated overfishing in this lake [48]. Environmental pollution also impacts fish and people who consume them; and high levels of heavy metals have been detected in fish from the area [49]. The Ethiope River flows through a densely populated area of Nigeria’s Delta state, where pollution has impaired water quality in some stretches, with documented effects on macroinvertebrates [50]. Effects of degraded water quality in this river has not yet been shown to impair fish survival directly [51]. The Igbokoda River is an important fishing location in the Ilaje local government at Ondo State, a major oil-producing state. Unsafe levels of heavy metals have been found in water samples and fish in this river [52]. Consumption of fish is also considered unsafe in other rivers located in the Ilaje local government due to oil and industrial pollution [53]. Given the low quality of its riverine environment, and high levels of inbreeding detected in this study, urgent attention needs to be paid to the African bonytongue in Igbokoda River. A decade ago, Epe Lagoon was reported to have diverse and abundant fish stocks threatened by a rapidly growing population in the Lagos metropolitan area [54]. Anthropogenic activities appear to largely contribute to pollution in Epe Lagoon. Recent studies of Epe lagoon reported the presence of benzene, toluene, ethylbenzene and xylene (BTEX) and unsafe levels of Polycyclic Aromatic Hydrocarbons (PAH) in H. niloticus, high levels of heavy metals in sediment of the lake, and a high prevalence of an intestinal parasite in H. niloticus, which may be attributable to the water pollution [55,56].

Studies of large-scale patterns of genetic structure of other important capture fisheries are limited in West Africa. A recent study of the Nile tilapia Oreochromis niloticus analyzed samples collected from 23 localities across eight West African countries, representing the major catchments of the Volta, Niger, Senegal and Gambia River basins [57]. That study found a pattern of IBD among all localities and significant spatial genetic structure that largely corresponds to major river basins and, to a lesser extent, sub-basins. Within the Volta Basin, a significant, yet much weaker relationship between genetic and geographic distances was observed, suggesting IBD was a relatively minor factor shaping genetic differentiation amongst populations. This is similar to our findings for the African bonytongue, for which IBD was significant at a large scale (i.e., when Kainji Lake samples were included), but not at smaller scales (i.e., when only samples of southern populations were included). Most of the pairwise F_ST values in the Nile tilapia study were significant, with the exception of pairwise comparisons among localities separated by some of the shortest distances. Significant genetic differentiation in the Nile tilapia was usually observed between populations separated by more than ~90 Km. Interestingly, Nile tilapia samples from the only two localities sampled in the Niger River, Malanville and Mopti, which are separated by ~1,400 Km, show high genetic similarity, suggesting high levels of gene flow. This is in contrast to what we observed for the African bonytongue in the Niger River, for which high genetic differentiation appears to occur at comparatively shorter distances within this river, i.e., between Malanville and Kainji Lake (~230 Km), and between Kainji Lake and the lower Niger portion (~700 Km).

An effect of floodplain connectivity and geographic scale has also been reported in the African bonytongue’s closest living relative, Arapaima gigas [58]. This species, which is distributed throughout the Amazon River Basin, is the only other living member of the Arapaiminae [59], and both species construct nests where they lay and protect eggs. The predatory arapaima protects free-swimming larvae much longer (several weeks) than the omnivorous African bonytongue (a few days) and also grows much larger (over 2 m and 100 Kg compared to Heterotis at 1 m and 10 Kg). At a fine scale (e.g. within the same floodplain system; < 25 Km in Arapaima and ~ 75 Km in Heterotis), the two species tend to exhibit genetic homogeneity [58] [20]. At a meso scale (e.g. in separate floodplain systems; ~100 Km in Arapaima and 69–400 Km in Heterotis), both species exhibit low but significant values of genetic differentiation. Finally, at the largest scale (e.g. >1300 Km in Arapaima and > 510 Km in Heterotis), the highest levels of genetic differentiation were observed in both species. Thus, despite occupying two separate continents these two sister taxa exhibit similar patterns of genetic differentiation.

A broad range of genetic diversity values are observed in African bonytongue populations of Benin and Nigeria. Globally, an average of Na = 9.1 and heterozygosity H = 0.54 is reported from freshwater fishes [60]. The sample from Kainji Lake had the highest heterozygosity (Ho = 0.73; He = 0.70) and second highest allelic diversity (Na = 8.33) among all African bonytongue populations examined to date. The sample from the Oueme–Sô river-floodplain system in Benin has the highest reported Na (9.25), but with a much larger sample examined (n = 184), and second highest heterozygosity (Ho = 0.60; He = 0.69). Genetic diversity values for fish from the southern Nigeria localities were lower (Na range = 4.89–5.67; Ho range = 0.44–0.50; He range = 0.47–0.54). The Igbokoda River sample revealed highest levels of inbreeding (F_IS = 0.18), and genetic diversity of fish from Malanville, Benin, was very low (Na = 3.50; Ho = 0.34; He = 0.43; F_IS = 0.20). Low genetic diversity at Malanville could be associated with intense fishing pressure [20]. Although we observed low allelic diversity (Na = 3.56) for the sample from the Oueme–Sô River at Porto Novo, Benin, Ho (0.57) was more similar to the value previously reported for the Oueme–Sô river-floodplain, and He (0.49) was similar to that observed for fish from southern Nigeria. The smaller Na and He values obtained for Porto Novo compared to the Oueme–Sô river-floodplain system are likely due to small (n = 6) sample size [43]. More similar values of Na to the one obtained in Porto Novo, however, were obtained for localities within the Oueme–Sô river floodplain system [20]. For a sample of n = 6 in the Ouemé River channel, values of Na = 3.38, Ho = 0.48 and He = 0.66 were obtained; whereas for a sample of n = 10 in the Sô River channel, values were Na = 3.62, Ho = 0.55 and He = 0.61.

Consistent with a previous study conducted in Benin, we found significant genetic structure among African bonytongue samples from locations in Nigeria. The highest values of genetic differentiation were observed between Kainji Lake and the southern localities examined. We found evidence of gene flow from Kainji Lake to the area just below the dam, but the large distance between Kainji Lake and southern populations probably restricts gene flow between the two regions. Relatively low but significant F_ST values were obtained for pairwise comparisons among fish from the southern localities. Based on these results, we suggest that populations of the African bonytongue in each of the Nigerian localities examined in this study be treated separately for conservation and management. Given the observed levels of genetic differentiation, understanding patterns of genetic diversity in the African bonytongue throughout the rest of its range is needed for adequate delineation of management stocks. We recommend that populations located in different rivers that are not part of the same floodplain system and are separated by >70 Km, or populations within a major river separated by >200 Km, be treated as separate management units pending genetic evidence.

ACKNOWLEDGEMENTS

We would like to express our appreciation to Kika E. Ogbuebunu, who assisted TO in obtaining samples from Kainji Lake, Zannou Patrice, who helped AA in obtaining samples from Porto Novo, and Elizabeth Carrera for technical advice.

FUNDING

Financial support for fieldwork in Benin and reagents and genotypic expenses was provided to KW by the state of George and Carolyn Kelso via the International Sportfish Fund.

AUTHOR CONTRIBUTIONS

TEO, MOA conceived the study; TEO, AA obtained the samples; TEO conducted laboratory analyses; LAH supervised laboratory work and population genetics analyses; ICC and MM provided laboratory guidance; KOW, LAH, MM contributed reagents and genotyping expenses; ICC, LAH, TEO, MM performed data analyses, prepared figures, and interpreted results; LAH wrote original draft / led manuscript writing; TEO, KOW, MM, ICC critically reviewed and edited the manuscript; all authors reviewed and approved final draft.

Béné, C. & Heck, S. Fish and food security in Africa. NAGA, WorldFish Center Quarterly 28, 8-13 (2005).
Funge-Smith, S. J. Review of the state of world fishery resources: inland fisheries. FAO Fisheries and Aquaculture Circular (2018).
Funge-Smith, S. & Bennett, A. A fresh look at inland fisheries and their role in food security and livelihoods. Fish Fish. (Oxf.) 20, 1176-1195 (2019).
De Graaf, G. & Garibaldi, L. The value of African fisheries. FAO fisheries and aquaculture circular, I (2015).
FAO. FAO yearbook. Fishery and Aquaculture Statistics 2018/FAO annuaire. Statistiques des pêches et de l’aquaculture 2018/FAO anuario. Estadísticas de pesca y acuicultura 2018., (2020).
Olaosebikan, B. D. & Bankole, N. O. in In: Proceedings of the 19th annual conference of the fisheries society of Nigeria (FISON), 29 Nov - 03 Dec 2004 754– 762 (Ilorin, Nigeria, 2005).
Marshall, B. E. Inland fisheries of tropical Africa. J. F, Graig, Freshwater Fisheries Ecology (eds.). Chichester: John Wiley & Sons, Ltd, 349 (2016).
Dudgeon, D. et al. Freshwater biodiversity: importance, threats, status and conservation challenges. Biol. Rev. (Camb.) 81, 163-182 (2006).
United Nations-Department of Economic and Social Affairs-Population Division. World population prospects 2019: Highlights (st/esa/ser. A/423). (2019).
FAO, IFAD, UNICEF, WFP & WHO. The state of food security and nutrition in the world 2019: safeguarding against economic slowdowns and downturns 2019. (Rome, Italy: FAO, <http://www.fao.org/3/ca5162en/ca5162en.pdf> 2019).
Carvalho, G. R. & Hauser, L. Molecular genetics and the stock concept in fisheries, in Molecular Genetics in Fisheries (eds GR Carvalho & TJ Pitcher) 55-79 (Springer, 1995).
Abban, E. K. Considerations for the conservation of African fish genetic resources for their sustainable exploitation, in Towards policies for conservation and sustainable use of aquatic genetic resources. ICLARM Conf. Proc. 59, 277p. (eds R.S.V. Pullin, D.M. Bartley, & J. Kooiman) 95-100 (International Center for Living Aquatic Resources Management (ICLARM) and FAO).
FAO. Fishery Statistical Collections: Global Capture Production 1950-2018, <http://www.fao.org/fishery/statistics/global-capture-production/query/en> (2020).
Chan, C. Y. et al. Prospects and challenges of fish for food security in Africa. Glob. Food Sec. 20, 17-25 (2019).
Olopade, O. A., Taiwo, I. O. & Dienye, H. E. Management of Overfishing in the Inland Capture Fisheries in Nigeria. LimnoFish 3, 189-194 (2017).
Gbaguidi, A. S. & Pfeiffer, V. Stastistiques des peches continentals, Annees 1987–1995. Cotonou, Benin: GTZ-GmbH, Benin Direction des Pêches (1996).
Monentcham, S.-E., Kouam, J., Pouomogne, V. & Kestemont, P. Biology and prospect for aquaculture of African bonytongue, Heterotis niloticus (Cuvier, 1829): A review. Aquaculture 289, 191-198 (2009).
FAO. The State of the World’s Aquatic Genetic Resources for Food and Agriculture., (Rome, 2019).
Mustapha, M. K. Heterotis niloticus (Cuvier, 1829) a threatened fish species in Oyun reservoir, Offa, Nigeria; the need for its conservation. Asian J. Exp. Biol. Sci. 1, 1-7 (2010).
Hurtado, L. A., Carrera, E., Adite, A. & Winemiller, K. O. Genetic differentiation of a primitive teleost, the African bonytongue Heterotis niloticus, among river basins and within a floodplain river system in Benin, West Africa. J. Fish Biol. 83, 682-690 (2013).
Hauber, M. E., Bierbach, D. & Linsenmair, K. E. A description of teleost fish diversity in floodplain pools ('Whedos') and the Middle-Niger at Malanville (north-eastern Benin). J. Appl. Ichthyol. 27, 1095-1099 (2011).
Carrera, E., Renshaw, M. A., Winemiller, K. O. & Hurtado, L. A. Isolation and characterization of nuclear-encoded microsatellite DNA primers for the African bonytongue, Heterotis niloticus. Conserv. Genet. Resour. 3, 537-539 (2011).
Lischer, H. E. L. & Excoffier, L. PGDSpider: an automated data conversion tool for connecting population genetics and genomics programs. Bioinformatics 28, 298-299 (2012).
Raymond, M. & Rousset, F. GENEPOP (version-1.2) - Population genetics software for exact tests and ecumenicism. J. Hered. 86, 248-249 (1995).
Rousset, F. Genepop'007: a complete re-implementation of the genepop software for Windows and Linux. Mol. Ecol. Resour. 8, 103-106 (2008).
Benjamini, Y. & Hochberg, Y. Controlling the False Discovery Rate - a practical and powerful approach to multiple testing. J. R. Stat. Soc. Series B Stat. Methodol. 57, 289-300 (1995).
Peakall, R. & Smouse, P. E. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol. Ecol. Notes 6, 288-295, doi:10.1111/j.1471-8286.2005.01155.x (2006).
Peakall, R. & Smouse, P. E. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics 28, 2537-2539 (2012).
Van Oosterhout, C., Hutchinson, W. F., Wills, D. P. M. & Shipley, P. MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol. Ecol. Notes 4, 535-538 (2004).
Chapuis, M. P. & Estoup, A. Microsatellite null alleles and estimation of population differentiation. Mol. Biol. Evol. 24, 621-631 (2007).
Excoffier, L. & Lischer, H. E. L. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 10, 564-567, doi:10.1111/j.1755-0998.2010.02847.x (2010).
Jombart, T. adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics 24, 1403-1405 (2008).
Jombart, T. & Ahmed, I. adegenet 1.3-1: new tools for the analysis of genome-wide SNP data. Bioinformatics 27, 3070-3071 (2011).
Jombart, T., Devillard, S. & Balloux, F. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 11, 94 (2010).
Miller, J. M., Cullingham, C. I. & Peery, R. M. The influence of a priori grouping on inference of genetic clusters: simulation study and literature review of the DAPC method. Heredity, doi:10.1038/s41437-020-0348-2 (2020).
Pritchard, J. K., Stephens, M. & Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 155, 945-959 (2000).
Mantel, N. The detection of disease clustering and a Generalized Regression Approach. Cancer Res. 27, 209-220 (1967).
Allendorf, F. W., Ryman, N. & Utter, F. M. Genetics and fishery management: Past, present, and future, in Population Genetics and Fishery Management (eds N Ryman & F. Utter) 1-19 (Washington Sea Grant Publications/ University of Washington Press, 1987).
Otobo, F. O. The commercial fishery of the middle River Niger, Nigeria, in Symposium on River and Floodplain Fisheries in Africa, Bujumbura, Burundi, 21-23 November 1977, Review and Experience Papers Vol. CIFA TECHNICAL PAPER No. 5 (ed R. L. Welcomme) (Committe for Inland Fisheries of Africa, FAO, 1978).
Lelek, A. & El-Zarka, A. Ecological comparison of the preimpoundment and postimpoundment fish faunas of the River Niger and Kainji Lake, Nigeria. Geophys. Monogr. Ser. 17, 655-660 (1973).
Morin, P. A., Manaster, C., Mesnick, S. L. & Holland, R. Normalization and binning of historical and multi-source microsatellite data: overcoming the problems of allele size shift with allelogram. Mol. Ecol. Resour. 9, 1451-1455 (2009).
Pruett, C. L. & Winker, K. The effects of sample size on population genetic diversity estimates in song sparrows Melospiza melodia. J. Avian Biol. 39, 252-256, doi: 10.1111/j.2008.0908-8857.0409 (2008).
Hale, M. L., Burg, T. M. & Steeves, T. E. Sampling for microsatellite-based population genetic studies: 25 to 30 individuals per population is enough to accurately estimate allele frequencies. PLoS ONE 7, e45170 (2012).
Macedo, D. et al. Population genetics and historical demographic inferences of the blue crab Callinectes sapidus in the US based on microsatellites. PeerJ 7, e7780 (2019).
Latch, E. K., Dharmarajan, G., Glaubitz, J. C. & Rhodes, O. E. Relative performance of Bayesian clustering software for inferring population substructure and individual assignment at low levels of population differentiation. Conserv. Genet. 7, 295-302 (2006).
Abiodun, J. A. Fisheries Statistical Bulletin Kainji Lake, Nigeria, 2001. 25p (2002).
Yem, I. Y., Sani, A. O., Bankole, N. O., Onimisi, H. U. & Musa, Y. M. in 21st Annual Conference of the Fisheries Society of Nigeria (FISON) 79-85 (Calabar, Nigeria, 2007).
Mshelia, M. B. et al. in 19th Annual Conference of the Fisheries Society of Nigeria (FISON) 597-604 ( Ilorin, Nigeria, 2005).
Adelakun, K. M. & Kehinde, A. S. Heavy metals bioaccumulations in Chrysichthys nigrodigitatus (Silver catfish) from River Oli, Kainji Lake National Park, Nigeria. Egyptian Journal of Aquatic Biology and Fisheries 23, 253-259 (2019).
Ikomi, R. B. & Arimoro, F. O. Effects of recreational activities on the littoral macroinvertebrates of Ethiope River, Niger Delta, Nigeria. J. Aquat. Sci. 29, 155-170 (2014).
Ushurhe, O., Origho, T. & Ewhuwhe-Ezo, J. Determinant of water quality and suitability of River Ethiope for fish survival in Southern Nigeria. Can. J. Agr. Crop. 1, 11-18 (2016).
Arojojoye, O. A., Oyagbemi, A. A. & Afolabi, J. M. Toxicological assessment of heavy metal bioaccumulation and oxidative stress biomarkers in Clarias gariepinus from Igbokoda River of South Western Nigeria. Bull. Environ. Contam. Toxicol. 100, 765-771 (2018).
Arojojoye, O. A. et al. Assessment of water quality of selected rivers in the Niger Delta region of Nigeria using biomarkers in Clarias gariepinus. Environ. Sci. Pollut. Res. 28, 22936-22943 (2021).
Olukolajo, S., Olufemi & Hillary, E., Chikezie. Species diversity and growth pattern of the fish fauna of Epe Lagoon, Nigeria. J. Fish. Aquat. Sci. 7, 392-401 (2012).
Akinsanya, B., Ayanda, I. O., Fadipe, A. O., Onwuka, B. & Saliu, J. K. Heavy metals, parasitologic and oxidative stress biomarker investigations in Heterotis niloticus from Lekki Lagoon, Lagos, Nigeria. Toxicol. Rep. 7, 1075-1082 (2020).
Akinsanya, B., Ayanda, I. O., Onwuka, B. & Saliu, J. K. Bioaccumulation of BTEX and PAHs in Heterotis niloticus (Actinopterygii) from the Epe Lagoon, Lagos, Nigeria. Heliyon 6, e03272, doi:https://doi.org/10.1016/j.heliyon.2020.e03272 (2020).
Lind, C. E. et al. Genetic diversity of Nile tilapia (Oreochromis niloticus) throughout West Africa. Sci. Rep. 9, 1-12 (2019).
Araripe, J., do Rêgo, P. S., Queiroz, H., Sampaio, I. & Schneider, H. Dispersal capacity and genetic structure of Arapaima gigas on different geographic scales using microsatellite markers. PLoS ONE 8, e54470 (2013).
Hilton, E. J. & Lavoué, S. A review of the systematic biology of fossil and living bony-tongue fishes, Osteoglossomorpha (Actinopterygii: Teleostei). Neotrop. Ichthyol. 16 (2018).
DeWoody, J. A. & Avise, J. C. Microsatellite variation in marine, freshwater and anadromous fishes compared with other animals. J. Fish Biol. 56, 461-473 (2000).

No competing interests reported.

Download PDF

Editorial decision: Major revision
28 Mar, 2022
Reviews received at journal
15 Mar, 2022
Reviews received at journal
12 Dec, 2021
Reviewers agreed at journal
24 Nov, 2021
Reviewers invited by journal
24 Nov, 2021
Editor assigned by journal
24 Nov, 2021
Editor invited by journal
30 Sep, 2021
Submission checks completed at journal
30 Sep, 2021
First submitted to journal
27 Sep, 2021

You are reading this latest preprint version

Genetic Identification and Diversity of Stocks of The African Bonytongue, Heterotis Niloticus (Osteoglossiformes: Arapaiminae), in Nigeria, West Africa

Status:

Version 1

Abstract

Figures

Introduction

Methods

Sampling and DNA extractions

Microsatellites amplification and genotyping

Genetic diversity and population differentiation analyses

Results

Genetic diversity

Population genetic structure

Discussion

Conclusion

Declarations

Literature Cited

Additional Declarations

Supplementary Files

Status:

Version 1