Population genetics of the endangered South American freshwater turtle, Podocnemis unifilis, inferred from microsatellite DNA data

The Impact of Weir Construction in Korea's Nakdong River on the Population Genetic Variability of the Endangered Fish Species, Rapid Small Gudgeon (Microphysogobio rapidus)

Microphysogobio rapidus, an endemic cyprinid fish species found exclusively in Korea, has been identified in only two tributaries of the Nakdong River. The species predominantly occupies the near-gravel bottom waters within shallow sections of the middle and lower reaches of the river, characterized by swift currents. M. rapidus is currently recognized as a critically endangered species due to its distinct habitat preference, as well as the negative impacts of stream dam development and water environment pollution. In this study, we used 10 microsatellite markers to examine the genetic diversity of M. rapidus in the upper Nam (UN), lower Nam (LN), and Deokcheon Rivers (DC) in Korea, with a specific focus on assessment of the impact of dam development. Fish sampled from the UN and LN showed a greater average number of alleles and allelic richness (A = 18.3-18.4, AR = 13.8) compared to those from DC (A = 11.8, AR = 11.5). The observed heterozygosity among the fish examined ranged from HO = 0.748 (LN) to 0.766 (DC). All three fish groups exhibited a significant departure from Hardy-Weinberg equilibrium (HWE) (p < 0.05). Despite having the largest effective population size (Ne = 175 and 157, respectively), the fish sampled from UN and LN showed the highest inbreeding coefficients (FIS = 0.056-0.053, respectively), which were highly significant (p < 0.01). In contrast, the fish sampled from DC exhibited the smallest effective population size (Ne = 61) and showed an inbreeding coefficient close to zero (p > 0.05). BOTTLENECK analysis and estimated M-ratio values (0.341-0.372) revealed indications of past population size reduction in all fish groups examined. No significant genetic differentiation (FST < 0.05) was detected using the DAPC, STRUCTURE, and AMOVA among the fish studied. However, pairwise comparisons of FST between fish sampled from the Nam and Deokcheon Rivers revealed significant values (p < 0.001) ranging from 0.013 to 0.014. In addition, the closest genetic distance (0.026) was observed between UN and LN, while the greatest distance (0.087) was found between UN and DC. Analysis of gene flow rates among the fish examined indicated asymmetrical gene exchange within the Nam River, which was 31.51% in the downstream direction (from UN to LN), with a minimal gene flow rate (0.41%) in the upstream (from LN to UN) direction. The opposite trend was recorded between DC and LN, with a higher gene flow rate (29.74%) in the upstream direction compared to the downstream direction (0.12%). Our study highlighted the importance of implementing long-term conservation efforts focused on maintaining river integrity by removing water barriers such as weirs that impede fish migration and implementing active protection measures, such as aquaculture breeding and reasonable stocking practices, to preserve M. rapidus in the study area.

Phylogeography and population genetic structure of the cardinal tetra (Paracheirodon axelrodi) in the Orinoco basin and Negro River (Amazon basin): evaluating connectivity and historical patterns of diversification

The Neotropics contain one of the most diverse assemblages of freshwater fishes worldwide. Part of this diversity is shared between the Orinoco and Amazon basins. These basins have been separated for a long time due to the Vaupes Arch, rising between 10-11 Ma. Today, there is only one permanent connection between the Orinoco and Negro (Amazon) basins, known as the Casiquiare Canal. However, alternative corridors allowing fish dispersion between both basins have been proposed. The cardinal tetra (Paracheirodon axelrodi), the most important fish in the ornamental world market, is distributed in both basins. Here we investigated P. axelrodi phylogeography, population structure, and potential routes of migration and connectivity between the two basins. A total of 468 bp of the mitochondrial gene (COI), 555 bp of the nuclear gene fragment (MYH6), and eight microsatellite loci were analyzed. As a result, we found two major genetic clusters as the most likely scenario (K = 2), but they were not discreetly distributed between basins. A gradient of genetic admixture was observed in Cucui and São Gabriel da Cachoeira, between the upper Negro River and the upper Orinoco. Samples from the middle-lower Negro River were highly structured. Cucui (Negro basin) was more similar to the Orinoco than to the rest of the Negro basin populations. However, substructure was also observed by the discriminant analysis, fixation indices and other hierarchichal structure analyses (K = 3 - 6), showing three major geographic clusters: Orinoco, Cucui, and the remaining Negro basin. Unidirectional migration patterns were detected between basins: via Cucui toward Orinoco and via the remaining of the Negro basin toward Orinoco. Results from the Relaxed Random Walk analysis support a very recent origin of this species in the headwater Orinoco basin (Western Guiana Shield, at late Pleistocene) with a later rapid colonization of the remaining Orinoco basin and almost simultaneously the Negro River via Cucui, between 0.115 until about 0.001 Ma. Historical biogeography and population genetic patterns observed here for Cardinal tetra, seem to be better explained by river capture, physical, or ecological barriers than due to the geographic distance.

