Artificial recombination in forensic mtDNA population databases

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.

Whole Mitochondrial DNA Sequencing Analysis in 47 Han Populations in Southwest China

Background

Mitochondrial DNA (mtDNA) sequencing has been used in many areas, including forensic genetics. Due to the rapid development of sequencing technology, whole mtDNA sequencing is now possible and may be used in epidemiological and forensic studies. This study aimed to use whole mtDNA sequencing to investigate 47 Chongqing Han populations in southwest China and the diversity in the mtGenome reference data.

Material/Methods

The mtDNA of 47 Chongqing Han populations was generated using the Ion Torrent Personal Genome Machine (PGM) system. The extent of the effects of the mtDNA on the subpopulations was investigated and compared with six other populations from published studies. Pairwise fixation index (FST), a measure of population differentiation due to genetic structure, were calculated. Analysis of molecular variance (AMOVA) was performed, and 1257 hypervariable region data sets were added to the principal component analysis (PCA).

Results

The whole mtDNA sequencing data of 47 southwest Chinese Han populations were successfully recovered. Expanding the sequencing rage increased the discrimination power of mtDNA from three-times to 25-times based on different populations. The subpopulation effects showed 20 times the differences in match probability when compared with south China regions.

Conclusions

Whole mtDNA sequencing distinguished between individuals from 47 Chongqing Han populations in southwest China and has potential applications that include high-quality forensic identification.

Concordance of mitochondrial DNA sequencing methods on bloodstains using Ion PGM™

In this study, the complete mitochondrial genome (mtGenome) of six samples from three forensic cases was sequenced using the Ion Torrent Personal Genome Machine (PGM). The analyzed samples from forensic cases included bloodstains from several materials, such as gauze, Flinder's Technology Associates (FTA) cards and swabs. The age of the samples ranged from two months to twelve years. The complete mtGenomes were amplified using the tiling sequencing strategy which divided the whole mtGenome into 162 amplicons. All amplicons were successfully recovered. A phylogenetic analysis was performed to determine the accuracy of the PGM data, and which were compared to partial Sanger-based sequencing data. The average coverage of the PGM data were above 4000× in all case samples, and 99.86% concordance was observed using both sequencing methods. In conclusion, we demonstrate the ability to recover the complete mtGenome from bloodstains with relatively poor DNA quality by PGM. Moreover, the results are concordant with Sanger sequencing data. This new method has potential use in forensic practice.

Optimized mtDNA Control Region Primer Extension Capture Analysis for Forensically Relevant Samples and Highly Compromised mtDNA of Different Age and Origin

The analysis of mitochondrial DNA (mtDNA) has proven useful in forensic genetics and ancient DNA (aDNA) studies, where specimens are often highly compromised and DNA quality and quantity are low. In forensic genetics, the mtDNA control region (CR) is commonly sequenced using established Sanger-type Sequencing (STS) protocols involving fragment sizes down to approximately 150 base pairs (bp). Recent developments include Massively Parallel Sequencing (MPS) of (multiplex) PCR-generated libraries using the same amplicon sizes. Molecular genetic studies on archaeological remains that harbor more degraded aDNA have pioneered alternative approaches to target mtDNA, such as capture hybridization and primer extension capture (PEC) methods followed by MPS. These assays target smaller mtDNA fragment sizes (down to 50 bp or less), and have proven to be substantially more successful in obtaining useful mtDNA sequences from these samples compared to electrophoretic methods. Here, we present the modification and optimization of a PEC method, earlier developed for sequencing the Neanderthal mitochondrial genome, with forensic applications in mind. Our approach was designed for a more sensitive enrichment of the mtDNA CR in a single tube assay and short laboratory turnaround times, thus complying with forensic practices. We characterized the method using sheared, high quantity mtDNA (six samples), and tested challenging forensic samples (n = 2) as well as compromised solid tissue samples (n = 15) up to 8 kyrs of age. The PEC MPS method produced reliable and plausible mtDNA haplotypes that were useful in the forensic context. It yielded plausible data in samples that did not provide results with STS and other MPS techniques. We addressed the issue of contamination by including four generations of negative controls, and discuss the results in the forensic context. We finally offer perspectives for future research to enable the validation and accreditation of the PEC MPS method for final implementation in forensic genetic laboratories.

HaploGrep 2: mitochondrial haplogroup classification in the era of high-throughput sequencing

Mitochondrial DNA (mtDNA) profiles can be classified into phylogenetic clusters (haplogroups), which is of great relevance for evolutionary, forensic and medical genetics. With the extensive growth of the underlying phylogenetic tree summarizing the published mtDNA sequences, the manual process of haplogroup classification would be too time-consuming. The previously published classification tool HaploGrep provided an automatic way to address this issue. Here, we present the completely updated version HaploGrep 2 offering several advanced features, including a generic rule-based system for immediate quality control (QC). This allows detecting artificial recombinants and missing variants as well as annotating rare and phantom mutations. Furthermore, the handling of high-throughput data in form of VCF files is now directly supported. For data output, several graphical reports are generated in real time, such as a multiple sequence alignment format, a VCF format and extended haplogroup QC reports, all viewable directly within the application. In addition, HaploGrep 2 generates a publication-ready phylogenetic tree of all input samples encoded relative to the revised Cambridge Reference Sequence. Finally, new distance measures and optimizations of the algorithm increase accuracy and speed-up the application. HaploGrep 2 can be accessed freely and without any registration at http://haplogrep.uibk.ac.at.

The complete mitogenome of a 500-year-old Inca child mummy

