Assessing heteroplasmic variant drift in the mtDNA control region of human hairs using an MPS approach

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

Forensic Investigative Genetic Genealogy, a recent sub discipline of forensic genomics, leverages the high throughput and sensitivity of detection of next generation sequencing and established genetic and genealogical approaches to support the identification of human remains from missing persons investigations and investigative lead generation in violent crimes. To facilitate forensic DNA evidence analysis, the ForenSeq® Kintelligence multiplex, consisting of 10,230 SNPs, was developed. Design of the ForenSeq Kintelligence Kit, the MiSeq FGx® Sequencing System and the ForenSeq Universal Analysis Software is described. Developmental validation in accordance with SWGDAM guidelines and forensic quality assurance standards, using single source samples, is reported for the end-to-end workflow from library preparation to data interpretation. Performance metrics support the conclusion that more genetic information can be obtained from challenging samples compared to other commercially available forensic targeted DNA assays developed for capillary electrophoresis (CE) or other current next generation sequencing (NGS) kits due to the higher number of markers, the overall shorter amplicon sizes (97.8% <150 bp), and kit design. Data indicate that the multiplex is robust and fit for purpose for a wide range of quantity and quality samples. The ForenSeq Kintelligence Kit and the Universal Analysis Software allow transfer of the genetic component of forensic investigative genetic genealogy to the operational forensic laboratory.

A New Tool for Probabilistic Assessment of MPS Data Associated with mtDNA Mixtures

Mitochondrial (mt) DNA plays an important role in the fields of forensic and clinical genetics, molecular anthropology, and population genetics, with mixture interpretation being of particular interest in medical and forensic genetics. The high copy number, haploid state (only a single haplotype contributed per individual), high mutation rate, and well-known phylogeny of mtDNA, makes it an attractive marker for mixture deconvolution in damaged and low quantity samples of all types. Given the desire to deconvolute mtDNA mixtures, the goals of this study were to (1) create a new software, MixtureAceMT™, to deconvolute mtDNA mixtures by assessing and combining two existing software tools, MixtureAce™ and Mixemt, (2) create a dataset of in-silico MPS mixtures from whole mitogenome haplotypes representing a diverse set of population groups, and consisting of two and three contributors at different dilution ratios, and (3) since amplicon targeted sequencing is desirable, and is a commonly used approach in forensic laboratories, create biological mixture data associated with two amplification kits: PowerSeq™ Whole Genome Mito (Promega™, Madison, WI, USA) and Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific by AB™, Waltham, MA, USA) to further validate the software for use in forensic laboratories. MixtureAceMT™ provides a user-friendly interface while reducing confounding features such as NUMTs and noise, reducing traditionally prohibitive processing times. The new software was able to detect the correct contributing haplogroups and closely estimate contributor proportions in sequencing data generated from small amplicons for mixtures with minor contributions of ≥5%. A challenge of mixture deconvolution using small amplicon sequencing is the potential generation of spurious haplogroups resulting from private mutations that differ from Phylotree. MixtureAceMT™ was able to resolve these additional haplogroups by including known haplotype/s in the evaluation. In addition, for some samples, the inclusion of known haplotypes was also able to resolve trace contributors (minor contribution 1-2%), which remain challenging to resolve even with deep sequencing.

Assessment of mitochondrial DNA copy number variation relative to nuclear DNA quantity between different tissues

Mitochondrial DNA is a widely tested genetic marker in various fields of research and diagnostics. Nonetheless, there is still little understanding on its abundance and quality within different tissues. Aiming to obtain deeper knowledge about the content and quality of mtDNA, we investigated nine tissues including blood, bone, brain, hair (root and shaft), cardiac muscle, liver, lung, skeletal muscle, and buccal mucosa of 32 deceased individuals using two real-time quantitative PCR-based assays with differently sized mtDNA and nDNA targets. The results revealed that the quantity of nDNA is a weak surrogate to estimate mtDNA quantities among tissues of an individual, as well as tissues across individuals. Especially hair showed extreme variation, depicting a range of multiple magnitudes of mtDNA molecules per hair fragment. Furthermore, degradation can lead to fewer fragments being available for PCR. The results call for parallel determination of the quantity and quality of mtDNA prior to downstream genotyping assays.

CmVCall: An automated and adjustable nanopore analysis pipeline for heteroplasmy detection of the control region in human mitochondrial genome

