Whole-genome SNP genotyping unveils ancestral and recent introgression in wild and domestic goats.

After the domestication of goats around 10,000 years before the present (BP), humans transported goats far beyond the range of their wild ancestor, the bezoar goat. This brought domestic goats into contact with many wild goat species such as ibex and markhor, enabling introgression between domestic and wild goats. To investigate this, while shedding light on the taxonomic status of wild and domestic goats, we analysed genome-wide SNP data of 613 specimens from 14 taxonomic units, including Capra hircus, C. pyrenaica, C. ibex (from Switzerland, Austria, Germany and Slovenia), C. aegagrus aegagrus, C. a. cretica, C. h. dorcas, C. caucasica caucasica, C. c. severtzovi, C. c. cylindricornis, C. falconeri, C. sibirica sibirica, C. s. alaiana and C. nubiana, as well as Oreamnos americanus (mountain goat) as an outgroup. To trace gene flow between domestic and wild goats, we integrated genotype data of local goat breeds from the Alps as well as from countries such as Spain, Greece, Türkiye, Egypt, Sudan, Iran, Russia (Caucasus and Altai) and Pakistan. Our phylogenetic analyses displayed a clear separation between bezoar-type and ibex-type clades with wild goats from the Greek islands of Crete and Youra clustered within domestic goats, confirming their feral origin. Our analyses also revealed gene flow between the lineages of Caucasian tur and domestic goats that most likely occurred before or during early domestication. Within the clade of domestic goats, analyses inferred gene flow between African and Iberian goats. The detected events of introgression were consistent with previous reports and offered interesting insights into the historical relationships among domestic and wild goats.

Cryptic diversity and phylogeographic patterns of Mediodactylus species in the Eastern Mediterranean region.

Cryptic diversity poses a great obstacle in our attempts to assess the current biodiversity crisis and may hamper conservation efforts. The gekkonid genus Mediodactylus, a well-known case of hidden species and genetic diversity, has been taxonomically reclassified several times during the last decade. Focusing on the Mediterranean populations, a recent study within the M. kotschyi species complex using classic mtDNA/nuDNA markers suggested the existence of five distinct species, some being endemic and some possibly threatened, yet their relationships have not been fully resolved. Here, we generated genome-wide SNPs (using ddRADseq) and applied molecular species delimitation approaches and population genomic analyses to further disentangle these relationships. Τhe most extensive nuclear dataset, so far, encompassing 2,360 loci and âˆ¼ 699,000 bp from across the genome of Mediodactylus gecko, enabled us to resolve previously obscure phylogenetic relationships among the five, recently elevated, Mediodactylus species and to support the hypothesis that the taxon includes several new, undescribed species. Population genomic analyses within each of the proposed species showed strong genetic structure and high levels of genetic differentiation among populations.

Population Genomics of Wall Lizards Reflects the Dynamic History of the Mediterranean Basin.

The Mediterranean Basin has experienced extensive change in geology and climate over the past six million years. Yet, the relative importance of key geological events for the distribution and genetic structure of the Mediterranean fauna remains poorly understood. Here, we use population genomic and phylogenomic analyses to establish the evolutionary history and genetic structure of common wall lizards (Podarcis muralis). This species is particularly informative because, in contrast to other Mediterranean lizards, it is widespread across the Iberian, Italian, and Balkan Peninsulas, and in extra-Mediterranean regions. We found strong support for six major lineages within P. muralis, which were largely discordant with the phylogenetic relationship of mitochondrial DNA. The most recent common ancestor of extant P. muralis was likely distributed in the Italian Peninsula, and experienced an "Out-of-Italy" expansion following the Messinian salinity crisis (∼5 Mya), resulting in the differentiation into the extant lineages on the Iberian, Italian, and Balkan Peninsulas. Introgression analysis revealed that both inter- and intraspecific gene flows have been pervasive throughout the evolutionary history of P. muralis. For example, the Southern Italy lineage has a hybrid origin, formed through admixture between the Central Italy lineage and an ancient lineage that was the sister to all other P. muralis. More recent genetic differentiation is associated with the onset of the Quaternary glaciations, which influenced population dynamics and genetic diversity of contemporary lineages. These results demonstrate the pervasive role of Mediterranean geology and climate for the evolutionary history and population genetic structure of extant species.

