Phylogeny of Philippine slender skinks (Scincidae: Brachymeles) reveals underestimated species diversity, complex biogeographical relationships, and cryptic patterns of lineage diversification

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.

Hidden limbs in the "limbless skink" Brachymeles lukbani: Developmental observations

Reduced limbs and limblessness have evolved independently in many lizard clades. Scincidae exhibit a wide range of limb-reduced morphologies, but only some species have been used to study the embryology of limb reduction (e.g., digit reduction in Chalcides and limb reduction in Scelotes). The genus Brachymeles, a Southeast Asian clade of skinks, includes species with a range of limb morphologies, from pentadactyl to functionally and structurally limbless species. Adults of the small, snake-like species Brachymeles lukbani show no sign of external limbs in the adult except for small depressions where they might be expected to occur. Here, we show that embryos of B. lukbani in early stages of development, on the other hand, show a truncated but well-developed limb with a stylopod and a zeugopod, but no signs of an autopod. As development proceeds, the limb's small size persists even while the embryo elongates. These observations are made based on external morphology. We used florescent whole-mount immunofluorescence to visualize the morphology of skeletal elements and muscles within the embryonic limb of B. lukabni. Early stages have a humerus and separated ulna and radius cartilages; associated with these structures are dorsal and ventral muscle masses as those found in the embryos of other limbed species. While the limb remains small, the pectoral girdle grows in proportion to the rest of the body, with well-developed skeletal elements and their associated muscles. In later stages of development, we find the small limb is still present under the skin, but there are few indications of its presence, save for the morphology of the scale covering it. By use of CT scanning, we find that the adult morphology consists of a well-developed pectoral girdle, small humerus, extremely reduced ulna and radius, and well-developed limb musculature connected to the pectoral girdle. These muscles form in association with a developing limb during embryonic stages, a hint that "limbless" lizards that possess these muscles may have or have had at least transient developing limbs, as we find in B. lukbani. Overall, this newly observed pattern of ontogenetic reduction leads to an externally limbless adult in which a limb rudiment is hidden and covered under the trunk skin, a situation called cryptomelia. The results of this work add to our growing understanding of clade-specific patterns of limb reduction and the convergent evolution of limbless phenotypes through different developmental processes.

Gene Flow Increases Phylogenetic Structure and Inflates Cryptic Species Estimations: A Case Study on Widespread Philippine Puddle Frogs (Occidozyga laevis)

In cryptic amphibian complexes, there is a growing trend to equate high levels of genetic structure with hidden cryptic species diversity. Typically, phylogenetic structure and distance-based approaches are used to demonstrate the distinctness of clades and justify the recognition of new cryptic species. However, this approach does not account for gene flow, spatial, and environmental processes that can obfuscate phylogenetic inference and bias species delimitation. As a case study, we sequenced genome-wide exons and introns to evince the processes that underlie the diversification of Philippine Puddle Frogs-a group that is widespread, phenotypically conserved, and exhibits high levels of geographically-based genetic structure. We showed that widely adopted tree- and distance-based approaches inferred up to 20 species, compared to genomic analyses that inferred an optimal number of five distinct genetic groups. Using a suite of clustering, admixture, and phylogenetic network analyses, we demonstrate extensive admixture among the five groups and elucidate two specific ways in which gene flow can cause overestimations of species diversity: (1) admixed populations can be inferred as distinct lineages characterized by long branches in phylograms; and (2) admixed lineages can appear to be genetically divergent, even from their parental populations when simple measures of genetic distance are used. We demonstrate that the relationship between mitochondrial and genome-wide nuclear p-distances is decoupled in admixed clades, leading to erroneous estimates of genetic distances and, consequently, species diversity. Additionally, genetic distance was also biased by spatial and environmental processes. Overall, we showed that high levels of genetic diversity in Philippine Puddle Frogs predominantly comprise metapopulation lineages that arose through complex patterns of admixture, isolation-by-distance, and isolation-by-environment as opposed to species divergence. Our findings suggest that speciation may not be the major process underlying the high levels of hidden diversity observed in many taxonomic groups and that widely-adopted tree- and distance-based methods overestimate species diversity in the presence of gene flow.

Locomotion and palaeoclimate explain the re-evolution of quadrupedal body form in Brachymeles lizards

