Molecular systematics of the Philippine forest skinks (Squamata: Scincidae: Sphenomorphus): testing morphological hypotheses of interspecific relationships

Philippine herpetology (Amphibia, Reptilia), 20 years on: two decades of progress towards an increasingly collaborative, equitable, and inclusive approach to the study of the archipelago's amphibians and reptiles

A first review of the history, status, and prospects for Philippine herpetology conducted more than two decades ago (2002) summarized the diverse topics studied and highlighted the development and achievements in research up to the year 2000. This study revisits and re-assesses what Philippine herpetology has accomplished, both as a discipline and a community, during the last two decades (2002-2022). A total of 423 herpetological publications was collated, revealing a substantial increase in annual publications, rising from approximately four per year during 2002-2008 to around 28 per year in 2009-2022. Half of the published studies focused on squamate reptiles (lizards 30.5%, snakes 21%) and 28.4% on amphibians, 5.9% on turtles, and 2.6% on crocodiles. The remaining 11.6% of studies focused simultaneously on multiple taxa (i.e., faunal inventories). Diversity and distribution (35.2%) and ecological (26.5%) studies remained popular, while studies on taxonomy (14.9%), phylogenetics and biogeography (11.8%), and conservation (11.6%) all increased. However, geographical gaps persist urging immediate surveys in many understudied regions of the country. Finally, we found a balanced representation between Filipino and foreign first authors (1.0:1.1), yet a substantial gender gap exists between male and female first authors (7.1:1.0). Nonetheless, the steep increase in publications and the diversity of people engaged in Philippine herpetology is a remarkable positive finding compared to the 20 years preceding the last review (1980-2000). Our hope is that the next decades will bring increasingly equitable, internationally collaborative, and broadly inclusive engagement in the study of amphibians and reptiles in the Philippines.

Highly-multiplexed and efficient long-amplicon PacBio and Nanopore sequencing of hundreds of full mitochondrial genomes

BACKGROUND

Mitochondrial genome sequences have become critical to the study of biodiversity. Genome skimming and other short-read based methods are the most common approaches, but they are not well-suited to scale up to multiplexing hundreds of samples. Here, we report on a new approach to sequence hundreds to thousands of complete mitochondrial genomes in parallel using long-amplicon sequencing. We amplified the mitochondrial genome of 677 specimens in two partially overlapping amplicons and implemented an asymmetric PCR-based indexing approach to multiplex 1,159 long amplicons together on a single PacBio SMRT Sequel II cell. We also tested this method on Oxford Nanopore Technologies (ONT) MinION R9.4 to assess if this method could be applied to other long-read technologies. We implemented several optimizations that make this method significantly more efficient than alternative mitochondrial genome sequencing methods.

RESULTS

With the PacBio sequencing data we recovered at least one of the two fragments for 96% of samples (~ 80-90%) with mean coverage ~ 1,500x. The ONT data recovered less than 50% of input fragments likely due to low throughput and the design of the Barcoded Universal Primers which were optimized for PacBio sequencing. We compared a single mitochondrial gene alignment to half and full mitochondrial genomes and found, as expected, increased tree support with longer alignments, though whole mitochondrial genomes were not significantly better than half mitochondrial genomes.

CONCLUSIONS

This method can effectively capture thousands of long amplicons in a single run and be used to build more robust phylogenies quickly and effectively. We provide several recommendations for future users depending on the evolutionary scale of their system. A natural extension of this method is to collect multi-locus datasets consisting of mitochondrial genomes and several long nuclear loci at once.

Generalizing Bayesian phylogenetics to infer shared evolutionary events

Phylogenetic models have long assumed that lineages diverge independently. Processes of diversification that are of interest in biogeography, epidemiology, and genome evolution violate this assumption by affecting multiple evolutionary lineages. To relax the assumption of independent divergences and infer patterns of divergences predicted by such processes, we introduce a way of conceptualizing, modeling, and inferring phylogenetic trees. We apply the approach to genomic data from geckos distributed across the Philippines and find support for patterns of shared divergences predicted by repeated fragmentation of the archipelago by interglacial rises in sea level.