High-throughput degraded DNA sequencing of subfossil shells of a critically endangered stenoendemic land snail in the Aegean.

High-throughput sequencing has enabled the comprehensive genetic exploration of biological diversity, especially by using natural history collections to study hard-to-find, threatened or even extinct-in-the-wild taxa. Mollusk shells are under-exploited as a source for DNA-based approaches, despite their apparent advantages in the field of conservation genetics. More particularly, degraded DNA techniques combined with high-throughput sequencing have never been used to gain insights about the DNA preservation in land snail subfossil or historical shells. Here, we applied degraded DNA analysis on two historical shells of Levantina rechingeri, a stenoendemic Critically Endangered species that has never been found alive, in order to explore the patterns of DNA preservation on land snail shells originating from the eastern Mediterranean, as well as to infer its molecular phylogenetic placement. Our results showed that centuries to decades-old DNA from an empty shell collected in an Aegean island exhibits characteristic post-mortem damage patterns similar to those observed in ancient DNA from eastern Mediterranean terrestrial animals, setting a precedent for future museomics studies on taxa distributed in areas with similar climate. Finally, genome skimming of the empty shell allowed high coverage of multiple nuclear and mitochondrial loci, enabling the phylogenetic placement of the focal taxon, the re-evaluation of its taxonomic classification, and the revealing of a new Aegean land snail lineage, Aristena genus novum. This approach is a non-invasive way to sample DNA from threatened land snail species and suitable to study the evolutionary history of taxa with cryptic ecology, stenoendemics, or extinct-in-the-wild, as well as old museum specimens.

Diversification within an oceanic Mediterranean island: Insights from a terrestrial isopod.

Understanding intra-island patterns of evolutionary divergence, including cases of cryptic diversity, is a crucial step towards deciphering speciation processes. Cyprus is an oceanic island isolated for at least 5.3 Mya from surrounding continental regions, while it remains unclear whether it was ever connected to the mainland, even during the Messinian Salinity Crisis. The terrestrial isopod species Armadillo officinalis, that is widespread across the Mediterranean, offers the opportunity to explore intra-island divergence patterns that might exhibit geographical structure related also to the region's known paleogeography. Genome-wide ddRADseq, as well as Sanger sequencing for four mitochondrial and three nuclear loci data were generated for this purpose. In total, 71 populations from Cyprus, neighbouring continental sites, i.e., Israel, Lebanon and Turkey, and other Mediterranean regions, i.e. Greece, Italy, and Tunisia, were included in the analysis. Phylogenetic reconstructions and population structure analyses support the existence of at least six genetically discrete groups across the study area. Five of these distinct genetic clades occur on Cyprus, four of which are endemic to the island and one is widely distributed along the circum-Mediterranean countries. The sixth clade is distributed in Israel. The closest evolutionary relationship of endemic Cypriot populations is with those from Israel, while the evolutionary clade that is present in countries all around the Mediterranean is very shallow. Cladochronological analyses date the origin of the species on the island at ∼6 Mya. Estimated f4 and D statistics as well as FST values indicate the genetic isolation between the populations sampled from Cyprus and surrounding continental areas, while there is evident gene flow among populations within the island. Species delimitation and population genetic metrics support the existence of three distinct taxonomic units across the study area, two of which occur on the island and correspond to the endemic clade and the widespread circum-Mediterranean one, respectively, while the third corresponds to Israel's clade. The islands' paleogeographic history and recent human activities seem to have shaped current patterns of genetic diversity in this group of species.

A new lineage of Galapagos giant tortoises identified from museum samples.

