Phylogenetic relationships and divergence dating in the Glass Lizards (Anguinae)

The Phylogenetic Relationships of Major Lizard Families Using Mitochondrial Genomes and Selection Pressure Analyses in Anguimorpha

Anguimorpha, within the order Squamata, represents a group with distinct morphological and behavioral characteristics in different ecological niches among lizards. Within Anguimorpha, there is a group characterized by limb loss, occupying lower ecological niches, concentrated within the subfamily Anguinae. Lizards with limbs and those without exhibit distinct locomotor abilities when adapting to their habitats, which in turn necessitate varying degrees of energy expenditure. Mitochondria, known as the metabolic powerhouses of cells, play a crucial role in providing approximately 95% of an organism's energy. Functionally, mitogenomes (mitochondrial genomes) can serve as a valuable tool for investigating potential adaptive evolutionary selection behind limb loss in reptiles. Due to the variation of mitogenome structures among each species, as well as its simple genetic structure, maternal inheritance, and high evolutionary rate, the mitogenome is increasingly utilized to reconstruct phylogenetic relationships of squamate animals. In this study, we sequenced the mitogenomes of two species within Anguimorpha as well as the mitogenomes of two species in Gekkota and four species in Scincoidea. We compared these data with the mitogenome content and evolutionary history of related species. Within Anguimorpha, between the mitogenomes of limbless and limbed lizards, a branch-site model analysis supported the presence of 10 positively selected sites: Cytb protein (at sites 183 and 187), ND2 protein (at sites 90, 155, and 198), ND3 protein (at site 21), ND5 protein (at sites 12 and 267), and ND6 protein (at sites 72 and 119). These findings suggested that positive selection of mitogenome in limbless lizards may be associated with the energy requirements for their locomotion. Additionally, we acquired data from 205 mitogenomes from the NCBI database. Bayesian inference (BI) and Maximum Likelihood (ML) trees were constructed using the 13 mitochondrial protein-coding genes (PCGs) and two rRNAs (12S rRNA and 16S rRNA) from 213 mitogenomes. Our phylogenetic tree and the divergence time estimates for Squamata based on mitogenome data are consistent with results from previous studies. Gekkota was placed at the root of Squamata in both BI and ML trees. However, within the Toxicofera clade, due to long-branch attraction, Anguimorpha and (Pleurodonta + (Serpentes + Acrodonta)) were closely related groupings, which might indicate errors and also demonstrate that mitogenome-based phylogenetic trees may not effectively resolve long-branch attraction issues. Additionally, we reviewed the origin and diversification of Squamata throughout the Mesozoic era, suggesting that Squamata originated in the Late Triassic (206.05 Mya), with the diversification of various superfamilies occurring during the Cretaceous period. Future improvements in constructing squamate phylogenetic relationships using mitogenomes will rely on identifying snake and acrodont species with slower evolutionary rates, ensuring comprehensive taxonomic coverage of squamate diversity, and increasing the number of genes analyzed.

Acknowledging more biodiversity without more species

Delimiting and naming biodiversity is a vital step toward wildlife conservation and research. However, species delimitation must be consistent across biota so that the limited resources available for nature protection can be spent effectively and objectively. To date, newly discovered lineages typically are either left undescribed and thus remain unprotected or are being erroneously proposed as new species despite mixed evidence for completed speciation, in turn contributing to the emerging problem of taxonomic inflation. Inspired by recent conceptual and methodological progress, we propose a standardized workflow for species delimitation that combines phylogenetic and hybrid zone analyses of genomic datasets ("genomic taxonomy"), in which phylogeographic lineages that do not freely admix are ranked as species, while those that have remained fully genetically compatible are ranked as subspecies. In both cases, we encourage their formal taxonomic naming, diagnosis, and description to promote social awareness toward biodiversity. The use of loci throughout the genome overcomes the unreliability of widely used barcoding genes when phylogeographic patterns are complex, while the evaluation of divergence and reproductive isolation unifies the long-opposed concepts of lineage species and biological species. We suggest that a shift in conservation assessments from a single level (species) toward a two-level hierarchy (species and subspecies) will lead to a more balanced perception of biodiversity in which both intraspecific and interspecific diversity are valued and more adequately protected.

