The effect of phylogeographic history on species boundaries: a comparative framework in Hyla tree frogs

Genomic, Phenotypic and Environmental Correlates of Speciation in the Midwife Toads (Alytes)

Speciation, i.e., the formation of new species, implies that diverging populations evolve genetic, phenotypic or ecological factors that promote reproductive isolation (RI), but the relative contributions of these factors remain elusive. Here we test which of genomic, bioacoustic, morphological, and environmental differences best predicts RI across a continuum of divergence in the midwife toads (genus Alytes), a group of Western Mediterranean amphibians, using a total evidence approach. We found that, without strong geographic barriers to dispersal, the extent of introgression across hybrid zones between phylogeographic lineages, which should reflect the strength of RI, predominantly covaries with genomic divergence. Overall phenotypic differentiation becomes substantial only between well established, fully isolated species. These results suggest that speciation in midwife toads initially involve cryptic lineages, which probably evolve RI through intrinsic (genetic) hybrid incompatibilities. As they continue to diverge, these nascent species eventually differentiate externally, which potentially enforces pre-mating barriers and facilitates sympatry. This speciation scenario has practical implications for species delimitation, notably when using hybrid zones and divergence thresholds as proxies for reproductive isolation.

Developing under artificial light is not so bright: Oxidative stress as a physiological response to light pollution across amphibian life stages

Light pollution has emerged as a new environmental challenge. Natural light has long been a selecting force in shaping physiological, morphological and behavioral traits across species. The introduction of artificial light at night (ALAN) disrupts the natural light-dark cycle, affecting a broad range of biological processes. In amphibians, ALAN has been suggested as a potential driver behind demographic declines in recent years. In this study, we investigated how exposure to cold LED light affects the development, body length, and oxidative status of the European tree frog (Hyla arborea). Individuals were exposed to ALAN from egg hatching (GS25) to two developmental endpoints: the larvae (GS38) and the juveniles. The experiment included three night-time light conditions: a natural dark control (<0.3 lx), LED light intensity at 20 lx, and LED light intensity at 90 lx. Larvae reared under artificial light reached metamorphosis in a shorter period and exhibited larger body sizes. Nocturnal light also induced changes in oxidative status across life stages. Tadpoles under both 20 lx and 90 lx levels experienced greater oxidative damage of lipids. Likewise, juveniles had higher oxidative damage of proteins and elevated activities of the first-line antioxidant defense enzymes (SOD, CAT and GSH-Px) in response to both light intensities. Our results show that exposure to ALAN shapes development and body size of amphibians, which appears to be accompanied by disturbances in redox balance. Overall, we highlight that ALAN can be a factor affecting amphibian populations, and thus it should be taken into consideration in future conservation strategies.

Genomic and bioacoustic variation in a midwife toad hybrid zone: A role for reinforcement?

Hybrid zones, i.e., geographic areas where diverging lineages meet, hybridize and eventually mix their genomes, offer opportunities to understand the mechanisms behind reproductive isolation and speciation. Hybrid zones are particularly well suited to study reinforcement, i.e., the process by which selection against hybridization increases reproductive barriers, which, in anuran amphibians, is typically expressed by increased divergence in advertisement calls-the main cue to assortative mating-in parapatric ranges. Using mitochondrial barcoding (16S sequences), population genomics (thousands of SNPs) and bioacoustic analyses (four call parameters), we examine the hybrid zone between two incipient species of midwife toads (Alytes obstetricans and A. almogavarii) in southern France, with the purposes of locating their transition, measuring genetic introgression, and documenting potential signatures of reinforcement. We map range boundaries in the Eastern Pyrenees and the southwestern foothills of the Massif Central, namely along the Ariège valley and the Montagne Noire area. Similarly to another transition between these species in Spain, we found the hybrid zone to be narrow, involving geographically restricted gene flow (~20 km wide allele frequency clines) and barrier loci (i.e., loci resisting introgression), both suggestive of partial post-zygotic isolation (hybrid incompatibilities). The calls of the species overlap less inside than outside the hybrid zone, due to a reduction of their standing variation rather than a shift towards distinctive variants. While neutral causes cannot be excluded, this pattern follows the general expectations of reinforcement, yet without reproductive character displacement. Our study highlights the potential of amphibian hybrid zones to assess the genetic and behavioral drivers of reproductive isolation in statu nascendi and under various evolutionary contexts.

