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.

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.

Unearthing who and Y at Harewood Cemetery and inference of George Washington's Y-chromosomal haplotype.

An excavation conducted at Harewood Cemetery to identify the unmarked grave of Samuel Washington resulted in the discovery of burials presumably belonging to George Washington's paternal grandnephews and their mother, Lucy Payne. To confirm their identities this study examined Y-chromosomal, mitochondrial, and autosomal DNA from the burials and a living Washington descendant. The burial's Y-STR profile was compared to FamilyTreeDNA's database, which resulted in a one-step difference from the living descendant and an exact match to another Washington. A more complete Y-STR and Y-SNP profile from the descendant was inferred to be the Washington Y profile. Kinship comparisons performed in relation to the descendant, who is a 4th and 5th degree relative of the putative individuals, resulted in >37,000 overlapping autosomal SNPs and strong statistical support with likelihood ratios exceeding one billion. This study highlights the benefits of a multi-marker approach for kinship prediction and DNA-assisted identification of historical remains.

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.

Impact of SNP microarray analysis of compromised DNA on kinship classification success in the context of investigative genetic genealogy.

Single nucleotide polymorphism (SNP) data generated with microarray technologies have been used to solve murder cases via investigative leads obtained from identifying relatives of the unknown perpetrator included in accessible genomic databases, an approach referred to as investigative genetic genealogy (IGG). However, SNP microarrays were developed for relatively high input DNA quantity and quality, while DNA typically obtainable from crime scene stains is of low DNA quantity and quality, and SNP microarray data obtained from compromised DNA are largely missing. By applying the Illumina Global Screening Array (GSA) to 264 DNA samples with systematically altered quantity and quality, we empirically tested the impact of SNP microarray analysis of compromised DNA on kinship classification success, as relevant in IGG. Reference data from manufacturer-recommended input DNA quality and quantity were used to estimate genotype accuracy in the compromised DNA samples and for simulating data of different degree relatives. Although stepwise decrease of input DNA amount from 200 ng to 6.25 pg led to decreased SNP call rates and increased genotyping errors, kinship classification success did not decrease down to 250 pg for siblings and 1st cousins, 1 ng for 2nd cousins, while at 25 pg and below kinship classification success was zero. Stepwise decrease of input DNA quality via increased DNA fragmentation resulted in the decrease of genotyping accuracy as well as kinship classification success, which went down to zero at the average DNA fragment size of 150 base pairs. Combining decreased DNA quantity and quality in mock casework and skeletal samples further highlighted possibilities and limitations. Overall, GSA analysis achieved maximal kinship classification success from 800 to 200 times lower input DNA quantities than manufacturer-recommended, although DNA quality plays a key role too, while compromised DNA produced false negative kinship classifications rather than false positive ones.

Challenges with co-amplification of microbial DNA in interpretation of STR profiles obtained from human skeletal remains.

STR-based DNA analysis is still the main tool for human DNA identification in most forensic DNA laboratories. DNA typing of aged human skeletal elements faces well-known interpretation challenges characteristic of degraded and low copy number DNA samples. Analyzing tens of thousands of human bone and teeth samples, we found that the occasional presence of artefactual peaks of presumed microbial origin adds another layer of complexity to the interpretation of STR profiles. In this paper, we present our approach and suggest guidelines for identifying and distinguishing non-human peaks, developed over the last 18 years. Additionally, we report a compendium of artefact peaks of presumed microbial origin recorded in human STR profiles obtained from bone and teeth samples, originating from Iraq, Chile, Maldives, Brazil and Western Balkans. Our experience has shown that these artefacts are not uncommon in bone STR testing, suggesting the possibility of occurrence in other forensic contexts, particularly trace DNA samples. Raising awareness among the forensic DNA community and accounting for this phenomenon is important for accurate STR interpretation.

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.

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.

Large scale DNA identification: The ICMP experience.

