Next Article in Journal
High Seroprevalence against SARS-CoV-2 among Dogs and Cats, Poland, 2021/2022
Previous Article in Journal
Mechanical Testing of the New Cage for Tibial Tuberosity Advancement with the Cranial Implant Fixation (TTA CF) Technique—Ex Vivo Study on Sheep Model
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 12
Issue 16
10.3390/ani12162015
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Complete Mitochondrial Genomes of Three Sphenomorphinae Species (Squamata: Scincidae) and the Selective Pressure Analysis on Mitochondrial Genomes of Limbless Isopachys gyldenstolpei

by
Lian Wu
1,
Yao Tong
1,
Sam Pedro Galilee Ayivi
1,
Kenneth B. Storey
2,
Jia-Yong Zhang
1,2 and
Dan-Na Yu
1,3,*
1
College of Chemistry and Life Science, Zhejiang Normal University, Jinhua 321004, China
2
Department of Biology, Carleton University, Ottawa, ON K1S5B6, Canada
3
Key Lab of Wildlife Biotechnology, Conservation and Utilization of Zhejiang Province, Zhejiang Normal University, Jinhua 321004, China
*
Author to whom correspondence should be addressed.
Animals 2022, 12(16), 2015; https://doi.org/10.3390/ani12162015
Submission received: 1 July 2022 / Revised: 27 July 2022 / Accepted: 7 August 2022 / Published: 9 August 2022
(This article belongs to the Section Animal Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Skinks are the most species-rich group of lizards and are widely distributed around the world. The family Scincidae (Reptilia: Lacertiformes) includes limbed and limbless representatives occupying diverse habitats and showing a range of morphologies. Both limbed and limbless skinks have unique locomotion patterns in their habitats. Locomotion is the process of energy consumption, of which different modes may have different energy demands. As the center of energy metabolism in organisms, mitochondria provide most of the energy for physiological and biochemical activities via oxidative phosphorylation. Here, we employed mitochondrial genomes to investigate potential selective pressures among limbless skinks. Isopachys gyldenstolpei, as a typical limbless skink, has a different locomotion pattern compared to a limbed skink. Thus, I. gyldenstolpei can be used to study whether limb loss has a positive selection on mitochondrial genes. Two typical limbed skinks, Sphenomorphus indicus and Tropidophorus hainanus, were included in this study to compare the selective pressure analysis on mitochondrial genomes. In addition, the phylogenetic relationships within Scincidae are also discussed.

Abstract

In order to adapt to diverse habitats, organisms often evolve corresponding adaptive mechanisms to cope with their survival needs. The species-rich family of Scincidae contains both limbed and limbless species, which differ fundamentally in their locomotor demands, such as relying on the movement of limbs or only body swing to move. Locomotion requires energy, and different types of locomotion have their own energy requirements. Mitochondria are the energy factories of living things, which provide a lot of energy for various physiological activities of organisms. Therefore, mitochondrial genomes could be tools to explore whether the limb loss of skinks are selected by adaptive evolution. Isopachys gyldenstolpei is a typical limbless skink. Here, we report the complete mitochondrial genomes of I. gyldenstolpei, Sphenomorphus indicus, and Tropidophorus hainanus. The latter two species were included as limbed comparator species to the limbless I. gyldenstolpei. The results showed that the full lengths of the mitochondrial genomes of I. gyldenstolpei, S. indicus, and T. hainanus were 17,210, 16,944, and 17,001 bp, respectively. Three mitochondrial genomes have typical circular double-stranded structures similar to other reptiles, including 13 protein-coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and the control region. Three mitochondrial genomes obtained in this study were combined with fifteen mitochondrially complete genomes of Scincidae in the NCBI database; the phylogenetic relationship between limbless I. gyldenstolpei and limbed skinks (S. indicus and T. hainanus) is discussed. Through BI and ML trees, Sphenomorphinae and Mabuyinae were monophyletic, while the paraphyly of Scincinae was also recovered. The limbless skink I. gyldenstolpei is closer to the species of Tropidophorus, which has formed a sister group with (T. hainanus + T. hangman). In the mitochondrial genome adaptations between limbless I. gyldenstolpei and limbed skinks, one positively selected site was found in the branch-site model analysis, which was located in ND2 (at position 28, BEB value = 0.907). Through analyzing the protein structure and function of the selected site, we found it was distributed in mitochondrial protein complex I. Positive selection of some mitochondrial genes in limbless skinks may be related to the requirement of energy to fit in their locomotion. Further research is still needed to confirm this conclusion though.

1. Introduction

Scincidae, as the most species-rich family in Squamata, has 1744 species currently recorded on the Reptile Database (http://www.reptile-database.org/, accessed on 22 March 2022) [1]. All live in tropical and temperate regions and can be found on all continents, except for Antarctica and most oceanic islands [2]. With such a wide distribution, it is not surprising that skinks exhibit diverse morphological characteristics to adapt to their various distribution environments [3,4]. Skinks range from fully pentadactyl forms to completely limbless body forms found in Scincidae [1]. Compared with the wide habitat range of limbed skinks, most limbless skinks are cryptozoic, mainly choosing to inhabit underground burrows or under rocks, with little exposure above ground [5,6]. At the same time, different movement abilities between limbed and limbless skinks have arisen in adapting to their habitats. Limbed body forms are often associated with high dispersal ability in Scincidae, whereas limbless body forms show reduced patterns of dispersal. The secretive living environments of the skink species and the scarcity of available samples have led to a poor understanding of the diversity, geographical distribution, genetic relationships, and phylogeny of these groups [7].
Isopachys gyldenstolpei belongs to the subfamily Sphenomorphinae and is mainly distributed in Thailand and Burma [8]. As an elongated, fossorial, and limbless lizard, I. gyldenstolpei burrows in dry soil covered with dead leaves and sticks [8,9]. Due to a scarcity of samples, no complete mitochondrial genome of I. gyldenstolpei has been reported to date. Compared with limbless skinks, the limbed skinks are not only more widely distributed, but also have diverse lifestyles. Thus, there is much interest in the typical morphological characteristics of limb loss and axial elongation in Scincidae. For example, the species-rich genus, Sphenomorphus, shows diverse lifestyles in a wide array of habitats, ranging from upland cloud forests to small, arid, virtually barren islands, to lowland and hill forests dominated by Dipterocarpus [10,11,12]. In contrast, lizards of the genus Tropidophorus are semi-aquatic and mainly dwell in lowlands near mainland forest streams [13,14,15]. The present study sequenced the mitochondrial genome of one limbless skink, I. gyldenstolpei, followed by the mitochondrial genomes of two limbed skinks Sphenomorphus indicus and Tropidophorus hainanus. Differences in the mitochondrial genome were compared to the morphological disparity between limbed and limbless skinks. At present, the monophyly of Scincidae has been widely recognized, and the taxonomic controversy is mainly focused on the monophyly of the subfamily Scincinae [16,17,18,19]. As an ideal molecular marker, mitochondrial genes have been widely used to solve the phylogeny of Scincidae [20]. Mitochondrial genomes were used to discuss the phylogenetic status of limbless I. gyldenstolpei and mitogenomes could similarly provide useful clues for the remaining disputes in the phylogenetic relationships of Scincinae.
Moving is an energy-expensive process, and different models of movement require different levels of energy consumption. Losing legs or damage to appendages associated with movement may impose additional energetic costs [21]. As the energy-producing center of organisms, up to 95% of the ATP used by eukaryotic cells is provided by mitochondria through oxidative phosphorylation (OXPHOS) [22,23,24,25]. The mitochondrial genomes of lizards have the same double-circular structures as those of other vertebrates, with typical lengths of 16–19 kb and conserved genetic compositions, including 13 protein-coding genes, 22 tRNAs, 2 rRNAs, and a major non-coding control region [22,26]. Due to its function in energy metabolism, mitochondrial DNA has attracted much attention for its selective role in adaptive evolution, and some studies have shown that the differences in locomotive ability are reflected in the different selection mechanisms of mitochondrial genomes [24,25,27,28,29,30]. For example, compared with flightless grasshoppers, flying grasshoppers have been positively selected for their mitochondrial genes in response to the energy requirements of flight [31]. Among fish, it has been found that the amount of swimming and energy consumption are different between migratory fishes and nonmigratory fishes, and genes related to the energy regulation in mitochondria have shown evidence for a direct response to selective pressure [24]. In birds, the mitochondrial genes of fast-flying birds have different effective energy metabolic requirements than those of weak-flying or flightless species, suggesting that evolution constrains the mitochondrial DNA of birds with different motor abilities, as has also been found in mammals [25]. Similarly, the mitochondrial genes of plateau birds were positively selected after adaptation to the plateau environment [32]. Research on the origins of flight in bats has shown that mitochondrial genes involved in energy metabolism were naturally selected to adapt to the huge changes in energy requirements during the origins of flight (as compared to other mammals) [28]. Selection pressures on mitochondrial genes have also been found in limb-reduced Dibamus bourreti and other limbless squamates, suggesting that selection pressures on mitochondrial genomes may play important roles in the energetic differences in locomotion between limbed and limbless squamates [33]. Therefore, considering that I. gyldenstolpei has morphological characteristics related to a complete loss of limbs, we speculate that the mitochondrial genes of limbless skinks may undergo positive selections compared to limbed skinks.

2. Material and Methods

2.1. Specimen Collection

The soak specimen of I. gyldenstolpei was housed in the Animal Herbarium of Zhejiang Normal University, which was collected from Myanmar in 2000. Specimens of S. indicus and T. hainanus were collected in Cixi, Zhejiang and Jinxiu, Guangxi, China, respectively. Our experimental procedures complied with current regulations on animal welfare and research in China. The Animal Research Ethics Committee of Zhejiang Normal University approved the experimental design (ZSDW2021057).

2.2. DNA Extraction, PCR Amplification, and Sequencing

Total genomic DNA of the three skinks was extracted from tail muscles using an Ezup Column Animal Genomic DNA Purification Kit (Sangon Biotech Company, Shanghai, China). Common primers [34] for lizards were used to amplify several partial segments, and species-specific primers were designed to complete the remaining gaps using Primer Premier 5.0 [35] based on previously obtained sequences (Table S1). The methods described in Zhang et al. [36] were used to amplify both PCR (product length < 3000 bp) and long-range-PCR (product length > 3000 bp). All PCR products were sequenced using the bi-directional primer-walking method by Sangon Biotech Company (Shanghai, China).

2.3. Mitochondrial Genome Annotation and Sequence Analyses

Sequences were checked and assembled using DNASTAR Package V.7.1 [37]. The tRNA genes of the sequences were annotated using the MITOS web server (http://mitos.bioinf.uni-leipzig.de/index.py, accessed on 15 December 2021) [38], and their secondary structures were predicted by the tRNAscan-SE online search server [39]. We identified the two rRNA genes (12S and 16S rRNA) and 13 protein-coding genes (PCGs) by comparing them with the complete mitochondrial genomes of other Scincidae in the Clustal W program of Mega 7.0 [40]. The CG View online server V 1.0 [41] was used to draw the maps of complete mitochondrial genomes (https://cgview.ca/, accessed on 15 January 2022). According to the formula: AT skew = (A − T) ÷ (A + T), GC skew = (G − C) ÷ (G + C), we calculated the CG and AT skews in PhyloSuite 1.2.2 [42]. Similarly, PhyloSuite 1.2.2 was also used to calculate the AT content, codon usage, and relative synonymous codon usage (RSCU) of protein-coding genes, and graphically drawn by Adobe Illustrator CS4 [43].

2.4. Phylogenetic Analyses

To explore the phylogenetic relationship of Scincidae, we prepared a dataset comprised of the complete mitochondrial genomes of the three skink species in this study and 15 skink mitochondrial genomes available in NCBI, including from Sphenomorphinae (9), Scincinae (5), and Mabuyinae (1) of Scincidae [44,45,46,47,48,49,50,51,52,53,54,55,56], as well as the two outgroups species belonging to Xantusiidae (1) and Cordyloidea (1) [57,58] (Table 1).
Table
Table 1. Information about the samples used in this study and the NCBI GenBank accession numbers.
Table 1. Information about the samples used in this study and the NCBI GenBank accession numbers.
FamilySubfamilySpeciesLengthAccession No.Reference
ScincidaeMabuyinaeEutropis multifasciata17,062 bpMN938934[44]
ScincinaeAteuchosaurus chinensis16,840 bpMW327509[45]
Plestiodon chinensis17,175 bpKT279358[46]
Plestiodon elegans17,304 bpKJ643142[47]
Plestiodon liui17,643 bpMT662111[48]
Plestiodon tunganus17,263 bpMK370739[49]
SphenomorphinaeAblepharus himalayanus17,304 bpMN885892[59]
Isopachys gyldenstolpei17,210 bpMH020638This study
Scincella huanrenensis17,212 bpKU507306[50]
Scincella modesta17,466 bpMN702771[51]
Scincella modesta17,511 bpMN786972[52]
Scincella reevesii15,424 bpMN832615[53]
Scincella vandenburghi17,103 bpKU646826[50]
Sphenomorphus indicus16,944 bpOM117611This study
Sphenomorphus indicus17,027 bpMK450438[55]
Sphenomorphus incognitus17,417 bpMH329292[54]
Tropidophorus hainanus17,001 bpOM117612This study
Tropidophorus hangnam16,777 bpMN977920[56]
CordyloideaCordylidaeSmaug warreni17,184 bpNC005962[57]
XantusiidaeXantusiinaeLepidophyma flavimaculatum16,158 bpNC008775[58]
Based on a dataset of the nucleotide sequences of 13 (PCGs), we constructed the Bayesian inference (BI) and maximum likelihood (ML) phylogenetic trees of Scincidae. MAFFT V. 7.475 [60] was used for the alignment of the nucleotide sequences of the 13 PCGs. Then the conservative region was detected by Gblock 0.91b [61] using the default setting. Because the third codon was saturated, our phylogenetic trees were constructed with the 1st and 2nd codons; the 1st and 2nd data were extracted and reserved using MEGA7.0 [40]. Based on the Bayesian information criterion (BIC), the program PartitionFinder 2.2.1 [62] was used to determine the best partitioning scheme and substitution model (Table S2). Considering that RAxML allows for only a single model of rate heterogeneity in partitioned analyses, the ML analysis was performed in RAxML 8.2.0 [63] using the GTR + I + G model, with branches of each node supporting the evaluation under 1000 ultrafast replications. The BI analysis was performed in MrBayes version 3.2 [64] using the partitioning results (Table S2) and was set for 10 million generations with sampling every 1000 generations. The first 25% of generations were discarded as burn-ins. When the phylogenetic tree converged (i.e., when the value of the average standard deviation of the split frequency was stable and the balance was less than 0.01), the tree was harvested.

2.5. Detecting Selective Pressure

The selective pressure analysis of genes was mainly obtained by calculating the replacement ratio of nonsynonymous codons and synonymous codons in the protein-coding sequence, namely, ω = dN/dS, where ω = 1 represents neutral mutation, ω < 1 means negative selection, and ω > 1 indicates positive selection [65]. As an interactive visualization tool, EasyCodeML can be used to detect the selection pressure of 13 PCGs in the Scincidae mitochondrial genomes [66]. In order to investigate whether the mitochondrial genes of I. gyldenstolpei were under positive selection, this limbless skink was chosen as the foreground branch and other limbed skinks as the background branch. Selection pressure analysis was conducted using 11,373 nucleotide sites (excluding the initiation and termination codons) and 3791 amino acid sites after alignment of the 13 PCGs from 18 species of Scincidae. Four analytical methods were conducted to explore the mitochondrial gene adaptative evolution of limbless skinks. The branch model mainly defined the heterogeneity of ω values of different lineages in the phylogenetic tree, and we observed whether the ω values of foreground branches and background branches were significantly different by comparing the one-ratio model and two-ratio model. The site model was used to compare the selection pressure of amino acid sites. The ω values of the branch-site model were different between the selected sites, and it was also assumed that the ω values of the branches were different, which was mainly used to detect the influence of positive selection on some sites in the foreground branch (model A vs. model A null). The clade model was mainly used to analyze multiple clades, simultaneously [65]. The posterior probability of these models and selected loci were evaluated using the likelihood ratio test (LRTs) and Bayesian empirical Bayes (BEB), respectively. The structural and functional information of selected genes were obtained in UniProt [67], and the three-dimensional structures of the corresponding proteins were constructed by the SWISS-MODEL workspace [68].

3. Results

3.1. Organization and Characteristics of Three Skink Mitochondrial Genomes

The lengths of the complete mitochondrial genomes were found to be I. gyldenstolpei 17,210 bp, S. indicus 16,944 bp, and T. hainanus 17,001 bp, and their corresponding GenBank accession numbers are MH020638, OM117611, and OM117612, respectively. With the typical circular structure, the mitochondrial genome arrangements of the three skink species are the same as those of other reptiles and include 13 PCGs (ND1-6, ND4L, COI-III, Cyt b, ATP6, and ATP8), 22 tRNAs, and two rRNA genes (12S rRNA and 16S rRNA), along with a non-coding region (D-loop) between tRNA-Pro and tRNA-Phe (Figure 1). Among the 37 genes, ND6 and 8 tRNAs (including tRNA-Gln, tRNA-Ala, tRNA-Asn, tRNA-Cys, tRNA-Tyr, tRNA-Ser, tRNA-Glu, and tRNA-Pro) are coded on the minor strand (N strand) and the remaining 28 genes are coded on the major strand (J strand). A + T bias of the whole mitochondrial genome was found in all three skinks, being 56% in I. gyldenstolpei, 57.9% in S. indicus, and 54.2% in T. hainanus. Positive AT skews and negative GC skews were found in all three mitochondrial genomes (Table 2).
Gene spacers and overlaps were found between adjacent genes in all three mitochondrial genomes. Nine gene-spacing regions were found in I. gyldenstolpei and S. indicus, and eight were found in T. hainanus, all with lengths from 1 to 14 bp. The longest length of the gene-spacing regions (14 bp) was between tRNA-Asn and tRNA-Cys in all three mitochondrial genomes. Eight, nine, and eleven base pair overlaps were found in the mitochondrial genomes of I. gyldenstolpei, S. indicus, and T. hainanus, respectively, with the lengths ranging from 1 to 10 bp. The longest overlap (10 bp) was between ATP6 and ATP8 in all three skinks.
Among the three complete mitochondrial genomes, the total lengths of the 13 PCGs were 11,382 bp in I. gyldenstolpei and S. indicus, and 11,379 bp in T. hainanus. After comparing the three sequences, two start codons (ATG and GTG) were used in the 13 PCGs. As the most common start codon, ATG was used in most protein-coding genes (11 genes in I. gyldenstolpei and T. hainanus, 12 genes in S. indicus), whereas GTG was used only in COI (I. gyldenstolpei, S. indicus, and T. hainanus), ND1 (T. hainanus), and ND3 (I. gyldenstolpei). TAN (TAA/TAG) was often found in the three sequences as the stop codon (eight protein-coding genes). In addition, AGG/AGA was used as the stop codon in some protein-coding genes: COI and ND6 in all three mitochondrial genomes, and COII only in I. gyldenstolpei. An incomplete stop codon T was found in three protein-coding genes: COIII and ND4 in all three mitochondrial genomes, and COII of S. indicus and T. hainanus (Tables S3–S5).
Figure 2 and Supplementary Table S6 show the relative synonymous codon usage (RSCU) and amino acid composition of the three mitochondrial genomes. Throughout the data, A and C nucleotides were the most abundant in the base composition of the three skinks and especially appeared more frequently at the third codon position of all three skinks. The utilization rate of the CUA (Leu1) codon was the highest in the three skinks, which was used 213 times in I. gyldenstolpei, 231 times in S. indicus, and 249 times in T. hainanus. In addition, the use ratio of AUC (Ile), AUA (Met), ACC (Thr), ACA (Thr), and GCC (Ala) was relatively high in all three skinks. At the same time, the utilization rates of some amino acids for the synonymous codon in these three skinks were different, such as UCA (Ser2) in S. indicus 2.45 > T. hainanus 2.04 > I. gyldenstolpei 1.95. Furthermore, A + T bias also existed in 13 PCGs of the three mitochondrial genomes (Table 2).
The two ribosomal RNAs (12S rRNA and 16S rRNA) were separated by tRNA-Val and were found on the H-strand. The full sizes of 12S rRNA were 949 bp (I. gyldenstolpei), 942 bp (S. indicus), and 958 bp (T. hainanus), respectively. Similarly, the total lengths of 16S rRNA for the three skinks were 1558 bp (I. gyldenstolpei), 1528 bp (S. indicus), and 1527 bp (T. hainanus). The sizes of all tRNA genes in the three skinks ranged from 64 to 76 bp. Compared with other tRNAs of the three sequences that showed typical cloverleaf structures, tRNA-Ser1 missed the dihydrouridine (DHU) arm (Figures S1–S3). The total lengths of tRNA in I. gyldenstolpei, S. indicus, and T. hainanus were between 1523 and 1542 bp.
With lengths of 1774 bp (I. gyldenstolpei), 1566 bp (S. indicus), and 1615 bp (T. hainanus), respectively, all control regions of the three skinks could be found between tRNA-Pro and tRNA-Phe. The A + T content of all three skink mitogenomes was more abundant, all greater than 60%. Positive AT skews and negative CG skews were also found in the control regions of the three skinks (Table 2).

3.2. Phylogenetic Relationships

Both BI and ML analyses recovered a consistent topology (Figure 3). In general, except for the outgroups, skinks in this study clustered mainly in three subfamilies: Sphenomorphinae, Scincinae, and Mabuyinae. Sphenomorphinae and Mabuyinae were recovered as monophyletic groups in both BI and ML trees. A sister-group relationship between Mabuyinae and the clade of (Sphenomorphinae + Scincinae) was supported in the present study. In ML and BI trees, the subfamily Scincinae is paraphyletic because Ateuchosaurus of Scincinae is nested into Sphenomorphinae. In the subfamily Sphenomorphinae, the limbless I. gyldenstolpei was a sister clade of (T. hainanus + T. hangnam) and then clustered in a clade with (Sp. indicus + Sp. incognitus) + (Sc. reevesii + (Sc. huanrenensis + (Sc. vandenburghi + Sc. Modesta)))). In ML and BI trees, Ab. himalayanus is considered the sister group to all other Sphenomorphinae species.

3.3. Positive Selection Analysis of 13 Protein-Coding Genes

The BI tree was used to analyze the selection pressure of 13 protein-coding genes. The results showed that LRT p-values were 0.859 and 0.704 in the branch and clade models, respectively. Thus, there were no significant selected sites found in the two models.
The analysis results of the branch-site model and the site model are shown in Table S7. Model A vs. model A null was significant (p < 0.05) in the branch-site model, where there was an amino acid selected site and the BEB value was 0.907 (amino acid residue 2017). Amino acid residue 2017 corresponds to amino acid position 28 of ND2 (Figure S4 and Table S8). At position 2017, amino acids were mutated from Leu (L) of the background branch to Thr (T) of the foreground branch. The positive selection site was found in mitochondrial complex I. Position 28 of the ND2 gene was found to be located in the domain of the proton-conducting membrane transporter through protein functional analysis. Two amino acid sites were positively selected (BEB > 0.90, positions 2416 and 3742) in the site model, with the LRT of M7-M8 showing high significance (p < 0.01). Amino acid residue 2416 corresponded to amino acid position 83 of the ND3 gene, whereas position 3742 corresponded to amino acid position 124 of the ND6 gene (Table S8). Positive selection sites of mitochondrial genes were found in the site model, which may suggest that some amino acid sites in some skinks have undergone adaptive selection.

4. Discussion

4.1. Phylogenetic Analyses

The monophyly of Scincidae was successfully recovered by both BI and ML trees, which is similar to previous studies [16,17,44,69,70]. The clade ((Sphenomorphinae + Scincinae) + Mabuyinae) of Scincidae was supported in both BI and ML trees, which is consistent with the results of Chen et al. (2021). However, this slightly contradicts the results from other studies that support the clade ((Mabuyinae + Scincinae) + Sphenomorphinae) [44,71]. Multiple studies suggest that Scincinae should be paraphyletic because Ateuchosaurus is nested into Lygosominae [16,18,70]. However, in the present study, although the paraphyly of Scincinae was recovered, the species of Ateuchosaurus was nested into Sphenomorphinae. This controversy may be related to the fact that Lygosominae was not included in our study. On the contrary, Scincinae was monophyletic and Ateuchosaurus was the sister to Lygosominae in a study using nuclear and mitochondrial genes by Pyron et al. [17]. When the genus Ateuchosaurus was not included in the phylogenetic tree, the monophyly of Scincinae was also restored [19,44,69]. In the future, samples of Ateuchosaurus and other skinks should be added to further explore the monophyly of Scincinae. In Mabuyinae, the phylogenetic trees constructed in this study and the results of Chen et al. (2021) based on 13 PCGs, all support that Mabuyinae was the sister of the clade (Scincinae + Sphenomorphinae). However, combined with morphological and molecular data, Mabuyinae and Sphenomorphinae showed the closest relationship [72]. Hence, the status of the Mabuyinae subfamily remains unclear due to the lack of mitochondrial genome data [44].
The species Isopachys clustered into a clade with the Tropidophorus genus and not with the Sphenomorphus genus, a result that is different from some previous conclusions [12,44,69,73]. Using 12S and 16S RNAs, the clades ((Isopachys + (Sphenomorphus + (Lipinia + Scincella))) + Tropidophorus) were supported in BI, ML, and maximum parsimony (MP) trees constructed by Honda et al. [73]. Although the same genes were used (12S and 16S RNA) in the phylogenetic analysis of Peninsular Malaysia skinks, Rizal et al. found that the clade of (((Tropidophorus + Sphenomorphus) + Lipinia) + Isopachys) was supported in both ML and BI trees [12]. These differences in topological structures may be related to the amount of data used to construct the phylogenetic tree. Sampling differences affecting the phylogenetic tree structure have also been found in other research [16,74,75]. The clade ((Sphenomorphus + Scincella) + Isopachys) was supported by using only the 13 PCGs in Chen et al. [69], whereas the tree that was constructed from 13 PCGs and 2 rRNAs (12S rRNA and 16S rRNA) showed ((Isopachys + Sphenomorphus) + Scincella) in Chen et al. [44]. At the same time, Tropidophorus was not included in the two references, and different genes used to construct the phylogenetic tree may also lead to inconsistent topological structure [16,44,69]. Thus, different samples and analytical methods will affect the results of the phylogenetic analysis. Hence, the relationship between Isopachys, Sphenomorphus, and Tropidophorus needs to be further verified.

4.2. Positive Selection Analyses

In general, positive selection acts on only a few amino acid sites over a short evolutionary period, so a signal of positive selection is usually swamped by subsequent successive negative selection at most sites in a gene sequence [76]. Similarly, after a short period of positive selection, there is usually a long purifying selection period, which neutralizes the results of the positive selection [28]. Therefore, positive selection in mitochondrial DNA is difficult to detect [28], especially in branch and clade models where ω varies only among branches and clades. On the contrary, the branch-site model and the site model are often used to distinguish the positive selection from the relaxed purification selection because they allow variations in selective pressure to occur at amino acid sites [28,76].
Consequently, analyses based on the branch-site model, suggesting a selective process that leads to the loss of limbs in I. gyldenstolpei was also recognizable in the ND2 gene of the mitochondrial complex I (NADH: ubiquinone oxidoreductase), which may indicate an adaptation of mitochondrial genes to the energy requirements of a limbless skink after limb loss occurs. Complex I is the first large protein complex in the respiratory chain, providing protonic power for ATP synthesis during the transfer of electrons to ubiquitin by NADH via transmembrane proton pumps. Indeed, the ND complex can drive over one-third of total energy (ATP) production in mitochondria, making it crucial in biological energy metabolism [77,78,79]. The smallest subsection of complex I is composed of the ND1-ND6 subunits that mainly form the transmembrane region core of the complex [80]. The main candidate genes for the proton pump in mitochondrial complex I are the ND2, ND4, and ND5 genes [81]. Mutations of these subunits may interfere with the efficiency of the proton pumping process [82]. In general, amino acid changes can cause inefficiencies in the electron transport chain system, which in turn affects mitochondrial energy production [82]. Certain research studies show a relationship between locomotion and energy requirements by using mitochondrial DNA. For example, it has been found that migratory fishes consume different amounts of energy compared with nonmigratory species due to differences in swimming that are reflected in the positive selections of ND1, ND3, and ND5 genes in the energy consumption of migratory fishes [24,83]. In squamates, comparisons can be made between different movement modes of limbed and limbless species. Locomotion was supported by the remaining limbs and body swing together in Dibamus bourreti, in which the ATP6 gene was selected compared with limbed species [33]. Furthermore, one of the major energy regulation genes in mitochondria of the limbless skink I. gyldenstolpei was selected (ND2 gene) during limb loss as compared with limbed skinks. This may be a relevant adaptation to the energetic requirements of I. gyldenstolpei. The results suggest that the difference in energy demand between limbless and limbed skinks may be related to the change of locomotion that promotes the evolution of mitochondrial DNA in different directions. However, considering the limited samples, this conclusion needs to be further verified. In future studies, we will consider adding more lizard samples to further explore whether the loss of limbs is associated with the adaptive evolution of energy-regulating genes.

5. Conclusions

In this study, the complete mitochondrial genome of a limbless skink, I. gyldenstolpei, was reported for the first time, representing the first complete mitochondrial genome of a limbless skink that has been recorded in Scincidae. At the same time, to compare with the mitochondrial genomes of limbless skinks, the complete mitochondrial genomes of two limbed skinks, S. indicus and T. hainanus, were also sequenced. Before this study, there was no complete mitochondrial genome of T. hainanus on NCBI, and only one complete mitochondrial genome of S. indicus. This study provides the complete mitochondrial genomes of S. indicus and T. hainanus to enrich the gene database of skinks. The mitochondrial genome structures of the three skinks were similar to the typical mitogenomes of reptiles. The clade of ((Sphenomorphinae + Scincinae) + Mabuyinae) in Scincidae was supported by both BI and ML analyses. The subfamily Scincinae was paraphyletic with respect to Sphenomorphinae and Mabuyinae. Our results support the conclusion that limbless Isopachys were more closely related to Tropidophorus, which as a result formed the clade (((Isopachys + Tropidophorus) + (Sphenomorphus + Scincella)) + Ablepharus). Compared with limbed skinks, one selective site was found in the mitochondrial genes of limbless I. gyldenstolpei. The selected site was located in the ND2 gene (position 2017, BEB = 0.907), and site 28 of the ND2 gene was located in the domain of the proton-conducting membrane transporter by the functional analysis.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani12162015/s1, Figure S1: The secondary structure of 22 tRNAs in the I. gyldenstolpei mitochondrial genome; Figure S2: The secondary structure of 22 tRNAs in the S. indicus mitochondrial genome; Figure S3: The secondary structure of 22 tRNAs in the T. hainanus mitochondrial genome; Figure S4: The protein structure of a selected limbless I. gyldenstolpei gene (ND2 gene) under Branch-site Model analysis and selected sites are also marked in the figure; Table S1: Specific primers used in this study to sequence the complete mitochondrial genomes of the three skink species; Table S2: The best partitioning scheme and best-fitting models were obtained by using the PartitionFinder program. The full names of all abbreviations are as follows: pos1: first codon; pos2: second codon; GTR: general time reversible; I: unchanged site proportion; G: Gamma distribution; Table S3: The localization characteristics of the mitochondrial genomes of I. gyldenstolpei; Table S4: The localization characteristics of the mitochondrial genomes of S. indicus; Table S5: The localization characteristics of the mitochondrial genomes of T. hainanus; Table S6: Codon numbers of the three skinks and relatively synonymous codon usage (RSCU); Table S7: The parameters and results analyzed by the branch-site model and the site model in this study; Table S8: The position and length of the nucleotides and amino acids from the 13 protein-coding genes that were used in this analysis.

Author Contributions

Conceptualization, L.W., Y.T., J.-Y.Z. and D.-N.Y.; data curation, L.W., Y.T., S.P.G.A., K.B.S., J.-Y.Z. and D.-N.Y.; formal analysis, L.W., Y.T., S.P.G.A. and J.-Y.Z.; funding acquisition, D.-N.Y.; investigation, L.W. and K.B.S.; methodology, L.W., Y.T., S.P.G.A., K.B.S., J.-Y.Z. and D.-N.Y.; project administration, S.P.G.A.; writing—original draft, L.W. and Y.T.; writing—review and editing, L.W., Y.T., S.P.G.A., K.B.S., J.-Y.Z. and D.-N.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This study was conducted under a grant from the National Natural Science Foundation of China (31801963). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Institutional Review Board Statement

All experimental designs and animal handling were approved by the Animal Research Ethics Committee of Zhejiang Normal University (ZSDW2021057).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data supporting the findings of this study are openly available from the National Center for Biotechnology Information at https://www.ncbi.nlm.nih.gov (accessed on 22 March 2022). Accession numbers are: MH020638, OM117611, and OM117612.

Acknowledgments

We are deeply grateful to Ziyi Zhang and Jian Chen for their help in revising the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Uetz, P.; Koo, M.S.; Aguilar, R.; Brings, E.; Catenazzi, A.; Chang, A.T.; Chaitanya, R.; Freed, P.; Gross, J.; Hammermann, M.; et al. A quarter century of reptile and amphibian databases. Herpetol. Rev. 2021, 52, 246–255. [Google Scholar]
  2. Grismer, L.L.; Wood, P.L.; Quah, E.S.; Anuar, S.; Poyarkov, N.A.; Thy, N.; Orlov, N.L.; Thammachoti, P.; Seiha, H. Integrative taxonomy of the Asian skinks Sphenomorphus stellatus (Boulenger, 1900) and S. praesignis (Boulenger, 1900) with the resurrection of S. annamiticus (Boettger, 1901) and the description of a new species from Cambodia. Zootaxa 2019, 4683, 381–411. [Google Scholar] [CrossRef] [PubMed]
  3. Greer, A.E. Limb reduction in squamates: Identification of the lineages and discussion of the trends. J. Herpetol. 1991, 95, 166–173. [Google Scholar] [CrossRef]
  4. Greer, A.E. The Biology and Evolution of Scincid Lizards; Surrey Beatty & Sons: Chipping Norton, NSW, Australia, 2007. [Google Scholar]
  5. Greer, A.E.; Wilson, G.D. Comments on the scincid lizard genus Ophiomorus, with a cladistic analysis of the species. Hamadryad 2001, 26, 261–271. [Google Scholar]
  6. Poulakakis, N.; Pakaki, V.; Mylonas, M.; Lymberakis, P. Molecular phylogeny of the Greek legless skink Ophiomorus punctatissimus (Squamata: Scincidae): The impact of the Mid-Aegean trench in its phylogeography. Mol. Phylogenet. Evol. 2008, 47, 396–402. [Google Scholar] [CrossRef]
  7. Schmitz, A.; Brandley, M.C.; Mausfeld, P.; Vences, M.; Glaw, F.; Nussbaum, R.A.; Reeder, T.W. Opening the black box: Phylogenetics and morphological evolution of the Malagasy fossorial lizards of the subfamily “Scincinae”. Mol. Phylogenet. Evol. 2005, 34, 118–133. [Google Scholar] [CrossRef]
  8. Wassersug, R.J.; Roberts, L.; Gimian, J.; Hughes, E.; Saunders, R.; Devison, D.; Woodbury, J.; O’Reilly, J.C. The behavioral responses of amphibians and reptiles to microgravity on parabolic flights. Zoology 2005, 108, 107–120. [Google Scholar] [CrossRef]
  9. Kunya, K.; Kamolnorranath, S.; Chan-ard, T.; Cota, M.; Makchai, S. Additional locality of an endemic legless skink, Isopachys roulei (Angel, 1920) (Squamata, Scincidae) from northeastern Thailand. Thailand Nat. Hist. Mus. J. 2011, 5, 89–93. [Google Scholar]
  10. Grismer, L.L.; Quah, E.S. The rediscovery of Sphenomorphus malayanus Doria, 1888 (Squamata: Scincidae) from the Titiwangsa Mountain range of Peninsular Malaysia and its re-description as S. senja sp. nov. Zootaxa 2015, 3931, 63–70. [Google Scholar] [CrossRef]
  11. Grismer, L.L. Lizards of Peninsular Malaysia, Singapore, and their Adjacent Archipelagos: Their Description, Distribution, and Natural History; Edition Chimaira: Frankfurt am Main, Germany, 2011. [Google Scholar]
  12. Sumarli, A.; Grismer, L.L.; Wood, P.L.; Ahmad, A.B.; Rizal, S.; Ismail, L.H.; Izam, N.A.M.; Ahmad, N.; Linkem, C. The first riparian skink (Genus: Sphenomorphus Strauch, 1887) from Peninsular Malaysia and its relationship to other Indochinese and Sundaic species. Zootaxa 2016, 4173, 029–044. [Google Scholar] [CrossRef]
  13. Barbour, T. Aquatic skincs and arboreal monitors. Copeia 1921, 97, 42–44. [Google Scholar] [CrossRef]
  14. Bauer, A.M.; Jackman, T. Global Diversity of Lizards in Freshwater (Reptilia: Lacertilia). In Freshwater Animal Diversity Assess; Springer: Dordrecht, NL, USA, 2007; pp. 581–586. [Google Scholar]
  15. Guo, K.J.; Shu, F.; Wu, N.F.; Chen, S.D.; Hou, M.; Shi, S.C.; Deng, X.J. Neotype designation and redescription of Tropidophorus guangxiensis Wen, 1992 (Squamata: Sauria: Scincidae), with description of a new subspecies from central South China. Zool. Res. 2021, 42, 606–613. [Google Scholar] [CrossRef]
  16. Zheng, Y.; Wiens, J.J. Combining phylogenomic and supermatrix approaches, and a time-calibrated phylogeny for squamate reptiles (lizards and snakes) based on 52 genes and 4162 species. Mol. Phylogenet. Evol. 2016, 94, 537–547. [Google Scholar] [CrossRef]
  17. Pyron, R.A.; Burbrink, F.T.; Wiens, J.J. A phylogeny and revised classification of Squamata, including 4161 species of lizards and snakes. BMC Evol. Biol. 2013, 13, 93. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Wiens, J.J.; Hutter, C.R.; Mulcahy, D.G.; Noonan, B.P.; Townsend, T.M.; Sites Jr, J.W.; Reeder, T.W. Resolving the phylogeny of lizards and snakes (Squamata) with extensive sampling of genes and species. Biol. Lett. 2012, 8, 1043–1046. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Lambert, S.M.; Reeder, T.W.; Wiens, J.J. When do species-tree and concatenated estimates disagree? An empirical analysis with higher-level scincid lizard phylogeny. Mol. Phylogenet. Evol. 2015, 82, 146–155. [Google Scholar] [CrossRef]
  20. Avise, J.C.; Arnold, J.; Ball, R.M.; Bermingham, E.; Lamb, T.; Neigel, J.E.; Reeb, C.A.; Saunders, N.C. Intraspecific phylogeography: The mitochondrial DNA bridge between population genetics and systematics. Annu. Rev. Ecol. Syst. 1987, 18, 489–522. [Google Scholar] [CrossRef]
  21. Escalante, I.; Ellis, V.R.; Elias, D.O. Leg loss decreases endurance and increases oxygen consumption during locomotion in harvestmen. J. Comp. Physiol. A 2021, 207, 257–268. [Google Scholar] [CrossRef]
  22. Boore, J.L. Animal mitochondrial genomes. Nucleic Acids Res. 1999, 27, 1767–1780. [Google Scholar] [CrossRef] [Green Version]
  23. Das, J. The role of mitochondrial respiration in physiological and evolutionary adaptation. BioEssays 2006, 28, 890–901. [Google Scholar] [CrossRef]
  24. Sun, Y.B.; Shen, Y.Y.; Irwin, D.M.; Zhang, Y.P. Evaluating the roles of energetic functional constraints on teleost mitochondrial-encoded protein evolution. Mol. Biol. Evol. 2011, 28, 39–44. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Shen, Y.Y.; Shi, P.; Sun, Y.B.; Zhang, Y.P. Relaxation of selective constraints on avian mitochondrial DNA following the degeneration of flight ability. Genome Res. 2009, 19, 1760–1765. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Pereira, S.L. Mitochondrial genome organization and vertebrate phylogenetics. Genet. Mol. Biol. 2000, 23, 745–752. [Google Scholar] [CrossRef] [Green Version]
  27. Galtier, N.; Jobson, R.W.; Nabholz, B.; Glémin, S.; Blier, P.U. Mitochondrial whims: Metabolic rate, longevity and the rate of molecular evolution. Biol. Lett. 2009, 5, 413–416. [Google Scholar] [CrossRef]
  28. Shen, Y.Y.; Liang, L.; Zhu, Z.H.; Zhou, W.P.; Irwin, D.M.; Zhang, Y.P. Adaptive evolution of energy metabolism genes and the origin of flight in bats. Proc. Natl. Acad. Sci. USA 2010, 107, 8666–8671. [Google Scholar] [CrossRef] [Green Version]
  29. Chong, R.A.; Mueller, R.L. Low metabolic rates in salamanders are correlated with weak selective constraints on mitochondrial genes. Evol. Int. J. Org. Evol. 2013, 67, 894–899. [Google Scholar] [CrossRef] [PubMed]
  30. Castoe, T.A.; Jiang, Z.J.; Gu, W.; Wang, Z.O.; Pollock, D.D. Adaptive evolution and functional redesign of core metabolic proteins in snakes. PLoS ONE 2008, 3, e2201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Li, X.D.; Jiang, G.F.; Yan, L.Y.; Li, R.; Mu, Y.; Deng, W.A. Positive selection drove the adaptation of mitochondrial genes to the demands of flight and high-altitude environments in grasshoppers. Front. Genet. 2018, 9, 605. [Google Scholar] [CrossRef] [Green Version]
  32. Zhou, X.; Xu, S.; Yang, Y.; Zhou, K.; Yang, G. Phylogenomic analyses and improved resolution of Cetartiodactyla. Mol. Phylogenet. Evol. 2011, 61, 255–264. [Google Scholar] [CrossRef]
  33. Wang, Z.; Wu, W.; Ren, J.; Peng, C.; Jiang, D.; Li, J. Evolution of phenotype and mitochondrial genome reveals limbless and body-elongated Squamates may change their energy basis for locomotion. Asian Herpetol. Res. 2021, 12, 213–220L. [Google Scholar]
  34. Kumazawa, Y.; Endo, H. Mitochondrial genome of the Komodo dragon: Efficient sequencing method with reptile-oriented primers and novel gene rearrangements. DNA Res. 2004, 11, 115–125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Lalitha, S. Primer Premier 5. Biotech Softw. Internet Rep. 2000, 1, 270–272. [Google Scholar] [CrossRef]
  36. Zhang, L.P.; Yu, D.N.; Storey, K.B.; Cheng, H.Y.; Zhang, J.Y. Higher tRNA gene duplication in mitogenomes of praying mantises (Dictyoptera, Mantodea) and the phylogeny within Mantodea. Int. J. Biol. Macromol. 2018, 111, 787–795. [Google Scholar] [CrossRef] [PubMed]
  37. Burland, T.G. DNASTAR’s lasergene sequence analysis software. Bioinf. Methods Protoc. 2000, 132, 71–91. [Google Scholar]
  38. Bernt, M.; Donath, A.; Jühling, F.; Externbrink, F.; Florentz, C.; Fritzsch, G.; Pütz, J.; Middendorf, M.; Stadler, P.F. MITOS: Improved de novo metazoan mitochondrial genome annotation. Mol. Phylogenet. Evol. 2013, 69, 313–319. [Google Scholar] [CrossRef]
  39. Lowe, T.M.; Chan, P.P. tRNAscan-SE On-line: Integrating search and context for analysis of transfer RNA genes. Nucleic Acids Res. 2016, 44, W54–W57. [Google Scholar] [CrossRef]
  40. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  41. Grant, J.R.; Stothard, P. The CGView Server: A comparative genomics tool for circular genomes. Nucleic Acids Res. 2008, 36, W181–W184. [Google Scholar] [CrossRef]
  42. Zhang, D.; Gao, F.; Jakovlić, I.; Zou, H.; Zhang, J.; Li, W.X.; Wang, G.T. PhyloSuite: An integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies. Mol. Ecol. Resour. 2020, 20, 348–355. [Google Scholar] [CrossRef]
  43. Golding, M. Real World Adobe Illustrator CS4; Peachpit Press: Berkeley, CA, USA, 2008. [Google Scholar]
  44. Chen, L.; Lin, Y.; Xiao, Q.; Lin, Y.; Du, Y.; Lin, C.; Ward-Fear, G.; Hu, C.; Qu, Y.; Li, H. Characterization of the complete mitochondrial genome of the many-lined sun skink (Eutropis multifasciata) and comparison with other Scincomorpha species. Genomics 2021, 113, 2526–2536. [Google Scholar] [CrossRef]
  45. Zhong, J.J.; Wu, Q.Q.; Wang, Y.M.; Guo, K.J.; Ding, G.H.; Luo, S.T. The first complete mitochondrial DNA of the Chinese short-limbed skink (Ateuchosaurus chinensis Gray, 1845) determined by next-generation sequencing. Mitochondrial DNA Part B 2021, 6, 995–996. [Google Scholar] [CrossRef] [PubMed]
  46. Zhang, C.; Sun, X.; Chen, L.; Xiao, W.; Zhu, X.; Xia, Y.; Chen, J.; Wang, H.; Zhang, B. The complete mitochondrial genome of Eumeces chinensis (Squamata: Scincidae) and implications for Scincidae taxonomy. Mitochondrial DNA Part A 2016, 27, 4691–4692. [Google Scholar] [CrossRef] [PubMed]
  47. Song, T.; Zhang, C.; Huang, X.; Zhang, B. Complete mitochondrial genome of Eumeces elegans (Squamata: Scincidae). Mitochondrial DNA Part A 2016, 27, 719–720. [Google Scholar] [CrossRef] [PubMed]
  48. Wu, N.; Cai, B.; Chen, M.; Guo, X. Next-generation sequencing yields a nearly complete mitochondrial genome of Plestiodon liui (Reptilia, Squamata, Scincidae) endemic to China. Mitochondrial DNA Part B 2020, 5, 3637–3638. [Google Scholar] [CrossRef] [PubMed]
  49. Chen, M.; Liu, J.; Chen, D.; Guo, X. The complete mitochondrial genome of a blue-tailed skink (Plestiodon tunganus) endemic to Sichuan Basin. Mitochondrial DNA Part B 2019, 4, 1109–1110. [Google Scholar] [CrossRef] [Green Version]
  50. Park, J.; Koo, K.S.; Kim, I.H.; Park, D. Complete mitochondrial genomes of Scincella vandenburghi and S. huanrenensis (Squamata: Scincidae). Mitochondrial DNA Part B 2016, 1, 237–238. [Google Scholar] [CrossRef] [Green Version]
  51. Zhong, J.; Ma, L.; Guo, K.; Du, Y. Complete mitochondrial genome of Scincella modesta (Squamata: Scincidae) and phylogenetic analysis. Mitochondrial DNA Part B 2020, 5, 373–374. [Google Scholar] [CrossRef] [Green Version]
  52. Chen, L.; Lin, Y.; Lin, Y.; Hu, Y.; Xiao, Q.; Zhou, X.; Liu, Y.; Li, H. Complete mitochondrial genome of Scincella modesta (Squamata: Scincidae). Mitochondrial DNA Part B 2020, 5, 676–677. [Google Scholar] [CrossRef] [Green Version]
  53. Zhong, J.; Guo, K.; Ma, L. First record of mitochondrial genome of Scincella reevesii (Squamata: Scincidae) and phylogenetic analysis. Mitochondrial DNA Part B 2021, 6, 564–565. [Google Scholar] [CrossRef]
  54. Tang, X.S.; Yang, D.C.; Huang, S. The complete mitochondrial genome of Sphenomorphus incognitus (Reptilia: Scincidae). Mitochondrial DNA Part B 2019, 4, 307–308. [Google Scholar] [CrossRef] [Green Version]
  55. Tang, X.S.; Yang, D.C.; Lin, Y.J.; Dai, L.L. The complete mitochondrial genome of Sphenomorphus indicus (Reptilia: Scincidae). Mitochondrial DNA Part B 2019, 4, 2727–2728. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Duangjai, S.; Srisodsuk, S.; Chuaynkern, C.; Chuaynkern, Y. Complete mitochondrial genome of Tropidophorus hangnam (Squamata: Scincidae) with phylogenetic analysis. Mitochondrial DNA Part B 2020, 5, 3683–3684. [Google Scholar] [CrossRef] [PubMed]
  57. Kumazawa, Y. Mitochondrial DNA sequences of five squamates: Phylogenetic affiliation of snakes. DNA Res. 2004, 11, 137–144. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Kumazawa, Y. Mitochondrial genomes from major lizard families suggest their phylogenetic relationships and ancient radiations. Gene 2007, 388, 19–26. [Google Scholar] [CrossRef]
  59. Weng, S.Y.; Huang, S.; Yang, L.; Peng, L.F.; Wei, C.; Li, J.C. The complete mitochondrial genome of Asymblepharus himalayanus (Scincidae; Sauria) and its phylogenetic position within Scincidae. Mitochondrial DNA Part B 2021, 6, 3031–3032. [Google Scholar] [CrossRef]
  60. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [Green Version]
  61. Castresana, J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol. 2000, 17, 540–552. [Google Scholar] [CrossRef] [Green Version]
  62. Lanfear, R.; Calcott, B.; Ho, S.Y.W.; Guindon, S. PartitionFinder: Combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol. Biol. Evol. 2012, 29, 1695–1701. [Google Scholar] [CrossRef] [Green Version]
  63. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef]
  64. Ronquist, F.; Teslenko, M.; Van Der Mark, P.; Ayres, D.L.; Darling, A.; Höhna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model apace. Syst. Biol. 2012, 61, 539–542. [Google Scholar] [CrossRef] [Green Version]
  65. Yang, Z. PAML 4: Phylogenetic analysis by Maximum Likelihood. Mol. Biol. Evol. 2007, 24, 1586–1591. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Gao, F.; Chen, C.; Arab, D.A.; Du, Z.; He, Y.; Ho, S.Y. EasyCodeML: A visual tool for analysis of selection using CodeML. Ecol. Evol. 2019, 9, 3891–3898. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Consortium, U. UniProt: The universal protein knowledgebase. Nucleic Acids Res. 2018, 46, D115–D119. [Google Scholar]
  68. Waterhouse, A.; Bertoni, M.; Bienert, S.; Studer, G.; Tauriello, G.; Gumienny, R.; Heer, F.T.; de Beer, T.A.; Rempfer, C.; Bordoli, L.; et al. SWISS-MODEL: Homology modelling of protein structures and complexes. Nucleic Acids Res. 2018, 46, W296–W303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Chen, M.; Liu, J.; Chen, D.; Guo, X. Complete mitochondrial genome of a blue-tailed skink Plestiodon capito (Reptilia, Squamata, Scincidae) and comparison with other Scincidae lizards. Genetica 2020, 148, 229–241. [Google Scholar] [CrossRef]
  70. Brandley, M.C.; OTA, H.; Hikida, T.; Nieto Montes De Oca, A.; Feria-Ortiz, M.; Guo, X.; Wang, Y. The phylogenetic systematics of blue-tailed skinks (Plestiodon) and the family Scincidae. J. Linn. Soc. Lond. Zool. 2012, 165, 163–189. [Google Scholar] [CrossRef] [Green Version]
  71. Pyron, R.A. Novel approaches for phylogenetic inference from morphological data and total-evidence dating in squamate reptiles (lizards, snakes, and amphisbaenians). Syst. Biol. 2017, 66, 38–56. [Google Scholar] [CrossRef]
  72. Hedges, S.B. The high-level classification of skinks (Reptilia, Squamata, Scincomorpha). Zootaxa 2014, 3765, 317–338. [Google Scholar] [CrossRef] [Green Version]
  73. Honda, M.; Ota, H.; Murphy, R.W.; Hikida, T. Phylogeny and biogeography of water skinks of the genus Tropidophorus (Reptilia: Scincidae): A molecular approach. Zool. Scr. 2006, 35, 85–95. [Google Scholar] [CrossRef]
  74. Cho, S.; Zwick, A.; Regier, J.C.; Mitter, C.; Cummings, M.P.; Yao, J.; Du, Z.; Zhao, H.; Kawahara, A.Y.; Weller, S.; et al. Can deliberately incomplete gene sample augmentation improve a phylogeny estimate for the advanced moths and butterflies (Hexapoda: Lepidoptera)? Syst. Biol. 2011, 60, 782–796. [Google Scholar] [CrossRef] [Green Version]
  75. Wiens, J.J.; Fetzner Jr, J.W.; Parkinson, C.L.; Reeder, T.W. Hylid frog phylogeny and sampling strategies for speciose clades. Syst. Biol. 2005, 54, 778–807. [Google Scholar] [CrossRef] [Green Version]
  76. Zhang, J.; Nielsen, R.; Yang, Z. Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level. Mol. Biol. Evol. 2005, 22, 2472–2479. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Dotson, E.M.; Beard, C.B. Sequence and organization of the mitochondrial genome of the Chagas disease vector, Triatoma dimidiata. Insect Mol. Biol. 2001, 10, 205–215. [Google Scholar] [CrossRef] [PubMed]
  78. Brandt, U. Energy converting NADH: Quinone oxidoreductase (complex I). Annu. Rev. Biochem. 2006, 75, 69–92. [Google Scholar] [CrossRef]
  79. Wirth, C.; Brandt, U.; Hunte, C.; Zickermann, V. Structure and function of mitochondrial complex I. Biochim. Biophys. Acta Biosyst. 2016, 1857, 902–914. [Google Scholar] [CrossRef] [PubMed]
  80. Bridges, H.R.; Birrell, J.A.; Hirst, J. The mitochondrial-encoded subunits of respiratory complex I (NADH: Ubiquinone oxidoreductase): Identifying residues important in mechanism and disease. Biochem. Soc. Trans. 2011, 39, 799–860. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Mathiesen, C.; Hägerhäll, C. Transmembrane topology of the NuoL, M and N subunits of NADH: Quinone oxidoreductase and their homologues among membrane-bound hydrogenases and bona fide antiporters. Biochim. Biophys. Acta Biosyst. 2002, 1556, 121–132. [Google Scholar] [CrossRef] [Green Version]
  82. Da Fonseca, R.R.; Johnson, W.E.; O’Brien, S.J.; Ramos, M.J.; Antunes, A. The adaptive evolution of the mammalian mitochondrial genome. BMC Genom. 2008, 9, 119. [Google Scholar] [CrossRef] [Green Version]
  83. Gross, M.R. Evolution of diadromy in fishes. Am. Fish. Soc. Symp. 1987, 1, 14–25. [Google Scholar]
Animals 12 02015 g001 550
Figure 1. Circular visualization maps of the complete mitochondrial genomes of I. gyldenstolpei, T. hainanus, and S. indicus (since the mitochondrial genomes of the three skinks are similar in structure, a map was used to represent the three skinks). The three circles from the outside to the inside show the gene map (PCGs, rRNAs, tRNAs, and the AT-rich region), the GC content, and GC skew, respectively. Among them, the genes outside the map are coded on the majority strand (J-strand), whereas the genes inside the map are coded on the minority strand (N-strand).
Figure 1. Circular visualization maps of the complete mitochondrial genomes of I. gyldenstolpei, T. hainanus, and S. indicus (since the mitochondrial genomes of the three skinks are similar in structure, a map was used to represent the three skinks). The three circles from the outside to the inside show the gene map (PCGs, rRNAs, tRNAs, and the AT-rich region), the GC content, and GC skew, respectively. Among them, the genes outside the map are coded on the majority strand (J-strand), whereas the genes inside the map are coded on the minority strand (N-strand).
Animals 12 02015 g001
Animals 12 02015 g002 550
Figure 2. The relative synonymous codon usage (RSCU) of the mitochondrial genomes of I. gyldenstolpei (A), S. indicus (B), and T. hainanus (C). Acronyms stand for different amino acids. The x-axis represents all codons used and different combinations of synonymous codons, and the RSCU values are listed on the y-axis.
Figure 2. The relative synonymous codon usage (RSCU) of the mitochondrial genomes of I. gyldenstolpei (A), S. indicus (B), and T. hainanus (C). Acronyms stand for different amino acids. The x-axis represents all codons used and different combinations of synonymous codons, and the RSCU values are listed on the y-axis.
Animals 12 02015 g002
Animals 12 02015 g003 550
Figure 3. Phylogenetic relationships among 18 species of Scincidae based on the nucleotide dataset of the 13 mitochondrial protein-coding genes. Smaug warreni and Lepidophyma flavimaculatum were used as the outgroups. The numbers above the branches specify bootstrap percentages from ML (left) and posterior probabilities as determined from BI (right). The GenBank accession numbers of all species are shown in the figure. Different colors represent different subfamilies.
Figure 3. Phylogenetic relationships among 18 species of Scincidae based on the nucleotide dataset of the 13 mitochondrial protein-coding genes. Smaug warreni and Lepidophyma flavimaculatum were used as the outgroups. The numbers above the branches specify bootstrap percentages from ML (left) and posterior probabilities as determined from BI (right). The GenBank accession numbers of all species are shown in the figure. Different colors represent different subfamilies.
Animals 12 02015 g003
Table
Table 2. Base composition of the mitochondrial genomes of the three skinks.
Table 2. Base composition of the mitochondrial genomes of the three skinks.
RegionI. gyldenstolpeiS. indicusT. hainanus
Length (bp)A + T (%)AT SkewGC SkewLength (bp)A + T (%)AT SkewGC SkewLength (bp)A + T (%)AT SkewGC Skew
mito 17,210560.121−0.35816,94457.90.11−0.30817,00154.20.118−0.33
PCGsJ11,38255.520.058−0.39211,38257.280.056−0.3511,37953.060.063−0.363
N−0.6510.602−0.510.508−0.5950.462
rRNAs 250754.20.313−0.207247057.40.29−0.147248553.70.296−0.191
A + T-rich region 177461.50.086−0.46156661.10.9−0.359161560.40.07−0.404
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wu, L.; Tong, Y.; Ayivi, S.P.G.; Storey, K.B.; Zhang, J.-Y.; Yu, D.-N. The Complete Mitochondrial Genomes of Three Sphenomorphinae Species (Squamata: Scincidae) and the Selective Pressure Analysis on Mitochondrial Genomes of Limbless Isopachys gyldenstolpei. Animals 2022, 12, 2015. https://doi.org/10.3390/ani12162015

AMA Style

Wu L, Tong Y, Ayivi SPG, Storey KB, Zhang J-Y, Yu D-N. The Complete Mitochondrial Genomes of Three Sphenomorphinae Species (Squamata: Scincidae) and the Selective Pressure Analysis on Mitochondrial Genomes of Limbless Isopachys gyldenstolpei. Animals. 2022; 12(16):2015. https://doi.org/10.3390/ani12162015

Chicago/Turabian Style

Wu, Lian, Yao Tong, Sam Pedro Galilee Ayivi, Kenneth B. Storey, Jia-Yong Zhang, and Dan-Na Yu. 2022. "The Complete Mitochondrial Genomes of Three Sphenomorphinae Species (Squamata: Scincidae) and the Selective Pressure Analysis on Mitochondrial Genomes of Limbless Isopachys gyldenstolpei" Animals 12, no. 16: 2015. https://doi.org/10.3390/ani12162015

APA Style

Wu, L., Tong, Y., Ayivi, S. P. G., Storey, K. B., Zhang, J. -Y., & Yu, D. -N. (2022). The Complete Mitochondrial Genomes of Three Sphenomorphinae Species (Squamata: Scincidae) and the Selective Pressure Analysis on Mitochondrial Genomes of Limbless Isopachys gyldenstolpei. Animals, 12(16), 2015. https://doi.org/10.3390/ani12162015

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop