Exploring Paleogene Tibet's warm temperate environments through target enrichment and phylogenetic niche modelling of Himalayan spiny frogs (Paini, Dicroglossidae).

The Cenozoic topographic development of the Himalaya-Tibet orogen (HTO) substantially affected the paleoenvironment and biodiversity patterns of High Asia. However, concepts on the evolution and paleoenvironmental history of the HTO differ massively in timing, elevational increase and sequence of surface uplift of the different elements of the orogen. Using target enrichment of a large set of transcriptome-derived markers, ancestral range estimation and paleoclimatic niche modelling, we assess a recently proposed concept of a warm temperate paleo-Tibet in Asian spiny frogs of the tribe Paini and reconstruct their historical biogeography. That concept was previously developed in invertebrates. Because of their early evolutionary origin, low dispersal capacity, high degree of local endemism, and strict dependence on temperature and humidity, the cladogenesis of spiny frogs may echo the evolution of the HTO paleoenvironment. We show that diversification of main lineages occurred during the early to Mid-Miocene, while the evolution of alpine taxa started during the late Miocene/early Pliocene. Our distribution and niche modelling results indicate range shifts and niche stability that may explain the modern disjunct distributions of spiny frogs. They probably maintained their (sub)tropical or (warm)temperate preferences and moved out of the ancestral paleo-Tibetan area into the Himalaya as the climate shifted, as opposed to adapting in situ. Based on ancestral range estimation, we assume the existence of low-elevation, climatically suitable corridors across paleo-Tibet during the Miocene along the Kunlun, Qiangtang and/or Gangdese Shan. Our results contribute to a deeper understanding of the mechanisms and processes of faunal evolution in the HTO.

Two new endemic species of Ameiva (Squamata: Teiidae) from the dry forest of northwestern Peru and additional information on Ameiva concolor Ruthven, 1924.

We describe two new species of Ameiva Meyer, 1795 from the dry forest of the Northern Peruvian Andes. The new species Ameiva nodam sp. nov. and Ameiva aggerecusans sp. nov. share a divided frontal plate and are differentiated from each other and from their congeners based on genetic (12S and 16S rRNA genes) and morphological characteristics. A. nodam sp. nov. has dilated postbrachials, a maximum known snout-vent length of 101 mm, 10 longitudinal rows of ventral plates, 86-113 midbody granules, 25-35 lamellae under the fourth toe, and a color pattern with 5 longitudinal yellow stripes on the dorsum. Ameiva aggerecusans sp. nov. has not or only hardly dilated postbrachials, a maximum known snout-vent length of 99.3 mm, 10-12 longitudinal rows of ventral plates, 73-92 midbody granules, 31-39 lamellae under the fourth toe, and the females and juveniles of the species normally exhibit a cream-colored vertebral stripe on a dark dorsum ground color. We provide information on the intraspecific variation and distribution of A. concolor. Furthermore, we provide information on the environmental niches of the taxa and test for niche conservatism. 

Phylogenetic and morphological influence on habitat choice in moisture-harvesting horned lizards (Phrynosoma spp.).

In previous studies, the superhydrophilic skin of moisture-harvesting lizards has been linked to the morphological traits of the lizards' integument, that is, the occurrence of honeycomb-shaped microstructures. Interestingly, these structures can also cover the skin of lizards inhabiting wet habitats. We therefore tested the influence of the microstructures' main features on the habitat choice and wettability in the genus Phrynosoma. The genus Phrynosoma comprises moisture-harvesting species as well as nonspecialists. Lizards of this genus inhabit large areas of North America with diverse climatic conditions. Remarkably, the differences in the manifestation of microstructures are just as versatile as their surroundings. The phylogeny of the lizards as well as the depth of their ventral microstructures, though independent of each other, correlated with the precipitation in their respective habitat. All other morphological traits, as well as the skin's wettability itself, could not predict the habitat of Phrynosoma species. Hence, it is unlikely that the microstructure influences the wettability, at least directly. Hence, we presume an indirect influence for the following reasons: (a) As the ventral side cannot get wet by rain, but the belly could easily interact with a wet surface, the microstructure might facilitate water absorption from wet soil following precipitation. (b) We found the number of dorsal microstructures to be linked to the occurrence of silt in the habitat. In our study, we observed scales being heavily contaminated, most likely with a mixture of dead skin (after shedding) and silt. As many lizards burrow themselves or even shovel sand onto their backs, deploying the substrate might be a mechanism to increase the skin's wettability.