Hemogregarine Diversity Infecting Brazilian Turtles with a Description of Six New Species of Haemogregarina (Apicomplexa: Adeleorina: Haemogregarinidae)

(1) Background: Hemoparasites of the genus Haemogregarina (Haemogregarinidae) are commonly reported in freshwater turtles. However, in Brazil, only three species have currently been characterised using molecular methods. This study aimed to bring new insights on the diversity of species of Haemogregarina infecting Brazilian freshwater turtles from Mato Grosso and Goiás states using molecular and morphological tools. (2) Methods: In total, 2 mL of blood was collected, with between two to five blood smears prepared, fixed with absolute methanol, and stained with a 10% Giemsa solution. Blood was stored at −20 °C for molecular analysis targeting the 18S rRNA gene. Fragments of the organs (liver, spleen, heart, and kidney) were separated and stained with hematoxylin-eosin. (3) Results: A total of 40 turtles were screened and hemogregarines were observed in 34 (85%) individuals (3 Podocnemis expansa and 31 Podocnemis unifilis). Haemogregarina embaubali and six new species of Haemogregarina were observed and formally described as follows: H. unifila n. sp., H. rubra n. sp., H. goianensis n. sp., H. araguaiensis n. sp., H. tigrina n. sp., and H. brasiliana n. sp. (4) Conclusions: This study contributes to the diversity and knowledge of Brazilian fresh-water turtle blood parasites, using integrative approaches for diagnosing and characterizing hemoparasites, with the identification of six undescribed species.

Genomic analyses of the scorpion mud turtle (Kinosternon scorpioides) (Linnaeus, 1766) in insular and continental Colombia: Evidence for multiple conservation and taxonomic units

The turtle genus Kinosternon is widespread with at least 25 species distributed from Mexico to northern Argentina. The taxonomy of this genus is controversial and requires a full revision using both morphological and molecular approaches. In this study, we did a genomic analysis on the species Kinosternon scorpioides distributed in insular and continental Colombia in order to define conservation units. Total DNA was extracted from 24 tissue samples and RADseq genotyping analysis was done. In addition, the intron R35 was amplified and sequenced for a subset of samples. A total of 35,507 SNPs combined with 1,047 bp of the intron were used for spatiotemporal colonization pattern reconstruction and phylogenetic analyses. In addition, SNPs were used for population structure inferences and allele frequency-based analyses. Reciprocal monophyly, significant differences in allele frequencies (Fst = 0.32 - 0.78), and evidence of reproductive isolation (no admixture/geneflow), indicate long-term divergence between groups (2-8 MYA), possibly due to geographical barriers. Four Evolutionarily Significant Units (ESUs) were defined within our samples. One ESU was represented by the insular subspecies K. scorpioides albogulare, found on San Andrés island, and three ESUs were defined for the subspecies K. s. scorpioides in continental Colombia: one trans-Andean, found in northwestern Colombia (Caribbean region) and two cis-Andean, found in eastern and southeastern Colombia in the Orinoco and Amazon regions, respectively. Colonization of this species occurred from an ancestral area on South of Central America region (~ 8.43 MYA), followed by an establishing of current populations on San Andrés Island and then, in the continent. First, in the Colombian Caribbean, next, in the Orinoco, and more recently, in the Amazon. We hypothesize that the emergence of the Panamá Isthmus, as well as the final uplift of the North Eastern Andes and Vaupes Arch, were key event leading to the differentiation of these ESUs. For management and conservation purposes, each of these ESUs should be considered as a separate management unit. A full revision of the taxonomy of the genus Kinosternon is warranted.

