Contamination detection in sequencing studies using the mitochondrial phylogeny.

Within-species contamination is a major issue in sequencing studies, especially for mitochondrial studies. Contamination can be detected by analyzing the nuclear genome or by inspecting polymorphic sites in the mitochondrial genome (mtDNA). Existing methods using the nuclear genome are computationally expensive, and no appropriate tool for detecting sample contamination in large-scale mtDNA data sets is available. Here we present haplocheck, a tool that requires only the mtDNA to detect contamination in both targeted mitochondrial and whole-genome sequencing studies. Our in silico simulations and amplicon mixture experiments indicate that haplocheck detects mtDNA contamination accurately and is independent of the phylogenetic distance within a sample mixture. By applying haplocheck to The 1000 Genomes Project Consortium data, we further evaluate the application of haplocheck as a fast proxy tool for nDNA-based contamination detection using the mtDNA and identify the mitochondrial copy number within a mixture as a critical component for the overall accuracy. The haplocheck tool is available both as a command-line tool and as a cloud web service producing interactive reports that facilitates the navigation through the phylogeny of contaminated samples.

Analyzing Low-Level mtDNA Heteroplasmy-Pitfalls and Challenges from Bench to Benchmarking.

Massive parallel sequencing technologies are promising a highly sensitive detection of low-level mutations, especially in mitochondrial DNA (mtDNA) studies. However, processes from DNA extraction and library construction to bioinformatic analysis include several varying tasks. Further, there is no validated recommendation for the comprehensive procedure. In this study, we examined potential pitfalls on the sequencing results based on two-person mtDNA mixtures. Therefore, we compared three DNA polymerases, six different variant callers in five mixtures between 50% and 0.5% variant allele frequencies generated with two different amplification protocols. In total, 48 samples were sequenced on Illumina MiSeq. Low-level variant calling at the 1% variant level and below was performed by comparing trimming and PCR duplicate removal as well as six different variant callers. The results indicate that sensitivity, specificity, and precision highly depend on the investigated polymerase but also vary based on the analysis tools. Our data highlight the advantage of prior standardization and validation of the individual laboratory setup with a DNA mixture model. Finally, we provide an artificial heteroplasmy benchmark dataset that can help improve somatic variant callers or pipelines, which may be of great interest for research related to cancer and aging.

Mitochondrial DNA copy number is associated with mortality and infections in a large cohort of patients with chronic kidney disease.

Damage of mitochondrial DNA (mtDNA) with reduction in copy number has been proposed as a biomarker for mitochondrial dysfunction and oxidative stress. Chronic kidney disease (CKD) is associated with increased mortality and risk of cardiovascular disease, but the underlying mechanisms remain incompletely understood. Here we investigated the prognostic role of mtDNA copy number for cause-specific mortality in 4812 patients from the German Chronic Kidney Disease study, an ongoing prospective observational national cohort study of patients with CKD stage G3 and A1-3 or G1-2 with overt proteinuria (A3) at enrollment. MtDNA was quantified in whole blood using a plasmid-normalized PCR-based assay. At baseline, 1235 patients had prevalent cardiovascular disease. These patients had a significantly lower mtDNA copy number than patients without cardiovascular disease (fully-adjusted model: odds ratio 1.03, 95% confidence interval [CI] 1.01-1.05 per 10 mtDNA copies decrease). After four years of follow-up, we observed a significant inverse association between mtDNA copy number and all-cause mortality, adjusted for kidney function and cardiovascular disease risk factors (hazard ratio 1.37, 95% CI 1.09-1.73 for quartile 1 compared to quartiles 2-4). When grouped by causes of death, estimates pointed in the same direction for all causes but in a fully-adjusted model decreased copy numbers were significantly lower only in infection-related death (hazard ratio 1.82, 95% CI 1.08-3.08). A similar association was observed for hospitalizations due to infections in 644 patients (hazard ratio 1.19, 95% CI 1.00-1.42 in the fully-adjusted model). Thus, our data support a role of mitochondrial dysfunction in increased cardiovascular disease and mortality risks as well as susceptibility to infections in patients with CKD.

Rapid coastal spread of First Americans: novel insights from South America's Southern Cone mitochondrial genomes.

