Delineating the dispersal of Y-chromosome sub-haplogroup O2a2b-P164 among Austronesian-speaking populations.

This article reports on an exploration of the Y-chromosome sub-haplogroup O2a2b-P164 in Austronesian-speaking populations. Moderate to high abundance of the P 164 mutation is seen in the West Pacific including the Amis of Formosa (36%) and the Filipinos of Mindanao (50%) as well as in the Kiritimati of Micronesia (70%), and Tonga and Samoa of West Polynesia (54% and 33%, respectively), and it drops to low frequencies in populations of East Polynesia. The communities of Polynesia and Micronesia exhibit considerable inter- and intra-population haplotype sharing suggesting extensive population affinity. The observed affinities, as well as the ages and diversity values within the P 164 sub-haplogroup among Austronesian-speaking populations signal an ancestral migration route and relationships that link the Amis of Taiwan with distant communities in West and East Polynesia, Micronesia, and the Maori of New Zealand. High resolution sequencing of the Austronesian Y chromosome indicate that the P 164 lineage originated about 19,000 ya and then split into three branches separating the Ami aborigines, Southeast Asian and Polynesian/Micronesian populations about 4700 ya, roughly coinciding with the initiation of the Austronesian diaspora. The Y-chromosomes of all the Polynesian and Micronesian population examined belong to the new FT 257096 haplogroup.

On the Y chromosome of Chennai, Tamil Nadu and the Indian subcontinent.

Several migratory waves from various origins along with cultural practices restricting marriages between people of different castes and tribes as well as continued endogamy have led to a complex and diverse society in the Indian subcontinent. Despite being widely represented in genetic studies, several interrogatives remain with regards to India's current genetic constituents and distributions, source populations and population relationships. To identify the forces that may have shaped Indian population's genetic relationships, we undertook a comprehensive comparative study of the Y-chromosomes across India utilizing Y-STR and Y-SNP chromosomal markers using the general population of Chennai in the state of Tamil Nadu as a point of reference. Our analyses identify differences in source populations for different regions within India, unique linguistic characteristics as well as demographic and cultural forces that may have shaped population structure.

Bidirectional dispersals during the peopling of the North American Arctic.

It is thought that Paleo-Inuit were the first people that settled the American Arctic about 5000 BP (before the present) from a migration that crossed Beringia from Northeast Asia. It is theorized that this group initially migrated to the North Slopes of Alaska and subsequently expanded eastward, eventually reaching Greenland. A second circumpolar dispersal of Neo-Inuit from the North Slopes associated with the Thule-Inuk culture has been postulated to have extended eastward around 800 BP, totally replacing the original Paleo-Inuit without admixing. Although generally accepted, this migration scenario is incompatible with previously reported indications of east to west gene flow across the American Arctic. Here we report on the Y-chromosome haplogroup and Y-STR diversity of the four circumpolar populations of the Tuva Republic (N = 24), Northeast Siberia (N = 9), Bethel, Alaska (N = 40), and Barrow, Alaska (N = 31). Four haplogroup lineages (Q-NWT01, Q-M3, Q-M346, and Q-M120) were detected, Q-NWT01 and Q-M3 being the most abundant at 11.11 and 66.67% in Northeast Siberia, 32.50 and 65.00% in Bethel, and 67.74 and 32.26% in Barrow, respectively. The same samples genotyped for Y-chromosome SNPs were typed for 17 Y-STYR loci using the AmpFlSTR Yfiler system. Age estimates and diversity values for the Q-NWT01 and Q-M3 mutations suggest extensive movement of male individuals along the entire longitudinal stretch of the American circumpolar region. Throughout the entire region, Q-M3 exhibits a west to east decreasing gradient in age and diversity while Q-NWT01 indicates the opposite with older TMRCA and higher diversity values running from east to west with the most recent estimates in Canada and Alaska. The high age and diversity values in Greenland are congruent with an origin of the Q-NWT01 mutation in the east of the circumpolar range about 2000-3000 ya. This scenario is incompatible with a complete biological replacement starting about 700 BP of Paleo-Inuit like the Dorset by the Thule-Inuit (Neo-Inuit), as is currently thought, and more parsimonious with gene flow carrying the NWT01 mutation from a pre-Thule population to the ancestors of the present-day Inuit.