Phylogenomics of Anguis and Pseudopus (Squamata, Anguidae) indicates Balkan-Apennine mitochondrial capture associated with the Messinian event

A dated phylogenetic hypothesis on the evolutionary history of the extant taxa of the Western Palearctic lizards Anguis and Pseudopus is revised using genome-wide nuclear DNA and mitogenomes. We found overall concordance between nuclear and mitochondrial DNA phylogenies, with one significant exception - the Apennine A. veronensis. In mitochondrial DNA, this species forms a common clade with the earliest diverging lineage, the southern Balkan endemic A. cephallonica, while it clusters together with A. fragilis in nuclear DNA. The nuclear phylogeny conforms to the morphology, which is relatively similar between A. veronensis and A. fragilis. The most plausible explanation for the mitonuclear discordance is ancient mitochondrial capture from the Balkan ancestor of A. cephallonica to the Apennine population of the A. fragilis-veronensis ancestor. We hypothesize that this capture occurred only in a geographically restricted population. The dating of this presumed mitochondrial introgression and capture coincides with the Messinian event, when the Balkan and Apennine Peninsulas were presumably largely connected. The dated nuclear phylogenomic reconstruction estimated the divergence of A. cephallonica around 12 Mya, while the sister clade representing the A. fragilis species complex consisting of the sister species A. fragilis-A. veronensis and A. colchica-A. graeca further diversified around 7 Mya. The depth of nuclear divergence among the evolutionary lineages of Pseudopus (0.5-1.2 Mya) supports their subspecies status.

Phylogenomics and historical biogeography of the cleptoparasitic bee genus Nomada (Hymenoptera: Apidae) using ultraconserved elements

The genus Nomada Scopoli (Hymenoptera: Apidae) is the largest genus of brood parasitic bees with nearly 800 species found across the globe and in nearly all biogeographic realms except Antarctica. There is no previous molecular phylogeny focused on Nomada despite their high species abundance nor is there an existing comprehensive biogeography for the genus. Using ultraconserved element (UCE) phylogenomic data, we constructed the first molecular phylogeny for the genus Nomada and tested the monophyly of 16 morphologically established species groups. We also estimated divergence dates using fossil calibration points and inferred the origin and of this genus around the globe. Our phylogeny recovered 14 of the 16 previously established species groups as monophyletic. The superba and ruficornis groups, however, were recovered as non-monophyletic and need to be re-evaluated using morphology. Divergence dating and historic biogeographic analyses performed on the phylogenetic reconstruction indicates that Nomada most likely originated in the Holarctic ∼65 Mya. Geodispersal into the southern hemisphere occurred three times; once during the Eocene into the Afrotropics, once during the Oligocene into the Neotropics, and once during the Miocene into Australasia. Geodispersal across the Holarctic was most frequent and occurred repeatedly throughout the Cenozoic era, using the De Geer, Thulean, and the Bering Land Bridges. This is the first instance of a bee using both the Thulean and De Geer land bridges and has implications of how early bee species dispersed throughout the Palearctic in the late Cretaceous and early Paleogene.

Phenotypic differentiation of the slow worm lizards (Squamata: Anguis) across their contact zone in Central Europe

Background

The application of molecular-phylogenetic approaches to taxonomy has had a dramatic effect on our understanding of the diversity of reptiles. These approaches have allowed researchers to reveal previously hidden lineages as well as taxonomic overestimation in morphologically plastic taxa. Slow worms, legless lizards of the genus Anguis (Squamata: Anguidae), were previously considered to comprise either one or two species, and morphology-based intraspecific taxonomy of Anguis fragilis remained controversial throughout the 20th century. After the discovery of deep genetic divergences within the genus, its taxonomy was reconsidered, and as a result, five extant species have been recognized. In order to better understand the patterns of their interspecific differentiation, here we studied phenotypic differences between the two most widespread of them—A. fragilis and A. colchica, and their putative hybrids across the contact zone of both species in Central Europe.

