Landscape genetics and the genetic legacy of Upper Paleolithic and Mesolithic hunter-gatherers in the modern Caucasus.

This study clarifies the role of refugia and landscape permeability in the formation of the current genetic structure of peoples of the Caucasus. We report novel genome-wide data for modern individuals from the Caucasus, and analyze them together with available Paleolithic and Mesolithic individuals from Eurasia and Africa in order (1) to link the current and ancient genetic structures via landscape permeability, and (2) thus to identify movement paths between the ancient refugial populations and the Caucasus. The ancient genetic ancestry is best explained by landscape permeability implying that human movement is impeded by terrain ruggedness, swamps, glaciers and desert. Major refugial source populations for the modern Caucasus are those of the Caucasus, Anatolia, the Balkans and Siberia. In Rugged areas new genetic signatures take a long time to form, but once they do so, they remain for a long time. These areas act as time capsules harboring genetic signatures of ancient source populations and making it possible to help reconstruct human history based on patterns of variation today.

Genotypic similarities among the parthenogenetic Darevskia rock lizards with different hybrid origins.

BACKGROUND: The majority of parthenogenetic vertebrates derive from hybridization between sexually reproducing species, but the exact number of hybridization events ancestral to currently extant clonal lineages is difficult to determine. Usually, we do not know whether the parental species are able to contribute their genes to the parthenogenetic vertebrate lineages after the initial hybridization. In this paper, we address the hypothesis, whether some genotypes of seven phenotypically distinct parthenogenetic rock lizards (genus Darevskia) could have resulted from back-crosses of parthenogens with their presumed parental species. We also tried to identify, as precise as possible, the ancestral populations of all seven parthenogens. RESULTS: We analysed partial mtDNA sequences and microsatellite genotypes of all seven parthenogens and their presumed ansectral species, sampled across the entire geographic range of parthenogenesis in this group. Our results confirm the previous designation of the parental species, but further specify the maternal populations that are likely ancestral to different parthenogenetic lineages. Contrary to the expectation of independent hybrid origins of the unisexual taxa, we found that genotypes at multiple loci were shared frequently between different parthenogenetic species. The highest proportions of shared genotypes were detected between (i) D. sapphirina and D. bendimahiensis and (ii) D. dahli and D. armeniaca, and less often between other parthenogens. In case (ii), genotypes at the remaining loci were notably distinct. CONCLUSIONS: We suggest that both observations (i-ii) can be explained by two parthenogenetic forms tracing their origin to a single initial hybridization event. In case (ii), however, occasional gene exchange between the unisexual and the parental bisexual species could have taken place after the onset of parthenogenetic reproduction. Indeed, backcrossed polyploid hybrids are relatively frequent in Darevskia, although no direct evidence of recent gene flow has been previously documented. Our results further suggest that parthenogens are losing heterozygosity as a result of allelic conversion, hence their fitness is expected to decline over time as genetic diversity declines. Backcrosses with the parental species could be a rescue mechanism which might prevent this decline, and therefore increase the persistance of unisexual forms.

Phylogeny of caucasian rock lizards (Darevskia) and other true lizards based on mitogenome analysis: Optimisation of the algorithms and gene selection.

We generated a phylogeny for Caucasian rock lizards (Darevskia), and included six other families of true lizards (Lacertini), based on complete mitochondrial genome analysis. Next-generation sequencing (NGS) of genomic DNA was used to obtain 16 new mitogenomes of Darevskia. These, along with 35 sequences downloaded from GenBank: genera Darevskia, Zootoca, Podarcis, Phoenicolacerta, Takydromus, Lacerta, and Eremias-were used in the analysis. All four analytical methods (Bayesian Inference, BI; Maximum Likelihood, ML; Maximum Parsimony, MP; and Neighbor-Joining, NJ) showed almost congruent intra-generic topologies for Darevskia and other lizard genera. However, ML and NJ methods on one side, and BI and MP methods on the other harvested conflicting phylogenies. The ML/NJ topology supports earlier published separation of Darevskia into three mitochondrial clades (Murphy, Fu, Macculloch, Darevsky, and Kupinova, 2000), but BI and MP topologies support that the basal branching occurred between D. parvula from the western Lesser Caucasus and the rest of Darevskia. All topologies altered the phylogenetic position of some individual species, including D. daghestanica, D. derjugini, and D. chlorogaster. Reanalysis after excluding four saturated genes from the data set, and excluding genus Eremias gives fully convergent topologies. The most basal branching for true lizards was between Far Eastern Takydromus and the Western Eurasian genera (BI). Comparing phylogenetic performance of individual genes relative to whole mitogenome data, concatenated 16S RNA (the least saturated gene in our analyses) and Cytochrome b genes generate a robust phylogeny that is fully congruent with that based on the complete mitogenome.

Partial genomic survival of cave bears in living brown bears.

