High-quality and high-throughput massively parallel sequencing of the human mitochondrial genome using the Illumina MiSeq

Mitochondrial genome study of Camellia oleifera revealed the tandem conserved gene cluster of nad5-nads in evolution

Camellia oleifera is a kind of high-quality oil supply species. Its seeds contain rich unsaturated fatty acids and antioxidant active ingredients, which is a kind of high-quality edible oil. In this study, we used bioinformatics methods to decipher a hexaploid Camellia oil tree's mitochondrial (mt) genome based on second-generation sequencing data. A 709,596 bp circular map of C. oleifera mt genome was found for the first time. And 74 genes were annotated in the whole genome. Mt genomes of C. oleifera and three Theaceae species had regions with high similarity, including gene composition and gene sequence. At the same time, five conserved gene pairs were found in 20 species. In all of the mt genomes, most of nad genes existed in tandem pairs. In addition, the species classification result, which, according to the gene differences in tandem with nad5 genes, was consistent with the phylogenetic tree. These initial results provide a valuable basis for the further researches of Camellia oleifera and a reference for the systematic evolution of plant mt genomes.

Complete Mitochondrial DNA Genome Variation in the Swedish Population

The development of complete mitochondrial genome (mitogenome) reference data for inclusion in publicly available population databases is currently underway, and the generation of more high-quality mitogenomes will only enhance the statistical power of this forensically useful locus. To characterize mitogenome variation in Sweden, the mitochondrial DNA (mtDNA) reads from the SweGen whole genome sequencing (WGS) dataset were analyzed. To overcome the interference from low-frequency nuclear mtDNA segments (NUMTs), a 10% variant frequency threshold was applied for the analysis. In total, 934 forensic-quality mitogenome haplotypes were characterized. Almost 45% of the SweGen haplotypes belonged to haplogroup H. Nearly all mitogenome haplotypes (99.1%) were assigned to European haplogroups, which was expected based on previous mtDNA studies of the Swedish population. There were signature northern Swedish and Finnish haplogroups observed in the dataset (e.g., U5b1, W1a), consistent with the nuclear DNA analyses of the SweGen data. The complete mitogenome analysis resulted in high haplotype diversity (0.9996) with a random match probability of 0.15%. Overall, the SweGen mitogenomes provide a large mtDNA reference dataset for the Swedish population and also contribute to the effort to estimate global mitogenome haplotype frequencies.

Exploring rare differences in mitochondrial genome between MZ twins using Ion Torrent semiconductor sequencing

Monozygotic (MZ) twins are considered to be genetically identical in that they have the same genomic DNA sequences in theory, and thus cannot be differentiated using forensic standard STR-based DNA profiling. However, a recent study employed deep sequencing to explore extremely rare mutations in the nuclear genome and reported that the mutation analysis could be applied to differentiate between MZ twins. Compared with the nuclear genome, the mitochondrial DNA (mtDNA) exhibits higher mutation rates due to fewer DNA repair mechanisms in the mitochondrial genome (mtGenome) and the lack of proofreading capability of the mtDNA polymerase. In a previous study, we used Illumina ultra-deep sequencing to describe point heteroplasmy (PHP) and nucleotide variant of the mtGenomes in venous blood samples of MZ twins. In the present study, we characterized minor differences of the mtGenomes in three tissue samples from seven sets of MZ twins using Ion Torrent semiconductor sequencing (Thermo Fisher Ion S5 XL system) and commercialized mtGenome sequencing kit (Precision ID mtDNA Whole Genome Panel). PHP was observed in blood samples from one set of MZ twins and in saliva samples from two sets of twins, but it presented in hair shaft samples from all seven sets of MZ twins. Overall, the coding region of the mtGenome exhibits more PHPs than the control region. The results of this study have further attested the competence of mtGenome sequencing in differentiating between MZ twins, and that among the three kinds of samples tested, hair shaft is more likely to accumulate minor differences in the mtGenomes of MZ twins.

Whole Mitochondrial Genome Detection and Analysis of Two- to Four-Generation Maternal Pedigrees Using a New Massively Parallel Sequencing Panel

Mitochondrial DNA (mtDNA) is an effective genetic marker in forensic practice, especially for aged bones and hair shafts. Detection of the whole mitochondrial genome (mtGenome) using traditional Sanger-type sequencing is laborious and time-consuming. Additionally, its ability to distinguish point heteroplasmy (PHP) and length heteroplasmy (LHP) is limited. The application of massively parallel sequencing in mtDNA detection helps researchers to study the mtGenome in-depth. The ForenSeq mtDNA Whole Genome Kit, which contains a total of 245 short amplicons, is one of the multiplex library preparation kits for the mtGenome. We used this system to detect the mtGenome in the blood samples and hair shafts of thirty-three individuals from eight two-generation pedigrees, one three-generation pedigree, and one four-generation pedigree. High-quality sequencing results were obtained. Ten unique mtGenome haplotypes were observed in the mothers from the ten pedigrees. A total of 26 PHPs were observed using the interpretation threshold of 6%. Eleven types of LHPs in six regions were evaluated in detail. When considering homoplasmic variants only, consistent mtGenome haplotypes were observed between the twice-sequenced libraries and between the blood and hair shafts from the same individual and among maternal relatives in the pedigrees. Four inherited PHPs were observed, and the remainder were de novo/disappearing PHPs in the pedigrees. Our results demonstrate the effective capability of the ForenSeq mtDNA Whole Genome Kit to generate the complete mtGenome in blood and hair shafts, as well as the complexity of mtDNA haplotype comparisons between different types of maternal relatives when heteroplasmy is considered.

Developmental validation of the ForenSeq MainstAY kit, MiSeq FGx sequencing system and ForenSeq Universal Analysis Software

For human identification purposes, forensic genetics has primarily relied upon a core set of autosomal (and to a lesser extent Y chromosome) short tandem repeat (STR) markers that are enriched by amplification using the polymerase chain reaction (PCR) that are subsequently separated and detected using capillary electrophoresis (CE). While STR typing conducted in this manner is well-developed and robust, advances in molecular biology that have occurred over the last 15 years, in particular massively parallel sequencing (MPS) [1-7], offer certain advantages as compared to CE-based typing. First and foremost is the high throughput capacity of MPS. Current bench top high throughput sequencers enable larger batteries of markers to be multiplexed and multiple samples to be sequenced simultaneously (e.g., millions to billions of nucleotides can be sequenced in one run). Second, compared to the length-based CE approach, sequencing STRs increases discrimination power, enhances sensitivity of detection, reduces noise due to instrumentation, and improves mixture interpretation [4,8-23]. Third, since detection of STRs is based on sequence and not fluorescence, amplicons can be designed that are shorter in length and of similar lengths among loci, where possible, which can improve amplification efficiency and analysis of degraded samples. Lastly, MPS offers a single format approach that can be applied to analysis of a wide variety of genetic markers of forensic interest (e.g., STRs, mitochondrial DNA, single nucleotide polymorphisms, insertion/deletions). These features make MPS a desirable technology for casework [14,15,24,25-48]. The developmental validation of the ForenSeq MainstAY library preparation kit with the MiSeq FGx Sequencing System and ForenSeq Universal Software is reported here to assist with validation of this MPS system for casework [49]. The results show that the system is sensitive, accurate and precise, specific, and performs well with mixtures and mock case-type samples.

