The taxonomy of the Tarentola mauritanica species complex (Gekkota: Phyllodactylidae): Bayesian species delimitation supports six candidate species

Lack of country-wide systematic herpetology collections in Portugal jeopardizes future research and conservation

Natural History Collections (NHCs) represent the world's largest repositories of long-term biodiversity datasets. Specimen collection and voucher deposition has been the backbone of NHCs since their inception, but recent decades have seen a drastic decline in rates of growth via active collecting. Amphibians and reptiles are amongst the most threatened zoological groups on the planet and are historically underrepresented in most worldwide NHCs. As part of an ongoing project to review the Portuguese zoological collections in the country's NHCs, herpetological data from its three major museums and smaller collections was gathered and used to examine the coverage and representation of the different taxa extant in Portugal. These collections are not taxonomically, geographically, or temporally complete. Approximately 90% of the Portuguese herpetological taxa are represented in the country's NHCs, and around half of the taxa are represented by less than 50 specimens. Geographically, the collections cover less than 30% of the country's territory and almost all of the occurring taxa have less than 10% of their known distribution represented in the collections. A discussion on the implications for science of such incomplete collections and a review of the current status of Portuguese NHCs is presented.

Colonization routes uncovered in a widely introduced Mediterranean gecko, Tarentola mauritanica

In this study, we aimed to understand the contemporary and ancient colonization routes of the Moorish gecko, Tarentola mauritanica, using simple sequence repeats. By analyzing the genetic diversity of populations in different regions, we found that Morocco is the genetic diversity hotspot for the species, followed by the Iberian Peninsula. However, historical gene flow estimates identified the Iberian Peninsula, not Morocco, as the primary contributor of colonizing individuals, along with continental Italy to a lesser extent. Currently, mainland Italy is the main source of introduced individuals, likely due to the plant nursery trade. The study suggests that human-facilitated introductions from various geographical origins, with numerous regions colonized through continental Italy during two distinct periods, are responsible for the recurrent entry of individuals belonging to the European lineage of T. mauritanica into the Mediterranean and Macaronesia. These findings can inform better monitoring surveys and conservation programs by identifying putative current colonization routes of alien species.

Trophic interactions of an invasive gecko in an endemic-rich oceanic island: Insights using DNA metabarcoding

Understanding the trophic interactions of introduced predators is key for evidence-based management of biological invasions. This is particularly important in oceanic islands, where predator-prey networks often include numerous endemic and range-restricted species. Geckos are successful island colonizers and in recent years numerous species have established populations in a wide array of oceanic islands. One such species is the Moorish gecko (Tarentola mauritanica), which has colonized multiple islands across the Mediterranean basin, Caribbean and Macaronesia. The species was first reported in Madeira Island in 1993 and over the last 30 years has colonized most of the islands' southern coast and expanded to the nearby island of Porto Santo. Here, we used DNA metabarcoding to provide the first insights into the diet of this successful colonizer in its introduced range. The species' diet was mainly composed of ground-dwelling arthropods belonging to the families Porcellionidae (Isopoda), Julidae (Diplopoda) and Formicidae (Hymenoptera). The diet richness and composition were not affected by neither sex nor size of adult geckos, instead they both change across populations. However, trophic niche-width differed among size classes, with smaller geckos feeding on a wider range of prey. We identified over 160 different Operational Taxonomic Units in the diet of T. mauritanica, with 21.6% of them belonging to introduced invertebrates and 13.6% to native species. Native prey taxa included the endemic Madeira wall lizard (Teira dugesii), the sole native reptile to Madeira. We also detected several agricultural pests and disease vectors in the diet of this exotic predator, and 19 taxa identified as prey had not yet been recorded to Madeira. Of these, several are serious agricultural pests, highlighting how this introduced gecko can be used as a natural sampler, in particular for the early detection of invasive arthropod pests. This study emphasizes the importance of trophic studies for monitoring the impacts of introduced predators in fragile insular systems.

Reptile biodiversity in Souss-Massa National Park: an internationally important hotspot in the Mediterranean region