Are genomic updates of well-studied species worth the investment for conservation? A case study of the Critically Endangered Magdalena River turtle

As genomic-scale data sets become economically feasible for most organisms, a key question for conservation biology is whether the increased resolution offered by new genomic approaches justifies repeating earlier studies based on traditional markers, rather than investing those same time and monetary resources in less-known species. Genomic studies offer clear advantages when the objective is to identify adaptive loci that may be critical to conservation policy-makers. However, the answer is far less certain for the population and landscape studies based on neutral loci that dominate the conservation genetics research agenda. We used RADseq to revisit earlier molecular studies of the IUCN Critically Endangered Magdalena River turtle (Podocnemis lewyana), documenting the conservation insights gained by increasing the number of neutral markers by several orders of magnitude. Earlier research indicated that P. lewyana has the lowest genetic diversity known for any chelonian, and little or no population differentiation among independent rivers. In contrast, the RADseq data revealed discrete population structure with isolation-by-distance within river segments and identified precise population breaks clearly delineating management units. It also confirmed that the species does not have extremely low heterozygosity and that effective population sizes are probably sufficient to maintain long-term evolutionary potential. Contrary to earlier inferences from more limited population genetic markers, our genomic data suggest that management strategies should shift from active genetic rescue to more passive protection without extreme interventions. We conclude with a list of examples of conservation studies in other vertebrates indicating that for many systems a genomic update is worth the investment.

Conservation genetics of the endangered Spotted Turtle (Clemmys guttata) illustrate the risks of “bottleneck tests”

Studies of population genetics in turtles have suggested that turtles do not experience genetic impacts of bottlenecks as strongly as expected. However, recent studies cast doubt on two commonly used tests implemented in the program BOTTLENECK, suggesting that these findings should be re-evaluated. The Spotted Turtle (Clemmys guttata (Schneider, 1792)) is endangered both globally and within Canada, but genetic data required to develop effective recovery strategies are unavailable. Here, we conducted the first study of population genetic structure in C. guttata. We then used multiple small populations of C. guttata as replicates to test whether the commonly used program BOTTLENECK could detect the genetic signature of bottlenecks in our study populations, which are all thought to have experienced significant declines in the past 2–3 generations (75 years). Turtles (n = 256) were genotyped at 11 microsatellite loci. A suite of Bayesian population genetics analyses and a principal coordinates analysis identified a minimum of 6 distinct genetic populations and a maximum of 10 differentiated subpopulations across the sampled Canadian range of C. guttata, which corresponded to demographically independent units. BOTTLENECK failed to detect population declines. A literature review found that bottleneck tests in 17 of 18 previous genetic studies of tortoises and freshwater turtles were based on suboptimal sampling, potentially confounding their results. High retention of genetic diversity (allelic richness and heterozygosity) in isolated populations of C. guttata and other turtle species is encouraging for species recovery, but conclusions about the prevalence of genetic bottlenecks in such populations should be re-examined.

Cryptic population structuring and the role of the Isthmus of Tehuantepec as a gene flow barrier in the critically endangered Central American River Turtle