It is now widely agreed that the Native American founders originated from a Beringian source population ~15-18 thousand years ago (kya) and rapidly populated all of the New World, probably mainly following the Pacific coastal route. However, details about the migration into the Americas and the routes pursued on the continent still remain unresolved, despite numerous genetic, archaeological, and linguistic investigations. To examine the pioneering peopling phase of the South American continent, we screened literature and mtDNA databases and identified two novel mitochondrial DNA (mtDNA) clades, here named D1g and D1j, within the pan-American haplogroup D1. They both show overall rare occurrences but local high frequencies, and are essentially restricted to populations from the Southern Cone of South America (Chile and Argentina). We selected and completely sequenced 43 D1g and D1j mtDNA genomes applying highest quality standards. Molecular and phylogeographic analyses revealed extensive variation within each of the two clades and possibly distinct dispersal patterns. Their age estimates agree with the dating of the earliest archaeological sites in South America and indicate that the Paleo-Indian spread along the entire longitude of the American double continent might have taken even <2000 yr. This study confirms that major sampling and sequencing efforts are mandatory for uncovering all of the most basal variation in the Native American mtDNA haplogroups and for clarification of Paleo-Indian migrations, by targeting, if possible, both the general mixed population of national states and autochthonous Native American groups, especially in South America.

Forensic and phylogeographic characterisation of mtDNA lineages from Somalia.

The African mitochondrial (mt) phylogeny is coarsely resolved but the majority of population data generated so far is limited to the analysis of the first hypervariable segment (HVS-1) of the control region (CR). Therefore, this study aimed on the investigation of the entire CR of 190 unrelated Somali individuals to enrich the severely underrepresented African mtDNA pool. The majority (60.5 %) of the haplotypes were of sub-Saharan origin with L0a1d, L2a1h and L3f being the most frequently observed haplogroups. This is in sharp contrast to previous data reported from the Y-chromosome, where only about 5 % of the observed haplogroups were of sub-Saharan provenance. We compared the genetic distances based on population pairwise F (st) values between 11 published East, Central and North African as well as western Asian populations and the Somali sequences and displayed them in a multi-dimensional scaling plot. Genetic proximity evidenced by clustering roughly reflected the relative geographic location of the populations. The sequences will be included in the EMPOP database ( www.empop.org ) under accession number EMP00397 upon publication (Parson and Dür Forensic Sci Int Genet 1:88-92, 2007).

Accumulation of mutations over the entire mitochondrial genome of breast cancer cells obtained by tissue microdissection.

The occurrence of heteroplasmy and mixtures is technically challenging for the analysis of mitochondrial DNA. More than that, observed mutations need to be carefully interpreted in the light of the phylogeny as mitochondrial DNA is a uniparental marker reflecting human evolution. Earlier attempts to explain the role of mtDNA in cancerous tissues led to substantial confusion in medical genetics mainly due to the presentation of low sequence data quality and misinterpretation of mutations representing a particular haplogroup background rather than being cancer-specific. The focus of this study is to characterize the extent and level of mutations in breast cancer samples obtained by tissue microdissection by application of an evaluated full mtDNA genome sequencing protocol. We amplified and sequenced the complete mitochondrial genomes of microdissected breast cancer cells of 15 patients and compared the results to those obtained from paired non-cancerous breast tissue derived from the same patients. We observed differences in the heteroplasmic states of substitutions between cancerous and normal cells, one of which was affecting a position that has been previously reported in lung cancer and another one that has been identified in 16 epithelial ovarian tumors, possibly indicating functional relevance. In the coding region, we found full transitions in two cancerous mitochondrial genomes and 12 heteroplasmic substitutions as compared to the non-cancerous breast cells. We identified somatic mutations over the entire mtDNA of human breast cancer cells potentially impairing the mitochondrial OXPHOS system.

Implications of Standardized Uptake Values of Oral Squamous Cell Carcinoma in PET-CT on Prognosis, Tumor Characteristics and Mitochondrial DNA Heteroplasmy.