The International Commission on Missing Persons (ICMP) is a treaty-based international organization with a global mandate to address the issue of missing persons. It works with governments, civil society organizations, and others, and utilizes data systems and technical assistance in forensic science. ICMP's initial work focused on the ∼40,000 people missing in the Western Balkans from the conflicts of the 1990s. A "DNA-led" approach to large-scale DNA identification of the missing was developed, based on high-throughput autosomal STR testing of skeletal remains from mass graves and other sites, and the establishment of a regional database of DNA profiles from family members of the missing. Database pairwise and pedigree kinship searching is conducted using in-house DNA matching software, the Identification Data Management System (iDMS), providing high-certainty DNA matches that are integrated in a multi-disciplinary identification process. Anthropological guidelines for sampling skeletal remains for DNA testing are based on tens of thousands of tests from a wide range of skeletal elements, allowing for prioritization based on DNA preservation. Large-scale collection of family reference samples has been conducted, resulting in a database of more than 100,000 family reference DNA profiles across all projects and delivering family DNA match reports for more than 20,000 individuals. From the 1995 Srebrenica event, ICMP provided DNA matches for 6887 of the ∼8000 missing from that event. In assistance to justice, ICMP has provided extensive evidence and expert testimony in multiple war crimes trials, including those conducted at the ICTY. This article provides an overview of ICMP's technical involvement over the last 17 years in areas of DNA testing and database matching, and training and capacity building projects with partners. It also touches on the development of massively parallel sequencing (MPS) strategies specifically tailored to missing persons applications.

An economical and efficient method for postmortem DNA sampling in mass fatalities.

In mass fatality events, the need to identify large numbers of deceased persons using DNA can be a significant drain on already overburdened forensic practitioners, both in the field setting and the laboratory. The laboratory may be required to extract DNA from a variety of postmortem sample types, family or direct reference samples related to the missing, and perform matching of these results in a short period of time. While most forensic institutions are well equipped to handle both family and direct reference samples, postmortem samples such as bone or heterogeneous tissue samples can be difficult for labs to analyze. We have devised an easily deployable, efficient, and inexpensive method for collecting postmortem DNA samples on commercially available DNA preservation cards ("FTA®" cards). FTA® cards are already widely used in forensic labs and are convenient for shipping due to their small volume and stability at room temperature. We evaluated the suitability of a protocol involving swabbing of incisions made on cadavers and sample deposition onto FTA® cards over various postmortem intervals and under different environmental conditions. Each trial took place during a different point in the calendar year to evaluate the effects of seasonal weather patterns and temperature on decomposition, DNA yield, and rates of DNA degradation. To further account for the effects of seasonality (temperature and humidity), the progression of body decomposition was recorded following the Total Body Score (TBS) method [1]. DNA degradation was assessed either through STR amplification of 1.2 mm FTA punches or DNA extraction from 3.0 mm punches followed by real-time PCR quantification and STR amplification and genotyping. No consistent relationship was observed between postmortem interval and DNA degradation. Instead, the TBS score, which captures the stage of body decomposition, was shown to correlate well with DNA quantity. A TBS of 15 and below consistently yielded strong partial or full profiles (20 STR loci and Amelogenin using the PowerPlex 21 System) from all individuals from either 1.2 mm or 3.0 mm punches. Transfer of sample swabs to FTA cards is shown to be a simple and effective method for both field and laboratory operations over a range of conditions that can be evaluated by field forensic practitioners based on a body decomposition score. The approach could be beneficially integrated into mass fatality response plans.

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.

Highly effective DNA extraction method for nuclear short tandem repeat testing of skeletal remains from mass graves.

AIM: To quantitatively compare a silica extraction method with a commonly used phenol/chloroform extraction method for DNA analysis of specimens exhumed from mass graves. METHODS: DNA was extracted from twenty randomly chosen femur samples, using the International Commission on Missing Persons (ICMP) silica method, based on Qiagen Blood Maxi Kit, and compared with the DNA extracted by the standard phenol/chloroform-based method. The efficacy of extraction methods was compared by real time polymerase chain reaction (PCR) to measure DNA quantity and the presence of inhibitors and by amplification with the PowerPlex 16 (PP16) multiplex nuclear short tandem repeat (STR) kit. RESULTS: DNA quantification results showed that the silica-based method extracted on average 1.94 ng of DNA per gram of bone (range 0.25-9.58 ng/g), compared with only 0.68 ng/g by the organic method extracted (range 0.0016-4.4880 ng/g). Inhibition tests showed that there were on average significantly lower levels of PCR inhibitors in DNA isolated by the organic method. When amplified with PP16, all samples extracted by silica-based method produced 16 full loci profiles, while only 75% of the DNA extracts obtained by organic technique amplified 16 loci profiles. CONCLUSIONS: The silica-based extraction method showed better results in nuclear STR typing from degraded bone samples than a commonly used phenol/chloroform method.