The critically endangered Central American River Turtle (Dermatemys mawii) is the only remaining member of the Dermatemydidae family, yet little is known about its population structuring. In a previous study of mitochondrial (mt) DNA in the species, three main lineages were described. One lineage (Central) was dominant across most of the range, while two other lineages were restricted to Papaloapan (PAP; isolated by the Isthmus of Tehuantepec and the Sierra de Santa Marta) or the south-eastern part of the range (1D). Here we provide data from seven polymorphic microsatellite loci and the R35 intron to re-evaluate these findings using DNA from the nuclear genome. Based on a slightly expanded data set of a total of 253 samples from the same localities, we find that mtDNA and nuclear DNA markers yield a highly congruent picture of the evolutionary history and population structuring of D. mawii. While resolution provided by the R35 intron (sequenced for a subset of the samples) was very limited, the microsatellite data revealed pronounced population structuring. Within the Grijalva-Usumacinta drainage basin, however, many populations separated by more than 300 kilometers showed signals of high gene flow. Across the entire range, neither mitochondrial nor nuclear DNA show a significant isolation-by-distance pattern, but both genomes highlight that the D. mawii population in the Papaloapan basin is genetically distinctive. Further, both marker systems detect unique genomic signals in four individuals with mtDNA clade 1D sampled on the southeast edge of the Grijalva-Usumacinta basin. These individuals may represent a separate cryptic taxon that is likely impacted by recent admixture.

Phylogeographical patterns shed light on evolutionary process in South America

The South American continent is composed of several biogeographical regions harbouring the highest biodiversity on the globe, encompassing five of the world's biodiversity 'hot spots'. Nonetheless, the patterns and processes responsible for shaping its astonishing species diversity are largely unknown. Here, we present a review of current South American phylogeographical knowledge based on published articles on this topic. An appraisal of the literature reveals emerging phylogeographical patterns in the biota of South America. The striking phylogeographical divergence observed among organism lineages in South American studies is suggestive of high levels of undocumented species diversity. The interplay between Pleistocene climatic oscillations and Pliocene/Miocene orogenic events has contributed to shaping the current diversity and distribution of modern lineages in both the tropical and temperate regions of South America. Although older divergence times were observed for a range of species, most herpetofauna underwent an intraspecific lineage split much earlier than other organisms. The geographical ranges of species associated with forest habitats were reduced mainly during glacial cycles, whereas species associated with open vegetation domains have shown variable responses to climatic oscillations. The results suggest a highly complex mosaic of phylogeographical patterns in South America. We suggest future research directions to promote a better understanding of the origin and maintenance of the South American biota.

A trophic study of the sympatric Amazonian freshwater turtles Podocnemis unifilis and Podocnemis expansa (Testudines, Podocnemidae) using carbon and nitrogen stable isotope analyses

The Yellow-spotted River Turtle ( Podocnemis unifilis Troschel, 1848) and the South American River Turtle (Podocnemis expansa (Schweigger, 1812)) are two turtles species that are widely distributed and have ecological, economic, and cultural importance in the Amazon basin. Although sympatric regarding most of their distribution, few studies have addressed the coexistence of these two species. To examine this, we analyzed the trophic level and the primary carbon source from the diets of both species in Baixo Araguaia, Tocantins, Brazil, using stable isotope analyses of carbon (δ13C) and nitrogen (δ15N). We also verified possible intraspecific variations (related to sex and body mass) in the trophic levels and primary carbon sources of their diets. Podocnemis unifilis had higher values of δ15N than P. expansa, averaging 7.59‰ and 5.06‰, respectively, a difference which may indicate a possible trophic change owing to exploiting different food resources. No differences were found between the two species in relation to δ13C (mean values of –26.2‰ and –26.1‰, respectively). The similarity between δ13C values suggests that the sources of their basal feeding are the same, consisting mainly of C3 plants. There was no intraspecific variation in the values of δ13C and δ15N.

Nine microsatellite markers for the Australian side-necked turtle Chelodina expansa (Chelidae) and cross species amplifications

Nine microsatellite DNA loci for the Australian broad-shelled freshwater turtle (Chelodina expansa) are presented. Markers were tailed with 20-mer oligonucleotides for use in four-colour fluorescent multiplex PCRs. The markers show high allelic richness (mean N(A)=10.9, range 2-38) and expected heretozygosity (mean H(E)=0.643; range 0.161-0.963) indicating that they will be valuable for population genetics studies in C. expansa. Cross-species amplification in three Australian freshwater turtle species further highlights the potential utility of these markers, particularly in the side-neck species C. longicollis and C. rugosa.