# Mitochondrial Genomes of Six Snakes (Lycodon) and Implications for Their Phylogeny

**Authors:** Fei Zhu, Anqiong Lu, Ke Sun

PMC · DOI: 10.3390/genes16050493 · Genes · 2025-04-26

## TL;DR

This study sequenced the mitochondrial genomes of six snake species to better understand their evolutionary relationships.

## Contribution

The study provides the first complete mitochondrial genome sequences for six Lycodon snake species and analyzes their phylogenetic relationships.

## Key findings

- The six snake species had mitochondrial genomes ranging from 17,143 to 17,298 bp in length.
- The atp8 gene was the fastest evolving, while cox1 was the most conserved.
- Phylogenetic analyses grouped the six species together with their congeneric relatives.

## Abstract

Background: Colubridae, known to be one of the most species-rich snake families, remains relatively understudied in termshe context of complete mitochondrial genome research. This study provide the first systematic characterization of the mitochondrial genomes of six colubrid species: Lycodon subcinctus, Lycodon rosozonatus, Lycodon fasciatus, Lycodon gongshan, Lycodon futsingensis, and Lycodon aulicus. Method: In this study, mitochondrial genomes were sequenced using Sanger sequencing. The raw data were subjected to quality- filtered withing using Fastp and subsequently assembled into complete mitochondrial genomes via SPAdes. Gene annotation was performed by Tblastn, Genewise (for CDS coding sequences), MiTFi (for transfer RNAs), and Rfam (for ribosomal RNAs). Sequence analyses were conducted with various tools, including MEGA, tRNAscan-SE, DnaSP, MISA, and REPuter. Finally, phylogenetic trees were reconstructed based on 13 protein-coding genes from 14 species. Results:The mitogenomes of these six species ranged from 17,143 to 17,298 bp in length and con-sisted of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and 2 control regions. The nucleotide composition of the Colu-bridae mitogenomes was comparable with an A + T composition ranging from 52.1% to 58.8% except for the trnS1 and trnC. All the tRNAs could fold into a stable secondary structure. The Pi and Ka/Ks values suggested that atp8 was the fastest-evolving gene, while cox1 was the most conserved gene. Bayesian inference and maximum likelihood phylogenetic analyses yielded consistent results, with the six sequenced species clus-tering together with their congeneric species. These findings will provide valuable references for further research on the phylogeny of Colubridae.

## Linked entities

- **Genes:** ATP8 (ATP synthase F0 subunit 8) [NCBI Gene 4509], COX1 (cytochrome c oxidase subunit I) [NCBI Gene 4512]
- **Species:** Lycodon subcinctus (taxon 1301316), Lycodon rosozonatus (taxon 2483678), Lycodon fasciatus (taxon 591076), Lycodon gongshan (taxon 1677021), Lycodon futsingensis (taxon 591077), Lycodon aulicus (taxon 1125755)

## Full-text entities

- **Genes:** COX1 (cytochrome c oxidase subunit I) [NCBI Gene 4512] {aka COI, MTCO1}, ATP8 (ATP synthase F0 subunit 8) [NCBI Gene 4509] {aka ATPase8, MTATP8}
- **Species:** Lycodon subcinctus (Malayan banded wolf snake, species) [taxon 1301316], Lycodon aulicus (Indian wolf snake, species) [taxon 1125755], Lycodon gongshan (Gongshan wolf snake, species) [taxon 1677021], Lycodon rosozonatus (species) [taxon 2483678], Lycodon fasciatus (species) [taxon 591076], Lycodon futsingensis (species) [taxon 591077]

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12110799/full.md

## Figures

4 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12110799/full.md

## References

70 references — full list in the complete paper: https://tomesphere.com/paper/PMC12110799/full.md

---
Source: https://tomesphere.com/paper/PMC12110799