The Galapagos Archipelago is recognized as a natural laboratory for studying evolutionary processes. San Cristóbal was one of the first islands colonized by tortoises, which radiated from there across the archipelago to inhabit 10 islands. Here, we sequenced the mitochondrial control region from six historical giant tortoises from San Cristóbal (five long deceased individuals found in a cave and one found alive during an expedition in 1906) and discovered that the five from the cave are from a clade that is distinct among known Galapagos giant tortoises but closely related to the species from Española and Pinta Islands. The haplotype of the individual collected alive in 1906 is in the same clade as the haplotype in the contemporary population. To search for traces of a second lineage in the contemporary population on San Cristóbal, we closely examined the population by sequencing the mitochondrial control region for 129 individuals and genotyping 70 of these for both 21 microsatellite loci and >12,000 genome-wide single nucleotide polymorphisms [SNPs]. Only a single mitochondrial haplotype was found, with no evidence to suggest substructure based on the nuclear markers. Given the geographic and temporal proximity of the two deeply divergent mitochondrial lineages in the historical samples, they were likely sympatric, raising the possibility that the lineages coexisted. Without the museum samples, this important discovery of an additional lineage of Galapagos giant tortoise would not have been possible, underscoring the value of such collections and providing insights into the early evolution of this iconic radiation.

The wall lizards of the Balkan peninsula: Tackling questions at the interface of phylogenomics and population genomics.

Wall lizards of the genus Podarcis (Sauria, Lacertidae) are the predominant reptile group in southern Europe, including 24 recognized species. Mitochondrial DNA data have shown that, with the exception of P. muralis, the Podarcis species distributed in the Balkan peninsula form a species group that is further sub-divided into two subgroups: the one of "P. tauricus" consisting of P. tauricus, P. milensis, P. gaigeae, and P. melisellensis, and the other of "P. erhardii" comprising P. erhardii, P. levendis, P. cretensis, and P. peloponnesiacus. In an attempt to explore the Balkan Podarcis phylogenomic relationships, assess the levels of genetic structure and to re-evaluate the number of extant species, we employed phylogenomic and admixture approaches on ddRADseq (double digested Restriction site Associated DNA sequencing) genomic data. With this efficient Next Generation Sequencing approach, we were able to obtain a large number of genomic loci randomly distributed throughout the genome and use them to resolve the previously obscure phylogenetic relationships among the different Podarcis species distributed in the Balkans. The obtained phylogenomic relationships support the monophyly of both aforementioned subgroups and revealed several divergent lineages within each subgroup, stressing the need for taxonomic re-evaluation of Podarcis' species in Balkans. The phylogenomic trees and the species delimitation analyses confirmed all recently recognized species (P. levendis, P. cretensis, and P. ionicus) and showed the presence of at least two more species, one in P. erhardii and the other in P. peloponnesiacus.

Multigene phylogeny, phylogeography and population structure of Podarcis cretensis species group in south Balkans.

The evolutionary history of taxa with limited overseas dispersal abilities is considered to be majorly influenced by vicariant events constituting them as model organisms for the interpretation of evolutionary processes. An excellent candidate are the wall lizards of the genus Podarcis exhibiting an impressive level of genetic and morphological diversification and harboring several cases of recently discovered cryptic diversity. In this study, we investigated the effect of palaeogeographic events on the wall lizards' biodiversity patterns in the Aegean (Greece) as well as the evolutionary processes that acted both in space and time. To accomplish that we studied a group of three endemic Podarcis species (i.e., P. cretensis, P. levendis, and P. peloponnesiacus) both at the intra and interspecific levels employing mitochondrial and nuclear DNA sequence data as well as microsatellites. Furthermore, presence information coupled with bioclimatic data (i.e., species distribution modeling, and niche similarity analyses) shed light on the necessary ecological factors for the species' occurrence. These approaches revealed yet another case of cryptic diversity for this group of lizards, with the existence of two slightly overlapping lineages within P. peloponnesiacus and highly structured populations within P. cretensis. Species diversification occurred during the Pliocene with P. peloponnesiacus divergence into the two lineages dating back to 1.86 Mya. Furthermore, temperature and precipitation related environmental parameters were the most important ones regarding the current distribution of the studied species. Based on the results, we propose a more detailed phylogeographic scenario where both the paleogeography of the area and several environmental parameters have shaped the genetic diversity and the current distribution pattern of this species group.

