The first complete mitochondrial genome of the Indian Tent Turtle, Pangshura tentoria (Testudines: Geoemydidae): Characterization and comparative analysis

A Genomic-Based Workflow for eDNA Assay Development for a Critically Endangered Turtle, Myuchelys georgesi

Environmental DNA (eDNA) analysis has become a popular conservation tool for detecting rare and elusive species. eDNA assays typically target mitochondrial DNA (mtDNA) due to its high copy number per cell and its ability to persist in the environment longer than nuclear DNA. Consequently, the development of eDNA assays has relied on mitochondrial reference sequences available in online databases, or in cases where such data are unavailable, de novo DNA extraction and sequencing of mtDNA. In this study, we designed eDNA primers for the critically endangered Bellinger River turtle (Myuchelys georgesi) using a bioinformatically assembled mitochondrial genome (mitogenome) derived from a reference genome. We confirmed the accuracy of this assembled mitogenome by comparing it to a Sanger-sequenced mitogenome of the same species, and no base pair mismatches were detected. Using the bioinformatically extracted mitogenome, we designed two 20 bp primers that target a 152-base-pair-long fragment of the cytochrome oxidase 1 (CO1) gene and a 186-base-pair-long fragment of the cytochrome B (CytB) gene. Both primers were successfully validated in silico, in vitro, and in situ.

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.

Complete mitogenome and intra-family comparative mitogenomics showed distinct position of Pama Croaker Otolithoides pama

The Pama Croaker, Otolithoides pama, is an economically important fish species in Bangladesh. Intra-family similarities in morphology and typical barcode sequences of cox1 create ambiguities in its identification. Therefore, morphology and the complete mitochondrial genome of O. pama, and comparative mitogenomics within the family Sciaenidae have been studied. Extracted genomic DNA was subjected to Illumina-based short read sequencing for De-Novo mitogenome assembly. The complete mitogenome of O. pama (Accession: OQ784575.1) was 16,513 bp, with strong AC biasness and strand asymmetry. Relative synonymous codon usage (RSCU) among 13 protein-coding genes (PCGs) of O. pama was also analyzed. The studied mitogenomes including O. pama exhibited consistent sizes and gene orders, except for the genus Johnius which possessed notably longer mitogenomes with unique gene rearrangements. Different genetic distance metrics across 30 species of Sciaenidae family demonstrated 12S rRNA and the control region (CR) as the most conserved and variable regions, respectively, while most of the PCGs undergone a purifying selection. Different phylogenetic trees were congruent with one another, where O. pama was distinctly placed. This study would contribute to distinguishing closely related fish species of Sciaenidae family and can be instrumental in conserving the genetic diversity of O. pama.

Phylogenetic relationships and adaptation in deep-sea carideans revealed by mitogenomes

The deep-sea environment is characterized by extreme and inhospitable conditions, including oxygen depletion, low temperatures, high pressure, absence of light, and limited food availability. Mitochondria and mitogenomes play a crudial role in aerobic respiration to generate energy for eukaryotes. Here, using the Illumina Hiseq 4000 platform, we performed mitogenome sequencing for five deep-sea caridean species: Lebbeus shinkaiae, Lebbeus Formosus, Glyphocrangon regalis, Heterocarpus dorsalis, and Heterocarpus laevigatus, and five deep-sea caridean mitogenomes were assembled and identified. Each of the five mitogenomes contained 13 protein-coding genes, 2 rRNAs and 22 tRNAs. Specific elements, such as tandem repeats and AT-rich sequences, were observed in the control regions of Lebbeus formosus and Lebbeus shinkaiae, potentially take a role in regulating mitochondrial genome replication and transcription. The gene order of all obtained mitogenomes follows caridean ancestral type organization. Phylogenetic analysis shows a robustly supported phylogenetic tree for the infraorder Caridea. The monophyly of the families included in this study was strongly supported. This study supports the monophyly of Oplophoroidea, but rejects the monophyletic status of Nematocarcinoidea, Crangonoidea, and Alpheoidea. At the genus level, Plesionika is polyphyletic and Rimicaris is paraphyletic in our analysis. Furthermore, Paralebbeus may be considered invalid and synonymous with Lebbeus. Positive selection analysis reveals evidence for adaptive changes in the mitogenome of different deep-sea caridean lineages. Nine residues located in cox1, cox3, atp6, nad1, nad2, nad4, nad5, nad6 and cytb were determined to have undergone positive selection. Mitogenome of different deep-sea lineages experienced different positive selection, and the lineage represented by Alvinocarididae living in deep-sea hydrothermal vents experienced the strongest positive selection. This study provides valuable insights into the adaptive evolution of deep-sea shrimps at the mitochondrial, highlighting the mitogenomic strategy that contribute to their unique adaptations in the deep-sea environment.

Matrilineal phylogeny and habitat suitability of the endangered spotted pond turtle (Geoclemys hamiltonii; Testudines: Geoemydidae): a two-dimensional approach to forecasting future conservation consequences