DNA and protein analyses of hair in forensic genetics

Hair is one of the most common pieces of biological evidence found at a crime scene and plays an essential role in forensic investigation. Hairs, especially non-follicular hairs, are usually found at various crime scenes, either by natural shedding or by forcible shedding. However, the genetic material in hairs is usually highly degraded, which makes forensic analysis difficult. As a result, the value of hair has not been fully exploited in forensic investigations and trials. In recent years, with advances in molecular biology, forensic analysis of hair has achieved remarkable strides and provided crucial clues in numerous cases. This article reviews recent developments in DNA and protein analysis of hair and attempts to provide a comprehensive solution to improve forensic hair analysis.

Routine Mitogenome MPS Analysis from 1 and 5 mm of Rootless Human Hair

While hair shafts are a common evidence type in forensic cases, they are often excluded from DNA analysis due to their limited DNA quantity and quality. Mitochondrial (mt) DNA sequencing is the method of choice when working with rootless hair shaft fragments due to the elevated copy number of mtDNA and the highly degraded nature of nuclear (n) DNA. Using massively parallel sequencing (MPS) of the mitochondrial (mito) genome, we studied the impact of hair age (time since collection) and physical characteristics (hair diameter, medullary structure, and length of hair tested) on mtDNA recovery and MPS data quality. Hair shaft cuttings of 1 and 5 mm from hairs less than five years to 46 years of age from 60 donors were characterized microscopically. Mitogenome sequences were generated using the Promega PowerSeqTM Whole Mito System prototype kit and the Illumina MiSeq instrument. Reportable mitogenome sequences were obtained from all hairs up to 27 years of age (37 donors), with at least 98% of the mitogenome reported for more than 94% of the 74 hair samples analyzed; the minimum reported sequence was 88%. Furthermore, data from the 1 and 5 mm replicates gave concordant haplotypes. As expected, mtDNA yield decreased, mtDNA degradation increased, and mitogenome MPS data quality declined as the age of the hair increased. Hair diameter and medullary structure had minimal impact on yield and data quality. Our findings support that MPS is a robust and reliable method for routinely generating mitogenome sequences from 1 and 5 mm hair shaft samples up to 27 years of age, which is of interest to the forensic community, biological anthropologists, and medical geneticists.

Evaluating the suitability of current mitochondrial DNA interpretation guidelines for multigenerational whole mitochondrial genome comparisons

Sanger sequencing of the mitochondrial DNA (mtDNA) control region was previously the only method available for forensic casework involving degraded samples from skeletal remains. The introduction of Next Generation Sequencing (NGS) has transformed genetic data generation and human identification using mtDNA. Whole mitochondrial genome (mtGenome) analysis is now being introduced into forensic laboratories around the world to analyze historical remains. Research into large pedigrees using the mtGenome is critical to evaluate currently available interpretation guidelines for mtDNA analysis, which were developed for comparisons using the control region. This study included mtGenomes from 225 individuals from the last four generations of the Norfolk Island (NI) genetic isolate pedigree consisting of 49 distinct maternal lineages. The data from these individuals were arranged into 2339 maternally related pairs separated by up to 18 meioses. Our results show that 97.3% of maternally related pairs were concordant at all nucleotide positions, resulting in the correct interpretation of “Cannot Exclude”; 2.7% of pairs produced an “Inconclusive” result, and there were no instances of false exclusion. While these results indicate that existing guidelines are suitable for multigenerational whole mtGenome analysis, we recommend caution be taken when classifying heteroplasmic changes as differences for human identification. Our data showed the classification of heteroplasmic changes as differences increases the prevalence of inconclusive identification by 6%, with false exclusions observed in 0.34% of pairs examined. Further studies of multigenerational pedigrees, however, are needed to validate mtGenome interpretation guidelines for historical case work to more fully utilize emerging advancements.

Next-Generation Bisulfite Sequencing for Targeted DNA Methylation Analysis

Bisulfite sequencing is the "gold-standard" technique for DNA methylation analysis. By combining bisulfite sequencing with high-throughput, next-generation sequencing technology, we can document methylation from many thousands of individual reads (equivalent to alleles or "cells"), for multiple target regions and from many samples simultaneously. Here, we describe a next-generation bisulfite-sequencing assay for targeted DNA methylation analysis which offers scope for the simultaneous interrogation of multiple genomic loci across numerous samples.

Forensic nanopore sequencing of microhaplotype markers using QitanTech's QNome

In recent years, extraordinary progress has been made in genome sequencing technologies, which has led to a decrease in cost and an increase in the diversity of sequenced genomes. Nanopore sequencing is one of the latest genome sequencing technologies. It aims to sequence longer contiguous pieces of DNA, which are essential for resolving structurally complex regions, and provides a new approach for forensic genetics to detect longer markers in real time. To date, multiple studies have been conducted to sequence forensic markers using MinION from Oxford Nanopore Technologies (ONT), and the results indicate that nanopore sequencing holds promise for forensic applications. Qitan Technology (QitanTech) recently launched its first commercial nanopore genome sequencer, QNome. It could achieve a read length of more than 150 kbp, and could generate approximately 500 Mb of data in 8 h. In this pilot study, we explored and validated this alternative nanopore sequencing device for microhaplotype (MH) profiling using a custom set of 15 MH loci. Seventy single-contributor samples were divided into 7 batches, each of which included 10 samples and control DNA 9947A and was sequenced by QNome. MH genotypes generated from QNome were compared to those from Ion Torrent sequencing (Ion S5XL system) to evaluate the accuracy and stability. Twelve samples randomly selected from the last three batches and Control DNA 9947A were also subjected to ONT MinION sequencing (with R9.4 flow cell) for parallel comparison. Based on MHtyper, a bioinformatics workflow developed for automated MH designation, all MH loci can be genotyped and reliably phased using the QNome data, with an overall accuracy of 99.83% (4 errors among 2310 genotypes). Three occurred near or in the region of homopolymer sequences, and one existed within 50 bp of the start of the sequencing reaction. In the last 15 samples (12 individual samples and 3 replicates of control DNA 9947A), two SNPs located at 4-mer homopolymers failed to obtain reliable genotypes on the MinION data. This study shows the potential of state-of-the-art nanopore sequencing methods to analyze forensic MH markers. Given the rapid pace of change, sporadic and nonrepetitive errors presented in this study are expected to be resolved by further developments of nanopore technologies and analysis tools.