Evolutionary reversals, including re-evolution of lost structures, are commonly found in phylogenetic studies. However, we lack an understanding of how these reversals happen mechanistically. A snake-like body form has evolved many times in vertebrates, and occasionally a quadrupedal form has re-evolved, including in Brachymeles lizards. We use body form and locomotion data for species ranging from snake-like to quadrupedal to address how a quadrupedal form could re-evolve. We show that large, quadrupedal species are faster at burying and surface locomotion than snake-like species, indicating a lack of expected performance trade-off between these modes of locomotion. Species with limbs use them while burying, suggesting that limbs are useful for burying in wet, packed substrates. Palaeoclimatological data suggest that Brachymeles originally evolved a snake-like form under a drier climate probably with looser soil in which it was easier to dig. The quadrupedal clade evolved as the climate became humid, where limbs and large size facilitated fossorial locomotion in packed soils.

Using 3D-digital photogrammetry to examine scaling of the body axis in burrowing skinks

Three-dimensional (3D) modeling techniques have been increasingly utilized across disciplines for the visualization and analysis of complex structures. We employ 3D-digital photogrammetry for understanding the scaling of the body axis of 12 species of scincid lizards in the genus Brachymeles. These skinks represent a diverse radiation which shows tremendous variation in body size and degree of axial elongation. Because of the complex nature of the body axis, 3D-methods are important for understanding how the body axis evolves. 3D-digital photogrammetry presents a flexible, inexpensive, and portable system for the reconstruction of biological forms. As body size increased among species, the cross-sectional area and circumference of the head and other portions of the body axis increased isometrically, which indicates that species of differing sizes possess proportionally similar head and body shapes. These results suggest that there are no substantial head and body shape changes with body size among the sampled species, but further comparative studies with larger sample sizes and functional studies of size and morphology effects on burrowing or above-ground locomotion are needed.

Phylogenetic relationships and biogeographic range evolution in cat-eyed snakes, Boiga (Serpentes: Colubridae)

The genus Boiga includes 35, primarily arboreal snake species distributed from the Middle East to Australia and many islands in the western Pacific, with particularly high species diversity in South-East Asia. Despite including the iconic mangrove snakes (Boiga dendrophila complex) and the brown tree snake (Boiga irregularis; infamous for avian extinctions on small islands of the Pacific), species-level phylogenetic relationships and the biogeographic history of this ecologically and morphologically distinct clade are poorly understood. In this study, we sequenced mitochondrial and nuclear DNA for 24 Boiga species and used these data to estimate a robust phylogenetic inference, in order to (1) test the hypothesis that Boiga is monophyletic, (2) evaluate the validity of current species-level taxonomy and (3) examine whether geographic range evolution in Boiga is consistent with expectations concerning dispersal and colonization of vertebrates between continents and islands. Our results support the prevailing view that most dispersal events are downstream – from continents to oceanic islands – but we also identify a role for upstream dispersal from oceanic islands to continents. Additionally, the novel phylogeny of Boiga presented here is informative for updating species-level taxonomy within the genus.

Convergent Evolution of Elongate Forms in Craniates and of Locomotion in Elongate Squamate Reptiles

Synopsis Elongate, snake- or eel-like, body forms have evolved convergently many times in most major lineages of vertebrates. Despite studies of various clades with elongate species, we still lack an understanding of their evolutionary dynamics and distribution on the vertebrate tree of life. We also do not know whether this convergence in body form coincides with convergence at other biological levels. Here, we present the first craniate-wide analysis of how many times elongate body forms have evolved, as well as rates of its evolution and reversion to a non-elongate form. We then focus on five convergently elongate squamate species and test if they converged in vertebral number and shape, as well as their locomotor performance and kinematics. We compared each elongate species to closely related quadrupedal species and determined whether the direction of vertebral or locomotor change matched in each case. The five lineages examined are obscure species from remote locations, providing a valuable glimpse into their biology. They are the skink lizards Brachymeles lukbani, Lerista praepedita, and Isopachys anguinoides, the basal squamate Dibamus novaeguineae, and the basal snake Malayotyphlops cf. ruficaudus. Our results support convergence among these species in the number of trunk and caudal vertebrae, but not vertebral shape. We also find that the elongate species are relatively slower than their limbed counterparts and move with lower frequency and higher amplitude body undulations, with the exception of Isopachys. This is among the first evidence of locomotor convergence across distantly related, elongate species.

Multilocus phylogeny of Bornean Bent-Toed geckos (Gekkonidae: Cyrtodactylus) reveals hidden diversity, taxonomic disarray, and novel biogeographic patterns