Spatial genetic structure in European marbled newts revealed with target enrichment by sequence capture

European marbled newts come in two species that have abutting ranges. The northern species, Triturus marmoratus, is found in France and the northern part of the Iberian Peninsula, whereas the southern species, T. pygmaeus, is found in the southwestern corner of the Iberian Peninsula. We study the intraspecific genetic differentiation of the group because morphological data show geographical variation and because the Iberian Peninsula is a recognized center of speciation and intraspecific genetic diversity for all kinds of organisms, amphibians included. We use target enrichment by sequence capture to generate c. 7 k nuclear DNA markers. We observe limited genetic exchange between the species, which confirms their distinctiveness. Both species show substantial genetic structuring that is only in part mirrored by morphological variation. Genetically differentiated groups are found in the south (T. marmoratus) and west (T. pygmaeus) of the species ranges. Our observations highlight the position of the Iberian Peninsula as a hotspot for genetic differentiation.

The phylogeography of Middle Eastern tree frogs in Israel

Western Palearctic treefrogs of the genus Hyla provide an example of a morphologically and ecologically cryptic group. Up to three distinct Hyla species have been proposed as resident in Israel and this number has consistently been subject to taxonomical debates. Here, we analyzed 16S rRNA and COI gene fragments of 658 individuals sampled at 47 pools in nine regions across Israel and the West Bank in order to resolve the taxonomic status of Hyla frogs. We generated both Bayesian and Maximum Likelihood phylogenies, and constructed time-calibrated trees to provide an evolutionary and historical context of sequence variations. We further applied SAMOVA as well as Monmonier's maximum-difference algorithm to study the genetic structure among populations and to identify potential zones acting as barriers to gene flow across locations. Our results revealed two distinct haplogroups for each gene fragment, with 95% CI divergence times dated from 8.9-17.1 Mya (16S) and 7.1-23.6 Mya (COI), respectively. SAMOVA and barrier analyses partitioned the populations into three groups. Our results highlight that, while there are probably only two Hyla species in Israel, one population of one of the species might qualify as a separate evolutionarily significant unit. Our findings elucidate the taxonomic status of Hyla frogs in Israel and provide the basis for determining appropriate management and conservation priorities.

Acknowledging more biodiversity without more species

Delimiting and naming biodiversity is a vital step toward wildlife conservation and research. However, species delimitation must be consistent across biota so that the limited resources available for nature protection can be spent effectively and objectively. To date, newly discovered lineages typically are either left undescribed and thus remain unprotected or are being erroneously proposed as new species despite mixed evidence for completed speciation, in turn contributing to the emerging problem of taxonomic inflation. Inspired by recent conceptual and methodological progress, we propose a standardized workflow for species delimitation that combines phylogenetic and hybrid zone analyses of genomic datasets ("genomic taxonomy"), in which phylogeographic lineages that do not freely admix are ranked as species, while those that have remained fully genetically compatible are ranked as subspecies. In both cases, we encourage their formal taxonomic naming, diagnosis, and description to promote social awareness toward biodiversity. The use of loci throughout the genome overcomes the unreliability of widely used barcoding genes when phylogeographic patterns are complex, while the evaluation of divergence and reproductive isolation unifies the long-opposed concepts of lineage species and biological species. We suggest that a shift in conservation assessments from a single level (species) toward a two-level hierarchy (species and subspecies) will lead to a more balanced perception of biodiversity in which both intraspecific and interspecific diversity are valued and more adequately protected.