Genetic associations between human mitochondrial DNA (mtDNA) heteroplasmy and mitochondrial diseases, aging, and cancer have been elaborated, contributing a lot to the further understanding of mtDNA polymorphic spectrum in anthropology, population, and forensic genetics. In the past decade, heteroplasmy detection using Sanger sequencing and next generation sequencing (NGS) was hampered by the former's inefficiency and the latter's inherent bias due to amplification and mapping of short reads, respectively. Nanopore sequencing stands out for its ability to yield long contiguous segments of DNA, providing a new insight into heterogeneity authentication. In addition to MinION from Oxford Nanopore Technologies, an alternative nanopore sequencer QNome (Qitan Technology) has also been applied to various biological research and the forensic applicability of this platform has been proved recently. In this study, we evaluated the performance of four commonly used variant callers in the heterogeneity authentication of the control region of human mtDNA based on simulations of different ratios generated by mixing QNome nanopore sequencing reads of two synthetic sequences. Then, an open-source and python-based nanopore analytics pipeline, CmVCall was developed and incorporated multiple programs including reads filtering, removal of nuclear mitochondrial sequences (NUMTs), alignment, optional 'Correction' mode, and heterogeneity identification. CmVCall can achieve high precision, accuracy, and recall of 100%, 99.9%, and 92.3% with a 5% heteroplasmy level in 'Correction' mode. Moreover, blood, saliva, and hair shaft samples from monozygotic (MZ) twins were used for heterogeneity evaluation and comparison with the NGS data. Results of MZ twin samples showed that CmVCall could identify more point heteroplasmy sites, revealing significant levels of inter- and intra-individual mtDNA polymorphism. In conclusion, we believe that this analysis pipeline will lay a solid foundation for the development of a comprehensive nanopore analysis pipeline targeting the whole mitochondrial genome.

Massively Parallel Sequencing of the Mitogenome from Human Hair Shafts in Forensic Investigations

This article highlights methods used to perform DNA extraction, mitochondrial DNA quantification, multiplex PCR amplification, amplicon-based massively parallel sequencing, and data analysis of the mitochondrial genome (mitogenome) from human hair shafts. The focus is on applications to forensic casework, but this set of protocols can be used for any purpose involving small cuttings (as small as 1 to 5 mm) of human hair shafts up to 40 years from the time of collection. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Extraction of mitochondrial DNA from human hair shafts Basic Protocol 2: Quantification of mitochondrial DNA (copies/μl) Basic Protocol 3: Multiplex amplification of the mitogenome Basic Protocol 4: Library preparation and sequencing of mitogenome amplicons Basic Protocol 5: Data analysis of mitogenome haplotypes.

Impact of storage conditions and time on DNA yield from ammunition cartridges

Recovery of suitable amounts of DNA from ammunition cartridges for short tandem repeat (STR) or mitochondrial (mt) DNA analysis has been a challenge for crime laboratories. The metal composition of cartridge cases and projectiles exposes the DNA to harmful ions that damage and ultimately degrade the DNA such that it cannot be effectively amplified. The current study assessed the impact of time and storage conditions on touch DNA deposited on cartridge components of varying metal content: aluminum, nickel, brass, and copper. Elevated humidity levels facilitated greater DNA degradation and loss compared to low humidity (or "dry") conditions, indicating that recovered cartridge component evidence should be stored in a low-humidity environment immediately after collection, preferably with a desiccant. As expected, a relationship was observed between the amount of time elapsed since the cartridge components were handled and the associated DNA yield. Interestingly, while yields dropped considerably in the first 48-96 h post-handling, regardless of the storage conditions, a layering effect was observed that helps maintain a relatively constant level of surface DNA over extended periods of time. An apparent layering effect was also observed on cartridge components following multiple surface depositions, where yields were two times higher than single deposition samples at similar timepoints. Overall, these findings suggest that storage conditions and a layering affect play an important role in the preservation of DNA on ammunition components.

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.

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.

Whole mitochondrial genome analysis of Tai-Kadai-speaking populations in Southwest China