The gekkonid genus Cyrtodactylus is a highly diverse group of lizards (280 + species), which covers an expansive geographic range. Although this genus has been the focus of many taxonomic and molecular systematic studies, species on the Southeast Asian island of Borneo have remained understudied, leading to an unclear evolutionary history with cascading effects on taxonomy and biogeographic inferences. We assembled the most comprehensive multilocus Bornean dataset (one mitochondrial and three nuclear loci) that included 129 novel sequences and representatives from each known Cyrtodactylus species on the island to validate taxonomic status, assess species diversity, and elucidate biogeographic patterns. Our results uncovered a high proportion of cryptic diversity and revealed numerous taxonomic complications, especially within the C. consobrinus, C. malayanus, and C. pubisulcus groups. Comparisons of pairwise genetic distances and a preliminary species delimitation analysis using the Automatic Barcode Gap Discovery (ABGD) method demonstrated that some wide-ranging species on Borneo likely comprise multiple distinct and deeply divergent lineages, each with more restricted distributional ranges. We also tested the prevailing biogeographic hypothesis of a single invasion from Borneo into the Philippines. Our analyses revealed that Philippine taxa were not monophyletic, but were likely derived from multiple separate invasions into the geopolitical areas comprising the Philippines. Although our investigation of Bornean Cyrtodactylus is the most comprehensive to-date, it highlights the need for expanded taxonomic sampling and suggests that our knowledge of the evolutionary history, systematics, and biogeography of Bornean Cyrtodactylus is far from complete.

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.

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.

Colonization and diversification of the white-browed shortwing (Aves: Muscicapidae: Brachypteryx montana) in the Philippines

Molecular phylogenetic approaches have greatly improved our knowledge of the pattern and process of biological diversification across the globe; however, many regions remain poorly documented, even for well-studied vertebrate taxa. The Philippine archipelago, one of the least-studied 'biodiversity hotspots', is an ideal natural laboratory for investigating the factors driving diversification in an insular and geologically dynamic setting. We investigated the history and geography of diversification of the Philippine populations of a widespread montane bird, the White-browed Shortwing (Brachypteryx montana). Leveraging dense archipelago-wide sampling, we generated a multi-locus genetic dataset (one nuclear and two mtDNA markers), which we analyzed using phylogenetic, population genetic, and coalescent-based methods. Our results demonstrate that Philippine shortwings (1) likely colonized the Philippines from the Sunda Shelf to Mindanao in the late Miocene or Pliocene, (2) diversified across inter-island barriers into three divergent lineages during the Pliocene and early Pleistocene, (3) have not diversified within the largest island, Luzon, contrary to patterns observed in other montane taxa, and (4) colonized Palawan from the oceanic Philippines rather than from Borneo, challenging the assumption of Palawan functioning exclusively as a biogeographic extension of the Sunda Shelf. Additionally, our finding that divergent (c. 4.0 mya) lineages are coexisting in secondary sympatry on Mindanao without apparent gene flow suggests that the speciation process is likely complete for these shortwing lineages. Overall, these investigations provide insight into how topography and island boundaries influence diversification within remote oceanic archipelagos and echo the results of many other studies in demonstrating that taxonomic diversity continues to be underestimated in the Philippines.

Ixora (Rubiaceae) on the Philippines - crossroad or cradle?

BACKGROUND

The Philippine archipelago is globally one of the most important model island systems for studying evolutionary processes. However, most plant species on this archipelago have not yet been studied in sufficient detail. The main aim of this study is to unravel the evolutionary history and biogeographic relationships of the Philippine members of the pantropical genus Ixora.

RESULTS

The complex plastid and nuclear divergence patterns in Philippine Ixora, documented using tree and network approaches, reveal a highly dynamic evolution in Ixora, involving several phases of radiation and recolonization. Philippine Ixora comprises at least five lineages, of which one is most closely related to species from Wallacea, and the remaining four to species from Asia.

CONCLUSIONS

Our study highlights the importance of Philippine species for understanding phytogeographic patterns in the Indomalayan-Australasian eco-region. The overall genetic differentiation, as well as the incongruence between genealogies based on the biparentally inherited nucleome and the maternally inherited plastome in Ixora, reflect the complex tectonic history of the Philippine archipelago. The Ixora lineage related to Wallacean species supports the delimitation of different ecozones along Huxley's line, because it is absent from Palawan. The remaining four lineages are all allied with Asian taxa, reflecting several waves of colonization. Close relationships between some widespread Philippine species and locally adapted narrow endemics suggest that the widespread, genetically diverse species act as pools for the formation of new species in a process of ongoing speciation. Our results suggest that the species concepts of some of the more widespread taxa need to be revised.

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.

