Predicting speciation probability from replicated population histories

In this issue of Molecular Ecology, Yamasaki et al. (2020) use genetic data from extensive sampling of Rhinogobius goby fish across the Ryukyu Archipelago in Japan to demonstrate the parallel speciation of a freshwater form from an ancestral amphidromous form. They then show that ecosystem size strongly predicts the probability of speciation between the two forms across islands. In doing so, this study connects population‐level processes (microevolution) to broad‐scale biodiversity patterns (macroevolution), an important but understudied link in evolutionary biology. Moving forward, we can build on this research to (a) more directly determine how geographic, ecological and historical factors influence the different stages of the speciation process, and (b) understand whether mechanisms inferred from insular radiations extend to those on continents, where both demographic histories and environmental regimes are likely more complex.     

The evolution of microendemism in a reef fish (Hypoplectrus maya)

Marine species tend to have extensive distributions, which are commonly attributed to the dispersal potential provided by planktonic larvae and the rarity of absolute barriers to dispersal in the ocean. Under this paradigm, the occurrence of marine microendemism without geographic isolation in species with planktonic larvae poses a dilemma. The recently described Maya hamlet (Hypoplectrus maya, Serranidae) is exactly such a case, being endemic to a 50‐km segment of the Mesoamerican Barrier Reef System (MBRS). We use whole‐genome analysis to infer the demographic history of the Maya hamlet and contrast it with the sympatric and pan‐Caribbean black (H. nigricans), barred (H. puella) and butter (H. unicolor) hamlets, as well as the allopatric but phenotypically similar blue hamlet (H. gemma). We show that H. maya is indeed a distinct evolutionary lineage, with genomic signatures of inbreeding and a unique demographic history of continuous decrease in effective population size since it diverged from congeners just ~3,000 generations ago. We suggest that this case of microendemism may be driven by the combination of a narrow ecological niche and restrictive oceanographic conditions in the southern MBRS, which is consistent with the occurrence of an unusually high number of marine microendemics in this region. The restricted distribution of the Maya hamlet, its decline in both census and effective population sizes, and the degradation of its habitat place it at risk of extinction. We conclude that the evolution of marine microendemism can be a fast and dynamic process, with extinction possibly occurring before speciation is complete.     

Influence of gene flow on divergence dating – implications for the speciation history of Takydromus grass lizards

Dating the time of divergence and understanding speciation processes are central to the study of the evolutionary history of organisms but are notoriously difficult. The difficulty is largely rooted in variations in the ancestral population size or in the genealogy variation across loci. To depict the speciation processes and divergence histories of three monophyletic Takydromus species endemic to Taiwan, we sequenced 20 nuclear loci and combined with one mitochondrial locus published in GenBank. They were analysed by a multispecies coalescent approach within a Bayesian framework. Divergence dating based on the gene tree approach showed high variation among loci, and the divergence was estimated at an earlier date than when derived by the species‐tree approach. To test whether variations in the ancestral population size accounted for the majority of this variation, we conducted computer inferences using isolation‐with‐migration (IM) and approximate Bayesian computation (ABC) frameworks. The results revealed that gene flow during the early stage of speciation was strongly favoured over the isolation model, and the initiation of the speciation process was far earlier than the dates estimated by gene‐ and species‐based divergence dating. Due to their limited dispersal ability, it is suggested that geographical isolation may have played a major role in the divergence of these Takydromus species. Nevertheless, this study reveals a more complex situation and demonstrates that gene flow during the speciation process cannot be overlooked and may have a great impact on divergence dating. By using multilocus data and incorporating Bayesian coalescence approaches, we provide a more biologically realistic framework for delineating the divergence history of Takydromus.     

Nuclear locus divergence at the early stages of speciation in the Orchard Oriole complex

