Detecting errors in mtDNA data by phylogenetic analysis

Guidelines for the Analysis of DNA Barcoding/Metabarcoding Sequencing Data and Interpretation of Publicly Available Databases

This chapter describes procedures for the use of DNA sequence data to obtain and compare taxonomic identification using the public databases GenBank and Barcode of Life Data System (BOLD). The chapter begins by describing procedures used to prepare quality sequences for uploading into GenBank and BOLD. Next, steps used to query the DNA sequences against the public databases are described using GenBank BLAST and BOLD identification engines. Interpretation guidelines for the taxonomic identification assignments are presented. Finally, a procedure for evaluating the accuracy and reliability of sequences from GenBank and BOLD is provided.

DNA Barcoding and Metabarcoding Protocols for Species Identification

The article presents the several steps to be performed on a plant, fungal, insect, or soil sample to obtain DNA sequences for DNA barcode analysis. The chapter begins with a description of sample preparation including procedures for cleaning and proceeds to DNA extraction with methods adapted for the specific type of sample. Next, DNA quantification is described so the proper amount is used for the amplification of the selected barcode regions. Information is provided for reaction mixes and amplification conditions for several referenced barcode primer pairs tuned for the individual sample of interest. This is followed by a description of procedures to access the success of amplification, cleanup, and quantification of the product ready for either Sanger sequencing or library preparation for massive parallel sequencing (MPS). Finally, procedures are provided for Sanger sequencing, library preparation, and MPS sequencing. The chapter provides several references of barcode regions for different sample types.

Effect of fuchsin fixation of pollen on DNA barcode recovery

Pollen grains attached to insects are a valuable source of ecological information which can be used to reconstruct visitation networks. Morphological pollen identification relies on light microscopy with pollen usually stained and mounted in fuchsin jelly, which is also used to remove pollen from the bodies of insects. Pollen embedded in fuchsin jelly could potentially be used for DNA barcoding and metabarcoding (large-scale taxonomic identification of complex mixed samples) and thus provide additional information for pollination networks. In this study, we determine whether fuchsin-embedded pollen can be used for downstream molecular applications. We evaluate the quality of plant barcode (ITS) sequences amplified from DNA extracted from both fresh (untreated) pollen, and pollen which had been embedded in fuchsin jelly. We show that the addition of fuchsin to DNA extraction does not impact DNA barcode sequence quality during short-term storage. DNA extractions from both untreated and fuchsin-treated pollen produced reliable barcode sequences of high quality. Our findings suggest that pollen which has been collected, stained, and embedded in fuchsin jelly for preliminary microscopy work can be used within several days for downstream genetic analysis, though the quality of DNA from pollen stored in fuchsin jelly for extended periods is yet to be established.

Sequencing errors in Native American mitogenomes: impact on clade definitions, haplogroup assignation, and beyond

Available evidence allows the interpretation that some cases of absence of otherwise expected variation, based on phylogenetic expectations in mitogenomes of Native American origin, are due to artificial recombination rather than to homoplasy, while other more complex scenarios involving combination of original Cambridge Reference Sequence mistakes plus incomplete or incorrect scoring of variation are also showed. Several instances of mismatched control and coding regions as well as partially duplicated HV2 are observed in Peruvians, while intra-haplogroup chimaeras of different D1 subhaplogroups are referred to in Mexican Native Americans. A revised definition for haplogroup B2h is proposed, and preventive quality control measures are suggested.

VLF: An R package for the analysis of very low frequency variants in DNA sequences

Here, we introduce VLF, an R package to determine the distribution of very low frequency variants (VLFs) in nucleotide and amino acid sequences for the analysis of errors in DNA sequence records. The package allows users to assess VLFs in aligned and trimmed protein-coding sequences by automatically calculating the frequency of nucleotides or amino acids in each sequence position and outputting those that occur under a user-specified frequency (default of p = 0.001). These results can then be used to explore fundamental population genetic and phylogeographic patterns, mechanisms and processes at the microevolutionary level, such as nucleotide and amino acid sequence conservation. Our package extends earlier work pertaining to an implementation of VLF analysis in Microsoft Excel, which was found to be both computationally slow and error prone. We compare those results to our own herein. Results between the two implementations are found to be highly consistent for a large DNA barcode dataset of bird species. Differences in results are readily explained by both manual human error and inadequate Linnean taxonomy (specifically, species synonymy). Here, VLF is also applied to a subset of avian barcodes to assess the extent of biological artifacts at the species level for Canada goose (Branta canadensis), as well as within a large dataset of DNA barcodes for fishes of forensic and regulatory importance. The novelty of VLF and its benefit over the previous implementation include its high level of automation, speed, scalability and ease-of-use, each desirable characteristics which will be extremely valuable as more sequence data are rapidly accumulated in popular reference databases, such as BOLD and GenBank.

