Evidence for multi-copy Mega-NUMTs in the human genome.

The maternal mode of mitochondrial DNA (mtDNA) inheritance is central to human genetics. Recently, evidence for bi-parental inheritance of mtDNA was claimed for individuals of three pedigrees that suffered mitochondrial disorders. We sequenced mtDNA using both direct Sanger and Massively Parallel Sequencing in several tissues of eleven maternally related and other affiliated healthy individuals of a family pedigree and observed mixed mitotypes in eight individuals. Cells without nuclear DNA, i.e. thrombocytes and hair shafts, only showed the mitotype of haplogroup (hg) V. Skin biopsies were prepared to generate ρ° cells void of mtDNA, sequencing of which resulted in a hg U4c1 mitotype. The position of the Mega-NUMT sequence was determined by fluorescence in situ hybridization and two different quantitative PCR assays were used to determine the number of contributing mtDNA copies. Thus, evidence for the presence of repetitive, full mitogenome Mega-NUMTs matching haplogroup U4c1 in various tissues of eight maternally related individuals was provided. Multi-copy Mega-NUMTs mimic mixtures of mtDNA that cannot be experimentally avoided and thus may appear in diverse fields of mtDNA research and diagnostics. We demonstrate that hair shaft mtDNA sequencing provides a simple but reliable approach to exclude NUMTs as source of misleading results.

Mitochondrial control region variation in a Korean population sample.

Mitochondrial DNA (mtDNA) control region (16024-576) sequences were generated from 281 individuals from South Korea. Robotic liquid handling, a redundant sequencing strategy, and a series of quality control checks were implemented to ensure the high quality of the dataset. This population sample showed a low random match probability (0.25 %) and high genetic diversity (0.9933). The haplogroup breakdown was consistent with previous studies describing Korean mtDNA variation. The 224 unique haplotypes (33 shared) presented will supplement the data already publically available.

A high-throughput Sanger strategy for human mitochondrial genome sequencing.

BACKGROUND: A population reference database of complete human mitochondrial genome (mtGenome) sequences is needed to enable the use of mitochondrial DNA (mtDNA) coding region data in forensic casework applications. However, the development of entire mtGenome haplotypes to forensic data quality standards is difficult and laborious. A Sanger-based amplification and sequencing strategy that is designed for automated processing, yet routinely produces high quality sequences, is needed to facilitate high-volume production of these mtGenome data sets. RESULTS: We developed a robust 8-amplicon Sanger sequencing strategy that regularly produces complete, forensic-quality mtGenome haplotypes in the first pass of data generation. The protocol works equally well on samples representing diverse mtDNA haplogroups and DNA input quantities ranging from 50 pg to 1 ng, and can be applied to specimens of varying DNA quality. The complete workflow was specifically designed for implementation on robotic instrumentation, which increases throughput and reduces both the opportunities for error inherent to manual processing and the cost of generating full mtGenome sequences. CONCLUSIONS: The described strategy will assist efforts to generate complete mtGenome haplotypes which meet the highest data quality expectations for forensic genetic and other applications. Additionally, high-quality data produced using this protocol can be used to assess mtDNA data developed using newer technologies and chemistries. Further, the amplification strategy can be used to enrich for mtDNA as a first step in sample preparation for targeted next-generation sequencing.

Characterization of mitochondrial DNA control region lineages in Iraq.

To evaluate the utility of mtDNA control region data for the purposes of forensic DNA testing in Iraq, a sample of 182 subjects (128 Arab Muslims, 15 Kurd Muslims, 22 Assyrian Christians and 17 Mandaean Arabs) was tested. High numbers of singleton haplotypes were observed among Arabs, Kurds and Assyrians, but fewer were found in Mandaeans. High molecular diversity and low random match probabilities confirmed the value of control region data in the investigation of maternal genetic lineages among the Iraqi population.

Mitochondrial DNA control region variation from samples of the Moroccan population.