Mitochondrial genome sequencing with short overlapping amplicons on MiSeq FGx system

With the development and maturation of massively parallel sequencing (MPS) technology, the mitochondrial genome (mitogenome) sequencing is increasingly applied in the forensic field. In this study, we employed the strategy of short overlapping amplicons for the whole mitogenome, library preparation with tagmentation using the Nextera® XT DNA Library Preparation Kit, sequencing on the MiSeq FGxTM Forensic Genomics System and analyzing data using the mitochondrial(mtDNA) MSR Plug-in and the mtDNA Variant Analyzer. A total of 27 libraries and 56 libraries were sequenced in a run using MiSeq Reagent Kit v2 and v3, respectively. Results showed more than 1800 × of averaged depth of coverage (DoC) at each position. Concordant haplotypes of 9947 A and 2800 M were obtained at 32 variants. Cross-reactivity was observed with 1 ng primate DNA and 10 ng non-primate DNA but could be easily distinguished. Full and accurate variants were obtained from at least 50 pg input DNA and from minor contributors between 19:1 and 1:19 mixed ratios with known reference profiles. More than 86% variants were detected from ≥200-bp degraded samples but its haplotype was assigned to more ancestral haplogroup. Further, a total of 3 962 variants were observed at 613 nucleotide positions from 103 Xibe mitogenomes with 25:1 ratio of transitions to transversions. Two new transversions (C13735A and A14755C) and two tri-alleles at nps 9824 and 16092 were identified. There were 103 unique mitogenome haplotypes from 103 Chinese Xibe that were assigned to 79 haplogroups. Haplogroup D was the preponderant top-level haplogroup in Xibe followed by F, B, M, A, N, G, C, Z, Y, HV and J. Random match probability (RMP) and haplotype diversity (HD) of the whole mitogenome was calculated as 0.0097 and 1.0000, respectively. Compared with HVS-I only, RMP decreased 33.56%, while the number of haplotypes and HD increased 15.73% and 0.49%, respectively. Principal component analysis (PCA) showed that Xibe was clustered to East and Southeast Asian. As a whole, this MPS strategy is suitable for the whole mitogenome sequencing especially for degraded samples and can facilitate generating mitogenome data to support the routine application in forensic sciences. EMP00726 is the first whole mitogenome dataset from Xibe contributed to the EMPOP. Supplemental data for this article are available online at.

MMDIT: A tool for the deconvolution and interpretation of mitochondrial DNA mixtures

Short tandem repeats of the nuclear genome have been the preferred markers for analyzing forensic DNA mixtures. However, when nuclear DNA in a sample is degraded or limited, mitochondrial DNA (mtDNA) markers provide a powerful alternative. Though historically considered challenging, the interpretation and analysis of mtDNA mixtures have recently seen renewed interest with the advent of massively parallel sequencing. However, there are only a few software tools available for mtDNA mixture interpretation. To address this gap, the Mitochondrial Mixture Deconvolution and Interpretation Tool (MMDIT) was developed. MMDIT is an interactive application complete with a graphical user interface that allows users to deconvolve mtDNA (whole or partial genomes) mixtures into constituent donor haplotypes and estimate random match probabilities on these resultant haplotypes. In cases where deconvolution might not be feasible, the software allows mixture analysis directly within a binary framework (i.e. qualitatively, only using data on allele presence/absence). This paper explains the functionality of MMDIT, using an example of an in vitro two-person mtDNA mixture with a ratio of 1:4. The uniqueness of MMDIT lies in its ability to resolve mixtures into complete donor haplotypes using a statistical phasing framework before mixture analysis and evaluating statistical weights employing a novel graph algorithm approach. MMDIT is the first available open-source software that can automate mtDNA mixture deconvolution and analysis. The MMDIT web application can be accessed online at https://www.unthsc.edu/mmdit/. The source code is available at https://github.com/SammedMandape/MMDIT_UI and archived on zenodo (https://doi.org/10.5281/zenodo.4770184).

Mitochondrial DNA in forensic use

Genetic analysis of mitochondrial DNA (mtDNA) has always been a useful tool for forensic geneticists, mainly because of its ubiquitous presence in biological material, even in the absence of nuclear DNA. Sequencing, however, is not a skill that is part of the routine forensic analysis because of the relative rarity of requests, and the need for retention of necessary skill sets and associated accreditation issues. While standard Sanger sequencing may be relatively simple, many requests are made in the face of compromised biological samples. Newer technologies, provided through massively parallel sequencing (MPS), will increase the opportunity for scientists to include this tool in their routine, particularly for missing person investigations. MPS has also enabled a different approach to sequencing that can increase sensitivity in a more targeted approach. In these circumstances it is likely that only a laboratory that specialises in undertaking forensic mtDNA analysis will be able to take these difficult cases forward, more so because reviews of the literature have revealed significantly high levels of typing errors in publications reporting mtDNA sequences. The forensic community has set out important guidelines, not only in the practical aspects of analysis, but also in the interpretation of that sequence to ensure that accurate comparisons can be made. Analysis of low-level, compromised and ancient DNA is not easy, however, as contamination is extremely difficult to eliminate and circumstances leading to sequencing errors are all too easily introduced. These problems, and solutions, are discussed in the article in relation to several historic cases.

Genetic and phylogeographic evidence for Jewish Holocaust victims at the Sobibór death camp

Six million Jews were killed by Nazi Germany and its collaborators during World War II. Archaeological excavations in the area of the death camp in Sobibór, Poland, revealed ten sets of human skeletal remains presumptively assigned to Polish victims of the totalitarian regimes. However, their genetic analyses indicate that the remains are of Ashkenazi Jews murdered as part of the mass extermination of European Jews by the Nazi regime and not of otherwise hypothesised non-Jewish partisan combatants. In accordance with traditional Jewish rite, the remains were reburied in the presence of a Rabbi at the place of their discovery.

Massively parallel sequencing of human skeletal remains in Vietnam using the precision ID mtDNA control region panel on the Ion S5™ system