Under aerobic conditions, some cancers switch to glycolysis to cover their energy requirements. Taking advantage of this process, functional imaging techniques such as PET-CT can be used to detect and assess tumorous tissues. The aim of this study was to investigate standardized uptake values and mitochondrial DNA mutations in oral squamous cell carcinoma. A cohort of 57 patients underwent 18[F]FDG-PET-CT and standardized uptake values were collected. In 15 patients, data on mitochondrial DNA mutations of the tumor were available. Kaplan-Meier curves were calculated, and correlation analyses as well as univariate Cox proportional hazard models were performed. Using ROC analysis to determine a statistical threshold for SUVmax in PET investigations, a cut-off value was determined at 9.765 MB/mL. Survival analysis for SUVmax in these groups showed a Hazard Ratio of 4 (95% CI 1.7-9) in the high SUVmax group with 5-year survival rates of 23.5% (p = 0.00042). For SUVmax and clinicopathological tumor features, significant correlations were found. A tendency towards higher mtDNA heteroplasmy levels in high SUVmax groups could be observed. We were able to confirm the prognostic value of SUVmax in OSCC, showing higher survival rates at lower SUVmax levels. Correlations between SUVmax and distinct tumor characteristics were highly significant, providing evidence that SUVmax may act as a reliable diagnostic parameter. Correlation analysis of mtDNA mutations suggests an influence on metabolic activity in OSCC.

OXPHOS remodeling in high-grade prostate cancer involves mtDNA mutations and increased succinate oxidation.

Rewiring of energy metabolism and adaptation of mitochondria are considered to impact on prostate cancer development and progression. Here, we report on mitochondrial respiration, DNA mutations and gene expression in paired benign/malignant human prostate tissue samples. Results reveal reduced respiratory capacities with NADH-pathway substrates glutamate and malate in malignant tissue and a significant metabolic shift towards higher succinate oxidation, particularly in high-grade tumors. The load of potentially deleterious mitochondrial-DNA mutations is higher in tumors and associated with unfavorable risk factors. High levels of potentially deleterious mutations in mitochondrial Complex I-encoding genes are associated with a 70% reduction in NADH-pathway capacity and compensation by increased succinate-pathway capacity. Structural analyses of these mutations reveal amino acid alterations leading to potentially deleterious effects on Complex I, supporting a causal relationship. A metagene signature extracted from the transcriptome of tumor samples exhibiting a severe mitochondrial phenotype enables identification of tumors with shorter survival times.

Profiling of Mitochondrial DNA Heteroplasmy in a Prospective Oral Squamous Cell Carcinoma Study.

While a shift in energy metabolism is essential to cancers, the knowledge about the involvement of the mitochondrial genome in tumorigenesis and progression in oral squamous cell carcinoma (OSCC) is still very limited. In this study, we evaluated 37 OSCC tumors and the corresponding benign mucosa tissue pairs by deep sequencing of the complete mitochondrial DNA (mtDNA). After extensive quality control, we identified 287 variants, 137 in tumor and 150 in benign samples exceeding the 1% threshold. Variant heteroplasmy levels were significantly increased in cancer compared to benign tissues (p = 0.0002). Furthermore, pairwise high heteroplasmy frequency difference variants (∆HF% > 20) with potential functional impact were increased in the cancer tissues (p = 0.024). Fourteen mutations were identified in the protein-coding region, out of which thirteen were detected in cancer and only one in benign tissue. After eight years of follow-up, the risk of mortality was higher for patients who harbored at least one ∆HF% > 20 variant in mtDNA protein-coding regions relative to those with no mutations (HR = 4.6, (95%CI = 1.3-17); p = 0.019 in primary tumor carriers). Haplogroup affiliation showed an impact on survival time, which however needs confirmation in a larger study. In conclusion, we observed a significantly higher accumulation of somatic mutations in the cancer tissues associated with a worse prognosis.

Sequencing strategy for the whole mitochondrial genome resulting in high quality sequences.

BACKGROUND: It has been demonstrated that a reliable and fail-safe sequencing strategy is mandatory for high-quality analysis of mitochondrial (mt) DNA, as the sequencing and base-calling process is prone to error. Here, we present a high quality, reliable and easy handling manual procedure for the sequencing of full mt genomes that is also appropriate for laboratories where fully automated processes are not available. RESULTS: We amplified whole mitochondrial genomes as two overlapping PCR-fragments comprising each about 8500 bases in length. We developed a set of 96 primers that can be applied to a (manual) 96 well-based technology, which resulted in at least double strand sequence coverage of the entire coding region (codR). CONCLUSION: This elaborated sequencing strategy is straightforward and allows for an unambiguous sequence analysis and interpretation including sometimes challenging phenomena such as point and length heteroplasmy that are relevant for the investigation of forensic and clinical samples.

Forensic and phylogeographic characterization of mtDNA lineages from northern Thailand (Chiang Mai).