Potential hotspots of amphibian roadkill risk in Spain

We test a forecasting strategy to identify potential hotspots of amphibian roadkill, combining the spatial distribution of amphibians, their relative risk of collision with vehicles and data on road density in Spain. We extracted a large dataset from studies reporting road casualties of 39 European amphibian species and then estimated the 'relative roadkill risk' of species as the frequency of occurrence of casualties for each amphibian and standardized by the range of distribution of the species in Europe. Using a map with the spatial distribution of Spanish amphibians at a spatial resolution of 10 × 10 Km squares, we estimated the 'cumulative relative risk of roadkill' for each amphibian assemblage as the sum of risk estimates previously calculated for each species. We also calculated the total length of roads in each square (road density). Finally, combining all layers of information, we elaborated a forecasting map highlighting the potential amphibian roadkill risk across Spain. Our findings are relevant to suggest areas that should be focused on at more detailed spatial scales. Additionally, we found that the frequency of roadkill was unrelated to the evolutionary distinctiveness score and conservation status of amphibian species, while was positively correlated with their distribution range.

A synoptic review of the Amphibians of Iran: bibliography, taxonomy, synonymy, distribution, conservation status, and identification key to the eggs, larvae, and adults

This study provides an illustrated account, a comprehensive update of the systematics, and a bibliography of the 15 species of anurans in five families, eight genera; and of the six species of urodeles in two families, four genera in Iran. Bufonidae, with eight species, is the most diverse family; Salamandridae has five species and Ranidae has four species. This study also presents updated identification keys for the eggs, larvae, and metamorphosed amphibians of Iran. We designated specimen NMW 19855.1 as neotype of Pelophylax persicus (Schneider, 1799) comb. nov.. Along with distribution maps obtained from all the reliable localities and museum specimens known at this time, the modelled habitat of species, and for the first time, the National Red List of amphibians based on the IUCN red list categories and criteria. Based on our evaluation we propose to categorize Bufo eichwaldi, Paradactylodon persicus, Neurergus derjugini, and N. kaiseri as Vulnerable at National Red List, and to move Bufotes (Calliopersa) luristanicus, B. (C.) surdus, Firouzophrynus olivaceus, and Rana pseudodalmatina from the category of Least Concern (LC) to Near Threatened (NT). The National Red List of amphibians that we propose has significant implications for endangered species management and conservation. Forty-one percent of amphibian species in Iran are endemic to the country, and more than forty percent of the Iranian amphibians are at risk of extinction. Zagros Mountain forest and Hyrcaniain forests have more than 80% (i.e. 18 species) of the diversity of Iranian amphibians. A considerable amount of scientific literature published on Iranian amphibians in Persian language is not easily accessible to researchers outside Iran. This monograph attempts to remedy the situation and provides broader access to international herpetology. We recognize that taxonomy is always in a state of flux, and the names and synonymies used here reflect our current view.

Dinarippiger gen. nov. (Tettigoniidae: Bradyporinae: Ephippigerini), a new saddle bush-cricket genus for Ephippiger discoidalis Fieber, 1853 from the Dinaric karst

A new genus of the tribe Ephippigerini, Dinarippiger Skejo, Kasalo, Fontana et Tvrtković gen. nov., is described based on the characters of occiput coloration, tegmina coloration, cerci and pronotum shape. The new genus is morphologically intermediate between the genera Ephippiger Berthold, 1827 and Uromenus Bolívar, 1878, and presently includes only Dalmatian Saddle Bush Cricket, Dinarippiger discoidalis (Fieber, 1853) comb. nov., hitherto known as Ephippiger discoidalis Fieber, 1853. The species inhabits NE Italy (mainly Carso Triestino), SW Slovenia, Croatia, Bosnia & Herzegovina, and Montenegro, i.e., islands and karst habitats along the eastern Adriatic coast, with isolated findings in Albania and Italy. Its prominent variation in size and coloration has already produced many synonyms (= limbata Fischer, 1853, = limbata var. major Krauss, 1879, = limbata var. minor Krauss, 1879, = selenophora Fieber, 1853, = sphacophila Krauss, 1879), which may suggest that what is currently regarded as a single species could represent a complex of distinct species with restricted distributions. This study also presents an annotated distribution map and a bioacoustic analysis of D. discoidalis comb. nov. Further research, especially adopting molecular methods, is necessary to assess possible cryptic diversity within the genus Dinarippiger gen. nov. and elucidate its evolutionary history.

