Two novel gene orders and the role of light-strand replication in rearrangement of the vertebrate mitochondrial genome

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.

The mitochondrial genome sequences of eleven leafhopper species of Batracomorphus (Hemiptera: Cicadellidae: Iassinae) reveal new gene rearrangements and phylogenetic implications

Batracomorphus is the most diverse and widely distributed genus of Iassinae. Nevertheless, there has been no systematic analysis of the genome structure and phylogenetic relationships of the genus. To determine the characteristics of the mitogenomes of Batracomorphus species as well as the phylogenetic relationships between them, we sequenced and compared the mitogenomes of 11 representative Batracomorphus species. The results revealed that the mitogenomes of the 11 Batracomorphus species exhibited highly similar gene and nucleotide composition, and codon usage compared with other reported mitogenomes of Iassinae. Of these 11 species, we found that the mitogenomes of four species were rearranged in the region from trnI-trnQ-trnM to trnQ-trnI-trnM, whereas the remaining species presented a typical gene order. The topologies of six phylogenetic trees were in agreement. Eurymelinae consistently formed paraphyletic groups. Ledrinae and Evacanthinae formed sister taxa within the same clade. Similarly, Typhlocybinae and Mileewinae consistently clustered together. All phylogenetic trees supported the monophyly of Iassinae, indicating its evolutionary distinctiveness while also revealing its sister relationship with Coelidiinae. Notably, the nodes for all species of the genus Batracomorphus were well supported and these taxa clustered into a large branch that indicated monophyly. Within this large branch, four Batracomorphus species with a gene rearrangement (trnQ-trnI-trnM) exhibited distinctive clustering, which divided the large branch into three minor branches. These findings expand our understanding of the taxonomy, evolution, genetics, and systematics of the genus Batracomorphus and broader Iassinae groups.

Integrative Taxonomy Revealed High Diversity of Hemiphyllodactylus Bleeker, 1860 (Squamata: Gekkonidae) and the Description of Three New Species from Yunnan Province, China

The karst landform in Yunnan Province, China, represents one of the most biodiverse regions for Hemiphyllodactylus. Previous research has revealed that the karst forests in this province host a greater diversity of Hemiphyllodactylus than previously acknowledged. However, substantial fundamental data essential for taxonomic and biogeographical studies are lacking. We conducted extensive surveys for Hemiphyllodactylus in the Yunnan Province that led to the discovery of three new species from Menglian Dai, Lahu, and Wa Autonomous County and Jinghong City based on morphological and genetic data. A phylogenetic reconstruction based on the ND2 gene (1038 bp) placed the three new species into clade 3 and clade 4 of Agung et al. The uncorrected genetic pairwise distance of the Menglian specimens were greater than 5.7%, and those of the Jinghong City specimens were greater than 5.2% and 8.5%, respectively. They could be distinguished from their congeners by body size, chin scales, internasal scales, ventral scales, dorsal scales, and the total number of femoral and precloacal pores. Furthermore, we update the distributional knowledge of the known species.

Full-Length Transcriptome Profiling of the Complete Mitochondrial Genome of Sericothrips houjii (Thysanoptera: Thripidae: Sericothripinae) Featuring Extensive Gene Rearrangement and Duplicated Control Regions

The mitochondrial genome (mitogenome) of Thysanoptera has extensive gene rearrangement, and some species have repeatable control regions. To investigate the characteristics of the gene expression, transcription and post-transcriptional processes in such extensively gene-rearranged mitogenomes, we sequenced the mitogenome and mitochondrial transcriptome of Sericothrips houjii to analyze. The mitogenome was 14,965 bp in length and included two CRs contains 140 bp repeats between COIII-trnN (CR1) and trnT-trnP (CR2). Unlike the putative ancestral arrangement of insects, S. houjii exhibited only six conserved gene blocks encompassing 14 genes (trnL2-COII, trnD-trnK, ND2-trnW, ATP8-ATP6, ND5-trnH-ND4-ND4L and trnV-lrRNA). A quantitative transcription map showed the gene with the highest relative expression in the mitogenome was ND4-ND4L. Based on analyses of polycistronic transcripts, non-coding RNAs (ncRNAs) and antisense transcripts, we proposed a transcriptional model of this mitogenome. Both CRs contained the transcription initiation sites (TISs) and transcription termination sites (TTSs) of both strands, and an additional TIS for the majority strand (J-strand) was found within antisense lrRNA. The post-transcriptional cleavage processes followed the "tRNA punctuation" model. After the cleavage of transfer RNAs (tRNAs), COI and ND3 matured as bicistronic mRNA COI/ND3 due to the translocation of intervening tRNAs, and the 3' untranslated region (UTR) remained in the mRNAs for COII, COIII, CYTB and ND5. Additionally, isoform RNAs of ND2, srRNA and lrRNA were identified. In summary, the relative mitochondrial gene expression levels, transcriptional model and post-transcriptional cleavage process of S. houjii are notably different from those insects with typical mitochondrial gene arrangements. In addition, the phylogenetic tree of Thripidae including S. houjii was reconstructed. Our study provides insights into the phylogenetic status of Sericothripinae and the transcriptional and post-transcriptional regulation processes of extensively gene-rearranged insect mitogenomes.

The description of the first rock-dwelling species of butterfly lizard Leiolepis Cuvier, 1829 (Squamata, Agamidae) from the Khorat Plateau in northeastern Thailand

A new species of rock-dwelling Leiolepis is described from the Khorat Plateau in northeastern Thailand. Leiolepisglaurung sp. nov. can be differentiated from all other sexual species of Leiolepis by a combination of having a black gular region with a wide medial yellow stripe, a yellow ventrum with black mottling, bright red to orange subcaudal coloration, having reduced to no expandable flanks, and having only one black transverse bar on the flanks. This is the first rocky habitat-adapted Leiolepis. Leiolepisglaurung sp. nov. demonstrates numerous ecological adaptations to survive in these rocky habitats. Leiolepis are known for their expandable flanks with bright display colors, however Leiolepisglaurung sp. nov. has reduced or no ability to expand its flanks. We hypothesize this is an adaptation to reduce their body diameter to better fit into smaller rocky burrows unlike the larger and deeper burrows constructed in looser soils by other Leiolepis species. This discovery increases the number of Leiolepis species in Thailand to six, and worldwide to 11.

Dragons in the tropics - Phylogeography and speciation in Diporiphora lizards and common geographic breaks in co-distributed taxa

Co-distributed taxa can respond both similarly or differently to the same climatic and geological events, resulting in a range of phylogeographic patterns across the region. Using a nested approach on a taxonomically diverse yet morphologically conservative group of agamid lizards, we first aimed to evaluate more precisely the extent of phylogeographic structuring within the genus. Then, focusing on four lineages within the more widespread species, we assessed the impact of biogeographic barriers on phylogeographic structuring and demographic history of species, comparing to patterns previously observed in co-distributed taxa. These species occur in the Australian Monsoonal Tropics, a vast tropical savanna system with high richness and endemism associated with environmental heterogeneity and past climate fluctuations. The employment of genomic data helped to determine the relationships between specific taxa that were previously difficult to place. We found a local influence of biogeographic and climatic breaks on population dynamics, analogous to other species. We detected high levels of population structure in the West Kimberley and Arnhem Plateau, which are already known for high endemism. However, we also highlighted unique lineages in areas that have been overlooked until recently, in the South Kimberley and West Top End. Climatic and geographical features in the Arnhem Plateau act as a soft barrier between populations in the east and west regions of the Top End. These observations reflect patterns observed for other vertebrates across this rich biome, indicating how climatic variation, species' ecology, and landscape features interact to shape regional diversity and endemism.

A new species of Hemiphyllodactylus (Squamata, Gekkonidae) from southwest Yunnan, China