As a single matrilineal gene, human mitochondrial DNA plays a very important role in the study of population genetics. The whole mitogenome sequences of 287 individuals of the Tai-Kadai-speaking population in Guizhou were obtained. It was discovered that there were 82, 104, and 94 haplotypes in 83 Bouyei individuals, 107 Dong individuals, and 97 Sui individuals, respectively; and the haplotype diversity in Bouyei, Dong, and Sui groups was 1.000 ± 0.02, 0.9993 ± 0.0015, and 0.999 ± 0.002, respectively. The result of neutrality tests of the Tai-Kadai-speaking population in Guizhou showed significant negative values, and the analysis of mismatch distribution showed an obvious unimodal distribution. The results implied that Guizhou Tai-Kadai-speaking populations had high genetic diversities and may have experienced recent population expansion. In addition, the primary haplogroups of studied populations were M*, F, B, D, and R*, implying that they may origin from Southern China. The matrilineal genetic structure of the Tai-Kadai-speaking populations in Guizhou was analyzed by merging the mitogenome data of 79 worldwide populations as reference data. The results showed that there were close relationships between studied populations and other Tai-Kadai as well as some Austronesian populations in East and Southeast Asia. Overall, the mitogenome data generated in this study will provide important data for the study of genetic structure of Tai-Kadai speaking populations.

Genetic Background of Kirgiz Ethnic Group From Northwest China Revealed by Mitochondrial DNA Control Region Sequences on Massively Parallel Sequencing

The mitochondrial DNA (mtDNA) has been used to trace population evolution and apply to forensic identification due to the characteristics including lack of recombination, higher copy number and matrilineal inheritance comparing with nuclear genome DNA. In this study, mtDNA control region sequences of 91 Kirgiz individuals from the Northwest region of China were sequenced to identify genetic polymorphisms and gain insight into the genetic background of the Kirgiz ethnic group. MtDNA control region sequences of Kirgiz individuals presented relatively high genetic polymorphisms. The 1,122 bp sequences of mtDNA control region could differ among unrelated Kirgiz individuals, which suggested the mtDNA control region sequences have a good maternal pedigree tracing capability among different Kirgiz individuals. The neutrality test, mismatch distribution, Bayesian phylogenetic inference, Bayesian skyline analysis, and the median network analyses showed that the Kirgiz group might occurred population expansion, and the expansion could be observed at about ∼53.41 kilo years ago (kya) when ancestries of modern humans began to thrive in Eurasia. The pairwise population comparisons, principal component analyses, and median network analyses were performed based on haplogroup frequencies or mtDNA control region sequences of 5,886 individuals from the Kirgiz group and the 48 reference populations all over the world. And the most homologous haplotypes were found between Kirgiz individuals and the East Asian individuals, which indicated that the Kirgiz group might have gene exchanges with the East Asian populations.

Exploring statistical weight estimates for mitochondrial DNA matches involving heteroplasmy

Massively parallel sequencing (MPS) of mitochondrial (mt) DNA allows forensic laboratories to report heteroplasmy on a routine basis. Statistical approaches will be needed to determine the relative frequency of observing an mtDNA haplotype when including the presence of a heteroplasmic site. Here, we examined 1301 control region (CR) sequences, collected from individuals in four major population groups (European, African, Asian, and Latino), and covering 24 geographically distributed haplogroups, to assess the rates of point heteroplasmy (PHP) on an individual and nucleotide position (np) basis. With a minor allele frequency (MAF) threshold of 2%, the data was similar across population groups, with an overall PHP rate of 37.7%, and the majority of heteroplasmic individuals (77.3%) having only one site of heteroplasmy. The majority (75.2%) of identified PHPs had an MAF of 2-10%, and were observed at 12.6% of the nps across the CR. Both the broad and phylogenetic testing suggested that in many cases the low number of observations of heteroplasmy at any one np results in a lack of statistical association. The posterior frequency estimates, which skew conservative to a degree depending on the sample size in a given haplogroup, had a mean of 0.152 (SD 0.134) and ranged from 0.031 to 0.83. As expected, posterior frequency estimates decreased in accordance with 1/n as the sample size (n) increased. This provides a proposed conservative statistical framework for assessing haplotype/heteroplasmy matches when applying an MPS technique in forensic cases and will allow for continual refinement as more data is generated, both within the CR and across the mitochondrial genome.

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.

Reducing the Number of Mismatches between Hairs and Buccal References When Analysing mtDNA Heteroplasmic Variation by Massively Parallel Sequencing

