Jelka Crnobrnja-Isailović

Phylogeography of two European newt species — discordance between mtDNA and morphology

The newts Triturus vulgaris and Triturus montandoni are sister species that exhibit contrasting levels of intraspecific morphological variation. Triturus vulgaris has a broad Eurasiatic distribution encompassing both formerly glaciated and unglaciated areas and shows substantial morphological differentiation in the southern part of its range, while T. montandoni, confined to the Carpathians, is morphologically uniform. We analysed sequence variation of two mtDNA fragments of the total length of c. 1850 bp in 285 individuals of both species collected from 103 localities. Phylogenetic analysis of 200 unique haplotypes defined 12 major clades, their age estimated at c. 4.5–1.0 million years (Myr). Most of the older clades were found in the southern part of the range, and also in central Europe, mainly in Romania. The distribution of mtDNA clades points to the existence of several glacial refugia, located in the Caucasus region, Anatolia, the Balkan Peninsula, Italy, and more to the north in central Europe. The concordance between mtDNA based phylogeny and the distribution of T. vulgaris subspecies was weak. Triturus montandoni haplotypes did not form a monophyletic group. Instead they were found in six clades, in five of them mixed with T. vulgaris haplotypes, most likely as a result of past or ongoing hybridization and multiple introgression of mtDNA from T. vulgaris to T. montandoni. Patterns of sequence variation within clades suggested long‐term demographic stability in the southern groups, moderate and relatively old demographic growth in the populations inhabiting central Europe, and high growth in some of the groups that colonized northern parts of Europe after the last glacial maximum.

Tracing glacial refugia of Triturus newts based on mitochondrial DNA phylogeography and species distribution modeling

IntroductionThe major climatic oscillations during the Quaternary Ice Age heavily influenced the distribution of species and left their mark on intraspecific genetic diversity. Past range shifts can be reconstructed with the aid of species distribution modeling and phylogeographical analyses. We test the responses of the different members of the genus Triturus (i.e. the marbled and crested newts) as the climate shifted from the previous glacial period (the Last Glacial Maximum, ~21 Ka) to the current interglacial.ResultsWe present the results of a dense mitochondrial DNA phylogeography (visualizing genetic diversity within and divergence among populations) and species distribution modeling (using two different climate simulations) for the nine Triturus species on composite maps.ConclusionsThe combined use of species distribution modeling and mitochondrial phylogeography provides insight in the glacial contraction and postglacial expansion of Triturus. The combined use of the two independent techniques yields a more complete understanding of the historical biogeography of Triturus than both approaches would on their own. Triturus newts generally conform to the ‘southern richness and northern purity’ paradigm, but we also find more intricate patterns, such as the absence of genetic variation and suitable area at the Last Glacial Maximum (T. dobrogicus), an ‘extra-Mediterranean’ refugium in the Carpathian Basin (T. cristatus), and areas where species displaced one another postglacially (e.g. T. macedonicus and western T. karelinii). We provide a biogeographical scenario for Triturus, showing the positions of glacial refugia, the regions that were postglacially colonized and the areas where species displaced one another as they shifted their ranges.

Multilocus species tree analyses resolve the radiation of the widespread Bufo bufo species group (Anura, Bufonidae)

New analytical methods are improving our ability to reconstruct robust species trees from multilocus datasets, despite difficulties in phylogenetic reconstruction associated with recent, rapid divergence, incomplete lineage sorting and/or introgression. In this study, we applied these methods to resolve the radiation of toads in the Bufo bufo (Anura, Bufonidae) species group, ranging from the Iberian Peninsula and North Africa to Siberia, based on sequences from two mitochondrial and four nuclear DNA regions (3490 base pairs). We obtained a fully-resolved topology, with the recently described Bufo eichwaldi from the Talysh Mountains in south Azerbaijan and Iran as the sister taxon to a clade including: (1) north African, Iberian, and most French populations, referred herein to Bufo spinosus based on the implied inclusion of populations from its type locality and (2) a second clade, sister to B. spinosus, including two sister subclades: one with all samples of Bufo verrucosissimus from the Caucasus and another one with samples of B. bufo from northern France to Russia, including the Apennine and Balkan peninsulas and most of Anatolia. Coalescent-based estimations of time to most recent common ancestors for each species and selected subclades allowed historical reconstruction of the diversification of the species group in the context of Mediterranean paleogeography and indicated a long evolutionary history in this region. Finally, we used our data to delimit the ranges of the four species, particularly the more widespread and historically confused B. spinosus and B. bufo, and identify potential contact zones, some of which show striking parallels with other co-distributed species.

Schoolchildren and One of the Most Unpopular Animals: Are They Ready to Protect Snakes?