The spotted pond turtle (Geoclemys hamiltonii) is a threatened and less explored species endemic to Bangladesh, India, Nepal, and Pakistan. To infer structural variation and matrilineal phylogenetic interpretation, the present research decoded the mitogenome of G. hamiltonii (16,509 bp) using next-generation sequencing technology. The mitogenome comprises 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one AT-rich control region (CR) with similar strand symmetry in vertebrates. The ATG was identified as a start codon in most of the PCGs except Cytochrome oxidase subunit 1 (cox1), which started with the GTG codon. The non-coding CR of G. hamiltonii was determined to have a unique structure and variation in different domains and stem-loop secondary structure as compared with other Batagurinae species. The PCGs-based Bayesian phylogeny inferred strong monophyletic support for all Batagurinae species and confirmed the sister relationship of G. hamiltonii with Pangshura and Batagur taxa. We recommend generating more mitogenomic data for other Batagurinae species to confirm their population structure and evolutionary relationships. In addition, the present study aims to infer the habitat suitability and habitat quality of G. hamiltonii in its global distribution, both in the present and future climatic scenarios. We identify that only 58,542 km2 (7.16%) of the total range extent (817,341 km2) is suitable for this species, along with the fragmented habitats in both the eastern and western ranges. Comparative habitat quality assessment suggests the level of patch shape in the western range is higher (71.3%) compared to the eastern range. Our results suggest a massive decline of approximately 65.73% to 70.31% and 70.53% to 75.30% under ssp245 and ssp585 future scenarios, respectively, for the years between 2021-2040 and 2061-2080 compared with the current distribution. The present study indicates that proper conservation management requires greater attention to the causes and solutions to the fragmented distribution and safeguarding of this endangered species in the Indus, Ganges, and Brahmaputra (IGB) river basins.

Characterization of the complete mitochondrial genome of the critically endangered Southern River Terrapin, Batagur affinis affinis (Reptilia: Geoemydidae)

In this study, we report the nearly complete mitochondrial sequence of Batagur affinis affinis. The assembled mitogenome consists of 13 PCGs, 22 tRNA genes, two rRNAs and one near-complete D-loop region. Of the annotated genes, the ND6 subunit gene and eight tRNA genes were encoded on the L-strand, while the remaining genes were dispersed on the H-strand. Except for CO1, which has a GTG start codon, all protein-coding genes begin with ATG. The mitogenome has been deposited in NCBI GenBank under the accession number OQ409915. Phylogenetic tree analysis based on publicly available mitogenomes indicate the sister grouping of B. affinis affinis with B. kachuga.

Low coverage sequencing provides insights into the key features of the nuclear and mitochondrial genomes of the Alligator Snapping Turtle Macrochelys temminckii

The alligator snapping turtle Macrochelys temminckii is a culturally, ecologically, and evolutionary relevant species of conservation concern. In this study, we conducted a genome survey of M. temminckii. Using a low-coverage short read sequencing strategy, this study estimated the genome size, repetitive genome content, annotated and quantified repetitive elements, assembled the 45S rRNA DNA operon, and characterized in detail the mitochondrial genome of M. temminckii. Using a k-mer strategy, the estimated haploid genome size varied between 3.77 and 3.19 Gbp, which is within the range previously reported for other representatives of the family Chelydridae. Repetitive genome content estimates using different k-mers (21 to 51) indicated that more than 75% of the genome of M. temminckii comprised repetitive elements. Taking into account only annotated repetitive elements, the most common repetitive elements were classified as Class I - Long Interspersed Nuclear Element (LINE) which were more abundant than Class I - Penelope and Class I - Long Terminal Repeat (LTR) Ty3-gypsy mobile elements. Less abundant repeat element families in the nuclear genome of M. temminckii included Class I - DIRS mobile elements and Satellite DNA. The nuclear ribosomal operon was partially assembled into two contigs, one encoding the complete ssrDNA and a second comprising the full lsrDNA. The AT-rich complete mitochondrial genome was 16,570 bp long. These new genomic resources are of utmost importance to aid in the development of conservation plans for this freshwater turtle.

Mitochondriomics of Clarias Fishes (Siluriformes: Clariidae) with a New Assembly of Clarias camerunensis: Insights into the Genetic Characterization and Diversification