The immigration of diverse ethnic groups over the past centuries from surrounding countries into Thailand left footprints in the genetic composition of Thai mitochondrial DNA (mtDNA) lineages. The entire mtDNA control region (1,122 bp) was typed in 190 unrelated male volunteers from the northern Thailand province of Chiang Mai following highest quality standards. For a more precise haplogroup classification, selected single nucleotide polymorphisms from the mtDNA coding region were genotyped. We found several new, so far undescribed mtDNA lineages. Quasi-median networks were constructed for visualisation of character conflicts. The data were put into population-genetic relationships with other Southeast Asian populations. Although the frequencies of the Thai haplogroups were characteristic for Southeast Asia in terms of haplotype composition and genetic structure, the Thai population was significantly different from other Southeast Asian populations. This necessitates establishing regional databases, especially for forensic applications. The population data have been submitted to the EMPOP database (www.empop.org) and will be available on publication.

Plasmid-normalized quantification of relative mitochondrial DNA copy number.

Alterations of mitochondrial DNA (mtDNA) copy number have been associated with a wide variety of phenotypes and diseases. Unfortunately, the literature provides scarce methodical information about duplex targeting of nuclear and mtDNA that meets the quality criteria for qPCR. Therefore, we established a method for mtDNA copy number quantification using a quantitative PCR assay that allows for simultaneous targeting of a single copy nuclear gene (beta-2-microglobulin) and the t-RNALeu gene on the mtDNA. We include a plasmid containing both targets in order to normalize against differences in emission intensities of the fluorescent dyes Yakima Yellow and FAM. Applying the plasmid calibrator on an internal control reduced the intra-assay variability from 21% (uncorrected) to 7% (plasmid-corrected). Moreover, we noted that DNA samples isolated with different methods revealed different numbers of mtDNA copies, thus highlighting an important influence of the pre-analytical procedures. In summary, we developed a precise assay for mitochondrial copy number detection relative to nuclear DNA. Our method is applicable to comparative mitochondrial DNA copy number studies since the use of the dual insert plasmid allows correcting for the unequal emission intensities of the different fluorescent labels of the two targets.

Validation of Next-Generation Sequencing of Entire Mitochondrial Genomes and the Diversity of Mitochondrial DNA Mutations in Oral Squamous Cell Carcinoma.

BACKGROUND: Oral squamous cell carcinoma (OSCC) is mainly caused by smoking and alcohol abuse and shows a five-year survival rate of ~50%. We aimed to explore the variation of somatic mitochondrial DNA (mtDNA) mutations in primary oral tumors, recurrences and metastases. METHODS: We performed an in-depth validation of mtDNA next-generation sequencing (NGS) on an Illumina HiSeq 2500 platform for its application to cancer tissues, with the goal to detect low-level heteroplasmies and to avoid artifacts. Therefore we genotyped the mitochondrial genome (16.6 kb) from 85 tissue samples (tumors, recurrences, resection edges, metastases and blood) collected from 28 prospectively recruited OSCC patients applying both Sanger sequencing and high-coverage NGS (~35,000 reads per base). RESULTS: We observed a strong correlation between Sanger sequencing and NGS in estimating the mixture ratio of heteroplasmies (r = 0.99; p<0.001). Non-synonymous heteroplasmic variants were enriched among cancerous tissues. The proportions of somatic and inherited variants in a given gene region were strongly correlated (r = 0.85; p<0.001). Half of the patients shared mutations between benign and cancerous tissue samples. Low level heteroplasmies (<10%) were more frequent in benign samples compared to tumor samples, where heteroplasmies >10% were predominant. Four out of six patients who developed a local tumor recurrence showed mutations in the recurrence that had also been observed in the primary tumor. Three out of five patients, who had tumor metastases in the lymph nodes of their necks, shared mtDNA mutations between primary tumors and lymph node metastases. The percentage of mutation heteroplasmy increased from the primary tumor to lymph node metastases. CONCLUSIONS: We conclude that Sanger sequencing is valid for heteroplasmy quantification for heteroplasmies ≥10% and that NGS is capable of reliably detecting and quantifying heteroplasmies down to the 1%-level. The finding of shared mutations between primary tumors, recurrences and metastasis indicates a clonal origin of malignant cells in oral cancer.

Full mtGenome reference data: development and characterization of 588 forensic-quality haplotypes representing three U.S. populations.