Genetic characterization of populations in the Marquesas Archipelago in the context of the Austronesian expansion.

Our exploration of the genetic constitution of Nuku Hiva (n = 51), Hiva Oa (n = 28) and Tahuata (n = 8) of the Marquesas Archipelago based on the analyses of genome-wide autosomal markers as well as high-resolution genotyping of paternal and maternal lineages provides us with information on the origins and settlement of these islands at the fringe of the Austronesian expansion. One widespread theme that emerges from this study is the genetic uniformity and relative isolation exhibited by the Marquesas and Society populations. This genetic homogeneity within East Polynesia groups is reflected in their limited average heterozygosity, uniformity of constituents in the Structure analyses, reiteration of complete mtDNA sequences, marked separation from Asian and other Oceanic populations in the PC analyses, limited differentiation in the PCAs and large number of IBD segments in common. Both the f3 and the Outgroup f3 results provide indications of intra-East Polynesian gene flow that may have promoted the observed intra-East Polynesia genetic homogeneity while ALDER analyses indicate that East Polynesia experienced two gene flow episodes, one relatively recent from Europe that coincides roughly with the European incursion into the region and an early one that may represent the original settlement of the islands by Austronesians. Median Network analysis based on high-resolution Y-STR loci under C2a-M208 generates a star-like topology with East Polynesian groups (especially from the Society Archipelago) in central stem positions and individuals from the different populations radiating out one mutational step away while several Samoan and outlier individuals occupy peripheral positions. This arrangement of populations is congruent with dispersals of C2a-M208 Y chromosomes from East Polynesia as a migration hub signaling dispersals in various directions. The equivalent ages of the C2a-M208 lineage of the populations in the Network corroborate an east to west flow of the most abundant Polynesian Y chromosome.

Post-Austronesian migrational wave of West Polynesians to Micronesia.

This study examines Y-chromosome and mtDNA markers in the population of the island of Kiritimati in the context of geographically targeted reference populations from the Pacific. Kiritimati derives its population from the atoll islands of the Gilbert Archipelago and representsa geographicaltransitional region between Micronesia, Polynesia and Melanesia that likely played a critical role during theAustronesian expansion. The large presence(84.1%)of individuals withO-M175, O2a-M324 and O2a2b-P164 sub-haplogroups, 69.9% being O2a2b-P164, the Y-STR homogeneity within O2a2b-P164 and the very recent age of the sub-haplogroup(363-548 years ago)inKiritimati suggestthe arrival ofa genetically homogenous population to the Gilberteses followed by a population expassion.The close Y-STR haplotype affinities with profiles from the Samoa and Tonga Archipelagos point to an unprecedented massive post-Austronesian expansionexodus from West Polynesia.Contrasting the abundance of AustronesianO2a2b-P164 sub-haplogroup, the most abundantMelanesian/Papuansub-haplogroup,C-M130is present at a frequency of 13.5%. Thenetwork topology suggests that C-M130 arrived to theKiribati Archipelago from West Polynesia, specifically from West Samoa, Tonga and/or Tutuila subsequent to the Austronesian expansion about 832-1408 years ago. The haplotype affinities withinO2a2b-P164 argue for anoriginal source in Taiwan and its dispersal to West Polynesia and then to Southeast Micronesia. The present investigation provides an understanding of the genetic composition and complex migration history of an understudied region of the Pacific and provides evidence for recent dispersals towards Micronesia from West Polynesia subsequent to the initial Austronesian expansion.

Association of LCT -13910C>T polymorphism and hip fracture in a cohort of older adult population from Northern Spain.

Hip fracture is a common health problem very frequent in the older adult population and is associated with significant morbidity, mortality, and societal costs. There are several factors that increase the risk of suffering a hip fracture, however, the effect of genetic lactase non-persistence is not clear-cut yet. For this reason, we investigated if the LCT -13910C>T polymorphism is a potential risk factor for osteoporotic hip fractures in older adult people from the Northern Spain population. A total of 740 individuals were included in this study. Of them, 364 belonged to the group of patients whit osteoporotic hip fracture while the control group consisted of 376 individuals without hip fracture. The genotypes for the LCT -13910C>T polymorphism were analyzed by using polymerase chain reaction and high resolution melting. The prevalence of the CC genotype, which is related to lactase non-persistence, did not differ significantly in both groups. Likewise, no differences were observed between groups when they were compared with regard to the C or the T allele, or when they were analyzed considering gender. Additionally, our results were compared with those obtained in a control group of 207 nonagenarian individuals originally from Northern Spain and no differences were observed. In conclusion, no significant association was observed between the LCT -13910C>T polymorphism and the risk for suffering hip fracture in the older adult population of Northern Spain.

