Genetic continuity in the Franco-Cantabrian region: new clues from autochthonous mitogenomes

Forensically relevant phylogeographic evaluation of mitogenome variation in the Basque Country

The Basque Country has been the focus of population (genetic) and evolutionary studies for decades, as it represents an interesting evolutionary feature: it is the only European country where a non-Indo-European language is still spoken today and, for which there are no known living or extinct relatives. Early studies that were based on anatomical and serological methods, along with subsequent molecular genetic investigations, contain controversial interpretations of their data. Additionally, the analysis of mitochondrial DNA, which is maternally inherited and thus suitable for the examination of the maternal phylogeny of the population, was the focus of some studies. Early mtDNA studies were however restricted to the information provided by the control region or its hypervariable segments only. These are known to harbour little phylogenetic information, particularly for haplogroup H that is dominant in Westeurasian populations including the Basques. Later studies analysed complete mitogenome sequences. Their information content is however limited, either because the number of samples was low, or because these studies only considered particular haplogroups. In this study we present the full mitogenome sequences of 178 autochthonous Basque individuals that were carefully selected based on their familial descent and discuss the observed phylogenetic signals in the light of earlier published findings. We confirm the presence of Basque-specific mtDNA lineages and extend the knowledge of these lineages by providing data on their distribution in comparison to other Basque and non-Basque populations. This dataset improves our understanding of the Basque mtDNA phylogeny and serves as a high-quality dataset that is provided via EMPOP for forensic genetic purposes.

Rare human mitochondrial HV lineages spread from the Near East and Caucasus during post-LGM and Neolithic expansions

Of particular significance to human population history in Eurasia are the migratory events that connected the Near East to Europe after the Last Glacial Maximum (LGM). Utilizing 315 HV*(xH,V) mitogenomes, including 27 contemporary lineages first reported here, we found the genetic signatures for distinctive movements out of the Near East and South Caucasus both westward into Europe and eastward into South Asia. The parallel phylogeographies of rare, yet widely distributed HV*(xH,V) subclades reveal a connection between the Italian Peninsula and South Caucasus, resulting from at least two (post-LGM, Neolithic) waves of migration. Many of these subclades originated in a population ancestral to contemporary Armenians and Assyrians. One such subclade, HV1b-152, supports a postexilic, northern Mesopotamian origin for the Ashkenazi HV1b2 lineages. In agreement with ancient DNA findings, our phylogenetic analysis of HV12 and HV14, the two exclusively Asian subclades of HV*(xH,V), point to the migration of lineages originating in Iran to South Asia before and during the Neolithic period. With HV12 being one of the oldest HV subclades, our results support an origin of HV haplogroup in the region defined by Western Iran, Mesopotamia, and the South Caucasus, where the highest prevalence of HV has been found.

The peopling of South America and the trans-Andean gene flow of the first settlers

Genetic and archaeological data indicate that the initial Paleoindian settlers of South America followed two entry routes separated by the Andes and the Amazon rainforest. The interactions between these paths and their impact on the peopling of South America remain unclear. Analysis of genetic variation in the Peruvian Andes and regions located south of the Amazon River might provide clues on this issue. We analyzed mitochondrial DNA variation at different Andean locations and >360,000 autosomal SNPs from 28 Native American ethnic groups to evaluate different trans-Andean demographic scenarios. Our data reveal that the Peruvian Altiplano was an important enclave for early Paleoindian expansions and point to a genetic continuity in the Andes until recent times, which was only marginally affected by gene flow from the Amazonian lowlands. Genomic variation shows a good fit with the archaeological evidence, indicating that the genetic interactions between the descendants of the settlers that followed the Pacific and Atlantic routes were extremely limited.

Ancient mitochondrial lineages support the prehistoric maternal root of Basques in Northern Iberian Peninsula

