The complete mitochondrial genome sequences of Chelodina rugosa and Chelus fimbriata (Pleurodira: Chelidae): implications of a common absence of initiation sites (OL) in pleurodiran turtles

Complete Mitochondrial Genomes of Four Pelodiscus sinensis Strains and Comparison with Other Trionychidae Species

The Chinese soft-shelled turtle (Pelodiscus sinensis) is an important aquaculture reptile with rich nutritional and medicinal values. In recent decades, the wild resources of P. sinensis have been depleting due to natural and artificial factors. Herein, we report the complete mitochondrial genome of four P. sinensis strains, including the Japanese (RB) strain, Qingxi Huabie (HB) strain, Jiangxi (JB) strain, and Qingxi Wubie (WB) strain. The nucleotide composition within the complete mitogenomes was biased towards A + T with a variable frequency ranging from 59.28% (cox3) to 70.31% (atp8). The mitogenomes of all four strains contained 13 protein-coding genes (PCGs), 22 tRNAs, 2 rRNAs, 1 control region, and a replication origin region of the L-strand replication (OL), which was consistent with most vertebrates. Additionally, the atp8, nad4l, nad6, and nad3 genes possessed high genetic variation and can be used as potential markers for the identification of these P. sinensis strains. Additionally, all PCGs genes were evolving primarily under purifying selection. Through comparative analysis, it was revealed that most of the tRNAs were structurally different in the TψC stem, DHU stem, and acceptor stem. The length of the tandem repeats in the control region was variable in the four P. sinensis strains, ranging from 2 bp to 50 bp. Phylogenetic analysis indicated that all P. sinensis strains clustered into one branch and were closely related to other Trionychinae species. Overall, this study provides mitochondrial genome information for different P. sinensis strains to support further species identification and germplasm resource conservation.

Structural Variation of the Turtle Mitochondrial Control Region

The present study describes the most comprehensive comparison of turtle mtD-loop regions to date. The primary structure was compared from DNA sequences accessed from GenBank from 48 species in 13 families of extant turtles, and secondary structures of the mtD-loop region were inferred from thermal stabilities, using the program Mfold, for each superfamiliy of turtles. Both primary and secondary structures were found to be highly variable across the order. The Cryptodira showed conservation in the primary structure at conserved sequence blocks (CSBs), but the Pleurodira displayed limited conservation of primary structural characters, other than the coreTAS, a binding site for the helicase TWINKLE, which was highly conserved in the Central and Right Domains across the order. No secondary structure was associated with a TAS, but an AT-rich fold (secondary structure) near the 3' terminus of the mtD-loop region was detected in all turtle superfamilies. Mapping of character states of structural features of the mtD-loop region revealed that most character states were autapomorphies and inferred a number of homoplasies. The Left Domain of turtles, containing no highly conserved structural elements, likely does not serve a functional role; therefore, the Central Domain in turtles is likely equivalent to the Left Domain of mammals. The AT-rich secondary structural element near the 3' terminus of the mtD-loop region may be conserved across turtles because of a functional role, perhaps containing the Light Strand Promotor, or perhaps interacting with the TWINKLE-coreTAS complex in the Central and Right Domains to regulate mtDNA replication and transcription.

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

The matamata is one of the most charismatic turtles on earth, widely distributed in northern South America. Debates have occurred over whether or not there should be two subspecies or species recognized due to its geographic variation in morphology. Even though the matamata is universally known, its natural history, conservation status and biogeography are largely unexplored. In this study we examined the phylogeographic differentiation of the matamata based on three mitochondrial DNA fragments (2168 bp of the control region, cytochrome oxidase subunit I, and the cytochrome b gene), one nuclear genomic DNA fragment (1068 bp of the R35 intron) and 1661 Single Nucleotide Polymorphisms (SNPs). Our molecular and morphological analyses revealed the existence of two distinct, genetically deeply divergent evolutionary lineages of matamatas that separated in the late Miocene (approximately 12.7 million years ago), corresponding well to the time when the Orinoco Basin was established. As a result of our analyses, we describe the genetically and morphologically highly distinct matamata from the Orinoco and Río Negro Basins and the Essequibo drainage as a species new to science (Chelus orinocensis sp. nov.). Chelus fimbriata sensu stricto is distributed in the Amazon Basin and the Mahury drainage. Additionally, the analyses revealed that each species displays phylogeographic differentiation. For C. orinocensis, there is moderate mitochondrial differentiation between the Orinoco and the Río Negro. For C. fimbriata, there is more pronounced differentiation matching different river systems. One mitochondrial clade was identified from the Amazon, Ucayali, and Mahury Rivers, and another one from the Madeira and Jaci Paraná Rivers. The C. orinocensis in the Essequibo and Branco Rivers have haplotypes that constitute a third clade clustering with C. fimbriata. Phylogenetic analyses of the R35 intron and SNP data link the matamatas from the Essequibo and Branco with the new species, suggesting past gene flow and old mitochondrial introgression. Chelus orinocensis is collected for the pet trade in Colombia and Venezuela. However, neither the extent of the harvest nor its impact are known. Hence, it is crucial to gather more information and to assess its exploitation throughout its distribution range to obtain a better understanding of its conservation status and to design appropriate conservation and management procedures. RESUMEN: La matamata es una de las tortugas más carismáticas del mundo, ampliamente distribuida en el norte de Sudamérica. Debido a su variación morfológica geográfica, se debate sobre el reconocimiento de dos subespecies o especies. A pesar de que la matamata es universalmente conocida, su historia natural, estado de conservación y biogeografía han sido muy poco estudiados. En este estudio examinamos la diferenciación filogeográfica de las matamatas en base ​​a tres fragmentos de ADN mitocondrial (2168 pb de la región de control, la subunidad I del citocromo oxidasa y el gen del citocromo b), un fragmento de ADN genómico nuclear (1068 pb del intrón R35) y 1661 polimorfismos de nucleótido único (SNPs). Nuestros análisis moleculares y morfológicos revelaron la existencia de dos linajes evolutivos distintos de matamatas, genéticamente divergentes que se separaron en el Mioceno tardio (hace aproximadamente 12.7 millones de años), correspondiendo al tiempo en que se estableció la cuenca del Orinoco. Como resultado de nuestros análisis, describimos las genéticamente y morfológicamente distintas matamatas de las cuencas del Orinoco, Río Negro y Essequibo como una especie nueva para la ciencia (Chelus orinocensis sp. nov.). Chelus fimbriata sensu stricto se distribuye en la cuenca del Amazonas y en el drenaje del Mahury. Adicionalmente, los análisis revelaron que cada especie muestra diferenciación filogeográfica. Para C. orinocensis, hay una moderada diferenciación mitocondrial entre el Orinoco y el Río Negro. Para C. fimbriata, hay una diferenciación más pronunciada, concordando con los diferentes sistemas fluviales. Se identificó un clado de los ríos Amazonas, Ucayali y Mahury y otro de los ríos Madeira y Jaci Paraná. Las C. orinocensis de los ríos Essequibo y Branco tienen haplotipos que constituyen un tercer clado que se agrupa con C. fimbriata. Los análisis filogenéticos del intrón R35 y los datos de SNP asocian las matamatas de Essequibo y Branco con la nueva especie, sugiriendo flujo de genes pasado ​​e introgresión mitocondrial antigua. Chelus orinocensis se colecta para el comercio de mascotas en Colombia y Venezuela. Sin embargo, ni se conoce el alcance de las colectas ni su impacto. Por lo tanto, es crucial recopilar más información y evaluar su explotación en todo su rango de distribución, comprender mejor su estado de conservación y para diseñar acciones apropiadas de conservación y manejo.

