Integrating hybrid zone analyses in species delimitation: lessons from two anuran radiations of the Western Mediterranean.

Molecular ecologists often rely on phylogenetic evidence for assessing the species-level systematics of newly discovered lineages. Alternatively, the extent of introgression at phylogeographic transitions can provide a more direct test to assign candidate taxa into subspecies or species categories. Here, we compared phylogenetic versus hybrid zone approaches of species delimitation in two groups of frogs from the Western Mediterranean region (Discoglossus and Pelodytes), by using genomic data (ddRAD). In both genera, coalescent analyses recovered almost all nominal taxa as "species". However, the least-diverged pairs D. g. galganoi/jeanneae and P. punctatus/hespericus admix over hundreds of kilometers, suggesting that they have not yet developed strong reproductive isolation and should be treated as conspecifics. In contrast, the comparatively older D. scovazzi/pictus and P. atlanticus/ibericus form narrow contact zones, consistent with species distinctiveness. Due to their complementarity, we recommend taxonomists to combine phylogenomics with hybrid zone analyses to scale the gray zone of speciation, i.e., the evolutionary window separating widely admixing lineages versus nascent reproductively isolated species. The radically different transitions documented here conform to the view that genetic incompatibilities accumulating with divergence generate a weak barrier to gene flow for long periods of time, until their effects multiply and the speciation process then advances rapidly. Given the variability of the gray zone among taxonomic groups, at least from our current abilities to measure it, we recommend to customize divergence thresholds within radiations to categorize lineages for which no direct test of speciation is possible.

Integrative systematic revision of the Montseny brook newt (Calotriton arnoldi), with the description of a new subspecies.

The Montseny brook newt (Calotriton arnoldi), a glacial relict endemic to a small, isolated massif in northeast Spain, is considered the only Critically Endangered urodele in Europe. Its restricted range is divided by a deep valley that acts as an impassable barrier to dispersal, separating two isolated metapopulations (Western and Eastern) that correspond to independent lineages with different evolutionary trajectories, based on genetic and genomic data. Here, we address the ecological differentiation between lineages and discuss its potential effect on the phenotypic distinctness of each lineage. Based on multiple lines of evidence, we formally describe the Western Montseny brook newt as a new subspecies: Calotriton arnoldi laietanus ssp. nov. Finally, our study underscores the importance of considering taxonomic progress in the conservation policies of endangered species, ensuring appropriate management and protection of the newly described taxa.

Genomic insights into the Montseny brook newt (Calotriton arnoldi), a Critically Endangered glacial relict.

The Montseny brook newt (Calotriton arnoldi), considered the most endangered amphibian in Europe, is a relict salamandrid species endemic to a small massif located in northeastern Spain. Although conservation efforts should always be guided by genomic studies, those are yet scarce among urodeles, hampered by the extreme sizes of their genomes. Here, we present the third available genome assembly for the order Caudata, and the first genomic study of the species and its sister taxon, the Pyrenean brook newt (Calotriton asper), combining whole-genome and ddRADseq data. Our results reveal significant demographic oscillations which accurately mirrored Europe's climatic history. Although severe bottlenecks have led to depauperate genomic diversity and long runs of homozygosity along a gigantic genome, inbreeding might have been avoided by assortative mating strategies. Other life history traits, however, seem to have been less advantageous, and the lack of land dispersal has driven to exceptional levels of population fragmentation.

Surviving on the edge: present and future effects of climate warming on the common frog (Rana temporaria) population in the Montseny massif (NE Iberia).

The Montseny massif shelters the southernmost western populations of common frogs (Rana temporaria) that live in a Mediterranean climate, one which poses a challenge for the species' persistence in a scenario of rising temperatures. We evaluated the effect of climate change at three levels. First, we analysed if there has been an advancement in the onset of spawning period due to the increase in temperatures. Second, we analysed the impact of climatic variables on the onset of the spawning period and, third, how the distribution of this species could vary according to the predictions with regard to rising temperatures for the end of this century. From 2009 to 2021, we found there had been an increase in temperatures of 0.439 °C/decade, more than the 0.1 °C indicated by estimates for the second half of the previous century. We found an advancement in the onset of the reproduction process of 26 days/decade for the period 2009-2022, a change that has been even more marked during the last eight years, when data were annually recorded. Minimum temperatures and the absence of frost days in the week prior to the onset of the spawning period determine the start of reproduction. Predictions on habitat availability for spawning provided by climatic niche analysis for the period 2021-2100 show a potential contraction of the species range in the Montseny and, remarkably, much isolation from the neighbouring populations.

Biogeographic Inferences on the Evolutionary History of the King Cobra (Ophiophagus hannah, Cantor 1836) Species Complex.