A new species of gekkonid, Hemiphyllodactylusgengmaensissp. nov., is described based on six specimens from Gengma Dai and Wa Autonomous County, Yunnan, China. The new species can be distinguished from its congeners by a significant genetic divergence of greater than 9.7% in the mitochondrial ND2 gene and a combination of the following characters: a maximum SVL of 43.24mm; 8 or 9 chin scales; six circumnasal scales; 2 or 3 internasal scales; 9-11 supralabial scales; 8 or 9 infralabial scales; 11-18 dorsal scales; 8-10 ventral scales; a manual lamellar formula of 5-5-5-4 or 5-6-5-4 and a pedal lamellar formula of 5-5-6-5; 20-25 precloacal and femoral pore-bearing scales contiguous in males; dark postorbital stripes or striping on body; dark dorsal transverse blotches present; and a brown postsacral mark bearing anteriorly projecting arms. The discovery of this new species brings the number of Hemiphyllodactylus species in China to 15.

Two new species of the Cnemaspisgalaxia complex (Squamata, Gekkonidae) from the eastern slopes of the southern Western Ghats

Two new species allied to Cnemaspisgalaxia are described from the eastern slopes of the south Western Ghats, Tamil Nadu, India. Both new species are members of the ornata subclade within the beddomei clade. The two new species can be easily distinguished from all other members of the beddomei clade and each other by a combination of nonoverlapping morphological characters such as small body size, distinct colouration of both sexes, the number of dorsal tubercles around the body, the number or arrangement of paravertebral tubercles, the number of midventral scales across the belly and longitudinal ventral scales from mental to cloaca, besides uncorrected pairwise ND2 and 16S sequence divergence of ≥ 7.4% and ≥ 2.7%. The two new species are distributed from low elevation, deciduous forests of Srivilliputhur, and add to the five previously known endemic vertebrates from Srivilliputhur-Megamalai Tiger Reserve.

A new species of the Cyrtodactyluschauquangensis species group (Squamata, Gekkonidae) from Lao Cai Province, Vietnam

We describe a new species of the genus Cyrtodactylus based on five adult specimens from Bac Ha District, Lao Cai Province, northern Vietnam. Cyrtodactyluslucisp. nov. is distinguished from the remaining Indochinese bent-toed geckos by a combination of the following morphological characteristics: medium size (SVL up to 89.5 mm); dorsal tubercles in 17-19 irregular transverse rows; ventral scales in 32-34 longitudinal rows at midbody; precloacal pores present in both sexes, 9 or 10 in males, 8 or 9 in females; 12-15 enlarged femoral scales on each thigh; femoral pores 9-12 in males, 5-10 in females; postcloacal tubercles 2-4; lamellae under toe IV 21-23; dorsal pattern consisting of 5 or 6 irregular dark bands, a thin neckband without V-shape or triangle shape in the middle, top of head with dark brown blotches; subcaudal scales transversely enlarged. Molecular phylogenetic analyses recovered the new species as the sister taxon to C.gulinqingensis from Yunnan Province, China, with strong support from all analyses and the two taxa are separated by approximately 8.87-9.22% genetic divergence based on a fragment of the mitochondrial ND2 gene. This is the first representative of Cyrtodactylus known from Lao Cai Province.

A comprehensive phylogenomic study unveils evolutionary patterns and challenges in the mitochondrial genomes of Carcharhiniformes: A focus on Triakidae

The complex evolutionary patterns in the mitochondrial genome (mitogenome) of the most species-rich shark order, the Carcharhiniformes (ground sharks) has led to challenges in the phylogenomic reconstruction of the families and genera belonging to the order, particularly the family Triakidae (houndsharks). The current state of Triakidae phylogeny remains controversial, with arguments for both monophyly and paraphyly within the family. We hypothesize that this variability is triggered by the selection of different a priori partitioning schemes to account for site and gene heterogeneity within the mitogenome. Here we used an extensive statistical framework to select the a priori partitioning scheme for inference of the mitochondrial phylogenomic relationships within Carcharhiniformes, tested site heterogeneous CAT + GTR + G4 models and incorporated the multi-species coalescent model (MSCM) into our analyses to account for the influence of gene tree discordance on species tree inference. We included five newly assembled houndshark mitogenomes to increase resolution of Triakidae. During the assembly procedure, we uncovered a 714 bp-duplication in the mitogenome of Galeorhinus galeus. Phylogenetic reconstruction confirmed monophyly within Triakidae and the existence of two distinct clades of the expanded Mustelus genus. The latter alludes to potential evolutionary reversal of reproductive mode from placental to aplacental, suggesting that reproductive mode has played a role in the trajectory of adaptive divergence. These new sequences have the potential to contribute to population genomic investigations, species phylogeography delineation, environmental DNA metabarcoding databases and, ultimately, improved conservation strategies for these ecologically and economically important species.

The extraordinary rearrangement of mitochondrial genome of the wheat pest, Aptinothrips stylifer and the mitochondrial phylogeny of Thripidae (Thysanoptera)

The mitochondrial gene order in Thysanoptera is notably distinct and highly rearranged, with each species exhibiting its own unique arrangement. To elucidate the relationship between gene rearrangements and phylogeny, the complete mitochondrial genome (mitogenome) of the wheat pest, Aptinothrips stylifer, was sequenced and assembled, spanning a total length of 16,033 bp. Compared with the ancestral arthropod mitogenome, significant rearrangement differences were evident in A. stylifer, whereas the gene order between A. stylifer and Anaphothrips obscurus was similar. Phylogenetic trees were reconstructed based on all 13 protein-coding gene sequences using Bayesian inference and maximum-likelihood methods, both yielding similar topological structures. Notably, A. stylifer was robustly clustered with A. obscurus, affirming its classification within Anaphothrips genus group. This exemplifies the potential correlation between gene rearrangements and phylogeny in the Thripidae family. Additionally, the mitogenome of A. stylifer exhibited several atypical features, including: (1) Three putative control regions (CRs) in close proximity, with CR2 and CR3 displaying partial similarity, and CR1 differing in base composition; (2) Two transfer RNAs (tRNAs), trnS1 and trnV, lacking the DHU arm; (3) Two ribosomal RNA (rRNA) genes inverted and positioned distant from each other; (4) Negative AT and GC skew (AT skew = -0.001, GC skew = -0.077); (5) One transposition (nad6), one inverse transposition (trnQ), four inversions (trnF, trnH, trnC, and gene block nad1-trnL1-rrnL-trnV-rrnS), and four tandem duplication random loss events; and (6) Two protein-coding genes, nad2 and atp8, terminated with an incomplete stop codon "T".

Two new species of Cnemaspis Strauch, 1887 (Squamata: Gekkonidae) from southern India

We describe two new species of Cnemaspis from Alagarkovil, Madurai District and Mekkarai, Shenkottai, Tenkasi, southern India, based on molecular and morphological data. The new species differ from their congeners by genetic divergences of 4.816.1% and 14.831.6% for the ND2 mitochondrial gene. Two new species belong to the Cnemaspis gracilis and Cnemaspis beddomei clades, respectively. Molecular data for the newly described Cnemaspis aaronbaueri from its type locality is also presented, and its phylogenetic position is established.

Its all in the name: Another new Cyrtodactylus Gray (Squamata: Gekkonidae) from northern Mizoram, North-east India

A new species of Cyrtodactylus is described from Vairengte town, situated in the Kolasib District of Mizoram State, north-eastern India. The new species is found to be a member of Indo-Burman Cyrtodactylus khasiensis clade based on ND2 gene sequences and morphological parameters, such as number of precloacal pores, mid-ventral scale rows, paravertebral tubercles on the trunk, dorsal tubercle rows, subdigital lamellae on pes and subcaudal scalation, making it the sixth endemic Cyrtodactylus from Mizoram and twenty second from north-east India. Moreover, phylogenetic evidence suggests the new species to be sister to the recently described C. aaronbaueri, and morphological analyses also reveal marginal separation between the two species based on the PCA of infralabials, lamellae on fingers and toes, paravertebral tubercles on the trunk, and dorsal tubercle rows.