In 1985, a frozen mummy was found in Cerro Aconcagua (Argentina). Archaeological studies identified the mummy as a seven-year-old Inca sacrifice victim who lived >500 years ago, at the time of the expansion of the Inca Empire towards the southern cone. The sequence of its entire mitogenome was obtained. After querying a large worldwide database of mitogenomes (>28,000) we found that the Inca haplotype belonged to a branch of haplogroup C1b (C1bi) that has not yet been identified in modern Native Americans. The expansion of C1b into the Americas, as estimated using 203 C1b mitogenomes, dates to the initial Paleoindian settlements (~18.3 thousand years ago [kya]); however, its internal variation differs between Mesoamerica and South America. By querying large databases of control region haplotypes (>150,000), we found only a few C1bi members in Peru and Bolivia (e.g. Aymaras), including one haplotype retrieved from ancient DNA of an individual belonging to the Wari Empire (Peruvian Andes). Overall, the results suggest that the profile of the mummy represents a very rare sub-clade that arose 14.3 (5–23.6) kya and could have been more frequent in the past. A Peruvian Inca origin for present-day C1bi haplotypes would satisfy both the genetic and paleo-anthropological findings.

Comprehensive Analysis of Pan-African Mitochondrial DNA Variation Provides New Insights into Continental Variation and Demography

Africa is the cradle of all human beings, and although it has been the focus of a number of genetic studies, there are many questions that remain unresolved. We have performed one of the largest and most comprehensive meta-analyses of mitochondrial DNA (mtDNA) lineages carried out in the African continent to date. We generated high-throughput mtDNA single nucleotide polymorphism (SNP) data (230 SNPs) from 2024 Africans, where more than 500 of them were additionally genotyped for the control region. These data were analyzed together with over 12,700 control region profiles collected from the literature, representing more than 300 population samples from Africa. Insights into the African homeland of humans are discussed. Phylogeographic patterns for the African continent are shown at a high phylogeographic resolution as well as at the population and regional levels. The deepest branch of the mtDNA tree, haplogroup L0, shows the highest sub-haplogroup diversity in Southeast and East Africa, suggesting this region as the homeland for modern humans. Several demographic estimates point to the coast as a facilitator of human migration in Africa, but the data indicate complex patterns, perhaps mirroring the effect of recent continental-scaled demographic events in re-shaping African mtDNA variability.

The Use and Effectiveness of Triple Multiplex System for Coding Region Single Nucleotide Polymorphism in Mitochondrial DNA Typing of Archaeologically Obtained Human Skeletons from Premodern Joseon Tombs of Korea

Previous study showed that East Asian mtDNA haplogroups, especially those of Koreans, could be successfully assigned by the coupled use of analyses on coding region SNP markers and control region mutation motifs. In this study, we tried to see if the same triple multiplex analysis for coding regions SNPs could be also applicable to ancient samples from East Asia as the complementation for sequence analysis of mtDNA control region. By the study on Joseon skeleton samples, we know that mtDNA haplogroup determined by coding region SNP markers successfully falls within the same haplogroup that sequence analysis on control region can assign. Considering that ancient samples in previous studies make no small number of errors in control region mtDNA sequencing, coding region SNP analysis can be used as good complimentary to the conventional haplogroup determination, especially of archaeological human bone samples buried underground over long periods.

Complete mitochondrial genome database and standardized classification system for Canis lupus familiaris

To contribute to the complete mitogenome database of the species Canis lupus familiaris and shed more light on its origin, we have sequenced mitochondrial genomes of 120 modern dogs from worldwide populations. Together with all the previously published mitogenome sequences of acceptable quality, we have reconstructed a global phylogenetic tree of 555 C. l. familiaris mitogenomes and standardized haplogroup nomenclature. The phylogenetic tree presented here and available online at http://clf.mtdna.tree.cm.umk.pl/ could be further used by forensic and evolutionary geneticists as well cynologists, for data quality control and unambiguous haplogroup classification. Our in-depth phylogeographic analysis of all C. l. familiaris mitogenomes confirmed that domestic dogs may have originated in East Asia during the Mesolithic and Upper Paleolithic time periods and started to expand to other parts of the world during Neolithic times.

Reviewing population studies for forensic purposes: Dog mitochondrial DNA

The identification of dog hair through mtDNA analysis has become increasingly important in the last 15 years, as it can provide associative evidence connecting victims and suspects. The evidential value of an mtDNA match between dog hair and its potential donor is determined by the random match probability of the haplotype. This probability is based on the haplotype’s population frequency estimate. Consequently, implementing a population study representative of the population relevant to the forensic case is vital to the correct evaluation of the evidence. This paper reviews numerous published dog mtDNA studies and shows that many of these studies vary widely in sampling strategies and data quality. Therefore, several features influencing the representativeness of a population sample are discussed. Moreover, recommendations are provided on how to set up a dog mtDNA population study and how to decide whether or not to include published data. This review emphasizes the need for improved dog mtDNA population data for forensic purposes, including targeting the entire mitochondrial genome. In particular, the creation of a publicly available database of qualitative dog mtDNA population studies would improve the genetic analysis of dog traces in forensic casework.

The genetic legacy of the pre-colonial period in contemporary Bolivians

Only a few genetic studies have been carried out to date in Bolivia. However, some of the most important (pre)historical enclaves of South America were located in these territories. Thus, the (sub)-Andean region of Bolivia was part of the Inca Empire, the largest state in Pre-Columbian America. We have genotyped the first hypervariable region (HVS-I) of 720 samples representing the main regions in Bolivia, and these data have been analyzed in the context of other pan-American samples (>19,000 HVS-I mtDNAs). Entire mtDNA genome sequencing was also undertaken on selected Native American lineages. Additionally, a panel of 46 Ancestry Informative Markers (AIMs) was genotyped in a sub-set of samples. The vast majority of the Bolivian mtDNAs (98.4%) were found to belong to the main Native American haplogroups (A: 14.3%, B: 52.6%, C: 21.9%, D: 9.6%), with little indication of sub-Saharan and/or European lineages; however, marked patterns of haplogroup frequencies between main regions exist (e.g. haplogroup B: Andean [71%], Sub-Andean [61%], Llanos [32%]). Analysis of entire genomes unraveled the phylogenetic characteristics of three Native haplogroups: the pan-American haplogroup B2b (originated ∼21.4 thousand years ago [kya]), A2ah (∼5.2 kya), and B2o (∼2.6 kya). The data suggest that B2b could have arisen in North California (an origin even in the north most region of the American continent cannot be disregarded), moved southward following the Pacific coastline and crossed Meso-America. Then, it most likely spread into South America following two routes: the Pacific path towards Peru and Bolivia (arriving here at about ∼15.2 kya), and the Amazonian route of Venezuela and Brazil southwards. In contrast to the mtDNA, Ancestry Informative Markers (AIMs) reveal a higher (although geographically variable) European introgression in Bolivians (25%). Bolivia shows a decreasing autosomal molecular diversity pattern along the longitudinal axis, from the Altiplano to the lowlands. Both autosomes and mtDNA revealed a low impact (1-2%) of a sub-Saharan component in Bolivians.

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.