In forensics, mitochondrial DNA (mtDNA) analysis is foremost applied to rootless hairs often lacking detectable nuclear DNA. Sanger sequencing is the routine mtDNA method in most forensic laboratories, even though interpretation of mixed samples and heteroplasmic sites can be challenging. Individuals may hold cells with low-level heteroplasmy variants below the detection threshold and other cells where this minor variant is the major one. This difference may be interpreted as a mismatch between reference and evidentiary trace samples, such as buccal specimens and rootless hairs. Such mismatches may be solved by Massively Parallel Sequencing (MPS), allowing more sensitive quantitative analysis for mixed positions than Sanger. The mtDNA control region was analysed in buccal reference samples from 26 individuals and 475 corresponding hairs by MPS and compared to Sanger sequencing data generated on the same samples. With MPS, mixed contributions down to 3% were regarded, leading to a substantial increase in the frequency of heteroplasmy. Our results demonstrate that previously reported mismatches between buccal reference and hair shaft samples by Sanger are detected as low-level heteroplasmy by MPS. A detailed overview of buccal and hair heteroplasmy is provided and implications for MPS-based mtDNA analysis in the context of forensic cases are discussed.

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.

Damage patterns observed in mtDNA control region MPS data for a range of template concentrations and when using different amplification approaches

Massively parallel sequencing (MPS) of mitochondrial (mt) DNA allows practitioners the ability to fully resolve heteroplasmic sites. In forensic DNA analysis, identifying heteroplasmy (a naturally occurring mixture of two mtDNA profiles) can provide additional mtDNA profile information which can lead to an increase in the discrimination potential of an mtDNA match between an evidentiary sample and reference source. Forensic samples such as hair and skeletal remains, especially older, more compromised samples, can often exhibit DNA damage. Because both damage and heteroplasmy can manifest as a mixture of two nucleotides, it is important to differentiate between the two conditions when interpreting mtDNA MPS data. In this study, DNA damage was applied under controlled conditions to samples containing a range of template concentrations, including some with identified heteroplasmy. Damage was applied via storage in water at room temperature on samples diluted before or after storage to mimic low template scenarios. Damage was assessed with respect to the following areas: mtDNA quantification and degradation ratios, MPS read depth, MPS profile results, overall damage rates, and the interpretation of heteroplasmy. Datasets were generated to assess and compare two different amplification and library preparation strategies: the Promega PowerSeq™ CRM Nested System kit and a 1.16 kb target amplicon of the entire mtDNA control region followed by a Nextera® XT library preparation. The results of this study provide an evaluation of the Promega 10-plex MPS procedure as an improved process to mitigate the impact of mtDNA damage on low template samples. Some of the negative effects of damage observed in this study were a decrease in mtDNA yield by 20-30% and lower quality MPS sequencing results. These effects were observed more frequently when samples were diluted prior to inducing damage, illustrating that low template samples are more susceptible to damage. The findings of this study will assist forensic laboratories in differentiating between damage and heteroplasmy, which is essential when developing robust mtDNA MPS interpretation guidelines such as setting appropriate reporting thresholds.

A Forensic Genomics Approach for the Identification of Sister Marija Crucifiksa Kozulić

Sister Marija Krucifiksa Kozulić (1852–1922) was a Croatian nun who is in consideration for beatification by the Vatican, which is facilitated by the identification of her 20th-century remains. Sister Marija was buried in a tomb in Rijeka, Croatia, along with other nuns including her biological sister, Tereza Kozulić (1861–1933). When the remains were exhumed in 2011, they were found in a deteriorated state and commingled with several other sets of remains. Thus, mitochondrial genome sequencing of the long bones was performed to sort the remains by mitochondrial haplotype. Two similar but unique haplotypes belonging to haplogroup H1bu were identified, and samples from these bones were subjected to autosomal short tandem repeat (STR) and single nucleotide polymorphism (SNP) sequencing. Although only partial profiles were obtained, the data were sufficient for kinship analysis with the profile of a paternal niece of Sister Marija (Fides Kozulić). The data indicate that it is 574,195-fold more likely that the two sets of skeletal remains represent 2nd-degree relatives of Fides than sisters who are unrelated to Fides. Although it is impossible to discern which set of remains belongs to Marija and which belongs to Tereza, forensic genomics methods have enabled identification of the sisters.

Characterization of background noise in MiSeq MPS data when sequencing human mitochondrial DNA from various sample sources and library preparation methods