In an effort to facilitate forensic mitochondrial DNA (mtDNA) testing in Morocco, high-quality control region sequences from 509 individuals were generated using a comprehensive processing and data review system. This large dataset of random samples from various Moroccan population groups (Arab speaking, Berber speaking, and Sahrawi speaking) exhibited a low random match probability (0.52 %) and a mean of pairwise comparisons of 13.24. The Moroccan mtDNA gene pool studied here was defined entirely by West Eurasian (58.15 %) and African haplogroups (41.85 %).

Improved DNA Extraction and Illumina Sequencing of DNA Recovered from Aged Rootless Hair Shafts Found in Relics Associated with the Romanov Family.

This study describes an optimized DNA extraction protocol targeting ultrashort DNA molecules from single rootless hairs. It was applied to the oldest samples available to us: locks of hairs that were found in relics associated with the Romanov family. Published mitochondrial DNA genome sequences of Tsar Nicholas II and his wife, Tsarina Alexandra, made these samples ideal to assess this DNA extraction protocol and evaluate the types of genetic information that can be recovered by sequencing ultrashort fragments. Using this method, the mtGenome of the Tsarina's lineage was identified in hairs that were concealed in a pendant made by Karl Fabergé for Alexandra Feodorovna Romanov. In addition, to determine if the lock originated from more than one individual, two hairs from the locket were extracted independently and converted into Illumina libraries for shotgun sequencing on a NextSeq 500 platform. From these data, autosomal SNPs were analyzed to assess relatedness. The results indicated that the two hairs came from a single individual. Genetic testing of hairs that were found in the second artifact, a framed photograph of Louise of Hesse-Kassel, Queen of Denmark and maternal grandmother of Tsar Nicholas II, revealed that the hair belonged to a woman who shared Tsar Nicholas' maternal lineage, including the well-known point heteroplasmy at position 16169.

The mtDNA composition of Uzbekistan: a microcosm of Central Asian patterns.

In order to better characterize and understand the mtDNA population genetics of Central Asia, the mtDNA control regions of over 1,500 individuals from Uzbekistan have been sequenced. Although all samples were obtained from individuals residing in Uzbekistan, individuals with direct ancestry from neighboring Central Asian countries are included. Individuals of Uzbek ancestry represent five distinct geographic regions of Uzbekistan: Fergana, Karakalpakstan, Khorezm, Qashkadarya, and Tashkent. Individuals with direct ancestry in nearby countries originate from Kazakhstan, Kyrgyzstan, Russia, Afghanistan, Turkmenistan, and Tajikistan. Our data reinforce the evidence of distinct clinal patterns that have been described among Central Asian populations with classical, mtDNA, and Y-chromosomal markers. Our data also reveal hallmarks of recent demographic events. Despite their current close geographic proximity, the populations with ancestry in neighboring countries show little sign of admixture and retain the primary mtDNA patterns of their source populations. The genetic distances and haplogroup distributions among the ethnic populations are more indicative of a broad east-west cline among their source populations than of their relatively small geographic distances from one another in Uzbekistan. Given the significant mtDNA heterogeneity detected, our results emphasize the need for heightened caution in the forensic interpretation of mtDNA data in regions as historically rich and genetically diverse as Central Asia.

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.

Investigation of heteroplasmy in the human mitochondrial DNA control region: a synthesis of observations from more than 5000 global population samples.

Instances of point and length heteroplasmy in the mitochondrial DNA control region were compiled and analyzed from over 5,000 global human population samples. These data represent observations from a large and broad population sample, representing nearly 20 global populations. As expected, length heteroplasmy was frequently observed in the HVI, HVII and HVIII C-stretches. Length heteroplasmy was also observed in the AC dinucleotide repeat region, as well as other locations. Point heteroplasmy was detected in approximately 6% of all samples, and while the vast majority of heteroplasmic samples comprised two molecules differing at a single position, samples exhibiting two and three mixed positions were also observed in this data set. In general, the sites at which heteroplasmy was most commonly observed correlated with reported control region mutational hotspots. However, for some sites, observations of heteroplasmy did not mirror established mutation rate data, suggesting the action of other mechanisms, both selective and neutral. Interestingly, these data indicate that the frequency of heteroplasmy differs between particular populations, perhaps reflecting variable mutation rates among different mtDNA lineages and/or artifacts of particular population groups. The results presented here contribute to our general understanding of mitochondrial DNA control region heteroplasmy and provide additional empirical information on the mechanisms contributing to mtDNA control region mutation and evolution.