Another step through the crux: a new microendemic rock-dwelling Paroedura (Squamata, Gekkonidae) from south-central Madagascar

Using an integrative taxonomic approach including genetic and morphological data, we formally describe a new microendemic gecko species belonging to the Paroedurabastardi clade, previously referred to as P.bastardi D. We name this taxon currently known from Anja Reserve and Tsaranoro Valley Forest (south-central Madagascar), as P.manongavatosp. nov. The new species differs from other species of the P.bastardi clade by ≥ 12.4% uncorrected p-distance at the mitochondrial 16S rRNA gene and it forms a monophyletic group in the COI mtDNA phylogenetic tree. It lacks haplotype sharing at the nuclear KIAA1239 and CMOS genes with the other species of the same complex, including the syntopic P.rennerae. Given its limited extent of occurrence and high levels of habitat fragmentation linked to forest clearances and fires, we propose the IUCN Red List Category of Critically Endangered, based on the B1ab(iii) criterion. The conservation value of Anja Reserve and Tsaranoro Valley Forest is remarkable. Preserving the remaining deciduous forest habitat is of paramount importance to protect these narrow-range reptile species.

Two-Fold ND5 Genes, Three-Fold Control Regions, lncRNA, and the "Missing" ATP8 Found in the Mitogenomes of Polypedates megacephalus (Rhacophridae: Polypedates)

In prior research on the mitochondrial genome (mitogenome) of Polypedates megacephalus, the one copy of ND5 gene was translocated to the control region (CR) and the ATP8 gene was not found. Gene loss is uncommon among vertebrates. However, in this study, we resequenced the mitogenomes of P. megacephalus from different regions using a "primer bridging" approach with Sanger sequencing technologies, which revealed the "missing" ATP8 gene in P. megacephalus as well as three other previously published Polypedates. The mitogenome of this species was found to contain two copies of the ND5 genes and three copies of the control regions. Furthermore, multiple tandem repeats were identified in the control regions. Notably, we observed that there was no correlation between genetic divergence and geographic distance. However, using the mitogenome, gene expression analysis was performed via RT-qPCR of liver samples and it was thus determined that COIII, ND2, ND4, and ND6 were reduced to 0.64 ± 0.24, 0.55 ± 0.34, 0.44 ± 0.21 and 0.65 ± 0.17, respectively, under low-temperature stress (8 °C) as compared with controls (p < 0.05). Remarkably, the transcript of long non-coding RNA (lncRNA) between positions 8029 and 8612 decreased significantly with exposure to low-temperature stress (8 °C). Antisense ND6 gene expression showed a downward trend, but this was not significant. These results reveal that modulations of protein-coding mitochondrial genes and lncRNAs of P. megacephalus play a crucial role in the molecular response to cold stress.

Comparative Mitogenomic Analyses of Darkling Beetles (Coleoptera: Tenebrionidae) Provide Evolutionary Insights into tRNA-like Sequences

Tenebrionidae is widely recognized owing to its species diversity and economic importance. Here, we determined the mitochondrial genomes (mitogenomes) of three Tenebrionidae species (Melanesthes exilidentata, Anatolica potanini, and Myladina unguiculina) and performed a comparative mitogenomic analysis to characterize the evolutionary characteristics of the family. The tenebrionid mitogenomes were highly conserved with respect to genome size, gene arrangement, base composition, and codon usage. All protein-coding genes evolved under purifying selection. The largest non-coding region (i.e., control region) showed several unusual features, including several conserved repetitive fragments (e.g., A+T-rich regions, G+C-rich regions, Poly-T tracts, TATA repeat units, and longer repetitive fragments) and tRNA-like structures. These tRNA-like structures can bind to the appropriate anticodon to form a cloverleaf structure, although base-pairing is not complete. We summarized the quantity, types, and conservation of tRNA-like sequences and performed functional and evolutionary analyses of tRNA-like sequences with various anticodons. Phylogenetic analyses based on three mitogenomic datasets and two tree inference methods largely supported the monophyly of each of the three subfamilies (Stenochiinae, Pimeliinae, and Lagriinae), whereas both Tenebrioninae and Diaperinae were consistently recovered as polyphyletic. We obtained a tenebrionid mitogenomic phylogeny: (Lagriinae, (Pimeliinae, ((Tenebrioninae + Diaperinae), Stenochiinae))). Our results provide insights into the evolution and function of tRNA-like sequences in tenebrionid mitogenomes and contribute to our general understanding of the evolution of Tenebrionidae.

A new species of the genus Hemiphyllodactylus Bleeker, 1860 (Squamata: Gekkonidae) from a lowland monsoon forest in Southern Vietnam

An integrative taxonomic analysis delimits and diagnoses the new species Hemiphyllodactylus cattien sp. nov. from Cat Tien National Park, Dong Nai Province, southern Vietnam-the eighth species of Hemiphyllodactylus to be described from Vietnam. The new species is recovered as the sister species to a clade comprised of H. indosobrinus from Laos and H. flaviventris from Thailand based on phylogenetic analyses using ND2 mtDNA gene sequences. Unlike the vast majority of Hemiphyllodactylus which are site-restricted upland endemics, Hemiphyllodactylus cattien sp. nov. is a lowland species collected from one of the most well-surveyed regions in Vietnam. Its discovery further emphasizes the need for continued re-sampling of survey sites already considered to be well-explored.

An expanded description of Cnemaspis gracilis (Beddome 1870) (Squamata: Gekkonidae) based on recent material

Cnemaspis gracilis is a poorly known species that has been reported from across southern India, with 11 named species within the clade. The species is known only from the type series, as many other records from areas outside the type locality have turned out to represent distinct species. We provide an expanded morphological description of the species based on topotypic material and other samples that were assigned using genetic data. Cnemaspis gracilis is now known with certainty from three localities in the vicinity of Palakkad, Palakkad District, Kerala, and Valparai town, Coimbatore district, Tamil Nadu, besides from an introduced population in the Timber Market, Kolhapur city, Kolhapur district, Maharashtra, India. The species is 7.6-15.8 % divergent on the ND2 gene from described members of the gracilis clade and can be diagnosed by a combination of morphological characters including body size, number of paravertebral tubercles between limb insertions, number of dorsal tubercle rows, number of ventral scale rows across the belly, number of femoral and precloacal pores and poreless scales separating these series in males, and the presence of a single central black ocellus on the neck and a smaller one on the occiput. The assignment of extant populations to the species, a detailed morphological description and genetic data will allow the description of many more species from within this diverse clade of diurnal geckos.

Hemidactylus paucifasciatus (Squamata: Gekkonidae), a new species of large-bodied, tuberculate gecko from Northern Odisha, India

A new Hemidactylus belonging to the H. prashadi group is described herein, based on a series of four specimens from the Chota Nagpur Plateau in northern Odisha. Hemidactylus paucifasciatus sp. nov. can be diagnosed from regional congeners by its large body size (SVL up to 117.7 mm), the number of dorsal tubercle rows at mid-body, three or four bands on the dorsum, the number of enlarged tubercles in paravertebral rows, the number of femoral pores and poreless scales separating the left and right series on the femoral-precloacal row in males, and the number of ventral scales across the belly at mid-body. A mitochondrial phylogeny demonstrates that the new species is nested within the H. triedrus complex of the northern clade within the H. prashadi group, with 16.2-16.5 % pairwise divergence in ND2 sequence data from members of the H. triedrus complex and > 21 % from other members of the prashadi group. Hemidactylus paucifasciatus sp. nov. is the first endemic gecko to be described from the Chota Nagpur Plateau.