King cobra (Ophiophagus hannah) is a snake widely distributed through southeastern tropical Asia, but in two separate subpopulations: one located in the Western Ghats (western Indian Peninsula) and the other much more extensive, ranging between the southern slopes of the Himalayas, Assam, Indochina to southeastern China. Similarly, it also appears in numerous tropical archipelagos such as Indonesia, the Philippines, and the Andaman Islands, but surprisingly it is absent from other large islands like Sri Lanka and Taiwan. In this study, we evaluated how climate could be shaping the distribution of this snake and estimated the future distribution of the species utilizing ecological niche modelling. To evaluate the effect of paleoclimatic conditions on the genetic structure of this species we performed Bayesian phylogenetic analysis under a molecular clock using mitochondrial DNA. Our analyses indicated that the current distribution of O. hannah is strongly influenced by the availability of humid climate conditions. King cobras have a long evolutionary history reflected in the appearance of four main mitochondrial lineages before the Pliocene (the Western Ghats, southeastern mainland Asia, Luzon, and Indonesia), congruently with paleoclimatic models that indicated the availability of suitable conditions for this species in these refugia during the glacial cycles. Climate history could explain the absence of O. hannah in Sri Lanka and Taiwan due to the absence of suitable climatic corridors when these islands were connected to the mainland (20,000 years ago). Future projections (2050‒2070) did not suggest significant range shifts in the region, even considering the worst global warming scenarios.

Jailed in the mountains: Genetic diversity and structure of an endemic newt species across the Pyrenees.

The Pyrenees represent a natural laboratory for biogeographic, evolutionary and ecological research of mountain fauna as a result of the high variety of habitats and the profound effect of the glacial and interglacial periods. There is a paucity of studies providing a detailed insight into genetic processes and better knowledge on the patterns of genetic diversity and how they are maintained under high altitude conditions. This is of particular interest when considering the course of past climate conditions and glaciations in a species which is considered site tenacious, with long generation times. Here we analyzed the genetic patterns of diversity and structure of the endemic Pyrenean brook newt (Calotriton asper) along its distribution range, with special emphasis on the distinct habitat types (caves, streams, and lakes), and the altitudinal and geographical ranges, using a total set of 900 individuals from 44 different localities across the Pyrenean mountain range genotyped for 19 microsatellite loci. We found evidence for a negative longitudinal and positive altitudinal gradient of genetic diversity in C. asper populations. The fact that genetic diversity was markedly higher westwards is in accordance with other Pyrenean species. However, the impact of altitudinal gradient on the genetic diversity seems to differ from other species, and mostly from other amphibians. We found that lower altitudes can act as a barrier probably because the lowlands do not provide a suitable habitat for C. asper. Regarding the distinct habitat types, caves had significantly lower values of genetic diversity compared to streams or lakes. The mean FST value was relatively high (0.304) with maximum values as high as 0.771, suggesting a highly structured total population. Indeed, populations were grouped into five subclusters, the eastern populations (cluster 1) remained grouped into two subclusters and the central-western Pyrenees (cluster 2) into three subclusters. The increase of isolation with geographical distance is consistent with the population structure detected. In conclusion, C. asper seems to be adapted to high altitude mountain habitats, and its genetic diversity is higher in the western Pyrenees. In terms of conservation priority, we consider more relevant the populations that represent a reservoir of genetic diversity.

Diversity, distribution and conservation of the terrestrial reptiles of Oman (Sauropsida, Squamata).

In the present work, we use an exceptional database including 5,359 records of 101 species of Oman's terrestrial reptiles together with spatial tools to infer the spatial patterns of species richness and endemicity, to infer the habitat preference of each species and to better define conservation priorities, with especial focus on the effectiveness of the protected areas in preserving this unique arid fauna. Our results indicate that the sampling effort is not only remarkable from a taxonomic point of view, with multiple observations for most species, but also for the spatial coverage achieved. The observations are distributed almost continuously across the two-dimensional climatic space of Oman defined by the mean annual temperature and the total annual precipitation and across the Principal Component Analysis (PCA) of the multivariate climatic space and are well represented within 17 out of the 20 climatic clusters grouping 10% of the explained climatic variance defined by PC1 and PC2. Species richness is highest in the Hajar and Dhofar Mountains, two of the most biodiverse areas of the Arabian Peninsula, and endemic species richness is greatest in the Jebel Akhdar, the highest part of the Hajar Mountains. Oman's 22 protected areas cover only 3.91% of the country, including within their limits 63.37% of terrestrial reptiles and 50% of all endemics. Our analyses show that large areas of the climatic space of Oman lie outside protected areas and that seven of the 20 climatic clusters are not protected at all. The results of the gap analysis indicate that most of the species are below the conservation target of 17% or even the less restrictive 12% of their total area within a protected area in order to be considered adequately protected. Therefore, an evaluation of the coverage of the current network of protected areas and the identification of priority protected areas for reptiles using reserve design algorithms are urgently needed. Our study also shows that more than half of the species are still pending of a definitive evaluation by the International Union for Conservation of Nature (IUCN).