The Value of Whole-Genome Sequencing for Mitochondrial DNA Population Studies: Strategies and Criteria for Extracting High-Quality Mitogenome Haplotypes

Whole-genome sequencing (WGS) data present a readily available resource for mitochondrial genome (mitogenome) haplotypes that can be utilized for genetics research including population studies. However, the reconstruction of the mitogenome is complicated by nuclear mitochondrial DNA (mtDNA) segments (NUMTs) that co-align with the mtDNA sequences and mimic authentic heteroplasmy. Two minimum variant detection thresholds, 5% and 10%, were assessed for the ability to produce authentic mitogenome haplotypes from a previously generated WGS dataset. Variants associated with NUMTs were detected in the mtDNA alignments for 91 of 917 (~8%) Swedish samples when the 5% frequency threshold was applied. The 413 observed NUMT variants were predominantly detected in two regions (nps 12,612–13,105 and 16,390–16,527), which were consistent with previously documented NUMTs. The number of NUMT variants was reduced by ~97% (400) using a 10% frequency threshold. Furthermore, the 5% frequency data were inconsistent with a platinum-quality mitogenome dataset with respect to observed heteroplasmy. These analyses illustrate that a 10% variant detection threshold may be necessary to ensure the generation of reliable mitogenome haplotypes from WGS data resources.

Forensically relevant phylogeographic evaluation of mitogenome variation in the Basque Country

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

Validation of NGS for mitochondrial DNA casework at the FBI Laboratory

As a first step towards integrating next generation sequencing (NGS) technology into the FBI Laboratory's operational casework, the PowerSeq™ CRM Nested System, an NGS-based mitochondrial DNA (mtDNA) control region assay, was developmentally and internally validated. The validation studies were conducted in accordance with the Scientific Working Group on DNA Analysis Methods (SWGDAM) Validation Guidelines for Forensic DNA Analysis Methods, and the FBI's Quality Assurance Standards (QAS) for Forensic DNA Testing Laboratories. The assay was shown to be highly reproducible, with variant frequencies across intra and inter-run replicates of the same sample differing, on average, by just 0.3% for substitutions and point heteroplasmies and 1.5% for insertions and deletions. The assay was also shown to be extremely sensitive, yielding complete control region sequence data from as few as 2000 copies of mtDNA. This is a more than 20-fold increase in sensitivity when compared to the FBI Laboratory's current Sanger sequencing-based protocols and, based on mtDNA quantitation values of samples routinely encountered in mtDNA casework, suggests that the percentage of questioned samples from which full control region data can be recovered will increase from our current 20% to approximately 90% success with NGS technology. In addition, the assay requires on average only 30% of the extract volume typically required to develop control region profiles from degraded samples via Sanger sequencing. Overall, these studies establish the reliability of the PowerSeq™ CRM Nested System for accurate mtDNA control region typing and can serve as a model for laboratories seeking to validate NGS protocols for forensic mtDNA analysis.

Analysis of the mitochondrial CYTB gene sequence in human populations of northeastern Bosnia

This study offers the first report on variation sequence of the mitochondrial cytochrome b (MTCYTB) gene in populations from Bosnia (northeastern Bosnia). This study was designed on the analysis of the genetic diversity of two populations of different cultural-anthropological and genetic origin, Roma population and native/non-Roma population. The main aim of our study was to estimate the usefulness of the CYTB sequence in the analysis of genetic categorization of different populations and intergroup diversity, as well as to provide some additional information on haplogroup-associated polymorphisms within the CYTB region in defining haplogroup status. Estimation of the genetic diversity was done using intra and intergroup genetic indices. The population-specific polymorphisms have been found in both categories of the populations. The results of the analysis of genetic differentiation show significant pairwise Fst differences between the Romani and native populations. Also, registered significant genetic differentiation is illustrated on the level of genetic variation between subpopulations of the Roma and non-Roma origin. The important result in our study is the confirmation of the significance of the triad of polymorphisms T14783C-G15043A-G15301A, indicating the influence of Asian component of the maternal gene pool on the genetic structure of the studied population of the Roma. Our data show that the haplogroup polymorphisms exist in the CYTB region and can provide useful information on the haplogroups that were defined only by the control region of the mtDNA. The results of this study indicate the region of CYTB gene can be a benefit in providing some additional information in the analysis of genetic structure of human populations and can be additionally applied in population studies.

