Phylogeography and species limits in the green woodpecker complex (Aves: Picidae): multiple Pleistocene refugia and range expansion across Europe and the Near East

Chromosome-level genome assembly of the European green woodpecker Picus viridis

The European green woodpecker, Picus viridis, is a widely distributed species found in the Western Palearctic region. Here, we assembled a highly contiguous genome assembly for this species using a combination of short- and long-read sequencing and scaffolded with chromatin conformation capture (Hi-C). The final genome assembly was 1.28 Gb and features a scaffold N50 of 37 Mb and a scaffold L50 of 39.165 Mb. The assembly incorporates 89.4% of the genes identified in birds in OrthoDB. Gene and repetitive content annotation on the assembly detected 15,805 genes and a ∼30.1% occurrence of repetitive elements, respectively. Analysis of synteny demonstrates the fragmented nature of the P. viridis genome when compared to the chicken (Gallus gallus). The assembly and annotations produced in this study will certainly help for further research into the genomics of P. viridis and the comparative evolution of woodpeckers. Five historical and seven contemporary samples have been resequenced and may give insights on the population history of this species.

Mapping endemic freshwater fish richness to identify high-priority areas for conservation: An ecoregion approach

Freshwater ecosystems are experiencing accelerating global biodiversity loss. Thus, knowing where these unique ecosystems' species richness reaches a peak can facilitate their conservation planning. By hosting more than 290 freshwater fishes, Iran is a major freshwater fish hotspot in the Middle East. Considering the accelerating rate of biodiversity loss, there is an urgent need to identify species-rich areas and understand the mechanisms driving biodiversity distribution. In this study, we gathered distribution records of all endemic freshwater fishes of Iran (85 species) to develop their richness map and determine the most critical drivers of their richness patterns from an ecoregion approach. We performed a generalized linear model (GLM) with quasi-Poisson distribution to identify contemporary and historical determinants of endemic freshwater fish richness. We also quantified endemic fish similarity among the 15 freshwater ecoregions of Iran. Results showed that endemic freshwater fish richness is highest in the Zagros Mountains while a moderate level of richness was observed between Zagros and Alborz Mountains. High, moderate, and low richness of endemic freshwater fish match with Upper Tigris & Euphrates, Namak, and Kavir & Lut Deserts ecoregions respectively. Kura - South Caspian Drainages and Caspian Highlands were the most similar ecoregions and Orumiyeh was the most unique ecoregion according to endemic fish presence. Precipitation and precipitation change velocity since the Last Glacial Maximum were the most important predictors of endemic freshwater fish richness. Areas identified to have the highest species richness have high priority for the conservation of freshwater fish in Iran, therefore, should be considered in future protected areas development.

Modeling habitat suitability of bats to identify high priority areas for field monitoring and conservation

Bats provide important ecosystem services but face severe threats due to land and climate changes. Although bats are an important component of mammal diversity in Iran, the ecology of many species remains virtually unstudied in the country. Here we applied the maximum entropy approach to model habitat suitability of bat species in Iran, identify the most important variables for their distribution, predict high priority areas for field monitoring and conservation, and estimate the coverage of the bat species' suitable habitats by the existing protected areas. We created a richness map for the twelve species to identify high priority areas for field monitoring and conservation. The results of species distribution modeling showed that Pipistrellus kuhlii (828,977.2 km²) and Miniopterus pallidus (646,581.9 km²) had the largest distribution ranges and Rhinopoma microphyllum (211,202.7 km²) and Rousettus aegyptiacus (218,278.6 km²) had the smallest distribution ranges in Iran. By averaging the importance of each ecological variable across the 12 species, we found that distance to forests (with a negative association) is the most important ecological driver of bat distribution in Iran. The Zagros Mountains were identified as a hotspot of bats based on the distribution of the 12 species. Our findings showed that small proportions of each species suitable habitats were covered by protected areas and protected suitable habitats varied from 3.2%for Pipistrellus kuhlii to15.9% for Tadarida teniotis. This study highlights the importance of forests for bat conservation showing that forest conservation is a high priority in the country. Areas which have the highest richness should be prioritized for field monitoring and conservation.

Pleistocene allopatric differentiation followed by recent range expansion explains the distribution and molecular diversity of two congeneric crustacean species in the Palaearctic