The Basque population inhabits the Franco-Cantabrian region in southwest Europe where Palaeolithic human groups took refuge during the Last Glacial Maximum. Basques have been an isolated population, largely considered as one of the most ancient European populations and it is possible that they maintained some pre-Neolithic genetic characteristics. This work shows the results of mitochondrial DNA analysis of seven ancient human remains from the Cave of Santimamiñe in the Basque Country dated from Mesolithic to the Late Roman period. In addition, we compared these data with those obtained from a modern sample of Basque population, 158 individuals that nowadays inhabits next to the cave. The results support the hypothesis that Iberians might have been less affected by the Neolithic mitochondrial lineages carried from the Near East than populations of Central Europe and revealed the unexpected presence of prehistoric maternal lineages such as U5a2a and U3a in the Basque region. Comparison between ancient and current population samples upholds the hypothesis of continuity of the maternal lineages in the area of the Franco-Cantabrian region.

Phylogenetic and population-based approaches to mitogenome variation do not support association with male infertility

Infertility has a complex multifactorial etiology and a high prevalence worldwide. Several studies have pointed to variation in the mitochondrial DNA (mtDNA) molecule as a factor responsible for the different disease phenotypes related to infertility. We analyzed 53 mitogenomes of infertile males from Galicia (northwest Spain), and these haplotypes were meta-analyzed phylogenetically with 43 previously reported from Portugal. Taking advantage of the large amount of information available, we additionally carried out association tests between patient mtDNA single-nucleotide polymorphisms (mtSNPs) and haplogroups against Iberian matched controls retrieved from The 1000 Genomes Project and the literature. Phylogenetic and association analyses did not reveal evidence of association between mtSNPs/haplogroups and infertility. Ratios and patterns in patients of nonsynonymous/synonymous changes, and variation at homoplasmic, heteroplasmic and private variants, fall within expected values for healthy individuals. Moreover, the haplogroup background of patients was variable and fits well with patterns typically observed in healthy western Europeans. We did not find evidence of association of mtSNPs or haplogroups pointing to a role for mtDNA in male infertility. A thorough review of the literature on mtDNA variation and infertility revealed contradictory findings and methodological and theoretical problems that overall undermine previous positive findings.

Meta-Analysis of Mitochondrial DNA Variation in the Iberian Peninsula

The Iberian Peninsula has been the focus of attention of numerous studies dealing with mitochondrial DNA (mtDNA) variation, most of them targeting the control region segment. In the present study we sequenced the control region of 3,024 Spanish individuals from areas where available data were still limited. We also compiled mtDNA haplotypes from the literature involving 4,588 sequences and 28 population groups or small regions. We meta-analyzed all these data in order to shed further light on patterns of geographic variation, taking advantage of the large sample size and geographic coverage, in contrast with the atomized sampling strategy of previous work. The results indicate that the main mtDNA haplogroups show primarily clinal geographic patterns across the Iberian geography, roughly along a North-South axis. Haplogroup HV0 (where haplogroup U is nested) is more prevalent in the Franco Cantabrian region, in good agreement with previous findings that identified this area as a climate refuge during the Last Glacial Maximum (LGM), prior to a subsequent demographic re-expansion towards Central Europe and the Mediterranean. Typical sub-Saharan and North African lineages are slightly more prevalent in South Iberia, although at low frequencies; this pattern has been shaped mainly by the transatlantic slave trade and the Arab invasion of the Iberian Peninsula. The results also indicate that summary statistics that aim to measure molecular variation, or AMOVA, have limited sensitivity to detect population substructure, in contrast to patterns revealed by phylogeographic analysis. Overall, the results suggest that mtDNA variation in Iberia is substantially stratified. These patterns might be relevant in biomedical studies given that stratification is a common cause of false positives in case-control mtDNA association studies, and should be also considered when weighting the DNA evidence in forensic casework, which is strongly dependent on haplotype frequencies.

Different Evolutionary History for Basque Diaspora Populations in USA and Argentina Unveiled by Mitochondrial DNA Analysis