Next generation database search algorithm for forensic mitogenome analyses

Mitochondrial DNA (mtDNA) variation is being reported relative to the corrected version of the first sequenced human mitochondrial genome. A review of the existing literature across disciplines that employ mtDNA demonstrates that insertions and deletions are not reported in a standardized way. This may lead to false exclusions of identical sequences, unidentified matches in missing persons mtDNA databases, biased mtDNA database frequency estimates and overestimation of the genetic evidence. Seven years ago we introduced alignment-free database search software (SAM) and implemented it into the mtDNA database EMPOP (https://empop.online) to produce reliable and conservative frequency estimates that are required in the forensic context. However, ambiguity remained in how laboratories have been reporting mitotypes, as often more than one single alignment of a given mtDNA sequence was feasible. In order to overcome this limitation we here describe a concept and provide software for producing stable, harmonized phylogenetic alignment of mtDNA sequences for database searches. The new software SAM 2 will be made available via EMPOP and provide the user with the already established conservative frequency estimates. In addition, SAM 2 offers the rCRS-coded haplotype of a given mtDNA sequence following the established and widely accepted phylogenetic alignment. This provides the user with feedback on how mitotypes are stored in EMPOP and how they should be reported in order to harmonize nomenclature. Finally, this approach does not only permit reliable mtDNA nomenclature in forensics but invites related disciplines to take advantage of a standardized way of reporting mtDNA variation, thus closing the ranks between different genetic fields and supporting dialogue and collaboration between mtDNA scholars from various disciplines.

A Comparison and Integration of MiSeq and MinION Platforms for Sequencing Single Source and Mixed Mitochondrial Genomes

Single source and multiple donor (mixed) samples of human mitochondrial DNA were analyzed and compared using the MinION and the MiSeq platforms. A generalized variant detection strategy was employed to provide a cursory framework for evaluating the reliability and accuracy of mitochondrial sequences produced by the MinION. The feasibility of long-read phasing was investigated to establish its efficacy in quantitatively distinguishing and deconvolving individuals in a mixture. Finally, a proof-of-concept was demonstrated by integrating both platforms in a hybrid assembly that leverages solely mixture data to accurately reconstruct full mitochondrial genomes.

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.

Mitochondrial DNA mutations in single human blood cells

Determination mitochondrial DNA (mtDNA) sequences from extremely small amounts of DNA extracted from tissue of limited amounts and/or degraded samples is frequently employed in medical, forensic, and anthropologic studies. Polymerase chain reaction (PCR) amplification followed by DNA cloning is a routine method, especially to examine heteroplasmy of mtDNA mutations. In this review, we compare the mtDNA mutation patterns detected by three different sequencing strategies. Cloning and sequencing methods that are based on PCR amplification of DNA extracted from either single cells or pooled cells yield a high frequency of mutations, partly due to the artifacts introduced by PCR and/or the DNA cloning process. Direct sequencing of PCR product which has been amplified from DNA in individual cells is able to detect the low levels of mtDNA mutations present within a cell. We further summarize the findings in our recent studies that utilized this single cell method to assay mtDNA mutation patterns in different human blood cells. Our data show that many somatic mutations observed in the end-stage differentiated cells are found in hematopoietic stem cells (HSCs) and progenitors within the CD34(+) cell compartment. Accumulation of mtDNA variations in the individual CD34+ cells is affected by both aging and family genetic background. Granulocytes harbor higher numbers of mutations compared with the other cells, such as CD34(+) cells and lymphocytes. Serial assessment of mtDNA mutations in a population of single CD34(+) cells obtained from the same donor over time suggests stability of some somatic mutations. CD34(+) cell clones from a donor marked by specific mtDNA somatic mutations can be found in the recipient after transplantation. The significance of these findings is discussed in terms of the lineage tracing of HSCs, aging effect on accumulation of mtDNA mutations and the usage of mtDNA sequence in forensic identification.

Mitochondrial DNA heteroplasmy in the emerging field of massively parallel sequencing

Long an important and useful tool in forensic genetic investigations, mitochondrial DNA (mtDNA) typing continues to mature. Research in the last few years has demonstrated both that data from the entire molecule will have practical benefits in forensic DNA casework, and that massively parallel sequencing (MPS) methods will make full mitochondrial genome (mtGenome) sequencing of forensic specimens feasible and cost-effective. A spate of recent studies has employed these new technologies to assess intraindividual mtDNA variation. However, in several instances, contamination and other sources of mixed mtDNA data have been erroneously identified as heteroplasmy. Well vetted mtGenome datasets based on both Sanger and MPS sequences have found authentic point heteroplasmy in approximately 25% of individuals when minor component detection thresholds are in the range of 10-20%, along with positional distribution patterns in the coding region that differ from patterns of point heteroplasmy in the well-studied control region. A few recent studies that examined very low-level heteroplasmy are concordant with these observations when the data are examined at a common level of resolution. In this review we provide an overview of considerations related to the use of MPS technologies to detect mtDNA heteroplasmy. In addition, we examine published reports on point heteroplasmy to characterize features of the data that will assist in the evaluation of future mtGenome data developed by any typing method.

How much information is needed to infer reticulate evolutionary histories?

Phylogenetic networks are a generalization of evolutionary trees and are an important tool for analyzing reticulate evolutionary histories. Recently, there has been great interest in developing new methods to construct rooted phylogenetic networks, that is, networks whose internal vertices correspond to hypothetical ancestors, whose leaves correspond to sampled taxa, and in which vertices with more than one parent correspond to taxa formed by reticulate evolutionary events such as recombination or hybridization. Several methods for constructing evolutionary trees use the strategy of building up a tree from simpler building blocks (such as triplets or clusters), and so it is natural to look for ways to construct networks from smaller networks. In this article, we shall demonstrate a fundamental issue with this approach. Namely, we show that even if we are given all of the subnetworks induced on all proper subsets of the leaves of some rooted phylogenetic network, we still do not have all of the information required to completely determine that network. This implies that even if all of the building blocks for some reticulate evolutionary history were to be taken as the input for any given network building method, the method might still output an incorrect history. We also discuss some potential consequences of this result for constructing phylogenetic networks.

Improved visibility of character conflicts in quasi-median networks with the EMPOP NETWORK software

Aim

To provide a valuable tool for graphical representation of mitochondrial DNA (mtDNA) data that enables visual emphasis on complex substructures within the network to highlight possible ambiguities and errors.

Method

We applied the new NETWORK graphical user interface, available via EMPOP (European DNA Profiling Group Mitochondrial DNA Population Database; www.empop.org) by means of two mtDNA data sets that were submitted for quality control.

Results

The quasi-median network torsi of the two data sets resulted in complex reticulations, suggesting ambiguous data. To check the corresponding raw data, accountable nodes and connecting branches of the network could be identified by highlighting induced subgraphs with concurrent dimming of their complements. This is achieved by accentuating the relevant substructures in the network: mouse clicking on a node displays a list of all mtDNA haplotypes included in that node; the selection of a branch specifies the mutation(s) connecting two nodes. It is indicated to evaluate these mutations by means of the raw data.

Conclusion

Inspection of the raw data confirmed the presence of phantom mutations due to suboptimal electrophoresis conditions and data misinterpretation. The network software proved to be a powerful tool to highlight problematic data and guide quality control of mtDNA data tables.

mitoSAVE: mitochondrial sequence analysis of variants in Excel

The mitochondrial genome (mtGenome) contains genetic information amenable to numerous applications such as medical research, population and evolutionary studies, and human identity testing. However, inconsistent nomenclature assignment makes haplotype comparison difficult and can lead to false exclusion of potentially useful profiles. Massively Parallel Sequencing (MPS) is a platform for sequencing large datasets and potentially whole populations with relative ease. However, the data generated are not easily parsed and interpreted. With this in mind, mitoSAVE has been developed to enable fast conversion of Variant Call Format (VCF) files. mitoSAVE is an Excel-based workbook that converts data within the VCF into mtDNA haplotypes using phylogenetically-established nomenclature as well as rule-based alignments consistent with current forensic standards. mitoSAVE is formatted for human mitochondrial genome; however, it can easily be adapted to support other reasonably small genomes.

Forensic informativity of ~3000 bp of coding sequence of domestic dog mtDNA

The discriminatory power of the noncoding control region (CR) of domestic dog mitochondrial DNA alone is relatively low. The extent to which the discriminatory power could be increased by analyzing additional highly variable coding regions of the mitochondrial genome (mtGenome) was therefore investigated. Genetic variability across the mtGenome was evaluated by phylogenetic analysis, and the three most variable ~1 kb coding regions identified. We then sampled 100 Swedish dogs to represent breeds in accordance with their frequency in the Swedish population. A previously published dataset of 59 dog mtGenomes collected in the United States was also analyzed. Inclusion of the three coding regions increased the exclusion capacity considerably for the Swedish sample, from 0.920 for the CR alone to 0.964 for all four regions. The number of mtDNA types among all 159 dogs increased from 41 to 72, the four most frequent CR haplotypes being resolved into 22 different haplotypes.

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.

Tracing Arab-Islamic inheritance in Madagascar: study of the Y-chromosome and mitochondrial DNA in the Antemoro

Madagascar is located at the crossroads of the Asian and African worlds and is therefore of particular interest for studies on human population migration. Within the large human diversity of the Great Island, we focused our study on a particular ethnic group, the Antemoro. Their culture presents an important Arab-Islamic influence, but the question of an Arab biological inheritance remains unresolved. We analyzed paternal (n=129) and maternal (n=135) lineages of this ethnic group. Although the majority of Antemoro genetic ancestry comes from sub-Saharan African and Southeast Asian gene pools, we observed in their paternal lineages two specific haplogroups (J1 and T1) linked to Middle Eastern origins. This inheritance was restricted to some Antemoro sub-groups. Statistical analyses tended to confirm significant Middle Eastern genetic contribution. This study gives a new perspective to the large human genetic diversity in Madagascar.

An update to MitoTool: using a new scoring system for faster mtDNA haplogroup determination

The determination of human mitochondrial DNA (mtDNA) haplogroups is not only crucial in anthropological and forensic studies, but is also helpful in the medical field to prevent establishment of wrong disease associations. In recent years, high-throughput technologies and the huge amounts of data they create, as well as the regular updates to the mtDNA phylogenetic tree, mean that there is a need for an automated approach which can make a speedier determination of haplogroups than can be made by using the traditional manual method. Here, we update the MitoTool (www.mitotool.org) by incorporating a novel scoring system for the determination of mtDNA into haplogroups, which has advantages on speed, accuracy and ease of implementation. In order to make the access to MitoTool easier, we also provide a stand-alone version of the program that will run on a local computer and this version is freely available at the MitoTool website.

The mitochondrial genome: a biosensor for early cancer detection?

Mutations in the mitochondrial genome have been reported as biomarkers for the detection of cancer. Hallmarks of cancer development include the accumulation of genetic alterations in the mitochondrial and nuclear genomes. Damage to mitochondria affects energy metabolism, generation of reactive oxygen species, apoptosis, cell growth and other processes that contribute to the neoplastic process. Furthermore, mitochondrial DNA mutations occur frequently in cancer. Little work has been done to link a pathway between mitochondrial mutations and cancer etiology. Volumes of work have been reported on the association of mitochondrial mutations and almost all types of cancer including the use of body fluids for early detection. This review examines the measurement of mitochondrial mutations for the application of detecting human tumor tissue.

Somatic mitochondrial DNA mutations in Chinese patients with osteosarcoma

Somatic mutations in mitochondrial DNA (mtDNA) have been long proposed to drive the pathogenesis and progression of human malignancies. Previous investigations have revealed a high frequency of somatic mutations in the D-loop control region of mtDNA in osteosarcoma. However, little is known with regard to whether or not somatic mutations also occur in the coding regions of mtDNA in osteosarcoma. To test this possibility, in the present study we screened somatic mutations over the full-length mitochondrial genome of 31 osteosarcoma tumour tissue samples, and corresponding peripheral blood samples from the same cohort of patients. We detected a sum of 11 somatic mutations in the mtDNA coding regions in our series. Nine of them were missense or frameshift mutations that have the potential to hamper mitochondrial respiratory function. In combination with our earlier observations on the D-loop fragment, 71.0% (22/31) of patients with osteosarcoma carried at least one somatic mtDNA mutation, and a total of 40 somatic mutations were identified. Amongst them, 29 (72.5%) were located in the D-loop region, two (5%) were in the sequences of the tRNA genes, two (5%) were in the mitochondrial ATP synthase subunit 6 gene and seven (17.5%) occurred in genes encoding components of the mitochondrial respiratory complexes. In addition, somatic mtDNA mutation was not closely associated with the clinicopathological characteristics of osteosarcoma. Together, these findings suggest that somatic mutations are highly prevalent events in both coding and non-coding regions of mtDNA in osteosarcoma. Some missense and frameshift mutations are putatively harmful to proper mitochondrial activity and might play vital roles in osteosarcoma carcinogenesis.

A "Copernican" reassessment of the human mitochondrial DNA tree from its root

Mutational events along the human mtDNA phylogeny are traditionally identified relative to the revised Cambridge Reference Sequence, a contemporary European sequence published in 1981. This historical choice is a continuous source of inconsistencies, misinterpretations, and errors in medical, forensic, and population genetic studies. Here, after having refined the human mtDNA phylogeny to an unprecedented level by adding information from 8,216 modern mitogenomes, we propose switching the reference to a Reconstructed Sapiens Reference Sequence, which was identified by considering all available mitogenomes from Homo neanderthalensis. This "Copernican" reassessment of the human mtDNA tree from its deepest root should resolve previous problems and will have a substantial practical and educational influence on the scientific and public perception of human evolution by clarifying the core principles of common ancestry for extant descendants.

Comparative evolution of S7 Intron 1 and ribosomal internal transcribed spacer in Coilia nasus (Clupeiformes: Engraulidae)

Coilia nasus is widely distributed in the Yangtze River, the coastal waters of China, Korea and the Ariake Sound of Japan. Several ecotypes exist and this provides a useful model for the study of comparative diversity between molecular markers. Here we analyze and compare the nucleotide sequences between single-copy ribosomal protein S7 gene intron 1 (rpS7) and multiple-copy ribosomal internal transcribed spacer 1 (ITS1) in this species to compare the phylogenetic signal of the two nuclear genes. Nucleotide substitutions among the two gene sequences and partial sequence of mitochondrial cytochrome c oxidase subunit I (COI) gene were also analyzed. A total of 115 clones for rpS7 and 122 clones for ITS1 were obtained from 37 specimens. The nucleotide sequence length is 741 to 743 bp for rpS7 and 334 to 348 bp for ITS1. Intra- and inter-specimen variation in rpS7 results from nucleotide substitution, while such variation in ITS1 is mainly due to different numbers of short base repeats. The content of G + C is lower in rpS7 (43.5%) than in ITS1 (68.2%). Our results indicate that the proportion of the sequence variable sites is higher in rpS7 (61) than in ITS1 (23); the informative parsimony of rpS7 is evidently higher than that of ITS1 (26 vs. 2); the overall ratio between transitions and transversions in ITS1 is slightly lower than in rpS7, but remarkably lower than in COI. These results suggest that rpS7 is more suitable than ITS1 as a marker for genetic divergence of this group. Furthermore, gene flow is observed between the different geographic populations of C. nasus from the phylogeny of this species based on rpS7, showing that rpS7 has more evolutionary characteristics for understanding the processes of genomic evolution at the intraspecific level.

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.

Capillary electrophoresis of Big-Dye terminator sequencing reactions for human mtDNA Control Region haplotyping in the identification of human remains

Cycle sequencing reaction with Big-Dye terminators provides the methodology to analyze mtDNA Control Region amplicons by means of capillary electrophoresis. DNA sequencing with ddNTPs or terminators was developed by (1). The progressive automation of the method by combining the use of fluorescent-dye terminators with cycle sequencing has made it possible to increase the sensibility and efficiency of the method and hence has allowed its introduction into the forensic field. PCR-generated mitochondrial DNA products are the templates for sequencing reactions. Different set of primers can be used to generate amplicons with different sizes according to the quality and quantity of the DNA extract providing sequence data for different ranges inside the Control Region.

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.

The landscape of mitochondrial DNA variation in human colorectal cancer on the background of phylogenetic knowledge

Recently, an increasing number of studies indicate that mutations in mitochondrial genome may contribute to cancer development or metastasis. Hence, it is important to determine whether the mitochondrial DNA might be a good, clinically applicable marker of cancer. This review describes hereditary as well as somatic mutations reported in mitochondrial DNA of colorectal cancer cells. We showed here that the entire mitochondrial genome mutational spectra are different in colorectal cancer and non-tumor cells. We also placed the described mutations on the phylogenetic context, which highlighted the recurrent problem of data quality. Therefore, the most important rules for adequately assessing the quality of mitochondrial DNA sequence analysis in cancer have been summarized. As follows from this review, neither the reliable spectrum of mtDNA somatic mutations nor the association between hereditary mutations and colorectal cancer risk have been resolved. This indicates that only high resolution studies on mtDNA variability, followed by a proper data interpretation employing phylogenetic knowledge may finally verify the utility of mtDNA sequence (if any) in clinical practice.

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).