Easternmost distribution of Bufo bufo (Linnaeus, 1758) in Türkiye: implications for the putative contact zone between B. bufo and B. verrucosissimus

The geographic range of a species is crucial for obtaining information on the exact distribution of the species. The geographic data are important for delimiting distinct species or exploring the degree of differentiation among different populations of a species. The local details of species boundaries facilitate the study of the importance of phylogeographic background, secondary contacts, and hybrid zones, along with the relations between the species and its extrinsic environmental factors. In the present study, the range boundaries of Bufo bufo and Bufo verrucosissimus in the north-eastern region of Türkiye were delineated using an integrative taxonomic approach that utilized a combination of molecular and morphological data. According to the mtDNA results of the present study, B. bufo inhabits a single distribution from İyidere town to Çayeli town in Rize, while B. verrucosissimus is distributed from Şavşat town of Artvin to Ardeşen town in Rize. In addition, the two species coexist in Pazar, Hemşin, and Çamlıhemşin towns in Rize. The demographic analyses indicated a distinct population expansion for the B. verrucosissimus species after the Last Glacial Maximum, while the same did not occur for B. bufo. The univariate and multivariate statistical analyses conducted for the morphological data of the two species corroborated the presence of a putative contact zone between B. bufo and B. verrucosissimus. In summary, the present study resolved the non-distinct geographic boundaries between B. bufo and B. verrucosissimus species and also revealed the easternmost distribution of B. bufo in Türkiye. In addition, important evidence on the putative contact zone between the two species was indicated using an integrative taxonomic approach.

Phylogeography of Kisaura Ross (Trichoptera: Philopotamidae) of the Japanese Archipelago and the character displacement evolution observed in a secondary contact area between genetically differentiated intra-specific lineages

In this study we present evidence of an ongoing speciation event in the Japanese Archipelago. The Kisaura nozakii species complex (Philopotamidae) is a group of a small trichopteran insects, and it has been reported that this species has large differences between populations in the copulatory organs (genitalia) of males. Therefore, we conducted molecular phylogenetic analyses based on the mtDNA COI region (658 bp) and genome-wide nDNA single nucleotide polymorphisms (SNPs) (16 254 loci) using the K. nozakii species complex and a sister-species, K. borealis. We also conducted principal component analyses (PCA) based on the forewing length and seven morphological characteristics of the male genitalia. It was revealed that the K. nozakii species complex is composed of multiple allopatrically differentiated genetic lineages. However, two genetic lineages were distributed sympatrically/parapatrically in the western area of the Japanese Archipelago, and the occurrence of ‘character displacement’ in the male genitalia was detected in this area. We consider that the evolutionary event in this species complex has high potential to become a model case for elucidating the processes of speciation and species diversification.

From Gondwana to the Yellow Sea, evolutionary diversifications of true toads Bufo sp. in the Eastern Palearctic and a revisit of species boundaries for Asian lineages