Detecting gene breakpoints in noisy genome sequences using position-annotated colored de-Bruijn graphs

BACKGROUND

Identifying the locations of gene breakpoints between species of different taxonomic groups can provide useful insights into the underlying evolutionary processes. Given the exact locations of their genes, the breakpoints can be computed without much effort. However, often, existing gene annotations are erroneous, or only nucleotide sequences are available. Especially in mitochondrial genomes, high variations in gene orders are usually accompanied by a high degree of sequence inconsistencies. This makes accurately locating breakpoints in mitogenomic nucleotide sequences a challenging task.

RESULTS

This contribution presents a novel method for detecting gene breakpoints in the nucleotide sequences of complete mitochondrial genomes, taking into account possible high substitution rates. The method is implemented in the software package DeBBI. DeBBI allows to analyze transposition- and inversion-based breakpoints independently and uses a parallel program design, allowing to make use of modern multi-processor systems. Extensive tests on synthetic data sets, covering a broad range of sequence dissimilarities and different numbers of introduced breakpoints, demonstrate DeBBI 's ability to produce accurate results. Case studies using species of various taxonomic groups further show DeBBI 's applicability to real-life data. While (some) multiple sequence alignment tools can also be used for the task at hand, we demonstrate that especially gene breaks between short, poorly conserved tRNA genes can be detected more frequently with the proposed approach.

CONCLUSION

The proposed method constructs a position-annotated de-Bruijn graph of the input sequences. Using a heuristic algorithm, this graph is searched for particular structures, called bulges, which may be associated with the breakpoint locations. Despite the large size of these structures, the algorithm only requires a small number of graph traversal steps.

The identity of Tokays (Squamata: Gekkonidae: Gekko) from Hong Kong

The Tokay gecko (Gekko gecko), found throughout Southeast Asia and India, is a heavily harvested species of high commercial value. Recent studies have supported the elevation of the black morph of Tokay gecko, found only in southern China and northern Vietnam, to the species Gekko reevesii. Previous genetic studies focused on specimens of G. reevesii from southwestern populations. Hong Kong, in southeastern China, has native populations of G. reevesii. To verify the identity of G. reevesii in Hong Kong, we employed three mitochondrial genes (COI, Cytb & ND2) and constructed a matrilineal genealogy using other specimens from Guangxi (southwestern China) and northern Vietnam, as well as G. gecko from a wide range of Southeast Asian countries. Our study confirmed that G. reevesii occurs naturally in Hong Kong, but one exotic population of G. gecko, likely a translocation from international trade, was also revealed. Our study did not reject the species ranking of G. reevesii. Moreover, like previous studies, we recovered a paraphyletic G. gecko, which may reflect a species complex, hybridization or incomplete lineage sorting. More extensive sampling of the two species over a broader range of their asserted distribution together with the use of both mitochondrial and nuclear DNA are required to better investigate their biogeography.

Two new species of Dixonius from Vietnam and Laos with a discussion of the taxonomy of Dixonius (Squamata, Gekkonidae)

Integrated analyses using maximum likelihood (ML), Bayesian inference (BI), principal component analysis (PCA), discriminate analysis of principal components (DAPC), multiple factor analysis (MFA), and analysis of variance (ANOVA) recovered two new diagnosable species of gekkonid lizards in the genus Dixonius, one from the Central Highlands, Gia Lai Province, Vietnam and another from the Vientiane Province, Laos. Phylogenetic analyses based on the mitochondrial NADH dehydrogenase subunit 2 gene (ND2) and adjacent tRNAs showed that Dixoniusgialaiensissp. nov. is the sister species of D.minhlei from Dong Nai Province, Vietnam and is nested within a clade that also includes the sister species D.siamensis and D.somchanhae. Dixoniusmuangfuangensissp. nov. is the sister species to D.lao from Khammouane Province, Laos and is embedded in a clade with D.vietnamensis, D.taoi, and undescribed species from Thailand. Multivariate (PCA, DAPC, and MFA) and univariate (ANOVA) analyses using combinations of 15 meristic (scale counts), six morphometric (measurements), and five categorical (color pattern and morphology) characters from 44 specimens encompassing all eight species of Dixonius from Vietnam and Laos clearly illustrate Dixoniusgialaiensissp. nov. and Dixoniusmuangfuangensissp. nov. are statistically different and discretely diagnosable from all closely related species of Dixonius. These integrative analyses also highlight additional taxonomic issues that remain unresolved within Dixonius and the need for additional studies. The discovery of these new species further emphasizes the underappreciated herpetological diversity of the genus Dixonius and illustrates the continued need for field work in these regions.

Rapid evolution, rearrangements and whole mitogenome duplication in the Australian stingless bees Tetragonula (Hymenoptera: Apidae): A steppingstone towards understanding mitochondrial function and evolution

The extreme conservation of mitochondrial genomes in metazoans poses a significant challenge to understanding mitogenome evolution. However, the presence of variation in gene order or genome structure, found in a small number of taxa, can provide unique insights into this evolution. Previous work on two stingless bees in the genus Tetragonula (T. carbonaria and T. hockingsi) revealed highly divergent CO1 regions between them and when compared to the bees from the same tribe (Meliponini), indicating rapid evolution. Using mtDNA isolation and Illumina sequencing, we elucidated the mitogenomes of both species. In both species, there has been a duplication of the whole mitogenome to give a total genome size of 30,666 bp in T. carbonaria; and 30,662 bp in T. hockingsi. These duplicated genomes present a circular structure with two identical and mirrored copies of all 13 protein coding genes and 22 tRNAs, with the exception of a few tRNAs that are present as single copies. In addition, the mitogenomes are characterized by rearrangements of two block of genes. We believe that rapid evolution is present in the whole Indo-Malay/Australasian group of Meliponini but is extraordinarily elevated in T. carbonaria and T. hockingsi, probably due to founder effect, low effective population size and the mitogenome duplication. All these features - rapid evolution, rearrangements, and duplication - deviate significantly from the vast majority of the mitogenomes described so far, making the mitogenomes of Tetragonula unique opportunities to address fundamental questions of mitogenome function and evolution.

A phylogenetic and taxonomic assessment of the Cnemaspis alwisi group (Reptilia: Gekkonidae) in Sri Lanka with a description of two new species from isolated misty-mountains

Sri Lanka is a local hotspot for Cnemaspis day geckos with 40 currently known species with 100% endemism. In this paper, we evaluate the phylogenetic relationships of Cnemaspis species belonging to the alwisi group of the podihuna clade and describe two additional new species of Cnemaspis from Sri Lanka; one from Galgiriya mountain, Kurunegala District, and another from Ethagala mountain, Ampara District. These new species were recorded from granite caves within forested areas in isolated mountains in the dry bioclimatic zone (point-endemics). Both new species are microhabitat specialists with narrow niches limited to humid, cool, canopy-shaded granite caves and old buildings associated with granite caves, where they are camouflaged by their cryptic morphology and body colouration. Furthermore, both species prefer narrow (~ 6–12 mm), long (~ 120–450 mm) and deep (~ 80–260 mm) crevices as refugia. The regions in which these habitats are located receive relatively low annual rainfall (1,000–1,500 mm). These new species are medium in size (28.5–36.8 mmSVL) and can be differentiated from all other Sri Lankan Cnemaspis by the presence of clearly enlarged, subhexagonal subcaudal scales and the absence of precloacal pores in males. Both species described here are categorised herein as Critically Endangered (CR) under the IUCN Red List criteria. The major threats for these new species are habitat loss due to expansion of commercial-scale agriculture, illicit forest encroachments, and forest fires. Therefore, we recommend that relevant authorities take immediate conservation action to ensure the protection of these forest areas with their buffer zones in the near future.