Many processes of biological diversification can simultaneously affect multiple evolutionary lineages. Examples include multiple members of a gene family diverging when a region of a chromosome is duplicated, multiple viral strains diverging at a “super-spreading” event, and a geological event fragmenting whole communities of species. It is difficult to test for patterns of shared divergences predicted by such processes because all phylogenetic methods assume that lineages diverge independently. We introduce a Bayesian phylogenetic approach to relax the assumption of independent, bifurcating divergences by expanding the space of topologies to include trees with shared and multifurcating divergences. This allows us to jointly infer phylogenetic relationships, divergence times, and patterns of divergences predicted by processes of diversification that affect multiple evolutionary lineages simultaneously or lead to more than two descendant lineages. Using simulations, we find that the method accurately infers shared and multifurcating divergence events when they occur and performs as well as current phylogenetic methods when divergences are independent and bifurcating. We apply our approach to genomic data from two genera of geckos from across the Philippines to test if past changes to the islands’ landscape caused bursts of speciation. Unlike previous analyses restricted to only pairs of gecko populations, we find evidence for patterns of shared divergences. By generalizing the space of phylogenetic trees in a way that is independent from the likelihood model, our approach opens many avenues for future research into processes of diversification across the life sciences.

Insights into the phylogenetic position and phylogeography of the monospecific skink-parasite genus

Neoentomelas asatoi Hasegawa, 1989 is a parasitic nematode that infests only the scincid lizard Ateuchosaurus pellopleurus (Hallowell, 1861) that inhabits the forest floor in the Northern and Central Ryukyu Archipelago, Japan. As a member of Rhabdiasidae, the reproductive mode of N. asatoi is characterised by the alternation of the protandrous hermaphroditic mode and gonochoristic mode throughout the life cycle. The intrafamily phylogenetic position and intraspecific diversity of this nematode species were inferred by molecular phylogenetic analyses. The results revealed the phylogenetic distinctiveness of Neoentomelas Hasegawa, 1989 in Rhabdiasidae that supports the unique generic status of Neoentomelas within the family. The intraspecific phylogenetic analyses of N. asatoi revealed a minor concordant phylogenetic pattern with the host and mosaic geographic arrangement of the major clades that was discordant with the host. The analyses and distribution pattern of subclades suggested that this geographic arrangement can be explained by at least three dispersal events and subsequent switching to indigenous host populations. Colonisation events might be promoted by the high establishment rate of new populations stemming from the parthenogenesis-like reproduction mode of N. asatoi. This study demonstrated that reproductive modes can affect the intraspecific genetic diversity of parasites.

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.

De novo Assembly, Annotation, and Analysis of Transcriptome Data of the Ladakh Ground Skink Provide Genetic Information on High-Altitude Adaptation

The Himalayan Arc is recognized as a global biodiversity hotspot. Among its numerous cryptic and undiscovered organisms, this composite high-mountain ecosystem harbors many taxa with adaptations to life in high elevations. However, evolutionary patterns and genomic features have been relatively rarely studied in Himalayan vertebrates. Here, we provide the first well-annotated transcriptome of a Greater Himalayan reptile species, the Ladakh Ground skink Asymblepharus ladacensis (Squamata: Scincidae). Based on tissues from the brain, an embryonic disc, and pooled organ material, using pair-end Illumina NextSeq 500 RNAseq, we assembled ~77,000 transcripts, which were annotated using seven functional databases. We tested ~1600 genes, known to be under positive selection in anurans and reptiles adapted to high elevations, and potentially detected positive selection for 114 of these genes in Asymblepharus. Even though the strength of these results is limited due to the single-animal approach, our transcriptome resource may be valuable data for further studies on squamate reptile evolution in the Himalayas as a hotspot of biodiversity.

Revision of the montane New Guinean skink genus Lobulia (Squamata: Scincidae), with the description of four new genera and nine new species

