Ancestral Gene Pool Research Articles

The genomic landscape of Nepalese Tibeto-Burmans reveals new insights into the recent peopling of Southern Himalayas

While much research attention has focused on demographic processes that enabled human diffusion on the Tibetan plateau, little is known about more recent colonization of Southern Himalayas. In particular, the history of migrations, admixture and/or isolation of populations speaking Tibeto-Burman languages, which is supposed to be quite complex and to have reshaped patterns of genetic variation on both sides of the Himalayan arc, remains only partially elucidated. We thus described the genomic landscape of previously unsurveyed Tibeto-Burman (i.e. Sherpa and Tamang) and Indo-Aryan communities from remote Nepalese valleys. Exploration of their genomic relationships with South/East Asian populations provided evidence for Tibetan admixture with low-altitude East Asians and for Sherpa isolation. We also showed that the other Southern Himalayan Tibeto-Burmans derived East Asian ancestry not from the Tibetan/Sherpa lineage, but from low-altitude ancestors who migrated from China plausibly across Northern India/Myanmar, having experienced extensive admixture that reshuffled the ancestral Tibeto-Burman gene pool. These findings improved the understanding of the impact of gene flow/drift on the evolution of high-altitude Himalayan peoples and shed light on migration events that drove colonization of the southern Himalayan slopes, as well as on the role played by different Tibeto-Burman groups in such a complex demographic scenario.

The genome of Chenopodium quinoa

Chenopodium quinoa (quinoa) is a highly nutritious grain identified as an important crop to improve world food security. Unfortunately, few resources are available to facilitate its genetic improvement. Here we report the assembly of a high-quality, chromosome-scale reference genome sequence for quinoa, which was produced using single-molecule real-time sequencing in combination with optical, chromosome-contact and genetic maps. We also report the sequencing of two diploids from the ancestral gene pools of quinoa, which enables the identification of sub-genomes in quinoa, and reduced-coverage genome sequences for 22 other samples of the allotetraploid goosefoot complex. The genome sequence facilitated the identification of the transcription factor likely to control the production of anti-nutritional triterpenoid saponins found in quinoa seeds, including a mutation that appears to cause alternative splicing and a premature stop codon in sweet quinoa strains. These genomic resources are an important first step towards the genetic improvement of quinoa.

Ancestral Origins and Genetic History of Tibetan Highlanders

The origin of Tibetans remains one of the most contentious puzzles in history, anthropology, and genetics. Analyses of deeply sequenced (30×-60×) genomes of 38 Tibetan highlanders and 39 Han Chinese lowlanders, together with available data on archaic and modern humans, allow us to comprehensively characterize the ancestral makeup of Tibetans and uncover their origins. Non-modern human sequences compose ∼6% of the Tibetan gene pool and form unique haplotypes in some genomic regions, where Denisovan-like, Neanderthal-like, ancient-Siberian-like, and unknown ancestries are entangled and elevated. The shared ancestry of Tibetan-enriched sequences dates back to ∼62,000-38,000 years ago, predating the Last Glacial Maximum (LGM) and representing early colonization of the plateau. Nonetheless, most of the Tibetan gene pool is of modern human origin and diverged from that of Han Chinese ∼15,000 to ∼9,000 years ago, which can be largely attributed to post-LGM arrivals. Analysis of ∼200 contemporary populations showed that Tibetans share ancestry with populations from East Asia (∼82%), Central Asia and Siberia (∼11%), South Asia (∼6%), and western Eurasia and Oceania (∼1%). Our results support that Tibetans arose from a mixture of multiple ancestral gene pools but that their origins are much more complicated and ancient than previously suspected. We provide compelling evidence of the co-existence of Paleolithic and Neolithic ancestries in the Tibetan gene pool, indicating a genetic continuity between pre-historical highland-foragers and present-day Tibetans. In particular, highly differentiated sequences harbored in highlanders' genomes were most likely inherited from pre-LGM settlers of multiple ancestral origins (SUNDer) and maintained in high frequency by natural selection.

Ecosystem size matters: the dimensionality of intralacustrine diversification in Icelandic stickleback is predicted by lake size.