Though investigations into the use of massively parallel sequencing technologies for the generation of complete mitochondrial genome (mtGenome) profiles from difficult forensic specimens are well underway in multiple laboratories, the high quality population reference data necessary to support full mtGenome typing in the forensic context are lacking. To address this deficiency, we have developed 588 complete mtGenome haplotypes, spanning three U.S. population groups (African American, Caucasian and Hispanic) from anonymized, randomly-sampled specimens. Data production utilized an 8-amplicon, 135 sequencing reaction Sanger-based protocol, performed in semi-automated fashion on robotic instrumentation. Data review followed an intensive multi-step strategy that included a minimum of three independent reviews of the raw data at two laboratories; repeat screenings of all insertions, deletions, heteroplasmies, transversions and any additional private mutations; and a check for phylogenetic feasibility. For all three populations, nearly complete resolution of the haplotypes was achieved with full mtGenome sequences: 90.3-98.8% of haplotypes were unique per population, an improvement of 7.7-29.2% over control region sequencing alone, and zero haplotypes overlapped between populations. Inferred maternal biogeographic ancestry frequencies for each population and heteroplasmy rates in the control region were generally consistent with published datasets. In the coding region, nearly 90% of individuals exhibited length heteroplasmy in the 12418-12425 adenine homopolymer; and despite a relatively high rate of point heteroplasmy (23.8% of individuals across the entire molecule), coding region point heteroplasmies shared by more than one individual were notably absent, and transversion-type heteroplasmies were extremely rare. The ratio of nonsynonymous to synonymous changes among point heteroplasmies in the protein-coding genes (1:1.3) and average pathogenicity scores in comparison to data reported for complete substitutions in previous studies seem to provide some additional support for the role of purifying selection in the evolution of the human mtGenome. Overall, these thoroughly vetted full mtGenome population reference data can serve as a standard against which the quality and features of future mtGenome datasets (especially those developed via massively parallel sequencing) may be evaluated, and will provide a solid foundation for the generation of complete mtGenome haplotype frequency estimates for forensic applications.

WITHDRAWN: Erratum to "Evaluation of next generation mtGenome sequencing using the Ion Torrent Personal Genome Machine (PGM)" [Forensic Sci. Int.: Genet. 7 (2013) 543-549].

The Publisher regrets that this article is an accidental duplication of an article that has already been published, http://dx.doi.org/10.1016/j.fsigen.2013.09.007. The duplicate article has therefore been withdrawn. The full Elsevier Policy on Article Withdrawal can be found at http://www.elsevier.com/locate/withdrawalpolicy.

Evaluation of next generation mtGenome sequencing using the Ion Torrent Personal Genome Machine (PGM).

Insights into the human mitochondrial phylogeny have been primarily achieved by sequencing full mitochondrial genomes (mtGenomes). In forensic genetics (partial) mtGenome information can be used to assign haplotypes to their phylogenetic backgrounds, which may, in turn, have characteristic geographic distributions that would offer useful information in a forensic case. In addition and perhaps even more relevant in the forensic context, haplogroup-specific patterns of mutations form the basis for quality control of mtDNA sequences. The current method for establishing (partial) mtDNA haplotypes is Sanger-type sequencing (STS), which is laborious, time-consuming, and expensive. With the emergence of Next Generation Sequencing (NGS) technologies, the body of available mtDNA data can potentially be extended much more quickly and cost-efficiently. Customized chemistries, laboratory workflows and data analysis packages could support the community and increase the utility of mtDNA analysis in forensics. We have evaluated the performance of mtGenome sequencing using the Personal Genome Machine (PGM) and compared the resulting haplotypes directly with conventional Sanger-type sequencing. A total of 64mtGenomes (>1 million bases) were established that yielded high concordance with the corresponding STS haplotypes (<0.02% differences). About two-thirds of the differences were observed in or around homopolymeric sequence stretches. In addition, the sequence alignment algorithm employed to align NGS reads played a significant role in the analysis of the data and the resulting mtDNA haplotypes. Further development of alignment software would be desirable to facilitate the application of NGS in mtDNA forensic genetics.

Reprint of: Evaluation of next generation mtGenome sequencing using the Ion Torrent Personal Genome Machine (PGM).