Mitochondrial DNA (mtDNA) analysis using Sanger sequencing has been a routine practice for the identification of human skeletal remains. However, this process is usually challenging since DNA from the remains is highly degraded and at low concentration. Recently, the advent and implementation of massively parallel sequencing (MPS) have been offered the ability to improve mtDNA sequence data for forensic analysis. To assess the utility of the Ion S5™ system - an MPS platform for mtDNA analysis in challenging samples, we sequenced the mitochondrial control region of 52 age-old skeletal remains. Using the Precision ID mtDNA Control Region Panel, 50 full and two partial control region haplotypes at relatively high mean coverage of 2494 × were achieved for variant calling. Further variant analysis at 10% threshold for point heteroplasmy showed high degradation degree in terms of DNA damage in our bone samples. A higher point heteroplasmy threshold of 20% was required to diminish most of background noise caused by the damage. The results from this study indicated the potential application of the Ion S5™ system in sequencing degraded samples in Vietnam and provided valuable data sources for forensic analyses in the future.

Resolution of mitochondrial DNA mixtures using a probe capture next generation sequencing system and phylogenetic-based software

Interpreting mixtures with nuclear genetic markers remains one of the persisting challenges in forensic DNA analysis, particularly when the DNA is degraded or present in trace amounts. In these scenarios, analyzing mitochondrial (mt) DNA can be useful due to the higher copy number per cell compared to nuclear DNA. However, until the emergence of Next-Generation Sequencing (NGS) with its clonal sequencing capability, analysis of mtDNA mixtures was very challenging. A mitochondrial genome probe capture Next-Generation Sequencing (NGS) system was used to sequence complex mtDNA mixtures and two different software programs to analyze the sequence data. Analysis of contrived mixtures of two contributors in 50:50 and 95:5 ratios as well as three-person mixtures ranging from near equal proportions (~33:33:33 ratio) to low amounts of the minor contributors (e.g., a 90:5:5 ratio) is reported. This system was also applied to the analysis of mtDNA mixtures from forensically relevant samples. For data analysis, both the variant frequency-based software program GeneMarker®HTS and the phylogenetic-based software program Mixemt was used to de-convolute the mixtures. Using the massively parallel, clonal features of NGS, one can bioinformatically separate and count the individual sequence reads to calculate the proportions of individual contributors using phylogenetically informative polymorphisms. GeneMarker®HTS allows us to detect all mutations, including "private" mutations (non-phylogenetically informative polymorphisms) and assign them to individual contributors based on the frequency of the sequence reads, provided that the proportions of the various contributors are sufficiently different. Using a probe capture NGS system and both GeneMarker®HTS and Mixemt software programs, the interpretation of complex mixtures of equal proportion contributors, trace amount contributors, and more than two contributors in contrived mixtures, as well as interpretation of challenging forensic specimens is demonstrated.

Human Mitochondrial Control Region and mtGenome: Design and Forensic Validation of NGS Multiplexes, Sequencing and Analytical Software

Forensic mitochondrial DNA (mtDNA) analysis conducted using next-generation sequencing (NGS), also known as massively parallel sequencing (MPS), as compared to Sanger-type sequencing brings modern advantages, such as deep coverage per base (herein referred to as read depth per base pair (bp)), simultaneous sequencing of multiple samples (libraries) and increased operational efficiencies. This report describes the design and developmental validation, according to forensic quality assurance standards, of end-to-end workflows for two multiplexes, comprised of ForenSeq mtDNA control region and mtDNA whole-genome kits the MiSeq FGx^TM instrument and ForenSeq universal analysis software (UAS) 2.0/2.1. Polymerase chain reaction (PCR) enrichment and a tiled amplicon approach target small, overlapping amplicons (60–150 bp and 60–209 bp for the control region and mtGenome, respectively). The system provides convenient access to data files that can be used outside of the UAS if desired. Studies assessed a range of environmental and situational variables, including but not limited to buccal samples, rootless hairs, dental and skeletal remains, concordance of control region typing between the two multiplexes and as compared to orthogonal data, assorted sensitivity studies, two-person DNA mixtures and PCR-based performance testing. Limitations of the system and implementation considerations are discussed. Data indicated that the two mtDNA multiplexes, MiSeq FGx and ForenSeq software, meet or exceed forensic DNA quality assurance (QA) guidelines with robust, reproducible performance on samples of various quantities and qualities.

Ancient and modern mitogenomes from Central Argentina: new insights into population continuity, temporal depth and migration in South America

The inverted triangle shape of South America places Argentina territory as a geographical crossroads between the two principal peopling streams that followed either the Pacific or the Atlantic coasts, which could have then merged in Central Argentina (CA). Although the genetic diversity from this region is therefore crucial to decipher past population movements in South America, its characterization has been overlooked so far. We report 92 modern and 22 ancient mitogenomes spanning a temporal range of 5000 years, which were compared with a large set of previously reported data. Leveraging this dataset representative of the mitochondrial diversity of the subcontinent, we investigate the maternal history of CA populations within a wider geographical context. We describe a large number of novel clades within the mitochondrial DNA tree, thus providing new phylogenetic interpretations for South America. We also identify several local clades of great temporal depth with continuity until the present time, which stem directly from the founder haplotypes, suggesting that they originated in the region and expanded from there. Moreover, the presence of lineages characteristic of other South American regions reveals the existence of gene flow to CA. Finally, we report some lineages with discontinuous distribution across the Americas, which suggest the persistence of relic lineages likely linked to the first population arrivals. The present study represents to date the most exhaustive attempt to elaborate a Native American genetic map from modern and ancient complete mitochondrial genomes in Argentina and provides relevant information about the general process of settlement in South America.

Interpreting NUMTs in forensic genetics: Seeing the forest for the trees

Nuclear mitochondrial DNA (mtDNA) segments (NUMTs) were discovered shortly after sequencing the first human mitochondrial genome. They have earlier been considered to represent archaic elements of ancient insertion events, but modern sequencing technologies and growing databases of mtDNA and NUMT sequences confirm that they are abundant and some of them phylogenetically young. Here, we build upon mtDNA/NUMT review articles published in the mid 2010 s and focus on the distinction of NUMTs and other artefacts that can be observed in aligned sequence reads, such as mixtures (contamination), point heteroplasmy, sequencing error and cytosine deamination. We show practical examples of the effect of the mtDNA enrichment method on the representation of NUMTs in the mapped sequence data and discuss methods to bioinformatically filter NUMTs from mtDNA reads.

Analysis of microbial diversity and functional differences in different types of high-temperature Daqu

