Genomic analyses reveal two species of the matamata (Testudines: Chelidae: Chelus spp.) and clarify their phylogeography

Diversity and biogeography of South American mud turtles elucidated by multilocus DNA sequencing (Testudines: Kinosternidae)

Kinosternon is the most speciose genus of extant turtles, with 22 currently recognized species, distributed across large parts of the Americas. Most species have small distributions, but K. leucostomum and K. scorpioides range from Mexico to South America. Previous studies have found discordance between mitochondrial and nuclear phylogenies in some kinosternid groups, with the current taxonomy following the nuclear-based results. Herein, based on extended molecular, geographic, and taxonomic sampling, we explore the phylogeographic structure and taxonomic limits for K. leucostomum and the K. scorpioides group and present a fossil-calibrated nuclear time tree for Kinosternon. Our results reveal contrasting differentiation patterns for the K. scorpioides group and K. leucostomum, despite overlapping distributions. Kinosternon leucostomum shows only shallow geographic divergence, whereas the K. scorpioides group is polyphyletic with up to 10 distinct taxa, some of them undescribed. We support the elevation of K. s. albogulare and K. s. cruentatum to species level. Given the deep divergence within the genus Kinosternon, we propose the recognition of three subgenera, Kinosternon, Cryptochelys and Thyrosternum, and the abandonment of the group-based classification, at least for the K. leucostomum and K. scorpioides groups. Our results show an initial split in Kinosternon that gave rise to two main radiations, one Nearctic and one mainly Neotropical. Most speciation events in Kinosternon occurred during the Quaternary and we hypothesize that they were mediated by both climatic and geological events. Additionally, our data imply that at least three South American colonizations occurred, two in the K. leucostomum group, and one in the K. scorpioides group. Additionally, we hypothesize that discordance between mitochondrial and nuclear phylogenetic signal is due to mitochondrial capture from an extinct kinosternine lineage.

Mercury biomagnification in the food chain of a piscivorous turtle species (Testudines: Chelidae: Chelus fimbriata) in the Central Amazon, Brazil

Due to their natural history and ecological attributes, turtles are excellent organisms for studies of heavy metal contamination. Turtles have a large geographical distribution, occupy different aquatic habitats, and occupy various trophic levels. The present study investigated mercury bioaccumulation in the carnivorous chelonian Chelus fimbriata (Matamata turtle) and Hg biomagnification in relation to its aquatic food chain in the middle Rio Negro, AM-Brazil. Tissue samples of muscle, carapace and claws were collected from 26 C. fimbriata individuals, as well as collections of autotrophic energy sources found in the turtle's aquatic habitat area. The samples were collected in February-March/2014 and analyzed for THg concentrations and carbon (δ13C) and nitrogen (δ15N) stable isotopes. The highest THg levels were found in claws (3780 ng.g-1), carapace (3622 ng.g-1) and muscle (403 ng.g-1), which were found to be significantly different [F(2.73) = 49.02 p < 0.01]. However, THg concentrations in muscle tissue were below the consumption threshold indicated by the WHO and Brazilian Health Ministry. The average δ13C and δ15N values in Matamata samples were -31.7‰ and 11.9‰, respectively. The principal energy source sustaining the food chain of C. fimbriata was found to be terrestrial shrubs, with smaller contributions from emergent aquatic herbaceous plants and algae, while δ15N values showed its trophic position to be two levels above the autotrophic energy sources. There was a positive correlation between THg and turtle size, while a significant relationship was found between THg and δ15N, showing strong biomagnification in the food chain of C. fimbriata: y = 0.21x + 0.46; r2 = 0.45; p < 0.001, for which the slope presented a value of 0.21.

Geographic patterns of distribution and ecological niche of the snake-necked turtle genus Hydromedusa