Cases of evolutionary diversification can be characterized along a continuum from weak to strong genetic and phenotypic differentiation. Several factors may facilitate or constrain the differentiation process. Comparative analyses of replicates of the same taxon at different stages of differentiation can be useful to identify these factors. We estimated the number of distinct phenotypic groups in three‐spine stickleback populations from nine lakes in Iceland and in one marine population. Using the inferred number of phenotypic groups in each lake, genetic divergence from the marine population, and physical lake and landscape variables, we tested whether ecosystem size, approximated by lake size and depth, or isolation from the ancestral marine gene pool predicts the occurrence and the extent of phenotypic and genetic diversification within lakes. We find intralacustrine phenotypic diversification to be the rule rather than the exception, occurring in all but the youngest lake population and being manifest in ecologically important phenotypic traits. Neutral genetic data further indicate nonrandom mating in four of nine studied lakes, and restricted gene flow between sympatric phenotypic groups in two. Although neither the phenotypic variation nor the number of intralacustrine phenotypic groups was associated with any of our environmental variables, the number of phenotypic traits that were differentiated was significantly positively related to lake size, and evidence for restricted gene flow between sympatric phenotypic groups was only found in the largest lakes where trait specific phenotypic differentiation was highest.

Rapid genome-wide evolution in Brassica rapa populations following drought revealed by sequencing of ancestral and descendant gene pools.