Bacterial communities that enrich in high-temperature Daqu are important for the Chinese maotai-flavor liquor brewing process. However, the bacterial communities in three different types of high-temperature Daqu (white Daqu, black Daqu, and yellow Daqu) are still undercharacterized. In this study, the bacterial diversity of three different types of high-temperature Daqu was investigated using Illumina MiSeq high-throughput sequencing. The bacterial community of high-temperature Daqu is mainly composed of thermophilic bacteria, and seven bacterial phyla along with 262 bacterial genera were identified in all 30 high-temperature Daqu samples. Firmicutes, Actinobacteria, Proteobacteria, and Acidobacteria were the dominant bacterial phyla in high-temperature Daqu samples, while Thermoactinomyces, Staphylococcus, Lentibacillus, Bacillus, Kroppenstedtia, Saccharopolyspora, Streptomyces, and Brevibacterium were the dominant bacterial genera. The bacterial community structure of three different types of high-temperature Daqu was significantly different (p < .05). In addition, the results of microbiome phenotype prediction by BugBase and bacterial functional potential prediction using PICRUSt show that bacteria from different types of high-temperature Daqu have similar functions as well as phenotypes, and bacteria in high-temperature Daqu have vigorous metabolism in the transport and decomposition of amino acids and carbohydrates. These results offer a reference for the comprehensive understanding of bacterial diversity of high-temperature Daqu.

Assessment of Illumina® Human mtDNA Genome assay: workflow evaluation with development of analysis and interpretation guidelines

Mitochondrial DNA (mtDNA) is a small but significant part of the human genome, whose applicability potential has gradually increased with the advent of massively parallel sequencing (MPS) technology. Knowledge of the particular workflow, equipment, and reagents used, along with extensive usage of negative controls to monitor all preparation steps constitute the prerequisites for confident reporting of results. In this study, we performed an assessment of Illumina® Human mtDNA Genome assay on MiSeq FGx™ instrument. Through analysis of several types of negative controls, as well as mtDNA positive controls, we established thresholds for data analysis and interpretation, consisting of several components: minimum read depth (220 reads), minimum quality score (41), percentage of minor allele sufficient for analysis (3.0%), percentage of minor allele sufficient for interpretation (6.0%), and percentage of major allele sufficient for homoplasmic variant call (97.0%). Based on these criteria, we defined internal guidelines for analysis and interpretation of mtDNA results obtained by MPS. Our study shows that the whole mtDNA assay on MiSeq FGx™ produces repeatable and reproducible results, independent of the analyst, which are also concordant with Sanger-type sequencing results for mtDNA control region, as well as with MPS results produced by NextSeq®. Overall, established thresholds and interpretation guidelines were successfully applied for the sequencing of complete mitochondrial genomes from high-quality samples. The underlying principles and proposed methodology on the definition of internal laboratory guidelines for analysis and interpretation of MPS results may be applicable to similar MPS workflows, e.g. targeting good-quality samples in forensic genetics and molecular diagnostics.

Graph Algorithms for Mixture Interpretation

The scale of genetic methods are presently being expanded: forensic genetic assays previously were limited to tens of loci, but now technologies allow for a transition to forensic genomic approaches that assess thousands to millions of loci. However, there are subtle distinctions between genetic assays and their genomic counterparts (especially in the context of forensics). For instance, forensic genetic approaches tend to describe a locus as a haplotype, be it a microhaplotype or a short tandem repeat with its accompanying flanking information. In contrast, genomic assays tend to provide not haplotypes but sequence variants or differences, variants which in turn describe how the alleles apparently differ from the reference sequence. By the given construction, mitochondrial genetic assays can be thought of as genomic as they often describe genetic differences in a similar way. The mitochondrial genetics literature makes clear that sequence differences, unlike the haplotypes they encode, are not comparable to each other. Different alignment algorithms and different variant calling conventions may cause the same haplotype to be encoded in multiple ways. This ambiguity can affect evidence and reference profile comparisons as well as how “match” statistics are computed. In this study, a graph algorithm is described (and implemented in the MMDIT (Mitochondrial Mixture Database and Interpretation Tool) R package) that permits the assessment of forensic match statistics on mitochondrial DNA mixtures in a way that is invariant to both the variant calling conventions followed and the alignment parameters considered. The algorithm described, given a few modest constraints, can be used to compute the “random man not excluded” statistic or the likelihood ratio. The performance of the approach is assessed in in silico mitochondrial DNA mixtures.

A Continuous Statistical Phasing Framework for the Analysis of Forensic Mitochondrial DNA Mixtures

Despite the benefits of quantitative data generated by massively parallel sequencing, resolving mitotypes from mixtures occurring in certain ratios remains challenging. In this study, a bioinformatic mixture deconvolution method centered on population-based phasing was developed and validated. The method was first tested on 270 in silico two-person mixtures varying in mixture proportions. An assortment of external reference panels containing information on haplotypic variation (from similar and different haplogroups) was leveraged to assess the effect of panel composition on phasing accuracy. Building on these simulations, mitochondrial genomes from the Human Mitochondrial DataBase were sourced to populate the panels and key parameter values were identified by deconvolving an additional 7290 in silico two-person mixtures. Finally, employing an optimized reference panel and phasing parameters, the approach was validated with in vitro two-person mixtures with differing proportions. Deconvolution was most accurate when the haplotypes in the mixture were similar to haplotypes present in the reference panel and when the mixture ratios were neither highly imbalanced nor subequal (e.g., 4:1). Overall, errors in haplotype estimation were largely bounded by the accuracy of the mixture's genotype results. The proposed framework is the first available approach that automates the reconstruction of complete individual mitotypes from mixtures, even in ratios that have traditionally been considered problematic.

An 18 Multi-InDels panel for analysis of highly degraded forensic biological samples

DNA genotyping from trace and highly degraded biological samples is one of the most significant challenges of forensic DNA identification. There is a lack of simple and effective methods for genotyping highly degraded samples. In this study, a multiple loci insertion/deletion polymorphisms (Multi-InDels) panel was designed for detecting 18 autosomal Multi-InDels through capillary electrophoresis (CE) with amplicon sizes no longer than 125 bp. Studies of sensitivity, degradation, and species specificity were performed and a population study was carried out using 192 samples from Han populations in Hunan province in the south of China. The combined random match probability (CMP) of these 18 Multi-InDels was 3.23 × 10^-12 and the cumulative probability of exclusion (CPE) was 0.9989, suggesting this panel could be used independently for human identification and could provide efficient supporting information for parentage testing. Complete profiles were obtained from as low as 62.5 pg of total input DNA after increasing the number of PCR cycles. Moreover, all alleles were detected from artificially highly degraded DNA after 80 min of boiling water bath treatment. This 18 Multi-InDels panel is simple, fast, and effective for the forensic analysis of highly degraded DNA.