Success rates of nuclear short tandem repeat typing from different skeletal elements.

AIM: To evaluate trends in DNA typing success rates of different skeletal elements from mass graves originating from conflicts that occurred in the former Yugoslavia (Bosnia and Herzegovina and Kosovo) during the 1990s, and to establish correlation between skeletal sample age and success of high throughput short tandem repeat (STR) typing in the large data set of the International Commission on Missing Persons. METHOD: DNA extraction and short tandem repeat (STR) typing have been attempted on over 25000 skeletal samples. The skeletal samples originated from different geographical locations where the conflicts occurred and from different time periods from 1992 to 1999. DNA preservation in these samples was highly variable, but was often significantly degraded and of limited quantity. For the purpose of this study, processed samples were categorized according to skeletal sample type, sample age since death, and success rates tabulated. RESULTS: Well-defined general trends in success rates of DNA analyses were observed with respect to the type of bone tested and sample age. The highest success rates were observed with samples from dense cortical bone of weight-bearing leg bones (femur 86.9%), whereas long bones of the arms showed significantly lower success (humerus 46.2%, radius 24.5%, ulna 22.8%). Intact teeth also exhibited high success rates (teeth 82.7%). DNA isolation from other skeletal elements differed considerably in success, making bone sample selection an important factor influencing success. CONCLUSION: The success of DNA typing is related to the type of skeletal sample. By carefully evaluating skeletal material available for forensic DNA testing with regard to sample age and type of skeletal element available, it is possible to increase the success and efficiency of forensic DNA testing.

Variant alleles, triallelic patterns, and point mutations observed in nuclear short tandem repeat typing of populations in Bosnia and Serbia.

AIM: To present a compendium of off-ladder alleles and other genotyping irregularities relating to rare/unexpected population genetic variation, observed in a large short tandem repeat (STR) database from Bosnia and Serbia. METHODS: DNA was extracted from blood stain cards relating to reference samples from a population of 32800 individuals from Bosnia and Serbia, and typed using Promega's PowerPlex16 STR kit. RESULTS: There were 31 distinct off-ladder alleles were observed in 10 of the 15 STR loci amplified from the PowerPlex16 STR kit. Of these 31, 3 have not been previously reported. Furthermore, 16 instances of triallelic patterns were observed in 9 of the 15 loci. Primer binding site mismatches that affected amplification were observed in two loci, D5S818 and D8S1179. CONCLUSION: Instances of deviations from manufacturer's allelic ladders should be expected and caution taken to properly designate the correct alleles in large DNA databases. Particular care should be taken in kinship matching or paternity cases as incorrect designation of any of these deviations from allelic ladders could lead to false exclusions.

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.

Identification of Missing Norwegian World War II Soldiers, in Karelia Russia.

This article presents the multidisciplinary effort in trying to identify the skeletal remains of 100 Norwegian soldiers serving in the German army, killed in Karelia Russia in 1944, from the recovery of the remains through the final identification using DNA. Of the 150 bone samples sent for DNA testing, 93 DNA profiles were obtained relating to 57 unique individuals. The relatives could not be directly contacted as the soldiers were considered as traitors to Norway; therefore, only 45 reference samples, relating to 42 cases of the missing, were donated. DNA matches for 14 soldiers and 12 additional body part re-associations for these individuals were found. Another 24 bone samples were re-associated with 16 individuals, but no familial match was found. More than six decades after the end of WWII, DNA analysis can significantly contribute to the identification of the remains.

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

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

Toward increased utility of mtDNA in forensic identifications.

The utility of mtDNA in forensic identifications is limited by its low power of discrimination and the absence of high quality mtDNA databases. Single nucleotide polymorphisms (SNPs) in the control region outside of hypervariable regions I and II (HVI/HVII), and in the coding region of the mtDNA genome, can provide additional discrimination in mtDNA testing. We have identified particularly useful SNP sites via high throughput sequencing of the entire mtDNA genome. We report here two cases in which an 11-plex SNP assay (panel "A") targeting the most common HVI/HVII type successfully resolved two cases in which identifications could not be made on the basis of HVI/HVII sequencing. Additionally, we established a database of 286 samples for SNP panel "A" generated with robotic protocols. We have addressed the need for high quality mtDNA control region (CR) databases by developing robotic protocols for lab processing, and a carefully devised electronic data review process. A large-scale databasing effort targeting several populations underrepresented in current mtDNA databases is underway.