The Basque Diaspora in Western USA and Argentina represents two populations which have maintained strong Basque cultural and social roots in a completely different geographic context. Hence, they provide an exceptional opportunity to study the maternal genetic legacy from the ancestral Basque population and assess the degree of genetic introgression from the host populations in two of the largest Basque communities outside the Basque Country. For this purpose, we analyzed the complete mitochondrial DNA control region of Basque descendants living in Western USA (n = 175) and in Argentina (n = 194). The Diaspora populations studied here displayed a genetic diversity in their European maternal input which was similar to that of the Basque source populations, indicating that not important founder effects would have occurred. Actually, the genetic legacy of the Basque population still prevailed in their present-day maternal pools, by means of a haplogroup distribution similar to the source population characterized by the presence of autochthonous Basque lineages, such as U5b1f1a and J1c5c1. However, introgression of non-Basque lineages, mostly Native American, has been observed in the Diaspora populations, particularly in Argentina, where the quick assimilation of the newcomers would have favored a wider admixture with host populations. In contrast, a longer isolation of the Diaspora groups in USA, because of language and cultural differences, would have limited the introgression of local lineages. This study reveals important differences in the maternal evolutionary histories of these Basque Diaspora populations, which have to be taken into consideration in forensic and medical genetic studies.

Fine Dissection of Human Mitochondrial DNA Haplogroup HV Lineages Reveals Paleolithic Signatures from European Glacial Refugia

Genetic signatures from the Paleolithic inhabitants of Eurasia can be traced from the early divergent mitochondrial DNA lineages still present in contemporary human populations. Previous studies already suggested a pre-Neolithic diffusion of mitochondrial haplogroup HV*(xH,V) lineages, a relatively rare class of mtDNA types that includes parallel branches mainly distributed across Europe and West Asia with a certain degree of structure. Up till now, variation within haplogroup HV was addressed mainly by analyzing sequence data from the mtDNA control region, except for specific sub-branches, such as HV4 or the widely distributed haplogroups H and V. In this study, we present a revised HV topology based on full mtDNA genome data, and we include a comprehensive dataset consisting of 316 complete mtDNA sequences including 60 new samples from the Italian peninsula, a previously underrepresented geographic area. We highlight points of instability in the particular topology of this haplogroup, reconstructed with BEAST-generated trees and networks. We also confirm a major lineage expansion that probably followed the Late Glacial Maximum and preceded Neolithic population movements. We finally observe that Italy harbors a reservoir of mtDNA diversity, with deep-rooting HV lineages often related to sequences present in the Caucasus and the Middle East. The resulting hypothesis of a glacial refugium in Southern Italy has implications for the understanding of late Paleolithic population movements and is discussed within the archaeological cultural shifts occurred over the entire continent.

No evidence of association between common European mitochondrial DNA variants in Alzheimer, Parkinson, and migraine in the Spanish population