The skink genus Lobulia is endemic to New Guinea, the largest and highest tropical island in the world. Lobulia and its related genera represent an important component of the montane herpetofauna of New Guinea, but it remains understudied and poorly known. We here provide the first, large-scale, systematic revision of Lobulia, using molecular phylogenetic and morphological comparisons to assess the monophyly of the genus and the diversity of species within it. We find that Lobulia, as currently defined, is polyphyletic. The eight species currently assigned to it form three clades. Furthermore, many specimens from New Guinea of unknown specific affinity are genetically and morphologically distinct from each other. Based on these data, we re-diagnose Lobulia and two of its closely related genera, Prasinohaema and Papuascincus. We erect four new genera (Alpinoscincus gen. nov., Nubeoscincus gen. nov., Ornithuroscincus gen. nov. and Palaia gen. nov.) to address the problem of polyphyly and describe nine new species Lobulia fortis sp. nov., Lobulia huonensis sp. nov., Loublia marmorata sp. nov., Lobulia vogelkopensis sp. nov., Ornithuroscincus bengaun sp. nov., Ornithuroscincus inornatus sp. nov., Ornithuroscincus pterophilus sp. nov., Ornithuroscincus shearmani sp. nov. and Ornithuroscincus viridis sp. nov. We supplement this taxonomic revision by investigating the biogeographic history of Lobulia s.l. and find evidence for a large radiation in the accreted terranes of New Guinea, with multiple independent colonizations of montane habitats and subsequent recolonization of lowland habitats. Our study reinforces the uniqueness and richness of the montane herpetofauna of New Guinea and the importance of mountains to biodiversity in the Tropics.

Selva Zoque, Mexico: an important Mesoamerican tropical region for reptile species diversity and conservation

The Selva Zoque region is characterized by a great variety of ecosystems for which there is little information about reptile species diversity and their conservation status. This study is the first assessment of the species richness, composition, and conservation status of reptiles of this region. Additionally, this information is compared with that of seven other tropical regions in northern Mesoamerica. In total, 141 native reptile species belonging to 81 genera and 29 families are recorded for the Selva Zoque region. Sixty species (42% of the total) recorded in Selva Zoque are in high-risk categories according to the Mexican Ministry of the Environment, the highest number for the Mexican regions of Mesoamerica. According to the IUCN, six species are in high-risk categories, seven species are in Data Deficient, and 23 (16%) have not been evaluated yet. According to the Environmental Vulnerability Scores approach, 28 species (20%) are in the high vulnerability category. The Selva Zoque species composition is most similar to Los Tuxtlas and Lacandona regions, and most dissimilar to Sian Ka´an Biosphere Reserve. The reptilian fauna of Selva Zoque has a distinctive composition, with the highest number (11 species) of endemic reptiles in the northern Mesoamerican, and species from two biogeographic provinces: the Gulf of Mexico and the Mexican Pacific Coast. These results indicate that the Selva Zoque is the most diverse region in native reptile species in northern Mesoamerica, highlighting it as extremely important for the conservation of the reptile fauna at local (southern Mexico) and regional levels (northern Mesoamerica).

Detection of cryptic diversity in lizards (Squamata) from two Biosphere Reserves in Mesoamerica

A combined approach based on karyology and DNA taxonomy allowed us to characterize the taxonomic peculiarities in 10 Mesoamerican lizard species, belonging to six genera and five families, inhabiting two Biosphere Reserve in Chiapas, Mexico: La Sepultura Biosphere Reserve, and Montes Azules Biosphere. The karyotypes of four species, Phyllodactylus sp. 3 (P. tuberculosus species group) (2n = 38), Holcosus festivus (Lichtenstein et von Martens, 1856) (2n = 50), Anolis lemurinus Cope, 1861 (2n = 40), and A. uniformis Cope, 1885 (2n = 29-30) are described for the first time, the last one showing a particular X₁X₁X₂X₂/X₁X₂Y condition. In Aspidoscelis deppii (Wiegmann, 1834) (2n = 50) and Anolis capito Peters, 1863 (2n = 42), we found a different karyotype from the ones previously reported for these species. Moreover, in A. capito, the cytogenetic observation is concurrent with a considerable genetic divergence (9%) at the studied mtDNA marker (MT-ND2), which is indicative of a putative new cryptic species. The skink Scincella cherriei (Cope, 1893), showed high values of genetic divergence (5.2% at 16S gene) between the specimens from Montes Azules and those from Costa Rica and Nicaragua, comparable to the values typical of sister species in skinks. A lower level of genetic divergence, compatible with an intraspecific phylogeographic structure, has been identified in Lepidophyma flavimaculatum Duméril, 1851. These new data identify taxa that urgently require more in-depth taxonomic studies especially in these areas where habitat alteration is proceeding at an alarming rate.

A taxonomic conundrum: Characterizing a cryptic radiation of Asian gracile skinks (Squamata: Scincidae: Riopa) in Myanmar