Continuous expansion of the geographic range linked to realized niche expansion in the invasive Mourning gecko Lepidodactylus lugubris (Duméril & Bibron, 1836).

Lepidodactylus lugubris is a parthenogenetic gecko which has been increasingly expanding its range during the last century. This invasive species has been reported from multiple tropical and subtropical countries in five continents, most of which were colonized in recent times. In order to understand how the realized niche of the species was affected by this dramatic geographic range expansion, we reconstructed the history of the geographic range expansion. We built models of the realized niche of the species at different points in time during the invasion process. This was achieved through the implementation of modern hypervolume construction methods, based on the Hutchinson's niche concept. The models were then compared to detect possible realized climatic niche expansion over time. Furthermore, we investigated possible pathways used by the species to spread. A progressive expansion of the realized niche was identified. As the species spread into new areas, we observed a tendency to colonize regions with warmer temperatures and higher precipitation rates. Finally, we found evidence for cargo shipping being the major pathway through which the species expands its range. Further studies on this topic should aim to investigate the role of biological interactions, and how they shape the distribution of L. lugubris on a global scale. A deeper understanding of this kind of processes will help us tackle the issue of invasive species, which has become a major challenge in conservation biology.

Applying n-dimensional hypervolumes for species delimitation: unexpected molecular, morphological, and ecological diversity in the Leaf-Toed Gecko Phyllodactylus reissii Peters, 1862 (Squamata: Phyllodactylidae) from northern Peru.

An integrative taxonomic approach based on morphology, molecular analyses, and climatic niche modeling was used to uncover cryptic diversity in the phyllodactylid gecko species Phyllodactylus reissii. At least three distinct species could be identified among the examined specimens from southern Ecuador and northern Peru. Phyllodactylus magister, described by Noble (1924) from arid Andean valleys of the Chinchipe and Marañón rivers in the Peruvian Department of Cajamarca and synonymized with P. reissii by Dixon & Huey (1970) is elevated from synonymy and a detailed redescription is provided. A new species of the genus Phyllodactylus from the Andean dry forest of the southern Marañón valley is identified and described herein. Phyllodactylus pachamama sp. nov. is differentiated from other South American congeners on the basis of mtDNA sequence divergence, morphological characters, and differences in the realized climatic niche. At least in Peru, P. reissii seems to primarily inhabit the northern coastal region west of the Andes, while the inter-Andean area along the Río Marañón and its tributaries seems to be inhabited mostly by other species of the genus, which are endemic to this area. The Andean valleys are underestimated in terms of biodiversity and lack thorough investigation and conservation actions.

Integrative Taxonomy of Southeast Asian Snail-Eating Turtles (Geoemydidae: Malayemys) Reveals a New Species and Mitochondrial Introgression.

