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.     

Mitochondrial sequence analysis for forensic identification using pyrosequencing technology.

Over recent years, requests for mtDNA analysis in the field of forensic medicine have notably increased, and the results of such analyses have proved to be very useful in forensic cases where nuclear DNA analysis cannot be performed. Traditionally, mtDNA has been analyzed by DNA sequencing of the two hypervariable regions, HVI and HVII, in the D-loop. DNA sequence analysis using the conventional Sanger sequencing is very robust but time consuming and labor intensive. By contrast, mtDNA analysis based on the pyrosequencing technology provides fast and accurate results from the human mtDNA present in many types of evidence materials in forensic casework. The assay has been developed to determine polymorphic sites in the mitochondrial D-loop as well as the coding region to further increase the discrimination power of mtDNA analysis. The pyrosequencing technology for analysis of mtDNA polymorphisms has been tested with regard to sensitivity, reproducibility, and success rate when applied to control samples and actual casework materials. The results show that the method is very accurate and sensitive; the results are easily interpreted and provide a high success rate on casework samples. The panel of pyrosequencing reactions for the mtDNA polymorphisms were chosen to result in an optimal discrimination power in relation to the number of bases determined.     

Comparison of the complete mtDNA genome sequences of human cell lines--HL-60 and GM10742A--from individuals with pro-myelocytic leukemia and leber hereditary optic neuropathy, respectively, and the inclusion of HL-60 in the NIST human mitochondrial DNA standard reference material--SRM 2392-I

Forensic and clinical laboratories benefit from DNA standard reference materials (SRMs) that provide the quality control and assurance that their results from sequencing unknown samples are correct. Therefore, the mitochondrial DNA (mtDNA) genome of HL-60, a promyelocytic leukemia cell line, has been completely sequenced by four laboratories and will be available to the forensic and medical communities in the spring of 2003; it will be called National Institute of Standards and Technology (NIST) SRM 2392-I. NIST human mtDNA SRM 2392 will continue to be available and includes the DNA from two apparently healthy individuals. Both SRM 2392 and 2392-I contain all the information (e.g. the sequences of 58 unique primer sets) needed to use these SRMs as positive controls for the amplification and sequencing any DNA. Compared to the templates in SRM 2392, the HL-60 mtDNA in SRM 2392-I has two tRNA differences and more polymorphisms resulting in amino acid changes. Four of these HL-60 mtDNA polymorphisms have been associated with Leber Hereditary Optic Neuropathy (LHON), one as an intermediate mutation and three as secondary mutations. The mtDNA from a cell line (GM10742A) from an individual with LHON was also completely sequenced for comparison and contained some of the same LHON mutations. The combination of these particular LHON associated mutations is also found in phylogenetic haplogroup J and its subset, J2, and may only be indicative that HL-60 belongs to haplogroup J, one of nine haplogroups that characterize Caucasian individuals of European descent or may mean that haplogroup J is more prone to LHON. Both these mtDNA SRMs will provide enhanced quality control in forensic identification, medical diagnosis, and single nucleotide polymorphism detection.     

Patterns of linkage disequilibrium in mitochondrial DNA of 16 ruminant populations

Mitochondrial DNA (mtDNA) is a widely employed molecular tool in phylogeography, in the inference of human evolutionary history, in dating the domestication of livestock and in forensic science. In humans and other vertebrates the popularity of mtDNA can be partially attributed to an assumption of strict maternal inheritance, such that there is no recombination between mitochondrial lineages. The recent demonstration that linkage disequilibrium (LD) declines as a function of distance between polymorphic sites in hominid mitochondrial genomes has been interpreted as evidence of recombination between mtDNA haplotypes, and hence nonclonal inheritance. However, critics of mtDNA recombination have suggested that this association is an artefact of an inappropriate measure of LD or of sequencing error, and subsequent studies of other populations have failed to replicate the initial finding. Here we report the analysis of 16 ruminant populations and present evidence that LD significantly declines with distance in five of them. A meta-analysis of the data indicates a nonsignificant trend of LD declining with distance. Most of the earlier criticisms of patterns between LD and distance in hominid mtDNA are not applicable to this data set. Our results suggest that either ruminant mtDNA is not strictly clonal or that compensatory selection has influenced patterns of variation at closely linked sites within the mitochondrial control region. The potential impact of these processes should be considered when using mtDNA as a tool in vertebrate population genetic, phylogenetic and forensic studies.     