Pleistocene glaciations had a tremendous impact on the biota across the Palaearctic, resulting in strong phylogeographic signals of range contraction and rapid postglacial recolonization of the deglaciated areas. Here, we explore the diversity patterns and history of two sibling species of passively dispersing taxa typical of temporary ponds, fairy shrimps (Anostraca). We combine mitochondrial (COI) and nuclear (ITS2 and 18S) markers to conduct a range-wide phylogeographic study including 56 populations of Branchinecta ferox and Branchinecta orientalis in the Palaearctic. Specifically, we investigate whether their largely overlapping ranges in Europe resulted from allopatric differentiation in separate glacial refugia followed by a secondary contact and reconstruct their postglacial recolonization from the inhabited refugia. Our results suggest the existence of distinct refugia for the two species, with genetic divergence among intraspecific lineages consistent with late Pleistocene glacial cycles. While B. ferox lineages originated from Mediterranean refugia, the origin of B. orientalis lineages was possibly located on the Pannonian Plain. We showed that most dispersal events predominantly happened within 100 km, coupled with several recent long-distance events (> 1000 km). Hence the regional habitat density of suitable habitats in Central Europe is possibly a key to the co-existence of the two species. Overall, our study illustrates how isolation in combination with stochastic effects linked to glacial periods are important drivers of the allopatric differentiation of Palaearctic taxa.

Common patterns in the molecular phylogeography of western palearctic birds: a comprehensive review

A plethora of studies have offered crucial insights in the phylogeographic status of Western Palearctic bird species. However, an overview integrating all this information and analyzing the combined results is still missing. In this study, we compiled all published peer-reviewed and grey literature available on the phylogeography of Western Palearctic bird species. Our literature review indicates a total number of 198 studies, with the overwhelming majority published as journal articles (n = 186). In total, these literature items offer information on 145 bird species. 85 of these species are characterized by low genetic differentiation, 46 species indicate genetic variation but no geographic structuring i.e. panmixia, while 14 species show geographically distinct lineages and haplotypes. Majority of bird species inhabiting the Western Palearctic display genetic admixture. The glaciation cycles in the past few million years were pivotal factors in shaping this situation: during warm periods many species expanded their distribution range to the north over wide areas of Eurasia; whereas, during ice ages most areas were no longer suitable and species retreated to refugia, where lineages mixed.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10336-021-01893-x.

Reptile species richness associated to ecological and historical variables in Iran

Spatial gradients of species richness can be shaped by the interplay between historical and ecological factors. They might interact in particularly complex ways in heterogeneous mountainous landscapes with strong climatic and geological contrasts. We mapped the distribution of 171 lizard species to investigate species richness patterns for all species (171), diurnal species (101), and nocturnal species (70) separately. We related species richness with the historical (past climate change, mountain uplifting) and ecological variables (climate, topography and vegetation). We found that assemblages in the Western Zagros Mountains, north eastern and north western parts of Central Iranian Plateau have the highest number of lizard species. Among the investigated variables, annual mean temperature explained the largest variance for all species (10%) and nocturnal species (31%). For diurnal species, temperature change velocity shows strongest explained variance in observed richness pattern (26%). Together, our results reveal that areas with annual temperature of 15–20 °C, which receive 400–600 mm precipitation and experienced moderate level of climate change since the Last Glacial Maximum (LGM) have highest number of species. Documented patterns of our study provide a baseline for understanding the potential effect of ongoing climate change on lizard diversity in Iran.

Finding complexity in complexes: Assessing the causes of mitonuclear discordance in a problematic species complex of Mesoamerican toads

Mitonuclear discordance is a frequently encountered pattern in phylogeographic studies and occurs when mitochondrial and nuclear DNA display conflicting signals. Discordance among these genetic markers can be caused by several factors including confounded taxonomies, gene flow, and incomplete lineage sorting. In this study, we present a strong case of mitonuclear discordance in a species complex of toads (Bufonidae: Incilius coccifer complex) found in the Chortís Block of Central America. To determine the cause of mitonuclear discordance in this complex, we used spatially explicit genetic data to test species limits and relationships, characterize demographic history, and quantify gene flow. We found extensive mitonuclear discordance among the three recognized species within this group, especially in populations within the Chortís Highlands of Honduras. Our data reveal nuclear introgression within the Chortís Highlands populations that was most probably driven by cyclical range expansions due to climatic fluctuations. Though we determined introgression occurred within the nuclear genome, our data suggest that it is not the key factor in driving mitonuclear discordance in the entire species complex. Rather, due to a lack of discernible geographic pattern between mitochondrial and nuclear DNA, as well as a relatively recent divergence time of this complex, we concluded that mitonuclear discordance has been caused by incomplete lineage sorting. Our study provides a framework to test sources of mitonuclear discordance and highlights the importance of using multiple marker types to test species boundaries in cryptic species.

Genetic characterization of the Turkish gray hamster (Cricetulus migratorius) based on mitochondrial cytochrome b and 12S rRNA sequences

