Comparison research and phylogenetic implications of mitochondrial control regions in four soft-shelled turtles of Trionychia (Reptilia, Testudinata)

Low depth sequencing reveals the critically endangered Batagur kachuga (Red-crowned roofed turtle) mitochondrial genome and its evolutionary implications

The Batagur kachuga (B. kachuga), commonly known as the Red-crowned roofed turtle, is a critically endangered species native to India and its neighboring countries like Bangladesh, and Nepal. The present study is the first report of the complete mitochondrial genome of B. kachuga (16,517 bp) construed via the next-generation sequencing (NGS) approach from eggshell DNA. There are 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), 13 protein-coding genes (PCGs), and one putative control region (CR/D-loop) in the mitogenome. The CR region from the current study reveals conserved TAS, CD, and CSB domains and two AT-rich tandem repeat regions. Most genes are encoded in the heavy strand except the NADH dehydrogenase subunit 6 (ND6) gene and seven tRNA genes. Most PCGs start with the initiation codon ATG, except the COI (Cytochrome Oxidase Subunit-I) gene, which starts with the GTG codon. The present investigation also predicts the distinctive cloverleaf structures of tRNAs except for tRNA-Ser1 and tRNA-Ser-2, which lack a DHU arm. The comparative analysis of Ka/Ks with other 33 species from Order Testudines, in relation to B. kachuga, revealed negative selection in most PCGs, indicating a process of preservation and purification that aids in eliminating undesirable or detrimental substitutes. Phylogenetic analysis of this species has been analysed using the complete mitogenome of 33 turtle species. The maximum likelihood phylogenetic tree strongly supports each family in different clades and also reveals a close relationship between the Pangashura and Batagur genera. Our study suggests the generation of genome-wide molecular data, in terms of mitogenomes, SNPs, and SSRs, is needed to improve the understanding of this species and their phylogenetics and evolutionary relationships, which will help to improve the conservation efforts of this species.

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.

Characterization of five complete Cyrtodactylus mitogenome structures reveals low structural diversity and conservation of repeated sequences in the lineage

Mitochondrial genomes (mitogenomes) of five Cyrtodactylus were determined. Their compositions and structures were similar to most of the available gecko lizard mitogenomes as 13 protein-coding, two rRNA and 22 tRNA genes. The non-coding control region (CR) of almost all Cyrtodactylus mitogenome structures contained a repeated sequence named the 75-bp box family, except for C. auribalteatus which contained the 225-bp box. Sequence similarities indicated that the 225-bp box resulted from the duplication event of 75-bp boxes, followed by homogenization and fixation in C. auribalteatus. The 75-bp box family was found in most gecko lizards with high conservation (55-75% similarities) and could form secondary structures, suggesting that this repeated sequence family played an important role under selective pressure and might involve mitogenome replication and the likelihood of rearrangements in CR. The 75-bp box family was acquired in the common ancestral genome of the gecko lizard, evolving gradually through each lineage by independent nucleotide mutation. Comparison of gecko lizard mitogenomes revealed low structural diversity with at least six types of mitochondrial gene rearrangements. Cyrtodactylus mitogenome structure showed the same gene rearrangement as found in most gecko lizards. Advanced mitogenome information will enable a better understanding of structure evolution mechanisms.

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.

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.

Phylogenetic relationships and divergence dates of softshell turtles (Testudines: Trionychidae) inferred from complete mitochondrial genomes

The softshell turtles (Trionychidae) are one of the most widely distributed reptile groups in the world, and fossils have been found on all continents except Antarctica. The phylogenetic relationships among members of this group have been previously studied; however, disagreements regarding its taxonomy, its phylogeography and divergence times are still poorly understood as well. Here, we present a comprehensive mitogenomic study of softshell turtles. We sequenced the complete mitochondrial genomes of 10 softshell turtles, in addition to the GenBank sequence of Dogania subplana, Lissemys punctata, Trionyx triunguis, which cover all extant genera within Trionychidae except for Cyclanorbis and Cycloderma. These data were combined with other mitogenomes of turtles for phylogenetic analyses. Divergence time calibration and ancestral reconstruction were calculated using BEAST and RASP software, respectively. Our phylogenetic analyses indicate that Trionychidae is the sister taxon of Carettochelyidae, and support the monophyly of Trionychinae and Cyclanorbinae, which is consistent with morphological data and molecular analysis. Our phylogenetic analyses have established a sister taxon relationship between the Asian Rafetus and the Asian Palea + Pelodiscus + Dogania + Nilssonia + Amyda, whereas a previous study grouped the Asian Rafetus with the American Apalone. The results of divergence time estimates and area ancestral reconstruction show that extant Trionychidae originated in Asia at around 108 million years ago (MA), and radiations mainly occurred during two warm periods, namely Late Cretaceous-Early Eocene and Oligocene. By combining the estimated divergence time and the reconstructed ancestral area of softshell turtles, we determined that the dispersal of softshell turtles out of Asia may have taken three routes. Furthermore, the times of dispersal seem to be in agreement with the time of the India-Asia collision and opening of the Bering Strait, which provide evidence for the accuracy of our estimation of divergence time. Overall, the mitogenomes of this group were used to explore the origin and dispersal route of Trionychidae and have provided new insights on the evolution of this group.