Improved resolution of massively parallel sequencing (MPS) allows for the characterization of mitochondrial (mt) DNA heteroplasmy to levels previously unattainable with traditional sequencing approaches. An essential criterion for the reporting of heteroplasmy is the ability of the MPS method to distinguish minor sequence variants (MSVs) from system noise, or error. Therefore, an assessment of the background noise in the MPS method is desirable to identify the point at which reliable data can be reported. Substitution and sequence specific error (SSE) was evaluated for a variety of sample types and two library preparations. Substitution error rates ranged from 0.18 to 0.49 per 100 nucleotides with C positions generally having the highest rate of misincorporation. Comparison of error rates across sample types indicated a significant increase for samples with damaged DNA. The positions of error were varied across datasets (pairwise concordance 0-68%), but had greater consistency within the damaged samples (80-96%). The most commonly observed motif preceding error in forward reads was CCG, while GGT was most common in reverse reads, both consistent with previous findings. The findings illustrate that for datasets containing samples with damaged DNA, reporting thresholds for heteroplasmy may have to be modified and individual sites with error levels exceeding thresholds should be scrutinized. Collectively, the shifting error profiles observed across the various sample types and library preparation methods demonstrates the need for an assessment of error under these varying circumstances. Characterization of the applicable background noise will help to ensure that thresholds are reliably set for detection of true MSVs.

Get it off, but keep it: Efficient cleaning of hair shafts with parallel DNA extraction of the surface stain

The analysis of hair samples is a common task in forensic investigations. Material transferred to the surface of a hair during a crime challenges the analysis as it has to be removed efficiently. However, the removal of the stain can also lead to a loss of information on stain contributors. DNA analysis of the stain itself might thus be helpful for the forensic investigation. The aim of this study was the examination of different methods to remove common biological surface stains completely from human hair shafts without hampering the parallel DNA extraction of the cleaned hair shaft and the isolated surface stain (blood, saliva, vaginal secretion, semen, and skin flocks). Four different methods of cleaning (water, lysis buffer, swabbing, NaClO) were compared to their cleaning efficiency as well as their success of mtDNA analysis of three hair donors and the original five stains on the hair. In order to test the suitability of this procedure for future analysis methods, a selection of samples were also sequenced with MPS. Additionally, nuclear DNA analysis of the stain DNA was performed using a screening STR assay to test the potential success for detection of a STR profile. The most efficient removal of the stain was achieved using NaClO, however compromising further analysis of the stain DNA. The best results for cleaning and parallel stain analysis were obtained using a swab moistened with 0.5 % SDS for surface cleaning. Especially water failed to remove stains efficiently, leading to a high amount of mixed mtDNA in the DNA extracts. MPS showed an increased sensitivity for detection of minute mixtures.

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.

Whole mitochondrial genome analysis of highland Tibetan ethnicity using massively parallel sequencing

Mitochondrial DNA (mtDNA) is a key player in numerous multifaceted and intricate biological processes and plays a pivotal role in dissecting the peopling of different populations, due to its maternally inherited property and comparatively high mutation rate. In this study, 119 Tibetan individuals from the Muli Tibetan Autonomous County of China (average altitude above 3,000 m) were employed in mitochondrial genome (mitogenome) sequencing by massively parallel sequencing (MPS) techniques using the Precision ID mtDNA Whole Genome Panel on an Ion S5XL system. The dataset presented 88 distinct haplotypes, resulting in the haplotype diversity of 0.9909. The majority of haplotypes were assigned to East Asian lineages and the distribution of haplogroups of Muli Tibetan significantly differed from reference Tibetan populations. The maximum parsimony phylogeny reconstructed by 119 newly generated mitogenomes revealed 12 major Muli Tibetan lineages. Intriguingly, a Sherpa-specific sub-haplogroup A15c1 with the lack of mutations at 4216 and 15,924 was discerned in our dataset, which suggested that the maternal gene pool of Sherpas may derive from Tibetan populations. The shared haplogroups between Muli Tibetan and lowland Han Chinese hinted that these lineages may derive from non-Tibetans and have already differentiated before their arrival on the Tibetan Plateau. Furthermore, extensive pairwise population comparisons displayed that Muli Tibetan had a closer genetic relationship with ethnically or linguistically close Nyingtri Tibetan, Nyingtri Lhoba and Chamdo Tibetan populations. Genetic affinity was also observed between the Muli Tibetan and North Han Chinese. Collectively, the results generated in this study enriched the existing forensic mtDNA database and raised additional interest in the application of whole mitogenome sequencing in forensic investigations.

Unraveling heteroplasmy patterns with NOVOPlasty