The Ami and Yami aborigines of Taiwan and their genetic relationship to East Asian and Pacific populations.

This article reports on the genetic characteristics of the Ami and Yami, two aboriginal populations of Taiwan. Y-SNP and mtDNA markers as well as autosomal SNPs were utilized to investigate the phylogenetic relationships to groups from MSEA (mainland Southeast Asia), ISEA (island Southeast Asia), and Oceania. Both the Ami and Yami have limited genetic diversity, with the Yami having even less diversity than the Ami. The partitioning of populations within the PCA plots based on autosomal SNPs, the profile constitution observed in the structure analyses demonstrating similar composition among specific populations, the average IBD (identical by descent) tract length gradients, the average total length of genome share among the populations, and the outgroup f3 results all indicate genetic affinities among populations that trace a geographical arc from Taiwan south into the Philippine Archipelago, Borneo, Indonesia, and Melanesia. Conversely, a more distant kinship between the Ami/Yami and MSEA based on all the markers examined, the total mtDNA sequences as well as the admixture f3 and f4 analyses argue against strong genetic contribution from MSEA to the Austronesian dispersal. The sharing of long IBD tracts, total genome length, and the large number of segments in common between the Ami/Yami and the Society Archipelago populations East Polynesia standout considering they are located about 10,700 km apart.

The Y chromosome of autochthonous Basque populations and the Bronze Age replacement.

Here we report on the Y haplogroup and Y-STR diversity of the three autochthonous Basque populations of Alava (n = 54), Guipuzcoa (n = 30) and Vizcaya (n = 61). The same samples genotyped for Y-chromosome SNPs were typed for 17 Y-STR loci (DYS19, DYS385a/b, DYS398I/II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635, Y-GATA H4) using the AmpFlSTR Yfiler system. Six major haplogroups (R, I, E, J, G, and DE) were detected, being R-S116 (P312) haplogroup the most abundant at 75.0% in Alava, 86.7% in Guipuzcoa and 87.3% in Vizcaya. Age estimates for the R-S116 mutation in the Basque Country are 3975 ± 303, 3680 ± 345 and 4553 ± 285 years for Alava, Guipuzcoa and Vizcaya, respectively. Pairwise Rst genetic distances demonstrated close Y-chromosome affinities among the three autochthonous Basque populations and between them and the male population of Ireland and Gascony. In a MDS plot, the population of Ireland segregates within the Basque cluster and closest to the population of Guipuzcoa, which plots closer to Ireland than to any of the other Basque populations. Overall, the results support the notion that during the Bronze Age a dispersal of individuals carrying the R-S116 mutation reached the Basque Country replacing the Paleolithic/Neolithic Y chromosome of the region.

Study of 17 X-STRs in Native American and Mestizo populations of Central America for forensic and population purposes.

In the present work, an extensive analysis of the X-chromosomal pool of Native American and Mestizo groups of Central America (Guatemala, El Salvador, Nicaragua, and Panama) has been carried out. Allele and haplotype frequency databases, as well as other forensic parameters for these populations, are presented. The admixture analysis supports the tri-hybrid composition in terms of ancestry in the Mestizo populations, with a predominant Native American contribution (54-69%), followed by European (19-28%) and African contributions (12-19%). Pairwise FST genetic distances highlight the genetic proximity between the northernmost Central American populations, especially among admixed populations. The unique and complex nature of this area, where populations from different origins intercrossed, as well as the informativity of X-STR data, highpoint the great interest of this genetic study. Furthermore, the X-chromosome databases for Central American populations here provided will be not only useful for forensic and population purposes not only in the target countries but also in the host countries.

Gene-Diet Interactions in Colorectal Cancer: Survey Design, Instruments, Participants and Descriptive Data of a Case-Control Study in the Basque Country.