Although many large mammal species went extinct at the end of the Pleistocene epoch, their DNA may persist due to past episodes of interspecies admixture. However, direct empirical evidence of the persistence of ancient alleles remains scarce. Here, we present multifold coverage genomic data from four Late Pleistocene cave bears (Ursus spelaeus complex) and show that cave bears hybridized with brown bears (Ursus arctos) during the Pleistocene. We develop an approach to assess both the directionality and relative timing of gene flow. We find that segments of cave bear DNA still persist in the genomes of living brown bears, with cave bears contributing 0.9 to 2.4% of the genomes of all brown bears investigated. Our results show that even though extinction is typically considered as absolute, following admixture, fragments of the gene pool of extinct species can survive for tens of thousands of years in the genomes of extant recipient species.

Landscape Complexity in the Caucasus Impedes Genetic Assimilation of Human Populations More Effectively than Language or Ethnicity.

The analyses of 15 autosomal and 23 Y-chromosome DNA single-tandem-repeat loci in five rural populations from the Caucasus (four ethnically Georgian and one ethnically Armenian) indicated that two Georgian populations, one from the west and the other from the east of the Greater Caucasus Mountains, were both patrilineally and autosomally most differentiated from each other, and the other populations of Georgians and Armenians held an intermediate position between those two. This pattern may be due to human dispersal from two distinct glacial refugia in the last glacial period and the early Holocene, followed by less gene flow among the populations from the Greater Caucasus than among those from the rest of the Caucasus, where the populations have undergone substantial admixture in historical time. This hypothesis is supported by a strong correlation between genetic differentiation among the populations and landscape permeability to human migrations as determined by terrain ruggedness, forest cover, and snow cover. Although geographic patterns of autosomal and Y-chromosome DNA are not fully concordant, both are influenced by landscape permeability and show a similar east-west gradient. Our results suggest that this permeability was a stronger factor limiting gene flow among human populations in the Caucasus than were ethnic or linguistic boundaries.

Human paternal lineages, languages, and environment in the Caucasus.

Publications that describe the composition of the human Y-DNA haplogroup in diffferent ethnic or linguistic groups and geographic regions provide no explicit explanation of the distribution of human paternal lineages in relation to specific ecological conditions. Our research attempts to address this topic for the Caucasus, a geographic region that encompasses a relatively small area but harbors high linguistic, ethnic, and Y-DNA haplogroup diversity. We genotyped 224 men that identified themselves as ethnic Georgian for 23 Y-chromosome short tandem-repeat markers and assigned them to their geographic places of origin. The genotyped data were supplemented with published data on haplogroup composition and location of other ethnic groups of the Caucasus. We used multivariate statistical methods to see if linguistics, climate, and landscape accounted for geographical diffferences in frequencies of the Y-DNA haplogroups G2, R1a, R1b, J1, and J2. The analysis showed significant associations of (1) G2 with wellforested mountains, (2) J2 with warm areas or poorly forested mountains, and (3) J1 with poorly forested mountains. R1b showed no association with environment. Haplogroups J1 and R1a were significantly associated with Daghestanian and Kipchak speakers, respectively, but the other haplogroups showed no such simple associations with languages. Climate and landscape in the context of competition over productive areas among diffferent paternal lineages, arriving in the Caucasus in diffferent times, have played an important role in shaping the present-day spatial distribution of patrilineages in the Caucasus. This spatial pattern had formed before linguistic subdivisions were finally shaped, probably in the Neolithic to Bronze Age. Later historical turmoil had little influence on the patrilineage composition and spatial distribution. Based on our results, the scenario of postglacial expansions of humans and their languages to the Caucasus from the Middle East, western Eurasia, and the East European Plain is plausible.

Geographic and genetic boundaries of brown bear (Ursus arctos) population in the Caucasus.

The taxonomic status of brown bears in the Caucasus remains unclear. Several morphs or subspecies have been identified from the morphological (craniological) data, but the status of each of these subspecies has never been verified by molecular genetic methods. We analysed mitochondrial DNA sequences (control region) to reveal phylogenetic relationships and infer divergence time between brown bear subpopulations in the Caucasus. We estimated migration and gene flow from both mitochondrial DNA and microsatellite allele frequencies, and identified possible barriers to gene flow among the subpopulations. Our suggestion is that all Caucasian bears belong to the nominal subspecies of Ursus arctos. Our results revealed two genetically and geographically distinct maternal haplogroups: one from the Lesser Caucasus and the other one from the Greater Caucasus. The genetic divergence between these haplogroups dates as far back as the beginning of human colonization of the Caucasus. Our analysis of the least-cost distances between the subpopulations suggests humans as a major barrier to gene flow. The low genetic differentiation inferred from microsatellite allele frequencies indicates that gene flow between the two populations in the Caucasus is maintained through the movements of male brown bears. The Likhi Ridge that connects the Greater and Lesser Caucasus mountains is the most likely corridor for this migration.