ABSTRACT For cultural reasons and due to the narrow vision of environmental policy makers, most conservation efforts focus on a few charismatic species and consequently neglect the majority of others under threat; many unpopular species are even killed in large numbers with little concern. Redressing this bias through educational programs is therefore important. Snakes are unpopular animals; they suffer from human harassment in most places and many populations have declined worldwide. Consequently, they provide suitable substrate to better improve conservation education in schoolchildren. Responses to a questionnaire administered to 2,570 schoolchildren (7–14 years old) from 10 countries showed that many children liked snakes and that most of the students wanted to see snakes protected. Such counterintuitive results were supported by the explanations given by the children, notably the reasons they liked or disliked snakes. Previous physical contact with snakes was often associated with snake likeability. We also found strong and expected differences between countries: for instance, where venomous snakes represented a health risk, most children declared to be afraid of snakes. Overall, our results negate the simplistic, and previously unverified, adult view that snakes are necessarily perceived as frightening animals, thereby justifying their persecution. This study provides an encouraging message, suggesting that it is not compulsory to focus on charismatic animals to convince children to protect wildlife.

Conservation phylogeography: does historical diversity contribute to regional vulnerability in European tree frogs (Hyla arborea)?

Documenting and preserving the genetic diversity of populations, which conditions their long‐term survival, have become a major issue in conservation biology. The loss of diversity often documented in declining populations is usually assumed to result from human disturbances; however, historical biogeographic events, otherwise known to strongly impact diversity, are rarely considered in this context. We apply a multilocus phylogeographic study to investigate the late‐Quaternary history of a tree frog (Hyla arborea) with declining populations in the northern and western part of its distribution range. Mitochondrial and nuclear polymorphisms reveal high genetic diversity in the Balkan Peninsula, with a spatial structure moulded by the last glaciations. While two of the main refugial lineages remained limited to the Balkans (Adriatic coast, southern Balkans), a third one expanded to recolonize Northern and Western Europe, loosing much of its diversity in the process. Our findings show that mobile and a priori homogeneous taxa may also display substructure within glacial refugia (‘refugia within refugia’) and emphasize the importance of the Balkans as a major European biodiversity centre. Moreover, the distribution of diversity roughly coincides with regional conservation situations, consistent with the idea that historically impoverished genetic diversity may interact with anthropogenic disturbances, and increase the vulnerability of populations. Phylogeographic models seem important to fully appreciate the risks of local declines and inform conservation strategies.

Lizards from urban areas are more asymmetric: using fluctuating asymmetry to evaluate environmental disturbance.

The increase in human activities that leads to wildlife decline and species extinction poses an urgent need for simple indicators of environmental stress in animal populations. Several studies have suggested that fluctuating asymmetry (FA) can be an easy, direct measure of developmental instability because it is associated to environmental stress and, as such, it can be a useful indicator of population disturbance. We examined three different morphological traits in urban and rural populations of the common wall lizard (Podarcis muralis) to test whether anthropogenic disturbance causes an increase in FA. Compared to rural populations, urban ones showed higher levels of FA in all analyzed traits, thus providing evidence that FA can respond to anthropogenic disturbance. However, we also found significant differences in FA among traits, where femoral pores and subdigital lamellae, traits with a functional relevance, were more stable developmentally compared to supracilliar granules which have no evident function. Unsigned FA [abs(right-left)] exhibited significant, but weak, positive correlations among traits, indicating that developmental noise does not have a uniform effect across characters and thus questioning the view of developmental stability as an organism-wide property. The degree of signed FA (right-left) was more similar between structurally associated traits, possibly as an outcome of morphological integration. In conclusion, our results demonstrate that FA can be a reliable indicator of disturbance provided that it is analyzed on multiple traits simultaneously and examined at the population level.

Effects of Environmental Disturbance on Phenotypic Variation: An Integrated Assessment of Canalization, Developmental Stability, Modularity, and Allometry in Lizard Head Shape

When populations experience suboptimal conditions, the mechanisms involved in the regulation of phenotypic variation can be challenged, resulting in increased phenotypic variance. This kind of disturbance can be diagnosed by using morphometric tools to study morphological patterns at different hierarchical levels and evaluate canalization, developmental stability, integration, modularity, and allometry. We assess the effect of urbanization on phenotypic variation in the common wall lizard (Podarcis muralis) by using geometric morphometrics to assess disturbance to head shape development. The head shapes of urban lizards were more variable and less symmetric, suggesting that urban living is more likely to disturb development. Head shape variation was congruent within and across individuals, which indicated that canalization and developmental stability are two related phenomena in these organisms. Furthermore, urban lizards exhibited smaller mean head sizes, divergent size-shape allometries, and increased deviation from within-group allometric lines. This suggests that mechanisms regulating head shape allometry may also be disrupted. The integrated evaluation of several measures of developmental instability at different hierarchical levels, which provided in this case congruent results, can be a powerful methodological guide for future studies, as it enhances the detection of environmental disturbances on phenotypic variation and aids biological interpretation of the results.

