Targeted enrichment of whole-genome SNPs from highly burned skeletal remains.

Genetic assessment of highly incinerated and/or degraded human skeletal material is a persistent challenge in forensic DNA analysis, including identifying victims of mass disasters. Few studies have investigated the impact of thermal degradation on whole-genome single-nucleotide polymorphism (SNP) quality and quantity using next-generation sequencing (NGS). We present whole-genome SNP data obtained from the bones and teeth of 27 fire victims using two DNA extraction techniques. Extracts were converted to double-stranded DNA libraries then enriched for whole-genome SNPs using unpublished biotinylated RNA baits and sequenced on an Illumina NextSeq 550 platform. Raw reads were processed using the EAGER (Efficient Ancient Genome Reconstruction) pipeline, and the SNPs filtered and called using FreeBayes and GATK (v. 3.8). Mixed-effects modeling of the data suggest that SNP variability and preservation is predominantly determined by skeletal element and burn category, and not by extraction type. Whole-genome SNP data suggest that selecting long bones, hand and foot bones, and teeth subjected to temperatures <350°C are the most likely sources for higher genomic DNA yields. Furthermore, we observed an inverse correlation between the number of captured SNPs and the extent to which samples were burned, as well as a significant decrease in the total number of SNPs measured for samples subjected to temperatures >350°C. Our data complement previous analyses of burned human remains that compare extraction methods for downstream forensic applications and support the idea of adopting a modified Dabney extraction technique when traditional forensic methods fail to produce DNA yields sufficient for genetic identification.

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.

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.

Development of a multiplex single base extension assay for mitochondrial DNA haplogroup typing.

AIM: To provide a screening tool to reduce time and sample consumption when attempting mitochondrial DNA (mtDNA) haplogroup typing. METHODS: A single base primer extension assay was developed to enable typing, in a single reaction, of twelve mtDNA haplogroup specific polymorphisms. For validation purposes a total of 147 samples were tested including 73 samples successfully haplogroup typed using mtDNA control region (CR) sequence data, 20 samples inconclusively haplogroup typed by CR sequence data, 21 samples previously haplogroup typed using RFLP analysis, and 31 samples of known ancestral origin without previous haplogroup typing. Additionally, two highly degraded human bones embalmed and buried in the early 1950s were analyzed using the SNP multiplex. RESULTS: When the SNP multiplex was used to type the 96 previously CR sequenced specimens, an increase in haplogroup or macrohaplogroup assignment relative to conventional CR sequence analysis was observed. The single base extension assay was also successfully used to assign a haplogroup to decades-old, embalmed skeletal remains dating to World War II. CONCLUSION: The SNP multiplex was successfully used to obtain haplogroup status of highly degraded human bones, and demonstrated the ability to eliminate possible contributors. The SNP multiplex provides a low-cost, high throughput method for typing of mtDNA haplogroups A, B, C, D, E, F, G, H, L1/L2, L3, M, and N that could be useful for screening purposes for human identification efforts and anthropological studies.

Complete Mitochondrial Genome Sequencing of a Burial from a Romano-Christian Cemetery in the Dakhleh Oasis, Egypt: Preliminary Indications.

The curse of ancient Egyptian DNA was lifted by a recent study which sequenced the mitochondrial genomes (mtGenome) of 90 ancient Egyptians from the archaeological site of Abusir el-Meleq. Surprisingly, these ancient inhabitants were more closely related to those from the Near East than to contemporary Egyptians. It has been accepted that the timeless highway of the Nile River seeded Egypt with African genetic influence, well before pre-Dynastic times. Here we report on the successful recovery and analysis of the complete mtGenome from a burial recovered from a remote Romano-Christian cemetery, Kellis 2 (K2). K2 serviced the ancient municipality of Kellis, a village located in the Dakhleh Oasis in the southwest desert in Egypt. The data were obtained by high throughput sequencing (HTS) performed independently at two ancient DNA facilities (Armed Forces DNA Identification Laboratory, Dover, DE, USA and Carl R. Woese Institute for Genomic Biology, University of Illinois Urbana-Champaign, Urbana, IL, USA). These efforts produced concordant haplotypes representing a U1a1a haplogroup lineage. This result indicates that Near Eastern maternal influence previously identified at Abusir el-Meleq was also present further south, in ancient Kellis during the Romano-Christian period.