Heteroplasmy, the existence of multiple mitochondrial haplotypes within an individual, has been studied across different scientific fields. Mitochondrial genome polymorphisms have been linked to multiple severe disorders and are of interest to evolutionary studies and forensic science. Before the development of massive parallel sequencing (MPS), most studies of mitochondrial genome variation were limited to short fragments and to heteroplasmic variants associated with a relatively high frequency (>10%). By utilizing ultra-deep sequencing, it has now become possible to uncover previously undiscovered patterns of intra-individual polymorphisms. Despite these technological advances, it is still challenging to determine the origin of the observed intra-individual polymorphisms. We therefore developed a new method that not only detects intra-individual polymorphisms within mitochondrial and chloroplast genomes more accurately, but also looks for linkage among polymorphic sites by assembling the sequence around each detected polymorphic site. Our benchmark study shows that this method is capable of detecting heteroplasmy more accurately than any method previously available and is the first tool that is able to completely or partially reconstruct the sequence for each mitochondrial haplotype (allele). The method is implemented in our open source software NOVOPlasty that can be downloaded at https://github.com/ndierckx/NOVOPlasty.

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.

Impact of DNA degradation on massively parallel sequencing-based autosomal STR, iiSNP, and mitochondrial DNA typing systems

Biological samples, including skeletal remains exposed to environmental insults for extended periods of time, exhibit increasing levels of DNA damage and fragmentation. Human forensic identification methods typically use a combination of mitochondrial (mt) DNA sequencing and short tandem repeat (STR) analysis, which target segments of DNA ranging from 80 to 500 base pairs (bps). Larger templates are often unavailable as skeletal samples age and the associated DNA degrades. Single-nucleotide polymorphism (SNP) loci target shorter templates and may serve as a solution to the problem. Recently developed assays for STR and SNP analysis using a massively parallel sequencing approach, such as the ForenSeq kit (Verogen, San Diego, CA), offer a means for generating results from degraded samples as they target templates down to 60 to 170 bps. We performed a modeling study that demonstrates that SNPs can increase the significance of an identification when analyzing DNA down to an average size of 100 bps for input amounts between 0.375 and 1 ng of nuclear DNA. Observations from this study were then compared with human skeletal material results (n = 14, ninth to eighteenth centuries), which further demonstrated the utility of the ForenSeq kit for degraded samples. The robustness of the Promega PowerSeq™ Mito System was also tested with human skeletal remains (n = 70, ninth to eighteenth centuries), resulting in successful coverage of 99.29% of the mtDNA control region at 50× coverage or more. This was accompanied by modifications to a mainstream DNA extraction technique for skeletal remains that improved recovery of shorter templates.

Discrimination of haplotype in mitochondrial DNA mixtures using LNA-mediated PCR clamping

Locked nucleic acid (LNA) has been widely used for various genetic analyses, and has many benefits, in terms of the specificity or sensitivity of amplification, because LNA-containing primers/probes form more stable duplexes with template DNA than probes lacking LNA. Here, we developed a new method for discriminating HV1 haplotypes from mitochondrial DNA (mtDNA) mixtures by applying PCR clamping using LNA. PCR clamping is based on the selective inhibition of amplification using LNA-containing probes, which can discriminate single-nucleotide differences. Before designing probes, we selected 171 sequences with single-nucleotide variations from the HV1 region, and evaluated the specificity of LNA-containing probes for them by predicting Tm values. The differences of Tm between mismatched and exactly matched probe-template duplexes depended markedly on the type of LNA nucleotides for discriminating single-nucleotide differences, and the cytosine LNA nucleotide at the site of variations in the probes was most effective to discriminate these differences. For mixture analysis, each probe targeted one or two variations (16209C, 16217C, 16257A/16261T, 16297C/16298C, 16304C, 16362C, or 16362T) that are particularly common in the Japanese population, and seven designed probes completely inhibited the amplification of exactly matched templates. We prepared mixed samples by mixing DNA from two individuals at a ratio of 1:9, 1:4, 1:1, 4:1, or 9:1, and then performed Sanger sequencing analysis after PCR clamping with each probe. Our method distinguished each haplotype at lower ratios from two-person mixtures, and enabled sensitive detection at 12 pg of total DNA including 600 copies of mtDNA. Moreover, we analyzed three-person mixtures with representative sequences, and detected the minor haplotype of one individual present at a rate of 10% by adding two selected probes. The ability to discriminate haplotypes in mixed samples by using LNA-mediated PCR clamping indicates the potential value of mtDNA analysis in criminal investigations.

Mitigating the effects of reference sequence bias in single-multiplex massively parallel sequencing of the mitochondrial DNA control region

•

mtDNA control region of 101 diverse samples amplified in a single reaction as 10 overlapping amplicons and sequenced via MPS.

•

Primers create reference bias, compromising ability to call variants or heteroplasmy in primer-binding regions.

•

Bioinformatic selection of overarching reads bypasses effects of proprietary primers and mitigates bias.

•