Taxa with vast distribution ranges often display unresolved phylogeographic structures and unclear taxonomic boundaries resulting in hidden diversity. This hypothesis-driven study reveals the evolutionary history of Bufonidae, covering the phylogeographic patterns found in Holarctic bufonids from the West Gondwana to the phylogenetic taxonomy of Asiatic true toads in the Eastern Palearctic. We used an integrative approach relying on fossilized birth-death calibrations, population dynamics, gene-flow, species distribution, and species delimitation modeling to resolve the biogeography of the clade and highlight cryptic lineages. We verified the near-simultaneous Miocene radiations within Western and Eastern Palearctic Bufo, c. 14.49-10.00 Mya, temporally matching with the maximum dust outflows in Central Asian deserts. Contrary to earlier studies, we demonstrated that the combined impacts of long dispersal and ice-age refugia equally contributed to the current genetic structure of Bufo in East Asia. Our findings reveal a climate-driven adaptation in septentrional Eastern Asian Bufo, explaining its range shifts toward northern latitudes. We resolve species boundaries within the Eastern Palearctic Bufo, and redefine the taxonomic and conservation units of the northeastern species: B. sachalinensis and its subspecies.

Peering into the Darkness: DNA Barcoding Reveals Surprisingly High Diversity of Unknown Species of Diptera (Insecta) in Germany

Determining the size of the German insect fauna requires better knowledge of several megadiverse families of Diptera and Hymenoptera that are taxonomically challenging. This study takes the first step in assessing these “dark taxa” families and provides species estimates for four challenging groups of Diptera (Cecidomyiidae, Chironomidae, Phoridae, and Sciaridae). These estimates are based on more than 48,000 DNA barcodes (COI) from Diptera collected by Malaise traps that were deployed in southern Germany. We assessed the fraction of German species belonging to 11 fly families with well-studied taxonomy in these samples. The resultant ratios were then used to estimate the species richness of the four “dark taxa” families (DT families hereafter). Our results suggest a surprisingly high proportion of undetected biodiversity in a supposedly well-investigated country: at least 1800–2200 species await discovery in Germany in these four families. As this estimate is based on collections from one region of Germany, the species count will likely increase with expanded geographic sampling.

Morphological diversification of Mediterranean anurans: the roles of evolutionary history and climate

Investigation of the ecological and evolutionary mechanisms governing the origin and diversification of species requires integrative approaches that often have to accommodate strong discordance among datasets. A common source of conflict is the combination of morphological and molecular characters with different evolutionary rates. Resolution of these discordances is crucial to assess the relative roles of different processes in generating and maintaining biodiversity. Anuran amphibians provide many examples of morphologically similar, genetically divergent lineages, posing questions about the relative roles of phylogeny and ecological factors in phenotypic evolution. We focused on three circum-Mediterranean anuran genera (Hyla, Alytes and Discoglossus), characterizing morphological and environmental disparity and comparing diversity patterns across biological levels of organization. Using a comparative phylogenetic framework, we tested how shared ancestry and climatic factors come together to shape phenotypic diversity. We found higher morphological differentiation within Hyla and Alytes than in Discoglossus. Body size and limb morphology contributed most to inter- and intraspecific morphological variation in Hyla and Alytes, but there was no strong phylogenetic signal, indicating that shared ancestry does not predict patterns of phenotypic divergence. In contrast, we uncovered a significant association between morphology and climatic descriptors, supporting the hypothesis that morphological disparity between species results from adaptive evolution.

Mass of genes rather than master genes underlie the genomic architecture of amphibian speciation

The genetic architecture of speciation, i.e., how intrinsic genomic incompatibilities promote reproductive isolation (RI) between diverging lineages, is one of the best-kept secrets of evolution. To directly assess whether incompatibilities arise in a limited set of large-effect speciation genes, or in a multitude of loci, we examined the geographic and genomic landscapes of introgression across the hybrid zones of 41 pairs of frog and toad lineages in the Western Palearctic region. As the divergence between lineages increases, phylogeographic transitions progressively become narrower, and larger parts of the genome resist introgression. This suggests that anuran speciation proceeds through a gradual accumulation of multiple barrier loci scattered across the genome, which ultimately deplete hybrid fitness by intrinsic postzygotic isolation, with behavioral isolation being achieved only at later stages. Moreover, these loci were disproportionately sex linked in one group (Hyla) but not in others (Rana and Bufotes), implying that large X-effects are not necessarily a rule of speciation with undifferentiated sex chromosomes. The highly polygenic nature of RI and the lack of hemizygous X/Z chromosomes could explain why the speciation clock ticks slower in amphibians compared to other vertebrates. The clock-like dynamics of speciation combined with the analytical focus on hybrid zones offer perspectives for more standardized practices of species delimitation.

