The emerging complexity of Open Science: assessing Intelligent Data Openness in Genomic Anthropology and Human Genomics

In recent decades, the scientific community has become aware of the importance of science being effectively open in order to speed up scientific and technological progress. In this context, the achievement of a robust, effective and responsible form of data sharing is now widely acknowledged as a fundamental part of the research process. The production and resolution of human genomic data has steadily increased in recent years, mainly due to technological advances and decreasing costs of DNA genotyping and sequencing. There is, however, a downside to this process due to the huge increase in the complexity of the data and related metadata. This means it is advisable to go beyond traditional forms of sharing analysis, which have focused on data availability only. Here we present a pilot study that aims to complement a survey on the availability of data related to peer-reviewed publications with an analysis of their findability, accessibility, useability and assessability (according to the "intelligent data openness" scheme). Sharing rates in genomic anthropology (73.0%) were found to be higher than human genomics (32.4%), but lower than closely related research fields (from 96.8% to 79.2% for paleogenetics and evolutionary genetics, respectively). We discuss the privacy and methodological issues that could be linked to this finding. Comparisons of sharing rates across a wide range of disciplines has suggested that the idea of human genomics as a forerunner for the open data movement should be questioned. Finally, both in genomic anthropology and human genomics, findability and useability were found to be compliant with the expectations of an intelligent data openness, whereas only a minor part of studies met the need to make the data completely assessable.

Population structure of modern-day Italians reveals patterns of ancient and archaic ancestries in Southern Europe

European populations display low genetic differentiation as the result of long-term blending of their ancient founding ancestries. However, it is unclear how the combination of ancient ancestries related to early foragers, Neolithic farmers, and Bronze Age nomadic pastoralists can explain the distribution of genetic variation across Europe. Populations in natural crossroads like the Italian peninsula are expected to recapitulate the continental diversity, but have been systematically understudied. Here, we characterize the ancestry profiles of Italian populations using a genome-wide dataset representative of modern and ancient samples from across Italy, Europe, and the rest of the world. Italian genomes capture several ancient signatures, including a non-steppe contribution derived ultimately from the Caucasus. Differences in ancestry composition, as the result of migration and admixture, have generated in Italy the largest degree of population structure detected so far in the continent, as well as shaping the amount of Neanderthal DNA in modern-day populations.

Genetic susceptibility in pharmacodynamic and pharmacokinetic pathways underlying drug-induced arrhythmia and sudden unexplained deaths

Drug-induced arrhythmia is an adverse drug reaction that can be potentially fatal since it is mostly related to drug-induced QT prolongation, a known risk factor for Torsade de Pointes and sudden cardiac death (SCD). Several risk factors have been described in association to these drug-induced events, such as preexistent cardiac disease and genetic variation. Our objective was to study the genetic susceptibility in pharmacodynamic and pharmacokinetic pathways underlying suspected drug-induced arrhythmias and sudden unexplained deaths in 32 patients. The genetic component in the pharmacodynamic pathway was studied by analysing 96 genes associated with higher risk of SCD through massive parallel sequencing. Pharmacokinetic-mediated genetic susceptibility was investigated by studying the genes encoding cytochrome P450 enzymes using medium-throughput genotyping. Pharmacodynamic analysis showed three probably pathogenic variants and 45 variants of uncertain significance in 28 patients, several of them previously described in relation to mild or late onset cardiomyopathies. These results suggest that genetic variants in cardiomyopathy genes, in addition to those related with channelopathies, could be relevant to drug-induced cardiotoxicity and contribute to the arrhythmogenic phenotype. Pharmacokinetic analysis showed three patients that could have an altered metabolism of the drugs they received involving CYP2C19 and/or CYP2D6, probably contributing to the arrhythmogenic phenotype. The study of genetic variants in both pharmacodynamic and pharmacokinetic pathways may be a useful strategy to understand the multifactorial mechanism of drug-induced events in both clinical practice and forensic field. However, it is necessary to comprehensively study and evaluate the contribution of the genetic susceptibility to drug-induced cardiotoxicity.

Iron Age Italic population genetics: the Piceni from Novilara (8th-7th century BC)

BACKGROUND

Archaeological data provide evidence that Italy, during the Iron Age, witnessed the appearance of the first communities with well defined cultural identities. To date, only a few studies report genetic data about these populations and, in particular, the Piceni have never been analysed.