Recognizing species-level diversity is important for studying evolutionary patterns across biological disciplines and is critical for conservation efforts. However, challenges remain in delimiting species-level diversity, especially in cryptic radiations where species are genetically divergent but show little morphological differentiation. Using multilocus molecular data, phylogenetic analyses, species delimitation analyses, and morphological data, we examine lineage diversification in a cryptic radiation of Riopa skinks in Myanmar. Four species of Riopa skinks are currently recognized from Myanmar based on morphological traits, but the boundaries between three of these species, R. anguina, R. lineolata, and R. popae, are not well-defined. We find high levels of genetic diversity within these three species. Our analyses suggest that they may comprise as many as 12 independently evolving lineages, highlighting the extent to which species diversity in the region is underestimated. However, quantitative trait data suggest that these lineages have not differentiated morphologically, possibly indicating that this cryptic radiation represents non-adaptive evolution, although additional data is needed to corroborate this.

Cryptic diversity and non-adaptive radiation of montane New Guinea skinks (Papuascincus; Scincidae)

New Guinea, the world's largest and highest tropical island, has a rich but poorly known biota. Papuascincus is a genus of skinks endemic to New Guinea's mountain regions, comprising two wide-ranging species and two species known only from their type series. The phylogeny of the genus has never been examined and the relationships among its species - as well as between it and closely related taxa - are hitherto unknown. We performed the first large-scale molecular-phylogenetic study of Papuascincus, including sampling across the genus' range in Papua New Guinea. We sequenced three mitochondrial and two nuclear markers from 65 specimens of Papuascincus and reconstructed their phylogenetic relationships. We also performed species-delimitation analyses, estimated divergence times and ancestral biogeography, and examined body-size evolution within the genus. Papuascincus was strongly supported as monophyletic. It began radiating during the mid-Miocene in the area now comprising the Central Cordillera of New Guinea, then dispersed eastward colonising the Papuan Peninsula. We found evidence of extensive cryptic diversity within the genus, with between nine and 20 supported genetic lineages. These were estimated using three methods of species delimitation and predominantly occur in allopatry. Distribution and body-size divergence patterns indicated that character displacement in size took place during the evolutionary history of Papuascincus. We conclude that the genus requires comprehensive taxonomic revision and likely represents a species-rich lineage of montane skinks.

Complete mitochondrial genome of Scincella modesta (Squamata: Scincidae)

The first complete mitochondrial genome sequence was determined for an oviparous lizard, Scincella modesta (Scincella, Scincidae). The total length of the complete mitochondrial genome was 17,511 bp, encodes 13 protein-coding genes, 22 tRNAs, 2 rRNA genes, and 2 non-coding regions. The overall base composition of S. modesta is A: 31.9%, T: 27.2%, G: 14.5%, and C: 26.5%. Most of the S. modesta mitochondrial genes are encoded on the H-strand except for the ND6 gene and eight tRNA genes, which are encoded on the L-strand. Mrbayes and ML tree based on 13 protein-coding genes indicated that S. vandenburghi is the sister group of the S. modesta within the genus Scincella. The complete mitogenome sequence of S. modesta provided fundamental data for resolving phylogenetic and genetic problems related to genus Scincella.

Hemipenial morphology of nine South American species of Mabuya (Scincidae: Lygosominae) with comments on the morphology of the family

Species of the genus Mabuya are barely recognizable because of their highly conservative external morphology, resulting in unstable systematic and taxonomic status. Considering that hemipenial morphology has been extensively used as a source of systematic information for several groups of Squamata, we described the hemipenial morphology of nine species of Mabuya, searching for diagnostic features useful to taxonomic and systematic studies. To have comparative parameters within the family, we also presented the description of the hemipenes of two other Scincidae species (i.e., Trachylepis atlantica and Tiliqua rugosa). This information, together with the available literature, was analyzed in a comparative framework, allowing us to define some variation for the genus and family. The hemipenial morphology of the species of Mabuya studied is very conservative; features that show subtle variation were the shape of the protuberances present on the asulcate face of the hemipenial body and the presence and type of a lateral fold, on the body. Also, within the family, some variation was found in the shape and length of lobes, length of branches of the sulcus spermaticus, and in the kind of folds and protuberances in both faces of the lobes and hemipenial body. Although these features were not tested in a systematic context, this information suggested variations in hemipenial morphology that could contribute to taxonomic diagnostics and phylogenetic characters. The importance of obtaining proper preparations for gathering accurate information is emphasized, and the use of fresh-killed specimens that shows fewer difficulties for the procedures of preparation is highly recommended.