Short tandem repeat typing on the 454 platform: strategies and considerations for targeted sequencing of common forensic markers.

To investigate the feasibility of next generation sequencing technology (NGS) for the multiplex detection and sequence production of short tandem repeats (STRs) from degraded and low DNA quantity samples, standard polymerase chain reaction amplification methods were used to enrich for commonly employed STR markers. Samples were amplified with two multiplexing strategies: a multiplex containing thirteen miniSTR markers and a series of multiplexes containing four miniSTR markers each. Each sample multiplex was barcoded with a sample-specific multiplex identifier for subsequent parallel tagged sequencing on the GS Junior System (454 Life Sciences, a Roche company, Branford, CT). Sequencing results from over fifty DNA extracts representing both pristine samples and low-quality evidentiary specimens reflected known genotypes and were consistent across multiple extracts and/or amplifications of the same sample. Furthermore, the NGS data revealed sequence information not available with standard capillary electrophoresis-based detection alone. For the population samples tested, a total of 152 single nucleotide polymorphisms or insertions/deletions were identified in over 935 recovered alleles, averaging one polymorphism for every six recovered alleles. For three of the loci, the sequence information doubled the number of alleles detected via traditional STR typing by fragment analysis. In addition, twenty-eight of these variants were only seen once within our dataset, highlighting the potential for discrimination among individuals. These additional data are likely to be particularly valuable in missing persons and disaster victim identification cases for which only partial profiles may be recovered and/or only distant kin are available as references. And, considering the opportunity to target only small amplicons with NGS, this type of STR typing will allow for greater information recovery from challenging casework samples. While our results highlight the potential of new technologies for recovering discriminatory genetic information from evidentiary specimens, our data also reveal the complexities of NGS-based STR typing, both in terms of the laboratory assays themselves as well as the downstream data processing and analysis.

Capillary electrophoresis of human mtDNA control region sequences from highly degraded samples using short mtDNA amplicons.

The forensic applications of mtDNA sequencing center primarily on samples that are either highly degraded or contain little or no nuclear DNA, since the testing of these sample types is often unsuccessful with more widely used nuclear STR profiling assays. In these cases, sequence data from the noncoding mtDNA control region are targeted due to its high variability. However, the ease of authentic DNA recovery and the strategy used for recovery depend strictly on the quality of the sample. In this chapter, we will cover mitochondrial DNA sequencing procedures for short mtDNA amplicons which range in size from 100 to 350 bp. Generally speaking, amplicons of this size are required only for the most degraded specimens, and the protocols described here have been specifically developed for recalcitrant human skeletal remains encountered during the course of a large-scale missing persons' identification effort. DNA templates from these types of specimens tend to exhibit various forms of intrastrand damage that, in turn, manifest as artifacts in the sequence data. Because these artifacts are not generally observed among sequence data from pristine templates, we address the particular data idiosyncrasies that warrant additional scrutiny. Although this chapter will primarily highlight this particular application, the basic experimental parameters and data considerations should easily extend to other applications and/or sample types. The protocols described here have been deliberately designed to produce raw sequence electropherograms and final mtDNA profiles that adhere to the strictest forensic guidelines in terms of overall data quality.

Characterization of a modified amplification approach for improved STR recovery from severely degraded skeletal elements.