MitoLSDB: A Comprehensive Resource to Study Genotype to Phenotype Correlations in Human Mitochondrial DNA Variations

Human mitochondrial DNA (mtDNA) encodes a set of 37 genes which are essential structural and functional components of the electron transport chain. Variations in these genes have been implicated in a broad spectrum of diseases and are extensively reported in literature and various databases. In this study, we describe MitoLSDB, an integrated platform to catalogue disease association studies on mtDNA (http://mitolsdb.igib.res.in). The main goal of MitoLSDB is to provide a central platform for direct submissions of novel variants that can be curated by the Mitochondrial Research Community. MitoLSDB provides access to standardized and annotated data from literature and databases encompassing information from 5231 individuals, 675 populations and 27 phenotypes. This platform is developed using the Leiden Open (source) Variation Database (LOVD) software. MitoLSDB houses information on all 37 genes in each population amounting to 132397 variants, 5147 unique variants. For each variant its genomic location as per the Revised Cambridge Reference Sequence, codon and amino acid change for variations in protein-coding regions, frequency, disease/phenotype, population, reference and remarks are also listed. MitoLSDB curators have also reported errors documented in literature which includes 94 phantom mutations, 10 NUMTs, six documentation errors and one artefactual recombination. MitoLSDB is the largest repository of mtDNA variants systematically standardized and presented using the LOVD platform. We believe that this is a good starting resource to curate mtDNA variants and will facilitate direct submissions enhancing data coverage, annotation in context of pathogenesis and quality control by ensuring non-redundancy in reporting novel disease associated variants.     

Capturing single-copy nuclear genes,organellar genomes,and nuclear ribosomal DNA from deep genome skimming data for plant phylogenetics: A case study in Vitaceae

With the decreasing cost and availability of many newly developed bioinformatics pipelines, next-generation sequencing (NGS) has revolutionized plant systematics in recent years. Genome skimming has been widely used to obtain high-copy fractions of the genomes, including plastomes, mitochondrial DNA (mtDNA), and nuclear ribosomal DNA (nrDNA). In this study, through simulations, we evaluated the optimal (minimum) sequencing depth and performance for recovering single-copy nuclear genes (SCNs) from genome skimming data, by subsampling genome resequencing data and generating 10 data sets with different sequencing coverage in silico. We tested the performance of four data sets (plastome, nrDNA, mtDNA, and SCNs) obtained from genome skimming based on phylogenetic analyses of the Vitis clade at the genus level and Vitaceae at the family level, respectively. Our results showed that optimal minimum sequencing depth for high-quality SCNs assembly via genome skimming was about 10× coverage. Without the steps of synthesizing baits and enrichment experiments, coupled with incredibly low sequencing costs, we showcase that deep genome skimming (DGS) is as effective for capturing large data sets of SCNs as the widely used Hyb-Seq approach, in addition to capturing plastomes, mtDNA, and entire nrDNA repeats. DGS may serve as an efficient and economical alternative and may be superior to the popular target enrichment/Hyb-Seq approach.     

人类线粒体DNA变异的检测方法和思路

基于线粒体DNA（mtDNA)的研究对于人群源流迁移、线粒体相关疾病病因的探讨和法医鉴定等具有重意义,就检测人线粒体突变的一些常用方法,如RFLP、SSO和控制区测序等作一小结和归纳,并重点介绍目前mtDNA突变的筛选方法和思路,另外,还总结了近年来对人mtDNA方面的研究结果,对世界人群中主要单倍型类群（haplogroup)特征变异位点和相应的酶切检测引物作了归纳。     

Recombination or mutation rate heterogeneity? Implications for Mitochondrial Eve

The study of mitochondrial DNA (mtDNA) has helped to demonstrate the African origin of our species and the relationship between living humans and the Neanderthals. mtDNA data have also been used to establish the time and route of major events in human history, such as the expansion of Neolithic farmers into Europe, and the settlement of the Pacific and the New World. However, it is becoming apparent that mtDNA evolution is more complex than previously believed. Anomalous mutation patterns perturb phylogenetic assumptions based on mtDNA data. Although they are frequently dismissed as sequencing errors or mutation hotspots, some of the anomalies have no satisfactory explanation. The mechanisms behind apparent mutation rate heterogeneity, or even possible mtDNA recombination, remain unknown. These issues need to be addressed, as they have profound consequences for the interpretation of mtDNA data.     