Evidence of introgressive hybridization between Stenella coeruleoalba and Delphinus delphis in the Greek Seas.

Natural interspecific hybridization might be more important for the evolutionary history and speciation of animals than previously thought, considering several demographic and life history traits as well as habitat disturbance as factors that promote it. In this aspect, cetaceans comprise an interesting case in which the occurrence of sympatric species in mixed associations provides excellent opportunities for interspecific sexual interaction and the potential for hybridization. Here, we present evidence of natural hybridization for two cetacean species commonly occurring in the Greek Seas (Stenella coeruleoalba and Delphinus delphis), which naturally overlap in the Gulf of Corinth by analyzing highly resolving microsatellite DNA markers and mitochondrial DNA sequences in skin samples from 45 individuals of S. coeruleoalba, 12 D. delphis and three intermediate morphs. Employing several phylogenetic and population genetic approaches, we found 15 individuals that are potential hybrids including the three intermediate morphs, verifying the occurrence of natural hybridization between species of different genera. Their hybrids are fertile and able to reproduce not only with the other hybrids but also with each of the two-parental species. However, current evidence does not allow firm conclusions whether hybridization might constitute a step towards the generation of a new species and/or the swan song of an already existing species (i.e., D. delphis). Given that the focal species form mixed pods in several areas of Mediterranean, this study is an excellent opportunity to understand the mechanisms leading to hybridization in the context of gene flow and urges for the evaluation of the genetic status of common dolphins in the Mediterranean.

Multilocus phylogeny and coalescent species delimitation in Kotschy's gecko, Mediodactylus kotschyi: Hidden diversity and cryptic species.

Kotschy's Gecko, Mediodactylus kotschyi, is a small gecko native to southeastern Europe and the Levant. It displays great morphological variation with a large number of morphologically recognized subspecies. However, it has been suggested that it constitutes a species complex of several yet unrecognized species. In this study, we used multilocus sequence data (three mitochondrial and three nuclear gene fragments) to estimate the phylogenetic relationships of 174 specimens from 129 sampling localities, covering a substantial part of the distribution range of the species. Our results revealed high genetic diversity of M. kotschyi populations and contributed to our knowledge about the phylogenetic relationships and the estimation of the divergence times between them. Diversification within M. kotschyi began approximately 15 million years ago (Mya) in the Middle Miocene, whereas the diversification within most of the major clades have been occurred in the last 5 Mya. Species delimitation analysis suggests there exists five species within the complex, and we propose to tentatively recognize the following taxa as full species: M. kotschyi (mainland Balkans, most of Aegean islands, and Italy), M. orientalis (Levant, Cyprus, southern Anatolia, and south-eastern Aegean islands), M. danilewskii (Black Sea region and south-western Anatolia), M. bartoni (Crete), and M. oertzeni (southern Dodecanese Islands). This newly recognized diversity underlines the complex biogeographical history of the Eastern Mediterranean region.

Resolving complex phylogeographic patterns in the Balkan Peninsula using closely related wall-lizard species as a model system.

The Balkan Peninsula constitutes a biodiversity hotspot with high levels of species richness and endemism. The complex geological history of the Balkans in conjunction with the climate evolution are hypothesized as the main drivers generating this biodiversity. We investigated the phylogeography, historical demography, and population structure of closely related wall-lizard species from the Balkan Peninsula and southeastern Europe to better understand diversification processes of species with limited dispersal ability, from Late Miocene to the Holocene. We used several analytical methods integrating genome-wide SNPs (ddRADseq), microsatellites, mitochondrial and nuclear DNA data, as well as species distribution modelling. Phylogenomic analysis resulted in a completely resolved species level phylogeny, population level analyses confirmed the existence of at least two cryptic evolutionary lineages and extensive within species genetic structuring. Divergence time estimations indicated that the Messinian Salinity Crisis played a key role in shaping patterns of species divergence, whereas intraspecific genetic structuring was mainly driven by Pliocene tectonic events and Quaternary climatic oscillations. The present work highlights the effectiveness of utilizing multiple methods and data types coupled with extensive geographic sampling to uncover the evolutionary processes that shaped the species over space and time.

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.

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.