Evaluation of modified Yfiler amplification strategy for compromised samples.

AIM: To characterize the data produced using a modified amplification protocol for the AmpFlSTR Yfiler PCR Amplification Kit (Applied Biosystems) and explore the potential of Y-chromosomal short tandem repeat (Y-STR) recovery from severely degraded skeletal remains encountered at the Armed Forces DNA Identification Laboratory. METHODS: Experiments were performed using two sets of Yfiler amplification parameters. One set of parameters reflected the manufacturer's recommendations. The second set of parameters included twice the recommended Taq concentration and 6 additional cycles. Recovery of authentic alleles and the incidence of drop-in alleles were assessed for 3 data sets: 8 different quantities of pristine DNA, 8 artificially-degraded samples, and 31 non-probative case samples. RESULTS: Samples tested with both protocols from all 3 data sets yielded twice as many authentic alleles under the modified parameters than under the standard parameters (62% vs 31%), with only a nominal associated increase in the occurrence of non-authentic alleles (1.36% of all alleles detected). When applied to a range of representative casework samples, the modified protocol leveraged 9 or more reproducible alleles from over half of the specimens tested. CONCLUSION: Reproducible and informative Y-STR profiles can be recovered from a broad range of degraded and inhibited skeletal remains extracts when a commercially available kit is employed under modified amplification parameters.

Use of the LUS in sequence allele designations to facilitate probabilistic genotyping of NGS-based STR typing results.

Some of the expected advantages of next generation sequencing (NGS) for short tandem repeat (STR) typing include enhanced mixture detection and genotype resolution via sequence variation among non-homologous alleles of the same length. However, at the same time that NGS methods for forensic DNA typing have advanced in recent years, many caseworking laboratories have implemented or are transitioning to probabilistic genotyping to assist the interpretation of complex autosomal STR typing results. Current probabilistic software programs are designed for length-based data, and were not intended to accommodate sequence strings as the product input. Yet to leverage the benefits of NGS for enhanced genotyping and mixture deconvolution, the sequence variation among same-length products must be utilized in some form. Here, we propose use of the longest uninterrupted stretch (LUS) in allele designations as a simple method to represent sequence variation within the STR repeat regions and facilitate - in the nearterm - probabilistic interpretation of NGS-based typing results. An examination of published population data indicated that a reference LUS region is straightforward to define for most autosomal STR loci, and that using repeat unit plus LUS length as the allele designator can represent greater than 80% of the alleles detected by sequencing. A proof of concept study performed using a freely available probabilistic software demonstrated that the LUS length can be used in allele designations when a program does not require alleles to be integers, and that utilizing sequence information improves interpretation of both single-source and mixed contributor STR typing results as compared to using repeat unit information alone. The LUS concept for allele designation maintains the repeat-based allele nomenclature that will permit backward compatibility to extant STR databases, and the LUS lengths themselves will be concordant regardless of the NGS assay or analysis tools employed. Further, these biologically based, easy-to-derive designations uphold clear relationships between parent alleles and their stutter products, enabling analysis in fully continuous probabilistic programs that model stutter while avoiding the algorithmic complexities that come with string based searches. Though using repeat unit plus LUS length as the allele designator does not capture variation that occurs outside of the core repeat regions, this straightforward approach would permit the large majority of known STR sequence variation to be used for mixture deconvolution and, in turn, result in more informative mixture statistics in the near term. Ultimately, the method could bridge the gap from current length-based probabilistic systems to facilitate broader adoption of NGS by forensic DNA testing laboratories.

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.

Biological Sexing of a 4000-Year-Old Egyptian Mummy Head to Assess the Potential of Nuclear DNA Recovery from the Most Damaged and Limited Forensic Specimens.