The comparative biogeography of Philippine geckos challenges predictions from a paradigm of climate-driven vicariant diversification across an island archipelago

A primary goal of biogeography is to understand how large-scale environmental processes, like climate change, affect diversification. One often-invoked but seldom tested process is the "species-pump" model, in which repeated bouts of cospeciation are driven by oscillating climate-induced habitat connectivity cycles. For example, over the past three million years, the landscape of the Philippine Islands has repeatedly coalesced and fragmented due to sea-level changes associated with glacial cycles. This repeated climate-driven vicariance has been proposed as a model of speciation across evolutionary lineages codistributed throughout the islands. This model predicts speciation times that are temporally clustered around the times when interglacial rises in sea level fragmented the islands. To test this prediction, we collected comparative genomic data from 16 pairs of insular gecko populations. We analyze these data in a full-likelihood, Bayesian model-choice framework to test for shared divergence times among the pairs. Our results provide support against the species-pump model prediction in favor of an alternative interpretation, namely that each pair of gecko populations diverged independently. These results suggest the repeated bouts of climate-driven landscape fragmentation have not been an important mechanism of speciation for gekkonid lizards across the Philippine Archipelago.

Multilocus phylogeny and a new classification for African, Asian and Indian supple and writhing skinks (Scincidae: Lygosominae)

The genera Lepidothyris, Lygosoma and Mochlus comprise the writhing or supple skinks, a group of semi-fossorial, elongate-bodied skinks distributed across the Old World Tropics. Due to their generalized morphology and lack of diagnostic characters, species- and clade-level relationships have long been debated. Recent molecular phylogenetic studies of the group have provided some clarification of species-level relationships, but a number of issues regarding higher level relationships among genera still remain. Here we present a phylogenetic estimate of relationships among species in Lygosoma, Mochlus and Lepidothyris generated by concatenated and species tree analyses of multilocus data using the most extensive taxonomic sampling of the group to date. We also use multivariate statistics to examine species and clade distributions in morpho space. Our results reject the monophyly of Lygosoma s.l., Lygosoma s.s. and Mochlus, which highlights the instability of the current taxonomic classification of the group. We, therefore, revise the taxonomy of the writhing skinks to better reflect the evolutionary history of Lygosoma s.l. by restricting Lygosoma for Southeast Asia, resurrecting the genus Riopa for a clade of Indian and Southeast Asian species, expanding the genus Mochlus to include all African species of writhing skinks and describing a new genus in Southeast Asia.

Amphibians and reptiles of the state of San Luis Potosí, Mexico, with comparisons with adjoining states

A summary of the species of amphibians and reptiles of the state has been compiled, including their geographic distributions, habitats, and conservation statuses. The herpetofauna of San Luis Potosí consists of 41 species of amphibians and 141 species of reptiles. San Luis Potosí shares the highest number of species with Hidalgo and Tamaulipas, and the least number of species with Nuevo León. In San Luis Potosí, there are several taxa of particular conservation concern including salamanders, emydid and trionychid turtles, anguid and xenosaurid lizards, and natricid and colubrid snakes.

The amphibians and reptiles of Mindanao Island, southern Philippines, II: the herpetofauna of northeast Mindanao and adjacent islands

We summarize all available amphibian and reptile species distribution data from the northeast Mindanao faunal region, including small islands associated with this subcenter of endemic vertebrate biodiversity. Together with all publicly available historical information from biodiversity repositories, we present new data from several major herpetological surveys, including recently conducted inventories on four major mountains of northeast Mindanao, and adjacent islands of Camiguin Sur, Dinagat, and Siargao. We present species accounts for all taxa, comment on unresolved taxonomic problems, and provide revisions to outdated IUCN conservation status assessments in cases where our new data significantly alter earlier classification status summaries. Together, our comprehensive analysis of this fauna suggests that the greater Mindanao faunal region possesses distinct subcenters of amphibian and reptile species diversity, and that until this area is revisited and its fauna and actually studied, with on-the-ground field work including targeted surveys of species distributions coupled to the study their natural history, our understanding of the diversity and conservation status of southern Philippine herpetological fauna will remain incomplete. Nevertheless, the northeast Mindanao geographical area (Caraga Region) appears to have the highest herpetological species diversity (at least 126 species) of any comparably-sized Philippine faunal subregion.