The mitogenome of an endemic catfish Clarias camerunensis was determined from the Cameroon water. This circular mitogenome was 16,511 bp in length and comprised 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a single AT-rich control region. The heavy strand accommodates 28 genes, whereas the light strand is constituted by ND6 and eight transfer RNA (tRNA) genes. The C. camerunensis mitochondrial genome is AT biased (56.89%), as showcased in other Clarias species. The comparative analyses revealed that most of the Clarias species have 6 overlapping and 11 intergenic spacer regions. Most of the PCGs were initiated and terminated with the ATG start codon and TAA stop codon, respectively. The tRNAs of C. camerunensis folded into the distinctive cloverleaf secondary structure, except trnS1. The placement of the conserved domains in the control region was similar in all the Clarias species with highly variable nucleotides in CSB-I. Both maximum likelihood and Bayesian-based matrilineal phylogenies distinctly separated all Clarias species into five clades on the basis of their known distributions (South China, Sundaland, Indochina, India, and Africa). The TimeTree analysis revealed that the two major clades (Indo-Africa and Asia) of Clarias species might have diverged during the Paleogene (≈28.66 MYA). Our findings revealed the separation of Indian species (C. dussumieri) and African species (C. camerunensis and Clarias gariepinus) took place during the Paleogene, as well as the South Chinese species (Clarias fuscus) and Sundaland species (Clarias batrachus) splits from the Indochinese species (Clarias macrocephalus) during the Neogene through independent colonization. This pattern of biotic relationships highlights the influence of topography and geological events in determining the evolutionary history of Clarias species. The enrichment of mitogenomic data and multiple nuclear loci from their native range or type locality will confirm the true diversification of Clarias species in African and Asian countries.

The First Complete Mitochondrial Genome of Eucrate crenata (Decapoda: Brachyura: Goneplacidae) and Phylogenetic Relationships within Infraorder Brachyura

Characterizing the complete mitochondrial genome (mitogenome) of an organism is useful for genomic studies in taxonomy and evolution. The mitogenomic characteristics of Eucrate crenata (Decapoda: Brachyura: Goneplacidae) have never been studied. The present study decodes the first mitogenome of E. crenata by high-throughput sequencing (HTS). The length of the mitogenome is 15,597 bp, and it contains 13 protein-coding genes, 2 ribosomal RNA genes (rrnS and rrnL), and 22 transfer RNA genes. There are 14 and 23 genes observed on the heavy and light strands, respectively. E. crenata possesses a trnH-cac translocation, with the trnH-cac shifted between trnE-gaa and trnF-ttc instead of the usual location between nad5 and nad4 in decapods. Phylogenetic analyses based on the current dataset of 33 Brachyuran mitogenomes indicate that E. crenata. is closely related to Ashtoret lunaris of Matutidae. The similar codon usage and rearrangements in the two species provide evidence for their close phylogenetic relationship. Positive selection analysis showed that one residue located in cox1 was identified as a positively selected site with high BEB value (>95%), indicating that this gene was under positive selection pressure. This study is the first complete mitogenome record for the family Goneplacidae, and the results obtained may improve the understanding of the phylogeny of Goneplacidae in Brachyura.

Characterization of the complete mitochondrial genome of the many-lined sun skink (Eutropis multifasciata) and comparison with other Scincomorpha species

Characterizating the complete mitochondrial genome (mitogenome) of an organism allows detailed genomic studies in systematics and evolution. The present study decodes the mitogenome (17,062 bp) of the many-lined sun skink, Eutropis multifasciata, using next-generation sequencing. To compare the diversity of mitogenomic structure and investigate intraspecific evolutionary relationships among the Asian Scincomorpha, the mitogenomes of 46 other species were examined concurrently. Within the group, the size of mitogenomes varied predominantly in the length at their control regions. The Ka/Ks ratios of 12 protein codon genes (PCGs) were lower than 1.00, demonstrating that they were under relaxed or moderate purifying selection. However, the ND5 had a Ka/Ks ratio >1, and was considered to be under positive selection. Currently there are two superfamilies in Scincomorpha (i.e. Scincoidea and Lacertoidea), but phylogenetic analysis using Bayesian Inference and Maximum-Likelihood Estimations produced phylogenetic trees with three clades in Scincomorpha ((Scincoidea + Lacertoidea (part)) + Gymnophthalmidae)).

The complete mitochondrial genome of the endangered Assam Roofed Turtle, Pangshura sylhetensis (Testudines: Geoemydidae): Genomic features and phylogeny

The Assam Roofed Turtle, Pangshura sylhetensis is an endangered and least studied species endemic to India and Bangladesh. The present study decodes the first complete mitochondrial genome of P. sylhetensis (16,568 bp) by using next-generation sequencing. The assembly encodes 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one control region (CR). Most of the genes were encoded on the majority strand, except NADH dehydrogenase subunit 6 (nad6) and eight tRNAs. All PCGs start with an ATG initiation codon, except for Cytochrome oxidase subunit 1 (cox1) and NADH dehydrogenase subunit 5 (nad5), which both start with GTG codon. The study also found the typical cloverleaf secondary structures in most of the predicted tRNA structures, except for serine (trnS1) which lacks of conventional DHU arm and loop. Both Bayesian and maximum-likelihood phylogenetic inference using 13 concatenated PCGs demonstrated strong support for the monophyly of all 52 Testudines species within their respective families and revealed Batagur trivittata as the nearest neighbor of P. sylhetensis. The mitogenomic phylogeny with other amniotes is congruent with previous research, supporting the sister relationship of Testudines and Archosaurians (birds and crocodilians). Additionally, the mitochondrial Gene Order (GO) analysis indicated plesiomorphy with the typical vertebrate GO in most of the Testudines species.