Results of the analysis of mtDNA in 15 individuals of El Aramo mine.
Dynamics of the Neolithic transition across Europe using ancient DNA datasets have established that Neolithic European populations received a limited amount of admixture from resident hunter-gatherers. However, the genetic diversity of Neolithic and Chalcolithic human populations was shaped predominantly by local processes. In the Iberian Peninsula, the Cantabrian fringe showed different proportions of local hunter-gatherers’ ancestry through time. The objective of this chapter is to analyze the mitochondrial variation of populations from the northern Iberian Peninsula from Neolithic to Chalcolithic time using new data from El Aramo mine (Asturias), in the context of the debate about the origin and dispersion of the Beaker culture in Europe.
- northern Iberia
- El Aramo mine
- Bell Beaker culture
1.1. State of the art
Ancient mitochondrial DNA (mtDNA) provides important insights into the movement and spread of human populations. In particular, European populations exhibit some remarkable changes after the end of the last glacial maximum (after 20,000 YBP). The changes in the early postglacial period are thought to be the result of the arrival of new human population groups to Europe . These new populations brought to Europe new mitochondrial DNAs that caused a change in the frequency of the indigenous mtDNA lineages. By reconstructing the variability of the mtDNA of past populations, it is possible to infer population movements that shaped the current genetic variability of our species.
One of the most studied population movements is Neolithization, the transition from a nomadic hunter-gatherer to an agro-pastoralist lifestyle. The debate on the mechanisms of the Neolithic transition has been framed within a dichotomy based on either a demic (DD) or a cultural diffusion (CD). According to the DD model, the migrating people bringing new knowledge experienced some gene flow with the local hunter-gatherer groups. On the other hand, the CD model postulates that the Neolithic transition was mediated mainly through the transmission of the agro-pastoralist system without substantial movement of people .
However, several DNA studies on different ancient European populations indicated a more complex pattern for the Neolithic transition. Unlike the initial proposal based on classical genetic markers that suggested a major migration wave , further studies have shown that the Neolithization process varied in different regions, occurring along several different routes into and across Europe, and having a different genetic impact on the various regions and at various times [4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17].
The mtDNA frequency distribution observed in hunter-gatherers and farmers from Europe provides support for a random dispersion model for Neolithic farmers, with different impacts on the various geographic regions (Central Europe, Mediterranean Europe, and the Cantabrian fringe) .
The transition from the Neolithic to the Chalcolithic period in Europe has been debated. Previous mitochondrial DNA analyses on ancient Europeans have suggested that the current distribution of haplogroup H was modeled by the expansion of the Bell Beaker culture (BBC) out of Iberia during the Chalcolithic period. In addition, it has been suggested that these groups with Bell Beaker (BB) culture in Central Europe represented a population movement from the Iberian Peninsula . However, according to the mtDNA variability in Chalcolithic groups from the Cantabrian fringe of Iberia, no genetic relationships have been detected between these Iberian and Central European groups . This suggestion has been confirmed by the recent study  about the Beaker phenomenon and genomic data of Europe.
1.2. Paleogenetics and paleogenomics
Massive sequencing has allowed important achievements in the field of human evolution, such as the “Neanderthal Genome” project, the discovery of new species (e.g., denisovans), and the recovery of the genome of very ancient humans (remains of La Sima de los Huesos—Atapuerca, Spain—dated to more than 400,000 years BP) [25, 26, 27, 28, 29, 30, 31].
The first paleogenomic studies about modern humans was the 7x coverage genome of the exceptionally well-preserved Tyrolean Ice man, Ötzi, dated to about 5300 years BP . Currently, there are complete genomes from over 90 humans that inhabited Eurasia between 50,000 and 5000 years BP (hunter-gatherers and Neolithic farmers), shedding light on the migratory movements that shaped the genetic variability of modern humans and validating hypotheses proposed from the inference of modern genomes or partial sequences of these individuals [11, 13, 14, 15, 23, 24, 33, 34, 35, 36, 37, 38, 39]. These paleogenomic studies will enlarge the possibilities of selective and demographic analyses of the European prehistoric populations. The genomic data from European hunter-gatherers and farmers show that there is no evidence that the first modern humans in Europe (~45,000–37,000 years ago) contributed to the genetic makeup of current Europeans; these data rather suggest that individuals between ~37,000 and ~14,000 years descended from a single-founder population that is part of the ancestry of today’s Europeans. During the period of greatest warming after ~14,000 years ago, a genetic component related to the inhabitants of the Middle East region became widespread in Europe. These results document how population rotation and migration have been recurring themes of European prehistory .
Recently, 400 European individuals ranging from the Neolithic period to the Bronze Age were analyzed using paleogenomic techniques, including 226 individuals associated to Beaker complex artifacts . Limited genetic affinity between BBC-associated individuals from the Iberian Peninsula and Central Europe was observed, and thus the authors excluded migration as an important mechanism of spread between these regions . This result rejects the hypothesis of the migratory movement of humans from the Iberian Peninsula to Central Europe in the Chalcolithic period accompanied by the BB culture .
In the debate about the biological influence of the dispersion of the Beaker culture in Europe, we have analyzed the mtDNA of remains recovered in El Aramo Mine in Asturias (Cantabrian fringe) from the Late Chalcolithic period that were not accompanied by BB cultural artifacts . This human group is contemporaneous to other Iberian Chalcolithic populations both without Beaker complex artifacts associated and with Beaker culture associated. Sites without BBC associated are those of Longar and San Jaun Ante Portam Latinam (SJAPL) in the Basque Country . Contemporaneous sites with Beaker complex culture associated are the central and southern Iberian and central European groups published by [18, 40]. The aim of this study is to contribute new mtDNA data variability of the Chalcolithic site from El Aramo Mine (Asturias) and to determine whether there is either a common genetic signal or a heterogeneous genetic landscape among Chalcolithic European groups (with and without BBC culture).
2. Material and methods
In this chapter, we have analyzed the human remains from El Aramo Mine discovered in 1888, a mine located in the Asturias region in the Cantabrian fringe of the Iberian Peninsula . The direct 14C analysis of the human remains from this mine indicated a dating between the Late Chalcolithic period and the Early Bronze Age. The anthropological remains from El Aramo Mine consists of 9 skulls and 12 skeletal remains. We have isolated DNA mainly from dental pieces (since it is the material that offers the greatest guarantees when recovering DNA). However, in some cases, we had to pulverize bone remains in order to perform DNA extraction, since it was the only anthropological material available.
In the case of teeth, we have selected those without caries or deep fissures that might extend into the dental pulp. The surface of the teeth was thoroughly cleaned with acids and ultraviolet (UV) irradiation to eliminate any possible DNA contaminants . In the case of bones, the surface was thoroughly cleaned by abrasion and pulverized using a Freezer miller. Then we extracted DNA from bone and dental tissue by means of the phenol/chloroform method with some modifications [20, 21, 22, 43].
The sequencing of a 399 bp (nps 16,000–16,399) segment of HVS-I and 394 bp (nps1–394) of HVS-II of the mtDNA as per  was conducted by amplifying 6 overlapping fragments of 93–133 bp in length. The protocol followed and the primers used are described in [9, 45]. Likewise, in order to verify the obtained mtDNA haplogroups, the nucleotide position of the coding region of mtDNA was determined by means of PCR-restriction fragment length polymorphisms (RFLPs) [43, 46].
The extraction of DNA and the preparation of samples for PCR were performed in a sterile chamber with positive pressure, free of modern DNA, in which no post-PCR process had ever been carried out. Ancient DNA results were validated through the application of the following criteria [47, 48]: (1) suitable clothing was used (disposable cap, gloves, mask and laboratory coat), (2) controls were applied to detect contamination during the extraction process and in each one of the amplifications, (3) Real-time PCR quantification of amplifiable DNA to quantify one mtDNA fragment of 113 bp was conducted [9, 49], (4) a duplicate analysis was performed on the greatest possible number of individuals, and (5) Cloning of PCR products was performed with subsequent sequencing of the clones. The cloning was carried out using TOPO TA Cloning® Kits (Invitrogen), following the supplier’s instructions.
The mitochondrial variability resulted from El Aramo Mine was compared with other ancient and present-day populations. With respect to hunter-gatherers, three groups were considered: one from Scandinavia, one from Central Europe [13, 14, 50, 51], and one from the Cantabrian fringe of the Iberian Peninsula [9, 17, 33, 52]. Regarding the Neolithic DNA, 14 populations were selected: 3 from the Near East , 4 from Central and Eastern Europe [16, 45], 5 from the Mediterranean area of Europe (Hungary, Romania, Catalonia and France) [6, 7, 10, 12], and 2 from northern Iberia [9, 11]. With regard to the Chalcolithic groups, we considered one from Central Europe with BB artifacts associated [16, 18], one from the Cantabrian fringe of Iberia without BC culture (Longar and SJAPL sites)  and another two from Iberia, one with BB culture, and another one without BB culture [BBC: Arroyal (Burgos), Camino de las Yeseras (Madrid), Humanajes (Madrid), La Magdalena (Madrid), and Paris Street (Barcelona). Without BBC: Camino del Molino (Murcia), Bolares (Extremadura), el Sotillo, chabola de la Hechicera (Alava), el Mirador (Burgos), La Mina, Trocs (Huesca), and El Portalón (Burgos)] . The Bronze Age period is represented by three groups from Siberia, Kazakhstan, and Romania [10, 53, 54]. The present-day populations database corresponds to that described in , to which the present-day population of Asturias, where El Aramo Mine is located, has been added [55, 56].
The FST distance matrix between present and ancient populations was calculated from the mtDNA haplogroup frequencies using Arlequin 3.11 . Relationships between populations were studied through Multidimensional Scaling analysis (MDS), based on the FST distance matrix, using SPSS 20 Software. Median Joining Network (MJN) for certain haplogroups was generated to infer phylogenetic relationships between the mitochondrial lineages from the Paleolithic period to the present day using Network software v22.214.171.124 (available at http://www.fluxus-engineering.com). Different mutation weights were applied in accordance with previous papers [58, 59, 60], and the point insertions and deletions were excluded from the analysis.
3. Results and discussion
Mitochondrial DNA variability was analyzed in 21 skeletal remains recovered from El Aramo Mine, in the Asturias region (Cantabrian fringe of the Iberian Peninsula). The quantification of the template mtDNA number of each of the samples showed values above 1000 molecules/μl in all the DNA extracts from teeth and values below 1000 molecules/μl in all the DNA extracts from bones (Table 1). These results indicate the greater efficiency of the DNA extraction from teeth when compared to that from bones. Furthermore, in order to authenticate the results, 5.26% of the samples analyzed were duplicated, and these results were consistent in all the samples (Table 1). Moreover, 24 PCR products were cloned, estimating an average of 7.46 mutations per cloned fragment (~100 bp). These mutations have been interpreted as artifacts produced by the postmortem damage of aDNA.
|Skull 1 (bone)||AR13||854.6||rCRS||73A||H|
|Skull 2 (bone)||AR15||551.4||rCRS||73A-263G-315.1C||H|
|Skull 3 (tooth)d||AR22||89190.2||183C-189C-270T||73G||U5b1b|
|Skull 4 (tooth)||AR21||33668.7||189C||73A||H|
|Skull 5 (tooth)||AR31||17194.3||rCRS||73A-150T-263-309.1C-315.1C-320C||H|
|Skull 6 (tooth)||AR32||18706.4||51G||73A||H1|
|Skull 7 (tooth)d||AR42||2078.3||rCRS||73A-153T-263G-315.1C||H|
|Skull 8 (tooth)d||AR41||5847.2||126C-292T-294T-296T||73G||T2c|
|Skull 9 (tooth)||AR5||178331.7||rCRS||73A-263G-315.1C||H|
|Right Hemimandible (tooth)||AR7-9||153133.6||86C-129A-223T||73G||I3|
|Left Hemimandible (tooth)||AR6||141554.8||69T-126C-320T-360T||73G||J1c|
|Right femur (bone)||ARH10||372.5||126C-163G-187T-189T-294T||73G||T1|
The mitochondrial variability obtained from the 21 human remains from El Aramo Mine showed 15 different haplotypes (genetic diversity: 0.9608 ± 0.0394). The nine skulls studied presented nine different mitochondrial haplotypes, which allow us to rule out the existence of maternal kinship among these individuals. The 12 postcranial remains analyzed showed 7 different mitochondrial haplotypes, which, compared with the haplotypes of the skulls, lead to reject possible coincidences, since the postcranial remains were not associated with the skulls. Finally, the minimum estimated number of individuals was 15, with 15 different mitochondrial haplotypes described (Table 1). The high genetic diversity obtained in El Aramo site allows us to indicate that it is a representative sample of the original population, with no evidence of kinship among these individuals.
The 15 mitochondrial haplotypes obtained from El Aramo Mine were classified into 5 different mitochondrial haplogroups (H, T, J, U5b, and I3), with a genetic diversity of 0.6381 ± 0.1288 and a heterogeneous distribution of their frequency values (60, 13, 13, 7, and 7%, respectively). Haplogroup H is the most frequent one in the population of El Aramo (60%), whose value is close to that shown by the current population of the Asturias region (56%), where El Aramo Mine is located, and much higher than the average value found in European (45%) and Near Eastern (16%) populations [54, 55, 61] (Figure 1). In El Aramo, haplogroup H is represented by nine different mitochondrial haplotypes, five of which correspond to haplogroup H-rCRS (33%), which is the most frequent one in the present-day European population (15–30%) and in the current population of Asturias ( ̴23.5%). The rest of the haplogroups detected in El Aramo mine have low frequencies in the extant populations of Europe and the Near East [54, 55, 61]. Haplogroup T has a frequency of 13% in El Aramo, which is close to that of the present-day population of Asturias (11.2%) and higher than that of the populations of Europe (6.8%) and the Near East (9.91%) (Figure 1). Haplogroup J shows similar frequencies in El Aramo and in the extant population of the Near East (13%), and lower frequencies in the present-day populations of Asturias and Europe (10.5 and 8.6%, respectively). Haplogroup U5b in El Aramo Mine has a frequency of 7%, which is close to that observed in the extant population of Europe (6%) and considerably higher than that detected in the current populations of Asturias and the Near East (0.7 and 0.05%, respectively) (Figure 1). Lastly, haplogroup I3 shows in El Aramo a value (7%) that is considerably higher than that observed in the current populations, among which the population of the Near East (3%) shows the highest frequency value for this haplogroup, which is absent in the current population of Asturias and very infrequent in Europe (1.8%) (Figure 1). To sum up, the main differences in the frequencies of the mitochondrial lineages of El Aramo Mine and the present-day populations lie in haplogroups H, U5b, and I3.
The frequency distribution of the mitochondrial haplogroups of El Aramo Mine was compared with that of other Chalcolithic groups of the Iberian Peninsula and Central Europe [9, 18]. For this analysis, five population groups were defined: (1) El Aramo Mine (Late Chalcolithic-Early Bronze Age) without BBC associated (ARAMO), (2) Chalcolithic from the Cantabrian fringe of the Iberian Peninsula without BBC (CA_Cantabrian), (3) Chalcolithic from Iberia without BBC (CA_Iberian), (4) Chalcolithic from Iberia associated with BBC (BBC_Iberian), and (5) Chalcolithic from Central Europe associated with BBC (BBC_CE) (Figure 2). All the populations included in this analysis share the presence of haplogroups H and T, although the distribution of their frequencies is different. The group of El Aramo shows the highest frequency value for haplogroup H (60%), followed by the two Chalcolithic groups of the Iberian Peninsula without BBC, CA_Cantabrian and CA_Iberian (39 and 31.2%, respectively); the lowest frequencies of haplogroup H were found in the BBC-associated groups of the Iberian Peninsula (BBC_Iberian, with 22%) and Central Europe (BBC_CE, with 29%). This frequency distribution of haplogroup H does not support the hypothesis of Brandt et al. , who suggest that the BBC was spread from the Iberian Peninsula toward Central Europe by haplogroup H carriers, since the highest frequency values were not found in population associated with the BB culture; moreover, the group of the Iberian Peninsula (BBC_Iberian) has the lowest value (22%) (Figure 2).
With respect to haplogroup T, there is also a heterogeneous distribution among the five population groups analyzed. The group of El Aramo shows the highest value for this haplogroup (13%); the Chalcolithic groups associated with BBC (BBC_Iberian and BBC_CE) have frequencies of 5.5 and 7.8%, respectively, which are lower than that of CA_Cantabrian (9.5%) and higher than that of CA_Iberian (2.08%). Given the fact that haplogroup T has been proposed as a marker of the diffusion of the Neolithic culture, the heterogeneity of the frequencies of this haplogroup in the Chalcolithic populations supports the model of random Neolithic cultural diffusion (Hervella et al. ). This behavior has also been observed in the MJN (Figure 3), where T haplotypes of prehistoric and present-day populations were shown. In MJN main nodes are shared by all the populations compared, both prehistoric and modern. The polymorphisms that define the two T haplotypes of El Aramo Mine are shared with Neolithic and Chalcolithic groups, indicating their relation. Furthermore, the polymorphisms of one of the lineages found in El Aramo Mine are shared with current samples of Europe and the Near East, showing its prevalence to the present time (Figure 3).
Haplogroup J shows high-frequency values in the groups that are not associated with BBC (Aramo (13%), CA_Cantabrian (13%), and CA_Iberian (23%)), with the sample of CA_Iberian showing the highest value. The BBC-associated groups have lower frequencies, with 7.8% in the sample of BBC_CE and 0% in BBC_Iberian, showing once again the difference between the Chalcolithic groups with BBC and without BBC. Haplogroup J has also been proposed as a marker of the Neolithic diffusion in Europe. The distribution of frequencies of this haplogroup in the populations analyzed (Figure 2) shows a greater Neolithic influence in the non-BBC Chalcolithic groups. This influence of Neolithic diffusion is also observed in the MJN analysis of haplogroup J (Figure 4), which shows that the haplotypes of the Chalcolithic samples non- BBC are included in the mitochondrial variability of the Near East. With regard to the BBC-associated populations, the J lineages correspond solely to the BBC_CE population, and not to the BBC_Iberian population (Figures 2 and 4), highlighting the differentiation between the two groups with BB culture (BBC_CE and BBC_Iberian). The frequency of haplogroup U5b in El Aramo (7%) indicates the persistence of Paleolithic lineages even after the Neolithic period. The prevalence of this lineage seems to be higher in the Chalcolithic groups of the Iberian Peninsula (16.6% in CA_Iberian and 27.8% in BBC_Iberian) when compared to those of Central Europe (6.3%) (Figure 2). The absence of haplogroup U5b in the Chalcolithic sample of the Cantabrian region (CA_Cantabrian) is due to the lack of differentiation in the U lineages published in this sample .
Lastly, haplogroup I3 appears only in the sample of El Aramo (7%) and in the BBC population of Central Europe (5.5%). The origin of this mitochondrial lineage has been thoroughly debated, and the consensus reached is that it began in the recent Neolithic period , with very low-frequency values found in the modern populations. In the MJN conducted for haplogroup I3, it was observed that El Aramo shares polymorphisms with Neolithic and Chalcolithic groups and with the current populations of Europe and the Near East, showing a low diversity for this lineage (Figure 5).
After the comparative analysis of the frequency distribution of the mitochondrial haplogroups of El Aramo and other prehistoric populations, it is possible to distinguish some general tendencies. On the one hand, there is a clear genetic differentiation between the Chalcolithic groups with BB culture (BBC_CE and BBC_Iberian) and those without BB culture (CA_Iberian, CA_Cantabrian and El Aramo); on the other hand, there is a differentiation between the BBC groups of Central Europe and those of the Iberian Peninsula (Figure 2). The population of El Aramo Mine shows a distribution of mitochondrial haplotypes that is closer to that of Chalcolithic populations without BBC than to that of chalcolithic populations with BBC, highlighting the specificity of El Aramo, the high frequency of haplogroup H (60%) and the presence of haplogroup I3, whose frequency is, in addition, one of the highest ones described to the present day.
The mitochondrial variability of the population of El Aramo was analyzed in the context of other prehistoric groups and the modern populations of Europe and the Near East through a Multidimensional Scaling analysis (MDS) (Figure 6). The MDS was done through a matrix of FST distances, calculated by the frequencies of the mitochondrial haplogroups of hunter-gatherer, farmer, Chalcolithic, Bronze Age, and present-day populations from the Cantabrian fringe, Europe, and the Near East. This analysis shows the differentiation between the two hunter-gatherer populations from Central Europe and Scandinavia (Figure 6). It has been proposed that the mtDNA variation of these groups indicates a genetic discontinuity between the hunter-gatherers and later populations in these two geographic regions [13, 14, 50, 51]. However, this suggested discontinuity is not so obvious in the case of the hunter-gatherers from the Cantabrian fringe who separated from those of Central Europe and Scandinavia in this analysis [9, 17] (Figure 6).
With regard to the European Neolithic populations, the heterogeneity observed in the mtDNA haplogroup frequency variation is revealed by their position on the two-dimensional plot of the MDS analysis (Figure 6). On the one hand, a group of populations (Near East, Central Europe, Hungary, and Eastern Pyrenees) with high-frequency values for haplogroup N is separated from the other Neolithic populations. On the other hand, heterogeneity is also apparent within the Mediterranean area, with a Neolithic population of Southern France being closer to present-day populations in the Near East, due to its high frequency for haplogroups J and U, whereas the Neolithic populations from the Iberian Peninsula (Catalonia and the Cantabrian fringe) show lower frequency for those haplogroups (J and U) (Figure 6). The genetic distances observed between the European Neolithic groups suggest a different genetic impact of the Neolithic farmers from the Near East on Central Europe, Mediterranean Europe, and the Cantabrian fringe. These data support a random dispersion model for Neolithic farmers, with different impact on the various geographic regions .
With respect to the Chalcolithic prehistoric groups included in the analysis (BBC_CE, BBC_Iberian, CA_Iberian, CA_Cantabrian and Aramo), those with BB culture (BBC_CE and BBC_Iberian) are differentiated from the rest. On the other hand, the distance between the BBC groups of Central Europe and the Iberian Peninsula is due to both the greater persistence of Paleolithic U5 lineages in the BBC_Iberian group and the higher frequency of Neolithic lineages T and J in the BBC_CE group. This differentiation suggests that the relation between these two Chalcolithic groups (BBC_CE and BBC_Iberian) is only cultural but not genetic, supporting the study  about the Beaker phenomenon and genomic data from Europe, who reject the hypothesis that the genetic substrate of the BBC_CE groups came from the BBC_Iberian groups .
Regarding the Chalcolithic populations without BB culture, CA_Iberian and CA_Cantabrian are genetically close to one another, with El Aramo being further from them due to the high frequencies of haplogroups H and I3. These Chalcolithic groups are not distant from their Neolithic ancestor populations, although they are distant from the present-day populations of these regions (Figure 6), which could be attributed to a post-Neolithic population restructuring .
In view of the results obtained, it can be inferred that the influence of the Neolithic period on the local groups was complex, and its result could generate the existence of Chalcolithic groups in the Iberian Peninsula with genetic and cultural differences, with the latter being mainly related to the Beaker phenomenon.
The human group of El Aramo Mine, without artifacts associated with BBC, shows some peculiarities. Its chronology expands from the Chalcolithic period to the Early Bronze Age, during which a change in the landscape and the subsistence conditions have been detected in this region . El Aramo consists of males who had a strong relationship with mining, and the burials reflect a ritual related to the exploitation of the mine . Therefore, it is a group with differentiated cultural characteristics within the Chalcolithic groups of the Iberian Peninsula. This cultural differentiation seems to be accompanied by a genetic differentiation, since in the MDS, El Aramo is distant from the rest of the Chalcolithic groups (Figure 6). These results indicate the existence of local processes in the Chalcolithic period that could be related to subsistence strategies linked to the characteristics of the environment. In the case of El Aramo (Asturias), these strategies are linked to mining activity, since the region of Asturias is one of the traditional mining areas since the Neolithic period .
This work was supported by the Basque Government to Research Groups of the Basque University System (IT1138-16) and grant from the Spanish Ministry of Science and Innovation (GCL2016-79093/P). We are grateful to Miguel Angel de Blas Cortina, Director of the archeological intervention in El Aramo Mine, for providing archeological data and fruitful discussion.
Conflict of interest
The authors have declared that no competing interests exist.
Revesz PZ. A spatio-temporal analysis of mitochondrial DNA Haplogroup I. MATEC Web Conference. 2016; 76:04048. DOI: 10.1051/matecconf/2016760
Zvelebil M, Zvelebil KV. Agricultural transition and indo-European dispersal. Antiquity. 1988; 62:574-583
Ammerman AJ, Cavalli-Sforza LL. The Neolithic Transition and the Genetics of Populations in Europe. Princeton: Princeton University Press; 1984
Haak W, Forster P, Bramanti B, Matsumura S, Brandt G, Tanzer M, Villems R, Renfrew C, Gronenborn D, Alt KW, Burger J. Ancient DNA from the first European farmers in 7500-year-old Neolithic sites. Science. 2005; 310:1016-1018. DOI: 10.1126/science.1118725
Haak W, Balanovsky O, Sanchez JJ, Koshel S, Zaporozhchenko V, Christina J, Adler C, Der Sarkissian SI, Brandt G, Schwarz C, Nicklisch N, Dresely V, Fritsch B, Balanovska E, Villems R, Meller H, Alt K, Cooper A. Genographic consortium. Ancient DNA from European early Neolithic farmers reveals their near eastern affinities. PLoS One Biology. 2010; 8:e1000536. DOI: 10.1371/journal.pbio.1000536
Sampietro ML, Lao O, Caramelli DL, Pou R, Marti M, Bertranpetit J, Lalueza-Fox C. Palaeogenetic evidence supports a dual model of Neolithic spreading into Europe. Philosophical Transactions of the Royal Society B: Biological Sciences. 2007; 274:2161-2167. DOI: 10.1098/rspb.2007.0465
Lacan M, Keyser C, Ricaut FX, Brucato N, Duranthon F, Duranthon F, Guilaine J, Eric Crubézy E, Ludes B. Ancient DNA reveals male diffusion through the Neolithic Mediterranean route. Proceedings of the National Academy of Science of the USA. 2011; 108:9788-9791. DOI: 10.1073/pnas.1100723108
Deguilloux MF, Soler L, Pemonge MH, Scarre C, Joussaume R, Laporte L. News from the west: Ancient DNA from a French megalithic burial chamber. American Journal of Physical Anthropology. 2011; 144:108-118. DOI: 10.1002/ajpa.21376
Hervella M, Izagirre N, Alonso S, Fregel R, Alonso A, Cabrera VM, de-la-Rua C. Ancient DNA from hunter-gatherer and farmer groups from northern Spain supports a random dispersion model for the Neolithic expansion into Europe. PLoS One. 2012; 7:e34417. DOI: 10.1371/journal.pone.0034417
Hervella M, Rotea M, Izagirre N, Constantinescu M, Alonso S, Ioana M, Lazar C, Ridichi F, Soficaru AD, Netea MG, de-la-Rua C. Ancient DNA from south-East Europe reveals different events during early and middle Neolithic influencing the European genetic heritage. PLoS One. 2015; 10:e0128810. DOI: 10.1371/journal.pone.0128810
Gamba C, Fernández E, Tirado M, Deguilloux MF, Pemonge MH, Utrilla P, Edo M, Molist M, Rasteiro R, Chikhi L, Arroyo-Pardo E. Ancient DNA from an early Neolithic Iberian population supports a pioneer colonization by first farmers. Molecular Ecology. 2012; 21:45-56. DOI: 10.1111/j.1365-294X.2011.05361.x
Gamba C, Jones ER, Teasdale MD, McLaughlin RL, Gonzalez-Fortes G, Mattiangeli V, Domboróczki L, Kővári I, Pap I, Anders A, Whittle A, János D, Raczky P, Higham TFG, Hofreiter M, Bradley DG, Pinhasi R. Genome flux and stasis in a five millennium transect of European prehistory. Nature Communications. 2014; 5:5257. DOI: 10.1038/ncomms6257
Skoglund P, Malmström H, Raghavan M, Stora J, Hall P, Willerslev E, Gilbert MT, Götherström A, Jakobsson M. Origins and genetic legacy of Neolithic farmers and hunter-gatherers in Europe. Science. 2012; 336:466-469. DOI: 10.1126/science.1216304
Skoglund P, Malmström H, Omrak A, Raghavan M, Valdiosera V, Günther T, Hall P, Tambets K, Parik J, Stören KG, Apel J, Hall P, Willerslev E, Stora J, Götherström A, Jakobsson M. Genomic diversity and admixture differs for stone-age Scandinavian foragers and farmers. Science. 2014; 344:747-750. DOI: 10.1126/science.1253448
Fernández E, Pérez-Pérez A, Gamba C, Prats E, Cuesta P, Cuesta P, Anfruns J, Molist M, Arroyo-Pardo E, Turbón D. Ancient DNA analysis of 8000 B.C. near eastern farmers supports an early Neolithic Pioneer maritime colonization of mainland Europe through Cyprus and the Aegean Islands. PLoS Genetics. 2014; 10:e1004401. DOI: 10.1371/journal.pgen.1004401
Brandt G, Haak W, Adler CJ, Roth C, Szécsényi-Nagy A, Karimnia S, Möller-Rieker S, Meller H, Ganslmeier R, Friederich S, Dresely V, Nicklisch N, Pickrell JK, Sirocko F, Reich D, Cooper A, Alt KW, Genographic Consortium. Ancient DNA reveals key stages in the formation of central European mitochondrial genetic diversity. Science. 2013; 342:257-261. DOI: 10.1126/science.1241844
de-la-Rua C, Izagirre N, Alonso S, Hervella M. Ancient DNA in the Cantabrian fringe populations: A mtDNA study from prehistory to late antiquity. Quaternary International. 2015; 364:306-311. DOI: 10.1016/j.quaint.2015.01.035
Olalde I, Brace S, Allentoft ME, Armit I, Kristiansen K, Booth T, Rohland N, Mallick S, Szécsényi-Nagy A, Mittnik A, Altena E, Lipson M, Lazaridis I, Harper TK, Patterson N, Broomandkhoshbacht N, Diekmann Y, Faltyskova Z, Fernandes D, Ferry M, Harney E, de Knijff P, Michel M, Oppenheimer J, Stewardson K, Barclay A, Alt KW, Liesau C, Ríos P, Blasco C, Miguel JV, García RM, Fernández AA, Bánffy E, Bernabò-Brea M, Billoin D, Bonsall C, Bonsall L, Allen T, Büster L, Carver S, Navarro LC, Craig OE, Cook GT, Cunliffe B, Denaire A, Dinwiddy KE, Dodwell N, Ernée M, Evans C, Kuchařík M, Farré JF, Fowler C, Gazenbeek M, Pena RG, Haber-Uriarte M, Haduch E, Hey G, Jowett N, Knowles T, Massy K, Pfrengle S, Lefranc P, Lemercier O, Lefebvre A, Martínez CH, Olmo VG, Ramírez AB, Maurandi JL, Majó T, JI MK, Mc Sweeney K, Mende BG, Mod A, Kulcsár G, Kiss V, Czene A, Patay R, Endrődi A, Köhler K, Hajdu T, Szeniczey T, Dani J, Bernert Z, Hoole M, Cheronet O, Keating D, Velemínský P, Dobeš M, Candilio F, Brown F, Fernández RF, Herrero-Corral AM, Tusa S, Carnieri E, Lentini L, Valenti A, Zanini A, Waddington C, Delibes G, Guerra-Doce E, Neil B, Brittain M, Luke M, Mortimer R, Desideri J, Besse M, Brücken G, Furmanek M, Hałuszko A, Mackiewicz M, Rapiński A, Leach S, Soriano I, Lillios KT, Cardoso JL, Pearson MP, Włodarczak P, Price TD, Prieto P, Rey PJ, Risch R, Rojo Guerra MA, Schmitt A, Serralongue J, Silva AM, Smrčka V, Vergnaud L, Zilhão J, Caramelli D, Higham T, Thomas MG, Kennett DJ, Fokkens H, Heyd V, Sheridan A, Sjögren KG, Stockhammer PW, Krause J, Pinhasi R, Haak W, Barnes I, Lalueza-Fox C, Reich D. The beaker phenomenon and the genomic transformation of Northwest Europe. Nature. 2018; 555:190-196. DOI: 10.1038/nature25738
Hervella M, Plantinga TS, Alonso S, Ferwerda B, Izagirre N, Fontecha L, Fregel R, van de Meer JWM, de-la-Rua C, Netea MG. The loss of functional Caspase-12 in Europe is a pre-Neolithic event. PLoS One. 2012; 7:e37022. DOI: 10.1371/journal.pone.0037022
Alzualde A, Izagirre N, Alonso S, Alonso A, de-la-Rua C. Temporal mitochondrial DNA variation in the Basque Country: Influence of post-Neolithic events. Annals of Human Genetics. 2005; 69:665-679
Alzualde A, Izagirre N, Alonso S, Alonso A, Albarrán C, Azcarate A, de-la-Rúa C. Insights into the “isolation” of the Basques: mtDNA lineages from the historical site of Aldaieta (6th-7th centuries AD). American Journal Physical Anthropology. 2006; 130:394-404. DOI: 10.1002/ajpa.20375
Alzualde A, Izagirre N, Alonso S, Alonso A, Albarrán C, Azcarate A, de-la-Rua C. Influence of the European kingdoms of late antiquity on the Basque Country: An ancient DNA study. Current Anthropology. 2007; 48:155-163. DOI: 10.1086/510464
Fu Q, Posth C, Hajdinjak M, Petr M, Mallick S, Fernandes D, Furtwängler A, Haak W, Meyer M, Mittnik A, Nickel B, Peltzer A, Rohland N, Slon V, Talamo S, Lazaridis I, Lipson M, Mathieson I, Schiffels S, Skoglund P, Derevianko AP, Drozdov N, Slavinsky V, Tsybankov A, Cremonesi RG, Mallegni F, Gély B, Vacca E, Morales MR, Straus LG, Neugebauer-Maresch C, Teschler-Nicola M, Constantin S, Moldovan OT, Benazzi S, Peresani M, Coppola D, Lari M, Ricci S, Ronchitelli A, Valentin F, Thevenet C, Wehrberger K, Grigorescu D, Rougier H, Crevecoeur I, Flas D, Semal P, Mannino MA, Cupillard C, Bocherens H, Conard NJ, Harvati K, Moiseyev V, Drucker DG, Svoboda J, Richards MP, Caramelli D, Pinhasi R, Kelso J, Patterson N, Krause J, Pääbo S, Reich D. The genetic history of Ice Age Europe. Nature. 2016; 534:200-205. DOI: 10.1038/nature17993
Lazaridis I, Nadel D, Rollefson G, Merrett DC, Rohland N, Mallick S, Fernandes D, Novak M, Gamarra B, Sirak K, Connell S, Stewardson K, Harney E, Fu Q, Gonzalez-Fortes G, Jones ER, Roodenberg SA, Lengyel G, Bocquentin F, Gasparian B, Monge JM, Gregg M, Eshed V, Mizrahi AS, Meiklejohn C, Gerritsen F, Bejenaru L, Blüher M, Campbell A, Cavalleri G, Comas D, Froguel P, Gilbert E, Kerr SM, Kovacs P, Krause J, McGettigan D, Merrigan M, Merriwether DA, O'Reilly S, Richards MB, Semino O, Shamoon-Pour M, Stefanescu G, Stumvoll M, Tönjes A, Torroni A, Wilson JF, Yengo L, Hovhannisyan NA, Patterson N, Pinhasi R, Reich D. Genomic insights into the origin of farming in the ancient Near East. Nature. 2016; 536:419-424. DOI: 10.1038/nature19310
Green RE, Krause J, Ptak SE, Briggs AW, Ronan MT, Simons JF, Du L, Egholm M, Rothberg JM, Paunovic M, Pääbo S. Analysis of one million base pairs of Neanderthal DNA. Nature. 2006; 444:330-336. DOI: 10.1038/nature05336
Green RE, Krause J, Briggs AW, Maricic T, Stenzel U, Kircher M, Patterson N, Li H, Zhai W, Fritz MH, Hansen NF, Durand EY, Malaspinas AS, Jensen JD, Marques-Bonet T, Alkan C, Prüfer K, Meyer M, Burbano HA, Good JM, Schultz R, Aximu-Petri A, Butthof A, Höber B, Höffner B, Siegemund M, Weihmann A, Nusbaum C, Lander ES, Russ C, Novod N, Affourtit J, Egholm M, Verna C, Rudan P, Brajkovic D, Kucan Ž, Gušic I, Doronichev VB, Golovanova LV, Lalueza-Fox C, de la Rasilla M, Fortea J, Rosas A, Schmitz RW, Johnson PLF, Eichler EE, Falush D, Birney E, Mullikin JC, Slatkin M, Nielsen R, Kelso J, Lachmann M, Reich D, Pääbo S. A draft sequence of the Neandertal genome. Science. 2010; 328:710-722. DOI: 10.1126/science.1188021
Noonan JP, Coop G, Kudaravalli S, Smith D, Krause J, Alessi J, Chen F, Platt D, Pääbo S, Pritchard JK, Rubin EM. Sequencing and analysis of Neanderthal genomic DNA. Science. 2006; 314:1113-1118. DOI: 10.1126/science.1131412
Prüfer K, Racimo F, Patterson N, Jay F, Sankararaman S, Sawyer S, Heinze A, Renaud G, Sudmant PH, de Filippo C, Li H, Mallick S, Dannemann M, Fu Q, Kircher M, Kuhlwilm M, Lachmann M, Meyer M, Ongyerth M, Siebauer M, Theunert C, Tandon A, Moorjani P, Pickrell J, Mullikin JC, Vohr SH, Green RE, Hellmann I, Johnson PL, Blanche H, Cann H, Kitzman JO, Shendure J, Eichler EE, Lein ES, Bakken TE, Golovanova LV, Doronichev VB, Shunkov MV, Derevianko AP, Viola B, Slatkin M, Reich D, Kelso J, Pääbo S. The complete genome sequence of a Neanderthal from the Altai Mountains. Nature. 2014; 505:43-49. DOI: 10.1038/nature12886
Reich D, Green RE, Kircher M, Krause J, Patterson N, Durand EY, Viola B, Briggs AW, Stenzel U, Johnson PL, Maricic T, Good JM, Marques-Bonet T, Alkan C, Fu Q, Mallick S, Li H, Meyer M, Eichler EE, Stoneking M, Richards M, Talamo S, Shunkov MV, Derevianko AP, Hublin JJ, Kelso J, Slatkin M, Pääbo S. Genetic history of an archaic hominin group from Denisova Cave in Siberia. Nature. 2010; 468:1053-1060. DOI: 10.1038/nature09710
Meyer M, Kircher M, Gansauge MT, Li H, Racimo F, Mallick S, Schraiber JG, Jay F, Prüfer K, de Filippo C, Sudmant PH, Alkan C, Fu Q, Do R, Rohland N, Tandon A, Siebauer M, Green RE, Bryc K, Briggs AW, Stenzel U, Dabney J, Shendure J, Kitzman J, Hammer MF, Shunkov MV, Derevianko AP, Patterson N, Andrés AM, Eichler EE, Slatkin M, Reich D, Kelso J, Pääbo S. A high-coverage genome sequence from an archaic Denisovan individual. Science. 2012; 338:222-226. DOI: 10.1126/science.1224344
Meyer M, Fu Q, Aximu-Petri A, Glocke I, Nickel B, Arsuaga JL, Martínez I, Gracia A, Bermúdez de Castro JM, Carbonell E, Pääbo S. A mitochondrial genome sequence of a hominin from Sima de los Huesos. Nature. 2014; 505:403-406. DOI: 10.1038/nature12788
Keller A, Graefen A, Ball M, Matzas M, Boisguerin V, Maixner F, Leidinger P, Backes C, Khairat R, Forster M, Stade B, Franke A, Mayer J, Spangler J, McLaughlin S, Shah M, Lee C, Harkins TT, Sartori A, Moreno-Estrada A, Henn B, Sikora M, Semino O, Chiaroni J, Rootsi S, Myres NM, Cabrera VM, Underhill PA, Bustamante CD, Vigl EE, Samadelli M, Cipollini G, Haas J, Katus H, O'Connor BD, Carlson MR, Meder B, Blin N, Meese E, Pusch CM, Zink A. New insights into the Tyrolean Iceman's origin and phenotype as inferred by whole-genome sequencing. Nature Communications. 2012; 3:698. DOI: 10.1038/ncomms1701
Sanchez-Quinto F, Schroeder H, Ramirez O, Avila-Arcos MC, Pybus M, Olalde I, Velazquez AM, Marcos ME, Encinas JM, Bertranpetit J, Orlando L, Gilbert MT, Lalueza-Fox C. Genomic affinities of two 7,000-year-old Iberian hunter-gatherers. Current Biology. 2012; 22:1494-1499. DOI: 10.1016/j.cub.2012.06.005
Brotherton P, Haak W, Templeton J, Brandt G, Soubrier J, Jane Adler C, Richards SM, Der Sarkissian C, Ganslmeier R, Friederich S, Dresely V, van Oven M, Kenyon R, Van der Hoek MB, Korlach J, Luong K, Ho SYW, Quintana-Murci L, Behar DM, Meller H, Alt KW, Cooper A, Genographic Consortium. Neolithic mitochondrial haplogroup H genomes and the genetic origins of Europeans. Nature Communications. 2013; 4:1764. DOI: 10.1038/ncomms2656
Benazzi S, Slon V, Talamo S, Negrino F, Peresani M, Bailey SE, Sawyer S, Panetta D, Vicino G, Starnini E, Mannino MA, Salvadori PA, Meyer M, Pääbo S, Hublin JJ. Archaeology. The makers of the Protoaurignacian and implications for Neandertal extinction. Science. 2015; 348:793-796. DOI: 10.1126/science.aaa2773
Lazaridis I, Patterson N, Mittnik A, Renaud G, Mallick S, Kirsanow K, Sudmant PH, Schraiber JG, Castellano S, Lipson M, Berger B, Economou C, Bollongino R, Fu Q, Bos KI, Nordenfelt S, Li H, de Filippo C, Prüfer K, Sawyer S, Posth C, Haak W, Hallgren F, Fornander E, Rohland N, Delsate D, Francken M, Guinet JM, Wahl J, Ayodo G, Babiker HA, Bailliet G, Balanovska E, Balanovsky O, Barrantes R, Bedoya G, Ben-Ami H, Bene J, Berrada F, Bravi CM, Brisighelli F, Busby GB, Cali F, Churnosov M, Cole DE, Corach D, Damba L, van Driem G, Dryomov S, Dugoujon JM, Fedorova SA, Gallego Romero I, Gubina M, Hammer M, Henn BM, Hervig T, Hodoglugil U, Jha AR, Karachanak-Yankova S, Khusainova R, Khusnutdinova E, Kittles R, Kivisild T, Klitz W, Kučinskas V, Kushniarevich A, Laredj L, Litvinov S, Loukidis T, Mahley RW, Melegh B, Metspalu E, Molina J, Mountain J, Näkkäläjärvi K, Nesheva D, Nyambo T, Osipova L, Parik J, Platonov F, Posukh O, Romano V, Rothhammer F, Rudan I, Ruizbakiev R, Sahakyan H, Sajantila A, Salas A, Starikovskaya EB, Tarekegn A, Toncheva D, Turdikulova S, Uktveryte I, Utevska O, Vasquez R, Villena M, Voevoda M, Winkler CA, Yepiskoposyan L, Zalloua P, Zemunik T, Cooper A, Capelli C, Thomas MG, Ruiz-Linares A, Tishkoff SA, Singh L, Thangaraj K, Villems R, Comas D, Sukernik R, Metspalu M, Meyer M, Eichler EE, Burger J, Slatkin M, Pääbo S, Kelso J, Reich D, Krause J. Ancient human genomes suggest three ancestral populations for present-day Europeans. Nature. 2014; 513:409-413. DOI: 10.1038/nature13673
Sikora M, Carpenter ML, Moreno-Estrada A, Henn BM, Underhill PA, Sánchez-Quinto F, Zara I, Pitzalis M, Sidore C, Busonero F, Maschio A, Angius A, Jones C, Mendoza-Revilla J, Nekhrizov G, Dimitrova D, Theodossiev N, Harkins TT, Keller A, Maixner F, Zink A, Abecasis G, Sanna S, Cucca F, Bustamante CD. Population genomic analysis of ancient and modern genomesyields new insights into the genetic ancestry of the Tyrolean Iceman and the genetic structure of Europe. PLoS Genetics. 2014; 10:e1004353. DOI: 10.1371/journal.pgen.1004353
Sikora M, Seguin-Orlando A, Sousa VC, Albrechtsen A, Korneliussen T, Ko A, Rasmussen S, Dupanloup I, Nigst PR, Bosch MD, Renaud G, Allentoft ME, Margaryan A, Vasilyev SV, Veselovskaya EV, Borutskaya SB, Deviese T, Comeskey D, Higham T, Manica A, Foley R, Meltzer DJ, Nielsen R, Excoffier L, Mirazon Lahr M, Orlando L, Willerslev E. Ancient genomes show social and reproductive behavior of early Upper Paleolithic foragers. Science. 2017; 358:659-662. DOI: 10.1126/science.aao1807
Günther T, Valdiosera C, Malmström H, Ureña I, Rodriguez-Varela R, Sverrisdóttir ÓO, Daskalaki EA, Skoglund P, Naidoo T, Svensson EM, Bermúdez de Castro JM, Carbonell E, Dunn M, Storå J, Iriarte E, Arsuaga JL, Carretero JM, Götherström A, Jakobsson M. Ancient genomes link early farmers from Atapuerca in Spain to modern-day Basques. Proceedings of the National Academy of Sciences of the United States of America. 2015; 112:11917-11922. DOI: 10.1073/pnas.1509851112
Gómez-Sánchez D, Olalde I, Pierini F, Matas-Lalueza L, Gigli E, Lari M, Civit S, Lozano M, Vergès JM, Caramelli D, Ramírez O, Lalueza-Fox C. Mitochondrial DNA from El Mirador cave (Atapuerca, Spain) reveals the heterogeneity of Chalcolithic populations. PLoS One. 2014; 9:e105105. DOI: 10.1371/journal.pone.0105105
de Blas Cortina MA. Producción e intercambio de metal: la singularidad de las minas de cobre prehistóricas del Aramo y El Milagro (Asturias). Studia Archaeologica. 1998; 88:71-103
Ginther C, Issel-Tarver L, King MC. Identifying individuals by sequencing mitochondrial DNA from teeth. Nature Genetics. 1992; 2:135-138. DOI: 10.1038/ng1092-135
Izagirre N, de-la-Rua C. An mtDNA analysis in ancient Basque populations: Implications for haplogroup V as a marker for a major Paleolithic expansion from southwestern Europe. American Journal of Human Genetics. 1999; 65:199-207. DOI: 10.1086/302442
Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turnbull DM, Howell N. Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nature Genetics. 1999; 23:147. DOI: 10.1038/13779
Alonso A, Martin P, Albarran C, García P, García O, de Simón LF, García-Hirschfel J, Sancho M, de-la-Rúa C, Fernández-Piqueras J. Multiplex-PCR of short amplicons for mtDNA sequencing from ancient DNA. International Congress Series. 2003; 1239:585-588
Torroni A, Huoponen K, Francalacci P, Petrozzi M, Morelli L, Scozzari R, Obinu D, Savontaus ML. Classification of European mtDNAs from an analysis of three European populations. Genetics. 1996; 144:1835-1850
Pääbo S, Poinar H, Serre D, Jaenicke-Despres V, Hebler J, Rohland N, Kuch M, Krause J, Vigilant L, Hofreiter M. Genetic analyses from ancient DNA. Annual Reviews in Genetics. 2004; 38:645-679. DOI: 10.1146/annurev.genet.37.110801.143214
Gilbert MT, Willerslev E. Authenticity in ancient DNA studies. Medicina Nei Scoli. 2006; 18:701-723
Handt O, Krings M, Ward RH, Pääbo S. The retrieval of ancient human DNA sequences. American Journal of Human Genetics. 1996; 59:368-376
Malmström H, Gilbert MT, Thomas MG, Brandstrom M, Stora J, Molnar P, Andersen PK, Bendixen C, Holmlund G, Gotherstrom A, Willerslev E. Ancient DNA reveals lack of continuity between Neolithic hunter-gatherers and contemporary Scandinavians. Current Biology. 2009; 19:1758-1762. DOI: 10.1016/j.cub.2009.09.017
Bramanti B, Thomas MG, Haak W, Unterlaender M, Jores P, Tambets K, Antanaitis-Jacobs I, Haidle MN, Jankauskas R, Kind CJ, Lueth F, Terberger T, Hiller J, Matsumura S, Forster P, Burger J. Genetic discontinuity between local hunter-gatherers and Central Europe's first farmers. Science. 2009; 326:137-140. DOI: 0.1126/science.1176869
Haak W, Lazaridis I, Patterson N, Rohland N, Mallick S, Llamas B, Brandt G, Nordenfelt S, Harney E, Stewardson K, Fu Q, Mittnik A, Bánffy E, Economou C, Francken M, Friederich S, Pena RG, Hallgren F, Khartanovich V, Khokhlov A, Kunst M, Kuznetsov P, Meller H, Mochalov O, Moiseyev V, Nicklisch N, Pichler SL, Risch R, Rojo Guerra MA, Roth C, Szécsényi-Nagy A, Wahl J, Meyer M, Krause J, Brown D, Anthony D, Cooper A, Alt KW, Reich D. Massive migration from the steppe was a source for indo-European languages in Europe. Nature. 2015; 522:207-211. DOI: 10.1038/nature14317
Keyser C, Bouakaze C, Crubézy E, Nikolaev VG, Montagnon D, Reis T, Ludes B. Ancient DNA provides new insights into the history of south Siberian kurgan people. Human Genetics. 2009; 126:395-410. DOI: 10.1007/s00439-009-0683-0
Lalueza-Fox C, Sampietro ML, Gilbert MT, Castri L, Facchini F, Pettener D, Bertranpetit J. Unravelling migrations in the steppe: Mitochondrial DNA sequences from ancient central Asians. Proceedings of the Biological Sciences. 2004; 271:941-947. DOI: 10.1098/rspb.2004.2698
Pardiñas AF, Roca A, García-Vazquez E, López B. Assessing the genetic influence of ancient sociopolitical structure: micro-differentiation patterns in the population of Asturias (Northern Spain). PLoS One. 2012; 7:e50206. DOI: 10.1371/journal.pone.0050206
Pardiñas AF, Martínez JL, Roca A, García-Vazquez E, López B. Over the sands and far away: interpreting an Iberian mitochondrial lineage with ancient Western African origins. American Journal Human Biology. 2014; 26:777-783. DOI: 10.1002/ajhb.22601
Schneider S, Excoffier L. Estimation of past demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: Application to human mitochondrial DNA. Genetics. 1999; 152:1079-1089
Meyer S, Weiss G, von Haeseler A. Pattern of nucleotide substitution and rate heterogeneity in the hypervariable regions I and II of human mtDNA. Genetics. 1999; 152:1103-1110
Allard MW, Miller K, Wilson M, Monson K, Budowle B. Characterization of the Caucasian haplogroups present in the SWGDAM forensic mtDNA dataset for 1771 human control region sequences. Scientific Working Group on DNA Analysis Methods. Journal of Forensic Sciences. 2002; 47:1215-1223
Santos C, Montiel R, Sierra B, Bettencourt C, Fernandez E, Alvarez L, Lima M, Abade A, Aluja MP. Understanding differences between phylogenetic and pedigree-derived mtDNA mutation rate: A model using families from the Azores Islands (Portugal). Molecular Biology and Evolution. 2005; 22:1490-1505. DOI: 10.1093/molbev/msi141
van Oven M, Kayser M. Updated comprehensive phylogenetic tree of global human mitochondrial DNA variation. Human Mutation. 2009; 30:386-394. DOI: 10.1002/humu.20921
Hervella M, Fernandez-Crespo T, Schulting RJ, Izagirre N, de Blas MA, de-la-Rua C. Anthropological, genetic and isotopic analysis of human remains recovered from the copper mine of the Aramo mountains (Asturias). Actas del VI Congreso del Neolítico en la Península Ibérica; Granada; 2018, in press
de Blas Cortina MA. La mina como ámbito infraterreno y el cadáver como ofrenda ritual. A propósito de los esqueletos humanos hallados en las explotaciones cupríferas del Aramo. In: Fernández Manzano J, Herrán Martínez JI, editors. Mineros y fundidores en el inicio de la Edad de los Metales. El Midí francés y el Norte de la Península Ibérica. León: Caja España y Fundación Las Médulas; 2003. pp. 32-48