The maintenance of mitochondrial genetic stability is crucial during the oncogenic process

The main energetic resources of the cell are the mitochondria. As such, these organelles control a number of processes related to the life and death of the cell and also have a prominent function in the maintenance of tumor cells. In the last years, several authors have proposed an active role for mitochondria in tumorigenesis, more specifically concerning somatic mutations in mitochondrial DNA (mtDNA). Here, we wanted to evaluate this hypothesis based on the conclusions obtained in a model of gliomagenesis with elevated levels of ROS (reactive oxygen species), a toxic by-product of tumor metabolism. According to our findings, none of the mtDNA variants were found relevant to the tumoral process or suggest the involvement of mitochondria in tumorigenesis beyond the metabolic requirements of the tumoral cell. We conclude that there is not enough evidence to support the claim that mitochondrial instability holds any relevant role in the tumoral process.

Haplogrouping mitochondrial DNA sequences in Legal Medicine/Forensic Genetics

Haplogrouping refers to the classification of (partial) mitochondrial DNA (mtDNA) sequences into haplogroups using the current knowledge of the worldwide mtDNA phylogeny. Haplogroup assignment of mtDNA control-region sequences assists in the focused comparison with closely related complete mtDNA sequences and thus serves two main goals in forensic genetics: first is the a posteriori quality analysis of sequencing results and second is the prediction of relevant coding-region sites for confirmation or further refinement of haplogroup status. The latter may be important in forensic casework where discrimination power needs to be as high as possible. However, most articles published in forensic genetics perform haplogrouping only in a rudimentary or incorrect way. The present study features PhyloTree as the key tool for assigning control-region sequences to haplogroups and elaborates on additional Web-based searches for finding near-matches with complete mtDNA genomes in the databases. In contrast, none of the automated haplogrouping tools available can yet compete with manual haplogrouping using PhyloTree plus additional Web-based searches, especially when confronted with artificial recombinants still present in forensic mtDNA datasets. We review and classify the various attempts at haplogrouping by using a multiplex approach or relying on automated haplogrouping. Furthermore, we re-examine a few articles in forensic journals providing mtDNA population data where appropriate haplogrouping following PhyloTree immediately highlights several kinds of sequence errors.

No evidence that major mtDNA European haplogroups confer risk to schizophrenia

Previous studies suggest that genetic factors could be involved in mitochondrial dysfunction observed in schizophrenia (SZ), some of them claiming a role of mtDNA common variants (mtSNPs) and/or haplogroups (hgs) in developing this disorder. These studies, however, have mainly been undertaken on relatively small cohorts of patients and control individuals and most have not yet been replicated. To further analyze the role of mtSNPs in SZ risk, we have carried out the largest genotyping effort to date using two Spanish case-control samples comprising a total of 942 schizophrenic patients and 1,231 unrelated controls: 454 patients and 616 controls from Santiago de Compostela (Galicia) and 488 patients and 615 controls from Reus (Catalonia). A set of 25 mtSNPs representing main branches of the European mtDNA phylogeny were genotyped in the Galician cohort and a subset of 16 out of these 25 mtSNPs was genotyped in the Catalan cohort. These 16 common variants characterize the most common European branches of the mtDNA phylogeny. We did not observe any positive association of mtSNPs and hgs with SZ. We discuss several deficiencies of previous studies that might explain the false positive nature of previous findings, including the confounding effect of population sub-structure and deficient statistical methodologies. It is unlikely that mtSNPs defining the most common European mtDNA haplogroups are related to SZ.

Interdisciplinary approach to the demography of Jamaica

Background

The trans-Atlantic slave trade dramatically changed the demographic makeup of the New World, with varying regions of the African coast exploited differently over roughly a 400 year period. When compared to the discrete mitochondrial haplotype distribution of historically appropriate source populations, the unique distribution within a specific source population can prove insightful in estimating the contribution of each population. Here, we analyzed the first hypervariable region of mitochondrial DNA in a sample from the Caribbean island of Jamaica and compared it to aggregated populations in Africa divided according to historiographically defined segments of the continent's coastline. The results from these admixture procedures were then compared to the wealth of historic knowledge surrounding the disembarkation of Africans on the island.

Results

In line with previous findings, the matriline of Jamaica is almost entirely of West African descent. Results from the admixture analyses suggest modern Jamaicans share a closer affinity with groups from the Gold Coast and Bight of Benin despite high mortality, low fecundity, and waning regional importation. The slaves from the Bight of Biafra and West-central Africa were imported in great numbers; however, the results suggest a deficit in expected maternal contribution from those regions.

Conclusions

When considering the demographic pressures imposed by chattel slavery on Jamaica during the slave era, the results seem incongruous. Ethnolinguistic and ethnographic evidence, however, may explain the apparent non-random levels of genetic perseverance. The application of genetics may prove useful in answering difficult demographic questions left by historically voiceless groups.

Interpretation guidelines of mtDNA control region sequence electropherograms in forensic genetics