Epidemiologic studies have revealed inconsistent evidence of gene-diet interaction in relation to colorectal cancer (CRC). The aim of this study was to analyze them in a sample of cases and controls from the population-based bowel cancer screening program of the Osakidetza/Basque Health Service. This study analyzed dietetic, genetic, demographic, socioeconomic factors and lifestyles. In the present manuscript, the survey design, sampling, instruments, measurements and related quality management were presented. Moreover, we analyze differences between cases and controls in some data, especially those related to diet. The participants were 308 cases and 308 age- and sex-matched subjects as controls. Cases were more likely than controls to have overweight/obesity (67.5% vs. 58.1%, p < 0.05), a lower intake of vitamin B2 (0.86 ± 0.23 vs. 0.92 ± 0.23 mg/1000 kcal, p < 0.01) and calcium:phosphorus ratio (0.62 ± 0.12 vs. 0.65 ± 0.13, p < 0.01). A higher proportion of cases than controls did not meet the Nutritional Objectives for saturated fatty acids (85.7% vs. 67.5%, p < 0.001) or cholesterol (35.4% vs. 25.0%, p < 0.01). In conclusion, the present study provides valuable data for analyzing the complexity of gene-diet interaction in relation to CRC. The results presented here suggest that overweight/obesity and a high intake of certain dietary components, especially saturated fatty acids and cholesterol, are more frequent in cases than in controls.

Evaluation of the VISAGE Basic Tool for Appearance and Ancestry Prediction Using PowerSeq Chemistry on the MiSeq FGx System.

The study of DNA to predict externally visible characteristics (EVCs) and the biogeographical ancestry (BGA) from unknown samples is gaining relevance in forensic genetics. Technical developments in Massively Parallel Sequencing (MPS) enable the simultaneous analysis of hundreds of DNA markers, which improves successful Forensic DNA Phenotyping (FDP). The EU-funded VISAGE (VISible Attributes through GEnomics) Consortium has developed various targeted MPS-based lab tools to apply FDP in routine forensic analyses. Here, we present an evaluation of the VISAGE Basic tool for appearance and ancestry prediction based on PowerSeq chemistry (Promega) on a MiSeq FGx System (Illumina). The panel consists of 153 single nucleotide polymorphisms (SNPs) that provide information about EVCs (41 SNPs for eye, hair and skin color from HIrisPlex-S) and continental BGA (115 SNPs; three overlap with the EVCs SNP set). The assay was evaluated for sensitivity, repeatability and genotyping concordance, as well as its performance with casework-type samples. This targeted MPS assay provided complete genotypes at all 153 SNPs down to 125 pg of input DNA and 99.67% correct genotypes at 50 pg. It was robust in terms of repeatability and concordance and provided useful results with casework-type samples. The results suggest that this MPS assay is a useful tool for basic appearance and ancestry prediction in forensic genetics for users interested in applying PowerSeq chemistry and MiSeq for this purpose.

Validation of a 52-mtSNP minisequencing panel for haplogroup classification of forensic DNA samples.

Mitochondrial DNA (mtDNA) is a useful tool in forensic investigation as it provides information about the matrilineal ancestry of individuals. In addition, mtDNA can be analyzed when the analysis of other nuclear markers is underperforming. Recently, we developed a minisequencing panel for the simultaneous analysis of 52 mtDNA SNPs to classify maternal lineages into the main haplogroups and their phylogeographic origin. In order to make this panel suitable for forensic genetics laboratories, a validation study has been performed in accordance with the Scientific Working Group on DNA Analysis Methods (SWGDAM) guidelines, including species specificity, reproducibility, sensitivity, and stability tests. The results demonstrate that the panel of 52 mtDNA SNPs is highly sensitive, since it enables to obtain complete genetic profiles of samples containing minimal amounts of DNA (1 pg). Furthermore, it provides sufficient genetic information to detect the matrilineal biogeographical origin of highly degraded samples, i.e., ancient dating skeletal remains, and samples with the presence of inhibitors, such as hematin and humic acid. In addition, this panel can detect mixtures in samples whose mtDNA haplogroups of contributors are different. Overall, the results of this study demonstrate the suitability of this minisequencing panel of 52 mtDNA SNPs to be used in forensic cases, with samples of low amount or degraded DNA.