Differentiating between monozygotic twins through next-generation mitochondrial genome sequencing

Monozygotic (MZ) twins, considered to be genetically identical, cannot be distinguished from one another by standard forensic DNA testing. A recent study employed whole genome sequencing to identify extremely rare mutations and reported that mutation analysis could be used to differentiate between MZ twins. Compared with nuclear DNA, mitochondrial DNA (mtDNA) has higher mutation rates; therefore, minor differences theoretically exist in MZ twins' mitochondrial genome (mtGenome). However, conventional Sanger-type sequencing (STS) is neither amenable to, nor feasible for, the detection of low-level sequence variants. The recent introduction of massively parallel sequencing (MPS) has the capability to sequence many targeted regions of multiple samples simultaneously with desirable depth of coverage. Thus, the aim of this study was to assess whether full mtGenome sequencing analysis can be used to differentiate between MZ twins. Ten sets of MZ twins provided blood samples that underwent extraction, quantification, mtDNA enrichment, library preparation, and ultra-deep sequencing. Point heteroplasmies were observed in eight sets of MZ twins, and a single nucleotide variant (nt15301) was detected in five sets of MZ twins. Thus, this study demonstrates that ultra-deep mtGenome sequencing could be used to differentiate between MZ twins.     

A Human Mitochondrial DNA Standard Reference Material for Quality Control in Forensic Identification,Medical Diagnosis,and Mutation Detection

A human mitochondrial DNA (mtDNA) standard reference material (SRM 2392) will provide quality control when mtDNA is sequenced for forensic identifications, medical diagnosis, or mutation detection. SRM 2392 includes DNA from two lymphoblast cell cultures (CHR and 9947A) and cloned DNA from the CHR HV1 region, which contains a C stretch and is difficult to sequence. The mtDNA sequence (but not the DNA) of a third human template GM03798 is provided for comparison. Fifty-eight unique primer sets allow any area or the entire mtDNA (16,569 bp) to be amplified and sequenced. While none of the differences in these three templates correspond to published mutations associated with specific diseases, some of these differences did result in animo acid changes compared with that published by S. Andersonet al.(1981,Nature290: 457–465). An interlaboratory evaluation of the amplification, sequencing, and data analysis of the CHR template was conducted by four laboratories. Corroboration of the SRM results will provide quality assurance that any unknown mtDNA is also being amplified and sequenced correctly.     

Mitochondrial DNA somatic mutation burden and heteroplasmy are associated with chronological age,smoking, and HIV infection

The gradual accumulation of mitochondrial DNA (mtDNA) mutations is implicated in aging and may contribute to the accelerated aging phenotype seen with tobacco smoking and HIV infection. mtDNA mutations are thought to arise from oxidative damage; however, recent reports implicate polymerase γ errors during mtDNA replication. Investigations of somatic mtDNA mutations have been hampered by technical challenges in measuring low‐frequency mutations. We use primer ID‐based next‐generation sequencing to quantify both somatic and heteroplasmic blood mtDNA point mutations within the D‐loop, in 164 women and girls aged 2–72 years, of whom 35% were smokers and 56% were HIV‐positive. Somatic mutations and the occurrence of heteroplasmic mutations increased with age. While transitions are theorized to result from polymerase γ errors, transversions are believed to arise from DNA oxidative damage. In our study, both transition and transversion mutations were associated with age. However, transition somatic mutations were more prevalent than transversions, and no heteroplasmic transversions were observed. We also measured elevated somatic mutations, but not heteroplasmy, in association with high peak HIV viremia. Conversely, heteroplasmy was higher among smokers, but somatic mutations were not, suggesting that smoking promotes the expansion of preexisting mutations rather than de novo mutations. Taken together, our results are consistent with blood mtDNA mutations increasing with age, inferring a greater contribution of polymerase γ errors in mtDNA mutagenesis. We further suggest that smoking and HIV infection both contribute to the accumulation of mtDNA mutations, though in different ways.     

Phylogeographic differentiation of mitochondrial DNA in Han Chinese