The mitochondrial genome of the Saunders's gull Chroicocephalus saundersi (Charadriiformes: Laridae) and a higher phylogeny of shorebirds (Charadriiformes)

The complete mitogenome of Chroicocephalus saundersi was characterized and compared with the 6 published Charadriiformes mitogenomes. The mitogenome of C. saundersi is a closed circular molecule 16,739 bp in size, and contains 37 genes and a control region. The AT and GC skews are positive and negative, respectively, and in agreement with those of the other Charadriiformes mitogenomes. The mitogenome of C. saundersi contains 3 start codons (ATG, GTG, and ATT), 4 stop codons (TAA, TAG, AGG, and AGA), and an incomplete stop codon (T-) in 13 PCGs. A codon usage analysis of all available Charadriiformes mitogenomes showed that the ATG (78%) and TAA (50.5%) were the most common start codon and stop codon, respectively. An unusual start codon, ATT, is commonly found in the ND3s of Charadriiformes mitogenomes, whereas the more common start codons, ATC and ATA, are rarely found. In all the Laridae species, one extra cytosine was inserted at position 174 in ND3. The control region of C. saundersi is 1180-bp long, with a nucleotide composition of 30.2% A, 28.6% T, 27.3% C, and 14.0% G. Variable numbers of tandem repeats (VNTRs) with nine copies of the 10 bp repeat sequence (AACAACAAAC) are found within the CSB domain of the control region. The ML/BI analyses, based on the amino acids of the 13 mitochondrial PCGs, strongly support the monophyly of the order Charadriiformes, with the suborder Lari considered sister to the Scolopaci, which is in turn a sister group to the suborder Charadrii.

Recombination and evolution of duplicate control regions in the mitochondrial genome of the Asian big-headed turtle, Platysternon megacephalum

Complete mitochondrial (mt) genome sequences with duplicate control regions (CRs) have been detected in various animal species. In Testudines, duplicate mtCRs have been reported in the mtDNA of the Asian big-headed turtle, Platysternon megacephalum, which has three living subspecies. However, the evolutionary pattern of these CRs remains unclear. In this study, we report the completed sequences of duplicate CRs from 20 individuals belonging to three subspecies of this turtle and discuss the micro-evolutionary analysis of the evolution of duplicate CRs. Genetic distances calculated with MEGA 4.1 using the complete duplicate CR sequences revealed that within turtle subspecies, genetic distances between orthologous copies from different individuals were 0.63% for CR1 and 1.2% for CR2app:addword:respectively, and the average distance between paralogous copies of CR1 and CR2 was 4.8%. Phylogenetic relationships were reconstructed from the CR sequences, excluding the variable number of tandem repeats (VNTRs) at the 3' end using three methods: neighbor-joining, maximum likelihood algorithm, and Bayesian inference. These data show that any two CRs within individuals were more genetically distant from orthologous genes in different individuals within the same subspecies. This suggests independent evolution of the two mtCRs within each P. megacephalum subspecies. Reconstruction of separate phylogenetic trees using different CR components (TAS, CD, CSB, and VNTRs) suggested the role of recombination in the evolution of duplicate CRs. Consequently, recombination events were detected using RDP software with break points at ≈290 bp and ≈1,080 bp. Based on these results, we hypothesize that duplicate CRs in P. megacephalum originated from heterological ancestral recombination of mtDNA. Subsequent recombination could have resulted in homogenization during independent evolutionary events, thus maintaining the functions of duplicate CRs in the mtDNA of P. megacephalum.