Unraveling the evolutionary history of the Chilostoma Fitzinger, 1833 (Mollusca, Gastropoda, Pulmonata) lineages in Greece.

The land snails of the genus Chilostoma Fitzinger, 1833 that includes, in Greece, the (sub)genera Cattania, Josephinella and Thiessea, are highly diversified and present high levels of endemism. However, their evolutionary history is unknown and their taxonomy is complex and continuously revised. The aim of this study is to investigate the phylogenetic relationships of the lineages of the genus Chilostoma distributed in Greece based on partial DNA sequences of two mitochondrial DNA (16S rRNA and COI) genes. Complete sequences of one nuclear gene (ITS1) representing the major mitochondrial lineages were also analyzed. The phylogenetic trees revealed three distinct major clades that correspond to the three (sub)genera. Several taxonomical incongruencies were made obvious, thus, raising questions about the "true" number of species in each clade, while rendering a taxonomic re-evaluation necessary. From a phylogeographic point of view, it seems that the three major phylogenetic clades were separated in the late Miocene. They started differentiating into distinct species during the Pliocene and Pleistocene through several vicariance and dispersal events.

Non-coding RNA gene families in the genomes of anopheline mosquitoes.

BACKGROUND: Only a small fraction of the mosquito species of the genus Anopheles are able to transmit malaria, one of the biggest killer diseases of poverty, which is mostly prevalent in the tropics. This diversity has genetic, yet unknown, causes. In a further attempt to contribute to the elucidation of these variances, the international "Anopheles Genomes Cluster Consortium" project (a.k.a. "16 Anopheles genomes project") was established, aiming at a comprehensive genomic analysis of several anopheline species, most of which are malaria vectors. In the frame of the international consortium carrying out this project our team studied the genes encoding families of non-coding RNAs (ncRNAs), concentrating on four classes: microRNA (miRNA), ribosomal RNA (rRNA), small nuclear RNA (snRNA), and in particular small nucleolar RNA (snoRNA) and, finally, transfer RNA (tRNA). RESULTS: Our analysis was carried out using, exclusively, computational approaches, and evaluating both the primary NGS reads as well as the respective genome assemblies produced by the consortium and stored in VectorBase; moreover, the results of RNAseq surveys in cases in which these were available and meaningful were also accessed in order to obtain supplementary data, as were "pre-genomic era" sequence data stored in nucleic acid databases. The investigation included the identification and analysis, in most species studied, of ncRNA genes belonging to several families, as well as the analysis of the evolutionary relations of some of those genes in cross-comparisons to other members of the genus Anopheles. CONCLUSIONS: Our study led to the identification of members of these gene families in the majority of twenty different anopheline taxa. A set of tools for the study of the evolution and molecular biology of important disease vectors has, thus, been obtained.

Molecular systematics and historical biogeography of the green lizards (Lacerta) in Greece: insights from mitochondrial and nuclear DNA.

The green lizards of the genus Lacerta (Sauria, Lacertidae) comprise nine recognized species, which in Europe are mainly restricted to the southern peninsulas. Four of them (L. trilineata, L. viridis, L. bilineata and L. agilis) occur in Greece. The uncertainty of morphological diversification renders the taxonomic assignment into species and subspecies problematic. In this study sequence data derived from two mitochondrial (cytochrome b and 16S rRNA) genes and one nuclear (NKTR) gene were used to (a) evaluate the taxonomic status of the genus Lacerta in Greece with emphasis on L. trilineata group and (b) investigate the evolutionary history of the genus through the application of phylogenetic and phylogeographic analyses, using Gallotia and Timon as outgroups. The phylogenetic analyses revealed the existence of four major clades. The first clade corresponds to L. trilineata group, the second to L. media, the third to L. agilis and the fourth to a complex of L. viridis and L. bilineata. However, the produced phylogenetic relationships are not congruent with the current taxonomy, especially in the first clade in which L. trilineata appeared to be paraphyletic in regard to L. pamphylica. Six distinct lineages were inferred within L. trilineata, despite the current recognition of nine morphological subspecies, the genetic differentiation of which exceeds that of other Lacerta species, imposing a thorough taxonomic revision of the species. Our results suggested a rapid diversification of L. trilineata group during the late Miocene. We believe that the present distribution of the genus in Greece is the result of several dispersal and vicariant events that took place during the late Miocene and early Pliocene.