High throughput sequencing (HTS) has been used for a number of years in the field of paleogenomics to facilitate the recovery of small DNA fragments from ancient specimens. Recently, these techniques have also been applied in forensics, where they have been used for the recovery of mitochondrial DNA sequences from samples where traditional PCR-based assays fail because of the very short length of endogenous DNA molecules. Here, we describe the biological sexing of a ~4000-year-old Egyptian mummy using shotgun sequencing and two established methods of biological sex determination (RX and RY), by way of mitochondrial genome analysis as a means of sequence data authentication. This particular case of historical interest increases the potential utility of HTS techniques for forensic purposes by demonstrating that data from the more discriminatory nuclear genome can be recovered from the most damaged specimens, even in cases where mitochondrial DNA cannot be recovered with current PCR-based forensic technologies. Although additional work remains to be done before nuclear DNA recovered via these methods can be used routinely in operational casework for individual identification purposes, these results indicate substantial promise for the retrieval of probative individually identifying DNA data from the most limited and degraded forensic specimens.

DNA identification of "Earthquake McGoon" 50 years postmortem.

This report describes the genetic identification of James "Earthquake McGoon" McGovern, a WWII fighter ace who perished in Laos while providing supplies to French troops during the French Indochina war. Because reference samples were unavailable for all of the potential casualties, testing of the entire mitochondrial genome, autosomal STRs and Y-chromosomal STRs was performed to increase the genetic information available for analysis. Kinship analyses performed on the evidentiary data and numerous indirect family references for McGovern excluded other possible casualties and definitively established McGovern's identity. This particular case demonstrates the practical utility of novel research technologies and aggressive genetic typing protocols in the identification of aged, degraded remains.

Application of low copy number STR typing to the identification of aged, degraded skeletal remains.

Low copy number (LCN) STR typing was successfully applied to four interesting cases during developmental validation of the approach for degraded skeletal remains. Specific questions were addressed in each case, with the acquisition of STR data largely serving as additional confirmatory or investigatory information in any specific situation, and not necessarily providing the definitive evidence to establish identity. The cases involve missing U.S. service members from World War I, World War II, and the Vietnam War. The variety of these cases, in terms of the questions addressed, the age of the remains, and the type of reference material available for comparison, demonstrates the broad utility of LCN STR typing in the identification of degraded skeletal remains from missing persons.

Performance and concordance of the ForenSeq™ system for autosomal and Y chromosome short tandem repeat sequencing of reference-type specimens.

Though the utility of next-generation sequencing (NGS) technologies for forensic short tandem repeat (STR) typing has been evident for several years, commercially available assays and software solutions developed specifically to meet forensic needs have only recently become available. One of these, the ForenSeq™ DNA Signature Prep Kit (Illumina, Inc.) sequences 27 autosomal STR (aSTR) and 24 Y chromosome STR (Y-STR) loci (concurrent with additional nuclear markers) per multiplexed sample, with automated secondary and tertiary analyses of the data accomplished via the associated ForenSeq™ Universal Analysis Software (UAS). In this study we investigated the performance of the ForenSeq system for aSTR and Y-STR typing by examination of 151 sample libraries developed from high quality DNAs amplified at the target 1ng template. Utilizing PCR Primer Mix B, greater than 99.5% of aSTR loci and 97.0% of Y-STR loci were recovered when 42 or fewer sample libraries were pooled for sequencing. A direct comparison of UAS developed fragment length results to capillary electrophoresis (CE) based data identified only two allele call discrepancies when no UAS quality flag was triggered. Review of the ForenSeq data indicated that most samples with total sequence read counts exceeding 40,000 could be interpreted to develop nearly complete aSTR genotypes or Y-STR haplotypes. However, markers D22S1045 and DYS392 produced poor or inconsistent results even when sample read counts were greater than 85,000. Excluding these two loci, analyst-interpreted aSTR and Y-STR ForenSeq profiles were 99.96% and 100% concordant, respectively, with CE data. In addition to demonstrating concordance on par with other CE kit to kit comparisons, the results from this study will assist laboratories seeking to develop workflows for high volume processing and analysis of aSTRs and Y-STRs from reference-type specimens using the ForenSeq system.