Data processing permits accurate calling of variants, and heteroplasmies down to 5% level.

Sequence analysis of the mitochondrial DNA (mtDNA) control region can provide forensically useful information, particularly in challenging samples where autosomal DNA profiling fails. Sub-division of the 1122-bp region into shorter PCR fragments improves data recovery, and such fragments can be analysed together via massively parallel sequencing (MPS). Here, we generate mtDNA data using the prototype PowerSeq™ Auto/Mito/Y System (Promega) MPS assay, in which a single PCR reaction amplifies ten overlapping amplicons of the control region, in a set of 101 highly diverse samples representing most major clades of the mtDNA phylogeny. The overlapping multiplex design leads to non-uniform coverage in the regions of overlap, where it is further increased by short amplicons generated alongside the intended products. Primer sequences in targeted amplification libraries are a potential source of reference sequence bias and thus should be removed, but the proprietary nature of the primers in commercial kits necessitates an alternative approach that minimises data loss: here, we introduce the bioinformatic selection of sequencing reads spanning putative primer sites (Overarching Read Enrichment Option, OREO). While OREO performs well in mitigating the effects of primer sequences at the ends of sequence reads, we still find evidence of the internalisation of primer-derived sequences by overlap extension, which may compromise the ability to call variants or to measure heteroplasmy in primer-binding regions. The commercially available PowerSeq™ CRM Nested System design prevents primer internalisation, as shown in a reanalysis of a subset of 57 samples that contain possible heteroplasmies. In combination with OREO, the CRM Nested kit mitigates reference sequence bias, allowing heteroplasmic variants to be estimated down to a 5% threshold. Provided appropriate steps are taken in data processing, single-reaction multiplex assays represent robust tools to analyse mtDNA control region variation. The OREO approach will allow users to bypass the effects of unknown primer sequences in any single-reaction tiled multiplex and eliminate primer-derived bias in overlapping amplicon sequencing studies, in both forensic and non-forensic settings.

Recovery of mtDNA from unfired metallic ammunition components with an assessment of sequence profile quality and DNA damage through MPS analysis

Recovery of suitable amounts of quality DNA from copper and brass surfaces, like those encountered in ammunition, has been a challenge for the forensic community. The ability of copper ions to rapidly facilitate oxidative damage leading to fragmentation of DNA significantly reduces the pool of templates for PCR amplification. We compared two methods for recovering mitochondrial (mt) DNA from the surface of unfired copper projectiles, brass casings, and aluminum casings, and found that using a cotton swab moistened with 0.5M EDTA was the favored approach, especially when the metallic surface was etched. Degradation was significantly higher for DNA samples recovered from copper and brass surfaces, when compared to aluminum. Massively parallel sequencing (MPS) of the control region, using the PowerSeq™ CRM Nested System kit and the Illumina MiSeq instrument, produced full haplotypes for aluminum samples regardless of the method used to deposit or collect DNA, while less than 60% of the copper and brass samples produced partial or full profile information. Touch DNA collected from copper and brass samples produced higher rates of partial or full MPS profile information (∼88-96%), while collection with 0.5M EDTA produced better results than when collection was performed with water; average of ∼70% versus ∼47%. While MPS data was not impacted by noise in the sequencing process, a higher than expected rate of noise was observed, potentially due to an increase in low-level damage lesions. Noise patterns were strikingly different when compared to control data, suggesting that noisy sites may be predictable when testing samples with high levels of oxidative damage. Library preparation was a poor predictor of MPS data quality, as a large percentage of reads did not align with the reference genome. This may impact the number of samples that can be run when a deep-coverage MPS approach is being considered for analysis of mtDNA heteroplasmy. Overall, when applying an MPS approach to the analysis of mtDNA recovered from ammunition, results are expected from touch DNA, will be limited for copper and brass components when the DNA is exposed to an aqueous environment, and DNA degradation will be accelerated when DNA comes in contact with copper or brass surfaces. Practitioners should consider collecting DNA from metallic surfaces with 0.5M EDTA, as this will maximize yield and mitigate degradation. The results of this study directly impact MPS analysis of minor mtDNA sequence variants from metallic surfaces, and are particularly relevant to forensic investigations.

Fragmented Nuclear DNA is the Predominant Genetic Material in Human Hair Shafts