Biotic and abiotic factors play a crucial role in determining the distribution of species. These factors dictate the conditions that must be met for a species to thrive in a particular area. Sister species that present some degree of niche overlap can shed light on how they are distributed and coexist in their environment. This study aims to investigate the geographical distribution and ecological niche of the sister species of snake-necked turtles Hydromedusa maximiliani and H. tectifera. By analyzing their niche overlap, we aim to obtain a better understanding of how these two species coexist and which variables are determining their occurences. We applied species distribution modeling and compared the niches using the niche equivalence and similarity tests. Our findings show that the distribution of H. maximiliani is most influenced by temperature seasonality and isothermality, while H. tectifera is most affected by the temperature seasonality, precipitation of warmest quarter and mean diurnal range. In addition, our results suggest that the niche expressed by H. maximiliani retained ecological characteristics that can accurately predict the H. tectifera distribution, but the inverse is not true. In this sense, differences are not solely due to the geographic availability of environmental conditions but can reflect niche restrictions, such as competition.

Amazonian birds in more dynamic habitats have less population genetic structure and higher gene flow

Understanding the factors that govern variation in genetic structure across species is key to the study of speciation and population genetics. Genetic structure has been linked to several aspects of life history, such as foraging strategy, habitat association, migration distance, and dispersal ability, all of which might influence dispersal and gene flow. Comparative studies of population genetic data from species with differing life histories provide opportunities to tease apart the role of dispersal in shaping gene flow and population genetic structure. Here, we examine population genetic data from sets of bird species specialized on a series of Amazonian habitat types hypothesized to filter for species with dramatically different dispersal abilities: stable upland forest, dynamic floodplain forest, and highly dynamic riverine islands. Using genome-wide markers, we show that habitat type has a significant effect on population genetic structure, with species in upland forest, floodplain forest, and riverine islands exhibiting progressively lower levels of structure. Although morphological traits used as proxies for individual-level dispersal ability did not explain this pattern, population genetic measures of gene flow are elevated in species from more dynamic riverine habitats. Our results suggest that the habitat in which a species occurs drives the degree of population genetic structuring via its impact on long-term fluctuations in levels of gene flow, with species in highly dynamic habitats having particularly elevated gene flow. These differences in genetic variation across taxa specialized in distinct habitats may lead to disparate responses to environmental change or habitat-specific diversification dynamics over evolutionary time scales.

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.

Matamatas Chelus spp. (Testudines, Chelidae) have a remarkable evolutionary history of sex chromosomes with a long-term stable XY microchromosome system

The genus Chelus, commonly known as Matamata is one of the most emblematic and remarkable species among the Neotropical chelids. It is an Amazonian species with an extensive distribution throughout Negro/Orinoco and Amazonas River basins. Currently, two species are formally recognized: Chelus orinocensis and Chelus fimbriata and although it is still classified as "Least Concern" in the IUCN, the Matamatas are very appreciated and illegally sold in the international pet trade. Regardless, little is known regarding many aspects of its natural history. Chromosomal features for Chelus, for instance, are meagre and practically restricted to the description of the diploid number (2n = 50) for Chelus fimbriata, and its sex determining strategies are yet to be fully investigated. Here, we examined the karyotype of Chelus fimbriata and the newly described Chelus orinocensis, applying an extensive conventional and molecular cytogenetic approach. This allowed us to identify a genetic sex determining mechanism with a micro XY sex chromosome system in both species, a system that was likely present in their most common recent ancestor Chelus colombiana. Furthermore, the XY system found in Chelus orinocensis and Chelus fimbriata, as seen in other chelid species, recruited several repeat motifs, possibly prior to the split of South America and Australasian lineages, indicating that such system indeed dates back to the earliest lineages of Chelid species.

A Pliocene-Pleistocene continental biota from Venezuela

The Pliocene-Pleistocene transition in the Neotropics is poorly understood despite the major climatic changes that occurred at the onset of the Quaternary. The San Gregorio Formation, the younger unit of the Urumaco Sequence, preserves a fauna that documents this critical transition. We report stingrays, freshwater bony fishes, amphibians, crocodiles, lizards, snakes, aquatic and terrestrial turtles, and mammals. A total of 49 taxa are reported from the Vergel Member (late Pliocene) and nine taxa from the Cocuiza Member (Early Pleistocene), with 28 and 18 taxa reported for the first time in the Urumaco sequence and Venezuela, respectively. Our findings include the first fossil record of the freshwater fishes Megaleporinus, Schizodon, Amblydoras, Scorpiodoras, and the pipesnake Anilius scytale, all from Pliocene strata. The late Pliocene and Early Pleistocene ages proposed here for the Vergel and Cocuiza members, respectively, are supported by their stratigraphic position, palynology, nannoplankton, and ⁸⁶Sr/⁸⁸Sr dating. Mammals from the Vergel Member are associated with the first major pulse of the Great American Biotic Interchange. In contrast to the dry conditions prevailing today, the San Gregorio Formation documents mixed open grassland/forest areas surrounding permanent freshwater systems, following the isolation of the northern South American basin from western Amazonia. These findings support the hypothesis that range contraction of many taxa to their current distribution in northern South America occurred rapidly during at least the last 1.5 million years.