Souss-Massa National Park (SMNP) is Morocco’s first coastal national park, created to preserve the high diversity of its continental and marine environments. Reptiles play an essential role in balancing SMNP ecosystems, yet little work has been done to study this fauna. The present work aims at providing the first reptile inventory of SMNP since its establishment in 1991. During the period 2019 to 2020, several field surveys were carried out at 30 sites using time-constrained visual encounter surveys (TCVES), with a total sampling effort of 300 person-hours. An inventory of 23 reptile species (including four endemic species) was obtained by combining TCVES results with additional data recorded during random encounters or provided by SMNP researchers. Based on TCVES data, both sampling effort and inventory completeness were evaluated by constructing sample-based accumulation curves and calculating non-parametric estimators (Chao 1, Chao 2, Jackknife 1 and Jackknife 2). These species richness estimators suggest that the current inventory is likely to be complete. Despite its small surface area, SMNP contains nearly 20% of all known Moroccan reptile species and constitutes an important biodiversity hotspot for reptiles in the Mediterranean Region. In terms of reptile conservation concern, five species in SMNP are classified as “vulnerable”, while two species are classified as “near threatened” on the IUCN Red List, underscoring the importance of protected areas for those species.

The Moorish gecko, Tarentola mauritanica Linnaeus, 1758 (Squamata, Phyllodactylidae), in Albania

We summarised the records of Tarentola mauritanica in Albania. The first record for the country was published in 1989 from Sazan Island; however, it was never confirmed again. Another population was found recently in an urban habitat of Vlorë City near the main city port, probably of human-mediated origin.

A review of the introduced herpetofauna of Mexico and Central America, with comments on the effects of invasive species and biosecurity methodology

Among the principal causes producing detrimental effects on global biodiversity are introductions of alien species. Very few attempts to control introduced amphibians and reptiles in Middle America (Mexico and Central America) can be identified, so listings are provided for 24 exotic species, 16 translocated species, and 11 species that were removed from the introduced species listing because of lack of substantiating evidence that they are from established populations. Biosecurity methods are also identified that can be applied for preventing, controlling, and managing introduced and especially invasive species.

Morphometric and population genomic evidence for species divergence in the Chimarrichthys fish complex of the Tibetan Plateau

The uplift of the Tibetan Plateau altered the environmental conditions of the local area substantially. Here, we conducted a comprehensive investigation based on morphometrics, population genomics, and climatic factors to evaluate phenotypic and genome-level variations in a radiation of Chimarrichthys catfish endemic to the Plateau. Discriminant function analysis showed phenotypic differences of Chimarrichthys between rivers with respect to elevation. Genetic structure analysis based on 6606 single nucleotide polymorphisms (SNPs) deduced genetic differences between rivers, and species delimitation indicated that the Chimarrichthys fish complex could be divided into three species. Restriction site-associated DNA tags were mapped to the gene sets of Glyptosternon maculatum, and matches were searched against databases for Gene Ontology annotation. Genomic regions exhibiting marked differences among localities represented a range of biological functions, including growth (gdf11), bone development (bmp8a), cellular response to light stimulus (opn3), regulation of the rhodopsin-mediated signalling pathway (grk1), immune response (rag1 and ung), reproductive process (antxr2), and regulation of intracellular iron levels (ireb2). The tag44126, where gene gdf11 is located, was identified as an outlier exhibiting divergence between rivers with altitude differences, and the SNP is thymine (T) in Dadu and Yalong River (~2700 m), but guanine (G) in Jinsha and Qingyi rivers (~2200 and ~ 684 m), suggesting a possible effect of altitude on its differentiation.

Evidence from Tarentola mauritanica (Gekkota: Phyllodactylidae) helps validate thermography as a tool to infer internal body temperatures of lizards

Infrared (IR) thermal imaging has become an increasingly popular tool to measure body temperature of animals. The high-resolution data it provides with short lag and minimum disturbance makes it an appealing tool when studying reptile thermal ecology. However, due to the common phenomenon of regional heterothermy and surface-to-core temperature gradients, it is essential to select the appropriate body part to measure and provide calibrations to accurately infer internal body temperatures. This work follows from a previous study on lacertid lizards to assess the reliability of thermography-measured body temperatures, from several body locations, as a proxy for internal body temperature in lizards. This study focuses on the Moorish gecko, Tarentola mauritanica, due to its distant phylogenetic relationship and its different ecology and morphology from the previously tested species. A total of 60 adult geckos of both sexes and of a range of sizes were tested in thermal gradients and subjected to a sequence of randomly assorted treatments of heating and cooling. The temperatures of the animals were periodically measured with a thermal camera at six different body parts and, immediately after, the cloacal temperature was then measured with a thermocouple probe. Body parts' temperatures, obtained thermographically, were regressed against cloacal temperature using OLS regression and the pairwise correlations were tested using Spearman coefficients. Relationships among all body parts and between all body parts and the cloaca were strong in all cases (R² > 0.87, Spearman Correlation > 0.95). The observed pattern was very similar to those previously obtained from lacertid lizards. Ultimately, the eye proved to provide the best overall proxy for internal temperature, when accounting for both the slope and intercept of the regression. Hence, this study provides further support for the establishment of the eye as the standard location to infer internal body temperatures of lizards through thermography.