Single nucleotide polymorphisms associated with susceptibility for development of colorectal cancer: Case-control study in a Basque population.

Given the significant population diversity in genetic variation, we aimed to investigate whether single nucleotide polymorphisms (SNPs) previously identified in studies of colorectal cancer (CRC) susceptibility were also relevant to the population of the Basque Country (North of Spain). We genotyped 230 CRC cases and 230 healthy controls for 48 previously reported CRC-susceptibility SNPs. Only the rs6687758 in DUPS10 exhibited a statistically significant association with CRC risk based on the crude analysis. The rs6687758 AG genotype conferred about 2.13-fold increased risk for CRC compared to the AA genotype. Moreover, we found significant associations in cases between smoking status, physical activity, and the rs6687758 SNP. The results of a Genetic Risk Score (GRS) showed that the risk alleles were more frequent in cases than controls and the score was associated with CRC in crude analysis. In conclusion, we have confirmed a CRC susceptibility locus and the existence of associations between modifiable factors and the rs6687758 SNP; moreover, the GRS was associated with CRC. However, further experimental validations are needed to establish the role of this SNP, the function of the gene identified, as well as the contribution of the interaction between environmental factors and this locusto the risk of CRC.

A unique common ancestor introduced P301L mutation in MAPT gene in frontotemporal dementia patients from Barcelona (Baix Llobregat, Spain).

The County of Baix Llobregat (Barcelona, Catalonia, Spain) presents a high prevalence of familial frontotemporal dementia (FTD) in the presence of P301L mutation in the MAPT gene. To evaluate a possible unique founder effect of P301L, and its age, the analysis of 20 single-nucleotide polymorphisms covering 50 kb and 12 single-nucleotide polymorphisms located along 30 Mb around the mutation was performed by developing 2 multiplex single-base extension reactions. In addition, families with affected and healthy individuals from France and Italy were analyzed. The FTD-affected individuals from Barcelona carried the same 50-kb haplotype linked to P301L mutation, suggesting a unique common ancestor, as opposed to French patients. Italian patients are also probably descendants of a unique ancestor, which would be different from that of Barcelona. Diversity of 30-Mb haplotypes found in Barcelona and the inference of the mutation age in these populations, among other reasons, suggest that prevalence of FTD linked to P301L MAPT mutation is the result of a locally originated mutation.

The Marquesans at the fringes of the Austronesian expansion.

In the present study, 87 unrelated individuals from the Marquesas Archipelago in French Polynesia were typed using mtDNA, Y-chromosome and autosomal (STRs) markers and compared to key target populations from Island South East Asia (ISEA), Taiwan, and West and East Polynesia to investigate their genetic relationships. The Marquesas, located at the eastern-most fringes of the Austronesian expansion, offer a unique opportunity to examine the effects of a protracted population expansion wave on population structure. We explore the contribution of Melanesian, Asian and European heritage to the Marquesan islands of Nuku-Hiva, Hiva-Oa and Tahuata. Overall, the Marquesas Islands are genetically homogeneous. In the Marquesan Archipelago all of the mtDNA haplogroups are of Austronesian origin belonging to the B4a1 subhaplogroup as the region marks the end of a west to east decreasing cline of Melanesian mtDNA starting with the West Polynesian population of Tonga. Genetic discrepancies are less pronounced between the Marquesan and Society islands, and among the Marquesan islands. Interestingly, a number of Melanesian, Polynesian and European Y-chromosome haplogroups exhibit very different distribution between the Marquesan islands of Nuku Hiva and Hiva Oa, likely resulting from drift, differential migration involving various source populations and/or unique trading routes.

Effective resolution of the Y chromosome sublineages of the Iberian haplogroup R1b-DF27 with forensic purposes.

