Further Evidence for Cryptic North-Western Refugia in Europe? Mitochondrial Phylogeography of the Sibling SpeciesPipistrellus pipistrellusandPipistrellus pygmaeus

Geographic Genetic Structure of Alectoris chukar in Türkiye: Post-LGM-Induced Hybridization and Human-Mediated Contaminations

Türkiye is considered an important evolutionary area for Chukar partridge (Alectoris chukar), since it is both a potential ancestral area and a diversification center for the species. Using 2 mitochondrial (Cty-b and D-loop) and 13 polymorphic microsatellite markers, we investigated the geographic genetic structure of A. chukar populations to determine how past climatic fluctuations and human activities have shaped the gene pool of this species in Türkiye. Our results indicate, firstly, that only A. chukar of the genus Alectoris is present in Türkiye (Anatolia and Thrace), with no natural or artificial gene flow from congenerics. Secondly, the geographic genetic structure of the species in Türkiye has been shaped by topographic heterogeneity, Pleistocene climatic fluctuations, and artificial transport by humans. Third, there appears to be three genetic clusters: Thracian, Eastern, and Western. Fourth, the post-LGM demographic expansion of the Eastern and Western populations has formed a hybrid zone in Central Anatolia (~8 kyBP). Fifth, the rate of China clade-B contamination in Türkiye is about 8% in mtDNA and about 12% in nuDNA, with the Southeastern Anatolian population having the highest contamination. Sixth, the Thracian population was the most genetically distinct, with the lowest genetic diversity and highest level of inbreeding and no China clad-B contamination. These results can contribute to the conservation regarding A. chukar populations, especially the Thracian population.

The genome sequence of the common pipistrelle, Pipistrellus pipistrellus Schreber 1774

We present a genome assembly from an individual female Pipistrellus pipistrellus (the common pipistrelle; Chordata; Mammalia; Chiroptera; Vespertilionidae). The genome sequence is 1.76 gigabases in span. The majority of the assembly is scaffolded into 21 chromosomal pseudomolecules, with the X sex chromosome assembled.

Red deer in Iberia: Molecular ecological studies in a southern refugium and inferences on European postglacial colonization history

The red deer (Cervus elaphus) is a widespread wild ungulate in Europe that has suffered strong anthropogenic impacts over their distribution during the last centuries, but also at the present time, due its economic importance as a game species. Here we focus on the evolutionary history of the red deer in Iberia, one of the three main southern refugial areas for temperate species in Europe, and addressed the hypothesis of a cryptic refugia at higher latitudes during the Last Glacial Maximum (LGM). A total of 911 individuals were sampled, genotyped for 34 microsatellites specifically developed for red deer and sequenced for a fragment of 670 bp of the mitochondrial (mtDNA) D-loop. The results were combined with published mtDNA sequences, and integrated with species distribution models and historical European paleo-distribution data, in order to further examine the alternative glacial refugial models and the influence of cryptic refugia on European postglacial colonization history. Clear genetic differentiation between Iberian and European contemporary populations was observed at nuclear and mtDNA levels, despite the mtDNA haplotypes central to the phylogenetic network are present across western Europe (including Iberia) suggesting a panmictic population in the past. Species distribution models, fossil records and genetic data support a timing of divergence between Iberian and European populations that overlap with the LGM. A notable population structure was also found within the Iberian Peninsula, although several populations displayed high levels of admixture as a consequence of recent red deer translocations. Five D-loop sub-lineages were found in Iberia that belong to the Western European mtDNA lineage, while there were four main clusters based on analysis of nuclear markers. Regarding glacial refugial models, our findings provide detailed support for the hypothesis that red deer may have persisted in cryptic northern refugia in western Europe during the LGM, most likely in southern France, southern Ireland, or in a region between them (continental shelf), and these regions were the source of individuals during the European re-colonization. This evidence heightens the importance of conserving the high mitochondrial and nuclear diversity currently observed in Iberian populations.

Historical explanation of genetic variation in the Mediterranean horseshoe bat Rhinolophus euryale (Chiroptera: Rhinolophidae) inferred from mitochondrial cytochrome-b and D-loop genes in Iran

Molecular phylogeography and species distribution modelling (SDM) suggest that late Quaternary glacial cycles have portrayed a significant role in structuring current population genetic structure and diversity. Based on phylogenetic relationships using Bayesian inference and maximum likelihood of 535 bp mtDNA (D-loop) and 745 bp mtDNA (Cytb) in 62 individuals of the Mediterranean Horseshoe Bat, Rhinolophus euryale, from 13 different localities in Iran we identified two subspecific populations with differing population genetic structure distributed in southern Zagros Mts. and northern Elburz Mts. Analysis of molecular variance (AMOVA) obtained from D-loop sequences indicates that 21.18% of sequence variation is distributed among populations and 10.84% within them. Moreover, a degree of genetic subdivision, mainly attributable to the existence of significant variance among the two regions is shown (θCT = 0.68, p = .005). The positive and significant correlation between geographic and genetic distances (R² = 0.28, r = 0.529, p = .000) is obtained following controlling for environmental distance. Spatial distribution of haplotypes indicates that marginal population of the species in southern part of the species range have occupied this section as a glacial refugia. However, this genetic variation, in conjunction with results of the SDM shows a massive postglacial range expansion for R. euryale towards higher latitudes in Iran.