Detecting the Anomaly Zone in Species Trees and Evidence for a Misleading Signal in Higher-Level Skink Phylogeny (Squamata: Scincidae)

The anomaly zone, defined by the presence of gene tree topologies that are more probable than the true species tree, presents a major challenge to the accurate resolution of many parts of the Tree of Life. This discrepancy can result from consecutive rapid speciation events in the species tree. Similar to the problem of long-branch attraction, including more data via loci concatenation will only reinforce the support for the incorrect species tree. Empirical phylogenetic studies often employ coalescent-based species tree methods to avoid the anomaly zone, but to this point these studies have not had a method for providing any direct evidence that the species tree is actually in the anomaly zone. In this study, we use 16 species of lizards in the family Scincidae to investigate whether nodes that are difficult to resolve place the species tree within the anomaly zone. We analyze new phylogenomic data (429 loci), using both concatenation and coalescent-based species tree estimation, to locate conflicting topological signal. We then use the unifying principle of the anomaly zone, together with estimates of ancestral population sizes and species persistence times, to determine whether the observed phylogenetic conflict is a result of the anomaly zone. We identify at least three regions of the Scincidae phylogeny that provide demographic signatures consistent with the anomaly zone, and this new information helps reconcile the phylogenetic conflict in previously published studies on these lizards. The anomaly zone presents a real problem in phylogenetics, and our new framework for identifying anomalous relationships will help empiricists leverage their resources appropriately for investigating and overcoming this challenge.

A new species of Eimeria Schneider, 1875 (Apicomplexa: Eimeriidae) from the Solomon ground skink, Sphenomorphus solomonis (Boulenger) (Sauria: Scincidae) from Papua New Guinea

Between September 1990 and November 1991, 19 Sphenomorphus spp. skinks, including nine S. jobiense, three S. simus, and seven Solomon ground skinks, S. solomonis (Boulenger), were collected from Madang and Morobe Provinces, Papua New Guinea (PNG), and examined for coccidia. A single S. solomonis was found to be infected with a new species of Eimeria Schneider, 1875. Oöcysts of Eimeria perkinsae n. sp. are ellipsoidal with a smooth, colourless, bi-layered wall, measure 18.6 × 14.7 μm, and have a length/width (L/W) ratio of 1.3; both micropyle and oöcyst residuum are absent, but a fragmented polar granule is present. Sporocysts are ovoidal, 8.9 × 6.4 μm, L/W 1.4; neither Stieda, sub-Stieda or para-Stieda bodies are present; a sporocyst residuum consisted of a loose cluster of granules dispersed between sporozoites. Sporozoites are comma-shaped with spheroidal anterior and posterior refractile bodies. This represents the first report of coccidia from this skink genus.

Phylogeny of Trachylepis sp. (Reptilia) from Turkey inferred from mtDNA sequences

The taxonomic status of the species included into the genus Trachylepis in Turkey are doubtful. So far, three morphological species have been attributed to this genus in Turkey; Trachylepis aurata. T. vittata, and T. septemtaeniata. Here, we investigated the taxonomy of the Turkish Trachylepis species by employing phylogenetic and phylogeographic approaches and using mitochondrial DNA (cytochrome b and 12 S rRNA). In total, 45 Trachylepis and 6 Mabuya specimens were used analyzed. Phylogenetic analyses were carried out using Maximum Likelihood (ML) and Bayesian Inference (BI) methods. The phylogenetic relationships and the genetic distances retrieved, revealed that the Turkish species, which currently recognized as Trachylepis, are highly diversified, forming a distinct clade that shows closer phylogenetic affinity with the species of the genus Mabuya rather than the other Trachylepis species. In this clade, the three Turkish species are monophyletic with T. vittata to branch off first in late Miocene (10.54 Mya). The other two species (T. septemtaeniata and T. aurata) seem to have sister group relationship that diverged at the end of Messinian Salinity Crisis (5.27 Mya). As a whole, the examination of mtDNA lineages in the Turkish lizards of the genus Trachylepis may contribute substantially to the refining of their taxonomic status, since the three species of Turkey, although monophyletic, represent a distinct radiation that would could probably recognized as a different genus in Mabuya sensu lato.