Transcriptomic and macroevolutionary evidence for phenotypic uncoupling between frog life history phases.

Anuran amphibians undergo major morphological transitions during development, but the contribution of their markedly different life-history phases to macroevolution has rarely been analysed. Here we generate testable predictions for coupling versus uncoupling of phenotypic evolution of tadpole and adult life-history phases, and for the underlying expression of genes related to morphological feature formation. We test these predictions by combining evidence from gene expression in two distantly related frogs, Xenopus laevis and Mantidactylus betsileanus, with patterns of morphological evolution in the entire radiation of Madagascan mantellid frogs. Genes linked to morphological structure formation are expressed in a highly phase-specific pattern, suggesting uncoupling of phenotypic evolution across life-history phases. This gene expression pattern agrees with uncoupled rates of trait evolution among life-history phases in the mantellids, which we show to have undergone an adaptive radiation. Our results validate a prevalence of uncoupling in the evolution of tadpole and adult phenotypes of frogs.

Phylogenetic relationships of Semaphore geckos (Squamata: Sphaerodactylidae: Pristurus) with an assessment of the taxonomy of Pristurus rupestris.

A molecular phylogeny of the sphaerodactylid geckos of the genus Pristurus is inferred based on an alignment of 1845 base pairs (bp) of concatenated mitochondrial (12S) and nuclear (acm4, cmos, rag1 and rag2) genes for 80 individuals, representing 18 of the 23-26 species, and the three subspecies of P. rupestris. The results indicate that P. rupestris is polyphyletic and includes two highly divergent clades: the eastern clade, found in coastal Iran and throughout the Hajar Mountain range in northern Oman and eastern UAE; and the western clade, distributed from central coastal Oman, through Yemen, Saudi Arabia and north to southern Jordan. Inferred haplotype networks for the four nuclear genes show that the eastern and western clades of "P. rupestris" are highly differentiated and do not share any alleles. Moreover, although the two clades are differentiated by a morphological multivariate analysis, no one character or set of characters was found to be diagnostic. Based on the molecular analysis of specimens from the type locality of P. rupestris rupestris, the name P. rupestris is applied to the eastern clade. The name that should apply to the western clade cannot be clarified until morphological and genetic data for "P. rupestris" is available from the vicinity of Bosaso, Somalia, and therefore we refer to it as Pristurus sp. 1. The phylogenetic tree of Pristurus supports the hypothesis that P. celerrimus is sister to all the other species in the analyses and that the Socotra Archipelago was independently colonized a minimum of two times.

Integrative phylogeography of Calotriton newts (Amphibia, Salamandridae), with special remarks on the conservation of the endangered Montseny brook newt (Calotriton arnoldi).

The genus Calotriton includes two species of newts highly adapted to live in cold and fast-flowing mountain springs. The Pyrenean brook newt (Calotriton asper), restricted to the Pyrenean region, and the Montseny brook newt (Calotriton arnoldi), endemic to the Montseny massif and one of the most endangered amphibian species in Europe. In the present manuscript, we use an integrative approach including species distribution modeling (SDM), molecular analyses of mitochondrial and nuclear DNA sequence data and morphology to unravel the historical processes that have contributed to shaping the biogeography and genetic structure of the genus Calotriton, with special emphasis on the conservation of C. arnoldi. The results of the molecular analyses confirm that, despite having originated recently, being ecologically similar and geographically very close, there is no signal of hybridization between C. asper and C. arnoldi. SDM results suggest that tough environmental conditions on mountains tops during glacial periods, together with subsequent warmer periods could have prevented the contact between the two species. Within the critically endangered C. arnoldi, a high genetic structure is revealed despite its extremely small distribution range compared to C. asper. Haplotype networks, AMOVA and SAMOVA analyses suggest that two distinct groups of populations can be clearly differentiated with absence of gene flow. This is in concordance with morphological differentiation and correlates with its geographical distribution, as the two groups are situated on the eastern and western sides of a river valley that acts as a barrier. The genetic and morphological results are highly important for the ongoing conservation program of C. arnoldi and strongly justify the management of this species into at least two independent evolutionary significant units (eastern and western sectors) to guarantee the long-term population viability.

Incidence of pigmented skin tumors in a population of wild Montseny brook newt (Calotriton arnoldi).

We report the presence of pigmented skin tumors in three populations of the endangered amphibian Montseny brook newt, Calotriton arnoldi, one of the European amphibian species with the smallest distribution range (40 km(2) in the Montseny Natural Park, Catalonia, Spain). Examination of one of the tumors by light microscopy was consistent with chromatophoroma and was most suggestive of a melanophoroma. Tumors were not found in juveniles. In adults, only two of three populations were affected. The proportions of males and females affected were not significantly different, but there was a positive correlation between body size and presence of tumors in both sexes. The etiology of chromatophoromas remains unknown but, in our study, they do not appear to have been caused by water quality or Ultraviolet B.