The complete mitochondrial genome of the Korean soft-shelled turtle Pelodiscus sinensis (Testudines, Trionychidae)

Complete mitochondrial genome of the southern painted turtle (Chrysemys dorsalis, Testudines: Emydidae) in Korea

The complete mitochondrial genome of Chrysemys dorsalis in Korea was sequenced and characterized. The mitochondrial genome is 17,258 bp in length and the GC content is 39%. It is constituted of 37 genes, 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and a noncoding region. Phylogenetic analysis reveals that C. dorsalis forms a monophyletic group with C. picta turtles but is distinctly separated from them, aligning with previous findings. In Korea, C. dorsalis forms a discrete clade, separate from both native and invasive turtle species. No evidence of genetic disturbance or intermingling is observed. This is the first case of a complete mitochondrial genome from C. dorsalis and provides crucial data for understanding C. dorsalis and managing invasive species effectively, emphasizing the need for continued mitochondrial genome data accumulation.

Complete mitochondrial genome of the razor-backed musk turtle (Sternotherus carinatus, testudines: emydidae) in Korea

Sternotherus carinatus has been considered as a potential invasive species in Korea. However, the mitochondrial genome information of S. carinatus which can be used to control its effect on ecosystem is lacking. In this study, the complete mitochondrial genome of S. carinatus in Korea was sequenced and characterized. The mitochondrial genome consists of 37 genes (13 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes) and a noncoding region. Phylogenetic analysis based on the mitochondrial genome sequences showed that S. carinatus from Korea is separated from other turtles which are the invasive species in Korea. Sequence divergence calculations indicated near-zero divergence between S. carinatus populations in Korea, the USA, and China, suggesting limited genetic differentiation. In the context of the broader issue of invasive species disrupting ecosystems, this research contributes to the identification of mitochondrial genomes for various freshwater turtle species, emphasizing the need for extended data collection to discern genetic mixing trends between native and non-native species. This study is a significant step toward managing S. carinatus as a potential invasive species in Korea.

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.

The complete mitochondrial genome of the Amur soft-shelled turtle (Pelodiscus maackii Brandt, 1858), from South Korea

In this study, we use a specimen from wild-caught individual to determine the complete mitochondrial genome of the Amur soft-shelled turtle (Pelodiscus maackii). The complete mitogenome of P. maackii has 16,258 bp in length and consists of 13 protein-coding genes (PCGs), 22 tRNAs, two rRNAs, and one control region. The arrangement of genes of P. maackii is identical with previously reported mitogenomes in the family Trionychoidea. According to our result, the ML tree for the phylogenetic reconstruction revealed that the individuals used in present study is closely related with the previously reported sequences of P. sinensis (AY962573 and MG431983) in p-distance 0.7% and 2.5%.

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.

Near-complete phylogeny of extant Crocodylia (Reptilia) using mitogenome-based data

Species of the order Crocodylia are mostly large, predatory and semi-aquatic reptiles. Crocodylia, the closest living relatives of birds, first appeared in the Late Cretaceous period. In the present study, the complete mitochondrial (mt) genomes of 19 Crocodylia species, including two species (Melanosuchus niger and Caiman yacare) that have not been previously sequenced for mitogenomes, were processed through Illumina sequencing to offer genetic resources and compare with the mitogenomes of Crocodylia species reported previously. In addition, a high-resolution phylogenetic tree of nearly all current recognized species of Crocodylia is constructed based on mitogenomic data. Phylogenetic analyses support monophyly of three families: Alligatoridae (four genera: Alligator, Caiman, Melanosuchus and Paleosuchus), Crocodylidae (three genera: Crocodylus, Mecistops and Osteolaemus) and Gavialidae (two genera: Gavialis and Tomistoma). The tree topology is generally similar to previous studies. Molecular dating suggests that the first split within Crocodylia date back to the Upper Cretaceous (approx. 86.75 Mya). The estimated time to the most recent common ancestor (TMRCA) of Alligatoridae is 53.33 Mya and that of Crocodylidae and Gavialidae is 50.13 Mya, which might be closely linked to climate changes during the Late Palaeocene and Early Eocene. Additionally, this study proves that the diversification rate within Crocodylia began to increase from the Late Eocene (about 36 Mya) and two diversification peak periods of Crocodylia (0–10 Mya and 10–20 Mya) are disclosed, which is roughly consistent with the estimated crocodylian species richness through time. Combining all these clues, we can suggest that climate fluctuation may have played a decisive role in the speciation of Crocodylia.