Single-nucleotide polymorphisms (SNPs) found within the non-recombining region of the Y chromosome (NRY) represent a powerful tool in forensic genetics for inferring the paternal ancestry of a vestige and complement the determination of biogeographical origin in combination with other markers like AIMs. In the present study, we introduce a panel of 15 Y-SNPs for a fine-resolution subtyping of the haplogroup R1b-DF27, in a single minisequencing reaction. This is the first minisequencing panel that allows a fine subtyping of R1b-DF27, which displays high frequencies in Iberian and Iberian-influenced populations. This panel includes subhaplogroups of DF27 that display moderate geographical differentiation, of interest to link a sample with a specific location of the Iberian Peninsula or with Iberian ancestry. Conversely, part of the intricacy of a new minisequencing panel is to have all the included variants available to test the effectiveness of the analysis method. We have overcome the absence of the least common variants through site-directed mutagenesis. Overall, the results show that our panel is a robust and effective method for subtyping R1b-DF27 lineages from a minimal amount of DNA, and its high resolution enables to improve male lineage discrimination in Iberian and Southwest European descent individuals. The small length of the amplicons and its reproducibility makes this assay suitable for forensic and population genetics purposes.

Development of a new screening method to determine the main 52 mitochondrial haplogroups through a single minisequencing reaction.

This work presents the design, development and optimization of a screening method based on single-base extension sequencing to simultaneously analyze a panel of 52 mitochondrial SNPs. This enables to recognize the main mitochondrial haplogroups and to discriminate even between lineages from the same phylogenetic branch that diverged in different continents. The unavailability of individuals harboring infrequent variants was a limitation to optimize the panel. To overcome this, we have modified DNA by site-directed mutagenesis to create the unavailable allelic variants. This allowed us to verify the reliability of this panel and its usefulness to be applied in biomedicine, forensic and population genetic studies.

Differentially methylated CpG regions analyzed by PCR-high resolution melting for monozygotic twin pair discrimination.

Discrimination between monozygotic (MZ) twins is a forensic limitation when using conventional DNA profiling techniques for human identification. Recent works based on epigenetics seem to open a new way to solve this issue due to methylation status of MZ twins change during their lifetime. Methylation analysis through BeadChip platforms allows the study up to 850 K CpG sites revealing that numerous differential methylation regions exist between MZ twins. However, this methodology is difficult to implement in forensic laboratories. On the contrary, PCR-HRM (High Resolution Melting) technology is one of the easiest methods for analyzing DNA methylation and it has been capable to discriminate between MZ twins. The purpose of this study is to contribute with new differential methylation regions in MZ twins to those that have been previously studied through PCR-HRM. Here, we have selected 6 CpG regions located at the ITGA2B, ASPA, PDE4C, ZIC5, USP11 and NOP14 loci that have shown methylation status variation during lifetime. The study has been carried out from saliva-derived DNA of 18 MZ twin pairs. The most discriminating regions were those located at ITGA2B, ASPA and ZIC5 loci showing significant within-pair differences in 44.4% of the cases. Non evidences of relation between age and significant differences between MZ twins were found, although the 50% of MZ twin pairs were discrimnated in the oldest age range (59-66 years old). These results support the use of these regions to increase the number of epigenetics age-related markers available to discriminate between MZ twins in a pair by PCR-HRM in forensic laboratories.

Frontotemporal Dementia Caused by the P301L Mutation in the MAPT Gene: Clinicopathological Features of 13 Cases from the Same Geographical Origin in Barcelona, Spain.

BACKGROUND/AIMS: We identified and studied 13 patients carrying the P301L mutation in the MAPT gene from the same area (Baix Llobregat County) in Barcelona, Spain. METHODS: The demographic and clinical features were reviewed retrospectively. Detailed neuropathological characterization was obtained in 9 subjects. To investigate the origin of the P301L mutation in these families, 20 single nucleotide polymorphisms (SNPs) in the MAPT gene were analyzed. RESULTS: The mean age at disease onset was 51 years and the mean disease duration was 7 years. The most common initial symptoms were behavioral changes (54%), followed by language disturbances (31%) and memory loss (15%). 46% developed parkinsonism. Neuropathology showed an extensive neuronal and glial 4-repeat (4R) tauopathy with "mini-Pick"-like bodies in the dentate gyrus as the characteristic underlying pathology in all cases. In 1 subject, additional 4R globular glial inclusions were observed. All the mutation carriers showed the same haplotype for the SNPs analyzed, suggesting a common ancestor. CONCLUSION: These findings suggest a relative homogeneous clinicopathological phenotype in P301L MAPT mutation carriers in our series. This phenotype might help in the differential diagnosis from other tauopathies and be a morphological hint for genetic testing. The haplotype analysis results suggest a founder effect of the P301L mutation in this area.

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.