Identification of the 18 World War II executed citizens of Adele, Rethymnon, Crete using an ancient DNA approach and low coverage genomes.

In the Battle of Crete during the World War II occupation of Greece, the German forces faced substantial civilian resistance. To retribute the numerous German losses, a series of mass executions took place in numerous places in Crete; a common practice reported from Greece and elsewhere. In Adele, a village in the regional unit of Rethymnon, 18 male civilians were executed and buried in a burial pit at the Sarakina site. In this study, the first one conducted for a conflict that occurred in Greece, we identified for humanitarian purposes the 18 skulls of the Sarakina victims, following a request from the local community of Adele. The molecular identification of historical human remains via ancient DNA approaches and low coverage whole genome sequencing has only recently been introduced. Here, we performed genome skimming on the living relatives of the victims, as well as high throughput historical DNA analysis on the skulls to infer the kinship degrees among the victims via genetic relatedness analyses. We also conducted targeted anthropological analysis to successfully complete the identification of all Sarakina victims. We demonstrate that our methodological approach constitutes a potentially highly informative forensic tool to identify war victims. It can hence be applied to analogous studies on degraded DNA, thus, paving the path for systematic war victim identification in Greece and beyond.

Molecular phylogeny and historical biogeography of the Anatolian lizard Apathya (Squamata, Lacertidae).

Apathya is a lacertid genus occurring mainly in south-east Turkey and its adjacent regions (part of Iran and Iraq). So far two morphological species have been attributed to the genus; A. cappadocica (with five subspecies, A. c.cappadocica, A. c.muhtari, A. c.schmidtlerorum, A. c. urmiana and A. c.wolteri) and A.yassujica. The first species occupies most of the genus' distribution range, while A. yassujica is endemic of the Zagros Mountains. Here, we explored Apathya's taxonomy and investigated the evolutionary history of the species by employing phylogenetic and phylogeographic approaches and using both mitochondrial (mtDNA) and nuclear markers. The phylogenetic relationships and the genetic distances retrieved, revealed that Apathya is a highly variable genus, which parallels its high morphological variation. Such levels of morphological and genetic differentiation often exceed those between species of other Lacertini genera that are already treated as full species, suggesting the necessity for a taxonomic revision of Apathya. The phylogeographical scenario emerging from the genetic data suggests that the present distribution of the genus was determined by a combination of dispersal and vicariance events between Anatolia and Southwest Asia dating back to the Miocene and continuing up to the Pleistocene. Key geological events for the understanding of the phylogeography of the genus are the movement of the Arabian plate that led to the configuration of Middle East (orogenesis of the mountain ranges of Turkey and Iran) and the formation of Anatolian Diagonal.

Spatial and temporal heterogeneity in human mobility patterns in Holocene Southwest Asia and the East Mediterranean.

We present a spatiotemporal picture of human genetic diversity in Anatolia, Iran, Levant, South Caucasus, and the Aegean, a broad region that experienced the earliest Neolithic transition and the emergence of complex hierarchical societies. Combining 35 new ancient shotgun genomes with 382 ancient and 23 present-day published genomes, we found that genetic diversity within each region steadily increased through the Holocene. We further observed that the inferred sources of gene flow shifted in time. In the first half of the Holocene, Southwest Asian and the East Mediterranean populations homogenized among themselves. Starting with the Bronze Age, however, regional populations diverged from each other, most likely driven by gene flow from external sources, which we term "the expanding mobility model." Interestingly, this increase in inter-regional divergence can be captured by outgroup-f3-based genetic distances, but not by the commonly used FST statistic, due to the sensitivity of FST, but not outgroup-f3, to within-population diversity. Finally, we report a temporal trend of increasing male bias in admixture events through the Holocene.