Aliens on Boats? The Eastern and Western Expansion of the African House Gecko

Invasive species disrupt relations between endemics and their ecosystem and are an increasing biodiversity conservation problem. The Hemidactylus genus comprises the most successful invasive reptile species, including the worldwide-distributed Hemidactylus mabouia. In this study, we used 12S and ND2 sequences to taxonomically identify and tentatively determine the diversity and origin of these invaders in Cabo Verde while also clarifying this for several Western Indian Ocean (WIO) populations. By comparing our sequences to recently published ones, we showed, for the first time, that Cabo Verde individuals belong to the H. mabouia sensu stricto lineage and that both of its sublineages (a and b) occur there. Both haplotypes are also in Madeira, which indicates a connection between these archipelagos, possibly related to the past Portuguese trading routes. Across the WIO, results clarified the identity of many island and coastal populations, showing that this likely invasive H. mabouia lineage is widespread in the region, including northern Madagascar, with important conservation implications. Colonisation origins were difficult to access due to the wide geographical spread of these haplotypes; thus, several possible scenarios were outlined. The introduction of this species throughout western and eastern Africa may threaten endemic taxa and needs to be closely monitored.

The complete mitochondrial genome and novel gene arrangement in Nesodiprion zhejiangensis Zhou & Xiao (Hymenoptera: Diprionidae)

The complete mitochondrial genome (mitogenome) of the sawfly, Nesodiprion zhejiangensis Zhou & Xiao, was sequenced, assembled, and deposited in GenBank (Accession Number: OM501121). The 15,660 bp N. zhejiangensis mitogenome encodes for 2 ribosomal RNAs (rrnL and rrnS), 22 transfer RNAs (tRNAs), 13 protein-coding genes (PCGs), and an AT-rich region of 450 bp in length. The nucleotide composition is biased toward adenine and thymine (A + T = 81.8%). Each PCG is initiated by an ATN codon, except for cox2, which starts with a TTG. Of 13 PCGs, 9 have a TAA termination codon, while the remainder terminate with a TAG or a single T. All tRNAs have the classic cloverleaf structure, except for the dihydrouridine (DHU) arm of tRNA^val, which forms a simple loop. There are 49 helices belonging to 6 domains in rrnL and 30 helices belonging to 4 domains in rrnS. In comparison to the ancestral architecture, N. zhejiangensis has the most rearranged mitogenome in Symphyta, in which rearrangement events of local inversion and transposition are identified in three gene clusters. Specifically, the main hotspot of gene rearrangement occurred between rrnS and trnY, and rearranged from rrnS-(AT-rich region)-I-Q-M-nd2-W-C-Y to rrnS-Q-W-C-nd2-I-M-(AT-rich region)-Y, involving a local inversion event of a large gene cluster and transposition events of some tRNAs. Transposition of trnA and trnR (rearranged from A-R to R-A) was observed at the nd3-nd5 gene junction while shuffling of trnP and trnT (rearranged from T-P to P-T) occurred at the nd4l-nd6 gene junction. While illegitimate inter-mtDNA recombination might explain the opposite orientations of transcription between rrnS and trnY, transposition events of tRNA in some gene blocks can be accounted for by the tandem duplication/random loss (TDRL) model. Our phylogenetic analysis suggests that N. zhejiangensis is closely related to congeneric species N. biremis and N. japonicus, which together form a sister lineage with the European pine sawfly, Neodiprion sertifer.

Scratching the surface: a new species of Bent-toed gecko (Squamata, Gekkonidae, Cyrtodactylus) from Timor-Leste of the darmandvillei group marks the potential for future discoveries

A new species of limestone-dwelling Bent-toed gecko (genus Cyrtodactylus) is described from Nino Konis Santana National Park in the far-east region of Timor-Leste. Both genetic and morphological data strongly support the evolutionary distinctness of the new species, which we describe herein as Cyrtodactylussantana sp. nov. Phylogenetic analysis based on the ND2 mitochondrial gene inferred the new species as part of the C.darmandvillei group with close genetic affinities to C.batucolus, C.seribuatensis, C.petani, C.sadleiri, and two undescribed lineages from the Moluccas in Indonesia. The new species represents the first species of Cyrtodactylus identified at the species level from Timor-Leste and fills an important gap in our understanding of the biogeography and evolutionary history of Cyrtodactylus especially in the Wallacean region. Our results strongly suggest that the diversity of Cyrtodactylus in Wallacea is still underestimated and many more unnamed species remain to be described.

Comparative analyses of the complete mitochondrial genomes of two southern African endemic guitarfish, Acroteriobatus annulatus and A. blochii

Shark-like rays (order Rhinopristiformes) are among the most threatened cartilaginous fish globally. Despite this, unresolved taxonomic issues still exist within the group. To date, no studies have used complete mitochondrial genomes to assess the phylogenetic placement of Acroteriobatus within the non-monophyletic family Rhinobatidae. The current study reports the first complete mitochondrial genomes for Acroteriobatus annulatus and A. blochii. Similar to other rhinopristiforms, the complete sequences of A. annulatus (16,773 bp) and A. blochii (16,771 bp) were circular molecules with gene organisations identical to that of the typical vertebrate mitogenome. The A + T content was higher than the G + C content, with a bias towards A and C nucleotides observed in all complete mitogenomes. The stem-and-loop secondary structures of the putative origin of light-strand replication were found to have highly conserved synthesis and stem regions, with all substitutions and indels restricted to the loop structure. The ratios of non-synonymous to synonymous substitution rates indicated that purifying selection has been the dominant driver of evolution in rhinopristiform mitogenomes. Phylogenetic reconstructions placed Acroteriobatus as a sister-group to Rhinobatos, confirming its affiliation with the family Rhinobatidae. However, based on its apparent polyphyly with the aforementioned genera, the familial assignment of Pseudobatos is not fully resolved and requires further investigation.

A multitude of spots! Five new microendemic species of the Cnemaspis gracilis group (Squamata: Gekkonidae) from massifs in the Shevaroy landscape, Tamil Nadu, India

South Asian Cnemaspis are one of the most diverse clades of gekkonids in South Asia with their highest diversity in the Western Ghats and Sri Lanka. These geckos include only a few nocturnal species and are largely diurnal or cathemeral and restricted to relatively cool habitats. One of the prominently diurnal subgroups in South Asian Cnemaspis is the bangara clade, which includes six species distributed in southern India on the eastern slopes of the Western Ghats, the southern Eastern Ghats and Palghat Gap. In this paper, we describe five more species of the bangara clade from the Shevaroyan landscape, including three from Kollimalai and one each from Yercaud and Pachaimalai, all in Tamil Nadu. These new species show 4.6–19.7 % uncorrected sequence divergence on the mitochondrial ND2 gene from each other and known species of the bangara clade and are morphologically diagnosable in body size, the number of paravertebral tubercles between limb insertions, the number of dorsal tubercle rows, the number of ventral scale rows across the belly, the number of femoral and precloacal pores and poreless scales separating these series, and aspects of colouration. The discovery of these five new species adds to the growing discoveries of cool-adapted species in southern India outside the Western Ghats and highlights the role of sky-islands in diversification. The Shevaroyan landscape shows high levels of microendemism with eight species distributed in an area of &lt; 2000 km2, and all these species restricted to much smaller areas of actual distribution. With an area of &lt; 500 km2 respectively, the massif of Pachaimalai has a single endemic and the massifs of Yercaud and Kollimalai have three endemic Cnemaspis species each.

Mitochondrial DNA Diversity of Mesocricetus auratus and Other Cricetinae Species among Cricetidae Family