Bacterial diversity and flavor profile of Zha-Chili, a traditional fermented food in China

Zha-chili is a traditional Chinese fermented food with special flavor, which is often used as an appetizer in condiments and an important energy source. The final quality of zha-chili is closely related to its microbial community structure. However, the differences of bacterial diversity in zha-chili from different regions and how bacterial species affect zha-chili fermentation process and flavor quality have not been reported. In this study, the bacterial diversity and flavor quality of zha-chili samples from different regions were analyzed using Illumina Miseq high-throughput sequencing, electronic nose and electronic tongue technology. Twenty-three bacterial phyla and 665 bacterial genera were identified in all zha-chili samples. Firmicutes, Proteobacteria and Actinobacteria were the dominant bacterial phyla in zha-chili samples, while Lactobacillus, Pseudomonas, Pediococcus, Weissella and Staphylococcus were the dominant bacterial genera. The bacterial community structure of zha-chili samples from different regions was significantly diverse (p < 0.05). The flavor of zha-chili samples also varied in different regions, and the discrepancy of taste was much greater than that of aroma. Moreover, there were significant correlations (p < 0.05) between 6 dominant bacterial genera and 8 flavor indicators (3 aroma indicators, 5 taste indicators). In addition, the results of microbiome phenotypes prediction by BugBase and bacterial functional potential prediction using PICRUSt showed that eight out of nine predicted phenotypic functions of zha-chili samples from different regions were significantly different (p < 0.05), bacterial metabolism was vigorous in the zha-chili samples, and Lactobacillus was the dominant bacterial genus involved in metabolism during fermentation.

Developmental Validation of a MPS Workflow with a PCR-Based Short Amplicon Whole Mitochondrial Genome Panel

For the adoption of massively parallel sequencing (MPS) systems by forensic laboratories, validation studies on specific workflows are needed to support the feasibility of implementation and the reliability of the data they produce. As such, the whole mitochondrial genome sequencing methodology-Precision ID mtDNA Whole Genome Panel, Ion Chef, Ion S5, and Converge-has been subjected to a variety of developmental validation studies. These validation studies were completed in accordance with the Scientific Working Group on DNA Analysis Methods (SWGDAM) validation guidelines and assessed reproducibility, repeatability, accuracy, sensitivity, specificity to human DNA, and ability to analyze challenging (e.g., mixed, degraded, or low quantity) samples. Intra- and inter-run replicates produced an average maximum pairwise difference in variant frequency of 1.2%. Concordance with data generated with traditional Sanger sequencing and an orthogonal MPS platform methodology was used to assess accuracy, and generation of complete and concordant haplotypes at DNA input levels as low as 37.5 pg of nuclear DNA or 187.5 mitochondrial genome copies illustrated the sensitivity of the system. Overall, data presented herein demonstrate that highly accurate and reproducible results were generated for a variety of sample qualities and quantities, supporting the reliability of this specific whole genome mitochondrial DNA MPS system for analysis of forensic biological evidence.

Mitochondrial Sequencing of Missing Persons DNA Casework by Implementing Thermo Fisher's Precision ID mtDNA Whole Genome Assay

The advent of massively parallel sequencing (MPS) in the past decade has opened the doors to mitochondrial whole-genome sequencing. Mitochondrial (mt) DNA is used in forensics due to its high copy number per cell and maternal mode of inheritance. Consequently, we have implemented the Thermo Fisher Precision ID mtDNA Whole Genome panel coupled with the Ion Chef™ and Ion S5™ for MPS analysis in the California Department of Justice, Missing Persons DNA Program. Thirty-one mostly challenging samples (degraded, inhibited, low template, or mixed) were evaluated for this study. The majority of these samples generated single source full or partial genome sequences with MPS, providing information in cases where previously there was none. The quantitative and sensitive nature of MPS analysis was beneficial, but also led to detection of low-level contaminants. In addition, we found Precision ID to be more susceptible to inhibition than our legacy Sanger assay. Overall, the success rate (full single source hypervariable regions I and II (HVI/HVII) for Sanger and control region for MPS result) for these challenging samples increased from 32.3% with Sanger sequencing to 74.2% with the Precision ID assay. Considering the increase in success rate, the simple workflow and the higher discriminating potential of whole genome data, the Precision ID platform is a significant improvement for the CA Department of Justice Missing Persons DNA Program.

Platinum-Quality Mitogenome Haplotypes from United States Populations

A total of 1327 platinum-quality mitochondrial DNA haplotypes from United States (U.S.) populations were generated using a robust, semi-automated next-generation sequencing (NGS) workflow with rigorous quality control (QC). The laboratory workflow involved long-range PCR to minimize the co-amplification of nuclear mitochondrial DNA segments (NUMTs), PCR-free library preparation to reduce amplification bias, and high-coverage Illumina MiSeq sequencing to produce an average per-sample read depth of 1000 × for low-frequency (5%) variant detection. Point heteroplasmies below 10% frequency were confirmed through replicate amplification, and length heteroplasmy was quantitatively assessed using a custom read count analysis tool. Data analysis involved a redundant, dual-analyst review to minimize errors in haplotype reporting with additional QC checks performed by EMPOP. Applying these methods, eight sample sets were processed from five U.S. metapopulations (African American, Caucasian, Hispanic, Asian American, and Native American) corresponding to self-reported identity at the time of sample collection. Population analyses (e.g., haplotype frequencies, random match probabilities, and genetic distance estimates) were performed to evaluate the eight datasets, with over 95% of haplotypes unique per dataset. The platinum-quality mitogenome haplotypes presented in this study will enable forensic statistical calculations and thereby support the usage of mitogenome sequencing in forensic laboratories.

Searching the undetected mtDNA variants in forensic MPS data

The efficiency of MPS in forensic mtDNA analysis has been thoroughly proven, although a reliable and well established data evaluation still remains a critical point. Numerous bioinformatics tools have been developed, but most of them require specific operating systems and high costs, while free open-source programs with user-friendly interfaces are few. In this study, 43 full mtGenomes were sequenced using the Ion Personal Genome Machine™ (PGM™) System and analyzed utilizing the plug-in Variant Caller (TVC) of the Ion Torrent Software Suite and the mtDNA-Server (mDS), a free web-based mitochondrial analysis tool for MPS data. The outcomes of these two different analysis tools were compared to variants noted after manual inspection of the aligned reads performed using Integrative Genomics Viewer (IGV). The comparison highlighted the presence of thirty-nine discordant variant calls, which were resolved by Sanger sequencing that confirmed the presence of all variants, except for 7 deletions. The combined adoption of IGV and Sanger type sequencing confirmatory steps, in addition of TVC and mDS analysis, resulted in a more accurate variants assignment with the detection of 32 additional true polymorphisms, which were noted in the final dataset. Regarding the heteroplasmy issue, out of a total of thirty heteroplasmic variants, twenty-eight were detected by the TVC, while the mDS detected twenty-two. Overall, none of the used bioinformatics tools were the perfect choice and a secondary analysis with an expert's opinion in complete mtGenome MPS data evaluation is still required in forensic genetic analysis.