The awakening of an advanced malignant cancer: an insult to the mitochondrial genome

BACKGROUND

In only months-to-years a primary cancer can progress to an advanced phenotype that is metastatic and resistant to clinical treatments. As early as the 1900s, it was discovered that the progression of a cancer to the advanced phenotype is often associated with a shift in the metabolic profile of the disease from a state of respiration to anaerobic fermentation - a phenomenon denoted as the Warburg Effect.

SCOPE OF REVIEW

Reports in the literature strongly suggest that the Warburg Effect is generated as a response to a loss in the integrity of the sequence and/or copy number of the mitochondrial genome content within a cancer.

MAJOR CONCLUSIONS

Multiple studies regarding the progression of cancer indicate that mutation, and/or, a flux in the copy number, of the mitochondrial genome content can support the early development of a cancer, until; the mutational load and/or the reduction-to-depletion of the copy number of the mitochondrial genome content induces the progression of the disease to an advanced phenotype.

GENERAL SIGNIFICANCE

Collectively, evidence has revealed that the human cell has incorporated the mitochondrial genome content into a cellular mechanism that, when pathologically actuated, can de(un)differentiate a cancer from the parental tissue of origin into an autonomous disease that disrupts the hierarchical structure-and-function of the human body. This article is part of a Special Issue entitled: Biochemistry of Mitochondria.