Certain mitochondrial DNA (mtDNA) variants and haplogroups have been found to be associated with neurological disorders. Several studies have suggested that mtDNA variation could have an etiologic role in these disorders by affecting the ATP production on high-energy demanding organs, such as the brain. We have analyzed 15 mtDNA SNPs (mtSNPs) in five cohorts of cases presenting Alzheimer disease (AD), Parkinson disease (PD), and migraine, and in controls, to evaluate the role mtDNA variation in disease risk. Association tests were undertaken both for mtSNPs and mitochondrial haplogroups. No significant association was detected for any mtSNP or haplogroup in AD and PD cohorts. Two mtSNPs were associated with one migraine cohort after correcting for multiple tests, namely, T4216C and G13708A and haplogroup J (FDR q-value = 0.02; Santiago's cohort). However, this association was not confirmed in a second replication migraine series. A review of the literature reveals the existence of inconsistent findings and methodological shortcomings affecting a large proportion of mtDNA association studies on AD, PD, and migraine. A detailed inspection of the literature highlights the need for performing more rigorous methodological and statistical standards in mtDNA genetic association studies aimed to avoid false positive results of association between mtDNA variants and neurological diseases.

Indian signatures in the westernmost edge of the European Romani diaspora: new insight from mitogenomes

In agreement with historical documentation, several genetic studies have revealed ancestral links between the European Romani and India. The entire mitochondrial DNA (mtDNA) of 27 Spanish Romani was sequenced in order to shed further light on the origins of this population. The data were analyzed together with a large published dataset (mainly hypervariable region I [HVS-I] haplotypes) of Romani (N=1,353) and non-Romani worldwide populations (N>150,000). Analysis of mitogenomes allowed the characterization of various Romani-specific clades. M5a1b1a1 is the most distinctive European Romani haplogroup; it is present in all Romani groups at variable frequencies (with only sporadic findings in non-Romani) and represents 18% of their mtDNA pool. Its phylogeographic features indicate that M5a1b1a1 originated 1.5 thousand years ago (kya; 95% CI: 1.3-1.8) in a proto-Romani population living in Northwest India. U3 represents the most characteristic Romani haplogroup of European/Near Eastern origin (12.4%); it appears at dissimilar frequencies across the continent (Iberia: ≈ 31%; Eastern/Central Europe: ≈ 13%). All U3 mitogenomes of our Iberian Romani sample fall within a new sub-clade, U3b1c, which can be dated to 0.5 kya (95% CI: 0.3-0.7); therefore, signaling a lower bound for the founder event that followed admixture in Europe/Near East. Other minor European/Near Eastern haplogroups (e.g. H24, H88a) were also assimilated into the Romani by introgression with neighboring populations during their diaspora into Europe; yet some show a differentiation from the phylogenetically closest non-Romani counterpart. The phylogeny of Romani mitogenomes shows clear signatures of low effective population sizes and founder effects. Overall, these results are in good agreement with historical documentation, suggesting that cultural identity and relative isolation have allowed the Romani to preserve a distinctive mtDNA heritage, with some features linking them unequivocally to their ancestral Indian homeland.

Mitogenomes from two uncommon haplogroups mark late glacial/postglacial expansions from the near east and neolithic dispersals within Europe

The current human mitochondrial (mtDNA) phylogeny does not equally represent all human populations but is biased in favour of representatives originally from north and central Europe. This especially affects the phylogeny of some uncommon West Eurasian haplogroups, including I and W, whose southern European and Near Eastern components are very poorly represented, suggesting that extensive hidden phylogenetic substructure remains to be uncovered. This study expanded and re-analysed the available datasets of I and W complete mtDNA genomes, reaching a comprehensive 419 mitogenomes, and searched for precise correlations between the ages and geographical distributions of their numerous newly identified subclades with events of human dispersal which contributed to the genetic formation of modern Europeans. Our results showed that haplogroups I (within N1a1b) and W originated in the Near East during the Last Glacial Maximum or pre-warming period (the period of gradual warming between the end of the LGM, ∼19 ky ago, and the beginning of the first main warming phase, ∼15 ky ago) and, like the much more common haplogroups J and T, may have been involved in Late Glacial expansions starting from the Near East. Thus our data contribute to a better definition of the Late and postglacial re-peopling of Europe, providing further evidence for the scenario that major population expansions started after the Last Glacial Maximum but before Neolithic times, but also evidencing traces of diffusion events in several I and W subclades dating to the European Neolithic and restricted to Europe.

The expanded mtDNA phylogeny of the Franco-Cantabrian region upholds the pre-neolithic genetic substrate of Basques

The European genetic landscape has been shaped by several human migrations occurred since Paleolithic times. The accumulation of archaeological records and the concordance of different lines of genetic evidence during the last two decades have triggered an interesting debate concerning the role of ancient settlers from the Franco-Cantabrian region in the postglacial resettlement of Europe. Among the Franco-Cantabrian populations, Basques are regarded as one of the oldest and more intriguing human groups of Europe. Recent data on complete mitochondrial DNA genomes focused on macrohaplogroup R0 revealed that Basques harbor some autochthonous lineages, suggesting a genetic continuity since pre-Neolithic times. However, excluding haplogroup H, the most representative lineage of macrohaplogroup R0, the majority of maternal lineages of this area remains virtually unexplored, so that further refinement of the mtDNA phylogeny based on analyses at the highest level of resolution is crucial for a better understanding of the European prehistory. We thus explored the maternal ancestry of 548 autochthonous individuals from various Franco-Cantabrian populations and sequenced 76 mitogenomes of the most representative lineages. Interestingly, we identified three mtDNA haplogroups, U5b1f, J1c5c1 and V22, that proved to be representative of Franco-Cantabria, notably of the Basque population. The seclusion and diversity of these female genetic lineages support a local origin in the Franco-Cantabrian area during the Mesolithic of southwestern Europe, ∼10,000 years before present (YBP), with signals of expansions at ∼3,500 YBP. These findings provide robust evidence of a partial genetic continuity between contemporary autochthonous populations from the Franco-Cantabrian region, specifically the Basques, and Paleolithic/Mesolithic hunter-gatherer groups. Furthermore, our results raise the current proportion (≈15%) of the Franco-Cantabrian maternal gene pool with a putative pre-Neolithic origin to ≈35%, further supporting the notion of a predominant Paleolithic genetic substrate in extant 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.

Assessing the genetic influence of ancient sociopolitical structure: micro-differentiation patterns in the population of Asturias (Northern Spain)

The human populations of the Iberian Peninsula are the varied result of a complex mixture of cultures throughout history, and are separated by clear social, cultural, linguistic or geographic barriers. The stronger genetic differences between closely related populations occur in the northern third of Spain, a phenomenon commonly known as "micro-differentiation". It has been argued and discussed how this form of genetic structuring can be related to both the rugged landscape and the ancient societies of Northern Iberia, but this is difficult to test in most regions due to the intense human mobility of previous centuries. Nevertheless, the Spanish autonomous community of Asturias shows a complex history which hints of a certain isolation of its population. This, joined together with a difficult terrain full of deep valleys and steep mountains, makes it suitable for performing a study of genetic structure, based on mitochondrial DNA and Y-Chromosome markers. Our analyses do not only show that there are micro-differentiation patterns inside the Asturian territory, but that these patterns are strikingly similar between both uniparental markers. The inference of barriers to gene flow also indicates that Asturian populations from the coastal north and the mountainous south seem to be relatively isolated from the rest of the territory. These findings are discussed in light of historic and geographic data and, coupled with previous evidence, show that the origin of the current genetic patterning might indeed lie in Roman and Pre-Roman sociopolitical divisions.

Multiple local and recent founder effects of TGM1 in Spanish families

BACKGROUND

Mutations in the TGM1 gene encoding transglutaminase 1 are a major cause of autosomal recessive congenital ichthyosis. In the Galician (NW Spain) population, three mutations, c.2278C>T, c.1223_1227delACAC and c.984+1G>A, were observed at high frequency, representing ~46%, ~21% and ~13% of all TGM1 gene mutations, respectively. Moreover, these mutations were reported only once outside of Galicia, pointing to the existence of historical episodes of local severe genetic drift in this region.

METHODOLOGY/PRINCIPAL FINDINGS

In order to determine whether these mutations were inherited from a common ancestor in the Galician population, and to estimate the number of generations since their initial appearance, we carried out a haplotype-based analysis by way of genotyping 21 SNPs within and flanking the TGM1 gene and 10 flanking polymorphic microsatellite markers spanning a region of 12 Mb. Two linkage disequilibrium based methods were used to estimate the time to the most recent common ancestor (TMRCA), while a Bayesian-based procedure was used to estimate the age of the two mutations. Haplotype reconstruction from unphased genotypes of all members of the affected pedigrees indicated that all carriers for each of the two mutations harbored the same haplotypes, indicating common ancestry.

CONCLUSIONS/SIGNIFICANCE

In good agreement with the documentation record and the census, both mutations arose between 2,800-2,900 years ago (y.a.), but their TMRCA was in the range 600-1,290 y.a., pointing to the existence of historical bottlenecks in the region followed by population growth. This demographic scenario finds further support on a Bayesian Coalescent Analysis based on TGM1 haplotypes that allowed estimating the occurrence of a dramatic reduction of effective population size around 900-4,500 y.a. (95% highest posterior density) followed by exponential growth.