Forensic mitochondrial DNA (mtDNA) analysis is a complementary technique to forensic nuclear DNA (nDNA) and trace evidence analysis. Its use has been accepted by the vast majority of courts of law around the world. However for the forensic community it is crucial to employ standardized methods and procedures to guaranty the quality of the results obtained in court. In this chapter, we describe the most important aspects regarding the interpretation and assessment of mtDNA analysis, and offer a simple guide which places particular emphasis on those aspects that can impact the final interpretation of the results. These include the criteria for authenticating a sequence excluding the contaminant origin, defining the quality of a sequence, editing procedure, alignment criteria for searching the databases, and the statistical evaluation of matches. It is not easy to establish a single guide to interpretation for mtDNA analysis; however, it is important to understand all variables that may in some way affect the final conclusion in the context of a forensic case. As a general rule, laboratories should be cautious before issuing the final conclusion of an mtDNA analysis, and consider any significant limitations regarding current understanding of specific aspects of the mtDNA molecule.

Mitochondria in anthropology and forensic medicine

Mitochondria's role in crucial metabolic pathways is probably the first answer which comes to our minds for the question: what do these tiny organelles serve for? However, specific features of their DNA made them extremely useful also in the field of anthropology and forensics. MtDNA analyses became a milestone in the complex task of unraveling earliest human migrations. Evidence provided by these experiments left no doubts on modern humans origins pointing to Africa being our cradle. It also contributed to interpretation of putative ways of our dispersal around Asia and Americas thousands years ago. On the other hand, analysis of mtDNA is well established and valuable tool in forensic genetics. When other definitely more popular markers give no answer on identity, it is the time to employ information carried by mitochondria. This chapter summarizes not only current reports on the role of mitochondria in forensics and reconstruction of modern humans phylogeny, but also calls one's attention to a broad range of difficulties and constraints associated with mtDNA analyses.

A statistical framework for the interpretation of mtDNA mixtures: forensic and medical applications

Background

Mitochondrial DNA (mtDNA) variation is commonly analyzed in a wide range of different biomedical applications. Cases where more than one individual contribute to a stain genotyped from some biological material give rise to a mixture. Most forensic mixture cases are analyzed using autosomal markers. In rape cases, Y-chromosome markers typically add useful information. However, there are important cases where autosomal and Y-chromosome markers fail to provide useful profiles. In some instances, usually involving small amounts or degraded DNA, mtDNA may be the only useful genetic evidence available. Mitochondrial DNA mixtures also arise in studies dealing with the role of mtDNA variation in tumorigenesis. Such mixtures may be generated by the tumor, but they could also originate in vitro due to inadvertent contamination or a sample mix-up.

Methods/Principal Findings

We present the statistical methods needed for mixture interpretation and emphasize the modifications required for the more well-known methods based on conventional markers to generalize to mtDNA mixtures. Two scenarios are considered. Firstly, only categorical mtDNA data is assumed available, that is, the variants contributing to the mixture. Secondly, quantitative data (peak heights or areas) on the allelic variants are also accessible. In cases where quantitative information is available in addition to allele designation, it is possible to extract more precise information by using regression models. More precisely, using quantitative information may lead to a unique solution in cases where the qualitative approach points to several possibilities. Importantly, these methods also apply to clinical cases where contamination is a potential alternative explanation for the data.

Conclusions/Significance

We argue that clinical and forensic scientists should give greater consideration to mtDNA for mixture interpretation. The results and examples show that the analysis of mtDNA mixtures contributes substantially to forensic casework and may also clarify erroneous claims made in clinical genetics regarding tumorigenesis.

Current genetic methodologies in the identification of disaster victims and in forensic analysis

This review presents the basic problems and currently available molecular techniques used for genetic profiling in disaster victim identification (DVI). The environmental conditions of a mass disaster often result in severe fragmentation, decomposition and intermixing of the remains of victims. In such cases, traditional identification based on the anthropological and physical characteristics of the victims is frequently inconclusive. This is the reason why DNA profiling became the gold standard for victim identification in mass-casualty incidents (MCIs) or any forensic cases where human remains are highly fragmented and/or degraded beyond recognition. The review provides general information about the sources of genetic material for DNA profiling, the genetic markers routinely used during genetic profiling (STR markers, mtDNA and single-nucleotide polymorphisms [SNP]) and the basic statistical approaches used in DNA-based disaster victim identification. Automated technological platforms that allow the simultaneous analysis of a multitude of genetic markers used in genetic identification (oligonucleotide microarray techniques and next-generation sequencing) are also presented. Forensic and population databases containing information on human variability, routinely used for statistical analyses, are discussed. The final part of this review is focused on recent developments, which offer particularly promising tools for forensic applications (mRNA analysis, transcriptome variation in individuals/populations and genetic profiling of specific cells separated from mixtures).

Current next generation sequencing technology may not meet forensic standards

In a Nature paper of 2010, the concern was raised that intra-individual mtDNA variation may be more pronounced than previously believed, in that heteroplasmies are common and vary markedly from tissue to tissue. This claim taken at face value would have considerable impact on forensic casework. It turns out however that the employed technology detected the germ-line variation relative to the reference sequence only incompletely: on average at least five mutations were missed per sample, as an in silico reassessment of the data reveals. Before one can really set out to access to entire mtDNA genome data with relative ease for forensic purposes, one needs careful calibration studies under strict forensic conditions-or might have to wait for another generation.

eCOMPAGT integrates mtDNA: import, validation and export of mitochondrial DNA profiles for population genetics, tumour dynamics and genotype-phenotype association studies

Background

Mitochondrial DNA (mtDNA) is widely being used for population genetics, forensic DNA fingerprinting and clinical disease association studies. The recent past has uncovered severe problems with mtDNA genotyping, not only due to the genotyping method itself, but mainly to the post-lab transcription, storage and report of mtDNA genotypes.

Description

eCOMPAGT, a system to store, administer and connect phenotype data to all kinds of genotype data is now enhanced by the possibility of storing mtDNA profiles and allowing their validation, linking to phenotypes and export as numerous formats. mtDNA profiles can be imported from different sequence evaluation programs, compared between evaluations and their haplogroup affiliations stored. Furthermore, eCOMPAGT has been improved in its sophisticated transparency (support of MySQL and Oracle), security aspects (by using database technology) and the option to import, manage and store genotypes derived from various genotyping methods (SNPlex, TaqMan, and STRs). It is a software solution designed for project management, laboratory work and the evaluation process all-in-one.