Drainage basins serve as multiple glacial refugia for alpine habitats in the Sierra Nevada Mountains, California

The evolutionary histories of alpine species are often directly associated with responses to glaciation. Deep divergence among populations and complex patterns of genetic variation have been inferred as consequences of persistence within glacier boundaries (i.e., on nunataks), while shallow divergence and limited genetic variation are assumed to result from expansion from large refugia at the edge of ice shields (i.e., massifs de refuge). However, for some species, dependence on specific microhabitats could profoundly influence their spatial and demographic response to glaciation, and such a simple dichotomy may obscure the localization of actual refugia. In this study, we use the Nebria ingens complex (Coleoptera: Carabidae), a water-affiliated ground beetle lineage, to test how drainage basins are linked to their observed population structure. By analysing mitochondrial COI gene sequences and genome-wide single nucleotide polymorphisms, we find that the major drainage systems of the Sierra Nevada Mountains in California best explain the population structure of the N. ingens complex. In addition, we find that an intermediate morphotype within the N. ingens complex is the product of historical hybridization of N. riversi and N. ingens in the San Joaquin basin during glaciation. This study highlights the importance of considering ecological preferences in how species respond to climate fluctuations and provides an explanation for discordances that are often observed in comparative phylogeographical studies.

An odyssey out of Africa: an integrative review of past and present invasions by the Mediterranean tree frog (Hyla meridionalis)

An illustration of the human footprint on biodiversity are the faunal movements that have accompanied commercial and cultural exchanges between civilizations throughout history. In this article, we provide an integrative review of biogeographical and archaeological knowledge to understand these processes for the Mediterranean tree frog (Hyla meridionalis), an African species that has reached several regions of south-western Europe and associated islands. By re-analysing molecular (multilocus mitochondrial and genome-wide markers) and bioclimatic (ecological niche modelling) data, we give a comprehensive picture of the genetic diversity, structure and environmental suitability of H. meridionalis across its natural and exotic ranges, which in turn offer specific clues to the putative routes of colonization and associated events. Long-term monitoring efforts suggest northwards shifts of local range margins, potentially due to global warming, and we further demonstrate that this species is the most frequent amphibian to travel via the food supply chain in Western Europe. High dispersal ability, ecological tolerance and proximity to human settlements have made H. meridionalis a recurrent witness to the complexity and diversity of the civilizations that ruled the Western Mediterranean.

Reconstruction of past distribution for the Mongolian toad, Strauchbufo raddei (Anura: Bufonidae) using environmental modeling

The use of ecological models enables determining the current distribution of species, but also their past distribution when matching climatic conditions are available. In specific cases, they can also be used to determine the likelihood of fossils to belong to the same species-under the hypothesis that all individuals of a species have the same ecological requirements. Here, using environmental modeling, we reconstructed the distribution of the Mongolian toad, Strauchbufo raddei, since the Last Glacial Maximum and thus covering the time period between the Late Pleistocene and the Holocene. We found the range of the species to have shifted over time, with the LGM population clustered around the current southern range of the species, before expanding east and north during the Pleistocene, and reaching the current range since the mid-Holocene. Finally, we determined that the ecological conditions during the life-time of the mid-Pleistocene fossils attributed to the species in Europe were too different from the one of the extant species or fossils occurring at the same period in Asia to belong to the same species.