Extra-Mediterranean glacial refuges in barred and common grass snakes (Natrix helvetica, N. natrix)

Extra-Mediterranean glacial refugia of thermophilic biota, in particular in northern latitudes, are controversial. In the present study we provide genetic evidence for extra-Mediterranean refugia in two species of grass snake. The refuge of a widely distributed western European lineage of the barred grass snake (Natrix helvetica) was most likely located in southern France, outside the classical refuges in the southern European peninsulas. One genetic lineage of the common grass snake (N. natrix), distributed in Scandinavia, Central Europe and the Balkan Peninsula, had two distinct glacial refuges. We show that one was located in the southern Balkan Peninsula. However, Central Europe and Scandinavia were not colonized from there, but from a second refuge in Central Europe. This refuge was located in between the northern ice sheet and the Alpine glaciers of the last glaciation and most likely in a permafrost region. Another co-distributed genetic lineage of N. natrix, now massively hybridizing with the aforementioned lineage, survived the last glaciation in a structured refuge in the southern Balkan Peninsula, according to the idea of 'refugia-within-refugia'. It reached Central Europe only very recently. This study reports for the first time the glacial survival of a thermophilic egg-laying reptile species in Central Europe.

Nuclear introns outperform mitochondrial DNA in inter-specific phylogenetic reconstruction: Lessons from horseshoe bats (Rhinolophidae: Chiroptera)

Despite many studies illustrating the perils of utilising mitochondrial DNA in phylogenetic studies, it remains one of the most widely used genetic markers for this purpose. Over the last decade, nuclear introns have been proposed as alternative markers for phylogenetic reconstruction. However, the resolution capabilities of mtDNA and nuclear introns have rarely been quantified and compared. In the current study we generated a novel ∼5kb dataset comprising six nuclear introns and a mtDNA fragment. We assessed the relative resolution capabilities of the six intronic fragments with respect to each other, when used in various combinations together, and when compared to the traditionally used mtDNA. We focused on a major clade in the horseshoe bat family (Afro-Palaearctic clade; Rhinolophidae) as our case study. This old, widely distributed and speciose group contains a high level of conserved morphology. This morphological stasis renders the reconstruction of the phylogeny of this group with traditional morphological characters complex. We sampled multiple individuals per species to represent their geographic distributions as best as possible (122 individuals, 24 species, 68 localities). We reconstructed the species phylogeny using several complementary methods (partitioned Maximum Likelihood and Bayesian and Bayesian multispecies-coalescent) and made inferences based on consensus across these methods. We computed pairwise comparisons based on Robinson-Foulds tree distance metric between all Bayesian topologies generated (27,000) for every gene(s) and visualised the tree space using multidimensional scaling (MDS) plots. Using our supported species phylogeny we estimated the ancestral state of key traits of interest within this group, e.g. echolocation peak frequency which has been implicated in speciation. Our results revealed many potential cryptic species within this group, even in taxa where this was not suspected a priori and also found evidence for mtDNA introgression. We demonstrated that by using just two introns one can recover a better supported species tree than when using the mtDNA alone, despite the shorter overall length of the combined introns. Additionally, when combining any single intron with mtDNA, we showed that the result is highly similar to the mtDNA gene tree and far from the true species tree and therefore this approach should be avoided. We caution against the indiscriminate use of mtDNA in phylogenetic studies and advocate for pilot studies to select nuclear introns. The selection of marker type and number is a crucial step that is best based on critical examination of preliminary or previously published data. Based on our findings and previous publications, we recommend the following markers to recover phylogenetic relationships between recently diverged taxa (<20 My) in bats and other mammals: ACOX2, COPS7A, BGN, ROGDI and STAT5A.

Two Mitochondrial Barcodes for one Biological Species: The Case of European Kuhl's Pipistrelles (Chiroptera)

The Kuhl's pipistrelle (Pipistrellus kuhlii) is a Western Palaearctic species of bat that exhibits several deeply divergent mitochondrial lineages across its range. These lineages could represent cryptic species or merely ancient polymorphism, but no nuclear markers have been studied so far to properly assess the taxonomic status of these lineages. We examined here two lineages occurring in Western Europe, and used both mitochondrial and nuclear markers to measure degrees of genetic isolation between bats carrying them. The sampling focused on an area of strict lineage sympatry in Switzerland but also included bats from further south, in North Africa. All individuals were barcoded for the COI gene to identify their mitochondrial lineages and five highly polymorphic microsatellite loci were used to cluster them according to their nuclear genotypes. Despite this low number of nuclear markers, all North African nuclear genotypes were grouped in a highly distinct subpopulation when compared with European samples sharing the same mitochondrial barcodes. The reverse situation prevailed in Switzerland where bats carrying distinct barcodes had similar nuclear genotypes. There was a weak east/west nuclear structure of populations, but this was independent of mitochondrial lineages as bats carrying either variant were completely admixed. Thus, the divergent mitochondrial barcodes present in Western Europe do not represent cryptic species, but are part of a single biological species. We argue that these distinct barcodes evolved in allopatry and came recently into secondary contact in an area of admixture north of the Alps. Historical records from this area and molecular dating support such a recent bipolar spatial expansion. These results also highlight the need for using appropriate markers before claiming the existence of cryptic species based on highly divergent barcodes.