AIMS

To provide new data about mitochondrial DNA (mtDNA) variability of an Iron Age Italic population, to understand the contribution of the Piceni in shaping the modern Italian gene pool and to ascertain the kinship between some individuals buried in the same grave within the Novilara necropolis.

SUBJECTS AND METHODS

In a first set of 10 individuals from Novilara, we performed deep sequencing of the HVS-I region of the mtDNA, combined with the genotyping of 22 SNPs in the coding region and the analysis of several autosomal markers.

RESULTS

The results show a low nucleotide diversity for the inhabitants of Novilara and highlight a genetic affinity of this ancient population with the current inhabitants of central Italy. No family relationship was observed between the individuals analysed here.

CONCLUSIONS

This study provides a preliminary characterisation of the mtDNA variability of the Piceni of Novilara, as well as a kinship assessment of two peculiar burials.

Charting the Y-chromosome ancestry of present-day Argentinean Mennonites

Old Order Mennonite communities initially arose in Northern Europe (centered in the Netherlands) and derived from the Anabaptist movement of the 16th century. Mennonites migrated to the New World in the early 18th century, first to North America, and more recently to Mesoamerica and South America. We analyzed Y-chromosome short tandem repeats (STRs) and single nucleotide polymorphisms in males from a community of Mennonites, 'La Nueva Esperanza', which arrived to Argentina in 1985 from colonies in Bolivia and Mexico. Molecular diversity indices coupled with demographic simulations show that Mennonites have a reduced variability when compared with local Argentinean populations and 69 European population samples. Mennonite Y-STR haplotypes were mainly observed in Central Europe. In agreement, multidimensional scaling analyses based on RST genetic distances indicate that Mennonite Y-chromosomes are closely related to Central/Northern Europeans (the Netherlands, Switzerland and Denmark). In addition, statistical inferences made on the most likely geographic origin of Y-chromosome haplotypes point more specifically to the Netherlands as the populations that best represent the majority of the Mennonite Y-chromosomes. Overall, Y-chromosome variation of Mennonites shows the signatures of moderate reduction of variability when compared with source populations, which is in good agreement with their lifestyle in small endogamous demes. These genetic singularities could also help to understand disease conditions that are more prevalent among Mennonites.

Whole mitochondrial DNA sequencing in Alpine populations and the genetic history of the Neolithic Tyrolean Iceman

The Tyrolean Iceman is an extraordinarily well-preserved natural mummy that lived south of the Alpine ridge ~5,200 years before present (ybp), during the Copper Age. Despite studies that have investigated his genetic profile, the relation of the Iceman´s maternal lineage with present-day mitochondrial variation remains elusive. Studies of the Iceman have shown that his mitochondrial DNA (mtDNA) belongs to a novel lineage of haplogroup K1 (K1f) not found in extant populations. We analyzed the complete mtDNA sequences of 42 haplogroup K bearing individuals from populations of the Eastern Italian Alps - putatively in genetic continuity with the Tyrolean Iceman-and compared his mitogenome with a large dataset of worldwide K1 sequences. Our results allow a re-definition of the K1 phylogeny, and indicate that the K1f haplogroup is absent or rare in present-day populations. We suggest that mtDNA Iceman´s lineage could have disappeared during demographic events starting in Europe from ~5,000 ybp. Based on the comparison of our results with published data, we propose a scenario that could explain the apparent contrast between the phylogeographic features of maternal and paternal lineages of the Tyrolean Iceman within the context of the demographic dynamics happening in Europe from 8,000 ybp.

Analysis of Y-chromosome STRs in Chile confirms an extensive introgression of European male lineages in urban populations