Degraded skeletal remains generally contain limited quantities of genetic material and thus DNA-based identification efforts often target the mitochondrial DNA (mtDNA) control region due to the relative abundance of intact mtDNA as compared to nuclear DNA. In many missing person cases, however, the discriminatory power of mtDNA is inadequate to permit identification when associated anthropological, odontological, or contextual evidence is also limited, and/or the event involves a large number of individuals. In situations such as these, more aggressive amplification protocols which can permit recovery of STR data are badly needed as they may represent the last hope for conclusive identification. We have previously demonstrated the potential of a modified Promega PowerPlex 16 amplification strategy for the recovery of autosomal STR data from severely degraded skeletal elements. Here, we further characterize the results obtained under these modified parameters on a variety of sample types including pristine control DNA and representative case work specimens. Not only is the amplification approach evaluated here sensitive to extremely low authentic DNA input quantities (6 pg), but when the method was applied to thirty-one challenging casework specimens, nine or more alleles were reproducibly recovered from 69% of the samples tested. Moreover, when we independently considered bone samples extracted with a protocol that includes complete demineralization of the bone matrix, the percentage of samples yielding nine or more reproducible alleles increased to 95% with the modified amplification parameters. Overall, direct comparisons between the modified amplification protocol and the standard amplification protocol demonstrated that allele recovery was significantly greater using the aggressive parameters, with only a minimal associated increase in artifactual data.

Automatable full demineralization DNA extraction procedure from degraded skeletal remains.

During the 7 year period from 2002 to 2009 a high volume, silica-binding DNA extraction protocol for bone, based on modified QIAGEN's Blood Maxi Kit protocol was highly successful permitting the DNA matching of >14,500 missing persons from former Yugoslavia. This method, however, requires large amount of bone material and large volumes of reagents. The logical evolution was to develop a more efficient extraction protocol for bone samples that uses significantly less starting material while increasing the success in obtaining DNA results from smaller, more challenging samples. In this study we compared the performance of ICMP's original protocol against an automatable full demineralization approach. In order to provide reliable results and to simulate a wide variety of cases, we analyzed 40 bone samples in a comparative study based on DNA concentrations and quality of resulting STR profiles. The new protocol results in the dissolution of the entire bone powder sample, thus eliminating the possibility that DNA is left behind, locked in remaining solid bone matrix. For the majority of samples tested, the DNA concentrations obtained from half a gram of fully digested bone material were equivalent to or greater than the ones obtained from 2g of partially demineralized bone powder. Furthermore, the full demineralization process significantly increases the proportion of full profiles reflecting the correlation with better DNA quality. This method has been adapted for the QIAcube robotic platform. The performance of this automated full demineralization protocol is similar to the manual version and increases overall lab throughput. It also simplifies the process by eliminating quality control procedures that are advisable in manual procedures, and overall reduces the chance of human error. Finally we described a simple and efficient post-extraction clean-up method that can be applied to DNA extracts obtained from different protocols. This protocol has also been adjusted for the QIAcube platform.

Assessing a novel room temperature DNA storage medium for forensic biological samples.

The ability to properly collect, analyze and preserve biological stains is important to preserving the integrity of forensic evidence. Stabilization of intact biological evidence in cells and the DNA extracts from them is particularly important since testing is generally not performed immediately following collection. Furthermore, retesting of stored DNA samples may be needed in casework for replicate testing, confirmation of results, and to accommodate future testing with new technologies. A novel room temperature DNA storage medium, SampleMatrix™ (SM; Biomatrica, Inc., San Diego, CA), was evaluated for stabilizing and protecting samples. Human genomic DNA samples at varying amounts (0.0625-200 ng) were stored dry in SM for 1 day to 1 year under varying conditions that included a typical ambient laboratory environment and also through successive freeze-thaw cycles (3 cycles). In addition, spiking of 1-4 × SM into samples prior to analysis was performed to determine any inhibitory effects of SM. Quantification of recovered DNA following storage was determined by quantitative PCR or by agarose gel electrophoresis, and evaluation of quantitative peak height results from multiplex short tandem repeat (STR) analyses were performed to assess the efficacy of SM for preserving DNA. Results indicate no substantial differences between the quality of samples stored frozen in liquid and those samples maintained dry at ambient temperatures protected in SM. For long-term storage and the storage of low concentration samples, SM provided a significant advantage over freezer storage through higher DNA recovery. No detectable inhibition of amplification was observed at the recommended SM concentration and complete profiles were obtained from genomic DNA samples even in the presence of higher than recommended concentrations of the SM storage medium. The ability to stabilize and protect DNA from degradation at ambient temperatures for extended time periods could have tremendous impact in simplifying and improving sample storage conditions and requirements. The current work focuses on forensics analysis; however this technology is applicable to all endeavors requiring storage of DNA.

