Next-generation sequencing yields a nearly complete mitochondrial genome of the Forsyth’s toad-headed agama, Phrynocephalus forsythii (Reptilia, Squamata, Agamidae)

Species divergence in valleys: the phylogeny of Phrynocephalus forsythii complex and description of a new species

Background

Geographic isolation caused by high-altitude valleys promotes the formation of geographic segregation of species, leading to species differentiation. The subgenus Oreosaura contains viviparous species from the Tibetan Plateau and the vicinity of the Tarim Basin, which can be divided into three species complexes according to their geographical distribution: Phrynocephalus vlangalii, Phrynocephalus theobaldi, and Phrynocephalus forsythii. However, molecular data for the P. forsythii complex are limited and the diversity of this species complex has been greatly underestimated. Therefore, this study aimed to explore the species diversity of Oreosaura and species differentiation within the P. forsythii complex.

Methods

We analysed the species diversity of Oreosaura by combining previous data, constructed a phylogenetic tree of the subgenus based on cytochrome c oxidase subunit I and 16S sequences, and estimated the divergence time.

Results

The results suggest significant genetic differences between the Tarim Basin populations and adjacent mountain valley populations of the P. forsythii complex and that the combination of deep valley landscapes in the high mountains and ice-age events have contributed to the differentiation of the viviparous toad-headed agama lizard, which is a key factor in the phylogenetics of the P. forsythii complex. Furthermore, we identified a population collected from Wuqia County, Xinjiang, as a new species, Phrynocephalus kangsuensis sp. nov. The results will provide data for phylogenetic studies following the P. forsythii complex and help demonstrate that valleys promote the formation of Phrynocephalus species.

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.

Next-generation sequencing yields a complete mitochondrial genome of the Tartar sand boa (Eryx tataricus) from Junggar Basin

The viviparous Tartar sand boa, Eryx tataricus, is a widespread species in arid Central Asia. A complete mitochondrial genome of one individual from Junggar Basin in Northwest China was determined by next-generation sequencing. The mitogenome is 17,537 bp in size, comprising 2 ribosomal RNA genes, 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and 2 control regions. The order and structure of the genes are similar to those of other Henophidia snakes. Phylogenetic analysis based on 13 concatenated PCGs recovered the monophyly of Boidae and indicated that E. tataricus is closely related to Boa constrictor plus Eunectes notaeus.