Timeframe of speciation inferred from secondary contact zones in the European tree frog radiation (Hyla arborea group)

BackgroundHybridization between incipient species is expected to become progressively limited as their genetic divergence increases and reproductive isolation proceeds. Amphibian radiations and their secondary contact zones are useful models to infer the timeframes of speciation, but empirical data from natural systems remains extremely scarce. Here we follow this approach in the European radiation of tree frogs (Hyla arborea group). We investigated a natural hybrid zone between two lineages (Hyla arborea and Hyla orientalis) of Mio-Pliocene divergence (~5 My) for comparison with other hybrid systems from this group.ResultsWe found concordant geographic distributions of nuclear and mitochondrial gene pools, and replicated narrow transitions (~30xa0km) across two independent transects, indicating an advanced state of reproductive isolation and potential local barriers to dispersal. This result parallels the situation between H. arborea and H. intermedia, which share the same amount of divergence with H. orientalis. In contrast, younger lineages show much stronger admixture at secondary contacts.ConclusionsOur findings corroborate the negative relationship between hybridizability and divergence time in European tree frogs, where 5 My are necessary to achieve almost complete reproductive isolation. Speciation seems to progress homogeneously in this radiation, and might thus be driven by gradual genome-wide changes rather than single speciation genes. However, the timescale differs greatly from that of other well-studied amphibians. General assumptions on the time necessary for speciation based on evidence from unrelated taxa may thus be unreliable. In contrast, comparative hybrid zone analyses within single radiations such as our case study are useful to appreciate the advance of speciation in space and time.

A reinvestigation of phylogeny and divergence times of the Ablepharus kitaibelii species complex (Sauria, Scincidae) based on mtDNA and nuDNA genes

Morphological and DNA data support that the East Mediterranean snake-eyed skink Ablepharus kitaibelii represents a species complex that includes four species A. kitaibelii, A. budaki, A. chernovi, and A. rueppellii, highlighting the need of its taxonomic reevaluation. Here, we used Bayesian and Maximum Likelihood methods to estimate the phylogenetic relationships of all members of the complex based on two mitochondrial (cyt b, 16S rRNA) and two nuclear markers (MC1R, and NKTR) and using Chalcides, Eumeces, and Eutropis as outgroups. The biogeographic history of the complex was also investigated through the application of several phylogeographic (BEAST) and biogeographic (BBM) analyses. Paleogeographic and paleoclimatic data were used to support the inferred phylogeographic patterns. The A. kitaibelli species complex exhibits high genetic diversity, revealing cases of hidden diversity and cases of non-monophyletic species such as A. kitaibelii and A. budaki. Our results indicate that A. pannonicus branches off first and a group that comprises specimens of A. kitaibelli and A. budaki from Kastelorizo Island group (southeast Greece) and southwest Turkey, respectively is differentiated from the rest A. kitaibelli and A. budaki populations and may represent a new species. The estimated divergence times place the origin of the complex in the Middle Miocene (∼16Mya) and the divergence of most currently recognized species in the Late Miocene. The inferred ancestral distribution suggests that the complex originated in Anatolia, supposing that several vicariance and dispersal events that are related with the formation of the Mid-Aegean Trench, the Anatolian Diagonal and the orogenesis of the mountain chains in southern and eastern Anatolia have led to current distribution pattern of A. kitaibelii species complex in the Balkans and Middle East.

Hidden diversity in the Podarcis tauricus (Sauria, Lacertidae) species subgroup in the light of multilocus phylogeny and species delimitation

The monophyletic species subgroup of Podarcis tauricus is distributed in the western and southern parts of the Balkans, and includes four species with unresolved and unstudied inter- and intra-specific phylogenetic relationships. Using sequence data from two mitochondrial and three nuclear genes and applying several phylogenetic methods and species delimitation approaches to an extensive dataset, we have reconstructed the phylogeny of the Podarcis wall lizards in the Balkans, and re-investigated the taxonomic status of the P. tauricus species subgroup. Multilocus analyses revealed that the aforementioned subgroup consists of five major clades, with P. melisellensis as its most basal taxon. Monophyly of P. tauricus sensu stricto is not supported, with one of the subspecies (P. t. ionicus) displaying great genetic diversity (hidden diversity or cryptic species). It comprises five, geographically distinct, subclades with genetic distances on the species level. Species delimitation approaches revealed nine species within the P. tauricus species subgroup (P. melisellensis, P. gaigeae, P. milensis, and six in the P. tauricus complex), underlining the necessity of taxonomic re-evaluation. We thus synonymize some previously recognized subspecies in this subgroup, elevate P. t. tauricus and P. g. gaigeae to the species level and suggest a distinct Albanian-Greek clade, provisionally named as the P. ionicus species complex. The latter clade comprises five unconfirmed candidate species that call for comprehensive studies in the future.