Methods

We used multivariate and univariate statistics and analyzed ten metric, eleven meristic, and six categorical phenotypic variables in material comprising a total of 326 individuals. We also genotyped individuals from the contact zone for one mitochondrial and two nuclear DNA fragments in order to delineate the distribution of individuals of hybrid and non-hybrid origin. The clines in morphological traits were studied using HZAR.

Results

We show that the two species are morphologically differentiated. Anguis fragilis has a less robust head, fewer scales covering the body, lower frequency of the external ear opening presence, lower frequency of separated prefrontal scales, higher frequency of prefrontal scales in contact with each other, and body coloration more similar to the juvenile coloration than A. colchica. Slow worms from the contact/hybrid zone are characterized by an intermediate morphology, with more similarities to A. fragilis than to A. colchica.

Discussion

None of the analyzed characters alone proved to be fully diagnostic, although more than 90% of all individuals could be successfully assigned to one or another species based on numbers of scales around the body. Our results indicate concordant, coincident, and steep clines in character states change. We present several hypotheses on the origin and evolutionary maintenance of the morphological divergence between both species and suggest that different evolutionary histories of the taxa rather than recently acting selection explain the observed morphological variation.

Comparative anatomy and ontogeny of appendicular skeleton of Pseudopus apodus (Pallas, 1775) (Anguimorpha, Anguidae) and a pattern of hindlimb loss in Anguinae

We present a detailed anatomy of the pectoral girdles, pelvic girdles, and hindlimbs of adult and juvenile specimens of Pseudopus apodus (Pallas, 1775). We compared the individual bones of the appendicular skeleton of P. apodus with those of Anguis fragilis and species of Ophisaurus living in North America, North Africa, and Southeast Asia. We found no anatomical features in P. apodus in common with the species of Ophisaurus living in only North America, North Africa, and Southeast Asia. Additionally, we present the prehatching ontogeny of the pelvic girdle of P. apodus and A. fragilis and the prehatching ontogeny of the hindlimb of P. apodus. In the ontogeny of the pelvic girdle of P. apodus, it is possible to distinguish the ossification centers of ilium, ischium, and pubis. In contrast, in the ontogeny of A. fragilis, no ossification centers of ilium, ischium, and pubis are present, and no hindlimb element was detected. In Stage 1 of ontogeny in Pseudopus, the femur and tibia are present; in Stage 2, the nodule representing the fibula appears; and in Stage 3, in addition to the femur, tibia, and fibula, four tarsal elements are present. This anatomical condition corresponds to the anatomical composition of the hindlimb of the adult O. koellikeri. In Stage 4, the involution of all tarsal elements and fibula begins, and in the last two prehatching stages, only femur and tibia remain; this condition is present not only in the adults of Pseudopus, but also in those of several other species of Ophisaurus.

Morphological and genetic differentiation in the anguid lizard Pseudopus apodus supports the existence of an endemic subspecies in the Levant

The Levant represents one of the most important reptile diversity hotspots and centers of endemism in the Western Palearctic. The region harbored numerous taxa in glacial refugia during the Pleistocene climatic oscillations. Due to the hostile arid conditions in the warmer periods they were not always able to spread or come into contact with populations from more distant regions. One large and conspicuous member of the Levantine herpetofauna is the legless anguid lizardPseudopus apodus. This species is distributed from the Balkans to Central Asia with a portion of its range running along the eastern Mediterranean coast. Mitochondrial and nuclear DNA sequences, microsatellite genotypes, and morphology show that populations in this region differ from the two named subspecies and presumably had a long independent evolutionary history during the Quaternary. Here we describe the Levantine population as a new subspecies and present biogeographic scenarios for its origin and diversification. The new subspecies is genetically highly diverse, and it forms a sister lineage toPseudopusfrom the remaining parts of the range according to mtDNA. It is the largest-bodied of the three subspecies, but occupies the smallest range.