Conclusions

The extended mtDNA version of eCOMPAGT was designed to enable error-free post-laboratory data handling of human mtDNA profiles. This software is suited for small to medium-sized human genetic, forensic and clinical genetic laboratories. The direct support of MySQL and the improved database security options render eCOMPAGT a powerful tool to build an automated workflow architecture for several genotyping methods. eCOMPAGT is freely available at http://dbis-informatik.uibk.ac.at/ecompagt.

Investigating the role of mitochondrial haplogroups in genetic predisposition to meningococcal disease

BACKGROUND AND AIMS

Meningococcal disease remains one of the most important infectious causes of death in industrialized countries. The highly diverse clinical presentation and prognosis of Neisseria meningitidis infections are the result of complex host genetics and environmental interactions. We investigated whether mitochondrial genetic background contributes to meningococcal disease (MD) susceptibility.

METHODOLOGY/PRINCIPAL FINDINGS

Prospective controlled study was performed through a national research network on MD that includes 41 Spanish hospitals. Cases were 307 paediatric patients with confirmed MD, representing the largest series of MD patients analysed to date. Two independent sets of ethnicity-matched control samples (CG1 [N = 917]), and CG2 [N = 616]) were used for comparison. Cases and controls underwent mtDNA haplotyping of a selected set of 25 mtDNA SNPs (mtSNPs), some of them defining major European branches of the mtDNA phylogeny. In addition, 34 ancestry informative markers (AIMs) were genotyped in cases and CG2 in order to monitor potential hidden population stratification. Samples of known African, Native American and European ancestry (N = 711) were used as classification sets for the determination of ancestral membership of our MD patients. A total of 39 individuals were eliminated from the main statistical analyses (including fourteen gypsies) on the basis of either non-Spanish self-reported ancestry or the results of AIMs indicating a European membership lower than 95%. Association analysis of the remaining 268 cases against CG1 suggested an overrepresentation of the synonym mtSNP G11719A variant (Pearson's chi-square test; adjusted P-value = 0.0188; OR [95% CI] = 1.63 [1.22-2.18]). When cases were compared with CG2, the positive association could not be replicated. No positive association has been observed between haplogroup (hg) status of cases and CG1/CG2 and hg status of cases and several clinical variants.

CONCLUSIONS

We did not find evidence of association between mtSNPs and mtDNA hgs with MD after carefully monitoring the confounding effect of population sub-structure. MtDNA variability is particularly stratified in human populations owing to its low effective population size in comparison with autosomal markers and therefore, special care should be taken in the interpretation of seeming signals of positive associations in mtDNA case-control association studies.

High mitochondrial DNA stability in B-cell chronic lymphocytic leukemia

BACKGROUND

Chronic Lymphocytic Leukemia (CLL) leads to progressive accumulation of lymphocytes in the blood, bone marrow, and lymphatic tissues. Previous findings have suggested that the mtDNA could play an important role in CLL.

METHODOLOGY/PRINCIPAL FINDINGS

The mitochondrial DNA (mtDNA) control-region was analyzed in lymphocyte cell DNA extracts and compared with their granulocyte counterpart extract of 146 patients suffering from B-Cell CLL; B-CLL (all recruited from the Basque country). Major efforts were undertaken to rule out methodological artefacts that would render a high false positive rate for mtDNA instabilities and thus lead to erroneous interpretation of sequence instabilities. Only twenty instabilities were finally confirmed, most of them affecting the homopolymeric stretch located in the second hypervariable segment (HVS-II) around position 310, which is well known to constitute an extreme mutational hotspot of length polymorphism, as these mutations are frequently observed in the general human population. A critical revision of the findings in previous studies indicates a lack of proper methodological standards, which eventually led to an overinterpretation of the role of the mtDNA in CLL tumorigenesis.

CONCLUSIONS/SIGNIFICANCE

Our results suggest that mtDNA instability is not the primary causal factor in B-CLL. A secondary role of mtDNA mutations cannot be fully ruled out under the hypothesis that the progressive accumulation of mtDNA instabilities could finally contribute to the tumoral process. Recommendations are given that would help to minimize erroneous interpretation of sequencing results in mtDNA studies in tumorigenesis.

Human mitochondrial DNA complete amplification and sequencing: a new validated primer set that prevents nuclear DNA sequences of mitochondrial origin co-amplification

To date, there are no published primers to amplify the entire mitochondrial DNA (mtDNA) that completely prevent the amplification of nuclear DNA (nDNA) sequences of mitochondrial origin. The main goal of this work was to design, validate and describe a set of primers, to specifically amplify and sequence the complete human mtDNA, allowing the correct interpretation of mtDNA heteroplasmy in healthy and pathological samples. Validation was performed using two different approaches: (i) Basic Local Alignment Search Tool and (ii) amplification using isolated nDNA obtained from sperm cells by differential lyses. During the validation process, two mtDNA regions, with high similarity with nDNA, represent the major problematic areas for primer design. One of these could represent a non-published nuclear DNA sequence of mitochondrial origin. For two of the initially designed fragments, the amplification results reveal PCR artifacts that can be attributed to the poor quality of the DNA. After the validation, nine overlapping primer pairs to perform mtDNA amplification and 22 additional internal primers for mtDNA sequencing were obtained. These primers could be a useful tool in future projects that deal with mtDNA complete sequencing and heteroplasmy detection, since they represent a set of primers that have been tested for the non-amplification of nDNA.

Testing the performance of mtSNP minisequencing in forensic samples

There is a growing interest among forensic geneticists in developing efficient protocols for genotyping coding region mitochondrial DNA (mtDNA) SNPs (mtSNPs). Minisequencing is becoming a popular method for SNP genotyping, but it is still used by few forensic laboratories. In part, this is due to the lack of studies testing its efficiency and reproducibility when applied to real and complex forensic samples. Here we tested a minisequencing design that consists of 71 mtSNPs (in three multiplexes) that are diagnostic of known branches of the R0 phylogeny, in real forensic samples, including degraded bones and teeth, hair shafts, and serial dilutions. The fact that amplicons are short coupled with the natural efficiency of the minisequencing technique allow these assays to perform well with all the samples tested either degraded and/or those containing low DNA amount. We did not observe phylogenetic inconsistencies in the 71 mtSNP haplotypes generated, indicating that the technique is robust against potential artefacts that could arise from unintended contamination and/or spurious amplification of nuclear mtDNA pseudogenes (NUMTs).

Contamination and sample mix-up can best explain some patterns of mtDNA instabilities in buccal cells and oral squamous cell carcinoma

The study of somatic DNA instabilities constitutes a debatable topic because different causes can lead to seeming DNA alteration patterns between different cells or tissues from the same individual. Carcinogenesis or the action of a particular toxic could generate such patterns, and this is in fact the leitmotif of a number of studies on mitochondrial DNA (mtDNA) instability. Patterns of seeming instabilities could also arise from technical errors at any stage of the analysis (DNA extraction, amplification, mutation screening/sequencing, and documentation). Specifically, inadvertent DNA contamination or sample mixing would yield mosaic variation that could be erroneously interpreted as real mutation differences (instabilities) between tissues from the same individual. From the very beginning, mtDNA studies comparing cancerous to non-cancerous tissues have suffered from such mosaic results. We demonstrate here that the phylogenetic linkage of whole arrays of mtDNA mutations provides strong evidence of artificial recombination in previous studies on buccal cells and oral squamous cell carcinoma.

New population and phylogenetic features of the internal variation within mitochondrial DNA macro-haplogroup R0

BACKGROUND

R0 embraces the most common mitochondrial DNA (mtDNA) lineage in West Eurasia, namely, haplogroup H (approximately 40%). R0 sub-lineages are badly defined in the control region and therefore, the analysis of diagnostic coding region polymorphisms is needed in order to gain resolution in population and medical studies.

METHODOLOGY/PRINCIPAL FINDINGS

We sequenced the first hypervariable segment (HVS-I) of 518 individuals from different North Iberian regions. The mtDNAs belonging to R0 (approximately 57%) were further genotyped for a set of 71 coding region SNPs characterizing major and minor branches of R0. We found that the North Iberian Peninsula shows moderate levels of population stratification; for instance, haplogroup V reaches the highest frequency in Cantabria (north-central Iberia), but lower in Galicia (northwest Iberia) and Catalonia (northeast Iberia). When compared to other European and Middle East populations, haplogroups H1, H3 and H5a show frequency peaks in the Franco-Cantabrian region, declining from West towards the East and South Europe. In addition, we have characterized, by way of complete genome sequencing, a new autochthonous clade of haplogroup H in the Basque country, named H2a5. Its coalescence age, 15.6+/-8 thousand years ago (kya), dates to the period immediately after the Last Glacial Maximum (LGM).

CONCLUSIONS/SIGNIFICANCE

In contrast to other H lineages that experienced re-expansion outside the Franco-Cantabrian refuge after the LGM (e.g. H1 and H3), H2a5 most likely remained confined to this area till present days.

Median network analysis of defectively sequenced entire mitochondrial genomes from early and contemporary disease studies

Sequence analysis of the mitochondrial genome has become a routine method in the study of mitochondrial diseases. Quite often, the sequencing efforts in the search of pathogenic or disease-associated mutations are affected by technical and interpretive problems, caused by sample mix-up, contamination, biochemical problems, incomplete sequencing, misdocumentation and insufficient reference to previously published data. To assess data quality in case studies of mitochondrial diseases, it is recommended to compare any mtDNA sequence under consideration to their phylogenetically closest lineages available in the Web. The median network method has proven useful for visualizing potential problems with the data. We contrast some early reports of complete mtDNA sequences to more recent total mtDNA sequencing efforts in studies of various mitochondrial diseases. We conclude that the quality of complete mtDNA sequences generated in the medical field in the past few years is somewhat unsatisfactory and may even fall behind that of pioneer manual sequencing in the early nineties. Our study provides a paradigm for an a posteriori evaluation of sequence quality and for detection of potential problems with inferring a pathogenic status of a particular mutation.

mtDNAmanager: a Web-based tool for the management and quality analysis of mitochondrial DNA control-region sequences

Background

For the past few years, scientific controversy has surrounded the large number of errors in forensic and literature mitochondrial DNA (mtDNA) data. However, recent research has shown that using mtDNA phylogeny and referring to known mtDNA haplotypes can be useful for checking the quality of sequence data.

Results

We developed a Web-based bioinformatics resource "mtDNAmanager" that offers a convenient interface supporting the management and quality analysis of mtDNA sequence data. The mtDNAmanager performs computations on mtDNA control-region sequences to estimate the most-probable mtDNA haplogroups and retrieves similar sequences from a selected database. By the phased designation of the most-probable haplogroups (both expected and estimated haplogroups), mtDNAmanager enables users to systematically detect errors whilst allowing for confirmation of the presence of clear key diagnostic mutations and accompanying mutations. The query tools of mtDNAmanager also facilitate database screening with two options of "match" and "include the queried nucleotide polymorphism". In addition, mtDNAmanager provides Web interfaces for users to manage and analyse their own data in batch mode.

Conclusion

The mtDNAmanager will provide systematic routines for mtDNA sequence data management and analysis via easily accessible Web interfaces, and thus should be very useful for population, medical and forensic studies that employ mtDNA analysis. mtDNAmanager can be accessed at .

The mtDNA ancestry of admixed Colombian populations

A total of 185 individuals from Colombia were sequenced for the first hypervariable region (HVS-I) of the mitochondrial DNA (mtDNA) genome, and a subset of these individuals were additionally genotyped for the second hypervariable segment (HVS-II). These individuals were collected according to their "self-reported ethnicity" in Colombia, comprising "Mestizos," "Mulatos," and "Afro-Colombians." We used databases containing more than 4,300 Native American lineages, 6,800 Africans, and 15,600 Europeans for population comparisons and phylogeographic inferences. We observe that Mulatos and Afro-Colombians have a dominant African mtDNA component, whereas Mestizos carry predominantly Native American haplotypes. All the populations analyzed have high diversity indices and there are no signatures of dramatic genetic drift episodes. Central and South America are the main candidate source populations of the Colombian Native American lineages, whereas west-central, southwest, and southeast Africa are the main original mtDNA sources for the African Colombian mtDNAs. We found that our results differ from those obtained in other studies for the same "population groups" in terms of haplogroup frequencies. This observation leads us to conclude that (i) self-reported ancestry is not a reliable proxy to indicate an individual's "ethnicity" in Colombia, (ii) our results do not support the use of outmoded race descriptions (Mestizos, Mulatos, etc.) mainly because these labels do not correspond to any genetically homogeneous population group, and (iii) studies relying on these terms to describe the population group of the individual, which then treat them as genetically homogeneous, carry a high risk of type I error (false positives) in medical studies in this country and of misinterpretation of the frequency of observed variation in forensic casework.

Distilling artificial recombinants from large sets of complete mtDNA genomes

Background

Large-scale genome sequencing poses enormous problems to the logistics of laboratory work and data handling. When numerous fragments of different genomes are PCR amplified and sequenced in a laboratory, there is a high immanent risk of sample confusion. For genetic markers, such as mitochondrial DNA (mtDNA), which are free of natural recombination, single instances of sample mix-up involving different branches of the mtDNA phylogeny would give rise to reticulate patterns and should therefore be detectable.

Methodology/Principal Findings

We have developed a strategy for comparing new complete mtDNA genomes, one by one, to a current skeleton of the worldwide mtDNA phylogeny. The mutations distinguishing the reference sequence from a putative recombinant sequence can then be allocated to two or more different branches of this phylogenetic skeleton. Thus, one would search for two (or three) near-matches in the total mtDNA database that together best explain the variation seen in the recombinants. The evolutionary pathway from the mtDNA tree connecting this pair together with the recombinant then generate a grid-like median network, from which one can read off the exchanged segments.

Conclusions

We have applied this procedure to a large collection of complete human mtDNA sequences, where several recombinants could be distilled by our method. All these recombinant sequences were subsequently corrected by de novo experiments – fully concordant with the predictions from our data-analytical approach.

Mitochondrial DNA control region polymorphism in the population of Alagoas state, north-eastern Brazil

The sequences of the two hypervariable (HV) segments of the mitochondrial DNA (mtDNA) control region were determined in 167 randomly selected, unrelated individuals living in the state of Alagoas, north-eastern Brazil. One hundred and forty-five different haplotypes, associated with 139 variable positions, were determined. More than 95% of the mtDNA sequences could be allocated to specific mtDNA haplogroups according to the mutational motifs. Length heteroplasmy in the C-stretch HV1 and HV2 regions was observed in 22 and 11%, respectively, of the population sample. The genetic diversity was estimated to be 0.9975 and the probability of two random individuals presenting identical mtDNA haplotypes was 0.0084. The most frequent haplotype was shared by six individuals. All sequences showed high-quality values and phantom mutations were not detected. The diversity revealed in the mitochondrial control region indicates the importance of this locus for forensic casework and population studies within Alagoas, Brazil.

Minisequencing mitochondrial DNA pathogenic mutations

BACKGROUND

There are a number of well-known mutations responsible of common mitochondrial DNA (mtDNA) diseases. In order to overcome technical problems related to the analysis of complete mtDNA genomes, a variety of different techniques have been proposed that allow the screening of coding region pathogenic mutations.

METHODS

We here propose a minisequencing assay for the analysis of mtDNA mutations. In a single reaction, we interrogate a total of 25 pathogenic mutations distributed all around the whole mtDNA genome in a sample of patients suspected for mtDNA disease.

RESULTS

We have detected 11 causal homoplasmic mutations in patients suspected for Leber disease, which were further confirmed by standard automatic sequencing. Mutations m.11778G>A and m.14484T>C occur at higher frequency than expected by change in the Galician (northwest Spain) patients carrying haplogroup J lineages (Fisher's Exact test, P-value < 0.01). The assay performs well in mixture experiments of wild:mutant DNAs that emulate heteroplasmic conditions in mtDNA diseases.

CONCLUSION

We here developed a minisequencing genotyping method for the screening of the most common pathogenic mtDNA mutations which is simple, fast, and low-cost. The technique is robust and reproducible and can easily be implemented in standard clinical laboratories.

The Genographic Project public participation mitochondrial DNA database

The Genographic Project is studying the genetic signatures of ancient human migrations and creating an open-source research database. It allows members of the public to participate in a real-time anthropological genetics study by submitting personal samples for analysis and donating the genetic results to the database. We report our experience from the first 18 months of public participation in the Genographic Project, during which we have created the largest standardized human mitochondrial DNA (mtDNA) database ever collected, comprising 78,590 genotypes. Here, we detail our genotyping and quality assurance protocols including direct sequencing of the mtDNA HVS-I, genotyping of 22 coding-region SNPs, and a series of computational quality checks based on phylogenetic principles. This database is very informative with respect to mtDNA phylogeny and mutational dynamics, and its size allows us to develop a nearest neighbor-based methodology for mtDNA haplogroup prediction based on HVS-I motifs that is superior to classic rule-based approaches. We make available to the scientific community and general public two new resources: a periodically updated database comprising all data donated by participants, and the nearest neighbor haplogroup prediction tool.

2006 GEP-ISFG collaborative exercise on mtDNA: reflections about interpretation, artefacts, and DNA mixtures

We report the results of the seventh edition of the GEP-ISFG mitochondrial DNA (mtDNA) collaborative exercise. The samples submitted to the participant laboratories were blood stains from a maternity case and simulated forensic samples, including a case of mixture. The success rate for the blood stains was moderate ( approximately 77%); even though four inexperienced laboratories concentrated about one-third of the total errors. A similar success was obtained for the analysis of mixed samples (78.8% for a hair-saliva mixture and 69.2% for a saliva-saliva mixture). Two laboratories also dissected the haplotypes contributing to the saliva-saliva mixture. Most of the errors were due to reading problems and misinterpretation of electropherograms, demonstrating once more that the lack of a solid devised experimental approach is the main cause of error in mtDNA testing.

Italian mitochondrial DNA database: results of a collaborative exercise and proficiency testing

This work is a review of a collaborative exercise on mtDNA analysis undertaken by the Italian working group (Ge.F.I.). A total of 593 samples from 11 forensic genetic laboratories were subjected to hypervariable region (HVS-I/HVS-II) sequence analysis. The raw lane data were sent to MtDNA Population Database (EMPOP) for an independent evaluation. For the inclusion of data for the Italian database, quality assurance procedures were applied to the control region profiles. Only eight laboratories with a final population sample of 395 subjects passed the quality conformance test. Control region haplogroup (hg) assignments were confirmed by restriction fragment length polymorphism (RFLP) typing of the most common European hg-diagnostic sites. A total of 306 unique haplotypes derived from the combined analysis of control and coding region polymorphisms were found; the most common haplotype--CRS, 263, 309.1C, 315.1C/ not7025 AluI--was shared by 20 subjects. The majority of mtDNAs detected in the Italian population fell into the most common west Eurasian hgs: R0a (0.76%), HV (4.81%), H (38.99%), HV0 (3.55%), J (7.85%), T (13.42%), U (11.65%), K (10.13%), I (1.52%), X (2.78%), and W (1.01%).

Specificity of mtDNA-directed PCR-influence of NUclear MTDNA insertion (NUMT) contamination in routine samples and techniques

Nuclear mitochondrial insertions (NUMTs) are sequences homologous to mtDNA, which are present throughout the human nuclear genome. The possibility that these sequences may be accidentally amplified in reactions directed to mtDNA has been raised and evaluated by different groups and by different means. Despite that, data is still missing on the specificity of PCRs in routine procedures in what concerns contamination with nuclear mtDNA insertions (NUMTs). In this work, we performed PCR sequencing reactions with primers directed either to mitochondrial or to NUMT DNA with different annealing temperatures and in different tissues. We observed that (a) contamination with NUMTs depends on the sample and tissue, and (b) employing routine techniques, there is no risk of co-amplification. Only when mtDNA is almost completely removed from the samples does the number of NUMT copies exceed mitochondrial sequences, i.e., only in samples with virtually no mtDNA, such as those resulting from preferential semen lysis, is there a risk of accidental amplification of NUMTs. We suggest that to evaluate a possible co-amplification of NUMT DNA, it is more relevant to take into account sample processing and original tissue of the samples, and consequently the relative proportions of NUMT and mtDNA, rather than the presence of NUMTs by itself, irrespectively of its proportion.

Utility of haplogroup determination for forensic mtDNA analysis in the Japanese population

Sequence analysis of the hypervariable regions (HVRs) of mitochondrial DNA (mtDNA) are routinely performed in forensic casework, however, there are still issues to be resolved, such as the existence of multiple errors in published databases or the limitations of individual discrimination in certain populations. Here, we analyzed the coding region of mtDNA in detail by examining 36 haplogroup (HG)-defining single nucleotide polymorphisms (SNPs) using amplified product-length polymorphisms (APLP) method in conjunction with sequence analysis of HVR1 and HVR2 to establish a methodology for forensically reliable and practical mtDNA testing. The mtDNAs from 217 unrelated Japanese were examined and could be classified into 27 haplogroups. By combining the data of the coding region with those of HVRs, genetic diversity was slightly increased from 0.9817 to 0.9888 for HVR1/HG and from 0.9967 to 0.9970 for HVR1/HVR2/HG, as compared to the results of HVRs only. Moreover, in most cases, reliability of the HVR data could be confirmed by haplogroup motif analysis. Our mtDNA profiling method can provide reliable data in a time and cost-saving way due to the rapid and economical nature of APLP analysis.

Mitochondrial diversity within modern human populations

With the recent increase in the available number of high-quality, full-length mitochondrial sequences, it is now possible to construct and analyze a comprehensive human mitochondrial consensus sequence. Using a data set of 827 carefully selected sequences, it is shown that modern humans contain extremely low levels of divergence from the mitochondrial consensus sequence, differing by a mere 21.6 nt sites on average. Fully 84.1% of the mitochondrial genome was found to be invariant and ‘private’ mutations accounted for 43.8% of the variable sites. Ninety eight percent of the variant sites had a primary nucleotide with an allele frequency of 0.90 or greater. Interestingly, the few truly ambiguous nucleotide sites could all be reliably assigned to either a purine or pyrimidine ancestral state. A comparison of this consensus sequence to several ancestral sequences derived from phylogenetic studies reveals a great deal of similarity, where, as expected, the most phylogenetically informative nucleotides in the ancestral studies tended to be the most variable nucleotides in the consensus. Allowing for this fact, the consensus approach provides variation data on the positions that do not contribute to phylogenetic reconstructions, and these data provide a baseline for measuring human mitochondrial variation in populations worldwide.

Development and expansion of high-quality control region databases to improve forensic mtDNA evidence interpretation

In an effort to increase the quantity, breadth and availability of mtDNA databases suitable for forensic comparisons, we have developed a high-throughput process to generate approximately 5000 control region sequences per year from regional US populations, global populations from which the current US population is derived and global populations currently under-represented in available forensic databases. The system utilizes robotic instrumentation for all laboratory steps from pre-extraction through sequence detection, and a rigorous eight-step, multi-laboratory data review process with entirely electronic data transfer. Over the past 3 years, nearly 10,000 control region sequences have been generated using this approach. These data are being made publicly available and should further address the need for consistent, high-quality mtDNA databases for forensic testing.