We analyzed the Y chromosome haplotypes (Yfiler) of 978 non-related Chilean males grouped in five sampling regions (Iquique, Santiago de Chile, Concepción, Temuco and Punta Arenas) covering main geo-political regions. Overall, 803 different haplotypes and 688 singletons were observed. Molecular diversity was moderately lower than in other neighboring countries (e.g. Argentina); and AMOVA analysis on Y-STR haplotypes showed that among variation within Chile accounted for only 0.25% of the total variation. Punta Arenas, in the southern cone, showed the lowest haplotype diversity, and discrimination capacity, and also the highest matching probability of the five Chilean samples, probably reflecting its more marked geographic isolation compared to the other regions. Multidimensional scaling (MDS) analysis based on RST genetic distances suggested a close proximity of Chilean Y-chromosome profiles to European ones. Consistently, haplogroups inferred from Y-STR profiles revealed that the Native American component constituted only 8% of all the haplotypes, and this component ranged from 5% in the Centre of the country to 9-10% in the South and 13% in the North, which is in good agreement with the distribution of Native American communities in these regions. AMOVA computed on inferred haplogroups confirmed the very low among variation observed in Chilean populations. The present project provides the first Chilean dataset to the international Y-chromosome STR Haplotype Reference Database (YHRD) and it is also the first reference database for Y-chromosome forensic casework of the country.

The Role of Recent Admixture in Forming the Contemporary West Eurasian Genomic Landscape

Over the past few years, studies of DNA isolated from human fossils and archaeological remains have generated considerable novel insight into the history of our species. Several landmark papers have described the genomes of ancient humans across West Eurasia, demonstrating the presence of large-scale, dynamic population movements over the last 10,000 years, such that ancestry across present-day populations is likely to be a mixture of several ancient groups [1, 2, 3, 4, 5, 6, 7]. While these efforts are bringing the details of West Eurasian prehistory into increasing focus, studies aimed at understanding the processes behind the generation of the current West Eurasian genetic landscape have been limited by the number of populations sampled or have been either too regional or global in their outlook [8, 9, 10, 11]. Here, using recently described haplotype-based techniques [11], we present the results of a systematic survey of recent admixture history across Western Eurasia and show that admixture is a universal property across almost all groups. Admixture in all regions except North Western Europe involved the influx of genetic material from outside of West Eurasia, which we date to specific time periods. Within Northern, Western, and Central Europe, admixture tended to occur between local groups during the period 300 to 1200 CE. Comparisons of the genetic profiles of West Eurasians before and after admixture show that population movements within the last 1,500 years are likely to have maintained differentiation among groups. Our analysis provides a timeline of the gene flow events that have generated the contemporary genetic landscape of West Eurasia.

•

Recent admixture events involved outside groups at the edges of West Eurasia

•

Admixture within Europe tended to fall within the European Migration Period

•

West Eurasian genetic structure today is likely to have been maintained by admixture

The Greeks in the West: genetic signatures of the Hellenic colonisation in southern Italy and Sicily

Greek colonisation of South Italy and Sicily (Magna Graecia) was a defining event in European cultural history, although the demographic processes and genetic impacts involved have not been systematically investigated. Here, we combine high-resolution surveys of the variability at the uni-parentally inherited Y chromosome and mitochondrial DNA in selected samples of putative source and recipient populations with forward-in-time simulations of alternative demographic models to detect signatures of that impact. Using a subset of haplotypes chosen to represent historical sources, we recover a clear signature of Greek ancestry in East Sicily compatible with the settlement from Euboea during the Archaic Period (eighth to fifth century BCE). We inferred moderate sex-bias in the numbers of individuals involved in the colonisation: a few thousand breeding men and a few hundred breeding women were the estimated number of migrants. Last, we demonstrate that studies aimed at quantifying Hellenic genetic flow by the proportion of specific lineages surviving in present-day populations may be misleading.

Shared language, diverging genetic histories: high-resolution analysis of Y-chromosome variability in Calabrian and Sicilian Arbereshe

The relationship between genetic and linguistic diversification in human populations has been often explored to interpret some specific issues in human history. The Albanian-speaking minorities of Sicily and Southern Italy (Arbereshe) constitute an important portion of the ethnolinguistic variability of Italy. Their linguistic isolation from neighboring Italian populations and their documented migration history, make such minorities particularly effective for investigating the interplay between cultural, geographic and historical factors. Nevertheless, the extent of Arbereshe genetic relationships with the Balkan homeland and the Italian recipient populations has been only partially investigated. In the present study we address the genetic history of Arbereshe people by combining highly resolved analyses of Y-chromosome lineages and extensive computer simulations. A large set of slow- and fast-evolving molecular markers was typed in different Arbereshe communities from Sicily and Southern Italy (Calabria), as well as in both the putative Balkan source and Italian sink populations. Our results revealed that the considered Arbereshe groups, despite speaking closely related languages and sharing common cultural features, actually experienced diverging genetic histories. The estimated proportions of genetic admixture confirm the tight relationship of Calabrian Arbereshe with modern Albanian populations, in accordance with linguistic hypotheses. On the other hand, population stratification and/or an increased permeability of linguistic and geographic barriers may be hypothesized for Sicilian groups, to account for their partial similarity with Greek populations and their higher levels of local admixture. These processes ultimately resulted in the differential acquisition or preservation of specific paternal lineages by the present-day Arbereshe communities.