Genetic variation for 15 autosomal STR loci (PowerPlex 16) in a population sample from northern Greece.

Fifteen autosomal STR loci were typed in a population sample of 318 unrelated individuals from northern Greece. Allele frequencies were determined, and Hardy-Weinberg equilibrium tested, for the following loci: D5S818, D7S820, TH01, TPOX, vWA, CSF1PO, D16S539, D13S317, D8S1179, D3S1358, FGA, PENTA D, PENTA E, D21S11, D18S51.

Concordance and reproducibility of a next generation mtGenome sequencing method for high-quality samples using the Illumina MiSeq.

Sanger-type sequencing (STS) of mitochondrial DNA (mtDNA), specifically the control region (CR), is routinely employed in forensics in human identification and missing persons scenarios. Yet next-generation sequencing (NGS) has the potential to overcome some of the major limitations of STS processing, permitting reasonable paths forward for full mitochondrial genome (mtGenome) sequencing, while also offering higher-throughput and higher sensitivity capabilities. To establish the accuracy and reproducibility of NGS for the development of mtDNA data, 90 DNA extracts that were previously used to generate forensic quality full mtGenomes using STS were sequenced using Nextera XT library preparation and the Illumina MiSeq. Using the same amplicon product, replicate library sets were generated and sequenced at different laboratories, and analysis was performed in replicate using the CLC Genomics Workbench. Both sequencing sets resulted in 99.998% of positions with greater than 10X coverage when 96 samples (including controls) were multiplexed. Overall, 99.9996% concordance was observed between the NGS data and the STS data for the full mtGenome. The only "discordant" calls involved low level point heteroplasmies, with the differences resulting from stochastic variation and/or the increased sensitivity of NGS. Higher sensitivity also allowed for the detection of a mixed sample previously not detected with STS. Additionally, variant calls were reproducible between sequencing sets and between software analysis versions with the variant frequency only differing by 0.23% and 0.01%, respectively. Further validation studies and specialized software functionality tailored to forensic practice should facilitate the incorporation of NGS processing into standard casework applications. The data herein comprise the largest, and likely most thoroughly examined, complete mtGenome STS-NGS concordance dataset available.

Massively parallel sequencing of forensic STRs: Considerations of the DNA commission of the International Society for Forensic Genetics (ISFG) on minimal nomenclature requirements.

The DNA Commission of the International Society for Forensic Genetics (ISFG) is reviewing factors that need to be considered ahead of the adoption by the forensic community of short tandem repeat (STR) genotyping by massively parallel sequencing (MPS) technologies. MPS produces sequence data that provide a precise description of the repeat allele structure of a STR marker and variants that may reside in the flanking areas of the repeat region. When a STR contains a complex arrangement of repeat motifs, the level of genetic polymorphism revealed by the sequence data can increase substantially. As repeat structures can be complex and include substitutions, insertions, deletions, variable tandem repeat arrangements of multiple nucleotide motifs, and flanking region SNPs, established capillary electrophoresis (CE) allele descriptions must be supplemented by a new system of STR allele nomenclature, which retains backward compatibility with the CE data that currently populate national DNA databases and that will continue to be produced for the coming years. Thus, there is a pressing need to produce a standardized framework for describing complex sequences that enable comparison with currently used repeat allele nomenclature derived from conventional CE systems. It is important to discern three levels of information in hierarchical order (i) the sequence, (ii) the alignment, and (iii) the nomenclature of STR sequence data. We propose a sequence (text) string format the minimal requirement of data storage that laboratories should follow when adopting MPS of STRs. We further discuss the variant annotation and sequence comparison framework necessary to maintain compatibility among established and future data. This system must be easy to use and interpret by the DNA specialist, based on a universally accessible genome assembly, and in place before the uptake of MPS by the general forensic community starts to generate sequence data on a large scale. While the established nomenclature for CE-based STR analysis will remain unchanged in the future, the nomenclature of sequence-based STR genotypes will need to follow updated rules and be generated by expert systems that translate MPS sequences to match CE conventions in order to guarantee compatibility between the different generations of STR data.