Figures
Abstract
Cyphocharax magdalenae, a Colombian freshwater fish species, plays a vital role in nutrients distribution and serves as a significant food source for other fish species and local fishing communities. Considered a short-distance migratory species, C. magdalenae populations face substantial extinction risk due to human activities impacting their habitats. To address the lack of knowledge on genetic diversity and population structure, this study used next-generation sequencing technology to develop species-specific microsatellite loci and conducted a population genetics analysis of C. magdalenae in the middle and lower sections of the Cauca River, Colombia. Out of 30 pairs of microsatellite primers evaluated in 324 individuals, 14 loci were found to be polymorphic, at linkage equilibrium and, in at least one population, their genotypic frequencies were in Hardy-Weinberg equilibrium. Results showed high genetic diversity levels compared to other neotropical Characiformes, with inbreeding coefficients similar to those reported for phylogenetically related species. Moreover, C. magdalenae exhibits seasonal population structure (rainy-dry) consisting of two genetic stocks showing bottleneck signals and high effective population sizes. This information is essential for understanding the current species genetics and developing future management programs for this fishery resource.
Citation: Ochoa-Aristizábal AM, Márquez EJ (2024) Genetic insights into Cyphocharax magdalenae (Characiformes: Curimatidae): Microsatellite loci development and population analysis in the Cauca River, Colombia. PLoS ONE 19(4): e0302273. https://doi.org/10.1371/journal.pone.0302273
Editor: Timothy Omara, BOKU: Universitat fur Bodenkultur Wien, AUSTRIA
Received: November 9, 2023; Accepted: March 29, 2024; Published: April 16, 2024
Copyright: © 2024 Ochoa-Aristizábal, Márquez. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the manuscript and its Supporting information files.
Funding: This study was supported by a grant framed under the Project “Variabilidad genética de un banco de peces de los sectores medio y bajo del río Cauca” (CT-2019-000661, Empresas Públicas de Medellín and Universidad Nacional de Colombia, Sede Medellín). Funders do not play any role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
The field of ichthyology is marked by a comprehensive fish diversity examination, with 36,775 valid species cataloged to date, 18,688 of which thrive exclusively in freshwater environments [1]. Within this diverse assemblage, the order Characiformes, a significant facet of Ostariophysi, stands out with a notable count of 2,334 valid species confined to freshwater habitats [2]. This taxonomically rich group spans the continents of Africa, South America, Central America, and southern North America, concentrating approximately 90% of its diversity in the Neotropical region [2]. Within Characiformes, the Curimatidae family encompasses 9 genera and 121 valid species distributed in the Neotropics [1]. Cyphocharax Fowler 1906 genus is one of the most abundant, having 47 valid species [1] and wide distribution ranging from the Pacific rivers in southern Costa Rica to the La Plata River and various coastal drainages in central Argentina [3, 4]. It has been found that Atlantic basins have a much more diverse group expanding from the Orinoco River, Amazonas River, rivers from the Guayanas, and eastern Brazil; on the contrary, a less diverse group of species habits rivers from northern Peru to Costa Rica, comprising the tributaries of Magdalena and Maracaibo [5].
Cyphocharax magdalenae (Steindachner 1878) commonly known as “viejito” or “Yalúa” (Fig 1) distributes from the Coto River in southeastern Cosa Rica [6] and Panamá [7, 8], to the Maracaibo lake and its tributaries in Venezuela [9] and Colombia [10]. In Colombia, C. magdalenae is the only species of the genus reported for the Atrato, Sinú, Canalete and Magdalena-Cauca basin [10, 11]. The Cauca River, main tributary of the Magdalena-Cauca basin [12, 13], has great economic importance due to its influence on most of the prioritized productive chains in Colombia, which include small exploitations, extensive industries, monocultures, energy generation, mineral extraction, agriculture and fishing [12, 14, 15]. This has resulted in fish being one of the most threatened biological resources, with few or unknown local priorities on rivers [16].
Courtesy of Jorge E. García-Melo and Luis J. García-Melo, Proyecto CaVFish Colombia.
Cyphocharax magdalenae constitutes 30% of the total population within the lower Magdalena-Cauca floodplains and plays a crucial role in energy, nutrients, coal, and minerals transportation and sediment resuspension in the Cauca basin [15]. This potamodromous species undergoes short-distance reproductive migrations (<100 km) [17], spawns twice a year [11, 15], and lays small (0.73 ± 0.04 mm) adhesive eggs with small perivitelline space [11, 18]. Moreover, this fish has a generation time of 2.6 years [19], lacks parental care for embryos, and does not build nests.
Additionally, the fishing importance of C. magdalenae has been increasing during the last years due to the replacement of traditionally captured species such as catfish Prochilodus magdalenae Steindachner 1879 and Blanquillo Sorubim cuspicaudus Littmann, Burr & Nass 2000 [15, 20]. Although classified as a Least Concern species by The International Union for Conservation of Nature (IUCN), its notable overfishing-caused mortality suggests a possible decrease in its populations [11, 21, 22]. Moreover, C. magdalenae lacks genetic studies that allow knowing the current species status and delving into its ecological processes, which, in turn, hinders the species protection and conservation plans development and application.
The scarce information for C. magdalenae contrasts with population genetic studies performed for Curimata mivartii [23] and species with fishing relevance of the phylogenetically close family Prochilodontidae like Prochilodus magdalenae [24–26], Prochilodus reticulatus [27] and Ichthyoelephas longirostris [28]. Likewise, it contrasts with studies on species of other families inhabiting the Magdalena-Cauca basin such as Pseudoplatystoma magdaleniatum [29, 30], Brycon henni [31–33], S. cuspicaudus, Ageneiosus pardalis, Pimelodus grosskopfii [34], Pimelodus yuma [35], Pseudopimelodus magnus and Pseudopimelodus atricaudus [36]. Among the different markers used in the cited studies, microsatellites stand out as the most used tool to solve questions at micro-evolutive scale in population genetics. In this area, the Single Nucleotide Polymorphisms have been utilized for species delimitation and for exploring genetic introgression in Pimelodus species [37].
To estimate the species genetic status in the lower section of the Cauca River, this study contrasted three hypotheses: (i) C. magdalenae shows signs of genetic deterioration with low genetic variability levels, high inbreeding levels and recent bottlenecks resulting from habitat alterations due to different anthropogenic activities; (ii) the species exhibit high gene flow as a result of the scarce or non-existent geographic barriers in this habitat, which is based on evidence found for other Curimatidae species [23]; and (iii) seasons influencing fish migration in the Cauca River shape the population structure of C. magdalenae. To test these hypotheses, this study identified de novo and characterized species-specific microsatellite loci to C. magdalenae. This approach was chosen instead of using heterologous loci to avoid potential pitfalls associated with cross-amplification [38], such as allele size homoplasy [39], unsuccessful amplification in phylogenetically distant species [40], low polymorphism, presence of null alleles [41], and amplification of non-orthologous loci [42].
Materials and methods
Biological material and study area
This study analyzed a total of 324 muscle tissues and/or fins of C. magdalenae preserved in ethanol 70%, provided between 2019 and 2021 by Universidad de Antioquia, Universidad de Córdoba, and Universidad Nacional de Colombia Sede Medellín, through scientific cooperation agreement CT-2019-000661, under environmental license # 0155 of January 30th, 2009, from Ministry of Environment, Housing and Territorial Development for the Ituango hydroelectric construction.
Samples from individuals were collected in 19 sites of the main channel, some swamps, and tributaries of the lower sections of the Cauca River. This zone, characterized by a large alluvial plain surrounded by mountains and flat and wavy surfaces where various swampy complexes are formed [43, 44], has been exposed to different anthropic activities including stockbreeding, fishing and extensive cultures [see 23]. Additionally, high mercury and sediment levels and elevated vegetation amounts are reported in the main courses of the river as a result of deforestation and flooding that may affect ichthyofauna [45]. According to their distribution, samples were grouped in five sections (S4-S8, Fig 2) previously described by Landínez-García & Márquez (2016) [28].
Self-made creation of the map based on contour lines scaled 1:100,000 from the Instituto Geográfico Agustín Codazzi source, 2019 (Available from: IGAC Geoportal, https://geoportal.igac.gov.co/contenido/datos-abiertos-cartografia-y-geografia).
Microsatellite loci identification and development of primers
Methodology previously described [23, 28] was followed for developing species-specific microsatellite loci for C. magdalenae. The genomic library was built from total DNA extraction from the fin (5,775 μg) of a C. magdalenae individual captured in the Cauca River (CYPHO11886) using the DNA extraction reagents and recommendations of the manufacturer PureLink genomic DNA Mini Kit (Invitrogen), then it was used for NGS surface sequencing (Whole Genome Sequencing), through Illumina Miseq (300 PE). After raw reads cleaning using PRINSEQ lite (removal of adapter sequences utilizing <Q30 quality bases), 100,000 high-quality extended reads were analyzed with PAL_FINDER v0.02.03 [46] for extracting reads containing tri-, tetra- and pentanucleotide microsatellite motifs. Then, Primer3 v2.0 [47] was used for the primer designs in the flanking sequences of the microsatellite loci and ultimately the correct evaluation of the primers was tested with ePCR [48]. Selection of the microsatellite loci set was performed according to already stablished features for validating new primers of microsatellite loci [49, 50].
For genotyping, 10 μl reaction mixtures were used with final concentrations of 0.3 pmol/μl of each forward primer tagged on the 5’ end with one of four adapters (tails A, B, C and D; [51]), 6 pmol/μl of each reverse primer, 0.5 pmol/μl fluorescently labeled adapter (6-FAM, VIC, NED and PET, Applied Biosystems), 1X Master Mix, 2.5% v/v GC Enhancer Platinum Multiplex PCR Master Mix (Applied Biosystems) and 3–5 μg/μl DNA isolated using the Purelink® purification kit (Thermo Scientific) with a modification in digestion time of 24 h when working with fins. Thermal profiles included an initial denaturalization step at 90°C for 35 s and an annealing step at 56°C for 35 s (with no final elongation). Subsequently, amplicons were separated by electrophoresis on an ABI 3730 XL automated sequencer using 600 LIZ as the internal molecular size. Then, GeneMarker® v3.0.0 was employed to denote allele fragments according to their molecular size (100 to 450 pb) and Micro-Checker v.2.2.3 [52] to detect potential genotyping errors.
Polymorphism, genetic diversity, outlier loci detection, and demographical events
Polymorphism information content (PIC) and average number of alleles per locus for each microsatellite marker were calculated using Cervus v3.0.7 [53] and GenAlEx v6.5.03 [54]. Estimation of observed (HO) and expected (HE) heterozygosities, inbreeding coefficient (FIS) and tests for departures from Hardy-Weinberg equilibrium were performed using Arlequin v3.5.2.2 [55], applying the sequential Bonferroni correction for the statistical significance in the multiple comparisons [56, 57].
For detecting outlier loci and determining evolutive forces acting over the microsatellite loci, BayeScan v2.1 [58] was used for performing an analysis employing the parameters for prior odds of 10:1 for the neutral model, 20 pilot runs, each having 5,000 iterations, followed by 100,000 iterations and a burn-in of 50,000. Critical values were based on the posterior probability of the Bayes factor, following the Jeffrey’s scale [59], which sets a probability of 0.76 as substantial evidence for selection.
To determine recent genetic bottlenecks, excess heterozygosity was tested under the mutation drift equilibrium assumption in three mutational models of the microsatellite loci (IAM: infinite alleles model, SMM: stepwise mutational model, TPM: two-phase model) through Wilcoxon signed rank test [60] included in Bottleneck v.1.2.02 [61]. Moreover, M ratio (mean ratio of the number of alleles to allele size range); [62] included in Arlequin v3.5.2.2 [55] was calculated following the criteria according to which values lower than 0.68 indicate recent and severe reductions in population size [62]. Furthermore, effective population size (Ne) of the species was determined in each of the assessed populations using the linkage disequilibrium method and a minimum allele frequency of 0.02 implemented in NeEstimator v2.1 [63]. Change in HE of the species in 10, 50 and 100 generations (t) was estimated based on the equation [64]. Genetic deterioration grade was stablished using the critical values of 25% Ht reduction in 10 (critically in danger), 50 (in danger) and 100 (vulnerable) generations [65].
Genetic structure
Bayesian analysis of population partitioning in STRUCTURE v.2.3.4 [66] was performed to determine samples grouping according to their co-ancestry coefficient. Parameters included 1,000,000 steps of the Markov chain Monte Carlo; 100,000 iterations were regarded as burn-in to estimate each K value (1–8). Each analysis was repeated 20 times. For better estimating the genetic stocks (K), STRUCTURESELECTOR [67] was used for calculating the ΔK ad hoc statistic [68], estimators MEDMEANK, MAXMEANK, MEDMEDK and MAXMEDK [69], and to graphically represent the results using the integrated CLUMPAK software [70].
Genetic differentiation of C. magdalenae among sections, genetic groups suggested by STRUCTURE, and seasons (rainy and dry) was calculated using the standardized statistics FST[71–73] and Jost’s DEST [74, 75], and Analysis of Molecular Variance (AMOVA) [73] using GenAlEx v6.503 software [54]. Additionally, a Discriminant Analysis of Principal Components (DAPC) was performed in polymorphic loci genotype resulting in 324 analyzed individuals using the R- package Adegenet [76]. Moreover, allele frequency distributions in genetic stocks were assessed through a G-test in Genepop v. 4.7.5. software [77].
Results
Identification of microsatellite loci and detection of outlier loci
Low-coverage genome sequencing raw data provided a total of 3,224,415,544 base pairs (bp), 10,712,344 reads, 41.28%GC and 83.72% Q30. After depuration, 100,000 reads (42,621,284 bp) were obtained, 37,119 out of which had compound repeat motifs, 1,054 interrupted repeat motifs and 71,376 contained microsatellite loci with direct repetitions (53,649 potentially amplifiable loci). Furthermore, 30,120 were validated using electronic PCR.
Out of 30 tested loci, 16 showed clearly defined peaks and absence of stutter bands in the electropherograms; these microsatellite loci include tri- (3) and tetra- (13) nucleotides motifs. Micro-Checker analysis did not indicate stutter-associated genotyping errors nor allele loss in the loci (dropout), and linkage disequilibrium tests were not significant. Departures of allele frequencies from Hardy-Weinberg equilibrium and heterozygote deficits were identified (Table 1). In general, microsatellite loci exhibited allele size between 106 and 414 bp with PIC values ranging from 0.733 to 0.954 and HO and HE between 0.408–0.916 and 0.762–0.957, respectively (Table 1); no outlier loci were detected by Bayescan. Nonetheless, null alleles were found in loci Cym10 and Cym16; thus, they were excluded in subsequent analyses.
Genetic diversity and demographic events
Genetic diversity estimators (Table 2) showed that S7 had the lower average number of alleles per locus (15.714 alleles/locus) while S5 exhibited the higher mean number (18.857 alleles/locus). Mean HO was lower in S7 (0.772) and higher in S8 (0.827) and mean HE displayed the lower value in S5 (0.880) and the higher in S4 (0.887). Only two (Cym8 and Cym14) out of the 14 selected loci departed from the Hardy-Weinberg equilibrium in four of the five evaluated sections and the remaining fulfilled the assumption in at least two sections. A similar result of genetic diversity was observed in the genetic stocks (Stock1 and Stock2, Table 2) found by the genetic structure Bayesian analysis, as shown below.
Moreover, the inbreeding coefficients across loci in each population showed statistical significance with greater impact in population S7 (FIS = 0.120, P = 0.000). These values, although remained significant, significantly decreased in the stocks (FISStock1 = 0.051, P = 0.000; FISStock2 = 0.071, P = 0.000).
It was found that assessed populations recently suffered a drastic reduction in population size since both the modified Garza-Williamson index (M-ratio: 0.238–0.256) and P values of the Wilcoxon signed rank test for the IAM and the TPM showed statistical significance (Table 3). Results for SMM, however, were not significant.
Statistical significance is noted in bold (P<0.05).
Furthermore, the effective size calculation in each sampled section showed the lower values in S4 and S7 and the higher in S6 and S8; nevertheless, it was not possible to obtain a value for S5, and estimation was obtained only for one of the genetic stocks (Stock1) (Table 3). Due to the latter, in subsequent analyses the lower limits of confidence intervals were used for S5 and Stock2.
For all evaluated generations, the species is classified as non-threatened as no reduction percentage reached 25% (critical value of the classification) (Table 4). S8 showed the lower reduction percentage in all generations (0.002, 0.011, 0.021); nevertheless, S7 and S4 exhibited the greater loss of heterozygosity percentage in 100 generations with 20.034 and 14.645, respectively. Moreover, in the evaluation on genetic stocks, no approximation surpassed the 5% reduction in all generations; in fact, similar values were found in each stock for each generation, although Stock1 showed slightly higher values.
Genetic structure
It was found that the most likely number of populations (ΔK) based on the Bayesian inference approach was 2. These two groups coexist in all sampled sections as confirmed by the co-ancestry histogram (Fig 3A). Additionally, it was observed that Stock1 predominates in sections S4 and S5 (pink), while Stock2 (purple) predominates in sections S6, S7 and S8. This was corroborated by the DAPC, showing that the two genetic stocks were not homogeneously distributed (Fig 3B and 3C); only just one similarity was found in the distribution of the tested individuals in sections S6 and S7 (Fig 3B). Differences in the frequencies of the stocks agree with results obtained in AMOVA (F´ST (4, 647) = 0.003, P = 0.001) and the pairwise comparisons of the standardized statistics F’ST, DEST (Table 5), which showed significant differences between S8 in relation to the remaining sections, and S6 in relation to S4 and S5.
Jost’s DEST (above diagonal) and F’ST (below diagonal).
Genetic structure comparison by collect season (rainy: April-June, October-December; dry: January-March, July-August; Serna et al., 2021) and sampling years, in addition to confirming the coexistence of two genetic stocks within the evaluated sections of the Cauca River (ΔK = 2, Mean LnP(K) = -22763.31000, MEDMEANK, MAXMEANK, MEDMEDK = 2; Fig 4A), revealed that during rainy seasons, the capture percentage was greater for Stock2 during the sampled years (2020: 91.7%; 2021: 94.3%), while Stock1 was primarily captured during dry season (2020: 67%; 2021: 58.8%) in comparison to the moderated captures of Stock2 (2020: 33%; 2021: 41.2%) (Fig 4B).
Rainy season (H): H_20: October 2020; H_21(1): April 2021; H_21(2): December 2021; Dry season (S): S_19: July-August-September 2019; S_20: February 2020; S_21: January 2021.
Discussion
To test three hypotheses related to the population genetics of C. magdalenae, this study developed a set of 16 polymorphic microsatellite loci capable of detecting diversity levels and genetic structure in the lower section of the Cauca River. The selected microsatellite loci were preferably those having tri- and tetranucleotide repetition motifs as these repetition patterns are the most recommended for their simplicity for genotyping and allele classification [46, 78, 79]. All microsatellite loci exhibited PIC values that allow describing them as highly informative according to Bostein et al., 1980 [80]; moreover, these values were higher than those reported for C. mivartii (0.549–0.946), I. longirostris (0.375–0.871) and P. magdalenae (0.399–0.949) [23, 26, 28].
Results obtained do not support the hypothesis of low genetic variability levels in C. magdalenae as it showed high alleles/locus and HE mean values in relation to those reported for neotropical Characiformes (Na: 10.92, HE: 0.675) [81]. As this is the first study performed for Cyphocharax genus, high diversity was found in comparison to C. mivartii (Na: 10.493, HO: 0.757, HE: 0.801), a Curimatidae species with which it shares its habitat [23]. Likewise, C. magdalenae displayed greater values over the two estimators with respect to I. longirostris (Na: 8.84–11. 05, HO: 0.701–0.767, HE: 0.771–0.798; [28]) and similar values to those reported for P. magdalenae (Na: 22.455, HO: 0.725, HE: 0.898; [26]) in the Cauca River. Even more, C. magdalenae showed greater values in alleles/locus and similar HE values compared to P. magdalenae, in rivers with great importance for Colombia, such as San Jorge (Na: 13.545, HO: 0.809; HE: 0.884), Sinú (Na: 15.273, HO: 0.767, HE: 0.882), Atrato (Na: 14.636, HO: 0.718, HE: 0.879), Nare (Na: 15.636, HO: 0.659; HE: 0.876), and Magdalena (Na: 19.455, HO: 0.758; HE: 0.896). C. magdalenae also exhibited much greater HE values than P. reticulatus (HE: 0.400) even though this estimator was calculated based on a different molecular marker [27].
Furthermore, the hypothesis of high inbreeding levels in C. magdalenae was not confirmed since the species showed inbreeding signals below 10% in each genetic stock (FIS Stock1 = 0.051; FIS Stock2 = 0.071). Although these values do not exceed the proposed limit [82, 83], it has been indicated that any inbreeding coefficient greater than cero has unfavorable effects on the fitness [84]. Considering these results and those reported for Characiformes in the Magdalena-Cauca basin using species specific microsatellites, C. magdalenae exhibited inbreeding levels relatively similar to those reported for C. mivartii (0.040–0.087; [23]), lower values to those described for P. magdalenae (0.125–0.255; [26]), and greater to those found in Brycon henni (-0.040 –-0.009; [33]), a phylogenetically farther species.
Additionally, this study corroborated the hypothesis of high gene flow despite having found significant differences among individuals of the lower section of the Cauca River as those differences are related to the coexistence of two genetic stocks with uneven distribution through the river. This indicates that there is no evidence of barriers separating the gene flow of C. magdalenae in the lower section of the Cauca River, which is consistent with data reported for C. mivartii [23], P. magdalenae [26], Cynopotamus magdalenae (Steindachner 1879), Megaleporinus muyscorum [85] and catfishes A. pardalis, P. grosskopfii, S. cuspicaudus [86], P. magdaleniatum [30], P. atricaudus and P. magnus [36]. Nonetheless, genetic structure resulting from seasonal variation was found since comparisons over the samples in rainy and dry seasons during the assessed years allowed detecting that one stock predominates in the dry season (Stock1) while the other predominates in the rainy season (Stock2). Sampling in both seasons for two years reflects an alternated cyclical behavior between both stocks and greater proliferative advantage of Stock2 over Stock1, which suggests a differential reproductive success leading to different cohorts (temporal Wahlund effect). These results are similar to those reported for Prochilodus lineatus (Valenciennes 1837) in the Mogi-Guaçu River in Brazil as individuals analyzed between 2005 and 2006 in rainy and dry seasons (January and August, respectively) showed evidence of temporal genetic structure [87]. Furthermore, it was found that Salminus brasiliensis (Cuvier 1816) in the Uruguay River in Brazil is formed by three genetic stocks resulting from a temporal genetic structure in which seasonal precipitations in the river are one of the main factors that could have originated genetic groups as the water level is fundamental in reproductive migration [88].
Furthermore, this study supports the hypothesis of a recent genetic bottleneck in C. magdalenae. These results may be a consequence of the joint action of various anthropogenic activities impacting the environment, directly affecting C. magdalenae individuals, or climatic changes. For instance, for The IUCN, C. magdalenae shows a decrease in its populations due to fishing pressure [11, 21, 22], which is consistent with the decline of this species reported in the Magdalena basin, with a diminishment in 887.7 tones unloaded for 2019 [21] in comparison to 436.1 tons for 2021 [22]. Moreover, there is no current regulation on a recommended or suggested minimum catch size for this species [89]; in fact, it was found that between 2018–2021 catch size ranges were lower than those described for size at maturity [22].
As it has been proposed for other species [86], contamination with mining-derived products could be affecting the viability of C. magdalenae, which is supported by evidence of methylmercury in fishes collected in the zone [45]. Another threat is the establishment of introduced exotic species [11, 90], such as Basa fish Pangasianodon hypophthalmus, a omnivorous species of fast size and weight growth, and having extension in Colombia, mainly in the Magdalena-Cauca basin [90], and the Nile Tilapia Oreochromis niloticus, which has the capacity to survive in all types of habitats, from salty, marine waters to estuarine and continental waters with temperatures between 8 and 42°C [90]. Despite the introduction of O. niloticus in Colombia was done with commercial purposes and control of Oreochromis mossambicus, an aggressive behavior against native species P. magdalenae has been displayed [91] and, after its introduction, there has been a reduction in the abundance of Triportheus magdalenae (Steindachner 1878) in the Guajaro reservoir in Magdalena [92].
In addition to the anthropogenic activities, it has been documented that the impact of climate factors has provoked changes in fish populations as they are responsible for abrupt changes in environments such as high temperatures, alterations of hydrological cycles, dissolved oxygen reduction, changes in mortality rates, growth, reproduction and distribution of fish populations [93]. This climate factor has also been used for explaining population genetics aspects of C. mivartii [23] and P. magdalenae [26].
Contrary to the bottleneck evidence, estimations performed over Ne (Stock1: 1464.6, Stock2: 2127) were higher than the critical value (Ne ≤ 1000), which allows deducting that the species has long term maintenance [84]. In conjunction with the previously mentioned high genetic diversity levels, this result suggests that C. magdalenae populations are large enough to sustain elevated genetic diversity and Ne, despite exhibiting signals of bottleneck. This idea is also supported by the results of the conservation status test as even in 100 generations (260 years: species with generational time of 2.6 years) [19] estimations did not show evidence of heterozygosity reductions greater than 25% (critical value of the classification), for which the species can be listed as non-threatened. This classification, however, should be interpreted carefully since it has been indicated that between 4–10% of freshwater species in South America are exposed to some extinction risk, mainly due to habitat loss or degradation [94].
To conclude, this study allowed determining that C. magdalenae in the lower section of the Cauca River exhibits seasonal population structure formed by two genetic stocks associated to the rainy and dry seasons, and which show high genetic diversity, low inbreeding coefficients, bottleneck signals and large Ne. Data obtained along with the microsatellite loci developed de novo in this study are a starting point for future research directed to monitoring the genetic diversity and population structure of this species to develop proper management plans.
Supporting information
S1 Table. Genotype data at 14 microsatellite loci included in the population genetic analysis of Cyphocharax magdalenae.
First row indicates respectively: Number of loci, Number of individuals, Number of sampling sections, Sample size for five sections. Third row indicates respectively: Sample ID, Sampling section, Genetic stock, Locus name for 14 loci.
https://doi.org/10.1371/journal.pone.0302273.s001
(XLSX)
S2 Table. Microsatellite loci not selected due to pitfalls in amplification or low levels of polymorphisms.
Na number of alleles per locus; HO Observed heterozygosity; HE Expected heterozygosity.
https://doi.org/10.1371/journal.pone.0302273.s002
(XLSX)
Acknowledgments
The authors express their gratitude to Universidad de Antioquia, Universidad de Córdoba, and Universidad Nacional de Colombia Sede Medellín for providing the fish tissue samples, and to the Proyecto CaVFish Colombia—Catálogo Visual de Peces de Agua Dulce de Colombia (cavfish.unibague.edu.co)—for providing the photograph of Fig 1.
References
- 1.
Fricke R, Eschmeyer WN. Fong JD. 2023 Species by family/subfamily. 2024 Feb 03. http://researcharchive.calacademy.org/research/ichthyology/catalog/SpeciesByFamily.asp.
- 2. Toledo-Piza M, Baena EG, Dagosta FCP, Menezes DA, Andrade M, Benine RC, et al. Checklist of the species of the Order Characiformes (Teleostei: Ostariophysi). Neotrop Ichthyol. 2024;22(1): e230086.
- 3. Maldonado-Ocampo JA, Vari RP, Usma JS. Checklist of the Freshwater Fishes of Colombia. Biota Colomb 9. 2008;9: 1–12.
- 4. Melo BF, Sidlauskas BL, Hoekzema K, Vari RP, Dillman CB, Oliveira C. Molecular phylogenetics of Neotropical detritivorous fishes of the family Curimatidae (Teleostei: Characiformes). Mol Phylogenet Evol. 2018;127: 800–812. pmid:29935301
- 5. Vari RP, Zanata AM, Camelier P. New species of Cyphocharax (Ostariophysi: Characiformes: Curimatidae) from the Rio de Contas Drainage, Bahia, Brazil. Copeia. 2010; 382–387.
- 6. Angulo A, Garita-Alvarado CA, Bussing WA, López MI. Annotated checklist of the freshwater fishes of continental and insular Costa Rica: Additions and nomenclatural revisions. Check List. 2013;9: 987–1019.
- 7. Smith SA, Bell G, Bermingham E. Cross-Cordillera exchange mediated by the Panama Canal increased the species richness of local freshwater fish assemblages. Proc R Soc B Biol Sci. 2004;271: 1889–1896. pmid:15347510
- 8. Vega AJ, Robles YA, Tuñón O, Barrera C. Fauna acuática del área centro occidental de Panamá. Tecnociencia. 2006;8: 87–100.
- 9.
Vari RP. Systematics of the Neotropical Genus Cyphocharax Fowler (Pisces: Ostariophysi). Smithsonian Contributions to Zoology. Washington, D.C; 1992. http://www.bioone.org/doi/abs/10.2994/SAJH-D-13-00022.1
- 10. DoNascimiento C, Agudelo-Zamora HD, Bogotá-Gregory JD, Méndez-López A, Ortega-Lara A, Lasso C, et al. Lista de especies de peces de agua dulce de Colombia / Checklist of the freshwater fishes of Colombia. v2.16. Asociación Colombiana de Ictiólogos. Dataset/Checklist. 2024. Available from:
- 11.
Lasso CA, Gutiérrez F de P, Morales-Betancourt MA, Agudelo E, Ramírez-Gil H, Ajiaco-Martínez RE, editors. II. Pesquerías continentales de Colombia: cuencas del Magdalena-Cauca, Sinú, Canalete, Atrato, Orinoco. Bogotá: Instituto de Investigación de los Recursos Biológicos Alexander von Humboldt; 2011. http://hdl.handle.net/20.500.12324/19395
- 12. González-López N, Carvajal-Escobar Y. Characterization of hydrological drought in the Cauca River high valley. Tecnol. cienc. agua. 2020;11: 235–264.
- 13. Cerón WL, Kayano MT, Ocampo-Marulanda C, Canchala T, Rivera IA, Ávila-Díaz A, et al. Spatio-temporal variability of hydroclimatology in the upper Cauca River basin in southwestern Colombia: Pre-and post-Salvajina dam perspective. Atmosphere (Basel). 2021;12.
- 14. Carvajal Y, Enciso A, Sandoval M. Hydrological analysis of historical floods in the upper valley of Cauca River. Ing. compet. 2016;18: 46.
- 15.
Jiménez-Segura L, Lasso CA, editors. XIX. Peces de la cuenca del río Magdalena, Colombia: diversidad, conservación y uso sostenible. Serie Editorial Recursos Hidrobiológicos y Pesqueros Continentales de Colombia. Instituto de Investigación de Recursos Biológicos Alexander von Humboldt; 2020. pp. 115–157. 2020. https://doi.org/10.21068/a2020rrhhxix
- 16. Hernández-Barrero S, Barco MV, Reyes CGB, Paramo J, Sierra LS, Stotz W. Is traditional fisheries management correctly addressing the possible causes of fish production decline? The relationship between environmental degradation and artisanal river fisheries in the Magdalena River basin, Colombia. Mar Freshw Res, 2022; 73(12):1475–1488.
- 17.
Usma Oviedo JS, Villa-Navarro F, Lasso CA, Castro F, Zúñiga-Upegui PT, Cipamocha CA, et al. Peces dulceacuícolas Migratorios de Colombia. In: Zapata LA, Usma JS. Guía de las especies migratorias de la biodiversidad en Colombia. Peces. Vol. 2. Bogotá: Ministerio de Ambiente y Desarrollo Sostenible, WWF-Colombia; 2013. p.215–441.
- 18.
Atencio-García VJ, Montes-Petro CD, Tapia-Pacheco CJ, Madariaga-Mendoza DL, Posada-González M, Ramírez-Guerra A, et al. Desarrollo embrionario y larvario de viejito Cyphocarax magdalenae. XVI Congreso Colombiano de Ictiología. ACICTIOS. 2022; Barranquilla, 27–29 de julio. https://acictios.org/congreso-colombiano-de-ictiologia-2022-memorias/
- 19.
Puentes V, Escobar FD, Polo CJ, Alonso JC, editors. Estados de los principales recursos pesqueros de Colombia 2014. Oficina de Generación del Conocimiento y la Información, Autoridad Nacional de Acuicultura y Pesca–AUNAP; 2014. http://hdl.handle.net/20.500.12324/34405
- 20. Blanco H V, Solipá J, Olaya-Nieto CW, Segura-Guevara FF, Brú-Cordero SB, Tordecilla-Petro G. Crecimiento y mortalidad de la Yalúa (Cyphocharax magdalenae Steindachner, 1878) en el río Sinú, Colombia. Rev MVZ Córdoba. 2005;10: 555–563. http://www.scielo.org.co/scielo.php?script=sci_arttext&pid=S0122-02682005000100005&lang=es%0Ahttp://www.scielo.org.co/pdf/mvz/v10n1/v10n1a05.pdf
- 21.
Duarte LO, Manjarrés–Martínez L, Reyes-Ardila H. Estadísticas de desembarco y esfuerzo de las pesquerías artesanales e industriales de Colombia entre febrero y diciembre de 2019. Autoridad Nacional de Acuicultura y Pesca. Bogotá; 2019. http://sepec.aunap.gov.co/Archivos/Boletines-2019/Boletin_Pesca_Consumo_artesanal_industrial_2019.pdf
- 22.
Duarte LO, Cuervo C, Vargas O, Gil-Manrique B, Tejeda K, De León G, et al. Estadísticas de desembarco y esfuerzo de las pesquerías artesanales de Colombia 2021. Santa Marta; 2021. http://sepec.aunap.gov.co/Archivos/Boletines-2021/SEPEC_Boletin_Desembarco_Artesanal_2021.pdf
- 23. Landínez-García RM, Marquez EJ. Microsatellite loci development and population genetics in neotropical fish Curimata mivartii (Characiformes: Curimatidae). PeerJ. 2018; 12–17. pmid:30479905
- 24. Aguirre-Pabón J, Narváez Barandica J, Castro García L. Mitochondrial DNA variation of the bocachico Prochilodus magdalenae (Characiformes, Prochilodontidae) in the Magdalena River Basin, Colombia. Aquat Conserv. 2013;23: 594–605.
- 25. Berdugo GO, Barandica JCN. Genetic diversity and population structure of bocachico Prochilodus magdalenae (Pisces, Prochilodontidae) in the Magdalena river basin and its tributaries, Colombia. Genet Mol Biol. 2014;37: 37–45. pmid:24688289
- 26. Landínez-García RM, Narváez JC, Márquez EJ. Population genetics of the freshwater fish Prochilodus magdalenae (Characiformes: Prochilodontidae), using species-specific microsatellite loci. PeerJ. 2020; 1–26. pmid:33240645
- 27. Macias JNL, Vallejo FG, Rincón ER, Giraldo AC, Cerón F. Diversidad genética del bocachico (Prochilodus reticulatus) de la cuenca alta del río Cauca (Colombia). Acta Biol Parana. 2009;38.
- 28. Landínez-García RM, Márquez EJ. Development and characterization of 24 polymorphic microsatellite loci for the freshwater fish Ichthyoelephas longirostris (Characiformes: Prochilodontidae). PeerJ. 2016;4: e2419. pmid:27635363
- 29. Gallo H, Díaz-sarmiento J. Variabilidad genética del bagre rayado Pseudoplatystoma fasciatum (Pisces: Pimelodidae) en el río magdalena. Rev Acad Colomb Cienc Exactas Fis Nat. 2003;105: 599–605.
- 30. García-Castro KL, Rangel-Medrano JD, Landínez-García RM, Márquez EJ. Population genetics of the endangered catfish Pseudoplatystoma magdaleniatum (Siluriformes: Pimelodidae) based on species-specific microsatellite loci. Neotrop Ichthyol. 2021;19: 1–15.
- 31. Pineda H, Arboleda L, Echeverry A, Urcuqui S, Pareja D, Oliveria M, et al. Caracterización de la diversidad genética en el pez Brycon henni (Characiformes: Characidae) en Colombia central por medio de marcadores RAPD. Rev Biol Trop. 2007;55: 1025–1035.
- 32. Hurtado-Alarcón JC, Mancera-Rodríguez NJ, Saldamando-Benjumea CI. Variabilidad genética de Brycon henni (Characiformes: Characidae) en la cuenca media de los ríos Nare y Guatapé, sistema Río Magdalena, Colombia. Rev Biol Trop. 2011;59: 269–282.
- 33. Landínez-García RM, Márquez EJ. Population genetics of the fish Brycon henni (Characiformes: Bryconidae) using species-specific polymorphic microsatellite loci. Rev Biol Trop. 2020;68: 847–861.
- 34. Restrepo-Escobar N, Márquez EJ. Microsatellite loci development for three catfish species from Northwestern South America. Neotrop Ichthyol. 2020;18: 1–16.
- 35. Joya CD, Landínez-García RM, Márquez EJ. Development of microsatellite loci and population genetics of the catfish Pimelodus yuma (Siluriformes: Pimelodidae). Neotrop Ichthyol. 2021;19: 1–15.
- 36. Rangel-Medrano JD, Márquez EJ. Development of microsatellite loci and population genetics in the bumblebee catfish species Pseudopimelodus atricaudus and Pseudopimelodus magnus (Siluriformes: Pseudopimelodidae). Neotrop Ichthyol. 2021;19: 1–15.
- 37. Martínez JG, Rangel-Medrano JD, Yepes-Acevedo AJ, Restrepo-Escobar N, Márquez EJ. Species limits and introgression in Pimelodus from the Magdalena-Cauca River basin. Mol Phylogenet Evol. 2022; 173:107517. pmid:35577289
- 38. Primmer CR, Painter JN, Koskinen MT, Palo JU, Merilä J. Factors affecting avian cross-species microsatellite amplification. J Avian Biol. 2005;36: 348–360.
- 39. Estoup A, Jarne P, Cornuet JM. Homoplasy and mutation model at microsatellite loci and their consequences for population genetics analysis. Mol Ecol. 2002; 11(9):1591–604. pmid:12207711
- 40. Barbará T, Palma-Silva C, Paggi GM, Bered F, Fay MF, Lexer C. Cross-species transfer of nuclear microsatellite markers: Potential and limitations. Mol Ecol. 2007;16: 3759–3767. pmid:17850543
- 41. Rutkowski R, Sielezniew M, Szostak A. Contrasting levels of polymorphism in cross-amplified microsatellites in two endangered xerothermophilous, obligatorily myrmecophilous, butterflies of the genus Phengaris (Maculinea)(Lepidoptera: Lycaenidae). Eur J Entomol. 2009; 106(4): 457–469.
- 42. Yue G-H, Kovacs B, Orban L. A new problem with cross-species amplification of microsatellites: generation of non-homologous products. Zool Res. 2010;31: 131–140. pmid:20545002
- 43. Mejía O, Betancur T, Londoño L. Aplicación de técnicas geoestadísticas en la hidrogeología del bajo Cauca antioqueño. Dyna. 2007;74 (152): 137–149.
- 44. Betancur T, Mejía O, Palacio C. Modelo hidrogeológico conceptual del Bajo Cauca antioqueño: un sistema acuífero tropical. Rev Fac Ing. 2009; 107–118.
- 45. Cruz-Esquivel Á, Marrugo-Negrete J. Concentraciones de metilmercurio en Prochilodus magdalenae (Teleostei: Curimatidae) y Hoplias malabaricus (Teleostei: Erythrinidae) en la cuenca baja del río Cauca-Magdalena, norte de Colombia. Acta Biol Colomb. 2022;27: 28–35.
- 46. Castoe TA, Poole AW, Gu W, Jason de Koning AP, Daza JM, Smith EN, et al. Rapid identification of thousands of copperhead snake (Agkistrodon contortrix) microsatellite loci from modest amounts of 454 shotgun genome sequence. Mol Ecol Resour. 2010;10: 341–347. pmid:21565030
- 47.
Rozen S, Skaletsky H. Primer3 on the World Wide Web for general users and for biologist programmers. In: Misener S, Krawetz S, editors. Bioinformatics Methods and Protocols. New Jersey: Humana Press; 2000. pp.365–386.
- 48. Rotmistrovsky K, Jang W, Schuler GD. A web server for performing electronic PCR. Nucleic Acids Res. 2004;32: 108–112. pmid:15215361
- 49. Fernández-Silva I, Whitney J, Wainwright B, Andrews KR, Ylitalo-Ward H, Bowen BW, et al. Microsatellites for Next-Generation Ecologists: A Post-Sequencing Bioinformatics Pipeline. PLoS One. 2013;8. pmid:23424642
- 50. Schoebel CN, Brodbeck S, Buehler D, Cornejo C, Gajurel J, Hartikainen H, et al. Lessons learned from microsatellite development for nonmodel organisms using 454 pyrosequencing. J Evol Biol. 2013;26: 600–611. pmid:23331991
- 51. Blacket MJ, Robin C, Good RT, Lee SF, Miller AD. Universal primers for fluorescent labelling of PCR fragments-an efficient and cost-effective approach to genotyping by fluorescence. Mol Ecol Resour. 2012;12: 456–463. pmid:22268566
- 52. Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. MICRO-CHECKER: Software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4: 535–538.
- 53. Marshall TC, Slate J, Kruuk LEB, Pemberton JM. Statistical confidence for likelihood-based paternity inference in natural populations. Mol Ecol. 1998;7: 639–655. pmid:9633105
- 54. Peakall R, Smouse PE. GENALEX 6: Genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6: 288–295.
- 55. Excoffier L, Laval G, Schneider S. Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol Bioinform Online. 2005;1: 117693430500100.
- 56. Holm S. A simple sequentially rejective multiple test procedure. Scand Stat Theory Appl. 1979;6: 65–70.
- 57. Rice WR. Analyzing Tables of Statistical Tests. Evolution. 1989;43: 223–225. pmid:28568501
- 58. Foll M, Gaggiotti O. A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: A Bayesian perspective. Genetics. 2008;180: 977–993. pmid:18780740
- 59. Nagin DS. Analyzing developmental trajectories: A semiparametric, group-based approach. Psychol Methods. 1999;4: 139–157.
- 60. Luikart G, Cornuet JM. Empirical evaluation of a test for identifying recently bottlenecked populations from allele frequency data. Conserv Biol. 1998;12: 228–237.
- 61. Piry S, Luikart G, Cornuet JM. BOTTLENECK: A computer program for detecting recent reductions in the effective population size using allele frequency data. J Hered. 1999;90: 502–503.
- 62. Garza JC, Williamson EG. Detection of reduction in population size using data from microsatellite loci. Mol Ecol. 2001;10: 305–318. pmid:11298947
- 63. Do C, Waples RS, Peel D, Macbeth GM, Tillett BJ, Ovenden JR. 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. pmid:23992227
- 64.
Crow JF, Kimura M. Introduction to population genetics theory. New York: Harper & Row; 1970.
- 65. Willoughby JR, Sundaram M, Wijayawardena BK, Kimble SJA, Ji Y, Fernandez NB, et al. The reduction of genetic diversity in threatened vertebrates and new recommendations regarding IUCN conservation rankings. Biol Conserv. 2015;191: 495–503.
- 66. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155: 945–959. pmid:10835412
- 67. Li Y-L, Liu J-X. StructureSelector: a web-based software to select and visualize the optimal number of clusters by using multiple methods. Mol Ecol Resour. 2018;18: 176–177. pmid:28921901
- 68. 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. pmid:15969739
- 69. Puechmaille SJ. The program structure does not reliably recover the correct population structure when sampling is uneven: Subsampling and new estimators alleviate the problem. Mol Ecol Resour. 2016;16: 608–627. pmid:26856252
- 70. Kopelman NM, Mayzel J, Jakobsson M, Rosenberg NA, Mayrose I. CLUMPAK: a program for identifying clustering modes and packaging population structure inferences across K. Mol Ecol Resour. 2015;15: 1179–1191. pmid:25684545
- 71. Wright S. Isolation by distance. Genetics. 1943;28: 114. pmid:17247074
- 72. Wright S. The interpretation of population structure by F-statistics with special regard to systems of mating. Evolution. 1965;19: 395–420.
- 73. Meirmans PG. Using the Amova framework to estimate a standardized genetic differentiation measure. Evolution. 2006;60: 2399–2402. pmid:17236430
- 74. Jost L. GST and its relatives do not measure differentiation. Mol Ecol. 2008;17: 4015–4026. pmid:19238703
- 75. Meirmans PG, Hedrick PW. Assessing population structure: FST and related measures. Mol Ecol Resour. 2011;11: 5–18. pmid:21429096
- 76. Jombart T. Adegenet: A R package for the multivariate analysis of genetic markers. Bioinformatics. 2008;24: 1403–1405. pmid:18397895
- 77. Rousset F. GENEPOP’007: A complete re-implementation of the GENEPOP software for Windows and Linux. Mol Ecol Resour. 2008;8: 103–106. pmid:21585727
- 78. Castoe TA, Poole AW, De Koning AJ, Jones KL, Tomback DF, Oyler-McCance , et al. Rapid microsatellite identification from Illumina paired-end genomic sequencing in two birds and a snake. PloS One. 2012;7(2): e30953. pmid:22348032
- 79. Guichoux E, Lagache L, Wagner S, Chaumeil P, Léger P, Lepais O, et al. Current trends in microsatellite genotyping. Mol Ecol Resour. 2011;11: 591–611. pmid:21565126
- 80. Bostein D, White RL, Skolnick M, Davis RW. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am J Hum Genet. 1980;32: 314–331. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1686077/ pmid:6247908
- 81.
Hilsdorf AWS, Hallerman EM. Genetic resources of neotropical fishes. Genetic Resources of Neotropical Fishes. Switzerland: Springer International Publishing; 2017.
- 82.
Franklin IR. Evolutionary changes in small populations. In: Soulé M, Wilcox B, editors. Conservation biology: an evolutionary-ecological prospective. Sunderland: Sinauer Associates; 1980. pp. 135–149.
- 83.
Soulé ME. Thresholds for Survival: Maintaining Fitness and Evolutionary Potential. In: Soulé M, Wilcox B, editors. Conservation Biology: An Evolutionary-Ecological Perspective. Sunderland: Sinauer Associates; 1980. pp. 151–169.
- 84. Frankham R, Bradshaw CJA, Brook BW. Genetics in conservation management: Revised recommendations for the 50/500 rules, Red List criteria and population viability analyses. Biol Conserv. 2014;170: 56–63.
- 85.
Márquez EJ., Restrepo-Escobar N, Yepes-Acevedo AJ, Narváez-Barandica JC. Diversidad y estructura genética de los peces de la cuenca del Río Magdalena, Colombia. In: Jiménez-Segura L, Lasso CA, editors. Peces de la cuenca del río Magdalena, Colombia: diversidad, conservación y uso sostenible. Serie Editorial Recursos Hidrobiológicos y Pesqueros Continentales de Colombia. Instituto de Investigación de Recursos Biológicos Alexander von Humboldt; 2020. pp. 115–157. https://doi.org/10.21068/a2020rrhhxix
- 86. Restrepo-Escobar N, Yepes-Acevedo AJ, Márquez EJ. Population genetics of three threatened catfish species in heterogeneous environments of the Cauca River, Colombia. Neotrop Ichthyol. 2021;19: 1–18.
- 87. Rosa IF, Oliveira DJ, Cruz VP, Ashikaga FY, Costa GO, Doretto LB, et al. Temporal genetic structure of a stock of Prochilodus lineatus (Characiformes: Prochilodontidae) in the Mogi-Guaçu River ecosystem, São Paulo, Brazil. Neotrop Ichthyol. 2022;20: 1–18.
- 88. Ribolli J, Hoeinghaus DJ, Johnson JA, Zaniboni-Filho E, de Freitas PD, Galetti PM. Isolation-by-time population structure in potamodromous Dourado Salminus brasiliensis in southern Brazil. Conserv Genet. 2017;18: 67–76.
- 89.
AUNAP-UNIMAGDALENA. Tallas mínimas de captura para el aprovechamiento sostenible de las principales especies de peces, crustáceos y moluscos comerciales de Colombia. 2013. http://sepec.unimagdalena.edu.co/Archivos/Cartilla - TALLAS MINIMAS_Digital (2).pdf
- 90.
Lasso CA, Escobar MD, Herrera J, Castellanos MC, Valencia-Rodríguez D, Campuzano J, et al. Peces introducidos en el río Magdalena y cuencas vecinas, Colombia. In: Jiménez-Segura L, Lasso CA, editors. XIX. Peces de la cuenca del río Magdalena, Colombia: diversidad, conservación y uso sostenible. Serie Editorial Recursos Hidrobiológicos y Pesqueros Continentales de Colombia. Instituto de Investigación de Recursos Biológicos Alexander von Humboldt; 2020.pp. 295–367.
- 91.
Rico G. Fauna exótica e invasora en Colombia. Instituto de Investigación de Recursos Biológicos Alexander von Humboldt. 2010. 71pp.
- 92. Caraballo P, Gandara M. El sector pesquero en el departamento de sucre, Colombia. Rev Colombiana cienc Anim RECIA. 2009;1: 273–278.
- 93. Serna-lópez JP, Cañón-barriga JE. Projecting the future of Ayapel Cienaga: A hydroecologic analysis under climate change scenarios. Rev Fac Ing Univ Antioquia. 2020; 73–87.
- 94. Reis RE, Albert JS, Di Dario F, Mincarone MM, Petry P, Rocha LA. Fish biodiversity and conservation in South America. J Fish Biol. 2016;89: 12–47. pmid:27312713