Insights into the human mitochondrial phylogeny have been primarily achieved by sequencing full mitochondrial genomes (mtGenomes). In forensic genetics (partial) mtGenome information can be used to assign haplotypes to their phylogenetic backgrounds, which may, in turn, have characteristic geographic distributions that would offer useful information in a forensic case. In addition and perhaps even more relevant in the forensic context, haplogroup-specific patterns of mutations form the basis for quality control of mtDNA sequences. The current method for establishing (partial) mtDNA haplotypes is Sanger-type sequencing (STS), which is laborious, time-consuming, and expensive. With the emergence of Next Generation Sequencing (NGS) technologies, the body of available mtDNA data can potentially be extended much more quickly and cost-efficiently. Customized chemistries, laboratory workflows and data analysis packages could support the community and increase the utility of mtDNA analysis in forensics. We have evaluated the performance of mtGenome sequencing using the Personal Genome Machine (PGM) and compared the resulting haplotypes directly with conventional Sanger-type sequencing. A total of 64mtGenomes (>1 million bases) were established that yielded high concordance with the corresponding STS haplotypes (<0.02% differences). About two-thirds of the differences were observed in or around homopolymeric sequence stretches. In addition, the sequence alignment algorithm employed to align NGS reads played a significant role in the analysis of the data and the resulting mtDNA haplotypes. Further development of alignment software would be desirable to facilitate the application of NGS in mtDNA forensic genetics.

Mitochondrial DNA control region data from indigenous Angolan Khoe-San lineages.

Here we provide 129 complete mitochondrial control region sequences of indigenous Khoe-San individuals from Angola to contribute to the still underrepresented pool of data from Africa. The dataset consists of exclusively African lineages with a majority of Sub-Saharan haplogroups. The probability of a random match was calculated as 0.09. The data set comprises 21 haplotypes occurring more than once and 17 unique haplotypes. Upon publication, haplotypes were incorporated in the EMPOP database (www.empop.org; EMP00069) [1].

MtDNA diversity of Ghana: a forensic and phylogeographic view.

West Africa is characterized by a migration history spanning more than 150,000 years. Climate changes but also political circumstances were responsible for several early but also recent population movements that shaped the West African mitochondrial landscape. The aim of the study was to establish a Ghanaian mtDNA dataset for forensic purposes and to investigate the diversity of the Ghanaian population sample with respect to surrounding populations. We sequenced full mitochondrial control regions of 193 Akan people from Ghana and excluded two apparently close maternally related individuals due to preceding kinship testing. The remaining dataset comprising 191 sequences was applied as etalon for quasi-median network analysis and was subsequently combined with 99 additional control region sequences from surrounding West African countries. All sequences were incorporated into the EMPOP database enriching the severely underrepresented African mtDNA pool. For phylogeographic considerations, the Ghanaian haplotypes were compared to those of 19 neighboring populations comprising a total number of 6198 HVS1 haplotypes. We found extensive genetic admixture between the Ghanaian lineages and those from adjacent populations diminishing with geographical distance. The extent of genetic admixture reflects the long but also recent history of migration waves within West Africa mainly caused by changing environmental conditions. Also, evidence for potential socio-economical influences such as trade routes is provided by the occurrence of U6b and U6d sequences found in Dubai but also in Tunisia leading to the African West Coast via Mauritania and Senegal but also via Niger, Nigeria to Cameroon.

Polymorphisms of mtDNA control region in Tunisian and Moroccan populations: an enrichment of forensic mtDNA databases with Northern Africa data.

Current forensic mitochondrial (mt)DNA databases are limited in representative population data of African origin. We investigated HVS-I/HVS-II sequences of 120 Tunisian and Moroccan healthy male donors applying stringent quality criteria to assure high quality of the data and phylogenetic alignment and notation of the sequences. Among 64 Tunisians, 56 different haplotypes were observed and the most common haplotype (16187T 16189C 16223T 16264T 16270T 16278T 16293G 16311C 73G 152C 182T 185T 195C 247A 263G 309.1C 315.1C; haplogroup (hg) L1b) was shared by four individuals. 56 Moroccans could be assigned to 52 different haplotypes where the most common haplotype was of West Eurasian origin with the hg H sequence motif 263G 315.1C and variations in the HVS-II polyC-stretch (309.1C 309.2C) shared by six samples. The majority of the observed haplotypes belong to the west Eurasian phylogeny (50% in Tunisians and 62.5% in Moroccans). Our data are consistent with the current phylogeographic knowledge displaying the occurrence of sub-Saharan haplogroup L sequences, found in 48.4% of Tunisians and 25% of Moroccans as well as the presence of the two re-migrated haplogroups U6 (7.8% and 1.8% in Tunisians and Moroccans, respectively) and M1 (1.6% in Tunisians and 8.9% in Moroccans).