While shed hairs are one of the most commonly encountered evidence types, they are among the most limited in terms of DNA quantity and quality. As a result, nuclear DNA short tandem repeat (STR) profiling is generally unsuccessful and DNA testing of shed hair is instead performed by targeting the mitochondrial DNA control region. Although the high copy number of mitochondrial DNA relative to nuclear DNA routinely permits the recovery of mitochondrial DNA (mtDNA) data in these cases, mtDNA profiles do not offer the discriminatory power of nuclear DNA profiles. In order to better understand the total content and degradation state of DNA in single shed hairs and assess the feasibility of recovering highly discriminatory nuclear DNA data from this common evidence type, high throughput shotgun sequencing was performed on both recently collected and aged (approximately 50-year-old) hair samples. The data reflect trends that have been demonstrated previously with other technologies, namely that mtDNA quantity and quality decrease along the length of the hair shaft. In addition, the shotgun data reveal that nuclear DNA is present in shed hair and surprisingly abundant relative to mitochondrial DNA, even in the most distal fragments. Nuclear DNA comprised, at minimum, 88% of the total human reads in any given sample, and generally more than 95%. Here, we characterize both the nuclear and mitochondrial DNA content of shed hairs and discuss the implications of these data for forensic investigations.

Massively parallel sequencing techniques for forensics: A review

DNA sequencing, starting with Sanger's chain termination method in 1977 and evolving into the next generation sequencing (NGS) techniques of today that employ massively parallel sequencing (MPS), has become essential in application areas such as biotechnology, virology, and medical diagnostics. Reflected by the growing number of articles published over the last 2-3 years, these techniques have also gained attention in the forensic field. This review contains a brief description of first, second, and third generation sequencing techniques, and focuses on the recent developments in human DNA analysis applicable in the forensic field. Relevance to the forensic analysis is that besides generation of standard STR-profiles, DNA repeats can also be sequenced to look for polymorphisms. Furthermore, additional SNPs can be sequenced to acquire information on ancestry, paternity or phenotype. The current MPS systems are also very helpful in cases where only a limited amount of DNA or highly degraded DNA has been secured from a crime scene. If enough autosomal DNA is not present, mitochondrial DNA can be sequenced for maternal lineage analysis. These developments clearly demonstrate that the use of NGS will grow into an indispensable tool for forensic science.

Evaluation of the precision ID whole MtDNA genome panel for forensic analyses

Mitochondrial DNA (mtDNA) amplification and Massively Parallel Sequencing (MPS) using an early access version of the Precision ID Whole MtDNA Genome Panel (Thermo Fisher Scientific) and the Ion Personal Genome Machine (PGM) were evaluated using 15 forensically relevant samples. Samples were selected to represent typical forensic specimens for mtDNA analysis including hairs, hair shafts, swabs and ancient solid tissue samples (bones and teeth) that were stored in the freezer for up to several years after having been typed with conventional Sanger-type Sequencing and Capillary Electrophoresis. The MPS haplotypes confirmed the earlier results in all samples and provided additional sequence information that improved discrimination power and haplogroup estimation. The results raised the appetite for further experiments to validate and apply the new technology in forensic practice.

Deep-Coverage MPS Analysis of Heteroplasmic Variants within the mtGenome Allows for Frequent Differentiation of Maternal Relatives

Distinguishing between maternal relatives through mitochondrial (mt) DNA sequence analysis has been a longstanding desire of the forensic community. Using a deep-coverage, massively parallel sequencing (DCMPS) approach, we studied the pattern of mtDNA heteroplasmy across the mtgenomes of 39 mother-child pairs of European decent; haplogroups H, J, K, R, T, U, and X. Both shared and differentiating heteroplasmy were observed on a frequent basis in these closely related maternal relatives, with the minor variant often presented as 2–10% of the sequencing reads. A total of 17 pairs exhibited differentiating heteroplasmy (44%), with the majority of sites (76%, 16 of 21) occurring in the coding region, further illustrating the value of conducting sequence analysis on the entire mtgenome. A number of the sites of differentiating heteroplasmy resulted in non-synonymous changes in protein sequence (5 of 21), and to changes in transfer or ribosomal RNA sequences (5 of 21), highlighting the potentially deleterious nature of these heteroplasmic states. Shared heteroplasmy was observed in 12 of the 39 mother-child pairs (31%), with no duplicate sites of either differentiating or shared heteroplasmy observed; a single nucleotide position (16093) was duplicated between the data sets. Finally, rates of heteroplasmy in blood and buccal cells were compared, as it is known that rates can vary across tissue types, with similar observations in the current study. Our data support the view that differentiating heteroplasmy across the mtgenome can be used to frequently distinguish maternal relatives, and could be of interest to both the medical genetics and forensic communities.