Pathogenic Variant Filtering for Mitochondrial Genome Haplotype Reporting

Given the enhanced discriminatory power of the mitochondrial DNA (mtDNA) genome (mitogenome) over the commonly sequenced control region (CR) portion, the scientific merit of mitogenome sequencing is generally accepted. However, many laboratories remain beholden to CR sequencing due to privacy policies and legal requirements restricting the use of disease information or coding region (codR) information. In this report, we present an approach to obviate the reporting of sensitive codR data in forensic haplotypes. We consulted the MitoMap database to identify 92 mtDNA codR variants with confirmed pathogenicity. We determined the frequencies of these pathogenic variants in literature-quality and forensic-quality databases to be very low, at 1.2% and 0.36%, respectively. The observed effect of pathogenic variant filtering on random match statistics in 2488 forensic-quality mitogenome haplotypes from four populations was nil. We propose that pathogenic variant filtering should be incorporated into variant calling algorithms for mitogenome haplotype reporting to maximize the discriminatory power of the locus while minimizing the reveal of sensitive genetic information.

Distinguishing mitochondrial DNA and NUMT sequences amplified with the precision ID mtDNA whole genome panel

Nuclear mitochondrial DNA segments (NUMTs) are generated via transfer of portions of the mitochondrial genome into the nuclear genome. Given their common origin, there is the possibility that both the mitochondrial and NUMT segments may co-amplify using the same set of primers. Thus, analysis of the variation of the mitochondrial genome must take into account this co-amplification of mitochondrial and NUMT sequences. The study herein builds on data from the study by Strobl et al. (Strobl et al., 2019), in which multiple point heteroplasmies were called with an "N" to prevent labeling NUMT sequences mimicking mitochondrial heteroplasmy and being interpreted as true mitochondrial in origin sequence variants. Each of these point heteroplasmies was studied in greater detail, both molecularly and bioinformatically, to determine whether NUMT or true mitochondrial DNA variation was present. The bioinformatic and molecular tools available to help distinguish between NUMT and mitochondrial DNA and the effect of NUMT sequences on interpretation were discussed.

Massive parallel sequencing in forensics: advantages, issues, technicalities, and prospects

In the last decade, next-generation sequencing (NGS) technology, alternatively massive parallel sequencing (MPS), was applied to all fields of biological research. Its introduction to the field of forensics was slower, mainly due to lack of accredited sequencers, kits, and relatively higher sequencing error rates as compared with standardized Sanger sequencing. Currently, a majority of the problematic issues have been solved, which is proven by the body of reports in the literature. Here, we discuss the utility of NGS sequencing in forensics, emphasizing the advantages, issues, the technical aspects of the experiments, commercial solutions, and the potentially interesting applications of MPS.

Complete mitogenome data for the Serbian population: the contribution to high-quality forensic databases

Mitochondrial genome (mtDNA) is a valuable resource in resolving various human forensic casework. The usage of variability of complete mtDNA genomes increases their discriminatory power to the maximum and enables ultimate resolution of distinct maternal lineages. However, their wider employment in forensic casework is nowadays limited by the lack of appropriate reference database. In order to fill in the gap in the reference data, which, considering Slavic-speaking populations, currently comprises only mitogenomes of East and West Slavs, we present mitogenome data for 226 Serbians, representatives of South Slavs from the Balkan Peninsula. We found 143 (sub)haplogroups among which West Eurasian ones were dominant. The percentage of unique haplotypes was 85%, and the random match probability was as low as 0.53%. We support previous findings on both high levels of genetic diversity in the Serbian population and patterns of genetic differentiation among this and ten studied European populations. However, our high-resolution data supported more pronounced genetic differentiation among Serbians and two Slavic populations (Russians and Poles) as well as expansion of the Serbian population after the Last Glacial Maximum and during the Migration period (fourth to ninth century A.D.), as inferred from the Bayesian skyline analysis. Phylogenetic analysis of haplotypes found in Serbians contributed towards the improvement of the worldwide mtDNA phylogeny, which is essential for the interpretation of the mtDNA casework.

Joint Genetic Analyses of Mitochondrial and Y-Chromosome Molecular Markers for a Population from Northwest China

The genetic markers on mitochondria DNA (mtDNA) and Y-chromosome can be applied as a powerful tool in population genetics. We present a study to reveal the genetic background of Kyrgyz group, a Chinese ethnic group living in northwest China, and genetic polymorphisms of 60 loci on maternal inherited mtDNA and 24 loci on paternal inherited Y-chromosome short tandem repeats (Y-STRs) were investigated. The relationship between the two systems was tested, and the result indicated that they were statistically independent from each other. The genetic distances between Kyrgyz group and 11 reference populations for mtDNA, and 13 reference populations for Y-STRs were also calculated, respectively. The present results demonstrated that the Kyrgyz group was genetically closer to East Asian populations than European populations based on the mtDNA loci but the other way around for the Y-STRs. The genetic analyses could largely strengthen the understanding for the genetic background of the Kyrgyz group.

The lot-to-lot variability in the mitochondrial genome of controls

Current research in the biomedical field has illustrated how cell lines used as reference standards can change over time and, more importantly, can affect research and diagnostic results obtained from these cell lines. With the use of increasingly sensitive and highly resolving technologies (e.g., massively parallel sequencing), forensic scientists must be aware of and account for potential variability in the cell lines used as controls in their validation studies and day-to-day casework. In this study, multiple lot numbers from four commonly-used control cell line DNAs were sequenced with massively parallel sequencing on the Ion S5. The variability among these different lots was evaluated, and the effect on forensic laboratory work discussed.

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.

Characterization of sequence variation at 30 autosomal STRs in Chinese Han and Tibetan populations