Forensic ancestry analysis with two capillary electrophoresis ancestry informative marker (AIM) panels: Results of a collaborative EDNAP exercise

There is increasing interest in forensic ancestry tests, which are part of a growing number of DNA analyses that can enhance routine profiling by obtaining additional genetic information about unidentified DNA donors. Nearly all ancestry tests use single nucleotide polymorphisms (SNPs), but these currently rely on SNaPshot single base extension chemistry that can fail to detect mixed DNA. Insertion-deletion polymorphism (Indel) tests have been developed using dye-labeled primers that allow direct capillary electrophoresis detection of PCR products (PCR-to-CE). PCR-to-CE maintains the direct relationship between input DNA and signal strength as each marker is detected with a single dye, so mixed DNA is more reliably detected. We report the results of a collaborative inter-laboratory exercise of 19 participants (15 from the EDNAP European DNA Profiling group) that assessed a 34-plex SNP test using SNaPshot and a 46-plex Indel test using PCR-to-CE. Laboratories were asked to type five samples with different ancestries and detect an additional mixed DNA sample. Statistical inference of ancestry was made by participants using the Snipper online Bayes analysis portal plus an optional PCA module that analyzes the genotype data alongside calculation of Bayes likelihood ratios. Exercise results indicated consistent genotyping performance from both tests, reaching a particularly high level of reliability for the Indel test. SNP genotyping gave 93.5% concordance (compared to the organizing laboratory's data) that rose to 97.3% excluding one laboratory with a large number of miscalled genotypes. Indel genotyping gave a higher concordance rate of 99.8% and a reduced no-call rate compared to SNP analysis. All participants detected the mixture from their Indel peak height data and successfully assigned the correct ancestry to the other samples using Snipper, with the exception of one laboratory with SNP miscalls that incorrectly assigned ancestry of two samples and did not obtain informative likelihood ratios for a third. Therefore, successful ancestry assignments were achieved by participants in 92 of 95 Snipper analyses. This exercise demonstrates that ancestry inference tests based on binary marker sets can be readily adopted by laboratories that already have well-established CE regimes in place. The Indel test proved to be easy to use and allowed all exercise participants to detect the DNA mixture as well as achieving complete and concordant profiles in nearly all cases. Lastly, two participants successfully ran parallel next-generation sequencing analyses (each using different systems) and achieved high levels of genotyping concordance using the exercise PCR primer mixes unmodified.

Ancient human genomes suggest three ancestral populations for present-day Europeans

We sequenced the genomes of a ~7,000 year old farmer from Germany and eight ~8,000 year old hunter-gatherers from Luxembourg and Sweden. We analyzed these and other ancient genomes^1–4 with 2,345 contemporary humans to show that most present Europeans derive from at least three highly differentiated populations: West European Hunter-Gatherers (WHG), who contributed ancestry to all Europeans but not to Near Easterners; Ancient North Eurasians (ANE) related to Upper Paleolithic Siberians³, who contributed to both Europeans and Near Easterners; and Early European Farmers (EEF), who were mainly of Near Eastern origin but also harbored WHG-related ancestry. We model these populations’ deep relationships and show that EEF had ~44% ancestry from a “Basal Eurasian” population that split prior to the diversification of other non-African lineages.

Static and moving frontiers: the genetic landscape of Southern African Bantu-speaking populations

A consensus on Bantu-speaking populations being genetically similar has emerged in the last few years, but the demographic scenarios associated with their dispersal are still a matter of debate. The frontier model proposed by archeologists postulates different degrees of interaction among incoming agropastoralist and resident foraging groups in the presence of "static" and "moving" frontiers. By combining mitochondrial DNA and Y chromosome data collected from several southern African populations, we show that Bantu-speaking populations from regions characterized by a moving frontier developing after a long-term static frontier have larger hunter-gatherer contributions than groups from areas where a static frontier was not followed by further spatial expansion. Differences in the female and male components suggest that the process of assimilation of the long-term resident groups into agropastoralist societies was gender biased. Our results show that the diffusion of Bantu languages and culture in Southern Africa was a process more complex than previously described and suggest that the admixture dynamics between farmers and foragers played an important role in shaping the current patterns of genetic diversity.

A global analysis of Y-chromosomal haplotype diversity for 23 STR loci

In a worldwide collaborative effort, 19,630 Y-chromosomes were sampled from 129 different populations in 51 countries. These chromosomes were typed for 23 short-tandem repeat (STR) loci (DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS385ab, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635, GATAH4, DYS481, DYS533, DYS549, DYS570, DYS576, and DYS643) and using the PowerPlex Y23 System (PPY23, Promega Corporation, Madison, WI). Locus-specific allelic spectra of these markers were determined and a consistently high level of allelic diversity was observed. A considerable number of null, duplicate and off-ladder alleles were revealed. Standard single-locus and haplotype-based parameters were calculated and compared between subsets of Y-STR markers established for forensic casework. The PPY23 marker set provides substantially stronger discriminatory power than other available kits but at the same time reveals the same general patterns of population structure as other marker sets. A strong correlation was observed between the number of Y-STRs included in a marker set and some of the forensic parameters under study. Interestingly a weak but consistent trend toward smaller genetic distances resulting from larger numbers of markers became apparent.

Reevaluating a model of gender-biased gene flow among Sub-Saharan Hunter-gatherers and farmers

In a previous study, we proposed a model for genetic admixture between African hunter-gatherers and food producers, in which we integrated demographic and genetic aspects together with ethnographic knowledge (Destro-Bisol et al. 2004b). In that study it was possible to test the model only using genetic information from widely dispersed and genetically heterogeneous populations. Here we reevaluate the congruence between the model and patterns of genetic variation using an anthropologically and geographically more homogeneous data set that includes Pygmies and farmers from Cameroon, Congo, and the Central African Republic. As implied by the model, the ratios of mtDNA to Y chromosome Nm estimates (effective population size, N, times the migration rate, m; 0.154 in Pygmies and 6.759 in farmers), support an asymmetric gene flow, with a higher Bantu-to-Pygmy gene flow for paternal than for maternal lineages, and vice versa for farmers. Analyses of intra- and interpopulation genetic variation further support the above observation, showing a prevailing effect of genetic drift on maternal lineages and gene flow on paternal lineages among Pygmies, and an opposite pattern among farmers. We also detected differences between patterns for classical and molecular measures of Y chromosome intrapopulation variation, which likely represent signatures of the introgression of Bantu lineages into the gene pool of Pygmy populations. On the whole, our results seem to reflect differences in the demographic history and the degree of patrilocality and polygyny between the two population groups, thus providing further support to our microevolutionary model in an anthropologically coherent framework.

Linguistic, geographic and genetic isolation: a collaborative study of Italian populations

The animal and plant biodiversity of the Italian territory is known to be one of the richest in the Mediterranean basin and Europe as a whole, but does the genetic diversity of extant human populations show a comparable pattern? According to a number of studies, the genetic structure of Italian populations retains the signatures of complex peopling processes which took place from the Paleolithic to modern era. Although the observed patterns highlight a remarkable degree of genetic heterogeneity, they do not, however, take into account an important source of variation. In fact, Italy is home to numerous ethnolinguistic minorities which have yet to be studied systematically. Due to their difference in geographical origin and demographic history, such groups not only signal the cultural and social diversity of our country, but they are also potential contributors to its bio-anthropological heterogeneity. To fill this gap, research groups from four Italian Universities (Bologna, Cagliari, Pisa and Roma Sapienza) started a collaborative study in 2007, which was funded by the Italian Ministry of Education, University and Research and received partial support by the Istituto Italiano di Antropologia. In this paper, we present an account of the results obtained in the course of this initiative. Four case-studies relative to linguistic minorities from the Eastern Alps, Sardinia, Apennines and Southern Italy are first described and discussed, focusing on their micro-evolutionary and anthropological implications. Thereafter, we present the results of a systematic analysis of the relations between linguistic, geographic and genetic isolation. Integrating the data obtained in the course of the long-term study with literature and unpublished results on Italian populations, we show that a combination of linguistic and geographic factors is probably responsible for the presence of the most robust signatures of genetic isolation. Finally, we evaluate the magnitude of the diversity of Italian populations in the European context. The human genetic diversity of our country was found to be greater than observed throughout the continent at short (0-200 km) and intermediate (700-800km) distances, and accounted for most of the highest values of genetic distances observed at all geographic ranges. Interestingly, an important contribution to this pattern comes from the "linguistic islands"( e.g. German speaking groups of Sappada and Luserna from the Eastern Italian Alps), further proof of the importance of considering social and cultural factors when studying human genetic variation.

Demographic histories, isolation and social factors as determinants of the genetic structure of Alpine linguistic groups

Great European mountain ranges have acted as barriers to gene flow for resident populations since prehistory and have offered a place for the settlement of small, and sometimes culturally diverse, communities. Therefore, the human groups that have settled in these areas are worth exploring as an important potential source of diversity in the genetic structure of European populations. In this study, we present new high resolution data concerning Y chromosomal variation in three distinct Alpine ethno-linguistic groups, Italian, Ladin and German. Combining unpublished and literature data on Y chromosome and mitochondrial variation, we were able to detect different genetic patterns. In fact, within and among population diversity values observed vary across linguistic groups, with German and Italian speakers at the two extremes, and seem to reflect their different demographic histories. Using simulations we inferred that the joint effect of continued genetic isolation and reduced founding group size may explain the apportionment of genetic diversity observed in all groups. Extending the analysis to other continental populations, we observed that the genetic differentiation of Ladins and German speakers from Europeans is comparable or even greater to that observed for well known outliers like Sardinian and Basques. Finally, we found that in south Tyroleans, the social practice of Geschlossener Hof, a hereditary norm which might have favored male dispersal, coincides with a significant intra-group diversity for mtDNA but not for Y chromosome, a genetic pattern which is opposite to those expected among patrilocal populations. Together with previous evidence regarding the possible effects of “local ethnicity” on the genetic structure of German speakers that have settled in the eastern Italian Alps, this finding suggests that taking socio-cultural factors into account together with geographical variables and linguistic diversity may help unveil some yet to be understood aspects of the genetic structure of European populations.

Uniparental markers of contemporary Italian population reveals details on its pre-Roman heritage

Background

According to archaeological records and historical documentation, Italy has been a melting point for populations of different geographical and ethnic matrices. Although Italy has been a favorite subject for numerous population genetic studies, genetic patterns have never been analyzed comprehensively, including uniparental and autosomal markers throughout the country.

Methods/Principal Findings

A total of 583 individuals were sampled from across the Italian Peninsula, from ten distant (if homogeneous by language) ethnic communities — and from two linguistic isolates (Ladins, Grecani Salentini). All samples were first typed for the mitochondrial DNA (mtDNA) control region and selected coding region SNPs (mtSNPs). This data was pooled for analysis with 3,778 mtDNA control-region profiles collected from the literature. Secondly, a set of Y-chromosome SNPs and STRs were also analyzed in 479 individuals together with a panel of autosomal ancestry informative markers (AIMs) from 441 samples. The resulting genetic record reveals clines of genetic frequencies laid according to the latitude slant along continental Italy – probably generated by demographical events dating back to the Neolithic. The Ladins showed distinctive, if more recent structure. The Neolithic contribution was estimated for the Y-chromosome as 14.5% and for mtDNA as 10.5%. Y-chromosome data showed larger differentiation between North, Center and South than mtDNA. AIMs detected a minor sub-Saharan component; this is however higher than for other European non-Mediterranean populations. The same signal of sub-Saharan heritage was also evident in uniparental markers.

Conclusions/Significance

Italy shows patterns of molecular variation mirroring other European countries, although some heterogeneity exists based on different analysis and molecular markers. From North to South, Italy shows clinal patterns that were most likely modulated during Neolithic times.

Reconstructing ancient mitochondrial DNA links between Africa and Europe

Mitochondrial DNA (mtDNA) lineages of macro-haplogroup L (excluding the derived L3 branches M and N) represent the majority of the typical sub-Saharan mtDNA variability. In Europe, these mtDNAs account for <1% of the total but, when analyzed at the level of control region, they show no signals of having evolved within the European continent, an observation that is compatible with a recent arrival from the African continent. To further evaluate this issue, we analyzed 69 mitochondrial genomes belonging to various L sublineages from a wide range of European populations. Phylogeographic analyses showed that ~65% of the European L lineages most likely arrived in rather recent historical times, including the Romanization period, the Arab conquest of the Iberian Peninsula and Sicily, and during the period of the Atlantic slave trade. However, the remaining 35% of L mtDNAs form European-specific subclades, revealing that there was gene flow from sub-Saharan Africa toward Europe as early as 11,000 yr ago.

The peopling of Europe and the cautionary tale of Y chromosome lineage R-M269

Recently, the debate on the origins of the major European Y chromosome haplogroup R1b1b2-M269 has reignited, and opinion has moved away from Palaeolithic origins to the notion of a younger Neolithic spread of these chromosomes from the Near East. Here, we address this debate by investigating frequency patterns and diversity in the largest collection of R1b1b2-M269 chromosomes yet assembled. Our analysis reveals no geographical trends in diversity, in contradiction to expectation under the Neolithic hypothesis, and suggests an alternative explanation for the apparent cline in diversity recently described. We further investigate the young, STR-based time to the most recent common ancestor estimates proposed so far for R-M269-related lineages and find evidence for an appreciable effect of microsatellite choice on age estimates. As a consequence, the existing data and tools are insufficient to make credible estimates for the age of this haplogroup, and conclusions about the timing of its origin and dispersal should be viewed with a large degree of caution.

Mitochondrial haplogroup U5b3: a distant echo of the epipaleolithic in Italy and the legacy of the early Sardinians

There are extensive data indicating that some glacial refuge zones of southern Europe (Franco-Cantabria, Balkans, and Ukraine) were major genetic sources for the human recolonization of the continent at the beginning of the Holocene. Intriguingly, there is no genetic evidence that the refuge area located in the Italian Peninsula contributed to this process. Here we show, through phylogeographic analyses of mitochondrial DNA (mtDNA) variation performed at the highest level of molecular resolution (52 entire mitochondrial genomes), that the most likely homeland for U5b3-a haplogroup present at a very low frequency across Europe-was the Italian Peninsula. In contrast to mtDNA haplogroups that expanded from other refugia, the Holocene expansion of haplogroup U5b3 toward the North was restricted by the Alps and occurred only along the Mediterranean coasts, mainly toward nearby Provence (southern France). From there, approximately 7,000-9,000 years ago, a subclade of this haplogroup moved to Sardinia, possibly as a result of the obsidian trade that linked the two regions, leaving a distinctive signature in the modern people of the island. This scenario strikingly matches the age, distribution, and postulated geographic source of a Sardinian Y chromosome haplogroup (I2a2-M26), a paradigmatic case in the European context of a founder event marking both female and male lineages.

Phylogenetic evidence for multiple independent duplication events at the DYS19 locus

Duplication events at Y chromosome STR loci have been repeatedly described in human populations. DYS19 is probably the best known example and it exhibits duplicate state in individuals from all continents. Despite the large amount of available data, evolutionary relationship between DYS19 duplication-bearing chromosomes has not been so far investigated. We address the genealogical correlation among such chromosomes by analysing newly identified DYS19 duplicated Y chromosomes by SNP genotyping and microsatellite-based network analysis. SNP and network analysis show that DYS19 duplicated Y chromosomes associate with different Y chromosome lineages. These results indicate that DYS19 duplication occurred more than once during human evolution.

A 9-loci Y chromosome haplotype in three Italian populations

Three geographic areas of Italy have been sampled and genotyped for 9 Y chromosome STRs: DYS19, DYS385, DYS388, DYS389 I and II, DYS390, DYS391, DYS392, DYS393. Sampling was focused on residents of small areas, well distant from major urban centres. Only individuals whose grandfather would live in the same area were included. A total of 210 unrelated individuals were collected. Distribution of genetic variation across the three samples and comparison with previously published Italian database indicated that so far Y chromosome diversity has been only partially explored in the Italian Peninsula.