There is increasing evidence that evolution can occur rapidly in response to selection. Recent advances in sequencing suggest the possibility of documenting genetic changes as they occur in populations, thus uncovering the genetic basis of evolution, particularly if samples are available from both before and after selection. Here, we had a unique opportunity to directly assess genetic changes in natural populations following an evolutionary response to a fluctuation in climate. We analysed genome‐wide differences between ancestors and descendants of natural populations of Brassica rapa plants from two locations that rapidly evolved changes in multiple phenotypic traits, including flowering time, following a multiyear late‐season drought in California. These ancestor‐descendant comparisons revealed evolutionary shifts in allele frequencies in many genes. Some genes showing evolutionary shifts have functions related to drought stress and flowering time, consistent with an adaptive response to selection. Loci differentiated between ancestors and descendants (F ST outliers) were generally different from those showing signatures of selection based on site frequency spectrum analysis (Tajima's D), indicating that the loci that evolved in response to the recent drought and those under historical selection were generally distinct. Very few genes showed similar evolutionary responses between two geographically distinct populations, suggesting independent genetic trajectories of evolution yielding parallel phenotypic changes. The results show that selection can result in rapid genome‐wide evolutionary shifts in allele frequencies in natural populations, and highlight the usefulness of combining resurrection experiments in natural populations with genomics for studying the genetic basis of adaptive evolution.

Population genomics reveals seahorses (Hippocampus erectus) of the western mid-Atlantic coast to be residents rather than vagrants.

Understanding population structure and areas of demographic persistence and transients is critical for effective species management. However, direct observational evidence to address the geographic scale and delineation of ephemeral or persistent populations for many marine fishes is limited. The Lined seahorse (Hippocampus erectus) can be commonly found in three western Atlantic zoogeographic provinces, though inhabitants of the temperate northern Virginia Province are often considered tropical vagrants that only arrive during warm seasons from the southern provinces and perish as temperatures decline. Although genetics can locate regions of historical population persistence and isolation, previous evidence of Virginia Province persistence is only provisional due to limited genetic sampling (i.e., mitochondrial DNA and five nuclear loci). To test alternative hypotheses of historical persistence versus the ephemerality of a northern Virginia Province population we used a RADseq generated dataset consisting of 11,708 single nucleotide polymorphisms (SNP) sampled from individuals collected from the eastern Gulf of Mexico to Long Island, NY. Concordant results from genomic analyses all infer three genetically divergent subpopulations, and strongly support Virginia Province inhabitants as a genetically diverged and a historically persistent ancestral gene pool. These results suggest that individuals that emerge in coastal areas during the warm season can be considered “local” and supports offshore migration during the colder months. This research demonstrates how a large number of genes sampled across a geographical range can capture the diversity of coalescent histories (across loci) while inferring population history. Moreover, these results clearly demonstrate the utility of population genomic data to infer peripheral subpopulation persistence in difficult-to-observe species.

Parallel evolutionary divergence in Arctic charSalvelinus alpinuscomplex from Transbaikalia: variation in differentiation degree and segregation of genetic diversity among sympatric forms

Despite the growing recognition of an important role of ecological speciation in evolution, the mechanisms driving sympatric divergence within a single population are not fully understood. We documented a series of parallel divergence events and tried to assess general regularities of the diversification of the ancestral gene pool. For that we analysed variation at microsatellite loci in populations of a highly polymorphic complex species Salvelinus alpinus from 14 Transbaikalian lakes, 10 of them hosting two or three forms: dwarf, small, and large. Our results suggest sympatric or parapatric origin of forms in all 10 lakes (in one lake, two forms out of three). Gene diversity and allelic richness of form populations are positively correlated with lake size and negatively correlated with modal length of mature fish, so that the dwarf form typically has the highest, and the large form, the lowest indices of genetic variation. The latter effect might be caused by differences in life histories, thus reflecting adaptive divergence. Sympatric forms have differently segregated gene pools (FST= 0.030–0.497, RST= 0.011–0.440) and restricted (m = 0.002–0.042) typically asymmetric long-term gene flow. The level of reproductive isolation among forms assessed using putatively neutral microsatellite loci is correlated with their differentiation in morphology, including trophic-related gill raker number, and with lake depth (i.e., with segregation of diets and with habitat diversity), which corresponds to “isolation-by-adaptation” pattern. Our data suggest that the advance in speciation stage apparently more depends upon ecological opportunities of lake ecosystems than upon their age.

Genetic patterns in <i>Pinus nigra</i> from the central Balkans inferred from plastid and mitochondrial data

J.F. Arnold, European black pine, is a typical component of Mediterranean and sub-Mediterranean coniferous forests with highly fragmentary distribution. Western Mediterranean populations of this species have been studied genetically to date, while eastern populations from the central Balkans, which are larger and more abundant, are still genetically understudied. We analyzed seven populations of representing all infraspecific taxa recognized within the central Balkans (subspecies with varieties and Đorđević; and subspecies (Lamb.) Holmboe with varieties and (Endl.) Georgescu et Ionescu), with three chloroplast microsatellites (cpDNA SSRs) and one mitochondrial (mtDNA) locus. Although our molecular data failed to support circumscription of studied infraspecific taxa, we found that genetic patterns at both genomes are in accordance with those found previously in westward populations of this species, that is – exceptionally high levels of genetic diversity ( = 0.949) and low genetic differentiation ( = 0.024) at the cpDNA level, and moderate levels of genetic diversity ( = 0.357) and genetic differentiation ( = 0.358) at the mtDNA level. Based on genealogical relations of mtDNA types currently present in Balkans’ and Iberian/African populations, we inferred that the ancestral gene pool of already harbored polymorphism at position 328 prior to the divergence to two lineages currently present in westward and eastward parts of the species range distribution. Subsequent occurrence of three mutations, which distinguish these two lineages, suggests their long-term isolation.Pinus nigraP. nigranigranigragocensis pallasianapallasianabanaticaHTGSTHTGSTP. nigra

The congruence between matrilineal genetic (mtDNA) and geographic diversity of Iranians and the territorial populations.

From the ancient era, emergence of Agriculture in the connecting region of Mesopotamia and the Iranian plateau at the foothills of the Zagros Mountains, made Iranian gene pool as an important source of populating the region. It has differentiated the population spread and different language groups. In order to trace the maternal genetic affinity between Iranians and other populations of the area and to establish the place of Iranians in a broad framework of ethnically and linguistically diverse groups of Middle Eastern and South Asian populations, a comparative study of territorial groups was designed and used in the population statistical analysis. Mix of 616 samples was sequenced for complete mtDNA or hyper variable regions in this study. A published dataset of neighboring populations was used as a comparison in the Iranian matrilineal lineage study based on mtDNA haplogroups. Statistical analyses data, demonstrate a close genetic structure of all Iranian populations, thus suggesting their origin from a common maternal ancestral gene pool and show that the diverse maternal genetic structure does not reflect population differentiation in the region in their language. In the aggregate of the eastward spreads of proto-Elamo-Dravidian language from the Southwest region of Iran, the Elam province, a reasonable degree of homogeneity has been observed among Iranians in this study. The approach will facilitate our perception of the more detailed relationship of the ethnic groups living in Iran with the other ancient peoples of the area, testing linguistic hypothesis and population movements.

Mycobacterial Pathogenomics and Evolution.

Most mycobacterial species are harmless saprophytes, often found in aquatic environments. A few species seem to have evolved from this pool of environmental mycobacteria into major human pathogens, such as Mycobacterium tuberculosis, the agent of tuberculosis, Mycobacterium leprae, the leprosy bacillus, and Mycobacterium ulcerans, the agent of Buruli ulcer. While the pathogenicity of M. ulcerans relates to the acquisition of a large plasmid encoding a polyketide-derived toxin, the molecular mechanisms by which M. leprae or M. tuberculosis have evolved to cause disease are complex and involve the interaction between the pathogen and the host. Here we focus on M. tuberculosis and closely related mycobacteria and discuss insights gained from recent genomic and functional studies. Comparison of M. tuberculosis genome data with sequences from nontuberculous mycobacteria, such as Mycobacterium marinum or Mycobacterium kansasii, provides a perception of the more distant evolution of M. tuberculosis, while the recently accomplished genome sequences of multiple tubercle bacilli with smooth colony morphology, named Mycobacterium canettii, have allowed the ancestral gene pool of tubercle bacilli to be estimated. The resulting findings are instrumental for our understanding of the pathogenomic evolution of tuberculosis-causing mycobacteria. Comparison of virulent and attenuated members of the M. tuberculosis complex has further contributed to identification of a specific secretion pathway, named ESX or Type VII secretion. The molecular machines involved are key elements for mycobacterial pathogenicity, strongly influencing the ability of M. tuberculosis to cope with the immune defense mounted by the host.

Unravelling the evolutionary history of Eucalyptus cordata (Myrtaceae) using molecular markers

We studied the evolutionary processes shaping the genetic diversity in the naturally fragmented Eucalyptus cordata, a rare homoblastic tree endemic to the island of Tasmania. A genome-wide scan showed that E. cordata and the endangered heteroblastic E. morrisbyi were closely related, suggesting a neotenous origin of E. cordata from an endemic heteroblastic ancestor. Bayesian cluster analysis based on nuclear microsatellites assayed in 567 E. cordata and E. morrisbyi individuals revealed five genetic clusters. Two clusters comprised populations that correspond to putative ancestral gene pools linking E. cordata and E. morrisbyi. Another cluster included populations that transgressed the drowned Derwent River valley, suggestive of a wider glacial distribution. However, the majority of individuals occurred in the two genetic clusters distributed in the south-west and north-east of the range of E. cordata. The elevated genetic diversity in populations comprising these clusters suggests that they represent two recently fragmented cores of the distribution. Genetic evidence suggests that the newly described, localised E. cordata subspecies quadrangulosa has been recently selected from within the morphologically diverse, south-western cluster. We argue that multiple phases of isolation and drift have led to the contemporary pattern of molecular variation and the scattering of relictual and more recently derived populations across the species distribution.

The complex history of the olive tree: from Late Quaternary diversification of Mediterranean lineages to primary domestication in the northern Levant

The location and timing of domestication of the olive tree, a key crop in Early Mediterranean societies, remain hotly debated. Here, we unravel the history of wild olives (oleasters), and then infer the primary origins of the domesticated olive. Phylogeography and Bayesian molecular dating analyses based on plastid genome profiling of 1263 oleasters and 534 cultivated genotypes reveal three main lineages of pre-Quaternary origin. Regional hotspots of plastid diversity, species distribution modelling and macrofossils support the existence of three long-term refugia; namely the Near East (including Cyprus), the Aegean area and the Strait of Gibraltar. These ancestral wild gene pools have provided the essential foundations for cultivated olive breeding. Comparison of the geographical pattern of plastid diversity between wild and cultivated olives indicates the cradle of first domestication in the northern Levant followed by dispersals across the Mediterranean basin in parallel with the expansion of civilizations and human exchanges in this part of the world.

Inconsistent genetic structure among members of a multitrophic system: did bruchid parasitoids (Horismenusspp.) escape the effects of bean domestication?

Anthropogenic range expansion and cultural practices have modified the distribution, abundance and genetic diversity of domesticated organisms, thereby altering multitrophic assemblages through space and time. The putative Mesoamerican domestication centre of the common bean, Phaseolus vulgaris L., in Mexico allows investigating the effects of plant domestication on the genetic structure of members of a multitrophic system. The aim of this study was to compare the evolutionary history of Horismenus parasitoids (Hymenoptera: Eulophidae) to those of their bruchid beetle hosts (Coleoptera: Bruchidae) and their domesticated host plant (P. vulgaris), in the context of traditional agriculture in Mexico. We analyzed the population genetic structure of four Horismenus species in Mexico using mitochondrial COI haplotype data. The two most abundant parasitoid species were Horismenus depressus and Horismenus missouriensis. Horismenus missouriensis were infected by Wolbachia endosymbionts and had little to no population differentiation (F(ST) = 0.06). We suspect the mitochondrial history of H. missouriensis to be blurred by Wolbachia, because differentiation among infected vs. non-infected individuals exists (F(ST) = 0.11). Populations of H. depressus were found to be highly differentiated (F(ST) = 0.34), but the genetic structuring could not be explained by tested spatial components. We then compared the genetic structure observed in this parasitoid species to previously published studies on bruchid beetles and their host plants. Despite extensive human-mediated migration and likely population homogenization of its two Acanthoscelides bruchid beetle hosts, H. depressus populations are structured like its host plant, by a recent dispersal from a diverse ancestral gene pool. Distinct evolutionary dynamics may explain inconsistent patterns among trophic levels. Parasitoids likely migrate from wild bean populations and are poorly adapted to bean storage conditions similar to their bruchid beetle hosts. Integrating several trophic levels to the study of evolutionary history has proven to be fruitful in detecting different ecological responses to human-mediated disturbances and host parasite interactions.

Species limits in Vachellia (Acacia) karroo (Mimosoideae: Leguminoseae): Evidence from automated ISSR DNA “fingerprinting”

Vachellia karroo (Hayne) Banfi & Galasso, previously known as Acacia karroo Hayne, is a very common woody species in South Africa, and displays a large amount of morphological and architectural variation. Previous morphological studies have provided evidence to support the recognition of a number of segregate species, but this has not been supported by early protein electrophoresis studies. Genetic variability of the species throughout South Africa was examined by means of the Inter-Simple Sequence Repeat (ISSR) DNA “fingerprinting” to determine whether there is any genetic structure that correlates to the morphological diversity. Based upon 33 samples from across the South African distribution of this species, we found no evidence for any genetic structuring, suggesting the species is one panmictic entity. The existence of morphologically discrete varieties in light of this genetic finding is difficult to explain, and perhaps the species ought to be considered as an ochlospecies. A scenario whereby the various morphotypes evolved in refugia during previous climates followed by subsequent expansion and introgression during the current interglacial may explain this conundrum. Alternatively, the evolution of the morphotypes has been recent and rapid, and the genetic variation observed here represents the ancestral gene pool that has not yet undergone lineage sorting as a consequence of isolation.

Evolution of HLA-DRB Genes

The HLA region shows diversity concerning the number and content of DRB genes present per haplotype. Similar observations are made for the equivalent regions in other primate species. To elucidate the evolutionary history of the various HLA-DRB genes, a large panel of intron sequences obtained from humans, chimpanzees, rhesus macaques, and common marmosets has been subjected to phylogenetic analyses. Special attention was paid to the presence and absence of particular transposable elements and/or to their segments. The sharing of different parts of the same long interspersed nuclear element-2 (LINE2, L2) and various Alu insertions by the species studied demonstrates that one precursor gene must have been duplicated several times before the Old World monkey (OWM) and hominid (HOM) divergence. At least four ancestral DRB gene families appear to have been present before the radiation of OWM and HOM, and one of these even predates the speciation of Old and New World primates. Two of these families represent the pseudogenes DRB6/DRB2 and DRB7, which have been locked in the genomes of various primate species over long evolutionary time spans. Furthermore, all phylogenies of different intron segments show consistently that, apart from the pseudogenes, only DRB5 genes are shared by OWM and HOM, and they demonstrate the common history of certain DRB genes/lineages of humans and chimpanzees. In contrast, the evolutionary history of some other DRB loci is difficult to decipher, thus illustrating the complex history of the evolution of DRB genes due to a combination of mutations and recombination-like events. The selected approach allowed us to shed light on the ancestral DRB gene pool in primates and on the evolutionary relationship of the various HLA-DRB genes.

Evolution and Association Analysis of Ghd7 in Rice

Plant height, heading date, and yield are the main targets for rice genetic improvement. Ghd7 is a pleiotropic gene that controls the aforementioned traits simultaneously. In this study, a rice germplasm collection of 104 accessions (Oryza sativa) and 3 wild rice varieties (O.rufipogon) was used to analyze the evolution and association of Ghd7 with plant height, heading date, and yield. Among the 104 accessions, 76 single nucleotide polymorphisms (SNPs) and six insertions and deletions were found within a 3932-bp DNA fragment of Ghd7. A higher pairwise π and θ in the promoter indicated a highly diversified promoter of Ghd7. Sixteen haplotypes and 8 types of Ghd7 protein were detected. SNP changes between haplotypes indicated that Ghd7 evolved from two distinct ancestral gene pools, and independent domestication processes were detected in indica and japonica varietals respectively. In addition to the previously reported premature stop mutation in the first exon of Ghd7, which caused phenotypic changes of multiple traits, we found another functional C/T mutation (SNP S_555) by structure-based association analysis. SNP S_555 is located in the promoter and was related to plant height probably by altering gene expression. Moreover, another seven SNP mutations in complete linkage were found to be associated with the number of spikelets per panicle, regardless of the photoperiod. These associations provide the potential for flexibility of Ghd7 application in rice breeding programs.

Complete mitochondrial DNA analysis of eastern Eurasian haplogroups rarely found in populations of northern Asia and eastern Europe.

With the aim of uncovering all of the most basal variation in the northern Asian mitochondrial DNA (mtDNA) haplogroups, we have analyzed mtDNA control region and coding region sequence variation in 98 Altaian Kazakhs from southern Siberia and 149 Barghuts from Inner Mongolia, China. Both populations exhibit the prevalence of eastern Eurasian lineages accounting for 91.9% in Barghuts and 60.2% in Altaian Kazakhs. The strong affinity of Altaian Kazakhs and populations of northern and central Asia has been revealed, reflecting both influences of central Asian inhabitants and essential genetic interaction with the Altai region indigenous populations. Statistical analyses data demonstrate a close positioning of all Mongolic-speaking populations (Mongolians, Buryats, Khamnigans, Kalmyks as well as Barghuts studied here) and Turkic-speaking Sojots, thus suggesting their origin from a common maternal ancestral gene pool. In order to achieve a thorough coverage of DNA lineages revealed in the northern Asian matrilineal gene pool, we have completely sequenced the mtDNA of 55 samples representing haplogroups R11b, B4, B5, F2, M9, M10, M11, M13, N9a and R9c1, which were pinpointed from a massive collection (over 5000 individuals) of northern and eastern Asian, as well as European control region mtDNA sequences. Applying the newly updated mtDNA tree to the previously reported northern Asian and eastern Asian mtDNA data sets has resolved the status of the poorly classified mtDNA types and allowed us to obtain the coalescence age estimates of the nodes of interest using different calibrated rates. Our findings confirm our previous conclusion that northern Asian maternal gene pool consists of predominantly post-LGM components of eastern Asian ancestry, though some genetic lineages may have a pre-LGM/LGM origin.

Fossil record of an archaeal HK97-like provirus

One of the outstanding questions in biology today is the origin of viruses. We have discovered a protein in the hyperthermophile Sulfolobus solfataricus while following proteome regulation during viral infection that led to the discovery of a fossil provirus. Characterization of the wild type and recombinant protein revealed that it assembled into virus-like particles with a diameter of ~ 32 nm. Sequence and structural analyses showed that the likely proviral capsid protein, Sso2749, is homologous to a protein from Pyrococcus furiosus that forms virus-like particles using the HK-97 major capsid protein fold. The SsP2-provirus appears mosaic and contains proteins with similarity to, among others, eukaryotic herpesviruses and tailed dsDNA bacteriophage families, reinforcing the hypothesis of a common ancestral gene pool across all three domains of life. This is the first description of the HK-97 fold in a crenarchaeal virus and the first direct genomic connection of linocin-like protein cages to a virus.

A Comparison of Y-Chromosome Variation in Sardinia and Anatolia Is More Consistent with Cultural Rather than Demic Diffusion of Agriculture

Two alternative models have been proposed to explain the spread of agriculture in Europe during the Neolithic period. The demic diffusion model postulates the spreading of farmers from the Middle East along a Southeast to Northeast axis. Conversely, the cultural diffusion model assumes transmission of agricultural techniques without substantial movements of people. Support for the demic model derives largely from the observation of frequency gradients among some genetic variants, in particular haplogroups defined by single nucleotide polymorphisms (SNPs) in the Y-chromosome. A recent network analysis of the R-M269 Y chromosome lineage has purportedly corroborated Neolithic expansion from Anatolia, the site of diffusion of agriculture. However, the data are still controversial and the analyses so far performed are prone to a number of biases. In the present study we show that the addition of a single marker, DYSA7.2, dramatically changes the shape of the R-M269 network into a topology showing a clear Western-Eastern dichotomy not consistent with a radial diffusion of people from the Middle East. We have also assessed other Y-chromosome haplogroups proposed to be markers of the Neolithic diffusion of farmers and compared their intra-lineage variation—defined by short tandem repeats (STRs)—in Anatolia and in Sardinia, the only Western population where these lineages are present at appreciable frequencies and where there is substantial archaeological and genetic evidence of pre-Neolithic human occupation. The data indicate that Sardinia does not contain a subset of the variability present in Anatolia and that the shared variability between these populations is best explained by an earlier, pre-Neolithic dispersal of haplogroups from a common ancestral gene pool. Overall, these results are consistent with the cultural diffusion and do not support the demic model of agriculture diffusion.

Genetic diversity, population structure and taxonomy of Calopteryx splendens (Odonata: Calopterygidae): An AFLP analysis

Calopteryx splendens is a widely distributed palaearctic damselfly with a remarkably uniform morphology. Variation in the size and shape of the pigmented spot on the wing is the main diagnostic character used to discriminate subspecies across its huge geographic range. Here, AFLP analysis was used to assess the genetic structure and diversity of nine populations representing 3 putative subspecies and evaluate the pigment spot as a taxonomic marker. Genetic diversity was high, with the number of polymorphic loci per population ranging from 141 to 280 out of a total of 333 variable sites (42.3-84.1%) and Nei's gene diversity from 0.160 to 0.283 (overall 0.299). Overall population genetic differentiation (FST = 0.2766) suggests limited gene flow and adaptation to local environments. Restricted gene flow and genetic differentiation among populations are supported by significant FST estimates. High levels of gene flow (Nm >1) were only recorded among three Asian populations (Russia - Kazakhstan - Turkey). The patterns of genotypic diversity suggest that a given wing spot size and shape may arise from the hybridization of a limited number, possibly not more than four, ancestral gene pools in different ways and at different times. Clearly, the sample analyzed was not sufficient to capture all of the complex history of C. splendens, but sufficient to indicate the taxa ancilla, waterstoni, and orientalis possibly represent three of the four ancestral gene pools, and originated in western Asia. The origin of the fourth, xanthostoma, is the western Mediterranean.