The complete mitochondrial genome of the Cyclemys fusca (Chelonia: Geoemydidae)

The complete mitochondrial genome (mitogenome) of Cyclemys fusca was obtained and characterized in this study. The circular molecule is 16,491 bp in length and contained 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and one non-coding region (control region). Its gene arrangement type is identical to the type of most vertebrate. All protein-coding genes initiate with ATG as start codon, except for COI started with GTG. Interestingly, COI and ND6 end up with AGG. The complete mitogenome of C. fusca provides the basic data to research molecular systematics of Geoemydidea.

Molecular phylogeny and divergence of the map turtles (Emydidae: Graptemys)

The map turtles (genus Graptemys) comprise a morphologically diverse clade that forms a major component of the southeastern US hotspot of chelonian diversity. Map turtles have experienced both recent and rapid diversification resulting in long-standing uncertainty regarding species boundaries and phylogenetic relationships within the genus as well as timing of their divergence. We present a phylogeny for the group that includes geographically representative sampling for all described species and subspecies. We make use of an empirical prior on rates of molecular evolution to estimate divergence times with a molecular clock under a coalescent framework. Together, the phylogeny and divergence time estimates suggest that diversification has been both more recent and more rapid than has so far been suspected. We provide a well-supported evolutionary framework for Graptemys that is necessary for understanding map turtle diversity, biogeography, and for conservation of this threatened clade of turtles.

Resolution of the enigmatic phylogenetic relationship of the critically endangered Western Swamp Tortoise Pseudemydura umbrina (Pleurodira: Chelidae) using a complete mitochondrial genome

Pseudemydura umbrina is one of the most endangered turtle species in the world, and the imperative for its conservation is its distinctive morphology and relict status among the Chelidae. We use Illumina sequencing to obtain the complete mitogenome for resolving its uncertain phylogenetic position. A novel nuclear paralogue confounded the assembly, and resolution of the authentic mitogenome required further Sanger sequencing. The P. umbrina mitogenome is 16,414bp comprising 37 genes organized in a conserved pattern for other vertebrates. The nuclear paralogue is 547bp, 97.8% identity to the corresponding mitochondrial sequence. Particular features of the mitogenome include an nd3 174+1A frameshift, loss of DHC loop in tRNA^Ser (AGN), and a light-strand replication initiation site in Wancy region that extends into an adjacent tRNA gene. Phylogenetic analysis showed that P. umbrina is the monotypic sister lineage to the remaining Australasian Chelidae, a lineage probably dating back to the Cretaceous.

Complete mitochondrial genome of the endangered Mary River turtle (Elusor macrurus) and low mtDNA variation across the species' range

Elusor macrurus is an endangered short-necked turtle restricted to the Mary River catchment in south-eastern Queensland. Shotgun sequencing of genomic DNA was used to generate a complete mitochondrial genome sequence for E. macrurus using the Illumina MiSeq platform. The mitogenome is 16 499 base pairs (bp) long with 37 genes arranged in the typical vertebrate order and a relatively short 918-bp control region, which does not feature extensive tandem repeats as observed in some turtles. Primers were designed to amplify a 1270-bp region that includes 81% of the typically hypervariable control region. Two haplotypes were detected in a sample of 22 wild-caught individuals from eight sites across its natural range. The Mary River turtle is a species with low mtDNA nucleotide variability relative to other Chelidae. The combination of a very restricted distribution and dramatic reduction in population size due to exploitation for the pet trade are the conditions likely to have led to very low mtDNA variability in this endangered species.