Plastic ingestion by freshwater turtles: a review and call to action

Plastic pollution, and especially plastic ingestion by animals, is a serious global issue. This problem is well documented in marine systems, but it is relatively understudied in freshwater systems. For turtles, it is unknown how plastic ingestion compares between marine and non-marine species. We review the relevant turtle dietary literature, and find that plastic ingestion is reported for all 7 marine turtle species, but only 5 of 352 non-marine turtle species. In the last 10 years, despite marine turtles representing just 2% of all turtle species, almost 50% of relevant turtle dietary studies involved only marine turtles. These results suggest that the potential threat of plastic ingestion is poorly studied in non-marine turtles. We also examine plastic ingestion frequency in a freshwater turtle population, finding that 7.7% of 65 turtles had ingested plastic. However, plastic-resembling organic material would have inflated our frequency results up to 40% higher were it not for verification using Raman spectroscopy. Additionally, we showcase how non-native turtles can be used as a proxy for understanding the potential for plastic ingestion by co-occurring native turtles of conservation concern. We conclude with recommendations for how scientists studying non-marine turtles can improve the implementation, quality, and discoverability of plastic ingestion research.

Cytogenetic Analysis of the Asian Box Turtles of the Genus Cuora (Testudines, Geoemydidae)

The Asian box turtle genus Cuora currently comprises 13 species with a wide distribution in Southeast Asia, including China and the islands of Indonesia and Philippines. The populations of these species are rapidly declining due to human pressure, including pollution, habitat loss, and harvesting for food consumption. Notably, the IUCN Red List identifies almost all species of the genus Cuora as Endangered (EN) or Critically Endangered (CR). In this study, we explore the karyotypes of 10 Cuora species with conventional (Giemsa staining, C-banding, karyogram reconstruction) and molecular cytogenetic methods (in situ hybridization with probes for rDNA loci and telomeric repeats). Our study reveals a diploid chromosome number of 2n = 52 chromosomes in all studied species, with karyotypes of similar chromosomal morphology. In all examined species, rDNA loci are detected at a single medium-sized chromosome pair and the telomeric repeats are restricted to the expected terminal position across all chromosomes. In contrast to a previous report, sex chromosomes are neither detected in Cuoragalbinifrons nor in any other species. Therefore, we assume that these turtles have either environmental sex determination or genotypic sex determination with poorly differentiated sex chromosomes. The conservation of genome organization could explain the numerous observed cases of interspecific hybridization both within the genus Cuora and across geoemydid turtles.

Interstitial Telomeric Repeats Are Rare in Turtles

Telomeres are nucleoprotein complexes protecting chromosome ends in most eukaryotic organisms. In addition to chromosome ends, telomeric-like motifs can be accumulated in centromeric, pericentromeric and intermediate (i.e., between centromeres and telomeres) positions as so-called interstitial telomeric repeats (ITRs). We mapped the distribution of (TTAGGG)_n repeats in the karyotypes of 30 species from nine families of turtles using fluorescence in situ hybridization. All examined species showed the expected terminal topology of telomeric motifs at the edges of chromosomes. We detected ITRs in only five species from three families. Combining our and literature data, we inferred seven independent origins of ITRs among turtles. ITRs occurred in turtles in centromeric positions, often in several chromosomal pairs, in a given species. Their distribution does not correspond directly to interchromosomal rearrangements. Our findings support that centromeres and non-recombining parts of sex chromosomes are very dynamic genomic regions, even in turtles, a group generally thought to be slowly evolving. However, in contrast to squamate reptiles (lizards and snakes), where ITRs were found in more than half of the examined species, and birds, the presence of ITRs is generally rare in turtles, which agrees with the expected low rates of chromosomal rearrangements and rather slow karyotype evolution in this group.