A morphometric assessment of species boundaries in a widespread anole lizard (Squamata: Dactyloidae)

Cryptic species – genetically distinct species that are morphologically difficult to distinguish – present challenges to systematists. Operationally, cryptic species are very difficult to identify and sole use of genetic data or morphological data can fail to recognize evolutionarily isolated lineages. We use morphometric data to test species boundaries hypothesized with genetic data in the North Caribbean bark anole (Anolis distichus), a suspected species complex. We use univariate and multivariate analyses to test if candidate species based on genetic data can be accurately diagnosed. We also test alternative species delimitation scenarios with a model fitting approach that evaluates normal mixture models capable of identifying morphological clusters. Our analyses reject the hypothesis that the candidate species are diagnosable. Neither uni- nor multivariate morphometric data distinguish candidate species. The best-supported model included two morphological clusters; however, these clusters were uneven and did not align with a plausible species divergence scenario. After removing two related traits driving this result, only one cluster was supported. Despite substantial differentiation revealed by genetic data, we recover no new evidence to delimit species and refrain from taxonomic revision. This study highlights the importance of considering other types of data along with molecular data when delimiting species.

Phylogenomics and molecular species delimitation reveals great cryptic diversity of leaf-toed geckos (Phyllodactylidae: Phyllodactylus), ancient origins, and diversification in Mexico

We utilize the efficient GBS technique to obtain thousands of nuclear loci and SNPs to reconstruct the evolutionary history of Mexican leaf-toed geckos (Phyllodactylus). Through the incorporation of unprecedented sampling for this group of geckos, in combination with genomic data analysis, we generate mostly consistent phylogenetic hypotheses using two approaches: supermatrix and coalescent-based inference. All topologies depict three, mutually exclusive major clades. Clade I comprises P. bordai and all species closer to P. bordai than to any other Phyllodactylus. Clade II comprises P. nocticolus and all species closer to P. nocticolus than to any other Phyllodactylus. Clade III comprises P. tuberculosus and all species closer to P. tuberculosus than to any other Phyllodactylus. Analyses estimate the age for the most recent common ancestor of Phyllodactylus in the Eocene (~43 mya), and the ancestors of each major clade date to the Eocene-Oligocene transition (32-36 mya). This group includes one late-Eocene lineage (P. bordai), Oligocene lineages (P. paucituberculatus, P. delcampi), but also topological patterns that indicate a recent radiation occurred during the Pleistocene on islands in the Gulf of California. The wide spatial and temporal scale indicates a complex and unique biogeographic history for each major clade. The 33 species delimited by BPP and stepping-stone BFD*coalescent based genomic approaches reflect this history. This diversity delimited for Mexican leaf-toed geckos demonstrates a vast underestimation in the number of species based on morphological data alone.

Untangling a mess of worms: Species delimitations reveal morphological crypsis and variability in Southeast Asian semi-aquatic earthworms (Almidae, Glyphidrilus)

Semi-aquatic freshwater earthworms in the genus Glyphidrilus from Southeast Asia are characterized by both an extreme morphological crypsis among divergent phylogenetic lineages and a high morphological variability within the same phylogenetic lineages. The present study provides a new taxonomic framework for this problematic genus in SE Asia by integrating DNA sequence and morphological data. When single-locus and multilocus multispecies coalescent-based (MSC) species delimitation methods were applied to DNA sequence data, they usually yielded highly incongruent results compared to morphology-based species identifications. This suggested the presence of several cryptic species and high levels of intraspecific morphological variation. Applying reciprocal monophyly to the cytochrome c oxidase subunit 1 (COI) gene tree allowed us to propose the existence of 33 monophyletic species. Yet, often substantially more molecular operational taxonomic units (MOTUs) were obtained when species delimitation was based on COI and 16S rRNA sequences. In contrast, the ITS1 and ITS2 sequences suggested fewer MOTUs and did not recover most of the monophyletic species from the Mekong basin. However, several of these latter taxa were better supported when MSC species delimitation methods were applied to the combined mtDNA and ITS datasets. The ITS2 secondary structure retrieved one unnamed Mekong basin species that was not uncovered by the other methods when applied to ITS2 sequences. In conclusion, based on an integrative taxonomic workflow, 26 Glyphidrilus candidate species were retained and two remained to be confirmed. As such, this study provides evidence to suggest nine species new to science and to synonymize 12 nominal morphospecies. It also illustrates that the uncritical use of COI as a universal DNA barcode may overestimate species diversity because COI may be unable to distinguish between divergent conspecific lineages and different candidate species.