Uncovering the profile of somatic mtDNA mutations in Chinese colorectal cancer patients

In the past decade, a high incidence of somatic mitochondrial DNA (mtDNA) mutations has been observed, mostly based on a fraction of the molecule, in various cancerous tissues; nevertheless, some of them were queried due to problems in data quality. Obviously, without a comprehensive understanding of mtDNA mutational profile in the cancerous tissue of a specific patient, it is unlikely to disclose the genuine relationship between somatic mtDNA mutations and tumorigenesis. To achieve this objective, the most straightforward way is to directly compare the whole mtDNA genome variation among three tissues (namely, cancerous tissue, para-cancerous tissue, and distant normal tissue) from the same patient. Considering the fact that most of the previous studies on the role of mtDNA in colorectal tumor focused merely on the D-loop or partial segment of the molecule, in the current study we have collected three tissues (cancerous, para-cancerous and normal tissues) respectively recruited from 20 patients with colorectal tumor and completely sequenced the mitochondrial genome of each tissue. Our results reveal a relatively lower incidence of somatic mutations in these patients; intriguingly, all somatic mutations are in heteroplasmic status. Surprisingly, the observed somatic mutations are not restricted to cancer tissues, for the para-cancer tissues and distant normal tissues also harbor somatic mtDNA mutations with a lower frequency than cancerous tissues but higher than that observed in the general population. Our results suggest that somatic mtDNA mutations in cancerous tissues could not be simply explained as a consequence of tumorigenesis; meanwhile, the somatic mtDNA mutations in normal tissues might reflect an altered physiological environment in cancer patients.