Structural analysis of leucine, lysine and tryptophan mitochondrial tRNA of nesting turtles Caretta caretta (Testudines: Chelonioidea) in the Colombian Caribbean

The understanding of the functional properties of mitochondrial transfer RNA (mt tRNAs) depend on the knowledge of its structure. tRNA acts as an interface between polynucleotides and polypeptides thus, they are key molecules in protein biosynthesis. The tRNA molecule has a functional design and, given its importance in the translation of mitochondrial genes, it is plausible that modifications of the structure can affect the synthesis of proteins and the functional properties of the mitochondria. In a previous work, the mitochondrial genome of an individual of the nesting Caretta caretta of the Colombian Caribbean was obtained, where specific mutations were identified in the only tRNALeu (CUN), tRNATrp and tRNALys genes. In order to analyze the effect of these mutations on these three mt tRNAs, the prediction of 2D and 3D structures was performed. Genes were sequenced in 11 nesting loggerhead turtles from the Colombian Caribbean. Two-dimensional structures were inferred using the ARWEN program, and three-dimensional structures were obtained with the RNA Composer 3D program. Two polymorphisms were identified in tRNATrp and another one was located in tRNALys, both specific to C. caretta. The thymine substitution in nucleotide position 14 of tRNATrp could constitute an endemic polymorphism of the nesting colony of the Colombian Caribbean. Two 2D and three 3D patterns were obtained for tRNATrp. In the case of tRNALys and tRNALeu 2D and 3D structures were obtained respectively, which showed compliance to canonical structures, with 4 bp in the D-arm, 4-5 bp in the T-arm, and 5 bp in the anticodon arm. Moderate deviations were found, such as a change in the number of nucleotides, elongation in loops or stems and non-Watson-Crick base pairing: adenine-adenine in stem D of tRNATrp, uracil-uracil and adenine-cytosine in the acceptor arm of the tRNALys and cytosine-cytosine in the anticodon stem of the tRNALeu. In addition, distortions or lack of typical interactions in 3D structures gave them unique characteristics. According to the size of the variable region (4-5 nt), the three analyzed tRNAs belong to class I. The interactions in the three studied tRNAs occur mainly between D loop-variable region, and between spacer bases-variable region, which classifies them as tRNA of typology II. The polymorphisms and structural changes described can, apparently, be post-transcriptionally stabilized. It will be crucial to perform studies at the population and functional levels to elucidate the synthetic pathways affected by these genes. This article analyses for the first time the 1D, 2D and 3D structures of the mitochondrial tRNALys, tRNATrp and tRNALeu in the loggerhead turtle.

Comparative mitogenomic and phylogenetic analysis of Apalone spinifera and Apalone ferox (Testudines: Trionychidae)

The soft-shell turtles Apalone spinifera (AS) and Apalone ferox (AF) are two important economic species. AF is found in the Yellow River of China, and is a confirmed member of the Trionychidae family. However, the classification of AS was in dispute. Mitochondrial genomes (mitogenomes) have been widely used for species identification, as well as population and phylogenetic analysis. In order to understand the phylogenetic and mitogenomic features of AS and AF, the complete mitogenomes were sequenced, annotated and analyzed in this study. The complete mitogenomes of AS and AF are 16,817 bp and 16,756 bp in length, respectively. Both mitogenomes contain 37 genes, seven short intergenic spacers and two long intergenic spacers. Comparative analysis showed that there are 1,137 variation sites (6.79%) between the two mitogenomes. AS and AF mitogenomes both show a usage preference in terms of nucleotides, codons and amino acids. In addition, the non-synonymous substitution rate/synonymous substitution rate indicates that all protein-coding genes (PCGs) have undergone a strong purifying selection. Phylogenetic trees constructed by 13 PCGs show a clear phylogenetic relationship of the soft-shell turtles and suggest that AS is a sister species to AF of the genus Apalone. The data could be useful for further research of species identification, population analysis and the mitogenomic features of soft-shell turtles.