Massively parallel sequencing (MPS) technologies have the ability to reveal sequence variations within STR alleles as well as their nominal allele lengths, which have traditionally been detected by CE instruments. Recently, Thermo Fisher Scientific has updated the MPS-STR panel, named the Precision ID GlobalFiler next-generation sequencing (NGS) STR Panel version 2, with primers redesigned to add two pentanucleotide tandem repeat loci and profile interpretation supported by the Converge software. Using the Ion Chef System, the Ion S5XL System, and the Converge software, genetic variations were characterized within STR repeat and flanking regions of 30 autosomal STR markers in 115 unrelated individuals from two Chinese population groups (58 Tibetans and 57 Hans). Nineteen STRs demonstrated a relative increase in diversity with the variant sequence alleles compared with those of traditional nominal length alleles. In total, 390 alleles were identified by their sequences compared with 258 alleles that were identified by length. Of these 92 sequence variants found within the STR repeat regions, 40 variants were located in STR flanking regions. Additionally, the agreement of the results with CE data was evaluated, as was the ability of this new MPS panel to analyze case-type (11 samples) and artificially degraded samples (seven samples in triplicate). The results generated from this study illustrate that extensive sequence variation exists in commonly used STR markers in the selected population samples and indicate that this NGS STR panel has the potential to be used as an effective tool for human forensics.

Comparison of whole mitochondrial genome variants between hair shafts and reference samples using massively parallel sequencing

Hair shafts are one of the most common types of evidence at crime scenes, and mitochondrial DNA (mtDNA) has been analyzed as a valuable genetic marker for hair shafts in forensic casework. However, the mtDNA analysis strategy may vary according to the quantity and quality of DNA extracted from a forensic sample and the available massively parallel sequencing (MPS) platform in laboratories. Forensic practitioners often have to interpret mtDNA sequences exhibiting point heteroplasmy (PHP) that are analyzed using different analytical methods. In the present study, the whole mitochondrial genome (mtGenome) variants of hair shaft samples obtained from 20 donors, which were sampled in duplicate and stored at room temperature for > 1 year, were analyzed using the Precision ID mtDNA Whole Genome Panel and Ion S5 system. The whole mtGenome variants of 20 blood and 20 buccal swab samples (reference samples) from the hair shaft donors were analyzed using the Nextera XT DNA Library Prep Kit and MiSeq System. A total of 20 unique mtGenome haplotypes were observed, and 56 PHP variants were identified across the 4 sets of tissue. When the major nucleotide of PHP was considered, 16 of 20 haplotypes of the hair shaft samples matched those of the corresponding blood and buccal swab samples. In four donors, the major nucleotide of PHP was inverted at one nucleotide position between the hair shaft and reference samples. However, the data obtained on MPS, showing high PHP resolution, provided substantial information to avoid false exclusion when comparing two haplotypes containing PHP with inverted major nucleotides. In conclusion, the present study demonstrates the utility of MPS in forensic casework in the comparative analysis of mtGenome variants containing PHP.

Estimation of the number of contributors to mixed samples of DNA by mitochondrial DNA analyses using massively parallel sequencing

We evaluated whether the number of contributors to mixed DNA samples can be estimated by analyzing the D-loop of mitochondrial DNA using massively parallel sequencing. The A- (positions 16,209-16,400) and B- (positions 30-284) amplicons in hypervariable regions 1 and 2, respectively, were sequenced using MiSeq with 2 × 251 cycles. Sequence extraction and trimming were performed using CLC Genomics Workbench 11 and the number of observed haplotypes was counted for each amplicon type using Microsoft Excel. The haplotype ratios were calculated by dividing the number of counted reads of the corresponding haplotype by the total number of sequence reads. Haplotypes that were over the threshold (5%) were defined as positive haplotypes. The number of larger positive haplotypes in either of the two amplicon types was defined as the number of contributors. Samples were collected from seven individuals. Seventeen mixed samples were prepared by mixing DNA from two to five contributors at various ratios. The number of contributors was correctly estimated from almost all of the mixed samples containing equal amounts of DNA from two to five people. In mixed samples of two or three people, the minor components were detected down to a ratio of 20:1 or 8:2:1. However, heteroplasmy, base deletions, and sharing of the same haplotypes caused incorrect estimations of the number of contributors. Although this method still has room for improvement, it may be useful for estimating the number of contributors in a mixed sample, as it does not rely on forensic mathematics.

Haplogroup Structure and Genetic Variation Analyses of 60 Mitochondrial DNA Markers in Southern Shaanxi Han Population

Mitochondrial DNA (mtDNA) has been widely employed as one tool for the studies of human migration and phylogenetic evolution owing to the characteristics of its lack of recombination and matrilineal inheritance. In this study, we analyze genetic distributions of 60 mtDNA markers in 126 unrelated individuals of Southern Shaanxi Han population and classify their haplogroups. Genetic distribution comparisons between Southern Shaanxi Han and other populations from different continents are conducted based on the same mtDNA markers. The majority of 60 mtDNA markers are polymorphic in Southern Shaanxi Han population. The most common haplogroups observed in Southern Shaanxi Han population are B5, followed by D5, A, D4e, and N9a1'3. Obtained matching probability for these 60 mtDNA markers indicates that the panel could be used as a valuable tool in forensic caseworks. Results of genetic distances (Fst) and multidimensional scaling analysis show that Southern Shaanxi Han population has relatively close genetic relationships with other Han populations in different regions. In conclusion, the panel comprising 60 mtDNA markers could be utilized for forensic applications in Southern Shaanxi Han population.

Rare human mitochondrial HV lineages spread from the Near East and Caucasus during post-LGM and Neolithic expansions

Of particular significance to human population history in Eurasia are the migratory events that connected the Near East to Europe after the Last Glacial Maximum (LGM). Utilizing 315 HV*(xH,V) mitogenomes, including 27 contemporary lineages first reported here, we found the genetic signatures for distinctive movements out of the Near East and South Caucasus both westward into Europe and eastward into South Asia. The parallel phylogeographies of rare, yet widely distributed HV*(xH,V) subclades reveal a connection between the Italian Peninsula and South Caucasus, resulting from at least two (post-LGM, Neolithic) waves of migration. Many of these subclades originated in a population ancestral to contemporary Armenians and Assyrians. One such subclade, HV1b-152, supports a postexilic, northern Mesopotamian origin for the Ashkenazi HV1b2 lineages. In agreement with ancient DNA findings, our phylogenetic analysis of HV12 and HV14, the two exclusively Asian subclades of HV*(xH,V), point to the migration of lineages originating in Iran to South Asia before and during the Neolithic period. With HV12 being one of the oldest HV subclades, our results support an origin of HV haplogroup in the region defined by Western Iran, Mesopotamia, and the South Caucasus, where the highest prevalence of HV has been found.