Unique anatomical and physiological features have made hamster species desirable research models. Comparative genomics and phylogenetic analysis of the hamster family members to clarify their evolution and genetic relationship, can provide a genetic basis for the comprehension of the variable research results obtained using different hamster models. The Syrian golden hamster (Mesocricetus auratus) is the most widely used species. In this study, we sequenced the complete mitochondrial genome (mitogenome) of M. auratus, compared it with the mitogenome of other Cricetinae subfamily species, and defined its phylogenetic position in the Cricetidae family. Our results show that the mitogenome organization, gene arrangement, base composition, and genetic analysis of the protein coding genes (PCGs) of M. auratus are similar to those observed in previous reports on Cricetinae species. Nonetheless, our analysis clarifies some striking differences of M. auratus relative to other subfamily members, namely distinct codon usage frequency of TAT (Tyr), AAT (Asn), and GAA (Glu) and the presence of the conserved sequence block 3 (CSB-3) in the control region of M. auratus mitogenome and other hamsters (not found in Arvicolinae). These results suggest the particularity of amino acid codon usage bias of M. auratus and special regulatory signals for the heavy strand replication in Cricetinae. Additionally, Bayesian inference/maximum likelihood (BI/ML) tree shows that Cricetinae and Arvicolinae are sister taxa sharing a common ancestor, and Neotominae split prior to the split between Cricetinae and Arvicolinae. Our results support taxonomy revisions in Cricetulus kamensis and Cricetulus migratorius, and further revision is needed within the other two subfamilies. Among the hamster research models, Cricetulus griseus is the species with highest sequence similarity and closer genetic relationship with M. auratus. Our results show mitochondrial DNA diversity of M. auratus and other Cricetinae species and provide genetic basis for judgement of different hamster models, promoting the development and usage of hamsters with regional characteristics.

Complete Mitogenome of Oreolalax Omeimontis Reveals Phylogenetic Status and Novel Gene Arrangement of Archaeobatrachia

Species of the genus Oreolalax displayed crucial morphological characteristics of vertebrates transitioning from aquatic to terrestrial habitats; thus, they can be regarded as a representative vertebrate genus for this landing phenomenon. But the present phylogenetic status of Oreolalax omeimontis has been controversial with morphological and molecular approaches, and specific gene rearrangements were discovered in all six published Oreolalax mitogenomes, which are rarely observed in Archaeobatrachia. Therefore, this study determined the complete mitogenome of O. omeimontis with the aim of identifying its precise phylogenetic position and novel gene arrangement in Archaeobatrachia. Phylogenetic analysis with Bayesian inference and maximum likelihood indicates O. omeimontis is a sister group to O. lichuanensis, which is consistent with previous phylogenetic analysis based on morphological characteristics, but contrasts with other studies using multiple gene fragments. Moreover, although the duplication of trnM occurred in all seven Oreolalax species, the translocation of trnQ and trnM occurred differently in O. omeimontis to the other six, and this unique rearrangement would happen after the speciation of O. omeimontis. In general, this study sheds new light on the phylogenetic relationships and gene rearrangements of Archaeobatrachia.

A new species of feather-tailed leaf-toed gecko, Kolekanos Heinicke, Daza, Greenbaum, Jackman, Bauer, 2014 (Squamata, Gekkonidae) from the poorly explored savannah of western Angola

We here describe a new species of feather-tailed leaf-toed gecko, Kolekanos, from southern Benguela Province, Angola, based on morphological and osteological evidence, supported by phylogenetic analysis of mitochondrial data. The new species adds to the rapidly growing and newly-recognised endemic biodiversity of Angola, doubling the number of Kolekanos species, breaking the pattern observed within other closely-related African members of a clade of circum-Indian Ocean leaf-toed geckos - Ramigekko, Cryptactites and Afrogecko - all of which are presently monotypic. The new species is easily distinguished from K.plumicaudus, based on spine-like (as opposed to feather-like) scales on the margins of the original tail. Phylogenetic analyses also recovered the new taxon as monophyletic, with a well-supported sister relationship to K.plumicaudus, from which it differs by a substantial 24.1% NADH-dehydrogenase subunit 2 mitochondrial gene uncorrected p-distance.

Comparative mitogenomes provide new insights into phylogeny and taxonomy of the subfamily Xenocyprinae (Cypriniformes: Cyprinidae)

Xenocyprinae is a cyprinid subfamily that not only has a discrete geographic distribution but also has a long history dating to the Early Miocene. However, it is controversial whether systematic classification and some species validity of Xenocyprinae exist, as well as its phylogenetic relationships and evolutionary history. In the present study, we first reviewed the description and taxonomic history of Xenocyprinae, and then the complete mitochondrial genome of Distoechodon compressus, an endemic and locally distributed species belonging to Xenocyprinae, was sequenced and annotated. Finally, all the mitogenomes of Xenocyprinae were compared to reconstruct the phylogenetic relationship and estimate the divergence time. The results showed that the mitogenomes are similar in organization and structure with 16618–16630 bp length from 12 mitogenomes of eight species. Phylogenetic analysis confirmed the monology of Xenocyprinae and illustrated three clades within the Xenocyprinae to consist of ambiguous generic classification. Plagiognathops is a valid genus located at the base of the phylogenetic tree. The genus Xenocypris was originally monophyletic, but X. fangi was excluded. Divergence time estimation revealed that the earliest divergence within Xenocyprinae occurred approximately 12.1 Mya when Plagiognathops separated from the primitive Xenocypris. The main two clades (Xenocypris and (Distoechodon + Pseudobrama + X. fangi)) diverged 10.0 Mya. The major divergence of Xenocyprinae species possibly occurred in the Middle to Late Miocene and Late Pliocene, suggesting that speciation and diversifications could be attributed to the Asian monsoon climate. This study clarifies some controversial issues of systematics and provides essential information on the taxonomy and phylogeny of the subfamily Xenocyprinae.

Novel Gene Rearrangement in the Mitochondrial Genome of Three Garra and Insights Into the Phylogenetic Relationships of Labeoninae

Complete mitochondrial genomes (mitogenomes) can provide valuable information for phylogenetic relationships, gene rearrangement, and molecular evolution. Here, we report the mitochondrial whole genomes of three Garra species and explore the mechanisms of rearrangements that occur in their mitochondrial genomes. The lengths of the mitogenomes’ sequences of Garra dengba, Garra tibetana, and Garra yajiangensis were 16,876, 16,861, and 16,835, respectively. They contained 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNA genes, and two identical control regions (CRs). The mitochondrial genomes of three Garra species were rearranged compared to other fish mitochondrial genomes. The tRNA-Thr, tRNA-Pro and CR (T-P-CR) genes undergo replication followed by random loss of the tRNA-Thr and tRNA-Pro genes to form tRNA-Thr, CR1, tRNA-Pro and CR2 (T-CR-P-CR). Tandem duplication and random loss best explain this mitochondrial gene rearrangement. These results provide a foundation for future characterization of the mitochondrial gene arrangement of Labeoninae and further phylogenetic studies.

Positive Correlation of the Gene Rearrangements and Evolutionary Rates in the Mitochondrial Genomes of Thrips (Insecta: Thysanoptera)

Simple Summary

Aeolothrips, commonly known as banded thrips, is the largest genus of the family Aeolothripidae (predatory thrips). In the current study, we sequenced the mitochondrial genome (mitogenome) of the banded thrip species Aeolothrips xinjiangensis. We found a novel gene arrangement in this mitogenome that has not been reported in Thysanoptera. By comparing the gene order and rearrangement patterns, we found seven identical gene blocks and three identical rearrangement events in two mitogenomes of banded thrips. There was marked variation in the mitochondrial gene order across thrip species, with only two conserved gene blocks shared by all 14 thrips. In addition, we found a positive correlation between the degree of gene rearrangement and evolutionary rate. Our results suggested that the mitogenomes of thrips have tended to be stable since their massive rearrangement.