Inferring ethnicity from mitochondrial DNA sequence

BACKGROUND

The assignment of DNA samples to coarse population groups can be a useful but difficult task. One such example is the inference of coarse ethnic groupings for forensic applications. Ethnicity plays an important role in forensic investigation and can be inferred with the help of genetic markers. Being maternally inherited, of high copy number, and robust persistence in degraded samples, mitochondrial DNA may be useful for inferring coarse ethnicity. In this study, we compare the performance of methods for inferring ethnicity from the sequence of the hypervariable region of the mitochondrial genome.

RESULTS

We present the results of comprehensive experiments conducted on datasets extracted from the mtDNA population database, showing that ethnicity inference based on support vector machines (SVM) achieves an overall accuracy of 80-90%, consistently outperforming nearest neighbor and discriminant analysis methods previously proposed in the literature. We also evaluate methods of handling missing data and characterize the most informative segments of the hypervariable region of the mitochondrial genome.

CONCLUSIONS

Support vector machines can be used to infer coarse ethnicity from a small region of mitochondrial DNA sequence with surprisingly high accuracy. In the presence of missing data, utilizing only the regions common to the training sequences and a test sequence proves to be the best strategy. Given these results, SVM algorithms are likely to also be useful in other DNA sequence classification applications.

A new strategy for the discrimination of mitochondrial DNA haplogroups in Han population