Molecular phylogeny of Panaspis and Afroablepharus skinks (Squamata: Scincidae) in the savannas of sub-Saharan Africa

African snake-eyed skinks are relatively small lizards of the genera Panaspis and Afroablepharus. Species allocation of these genera frequently changed during the 20th century based on morphology, ecology, and biogeography. Members of these genera occur primarily in savanna habitats throughout sub-Saharan Africa and include species whose highly conserved morphology poses challenges for taxonomic studies. We sequenced two mitochondrial (16S and cyt b) and two nuclear genes (PDC and RAG1) from 76 Panaspis and Afroablepharus samples from across eastern, central, and southern Africa. Concatenated gene-tree and divergence-dating analyses were conducted to infer phylogenies and biogeographic patterns. Molecular data sets revealed several cryptic lineages, with most radiations occurring during the mid-Miocene to Pliocene. We infer that rifting processes (including the formation of the East African Rift System) and climatic oscillations contributed to the expansion and contraction of savannas, and caused cladogenesis in snake-eyed skinks. Species in Panaspis and Afroablepharus used in this study, including type species for both genera, formed a monophyletic group. As a result, the latter genus should be synonymized with the former, which has priority. Conservatively, we continue to include the West African species P. breviceps and P. togoensis within an expanded Panaspis, but note that they occur in relatively divergent clades, and their taxonomic status may change with improved taxon sampling. Divergence estimates and cryptic speciation patterns of snake-eyed skinks were consistent with previous studies of other savanna vertebrate lineages from the same areas examined in this study.

Macroevolution of hyperdiverse flightless beetles reflects the complex geological history of the Sunda Arc

The Sunda Arc forms an almost continuous chain of islands and thus a potential dispersal corridor between mainland Southeast Asia and Melanesia. However, the Sunda Islands have rather different geological histories, which might have had an important impact on actual dispersal routes and community assembly. Here, we reveal the biogeographical history of hyperdiverse and flightless Trigonopterus weevils. Different approaches to ancestral area reconstruction suggest a complex east to west range expansion. Out of New Guinea, Trigonopterus repeatedly reached the Moluccas and Sulawesi transgressing Lydekker's Line. Sulawesi repeatedly acted as colonization hub for different segments of the Sunda Arc. West Java, East Java and Bali are recognized as distinct biogeographic areas. The timing and diversification of species largely coincides with the geological chronology of island emergence. Colonization was not inhibited by traditional biogeographical boundaries such as Wallace's Line. Rather, colonization patterns support distance dependent dispersal and island age limiting dispersal.

Appendicular skeletons of five Asian skink species of the genera Brachymeles and Ophiomorus, including species with vestigial appendicular structures

Vestigial biological structures provide an important line of evidence for macroevolution. They abound in the appendicular skeletons of limbless and reduced-limbed members of the skink subfamily Scincinae, which includes a predominantly Asian clade and a predominantly African clade. Reduced appendicular skeletons in the predominantly African clade have received much recent attention, but for most species in the predominantly Asian clade the appendicular skeleton has yet to be described. Here we provide descriptions of the appendicular skeletons of the reduced-limbed skinks Brachymeles bonitae and Ophiomorus blandfordi, the externally limbless skink Ophiomorus punctatissimus, and, for comparison, the pentadactyl skinks Brachymeles gracilis and B. talinis. We used x-ray radiographs to examine the skeletons of these species and to note similarities and differences in the previously-described appendicular skeletal morphology of related species. We found that in B. bonitae the pectoral and pelvic girdles are unreduced, the proximal limb elements are reduced, and the distal limb elements are vestigial. In O. punctatissimus vestigial pectoral and pelvic girdles are present. In O. blanfordi the fifth metatarsal is vestigial. The phylogenetic distribution of morphological features related to appendicular reduction shows that multiple, parallel reduction events have taken place within each of these two genera. In addition, the anatomical distribution of element reduction and loss in these genera shows that the bones are reduced and lost in the same sequence in the predominantly Asian scincine clade as they are in other squamate clades. This suggests a common evolutionary mechanism for appendicular reduction and loss across the Squamata.