Titanic's unknown child: the critical role of the mitochondrial DNA coding region in a re-identification effort.

This report describes a re-examination of the remains of a young male child recovered in the Northwest Atlantic following the loss of the Royal Mail Ship Titanic in 1912 and buried as an unknown in Halifax, Nova Scotia shortly thereafter. Following exhumation of the grave in 2001, mitochondrial DNA (mtDNA) hypervariable region 1 sequencing and odontological examination of the extremely limited skeletal remains resulted in the identification of the child as Eino Viljami Panula, a 13-month-old Finnish boy. This paper details recent and more extensive mitochondrial genome analyses that indicate the remains are instead most likely those of an English child, Sidney Leslie Goodwin. The case demonstrates the benefit of targeted mtDNA coding region typing in difficult forensic cases, and highlights the need for entire mtDNA sequence databases appropriate for forensic use.

Integrated DNA and fingerprint analyses in the identification of 60-year-old mummified human remains discovered in an Alaskan glacier.

This report describes the identification of a merchant mariner who perished in 1948 when Northwest Airlines Flight 4422, a DC-4 carrying 24 seamen and six crew members crashed into Mount Sanford, Alaska. Fifty-one years later, a human forearm and hand were found close by the wreckage of the plane, prompting identification efforts using DNA and fingerprints. There were significant challenges to both the fingerprint and DNA analyses. The hand was badly desiccated, making fingerprint friction-ridge detail almost invisible and the remains had been embalmed upon discovery, making DNA amplification difficult. We present the results of an interdisciplinary approach that successfully addressed these challenges and ultimately led to the identification of the remains. These efforts relied on efficient fingerprint rejuvenation and imaging techniques that improved print resolution, as well as new DNA extraction techniques optimized for aggressively embalmed remains.

Mystery solved: the identification of the two missing Romanov children using DNA analysis.

One of the greatest mysteries for most of the twentieth century was the fate of the Romanov family, the last Russian monarchy. Following the abdication of Tsar Nicholas II, he and his wife, Alexandra, and their five children were eventually exiled to the city of Yekaterinburg. The family, along with four loyal members of their staff, was held captive by members of the Ural Soviet. According to historical reports, in the early morning hours of July 17, 1918 the entire family along with four loyal members of their staff was executed by a firing squad. After a failed attempt to dispose of the remains in an abandoned mine shaft, the bodies were transported to an open field only a few kilometers from the mine shaft. Nine members of the group were buried in one mass grave while two of the children were buried in a separate grave. With the official discovery of the larger mass grave in 1991, and subsequent DNA testing to confirm the identities of the Tsar, the Tsarina, and three of their daughters--doubt persisted that these remains were in fact those of the Romanov family. In the summer of 2007, a group of amateur archeologists discovered a collection of remains from the second grave approximately 70 meters from the larger grave. We report forensic DNA testing on the remains discovered in 2007 using mitochondrial DNA (mtDNA), autosomal STR, and Y-STR testing. Combined with additional DNA testing of material from the 1991 grave, we have virtually irrefutable evidence that the two individuals recovered from the 2007 grave are the two missing children of the Romanov family: the Tsarevich Alexei and one of his sisters.

High efficiency DNA extraction from bone by total demineralization.