Abstract

Extensive gene rearrangement is characteristic in the mitogenomes of thrips (Thysanoptera), but the historical process giving rise to the contemporary gene rearrangement pattern remains unclear. To better understand the evolutionary processes of gene rearrangement in the mitogenomes of thrips, we sequenced the mitogenome of the banded thrip species Aeolothrips xinjiangensis. First, we found a novel mitochondrial gene order in this species. This mitogenome is 16,947 bp in length and encodes the typical 37 coding genes (13 protein-coding genes, 22 tRNA genes, and two rRNA genes) of insects. The gene arrangement was dramatically different from the putative ancestral mitogenome, with 26 genes being translocated, eight of which were inverted. Moreover, we found a novel, conserved gene block, trnC-trnY, which has not been previously reported in the mitogenomes of thrips. With this newly assembled mitogenome, we compared mitogenome sequences across Thysanoptera to assess the evolutionary processes giving rise to the current gene rearrangement pattern in thrips. Seven identical gene blocks were shared by two sequenced banded thrip mitogenomes, while the reversal of ND2 combined with TDRL events resulted in the different gene orders of these two species. In phylogenetic analysis, the monophyly of the suborders and families of Thysanoptera was well supported. Across the gene orders of 14 thrips, only two conserved gene blocks, ATP8-ATP6 and ND4-ND4L, could be found. Correlation analysis showed that the degree of gene rearrangement was positively correlated with the non-synonymous substitution rate in thrips. Our study suggests that the mitogenomes of thrips remain stable over long evolutionary timescales after massive rearrangement during early diversification.

Complete mitochondrial genome of Episymploce splendens (Blattodea: Ectobiidae): A large intergenic spacer and lacking of two tRNA genes

The complete mitochondrial genome of Episymploce splendens, 15,802 bp in length, was determined and annotated in this study. The mito-genome included 13 PCGs, 20 tRNAs and 2 rRNAs. Unlike most typical mito-genomes with conservative gene arrangement and exceptional economic organization, E. splendens mito-genome has two tRNAs (tRNA-Gln and tRNA-Met) absence and a long intergenic spacer sequence (93 bp) between tRNA-Val and srRNA, showing the diversified features of insect mito-genomes. This is the first report of the tRNAs deletion in blattarian mito-genomes and we supported the duplication/random loss model as the origin mechanism of the long intergenic spacer. Two Numts, Numt-1 (557 bp) and Numt-2 (975 bp) transferred to the nucleus at about 14.15 Ma to 22.34 Ma, and 19.19 Ma to 24.06 Ma respectively, were found in E. splendens. They can be used as molecular fossils in insect phylogenetic relationship inference. Our study provided useful data for further studies on the evolution of insect mito-genome.

Two new species of Hemidactylus Goldfuss, 1820 (Squamata: Gekkonidae) from the rocky outcrops of the Deccan plateau, peninsular India

The gekkonid diversity of the under-explored dry zones in Peninsular India is not yet completely inventoried, as evidenced by frequent new species descriptions. Here, we describe two new species of Hemidactylus from the dry rocky zones of Deccan Plateau viz. Hemidactylus mahonyi sp. nov. from Sandur Hills and Hemidactylus srikanthani sp. nov. from Devarayana Durga Hills, Karnataka State of Peninsular India. The new species belong to the morphologically cryptic H. murrayi Clade but differ from each other and the known congeners by a set of morphological characters and possess non-overlapping distribution ranges. Furthermore, the new species are also divergent (> 6.4% uncorrected p-distance) from their nearest related congeners, and from one another with respect to the ND2 mitochondrial gene. Recently, other new species of reptiles have been described from the two aforesaid type localities, indicating that further bio-inventories in the dry zones of peninsular India will yield more such novelties.  

How many more species are out there? Current taxonomy substantially underestimates the diversity of bent-toed geckos (Gekkonidae, Cyrtodactylus) in Laos and Vietnam

Cyrtodactylus is the most diverse genus of the family Gekkonidae and the world's third largest vertebrate genus. The number of species has increased more than fourfold over the last two decades. Indochina, especially Vietnam and Laos, has witnessed a surge in new species discoveries over the last three decades. The species number reported from Laos and Vietnam has remarkably increased from five in 1997 to 71 species in 2021. However, within the genus, several taxonomic issues have not yet been fully resolved. Based on recently collected samples from Laos and Vietnam, we conducted a comprehensive molecular review of Cyrtodactylus occurring in Laos and Vietnam. Our molecular analysis with support from morphological comparisons showed that C.thuongae is a junior synonym of C.dati and C.rufford is a junior synonym of C.lomyenensis. In total, 68 described species distributed in Laos and Vietnam are undisputed with strong support from both molecular and morphological evidence. On the other hand, the molecular analyses revealed that there are at least seven undescribed species in Vietnam and Laos, one in the C.angularis group, one in the C.chauquangensis, and five in the C.irregularis group. This number will likely increase significantly, as previous work suggested that the C.angularis and C.irregularis groups harbor three and six unnamed lineages, respectively. Based on survey gaps identified in our study, it is clear that additional new species will be discovered in poorly studied regions of central Vietnam and northern and southern Laos. As many species in the genus are facing high extinction risks, several undescribed populations might already be severely threatened by human activities in both countries. Therefore, urgent taxonomic research is needed before conservation assessments of newly discovered taxa can be undertaken to protect them from anthropogenic threats.

Phylogenetic position of Bopyroides hippolytes, with comments on the rearrangement of the mitochondrial genome in isopods (Isopoda: Epicaridea: Bopyridae)

Background

Classification of parasitic bopyrids has traditionally been based on morphological characteristics, but phylogenetic relationships have remained elusive due to limited information provided by morphological data and tendency for loss of morphological features as a result of parasitic lifestyle. Subfamily Argeiinae was separated from Bopyrinae based on morphological evidence, although the assignment of all genera has not been phylogenetically evaluated. Bopyroides hippolytes has been traditionally classified in Bopyrinae, but divergent morphological characters make this assignment questionable. To investigate the relationship of bopyrines, we sequenced the complete mitochondrial genome of B. hippolytes and four mitochondrial genes of two other Bopyrinae.

Results

The phylogenetic trees based on separate and combined cox1and 18S sequence data recovered Bopyridae as robustly monophyletic, but Bopyrinae as polyphyletic. Bopyroides hippolytes was a close sister to Argeia pugettensis, type species to Argeiinae. Mitochondrial phylogenomics also suggested that B. hippolytes was close to Argeiinae. We also found a novel gene order in B. hippolytes compared to other isopods.

Conclusions

Bopyroides hippolytes should be excluded from the Bopyrinae and has a close affinity with Argeia pugettensis based on molecular and morphological data. The conserved syntenic blocks of mitochondrial gene order have distinctive characteristics at a subordinal level and may be helpful for understanding the higher taxonomic level relationships of Isopoda.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08513-9.

Comparative genomic analysis of vertebrate mitochondrial reveals a differential of rearrangements rate between taxonomic class

Vertebrate mitochondrial genomes have been extensively studied for genetic and evolutionary purposes, these are normally believed to be extremely conserved, however, different cases of gene rearrangements have been reported. To verify the level of rearrangement and the mitogenome evolution, we performed a comparative genomic analysis of the 2831 vertebrate mitochondrial genomes representing 12 classes available in the NCBI database. Using a combination of bioinformatics methods, we determined there is a high number of errors in the annotation of mitochondrial genes, especially in tRNAs. We determined there is a large variation in the proportion of rearrangements per gene and per taxonomic class, with higher values observed in Actinopteri, Amphibia and Reptilia. We highlight that these are results for currently available vertebrate sequences, so an increase in sequence representativeness in some groups may alter the rearrangement rates, so in a few years it would be interesting to see if these rates are maintained or altered with the new mitogenome sequences. In addition, within each vertebrate class, different patterns in rearrangement proportion with distinct hotspots in the mitochondrial genome were found. We also determined that there are eleven convergence events in gene rearrangement, nine of which are new reports to the scientific community.