Hidden species diversity in Pachyhynobius: A multiple approaches species delimitation with mitogenomes

The lack of distinct morphological features of cryptic species is a hard problem for taxonomy, especially when the taxa are closely related with considerable amounts of ancestral polymorphism. Lately, intensive coalescent-based analyses involving multiple loci have become the preferred method to assess the extent of genetic distinctiveness in otherwise phenotypically similar populations. Previously, phylogenetic studies on Pachyhynobius shangchengensis uncovered five extremely deeply divergent clades, which suggested that this species may be a cryptic species complex. In this study, we used the complete mitochondrial genome data and samples from the entire range of stout salamander (Pachyhynobius), as well as publicly available mitochondrial genomes to assess species boundaries within this genus using a suite of diverse methodologies (e.g. general mixed Yule coalescent model, Automatic Barcode Gap Discovery). The phylogenetic relationships recovered two major groups within P. shangchengensis, with one group formed by four of the six extant populations and corresponding to the central and eastern range of the Dabie mountains, while the other group encompassed two other lineages in the north west of the Dabie mountain range. The species delimitation comparison within Pachyhynobius supported the presence of recognized species within the genus, and consensus was observed across methods for the existence of up to five cryptic species within what has been traditionally considered to be P. shangchengensis. While this implies the existence of four taxa in addition to the described P. shangchengensis species, morphological data and life history information are further required to contribute to the species definition. The observed pattern of genetic variation is likely the outcome of a discontinuous habitat combined with niche conservatism, which produced the sky-island effect observed in Pachyhynobius, and which led to formation of a hidden species diversity in this genus.

Integrating multilocus DNA data and 3D geometric morphometrics to elucidate species boundaries in the case of Pyrenaearia (Pulmonata: Hygromiidae)

To accurately delimit species the use of multiple character types is essential as all speciation processes are not equally reflected in different data (e.g. morphological, molecular or ecological characters). With the introduction of geometric morphometrics methods and advances in 3D technology, a comprehensive combination of molecular and morphological data has been enabled in groups where exhaustively quantifying and measuring morphological shape change was not possible before such as gastropod shells. In this study, we combined multilocus coalescent species delimitation methods with 3D geometric morphometrics of shell shape to delimit species within the land snail genus Pyrenaearia. A new taxonomic scheme was constructed for the genus identifying ten species. Two nominal species were synonymized and a hitherto unrecognized cryptic species was identified. Our findings support the importance of combining multiple lines of evidence as molecular and morphological data on their own do not yield the same information. Further, the integration of morphological and molecular data shows the importance of allometry in shell shape and suggests a combined effect of population history and selection in different environments on shells morphological variation. Our new taxonomy and phylogenetic reconstruction suggest that, besides the glacial cycles of the Pleistocene, passive dispersal and rock substrate complexity could also have been involved in the speciation of the genus.

Comparative phylogeography of amphibians and reptiles in Algeria suggests common causes for the east-west phylogeographic breaks in the Maghreb

A series of phylogeographic studies in the Maghreb identified a repeated pattern of deep genetic divergence between an eastern (Tunisia) and western (Morocco) lineage for several taxa but lack of sampling in Algeria made it difficult to know if the range limits between the eastern and western lineages were shared among taxa or not. To address this question, we designed a comparative phylogeographic study using 8 reptile and 3 amphibian species with wide distribution in the Maghreb as models. We selected species where previous studies had identified an East-West phylogeographic divide and collected sampled in Algeria to 1) examine whether the simple East-West divergence pattern still holds after filling the sampling gap in Algeria or if more complex diversity patterns emerge; 2) if the E-W pattern still holds, test whether the limits between the E and W clades are shared between species, suggesting that common historical process caused the E-W divergences; 3) if E-W limits are shared between species, use information on the age of the divergence to identify possible geological or climatic events that could have triggered these E-W differentiations. We found that the E-W pattern was generally maintained after additional sampling in Algeria and identified two common disjunction areas, one around the Algeria-Morocco border, the other one in Kabylia (central Algeria), suggesting that common historical mechanisms caused the E-W divergences in the Maghreb. Our estimates for the times to most common recent ancestors to the E and W clades span a wide range between the Messinian salinity crisis and the Plio-Pleistocene limit (except for one older split), suggesting different origins for the initial divergences and subsequent preservation of the E and W lineages in common climatic refugia in the west and the east of the Maghreb.