Hand/foot splitting and the 're-evolution' of mesopodial skeletal elements during the evolution and radiation of chameleons

Background

One of the most distinctive traits found within Chamaeleonidae is their split/cleft autopodia and the simplified and divergent morphology of the mesopodial skeleton. These anatomical characteristics have facilitated the adaptive radiation of chameleons to arboreal niches. To better understand the homology of chameleon carpal and tarsal elements, the process of syndactyly, cleft formation, and how modification of the mesopodial skeleton has played a role in the evolution and diversification of chameleons, we have studied the Veiled Chameleon (Chamaeleo calyptratus). We analysed limb patterning and morphogenesis through in situ hybridization, in vitro whole embryo culture and pharmacological perturbation, scoring for apoptosis, clefting, and skeletogenesis. Furthermore, we framed our data within a phylogenetic context by performing comparative skeletal analyses in 8 of the 12 currently recognized genera of extant chameleons.

Results

Our study uncovered a previously underappreciated degree of mesopodial skeletal diversity in chameleons. Phylogenetically derived chameleons exhibit a ‘typical’ outgroup complement of mesopodial elements (with the exception of centralia), with twice the number of currently recognized carpal and tarsal elements considered for this clade. In contrast to avians and rodents, mesenchymal clefting in chameleons commences in spite of the maintenance of a robust apical ectodermal ridge (AER). Furthermore, Bmp signaling appears to be important for cleft initiation but not for maintenance of apoptosis. Interdigital cell death therefore may be an ancestral characteristic of the autopodium, however syndactyly is an evolutionary novelty. In addition, we find that the pisiform segments from the ulnare and that chameleons lack an astragalus-calcaneum complex typical of amniotes and have evolved an ankle architecture convergent with amphibians in phylogenetically higher chameleons.

Conclusion

Our data underscores the importance of comparative and phylogenetic approaches when studying development. Body size may have played a role in the characteristic mesopodial skeletal architecture of chameleons by constraining deployment of the skeletogenic program in the smaller and earliest diverged and basal taxa. Our study challenges the ‘re-evolution’ of osteological features by showing that ‘re-evolving’ a ‘lost’ feature de novo (contrary to Dollo’s Law) may instead be due to so called ‘missing structures’ being present but underdeveloped and/or fused to other adjacent elements (cryptic features) whose independence may be re-established under changes in adaptive selective pressure.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0464-4) contains supplementary material, which is available to authorized users.

Distinct patterns of desynchronized limb regression in malagasy scincine lizards (squamata, scincidae)

Scincine lizards in Madagascar form an endemic clade of about 60 species exhibiting a variety of ecomorphological adaptations. Several subclades have adapted to burrowing and convergently regressed their limbs and eyes, resulting in a variety of partial and completely limbless morphologies among extant taxa. However, patterns of limb regression in these taxa have not been studied in detail. Here we fill this gap in knowledge by providing a phylogenetic analysis of DNA sequences of three mitochondrial and four nuclear gene fragments in an extended sampling of Malagasy skinks, and microtomographic analyses of osteology of various burrowing taxa adapted to sand substrate. Based on our data we propose to (i) consider Sirenoscincus Sakata & Hikida, 2003, as junior synonym of Voeltzkowia Boettger, 1893; (ii) resurrect the genus name Grandidierina Mocquard, 1894, for four species previously included in Voeltzkowia; and (iii) consider Androngo Brygoo, 1982, as junior synonym of Pygomeles Grandidier, 1867. By supporting the clade consisting of the limbless Voeltzkowia mira and the forelimb-only taxa V. mobydick and V. yamagishii, our data indicate that full regression of limbs and eyes occurred in parallel twice in the genus Voeltzkowia (as hitherto defined) that we consider as a sand-swimming ecomorph: in the Voeltzkowia clade sensu stricto the regression first affected the hindlimbs and subsequently the forelimbs, whereas the Grandidierina clade first regressed the forelimbs and subsequently the hindlimbs following the pattern prevalent in squamates. Timetree reconstructions for the Malagasy Scincidae contain a substantial amount of uncertainty due to the absence of suitable primary fossil calibrations. However, our preliminary reconstructions suggest rapid limb regression in Malagasy scincids with an estimated maximal duration of 6 MYr for a complete regression in Paracontias, and 4 and 8 MYr respectively for complete regression of forelimbs in Grandidierina and hindlimbs in Voeltzkowia.

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.