A phylogeny and revised classification of Squamata, including 4161 species of lizards and snakes

BACKGROUND

The extant squamates (>9400 known species of lizards and snakes) are one of the most diverse and conspicuous radiations of terrestrial vertebrates, but no studies have attempted to reconstruct a phylogeny for the group with large-scale taxon sampling. Such an estimate is invaluable for comparative evolutionary studies, and to address their classification. Here, we present the first large-scale phylogenetic estimate for Squamata.

RESULTS

The estimated phylogeny contains 4161 species, representing all currently recognized families and subfamilies. The analysis is based on up to 12896 base pairs of sequence data per species (average = 2497 bp) from 12 genes, including seven nuclear loci (BDNF, c-mos, NT3, PDC, R35, RAG-1, and RAG-2), and five mitochondrial genes (12S, 16S, cytochrome b, ND2, and ND4). The tree provides important confirmation for recent estimates of higher-level squamate phylogeny based on molecular data (but with more limited taxon sampling), estimates that are very different from previous morphology-based hypotheses. The tree also includes many relationships that differ from previous molecular estimates and many that differ from traditional taxonomy.

CONCLUSIONS

We present a new large-scale phylogeny of squamate reptiles that should be a valuable resource for future comparative studies. We also present a revised classification of squamates at the family and subfamily level to bring the taxonomy more in line with the new phylogenetic hypothesis. This classification includes new, resurrected, and modified subfamilies within gymnophthalmid and scincid lizards, and boid, colubrid, and lamprophiid snakes.

The amphibians and reptiles of Luzon Island, Philippines, VIII: the herpetofauna of Cagayan and Isabela Provinces, northern Sierra Madre Mountain Range

We provide the first report on the herpetological biodiversity (amphibians and reptiles) of the northern Sierra Madre Mountain Range (Cagayan and Isabela provinces), northeast Luzon Island, Philippines. New data from extensive previously unpublished surveys in the Municipalities of Gonzaga, Gattaran, Lasam, Santa Ana, and Baggao (Cagayan Province), as well as fieldwork in the Municipalities of Cabagan, San Mariano, and Palanan (Isabela Province), combined with all available historical museum records, suggest this region is quite diverse. Our new data indicate that at least 101 species are present (29 amphibians, 30 lizards, 35 snakes, two freshwater turtles, three marine turtles, and two crocodilians) and now represented with well-documented records and/or voucher specimens, confirmed in institutional biodiversity repositories. A high percentage of Philippine endemic species constitute the local fauna (approximately 70%). The results of this and other recent studies signify that the herpetological diversity of the northern Philippines is far more diverse than previously imagined. Thirty-eight percent of our recorded species are associated with unresolved taxonomic issues (suspected new species or species complexes in need of taxonomic partitioning). This suggests that despite past and present efforts to comprehensively characterize the fauna, the herpetological biodiversity of the northern Philippines is still substantially underestimated and warranting of further study.

Evidence for climate-driven diversification? A caution for interpreting ABC inferences of simultaneous historical events

Approximate Bayesian computation (ABC) is rapidly gaining popularity in population genetics. One example, msBayes, infers the distribution of divergence times among pairs of taxa, allowing phylogeographers to test hypotheses about historical causes of diversification in co-distributed groups of organisms. Using msBayes, we infer the distribution of divergence times among 22 pairs of populations of vertebrates distributed across the Philippine Archipelago. Our objective was to test whether sea-level oscillations during the Pleistocene caused diversification across the islands. To guide interpretation of our results, we perform a suite of simulation-based power analyses. Our empirical results strongly support a recent simultaneous divergence event for all 22 taxon pairs, consistent with the prediction of the Pleistocene-driven diversification hypothesis. However, our empirical estimates are sensitive to changes in prior distributions, and our simulations reveal low power of the method to detect random variation in divergence times and bias toward supporting clustered divergences. Our results demonstrate that analyses exploring power and prior sensitivity should accompany ABC model selection inferences. The problems we identify are potentially mitigable with uniform priors over divergence models (rather than classes of models) and more flexible prior distributions on demographic and divergence-time parameters.