Complete mitochondrial genome of Black Soft-shell Turtle (Nilssonia nigricans) and comparative analysis with other Trionychidae

The characterization of mitochondrial genome has been evidenced as an efficient field of study for phylogenetic and evolutionary analysis in vertebrates including turtles. The aim of this study was to distinguish the structure and variability of the Trionychidae species mitogenomes through comparative analysis. The complete mitogenome (16796 bp) of an endangered freshwater turtle, Nilssonia nigricans was sequenced and annotated. The mitogenome encoded for 37 genes and a major non-coding control region (CR). The mitogenome was A + T biased (62.16%) and included six overlapping and 19 intergenic spacer regions. The Relative synonymous codon usage (RSCU) value was consistent among all the Trionychidae species; with the exception of significant reduction of Serine (TCG) frequency in N. nigricans, N. formosa, and R. swinhoei. In N. nigricans, most of the transfer RNAs (tRNAs) were folded into classic clover-leaf secondary structures with Watson-Crick base pairing except for trnS1 (GCT). The comparative analysis revealed that most of the tRNAs were structurally different, except for trnE (TTC), trnQ (TTG), and trnM (CAT). The structural features of tRNAs resulted ≥ 10 mismatched or wobble base pairings in 12 tRNAs, which reflects the nucleotide composition in both H- and L-strands. The mitogenome of N. nigricans also revealed two unique tandem repeats (ATTAT)8, and (TATTA)20 in the CR. Further, the conserved motif 5'-GACATA-3' and stable stem-loop structure was detected in the CRs of all Trionychidae species, which play an significant role in regulating transcription and replication in the mitochondrial genome. Further, the comparative analysis of Ka/Ks indicated negative selection in most of the protein coding genes (PCGs). The constructed Maximum Likelihood (ML) phylogeny using all PCGs showed clustering of N. nigricans with N. formosa. The resulting phylogeny illustrated the similar topology as described previously and consistent with the taxonomic classification. However, more sampling from different taxonomic groups of Testudines and studies on their mitogenomics are desirable for better understanding of the phylogenetic and evolutionary relationships.

Millennium-old farm breeding of Chinese softshell turtles (Pelodiscus spp.) results in massive erosion of biodiversity

Chinese softshell turtles (Pelodiscus spp.) are widely distributed, ranging from the Amur and Ussuri Rivers in the Russian Far East through the Korean Peninsula, Japan, and eastern, central, and southern China to southern Vietnam. In East and Southeast Asia, Chinese softshell turtles are traditionally exploited for food and have been farm-bred in China since the Spring and Autumn Period, more than 2400 years ago. Currently, the annual production of Pelodiscus amounts to 340,000 t in China alone. Using mitochondrial DNA (2428 bp) and five nuclear loci (3704 bp), we examined broad sampling of wild and farm-bred Pelodiscus to infer genetic and taxonomic differentiation. We discovered four previously unknown mitochondrial lineages, all from China. One lineage from Jiangxi is deeply divergent and sister to the mitochondrial lineage of Pelodiscus axenaria. The nuclear loci supported species status for P. axenaria and the new lineage from Jiangxi. Pelodiscus maackii and P. parviformis, both harboring distinct mitochondrial lineages, were not differentiated from P. sinensis in the studied nuclear markers. The same is true for two new mitochondrial lineages from Zhejiang, China, represented by only one individual each, and another new lineage from Anhui, Guangdong, Jiangxi and Zhejiang, China. However, Vietnamese turtles yielding a mitochondrial lineage clustering within P. sinensis were distinct in nuclear markers, suggesting that these populations could represent another unknown species with introgressed mitochondria. Its species status is also supported by the syntopic occurrence with P. sinensis in northern Vietnam and by morphology. In addition, we confirmed sympatry of P. axenaria and P. parviformis in Guangxi, China, and found evidence for sympatry of P. sinensis and the new putative species from Jiangxi, China. We also discovered evidence for hybridization in turtle farms and for the occurrence of alien lineages in the wild (Zhejiang, China), highlighting the risk of genetic pollution of native stock. In the face of the large-scale breeding of Pelodiscus, we claim that the long-term survival of distinct genetic lineages and species can only be assured when an upscale market segment for pure-bred softshell turtles is established, making the breeding of pure lineages lucrative for turtle farms. Our findings underline that the diversity of Pelodiscus is currently underestimated and threatened by anthropogenic admixture. We recommend mass screening of genetic and morphological variation of Chinese softshell turtles as a first step to understand and preserve their diversity.