Biodiversity and the Species Concept-Lineages are not Enough

The nature and definition of species continue to be matters of debate. Current views of species often focus on their nature as lineages-maximal reproductive communities through time. Whereas many authors point to the Evolutionary Species Concept as optimal, in its original form it stressed the ecological role of species as well as their history as lineages, but most recent authors have ignored the role aspect of the concept, making it difficult to apply unambiguously in a time-extended way. This trend has been exacerbated by the application of methods and concepts emphasizing the notion of monophyly, originally applied only at higher levels, to the level of individuals, as well as by the current emphasis on molecular data. Hence, some current authors recognize units that are no more than probable exclusive lineages as species. We argue that biodiversity is inherently a phenotypic concept and that role, as manifested in the organismal extended phenotype, is a necessary component of the species concept. Viewing species as historically connected populations with unique role brings together the temporal and phenotypic natures of species, providing a clear way to view species both in a time-limited and time-extended way. Doing so alleviates perceived issues with "paraphyletic species" and returns the focus of species to units that are most relevant for biodiversity.

Molecular systematics, species delimitation and diversification patterns of the Phyllodactylus lanei complex (Gekkota: Phyllodactylidae) in Mexico

The description of cryptic gecko species worldwide has revealed both that many putative species are, in fact, conformed by a complex of morphologically conserved species that are genetically distinct and highly divergent, and that gecko species diversity could be underestimated. The taxonomy and species delimitation of geckos belonging to the genus Phyllodactylus is still controversial, 16 of which are distributed in Mexico and 13 are endemic. Although the large morphological variation shown by the Phyllodactylus species from Mexico has been amply documented, little is known about their genetic diversity and evolutionary relationships, and much less regarding cryptic speciation. Here, we included the most comprehensive sampling of populations and species of the Phyllodactylus lanei complex distributed in Mexico, and applied an analytical approach that included probabilistic phylogenetic analyses, jointly with species delimitation methods and Bayesian putative species validation analysis. Our results suggest the existence of 10 lineages within the complex, supporting the existence of cryptic species, and in great contrast with the current taxonomic proposal that includes only four subspecies. The most recent common ancestor (MRCA) for the P. lanei clade originated on the Early Eocene (∼54Mya), along the southern coasts of Mexico, followed by the highest diversification of the complex MRCA during the Eocene (34-56Mya). Lineages subsequently dispersed and diversified towards the northwest, and the diversification process ended with the most recent lineages inhabiting two islands on the coasts of Nayarit (Miocene; 5.5-23Mya). Our results highlight three vicariant events associated with the evolution of the lineages, two of them intimately related to the formation of the Sierra Madre del Sur and the Transmexican Volcanic Belt mountain ranges, main geographic barriers that isolated and facilitated the divergence and speciation in this group of geckos. Finally, we propose that there are 10 species in the P. lanei complex, from which four represent taxonomic changes and six are new species and require a formal description. We acknowledge that more analyses, including a detailed evaluation of morphological characters and use of more unlinked nuclear loci with enough variability, are needed to further support their taxonomic description.

Subterranean termite phylogeography reveals multiple postglacial colonization events in southwestern Europe

A long‐standing goal of evolutionary biology is to understand how paleoclimatic and geological events shape the geographical distribution and genetic structure within and among species. Using a diverse set of markers (cuticular hydrocarbons, mitochondrial and nuclear gene sequences, microsatellite loci), we studied Reticulitermes grassei and R. banyulensis, two closely related termite species in southwestern Europe. We sought to clarify the current genetic structure of populations that formed following postglacial dispersal from refugia in southern Spain and characterize the gene flow between the two lineages over the last several million years. Each marker type separately provided a fragmented picture of the evolutionary history at different timescales. Chemical analyses of cuticular hydrocarbons and phylogenetic analyses of mitochondrial and nuclear genes showed clear separation between the species, suggesting they diverged following vicariance events in the Late Miocene. However, the presence of intermediate chemical profiles and mtDNA introgression in some Spanish colonies suggests ongoing gene flow. The current genetic structure of Iberian populations is consistent with alternating isolation and dispersal events during Quaternary glacial periods. Analyses of population genetic structure revealed postglacial colonization routes from southern Spain to France, where populations underwent strong genetic bottlenecks after traversing the Pyrenees resulting in parapatric speciation.