To characterize the mitochondrial DNA (mtDNA) variation in Han Chinese from several provinces of China, we have sequenced the two hypervariable segments of the control region and the segment spanning nucleotide positions 10171-10659 of the coding region, and we have identified a number of specific coding-region mutations by direct sequencing or restriction-fragment-length-polymorphism tests. This allows us to define new haplogroups (clades of the mtDNA phylogeny) and to dissect the Han mtDNA pool on a phylogenetic basis, which is a prerequisite for any fine-grained phylogeographic analysis, the interpretation of ancient mtDNA, or future complete mtDNA sequencing efforts. Some of the haplogroups under study differ considerably in frequencies across different provinces. The southernmost provinces show more pronounced contrasts in their regional Han mtDNA pools than the central and northern provinces. These and other features of the geographical distribution of the mtDNA haplogroups observed in the Han Chinese make an initial Paleolithic colonization from south to north plausible but would suggest subsequent migration events in China that mainly proceeded from north to south and east to west. Lumping together all regional Han mtDNA pools into one fictive general mtDNA pool or choosing one or two regional Han populations to represent all Han Chinese is inappropriate for prehistoric considerations as well as for forensic purposes or medical disease studies.     

Low "penetrance" of phylogenetic knowledge in mitochondrial disease studies

An up-to-date view of the worldwide mitochondrial DNA (mtDNA) phylogeny together with an evaluation of the conservation of each site is a reliable tool for detecting errors in mtDNA studies and assessing the functional importance of alleged pathogenic mutations. However, most of the published studies on mitochondrial diseases make very little use of the phylogenetic knowledge that is currently available. This drawback has two inadvertent consequences: first, there is no sufficient a posteriori quality assessment of complete mtDNA sequencing efforts; and second, no feedback is provided for the general mtDNA database when apparently new mtDNA lineages are discovered. We demonstrate, by way of example, these issues by reanalysing three mtDNA sequencing attempts, two from Europe and another one from East Asia. To further validate our phylogenetic deductions, we completely sequenced two mtDNAs from healthy subjects that nearly match the mtDNAs of two patients, whose sequences gave problematic results.     

Detecting Heteroplasmy from High-Throughput Sequencing of Complete Human Mitochondrial DNA Genomes

Heteroplasmy, the existence of multiple mtDNA types within an individual, has been previously detected by using mostly indirect methods and focusing largely on just the hypervariable segments of the control region. Next-generation sequencing technologies should enable studies of heteroplasmy across the entire mtDNA genome at much higher resolution, because many independent reads are generated for each position. However, the higher error rate associated with these technologies must be taken into consideration to avoid false detection of heteroplasmy. We used simulations and phiX174 sequence data to design criteria for accurate detection of heteroplasmy with the Illumina Genome Analyzer platform, and we used artificial mixtures and replicate data to test and refine the criteria. We then applied these criteria to mtDNA sequence reads for 131 individuals from five Eurasian populations that had been generated via a parallel tagged approach. We identified 37 heteroplasmies at 10% frequency or higher at 34 sites in 32 individuals. The mutational spectrum does not differ between heteroplasmic mutations and polymorphisms in the same individuals, but the relative mutation rate at heteroplasmic mutations is significantly higher than that estimated for all mutable sites in the human mtDNA genome. Moreover, there is also a significant excess of nonsynonymous mutations observed among heteroplasmies, compared to polymorphism data from the same individuals. Both mutation-drift and negative selection influence the fate of heteroplasmies to determine the polymorphism spectrum in humans. With appropriate criteria for avoiding false positives due to sequencing errors, next-generation technologies can provide novel insights into genome-wide aspects of mtDNA heteroplasmy.     

人类线粒体DNA异质性研究进展及其法医学应用

线粒体DNA(mitochondrial DNA mtDNA)的异质性自从被发现以来,一直被遗传学、进化学、发育遗传学以及法医遗传学、分子生物学领域所重视。由于线粒体异质性的存在,使得很多涉及疾病、进化、系统发育线粒体基因组与核基因组的相互作用关系、线粒体DNA复制机制以及法医学运用线粒体DNA进行实际案件评估的问题变得复杂化。此外线粒体DNA异质性的发生原因以及对线粒体异质性的检测方法标准化问题还没有一个统一的答案。针对线粒体DNA异质性带来的种种问题,近年来国内外取得了不少研究进展。     

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.     

Base composition analysis of human mitochondrial DNA using electrospray ionization mass spectrometry: a novel tool for the identification and differentiation of humans