A new species of bent-toed geckos of the genus Cyrtodactylus Gray, 1827 from western Arunachal Pradesh, India

A new species of bent-toed geckos of the genus Cyrtodactylus is described from western Arunachal Pradesh. The new species is a member of the Indo-Burma clade and is embedded within the “peguensis” group, a relationship deduced contingent on a partial fragment of mitochondrial NADH subunit 2 gene. Cyrtodactylus kamengensissp. nov. is morphologically similar to C. himalayicus from which it differs in bearing a distinct ventrolateral fold and 6–8 basal lamellae on digit IV of pes. Genetically, the new species is sister to the Indian lineage of the “peguensis” group containing C. bhupathyi and C. gubernatoris. The Indian lineage of the “peguensis” group diverged from its Burmese relatives during the mid-Oligocene likely followed by the beginning of the Himalayan uplift, highlighting the role of the Himalayas in the diversification of biota.

Phylogenetic and multivariate analyses of Gekko smithii Gray, 1842 recover a new species from Peninsular Malaysia and support the resurrection of G. albomaculatus (Giebel, 1861) from Sumatra

Abstract

An integrative taxonomic analysis of Sundaic populations of Gekko smithii from the Thai-Malaya Peninsula, Sumatra, and Borneo recovered four deeply divergent mitochondrial lineages that are separated by major geographic barriers (mountains and seaways). Furthermore, they bear a number of concordant statistically significant differences in meristic and morphometric features, morphospatial separation in multivariate space, and discrete differences in color pattern. Gekko smithiisensu stricto is restricted to southern Thailand south of the Isthmus of Kra and Peninsular Malaysia west of the Banjaran (mountain range) Titiwangsa, being that the type locality is on Penang Island, Penang. Gekko hulksp. nov. is a new species from extreme southern Thailand and Peninsular Malaysia east of the Banjaran Titiwangsa and five east coast islands—the type locality being Pulau (island) Tioman, Pahang. Gekko cf. albofasciolatus is tentatively used to include Bornean populations west of the Iran Mountains in Sabah and Sarawak which, in the absence of molecular data, cannot unequivocally be separated morphologically from G. albofasciolatus from the type locality at Banjarmasin, Kalimantan, Indonesia east of the Iran Mountains. In the absence of molecular data, G. albomaculatus is resurrected to include mainland Sumatran, Nias Island, and Banyak Islands populations which, based on their morphology, cannot be separated from descriptions of G. albomaculatus from the type locality of Bangka Island, 15 km off the southeast coast of mainland Sumatra. Further integrative analyses of all Sumatran and Bornean populations are currently underway as well as the enigmatic Wallacean populations from Sulawesi. Data are presented that strongly suggest all references to G. smithii from Java stem from a 151 year-old misidentification of a specimen of G. gecko of unknown provenance. Additionally, there are no vouchered records of G. smithii from Myanmar. The phylogeographic patterns of Sundaic populations of the G. smithii complex are concordant with those of a plethora of other Sundaic lineages.

Comparative Analysis of Mitochondrial Genomes in Two Subspecies of the Sunwatcher Toad-Headed Agama (Phrynocephalus helioscopus): Prevalent Intraspecific Gene Rearrangements in Phrynocephalus

Intraspecific rearrangements of mitochondrial genomes are rarely reported in reptiles, even in vertebrates. The sunwatcher toad-headed agama, Phryncoephalus helioscopus, can serve as an excellent model for investigating the dynamic mitogenome structure at intraspecific level. To date, seven subspecies of P. helioscopus are well recognized, but little is known about the mitogenomic evolution among different subspecies. In this study, complete mitogenomes of subspecies P. helioscopus varius II and P. helioscopus cameranoi were determined by next-generation sequencing, and another P. helioscopus varius I retrieved from GenBank was compiled for comparative analysis. The nucleotide composition and the codon usage are similar to those previously published from toad-headed agamas. P. helioscopus varius II and P. helioscopus cameranoi have 23 tRNA genes, including standard 22 tRNA genes and one extra tRNA-Phe (tRNA-Phe duplication). Gene order and phylogenetic analyses in the genus Phrynocephalus support prevalent intraspecific gene rearrangement in P. helioscopus and other congener species including P. erythrurus, P. vlangalii, and P. forsythii. Six different mitochondrial gene arrangements are observed in Phrynocephalus. Overall, the occurrence of rearrangements may result from multiple independent structural dynamic events. The split of the two subspecies in P. helioscopus was dated at approximately 2.34 million years ago (Ma). Two types of gene rearrangements are found in the three mitogenomes of P. helioscopus, and this intraspecific rearrangement phenomenon can be explained by the tandem duplication/random loss (TDRL) model. Post duplication, the alternative loss types can occur in 0.23–0.72 Ma, suggesting that the duplication and fixation of these rearrangements can occur quite quickly. These findings highlight the need for more mitogenomes at the population level in order to better understand the potentially rampant intraspecific mitogenomic reorganization in Phrynocephalus.

The evolutionary history of an accidental model organism, the leopard gecko Eublepharis macularius (Squamata: Eublepharidae)

The leopard gecko, Eublepharis macularius, is a widely used model organism in laboratory and experimental studies. The high phenotypic diversity in the pet trade, the fact that the provenance of different breeding lines is unknown, and that distinct Eublepharis species are known to hybridize, implies that the continued use of E. macularius as a model requires clarity on the origin of the lineages in the pet trade. We combine multi-locus sequence data and the first range-wide sampling of the genus Eublepharis to reconstruct the evolutionary history of the Eublepharidae and Eublepharis, with an updated time-tree for the Eublepharidae. Our sampling includes five of the six recognized species and additional nominal taxa of uncertain status comprising 43 samples from 34 localities plus 48 pet-trade samples. The Eublepharidae began diversifying in the Cretaceous. Eublepharis split from its sister genera in Africa in the Palaeocene-Eocene, and began diversifying in the Oligocene-Miocene, with late Miocene-Pliocene cladogenesis giving rise to extant species. The current species diversity within this group is moderately underestimated. Our species delimitation suggests 10 species with four potentially unnamed divergent lineages in Iran, India and Pakistan. All 30 individuals of E. macularius that we sampled from the pet trade, which include diverse morphotypes, come from a few shallow E. macularius clades, confirming that lab and pet trade strains are part of a single taxon. One of the wild-caught haplotypes of E. macularius, from near Karachi, Pakistan, is identical to (10) pet-trade samples and all other captive populations are closely related to wild-caught animals from central/southern Pakistan (0.1-0.5 % minimum pairwise uncorrected ND2 sequence divergence).

A new ancient lineage of ablepharine skinks (Sauria: Scincidae) from eastern Himalayas with notes on origin and systematics of the group

The Himalayas represent a renowned biodiversity hotspot and an important biogeographic realm that has influenced origin and diversification of multiple taxa. A recent herpetological investigation of the eastern Himalayas of the Indian state of Arunachal Pradesh led to the discovery of a unique lineage of ablepharine skink, which is herein described as a new genus along with a new species. The findings are based an integrated taxonomic approach incorporating data from external morphology, microCT scans of the skull and molecular data. The molecular phylogeny of ablepharine skinks is also presented that suggests taxonomic amendments. Discovery of this unique lineage of skinks further highlights the biogeographic importance of the eastern Himalayas as a source for origin of several relic biota.