Mitochondrial DNA (mtDNA) haplogroup discrimination is interesting not only for phylogenetic and clinical but also for forensic studies. We discriminated the mtDNA haplogroups of 570 healthy unrelated Han people from Zhejiang Province, Southeast China, by comprehensive analysis mutations of the hypervariable segments-I sequence and diagnostic polymorphisms in mtDNA coding region using real-time polymerase chain reaction (RT-PCR), which was compared with the widely used PCR and restriction fragment length polymorphism (PCR-RFLP) method. The results showed that in superhaplogroup M, haplogroup D was the most common haplotype within this assay to 24.6%, and in the other superhaplotype N, haplogroup B and F were the most common groups. Samples re-identified by PCR-RFLP showed the consistent results that were got with RT-PCR. In conclusion, the RT-PCR strategy appears to be an accurate, reproducible, and sensitive technique for the discrimination of mtDNA haplogroups, especially for mass screenings quickly and economically.

The role of the mitochondrial genome in ageing and carcinogenesis

Mitochondrial DNA mutations and polymorphisms have been the focus of intensive investigations for well over a decade in an attempt to understand how they affect fundamental processes such as cancer and aging. Initial interest in mutations occurring in mitochondrial DNA of cancer cells diminished when most were found to be the same mutations which occurred during the evolution of human mitochondrial haplogroups. However, increasingly correlations are being found between various mitochondrial haplogroups and susceptibility to cancer or diseases in some cases and successful aging in others.