In traditional approaches, mitochondrial DNA (mtDNA) variation is exploited for forensic identity testing by sequencing the two hypervariable regions of the human mtDNA control region. To reduce time and labor, single nucleotide polymorphism (SNP) assays are being sought to possibly replace sequencing. However, most SNP assays capture only a portion of the total variation within the desired regions, require a priori knowledge of the position of the SNP in the genome, and are generally not quantitative. Furthermore, with mtDNA, the clustering of SNPs complicates the design of SNP extension primers or hybridization probes. This article describes an automated electrospray ionization mass spectrometry method that can detect a number of clustered SNPs within an amplicon without a priori knowledge of specific SNP positions and can do so quantitatively. With this technique, the base composition of a PCR amplicon, less than 140 nucleotides in length, can be calculated. The difference in base composition between two samples indicates the presence of an SNP. Therefore, no post-PCR analytical construct needs to be developed to assess variation within a fragment. Of the 2754 different mtDNA sequences in the public forensic mtDNA database, nearly 90% could be resolved by the assay. The mass spectrometer is well suited to characterize and quantitate heteroplasmic samples or those containing mixtures. This makes possible the interpretation of mtDNA mixtures (as well as mixtures when assaying other SNPs). This assay can be expanded to assess genetic variation in the coding region of the mtDNA genome and can be automated to facilitate analysis of a large number of samples such as those encountered after a mass disaster.     

一种新的线粒体基因组DNA捕获探针的制备及初步应用

目的:建立新的线粒体基因组DNA杂交捕获探针制备方法并用进行初步应用。方法:通过PCR技术扩增特异线粒体DNA片段,并与生物素偶联,最后与标记磁珠的亲和素混合获得捕获探针。并自行制备的线粒体基因组DNA文库捕获探针与肝癌全基因组测序文库进行液相杂交。分离捕获产物后PCR扩增并进行测序分析。结果:成功建立了线粒体基因组杂交捕获探针制备方法并成功分离线粒体基因组DNA;对测序数据的分析显示:90%以上测序数据来自线粒体基因组DNA,且覆盖率达到100%,且均一性良好。检测到的同质性变异位点数量和异质性变异位点数量与全基因组测序数据产生的结果接近(P=0.9152,P=0.8409)。结论:新方法制备的线粒体基因组DNA杂交捕获探针可以从全基因组文库中高效捕获线粒体基因组DNA测序文库。     

Using mtDNA Sequences to Estimate SNP Parameters in ESTs

Discovery of single nucleotide polymorphisms (SNPs) requires analysis of redundant sequences such as those available in large public databases. The ability to detect SNPs, especially those of low frequency, is dependent on the depth and scale of the discovery effort. Large numbers of SNPs have been identified by mining large-scale EST surveys and whole genome sequencing projects. These surveys however are subject to ascertainment bias and the inherent errors in large-scale single pass sequencing efforts. For example, the number of steps involved in the construction and sequencing of cDNA libraries make ESTs highly error prone, resulting in an increased frequency of nonvalid SNPs obtained in these surveys. Sequences of mtDNA genes are often incorporated into cDNA libraries as an artifact of the library construction process and are typically either subtracted from cDNA libraries or are considered superfluous when evaluating the information content of EST datasets. Sequences of mtDNA genes provide a unique resource for the analysis of SNP parameters in EST projects. This study uses sequences from four turkey muscle cDNA libraries to demonstrate how mtDNA sequences gleaned from collections of ESTs can be used to estimate SNP parameters and thus help predict the validity of SNPs.     

Using mtDNA sequences to estimate SNP parameters in ESTs

Discovery of single nucleotide polymorphisms (SNPs) requires analysis of redundant sequences such as those available in large public databases. The ability to detect SNPs, especially those of low frequency, is dependent on the depth and scale of the discovery effort. Large numbers of SNPs have been identified by mining large-scale EST surveys and whole genome sequencing projects. These surveys however are subject to ascertainment bias and the inherent errors in large-scale single pass sequencing efforts. For example, the number of steps involved in the construction and sequencing of cDNA libraries make ESTs highly error prone, resulting in an increased frequency of nonvalid SNPs obtained in these surveys. Sequences of mtDNA genes are often incorporated into cDNA libraries as an artifact of the library construction process and are typically either subtracted from cDNA libraries or are considered superfluous when evaluating the information content of EST datasets. Sequences of mtDNA genes provide a unique resource for the analysis of SNP parameters in EST projects. This study uses sequences from four turkey muscle cDNA libraries to demonstrate how mtDNA sequences gleaned from collections of ESTs can be used to estimate SNP parameters and thus help predict the validity of SNPs.