In historical cases, missing persons' identification, mass disasters, and ancient DNA investigations, bone and teeth samples are often the only, and almost always the best, biological material available for DNA typing. This is because of the physical and chemical barrier that the protein:mineral matrix of bone poses to environmental deterioration and biological attack. Most bone extraction protocols utilized in the forensic community involve an incubation period of bone powder in extraction buffer for proteinase digestion, followed by the collection of the supernatant, and the disposal of large quantities of undissolved bone powder. Here we present an extremely efficient protocol for recovery of DNA by complete demineralization, resulting in full physical dissolution of the bone sample. This is performed in a manner that retains and concentrates all the reagent volume, for complete DNA recovery. For our study, we selected 14 challenging bone samples. The bones were extracted side-by-side with our new demineralization protocol and the standard extraction protocol in use at AFDIL. A real-time quantification assay based on the amplification of a 143 bp mtDNA fragment showed that this new demineralization protocol significantly enhances the quantity of DNA that can be extracted and amplified from degraded skeletal remains. We have used this technique to successfully recover authentic DNA sequences from extremely challenging samples that failed repeatedly using the standard protocol.

Evolution of ascariasis in humans and pigs: a multi-disciplinary approach.

The nematode parasite Ascaris lumbricoides infects the digestive tracts of over 1.4 billion people worldwide, and its sister species, Ascaris suum, has infected a countless number of domesticated and feral pigs. It is generally thought that the putative ancestor to these worms infected either humans or pigs, but with the advent of domestication, they had ample opportunity to jump to a new host and subsequently specialize and evolve into a new species. While nuclear DNA makers decisively separate the two populations, mitochondrial sequences reveal that three major haplotypes are found in A. suum and in A. lumbricoides, indicating either occasional hybridization, causing introgression of gene trees, or retention of polymorphism dating back to the original ancestral species. This article provides an illustration of the combined contribution of parasitology, archaeoparasitology, genetics and paleogenetics to the history of ascariasis. We specifically investigate the molecular history of ascariasis in humans by sequencing DNA from the eggs of Ascaris found among ancient archeological remains. The findings of this paleogenetic survey will explain whether the three mitochondrial haplotypes result from recent hybridization and introgression, due to intensive human-pig interaction, or whether their co-occurrence predates pig husbandry, perhaps dating back to the common ancestor. We hope to show how human-pig interaction has shaped the recent evolutionary history of this disease, perhaps revealing the identify of the ancestral host.

Molecular phylogeny of the extinct cave lion Panthera leo spelaea.

To reconstruct the phylogenetic position of the extinct cave lion (Panthera leo spelaea), we sequenced 1 kb of the mitochondrial cytochrome b gene from two Pleistocene cave lion DNA samples (47 and 32 ky B.P.). Phylogenetic analysis shows that the ancient sequences form a clade that is most closely related to the extant lions from Africa and Asia; at the same time, cave lions appear to be highly distinct from their living relatives. Our data show that these cave lion sequences represent lineages that were isolated from lions in Africa and Asia since their dispersal over Europe about 600 ky B.P., as they are not found among our sample of extant populations. The cave lion lineages presented here went extinct without mitochondrial descendants on other continents. The high sequence divergence in the cytochrome b gene between cave and modern lions is notable.

Evolution of ascariasis in humans and pigs: a multi-disciplinary approach

The nematode parasite Ascaris lumbricoides infects the digestive tracts of over 1.4 billion people worldwide, and its sister species, Ascaris suum, has infected a countless number of domesticated and feral pigs. It is generally thought that the putative ancestor to these worms infected either humans or pigs, but with the advent of domestication, they had ample opportunity to jump to a new host and subsequently specialize and evolve into a new species. While nuclear DNA markers decisively separate the two populations, mitochondrial sequences reveal that three major haplotypes are found in A. suum and in A. lumbricoides, indicating either occasional hybridization, causing introgression of gene trees, or retention of polymorphism dating back to the original ancestral species. This article provides an illustration of the combined contribution of parasitology, archaeoparasitology, genetics and paleogenetics to the history of ascariasis. We specifically investigate the molecular history of ascariasis in humans by sequencing DNA from the eggs of Ascaris found among ancient archeological remains. The findings of this paleogenetic survey will explain whether the three mitochondrial haplotypes result from recent hybridization and introgression, due to intensive human-pig interaction, or whether their co-occurrence predates pig husbandry, perhaps dating back to the common ancestor. We hope to show how human-pig interaction has shaped the recent evolutionary history of this disease, perhaps revealing the identity of the ancestral host