Based on an integrative taxonomic approach, we examine the differentiation of Southeast Asian snail-eating turtles using information from 1863 bp of mitochondrial DNA, 12 microsatellite loci, morphology and a correlative species distribution model. Our analyses reveal three genetically distinct groups with limited mitochondrial introgression in one group. All three groups exhibit distinct nuclear gene pools and distinct morphology. Two of these groups correspond to the previously recognized species Malayemys macrocephala (Chao Phraya Basin) and M. subtrijuga (Lower Mekong Basin). The third and genetically most divergent group from the Khorat Basin represents a previously unrecognized species, which is described herein. Although Malayemys are extensively traded and used for religious release, only few studied turtles appear to be translocated by humans. Historic fluctuations in potential distributions were assessed using species distribution models (SDMs). The Last Glacial Maximum (LGM) projection of the predictive SDMs suggests two distinct glacial distribution ranges, implying that the divergence of M. macrocephala and M. subtrijuga occurred in allopatry and was triggered by Pleistocene climate fluctuations. Only the projection derived from the global circulation model MIROC reveals a distinct third glacial distribution range for the newly discovered Malayemys species.

Cryptic speciation patterns in Iranian rock lizards uncovered by integrative taxonomy.

While traditionally species recognition has been based solely on morphological differences either typological or quantitative, several newly developed methods can be used for a more objective and integrative approach on species delimitation. This may be especially relevant when dealing with cryptic species or species complexes, where high overall resemblance between species is coupled with comparatively high morphological variation within populations. Rock lizards, genus Darevskia, are such an example, as many of its members offer few diagnostic morphological features. Herein, we use a combination of genetic, morphological and ecological criteria to delimit cryptic species within two species complexes, D. chlorogaster and D. defilippii, both distributed in northern Iran. Our analyses are based on molecular information from two nuclear and two mitochondrial genes, morphological data (15 morphometric, 16 meristic and four categorical characters) and eleven newly calculated spatial environmental predictors. The phylogeny inferred for Darevskia confirmed monophyly of each species complex, with each of them comprising several highly divergent clades, especially when compared to other congeners. We identified seven candidate species within each complex, of which three and four species were supported by Bayesian species delimitation within D. chlorogaster and D. defilippii, respectively. Trained with genetically determined clades, Ecological Niche Modeling provided additional support for these cryptic species. Especially those within the D. defilippii-complex exhibit well-differentiated niches. Due to overall morphological resemblance, in a first approach PCA with mixed variables only showed the separation between the two complexes. However, MANCOVA and subsequent Discriminant Analysis performed separately for both complexes allowed for distinction of the species when sample size was large enough, namely within the D. chlorogaster-complex. In conclusion, the results support four new species, which are described herein.

Evaluating the significance of paleophylogeographic species distribution models in reconstructing quaternary range-shifts of nearctic chelonians.

The climatic cycles of the Quaternary, during which global mean annual temperatures have regularly changed by 5-10°C, provide a special opportunity for studying the rate, magnitude, and effects of geographic responses to changing climates. During the Quaternary, high- and mid-latitude species were extirpated from regions that were covered by ice or otherwise became unsuitable, persisting in refugial retreats where the environment was compatible with their tolerances. In this study we combine modern geographic range data, phylogeny, Pleistocene paleoclimatic models, and isotopic records of changes in global mean annual temperature, to produce a temporally continuous model of geographic changes in potential habitat for 59 species of North American turtles over the past 320 Ka (three full glacial-interglacial cycles). These paleophylogeographic models indicate the areas where past climates were compatible with the modern ranges of the species and serve as hypotheses for how their geographic ranges would have changed in response to Quaternary climate cycles. We test these hypotheses against physiological, genetic, taxonomic and fossil evidence, and we then use them to measure the effects of Quaternary climate cycles on species distributions. Patterns of range expansion, contraction, and fragmentation in the models are strongly congruent with (i) phylogeographic differentiation; (ii) morphological variation; (iii) physiological tolerances; and (iv) intraspecific genetic variability. Modern species with significant interspecific differentiation have geographic ranges that strongly fluctuated and repeatedly fragmented throughout the Quaternary. Modern species with low genetic diversity have geographic distributions that were highly variable and at times exceedingly small in the past. Our results reveal the potential for paleophylogeographic models to (i) reconstruct past geographic range modifications, (ii) identify geographic processes that result in genetic bottlenecks; and (iii) predict threats due to anthropogenic climate change in the future.