Although genetic diversity and phylogenetic status of the gray hamster (Cricetulus migratorius) have been investigated from different regions in previous studies, genetic data on this species from Turkey are still lacking, since previous data have been based on a limited number of gray hamsters sampled across the Anatolian part of Turkey. The aim of this study was to determine the genetic diversity of the Anatolian population and to reveal the phylogenetic relationships among the Anatolian population and conspecific populations of the gray hamster. The complete and partial fragments of mitochondrial Cyt b and 12S rRNA from the 20 Turkish samples were amplified and sequenced. Ten 12S rRNA (901 bp) and 15 Cyt b (1140 bp) haplotypes found in this work were not previously reported. Based on Bayesian, Maximum Likelihood, Neighbour-Joining and Median-Joining network analyses by using mitochondrial data under the name Cricetulus, the results of phylogenetic and network analyses indicated that there was a deep separation among the distinct lineages within the genus Cricetulus. When considering the position of the Turkish haplotypes in median joining network, the Anatolian part of Turkey may have hosted a source population of the gray hamster for expansion to adjacent regions in the past period. Additionally, the Anatolian population of gray hamster had relatively high haplotype diversity and the present study propounded the importance of data obtained from the Anatolian population of gray hamster to reveal the phylogenetic relationships among conspecific populations of the gray hamster.

Range instability leads to cytonuclear discordance in a morphologically cryptic ground squirrel species complex

The processes responsible for cytonuclear discordance frequently remain unclear. Here, we employed an exon capture data set and demographic methods to test hypotheses generated by species distribution models to examine how contrasting histories of range stability vs. fluctuation have caused cytonuclear concordance and discordance in ground squirrel lineages from the Otospermophilus beecheyi species complex. Previous studies in O. beecheyi revealed three morphologically cryptic and highly divergent mitochondrial DNA lineages (named the Northern, Central and Southern lineages based on geography) with only the Northern lineage exhibiting concordant divergence for nuclear genes. Here, we showed that these mtDNA lineages likely formed in allopatry during the Pleistocene, but responded differentially to climatic changes that occurred since the last interglacial (~120,000 years ago). We find that the Northern lineage maintained a stable range throughout this period, correlating with genetic distinctiveness among all genetic markers and low migration rates with the other lineages. In contrast, our results suggested that the Southern lineage expanded from Baja California Sur during the Late Pleistocene to overlap and potentially swamp a contracting Central lineage. High rates of intraspecific gene flow between Southern lineage individuals among expansion origin and expansion edge populations largely eroded Central ancestry from autosomal markers. However, male-biased dispersal in this system preserved signals of this past hybridization and introgression event in matrilineal-biased X-chromosome and mtDNA markers. Our results highlight the importance of range stability in maintaining the persistence of phylogeographic lineages, whereas unstable range dynamics can increase the tendency for lineages to merge upon secondary contact.

Wild gazelles of the southern Levant: genetic profiling defines new conservation priorities

The mountain gazelle (Gazella gazelle), Dorcas gazelle (Gazella Dorcas) and acacia gazelle (Gazella arabica acacia) were historically abundant in the southern Levant, and more specifically in Israel. Anthropogenic and natural changes have caused a rapid decline in gazelle populations, raising concerns about their conservation status and future survival. The genetic profile of 111 wild gazelles from Israel was determined based on three regions of mitochondrial DNA (control region, Cytochrome b and 12S ribosomal RNA) and nine nuclear microsatellite markers. Genetic analysis of the mountain gazelle population, the largest known population of this rare species, revealed adequate diversity levels and gene flow between subpopulations. Nevertheless, ongoing habitat degradation and other human effects, such as poaching, suggest the need for drastic measures to prevent species extinction. Dorcas gazelles in Israel displayed inbreeding within subpopulations while still maintaining considerable genetic diversity overall. This stable population, represented by a distinctive genetic profile, is fragmented and isolated from its relatives in neighboring localities. Based on the genetic profile of a newly sampled subpopulation in Israel, we provide an alternative hypothesis for the historic dispersal of Dorcas gazelle, from the Southern Levant to northern Africa. The small acacia gazelle population was closest to gazelles from the Farasan Islands of Saudi Arabia, based on mitochondrial markers. The two populations did not share haplotypes, suggesting that these two populations may be the last remnant wild gazelles of this species worldwide. Only a dozen acacia gazelles survive in Israel, and urgent steps are needed to ensure the survival of this genetically distinctive lineage. The genetic assessments of our study recognize new conservation priorities for each gazelle species in the Southern Levant.