Intraspecific variations in Cyt b and D-loop sequences of Testudine species, Lissemys punctata from south Karnataka

The freshwater Testudine species have gained importance in recent years, as most of their population is threatened due to exploitation for delicacy and pet trade. In this regard, Lissemys punctata, a freshwater terrapin, predominantly distributed in Asian countries has gained its significance for the study. A pilot study report on mitochondrial markers (Cyt b and D-loop) conducted on L. punctata species from southern Karnataka, India was presented in this investigation. A complete region spanning 1.14 kb and ∼1 kb was amplified by HotStart PCR and sequenced by Sanger sequencing. The Cyt b sequence revealed 85 substitution sites, no indels and 17 parsimony informative sites, whereas D-loop showed 189 variable sites, 51 parsimony informative sites with 5′ functional domains TAS, CSB-F, CSBs (1, 2, 3) preceding tandem repeat at 3′ end. Current data highlights the intraspecific variations in these target regions and variations validated using suitable evolutionary models points out that the overall point mutations observed in the region are transitions leading to no structural and functional alterations. The mitochondrial data generated uncover the genetic diversity within species and conservationist can utilize the data to estimate the effective population size or for forensic identification of animal or its seizures during unlawful trade activities.

Establishment and characterization of a cell line from the Chinese soft-shelled turtle Pelodiscus sinensis

The establishment and partial characterization of Pelodiscus sinensis continuous cell line is described here. A novel P. sinensis fibroblast cell line, designated PSF, was established from heart tissue by the semi-digestion explant culture technique. Since its initiation in July 2013, the cell line has been subcultured at 30°C in minimal essential medium (MEM) containing 15% (v/v) fetal bovine serum for more than 50 passages. The growth curve of the cell line revealed the population doubling time was 51.1 h. Karyotyping analysis indicated the modal chromosome number was 66, and no microbial contamination was detected. The PSF cell line produced significant fluorescent signals after transfection with plasmid pEGFP-C3. Analysis of mitochondrial cytochrome D-loop sequences revealed 96% identity among other Chinese turtle subspecies. Several cell line characterizations included morphological analysis and immunocytochemistry, which revealed the origin of the PSF cell line was fibroblast-like cells. Measurement of the isoenzymes lactic dehydrogenase and malic dehydrogenase showed no cross-contamination of this cell line with other species. This newly established cell line will be a valuable tool for transgenic and genetic manipulation studies and will act as an efficient instrument for studies of the viral diseases of the soft-shelled turtle.

The complete mitochondrial genome of the loggerhead turtle Caretta caretta (Testudines: Cheloniidae): genome description and phylogenetic considerations

The marine turtle Caretta caretta is a widely distributed species that is facing critical population decline, especially in the Mediterranean rookeries. Molecular markers, such as mitochondrial DNA (mtDNA) sequences, are of great importance for the description and monitoring of turtle migratory populations. The complete sequence of the C. caretta mitochondrial genome is presented here. The genome comprises 16,440 base pairs, containing 37 genes (13 protein-coding genes, 22 tRNA genes, and 2 rRNA genes), and a control region, all organized similar to the majority of vertebrate mitogenomes. MtDNA length polymorphism and heteroplasmy were observed among, and within, individuals due to the variable size of a microsatellite repeat residing at the 3' end of the control region. The use of the above repeat as a marker for individual fingerprinting is discussed. Furthermore, phylogenetic analyses among Testudines based on complete mitogenomes, as well as among marine turtles based on partial mtDNA sequences, are considered.

Genetic variation of the Nile soft-shelled turtle (Trionyx triunguis)