HaploGrep: a fast and reliable algorithm for automatic classification of mitochondrial DNA haplogroups

An ongoing source of controversy in mitochondrial DNA (mtDNA) research is based on the detection of numerous errors in mtDNA profiles that led to erroneous conclusions and false disease associations. Most of these controversies could be avoided if the samples' haplogroup status would be taken into consideration. Knowing the mtDNA haplogroup affiliation is a critical prerequisite for studying mechanisms of human evolution and discovering genes involved in complex diseases, and validating phylogenetic consistency using haplogroup classification is an important step in quality control. However, despite the availability of Phylotree, a regularly updated classification tree of global mtDNA variation, the process of haplogroup classification is still time-consuming and error-prone, as researchers have to manually compare the polymorphisms found in a population sample to those summarized in Phylotree, polymorphism by polymorphism, sample by sample. We present HaploGrep, a fast, reliable and straight-forward algorithm implemented in a Web application to determine the haplogroup affiliation of thousands of mtDNA profiles genotyped for the entire mtDNA or any part of it. HaploGrep uses the latest version of Phylotree and offers an all-in-one solution for quality assessment of mtDNA profiles in clinical genetics, population genetics and forensics. HaploGrep can be accessed freely at http://haplogrep.uibk.ac.at.

The GHEP-EMPOP collaboration on mtDNA population data--A new resource for forensic casework

Mitochondrial DNA (mtDNA) population data for forensic purposes are still scarce for some populations, which may limit the evaluation of forensic evidence especially when the rarity of a haplotype needs to be determined in a database search. In order to improve the collection of mtDNA lineages from the Iberian and South American subcontinents, we here report the results of a collaborative study involving nine laboratories from the Spanish and Portuguese Speaking Working Group of the International Society for Forensic Genetics (GHEP-ISFG) and EMPOP. The individual laboratories contributed population data that were generated throughout the past 10 years, but in the majority of cases have not been made available to the scientific community. A total of 1019 haplotypes from Iberia (Basque Country, 2 general Spanish populations, 2 North and 1 Central Portugal populations), and Latin America (3 populations from São Paulo) were collected, reviewed and harmonized according to defined EMPOP criteria. The majority of data ambiguities that were found during the reviewing process (41 in total) were transcription errors confirming that the documentation process is still the most error-prone stage in reporting mtDNA population data, especially when performed manually. This GHEP–EMPOP collaboration has significantly improved the quality of the individual mtDNA datasets and adds mtDNA population data as valuable resource to the EMPOP database (www.empop.org).