As two lineages diverge from one another, mitochondrial DNA should evolve fixed differences more rapidly than nuclear DNA due to its smaller effective population size and faster mutation rate. As a consequence, molecular systematists have focused on the criteria of reciprocal monophyly in mitochondrial DNA for delimiting species boundaries. However, mitochondrial gene trees do not necessarily reflect the evolutionary history of the taxa in question, and even mitochondrial loci are not expected to be reciprocally monophyletic when the speciation event happened very recently. The goal of this study was to examine mitochondrial paraphyly within the Orchard Oriole complex, which is composed of Icterus spurius (Orchard Oriole) and Icterus fuertesi (Fuertes' Oriole). We increased the geographic sampling, added four nuclear loci, and used a range of population genetic and coalescent methods to examine the divergence between the taxa. With increased taxon sampling, we found evidence of clear structure between the taxa for mitochondrial DNA. However, nuclear loci showed little evidence of population structure, indicating a very recent divergence between I. spurius and I. fuertesi. Another goal was to examine the genetic variation within each taxon to look for evidence of a past founder event within the I. fuertesi lineage. Based on the high amounts of genetic variation for all nuclear loci, we found no evidence of such an event – thus, we found no support for the possible founding of I. fuertesi through a change in migratory behavior, followed by peripheral isolates speciation. Our results demonstrate that these two taxa are in the earliest stages of speciation, at a point when they have fixed differences in plumage color that are not reflected in monophyly of the mitochondrial or nuclear DNA markers in this study. This very recent divergence makes them ideal for continued studies of species boundaries and the earliest stages of speciation.     

Molecular evidence for the coexistence of two sibling species in Pylaiella littoralis (Ectocarpales,Phaeophyceae) along the Brittany coast

The great phenotypic variability and the lack of diagnostic characters in the genus Pylaiella render the systematic study of this genus problematic. In this study, we investigated the diversity of Pylaiella littoralis along the Brittany (France) coast using a DNA barcoding multilocus approach with mitochondrial (cox1, nad1, and atp9) and chloroplastic (rbcL and atpB) markers associated with a population genetics approach using 10 microsatellite markers. In addition, spatio‐temporal sampling was conducted along the Brittany coast. We sampled 140 individuals from four sites located between Saint‐Malo and Concarneau (380 km) from April to October. Mitochondrial sequence data revealed the occurrence of two sibling species, with a minimum of 2.4% divergence between them. Microsatellite genotypic data congruently revealed two well‐supported clusters matching the two mitochondrial clades of Pylaiella. Although gene flow is limited between species, occurrence of genetic admixtures in some populations suggested that reproductive isolation is not complete. Our study highlighted the complementarity of barcoding and population genetics approaches to shed light on the evolutionary processes that lead to speciation.     

quaddRAD: a new high‐multiplexing and PCR duplicate removal ddRAD protocol produces novel evolutionary insights in a nonradiating cichlid lineage

The identification of thousands of variants across the genomes and their accurate genotyping are crucial for estimating the genetic parameters needed to address a host of molecular ecological and evolutionary questions. With rapid advances of massively parallel high‐throughput sequencing technologies, several methods have recently been developed to access genomewide data on population variation. One of the most successful and widely used techniques relies on the combination of restriction enzymes and sequencing‐by‐synthesis: restriction‐site‐associated DNA sequencing (RADSeq). We developed a new, more time‐ and cost‐efficient double‐digest RAD paired‐end protocol (quaddRAD) that simplifies and speeds up the identification of PCR duplicates and permits large‐scale multiplexing. Assessing its performance on a technical data set, we also applied the quaddRAD method on population samples of a Neotropical cichlid fish lineage (Archocentrus centrarchus) to assess its genetic structure and demographic history. While we identified allopatric interlake genetic divergence, most likely driven by drift, no signature of sympatric divergence was detected. This differs from what has been observed in the clade of Midas cichlids (Amphilophus citrinellus spp.), another cichlid lineage that inhabits the same lakes and shares a similar demographic history, but has evolved into small‐scale adaptive radiations via sympatric speciation. We demonstrate that quaddRAD is a robust and efficient method for genotyping a massive number and widely overlapping set of loci with high accuracy. Furthermore, the results on A. centrarchus open new research avenues providing an ideal system to investigate genome‐level mechanisms that could alter the speciation potential of different but closely related cichlid lineages.     

Genetic consequences of plant edaphic specialization to solfatara fields: Phylogenetic and population genetic analysis of Carex angustisquama (Cyperaceae)

Edaphic specialization is one of the main drivers of plant diversification and has multifaceted effects on population dynamics. Carex angustisquama is a sedge plant growing only on heavily acidified soil in solfatara fields, where only extremophytes can survive. Because of the lack of closely related species in similar habitats and its disjunct distribution, the species offers ideal settings to investigate the effects of adaptation to solfatara fields and of historical biogeography on the genetic consequences of plant edaphic specialization to solfatara fields. Here, genome‐wide single nucleotide polymorphisms were used to reveal the phylogenetic origin of C. angustisquama, and 16 expressed sequence tag–simple sequence repeat markers were employed to infer population demography of C angustisquama. Molecular phylogenetic analysis strongly indicated that C. angustisquama formed a monophyletic clade with Carex doenitzii, a species growing on nonacidified soil in the sympatric subalpine zone. The result of population genetic analysis showed that C. angustisquama has much lower genetic diversity than the sister species, and notably, all 16 loci were completely homozygous in most individuals of C. angustisquama. Approximate Bayesian computation analysis supported the model that assumed hierarchical declines of population size through its evolutionary sequence. We propose that the edaphic specialist in solfatara fields has newly attained the adaptation to solfatara fields in the process of speciation. Furthermore, we found evidence of a drastic reduction in genetic diversity in C. angustisquama, suggesting that the repeated founder effects associated with edaphic specialization and subsequent population demography lead to the loss of genetic diversity of this extremophyte in solfatara fields.     

Evidence for nonallopatric speciation among closely related sympatric Heliotropium species in the Atacama Desert

The genetic structure of populations of closely related, sympatric species may hold the signature of the geographical mode of the speciation process. In fully allopatric speciation, it is expected that genetic differentiation between species is homogeneously distributed across the genome. In nonallopatric speciation, the genomes may remain undifferentiated to a large extent. In this article, we analyzed the genetic structure of five sympatric species from the plant genus Heliotropium in the Atacama Desert. We used amplified fragment length polymorphisms (AFLPs) to characterize the genetic structure of these species and evaluate their genetic differentiation as well as the number of loci subject to positive selection using divergence outlier analysis (DOA). The five species form distinguishable groups in the genetic space, with zones of overlap, indicating that they are possibly not completely isolated. Among‐species differentiation accounts for 35% of the total genetic differentiation (F_ST = 0.35), and F_ST between species pairs is positively correlated with phylogenetic distance. DOA suggests that few loci are subject to positive selection, which is in line with a scenario of nonallopatric speciation. These results support the idea that sympatric species of Heliotropium sect. Cochranea are under an ongoing speciation process, characterized by a fluctuation of population ranges in response to pulses of arid and humid periods during Quaternary times.     

Occurrence and dominance of six Pacific Northwest conifer species

Questions: Can probability of occurrence and dominance be accurately estimated for six important conifer species with varying range sizes? Does range size impact the accuracy of species probability of occurrence models? Is species predicted probability of occurrence significantly related to observed dominance? Location: Pacific Northwest region, North America (60°–40°N, 140°–110°W). Methods: This study develops near range‐wide predictive distribution maps for six important conifer species (Pseudotsuga menziesii, Tsuga heterophylla, Pinus contorta, Thuja plicata, Larix occidentalis, and Picea glauca) using forest inventory data collected across the United States and Canada. Species model accuracies are compared with range size using a rank scoring system. A suite of climate and topographic predictor variables are used to investigate environmental constraints that limit species range and quantify relationships between species predicted probability of occurrence and dominance at both plot and landscape scales. Results: Evaluation statistics revealed significant and accurate probability of occurrence models were developed for all six species. Based on ranked evaluation statistics, Tsuga heterophylla had highest overall model accuracy (statistic rank score=5) and Pinus contorta the lowest (statistic rank score=17). Across species, ranked evaluation statistics also revealed a pattern of decreasing model accuracy with increasing range size. At plot level, correlations between dominance and probability of occurrence were weakly positive for all species with only half of the species having statistically significant correlations. Pseudotsuga menziesii had the highest correlation (r=0.36, P<0.001) and Thuja plicata lowest (r=0.038, P=0.799). At the 50‐km scale, correlations between dominance and probability of occurrence improved for all species except Pinus contorta. Pseudotsuga menziesii displayed the highest correlation (r=0.68, P<0.001) and Thuja plicata the lowest (r=0.07, P>0.709). Conclusions: Species probability of occurrence model accuracy decreased with increasing range size. The strength and significance of correlations between probability of occurrence and dominance varied considerably by species and across spatial scales. Apart from Pseudotsuga menziesii and L. occidentalis, the results suggest that probability of occurrence is not a consistently reliable surrogate for species dominance in Pacific Northwest forests. We demonstrate how the degree of correlation between species occurrence and dominance can be used as an indicator of how well predictions of occurrence characterize the optimal niche of a species.     

Late Pleistocene climate change promoted divergence between Picea asperata and P. crassifolia on the Qinghai–Tibet Plateau through recent bottlenecks

Divergence during the early stage of speciation can be driven by a population bottleneck via reduced gene flow and enhanced lineage sorting. In this study, we aimed to examine whether such bottlenecks occurred during the initial speciation of two closely related spruce species Picea asperata and P. crassifolia occurring on the Qinghai–Tibet Plateau (QTP). We analyzed sequences of three chloroplast, two mitochondrial DNA fragments and a further 13 nuclear loci from 216 individuals of the two species. Both species showed a low level of genetic diversity in contrast to other congeners occurring in the QTP and adjacent regions. The estimated population sizes of P. asperata and P. crassifolia are less than the ancestral population size before splitting. These results together with multiple statistical tests (Tajima's D, Fu and Li's D* and F*) suggest that these two species underwent recent bottlenecks. Based on approximate Bayesian computation (ABC), we also determined that the period of the population shrinkage was consistent with the interspecific divergence during the late Pleistocene. The reduced population sizes and the divergent selection may together have triggered the initial divergence under high gene flow between these two species. Our results therefore highlight the importance of climatic oscillations during the late Pleistocene in promoting speciation through changing demographic sizes of the ancestral species on the QTP and in adjacent regions.     

Distribution pattern of the epiphyte Neckera pennata on three spatial scales – importance of past landscape structure, connectivity and local conditions

We tested which factors explain the distribution pattern of the epiphytic moss Neckera pennata on three spatial scales using the framework of generalized linear models. First, we tested which factors explained its occurrence in forest stands in a 2500 ha landscape. At this scale, we also tested the effect of the historic landscape structure. We recorded its occurrence in all suitable stands. The occurrence probability increased with increasing present quantity of Acer platanoides, and with increasing present and past quantity of Fraxinus excelsior. The probability also increased with increasing connectivity to occupied stands. However, the connectivity to stands present in 1977 (recorded from infra‐red aerial photographs) explained more of the variation. This suggests that the regional metapopulation size of N. pennata has decreased during the past decades, and that its present distribution pattern reflects the age of the remaining stands, and the distribution of past dispersal sources in the landscape. Second, we tested which factors explained the occurrence and abundance on individual trees in three forest stands. Neckera pennata mainly occurred on Acer and Fraxinus stems. The most important variable in explaining occurrence probability was connectivity to surrounding occupied trees, which probably reflects the restricted dispersal range in this species. The abundance on occupied trees was also explained by this variable. The occurrence probability and abundance also increased with increasing tree diameter, probably reflecting the time that a tree has been available for colonization and the time since colonization, respectively. The occurrence probability and abundance furthermore decreased on strongly leaning (and deteriorating) trees. The occurrence probability increased with increasing bark roughness, probably reflecting increasing suitability regarding bark chemistry and moisture. Third, we tested its vertical distribution on occupied trees. The main distribution was below 1.6 m.     

Speciation in the marine crop Pyropia yezoensis (Bangiales,Rhodophyta)

In the marine red alga Pyropia yezoensis, commonly known in Japan as nori, sympatric occurrence of two cryptic species Pyropia sp. 2 and Pyropia sp. 3 on the same rock in a natural habitat has been confirmed by molecular analysis and detailed morphological observations. To confirm whether Pyropia sp. 2 and Pyropia sp. 3 were reproductively isolated in the sympatric population, 170 blades that had previously been studied using a maternally inherited plastid marker were examined with a nuclear gene marker. The results suggested that Pyropia sp. 2 and Pyropia sp. 3 with identical morphological features were reproductively isolated in the sympatric population and that they were different species based on the biological species concept. Although gametophytic blades of Pyropia were usually assumed to be haploid, 18 of 170 blades possessed both of the two genotypes derived from Pyropia sp. 2 and from Pyropia sp. 3. These results inferred that allodiploid blades were generated from the interspecific hybridization between these two cryptic species. The present findings provide insights for future studies on the speciation mechanism in seaweeds, particularly for genera that contain numerous species.     

Speciation in protists: Spatial and ecological divergence processes cause rapid species diversification in a freshwater chrysophyte

Although eukaryotic microorganisms are extremely numerous, diverse and essential to global ecosystem functioning, they are largely understudied by evolutionary biologists compared to multicellular macroscopic organisms. In particular, very little is known about the speciation mechanisms which may give rise to the diversity of microscopic eukaryotes. It was postulated that the enormous population sizes and ubiquitous distribution of these organisms could lead to a lack of population differentiation and therefore very low speciation rates. However, such assumptions have traditionally been based on morphospecies, which may not accurately reflect the true diversity, missing cryptic taxa. In this study, we aim to articulate the major diversification mechanisms leading to the contemporary molecular diversity by using a colonial freshwater flagellate, Synura sphagnicola, as an example. Phylogenetic analysis of five sequenced loci showed that S. sphagnicola differentiated into two morphologically distinct lineages approximately 15.4 million years ago, which further diverged into several evolutionarily recent haplotypes during the late Pleistocene. The most recent haplotypes are ecologically and biogeographically much more differentiated than the old lineages, presumably because of their persistent differentiation after the allopatric speciation events. Our study shows that in microbial eukaryotes, species diversification via the colonization of new geographical regions or ecological resources occurs much more readily than was previously thought. Consequently, divergence times of microorganisms in some lineages may be equivalent to the estimated times of speciation in plants and animals.     

Host population genetics and biogeography structure the microbiome of the sponge Cliona delitrix

Sponges occur across diverse marine biomes and host internal microbial communities that can provide critical ecological functions. While strong patterns of host specificity have been observed consistently in sponge microbiomes, the precise ecological relationships between hosts and their symbiotic microbial communities remain to be fully delineated. In the current study, we investigate the relative roles of host population genetics and biogeography in structuring the microbial communities hosted by the excavating sponge Cliona delitrix. A total of 53 samples, previously used to demarcate the population genetic structure of C. delitrix, were selected from two locations in the Caribbean Sea and from eight locations across the reefs of Florida and the Bahamas. Microbial community diversity and composition were measured using Illumina‐based high‐throughput sequencing of the 16S rRNA V4 region and related to host population structure and geographic distribution. Most operational taxonomic units (OTUs) specific to Cliona delitrix microbiomes were rare, while other OTUs were shared with congeneric hosts. Across a large regional scale (>1,000 km), geographic distance was associated with considerable variability of the sponge microbiome, suggesting a distance–decay relationship, but little impact over smaller spatial scales (<300 km) was observed. Host population structure had a moderate effect on the structure of these microbial communities, regardless of geographic distance. These results support the interplay between geographic, environmental, and host factors as forces determining the community structure of microbiomes associated with C. delitrix. Moreover, these data suggest that the mechanisms of host regulation can be observed at the population genetic scale, prior to the onset of speciation.     

Phylogeny and biogeography of the American live oaks (Quercus subsection Virentes): a genomic and population genetics approach

The nature and timing of evolution of niche differentiation among closely related species remains an important question in ecology and evolution. The American live oak clade, Virentes, which spans the unglaciated temperate and tropical regions of North America and Mesoamerica, provides an instructive system in which to examine speciation and niche evolution. We generated a fossil‐calibrated phylogeny of Virentes using RADseq data to estimate divergence times and used nuclear microsatellites, chloroplast sequences and an intron region of nitrate reductase (NIA‐i3) to examine genetic diversity within species, rates of gene flow among species and ancestral population size of disjunct sister species. Transitions in functional and morphological traits associated with ecological and climatic niche axes were examined across the phylogeny. We found the Virentes to be monophyletic with three subclades, including a southwest clade, a southeastern US clade and a Central American/Cuban clade. Despite high leaf morphological variation within species and transpecific chloroplast haplotypes, RADseq and nuclear SSR data showed genetic coherence of species. We estimated a crown date for Virentes of 11 Ma and implicated the formation of the Sea of Cortés in a speciation event ~5 Ma. Tree height at maturity, associated with fire tolerance, differs among the sympatric species, while freezing tolerance appears to have diverged repeatedly across the tropical–temperate divide. Sympatric species thus show evidence of ecological niche differentiation but share climatic niches, while allopatric and parapatric species conserve ecological niches, but diverge in climatic niches. The mode of speciation and/or degree of co‐occurrence may thus influence which niche axis plants diverge along.     

Conservation status assessment of the Sharp‐shinned Hawk,an endangered insular raptor in Puerto Rico

Sharp‐shinned Hawks (Accipiter striatus) are forest raptors that are widely distributed in the Americas. A subspecies endemic to Puerto Rico (A. s. venator) is listed as endangered and restricted to mature and old secondary montane forests and shade coffee plantations. However, recent information about the population status and distribution of Puerto Rican Sharp‐shinned Hawks is lacking. We developed a spatial geographic distribution model for Sharp‐shinned Hawks in Puerto Rico from 33 locations collected during four breeding seasons (2013–2016) using biologically relevant landscape variables (aspect, canopy closure, elevation, rainfall, slope, and terrain roughness). Elevation accounted for 89.8% of the model fit and predicted that the greatest probability of occurrence of Sharp‐shinned Hawks in Puerto Rico (> 60%) was at elevations above 900 m. Based on our model, an estimated 56.1 km² of habitat exists in Puerto Rico with a high probability of occurrence. This total represents ~0.6% of the island's area. Public lands included 43.8% of habitat with high probability of occurrence (24.6 km²), 96% of which was located within four protected areas. Our results suggest that Sharp‐shinned Hawks are rare in Puerto Rico and restricted to the higher elevations of the Cordillera Central. Additional research is needed to identify and address ecological limiting factors, and recovery actions are needed to avoid the extinction of this endemic island raptor.     

Leaf size variations in a dominant desert shrub,Reaumuria soongarica,adapted to heterogeneous environments

The climate in arid Central Asia (ACA) has changed rapidly in recent decades, but the ecological consequences of this are far from clear. To predict the impacts of climate change on ecosystem functioning, greater attention should be given to the relationships between leaf functional traits and environmental heterogeneity. As a dominant constructive shrub widely distributed in ACA, Reaumuria soongarica provided us with an ideal model to understand how leaf functional traits of desert ecosystems responded to the heterogeneous environments of ACA. Here, to determine the influences of genetic and ecological factors, we characterized species‐wide variations in leaf traits among 30 wild populations of R. soongarica and 16 populations grown in a common garden. We found that the leaf length, width, and leaf length to width ratio (L/W) of the northern lineage were significantly larger than those of other genetic lineages, and principal component analysis based on the in situ environmental factors distinguished the northern lineage from the other lineages studied. With increasing latitude, leaf length, width, and L/W in the wild populations increased significantly. Leaf length and L/W were negatively correlated with altitude, and first increased and then decreased with increasing mean annual temperature (MAT) and mean annual precipitation (MAP). Stepwise regression analyses further indicated that leaf length variation was mainly affected by latitude. However, leaf width was uncorrelated with altitude, MAT, or MAP. The common garden trial showed that leaf width variation among the eastern populations was caused by both local adaptation and phenotypic plasticity. Our findings suggest that R. soongarica preferentially changes leaf length to adjust leaf size to cope with environmental change. We also reveal phenotypic evidence for ecological speciation of R. soongarica. These results will help us better understand and predict the consequences of climate change for desert ecosystem functioning.     

Trapped within the city: integrating demography,time since isolation and population‐specific traits to assess the genetic effects of urbanization

Urbanization is a severe form of habitat fragmentation that can cause many species to be locally extirpated and many others to become trapped and isolated within an urban matrix. The role of drift in reducing genetic diversity and increasing genetic differentiation is well recognized in urban populations. However, explicit incorporation and analysis of the demographic and temporal factors promoting drift in urban environments are poorly studied. Here, we genotyped 15 microsatellites in 320 fire salamanders from the historical city of Oviedo (Est. 8th century) to assess the effects of time since isolation, demographic history (historical effective population size; N_e) and patch size on genetic diversity, population structure and contemporary N_e. Our results indicate that urban populations of fire salamanders are highly differentiated, most likely due to the recent N_e declines, as calculated in coalescence analyses, concomitant with the urban development of Oviedo. However, urbanization only caused a small loss of genetic diversity. Regression modelling showed that patch size was positively associated with contemporary N_e, while we found only moderate support for the effects of demographic history when excluding populations with unresolved history. This highlights the interplay between different factors in determining current genetic diversity and structure. Overall, the results of our study on urban populations of fire salamanders provide some of the very first insights into the mechanisms affecting changes in genetic diversity and population differentiation via drift in urban environments, a crucial subject in a world where increasing urbanization is forecasted.     

Phenotypic diversification and island evolution of pipistrelle bats (Pipistrellus pipistrellus group) in the Mediterranean region inferred from geometric morphometrics and molecular phylogenetics

Aim The Mediterranean Basin is a centre of radiation for numerous species groups. To increase our understanding of the mechanisms underlying speciation and radiation events in this region, we assessed the phenotypic variability within the Pipistrellus pipistrellus–pygmaeus–hanaki species complex. Although bats form the second largest mammalian order, studies of insular evolution in this group are scarce. We approached this problem from a microevolutionary perspective and tested for the recurrence of the insular syndrome. Location The Mediterranean Basin, with a special focus on isolated populations from Corsica, the Maghreb, Cyprus, Cyrenaica and Crete. Methods Phenotypic variability was assessed by cranial morphometrics using the coordinates of 41 3D landmarks and associated geometric‐morphometric methods. We analysed 125 specimens representing all of the lineages in the species complex. Differences between taxa and between insular and continental populations in cranial size, shape, form and allometries were tested using analyses of variance and visualized using boxplots and canonical variate analysis. Relationships between molecular data from a previous study (cytochrome b sequences) and morphometric data were tested with co‐inertia analyses (RV test) and multivariate regressions. Results The three species were relatively well differentiated in cranial size and shape, and each species showed a significant amount of inter‐population variability. Comparisons of pairs of insular versus continental populations revealed heterogeneities in cranial patterns among island phenotypes, suggesting no recurrent insular syndrome. Molecular and phenotypic traits were correlated, except for molecular and lateral cranium shape. Main conclusions The Pipistrellus pipistrellus – pygmaeus – hanaki species complex exhibits phenotypic variability as a result of the fragmentation of its distribution (especially on islands), its phylogenetic and phylogeographic history and, most probably, other evolutionary factors that were not investigated in this study. We found no recurrent pattern of evolution on islands, indicating that site‐specific factors play a prevailing role on Mediterranean islands. The correlation between molecular and phenotypic data is incomplete, suggesting that factors other than phylogenetic relationships, potentially connected with feeding ecology, have played a role in shaping cranial morphology in this species complex.