We studied the genetic structure of Trionyx triunguis populations from the Mediterranean and African continent based on mtDNA D-loop (776 bp) and nine microsatellite loci. A total of 102 polymorphic sites and 13 mtDNA haplotypes were described. Nucleotide diversity and haplotypes diversity were 0.047 and 0.974 respectively. Both mtDNA and nDNA supported the existence of two main management units as the Mediterranean and Africa. Based on the mtDNA results, the Mediterranean can be divided into two subunits; western Turkey and the eastern Mediterranean.

Comparing and combining distance-based and character-based approaches for barcoding turtles

Molecular barcoding can serve as a powerful tool in wildlife forensics and may prove to be a vital aid in conserving organisms that are threatened by illegal wildlife trade, such as turtles (Order Testudines). We produced cytochrome oxidase subunit one (COI) sequences (650 bp) for 174 turtle species and combined these with publicly available sequences for 50 species to produce a data set representative of the breadth of the order. Variability within the barcode region was assessed, and the utility of both distance-based and character-based methods for species identification was evaluated. For species in which genetic material from more than one individual was available (n = 69), intraspecific divergences were 1.3% on average, although divergences greater than the customary 2% barcode threshold occurred within 15 species. High intraspecific divergences could indicate species with a high degree of internal genetic structure or possibly even cryptic species, although introgression is also probable in some of these taxa. Divergences between species of the same genus were 6.4% on average; however, 49 species were <2% divergent from congeners. Low levels of interspecific divergence could be caused by recent evolutionary radiations coupled with the low rates of mtDNA evolution previously observed in turtles. Complementing distance-based barcoding with character-based methods for identifying diagnostic sets of nucleotides provided better resolution in several cases where distance-based methods failed to distinguish species. An online identification engine was created to provide character-based identifications. This study constitutes the first comprehensive barcoding effort for this seriously threatened order.

Complete mitochondrial genome of the bullhead torrent catfish, Liobagrus obesus (Siluriformes, Amblycipididae): Genome description and phylogenetic considerations inferred from the Cyt b and 16S rRNA genes

Mitochondrial DNA (mtDNA) from the bullhead torrent catfish, Liobagrus obesus, was isolated by long-polymerase chain reaction (Long-PCR) with universal primers and was fully sequenced by primer working using flanking sequences. The complete mtDNA from L. obesus was 16,531 bp in length and contained 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and a control region, demonstrating a structure very similar to that of other bony fish. An analysis of the protein-coding genes revealed a statistically substantiated bias in (T+C): (A+G) content, supporting earlier findings regarding this peculiarity. As indicated by a chi-square test, the observed scores for pyrimidine and purine content were different from those expected assuming a 50:50 ratio: chi(2)=41.63, d.f.=5, p<0.000001 for three categories, including the 1st, 2nd, and 3rd codon positions. Further, there was a difference in nucleotide content between ND6 and the other 12 protein-coding genes in L. obesus. The values of p-distances, as summarized for different scales of evolutionary history at the Cyt b gene, revealed a clear pattern of increased nucleotide diversity at four levels: (1) intraspecies, (2) intragenus, (3) intrafamily, and (4) intraorder. Scores of average p-distances of the four categories in catfish were (1) 1.59+/-0.54%, (2) 5.28+/-1.72% (3) 16.37+/-1.26%, and (4) 19.81+/-0.14%, respectively. These data support the hypothesis that speciation in the order Siluriformes, in most cases, follows a geographic mode through the accumulation of a numerous small genetic changes over a long time period. A phylogenetic tree for the bullhead torrent catfish and several other fish species belonging to the order Siluriformes was developed on the basis of respective Cyt b sequences (1138 bp); the analysis revealed a monophyletic origin for the five examined families. A species-specific clustering of sequences from single species was obtained, supporting additionally basic phylogenetic information for the catfish and the barcoding suitability of Cyt b sequence data. Lastly, one of the well-supported properties of our phylogenetic tree (99% repetition level in our analysis) was the monophyletic placement of all catfish (order Siluriformes) among other ray-finned fish of the class Actinopterigii. Also discussed herein are the aspects of phylogeny based on the 16S rRNA gene.