DNA barcoding reveals distinct population of Dopasia gracilis (Squamata: Anguidae) in Mizoram, Northeast India

The DNA barcode data of Asian Glass Lizard, Dopasia gracilis, is limited in the global database, especially from India. The present study aimed to generate a barcode sequence of morphologically identified D. gracilis from the Mizoram state in northeast India and compared with other Anguidae species. The studied species showed monophyletic clustering in the Bayesian analysis (BA) phylogeny with strong posterior probability support and also discriminated sufficient Kimura 2 parameter genetic distances. The barcode data of D. gracilis revealed high intra-species genetic variability and formed two clusters in BA phylogeny. The Templeton, Crandall, and Sing network also depicted four different haplotypes within the barcode sequences of D. gracilis. The DNA sequences generated from northeast India showed 6.5-6.6% and 7.3% genetic distances with the sequences generated from Yunnan Province and Tibetan Plateau, respectively. Considering the high genetic distances, multiple clustering, and distinct haplotypes, the present study assumed the presence of possible cryptic diversity of D. gracilis in the Indochina sub-region and a distinct population in northeast India. We recommended the generation of more DNA information from different localities to elucidate the actual diversity of D. gracilis within the known range distribution.

An assessment of vouchered records and field observations of the rare anguid, Dopasia buettikoferi (Lidth de Jeude, 1905) in Borneo

We report recent observations, two new locality records and an updated overview of the distribution of Dopasia buettikoferi (Lidth de Jeude, 1905) from Borneo, as well as photographic documentation of this species showing colouration and patterns in living individuals and the holotype. This note represents the first complete compilation of distributional knowledge of this rarely-encountered member of the family Anguidae.

Next-generation sequencing yields a complete mitochondrial genome of the Asian Glass Lizard (Dopasia gracillis) from the Yungui Plateau in Southwest China

The Asian Glass Lizard, Dopasia gracillis, has wide distribution in North India, Nepal, South China, and Indochina. In this study, a complete mitochondrial genome of D. gracillis from the Yungui Plateau in Southwest China was determined by next-generation sequencing. Similar to the typical mtDNA of vertebrates, the mitogenome was 17,133 bp in length and comprised the standard set of 13 protein-coding genes (PCGs), 22 tRNA genes, 2 rRNA genes, and 1 control region. The concatenated PCGs were used to conduct Bayesian phylogenetic analyses together with mitogenome data of Anguidae and related taxa in GenBank. The resulting phylogenetic tree confirmed the monophyly of Anguidae and Aguinae as well as D. gracillis, respectively. The mitogenome reported here will contribute to the examination of phylogeographic structure for D. gracillis and understanding of mitochondrial DNA evolution in Anguidae.

The complete mitochondrial genome of the Hainan Glass Lizard (Dopasia hainanensis) determined by next-generation sequencing

The complete mitochondrial genome of the Hainan Glass Lizard (Dopasia hainanensis) from its type locality (Diaoluo Mountain in Hainan Island, China) was determined using next-generation sequencing. The mitogenome was 17,000 bp in length and comprised the standard set of 1 control region, 2 rRNAs, 22 tRNAs, plus 13 protein-coding genes (PCGs). The PCGs were used to perform Bayesian phylogenetic analysis together with other Anguidae and Helodermatidae as well as Shinisauridae lizards with mitogenome in GenBank. The resulting phylogenetic tree recovered D. hainanensis as the sister-taxon to Dopasia harti. The mitogenome of D. hainanensis will provide a valuable resource for various study areas such as species delimitation, molecular evolution, and phylogenetic inference.