Next Article in Journal
Translational Advances of Hydrofection by Hydrodynamic Injection
Next Article in Special Issue
Dating Pupae of the Blow Fly Calliphora vicina Robineau–Desvoidy 1830 (Diptera: Calliphoridae) for Post Mortem Interval—Estimation: Validation of Molecular Age Markers
Previous Article in Journal
The Unforeseen Non-Coding RNAs in Head and Neck Cancer
Previous Article in Special Issue
Deep-Coverage MPS Analysis of Heteroplasmic Variants within the mtGenome Allows for Frequent Differentiation of Maternal Relatives
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 9
Issue 3
10.3390/genes9030135
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Odile Loreille
1,*,
Shashikala Ratnayake
2,
Adam L. Bazinet
2,
Timothy B. Stockwell
2,
Daniel D. Sommer
2,
Nadin Rohland
3,
Swapan Mallick
3,
Philip L.F. Johnson
4,
Pontus Skoglund
5,
Anthony J. Onorato
1,
Nicholas H. Bergman
2,
David Reich
3,6 and
Jodi A. Irwin
1
1
DNA Support Unit, FBI Laboratory, 2501 Investigation Parkway, Quantico, VA 22135, USA
2
National Biodefense Analysis and Countermeasures Center, 8300 Research Plaza, Fort Detrick, MD 21702, USA
3
Department of Genetics Harvard Medical School, 77 Avenue Louis Pasteur, Boston, MA 02115, USA
4
Department of Biology, University of Maryland, 1210 Biology-Psychology Building, 4094 Campus Drive, College Park, MD 20742, USA
5
The Francis Crick Institute, 1 Midland Rd, London NW1 1AT, UK
6
Broad Institute of MIT and Harvard, 415 Main Street, Cambridge, MA 02142, USA
*
Author to whom correspondence should be addressed.
Genes 2018, 9(3), 135; https://doi.org/10.3390/genes9030135
Submission received: 12 January 2018 / Revised: 6 February 2018 / Accepted: 6 February 2018 / Published: 1 March 2018
(This article belongs to the Special Issue Forensic Genomics)
Browse Figures
Versions Notes

Abstract

:
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.

1. Introduction

On 18 October 2009, the Museum of Fine Arts, Boston (MFA) opened an exhibition called The Secrets of Tomb 10A: Egypt 2000 BC. One of the most intriguing items in the collection was a mummified human head discovered over a century ago in the necropolis of Deir el-Bersha (also known as Dayr al-Barshā). The site is located on the east bank of the Nile River in close proximity to the city of Mallawi, approximately 250 km south of Cairo. Deir el-Bersha is known for tombs cut into cliffs of limestone that date back to the Old Kingdom (2686–2181 BC), the First Intermediate Period (about 2100–2040 BC) and the Middle Kingdom (MK; 2040–approx. 1640 BC). During the MK 11th and 12th Dynasties (2040–1783 BC), it served as the chief cemetery for the governors or regional lords (a.k.a. nomarchs) of the 15th Upper Egyptian Nome (a.k.a. the Hare Nome).
In 1915, Deir el-Bersha was excavated by a joint MFA-Harvard University team directed by George A. Reisner with the assistance of Hanford L. Story and Said Ahmed Said [1]. On 23 April, Reisner’s men began clearing the burial shaft from tomb number 10A. Six days later, at the bottom of a 30-foot pit, they discovered the burial chamber of an early MK governor named Djehutynakht. It is not yet clear whether he is Djehutynakht IV, son of Ahanakht I, or Djehutynakht V, son of Nehri I [1,2,3]. Both were nomarchs of the Hare Nome, and while they shared the same name (which means Thoth [the main local deity] is Strong), there is no evidence they were related. Tomb 10A contained a second occupant: the governor’s wife, who was also called Djehutynakht. Although the tomb had been plundered in antiquity and most of the valuable jewels stolen, many objects were left behind. In fact, the recovered contents of Tomb 10A are considered one of the largest burial assemblages of the MK ever discovered. The funerary equipment includes pottery, canopic jars, models showing men and women in different daily life activities, nearly 60 model boats, and a famous, exquisitely carved and painted processional group composed of a priest and four offering bearers, known as the “Bersha procession” [4]. Befitting their high status, the governor and his wife were both buried in finely decorated rectangular wooden coffins placed within larger coffins, all made of thick cedar of Lebanon boards. Upon discovery, the coffins of the governor were nearly intact with the exception of the head end that had been removed by tomb robbers. The intricate carvings and paintings on the governor’s outer coffin make it an unparalleled masterpiece of MK art (see pictures in [5]).
When the tomb was looted, the mummies of the governor and his wife were damaged by thieves in search of fine jewelry. A torso, originally attributed to Lady Djehutynakht, was found in the far corner of the burial chamber, but was recently argued to belong to the governor after re-excavation of Tomb 10A [6]. A mummified head, which could neither be attributed to the governor nor to his wife, was found atop the governor’s coffin (Figure 1).
In an effort to learn more, the head was analyzed with computerized tomography (CT) in 2005. The CT scanning of the head revealed extensive bilateral post-mortem alterations of the facial bones [7]. The absence of these bones, together with the lack of comparative data on ancient Egyptian skulls, preclude definitive morphological sex determination; however, the presence of large mastoid processes, robust occipital and temporal regions, and pronounced gonial flaring of the mandible, suggest that the skull more likely belonged to a man [8], (Figure S1). In order to unequivocally determine the biological sex of the individual, the MFA collaborated with the Federal Bureau of Investigation (FBI) Laboratory to perform DNA analysis.
At the time the FBI was contacted, the ancient DNA community had largely given up on the testing of ancient Egyptian human remains. Though DNA extraction and amplification from ancient Egyptian samples had been attempted in the early days of paleogenetics, these initial attempts either resulted in failure (4500-year-old human femurs [9]) or yielded data that turned out to be the product of modern DNA contamination (2400-year-old mummy [10] and 1600-year-old sacred monkey bones from the Saqqara Baboon Galleries [11]), both contaminated with modern human DNA). These early failures prompted studies on DNA survival [9,12]. Together with older research showing that DNA degradation (depurination in particular) is primarily influenced by temperature, pH, oxygen, and water [13], these analyses suggested that Egyptian environmental conditions likely cause DNA to degrade to fragments smaller than 100 base pairs (bp) in just a few hundred years [14,15].
Currently, forensic DNA testing at the FBI, as well as in nearly all global operational laboratories, is based on targeted PCR amplification of fragments ranging from 90–1200 bp in size, followed by size-based capillary electrophoresis for short tandem repeats (STRs) or Sanger sequencing for mitochondrial DNA (mtDNA). These approaches have served the forensic community well over the past twenty years. They present limitations, however, to both the quality and quantity of genetic data that can be recovered from the most challenging specimens. Recently, a number of commercial high throughput sequencing (HTS) assays, designed specifically for forensic applications, have become available [16,17]. While these assays overcome many of the limitations of traditional capillary electrophoresis-based forensic DNA analyses, they are still based on targeted amplification of defined genomic regions. As a result, their utility is limited to samples harboring DNA fragments large enough for PCR amplification.
One of the primary advantages of HTS is that DNA fragments of very short size can be recovered and sequenced, thus obviating the need for targeted PCR. Here, we exploit this feature and describe the use of shotgun sequencing to determine the biological sex of a 4000-year-old Egyptian mummy.

2. Materials and Methods

2.1. Tooth Extraction

In 2009, a second mandibular molar, in pristine condition, was extracted from the skull using a flexible fiber optic endoscope with grasping forceps [18]. Immediately following its extraction, the molar was placed in a sterile container and sent to a DNA laboratory where it was embedded in epoxy and cut in two parts (see photo in Figure 2). The block containing the crown and upper roots (~2/3 of the tooth) was sent to the FBI Laboratory and later to the Harvard Medical School (HMS) ancient DNA facility.

2.2. Strategy for Analyses

Given our understanding that the recovery of DNA from a 4000-year-old Egyptian mummy specimen would be challenging, if not impossible, the first goal was to determine whether any endogenous mitochondrial DNA could be recovered from the sample. On a practical level, the high copy number of the mitochondrial genome (mtGenome) per cell significantly improves the chances of recovering genetic material from ancient samples. Furthermore, because consistently higher HTS coverage is achieved with mtDNA than with nuclear DNA, the authenticity of the genetic data can be more easily addressed. Data authenticity was assessed in three ways. First, the average size of the human mtDNA sequences was determined bioinformatically. Given the ancient Egyptian origin of the sample, the recovered fragments were expected to be very small as a result of DNA degradation over time, while modern contamination would be more likely to contain longer fragments. Second, the occurrence and pattern of one particular type of DNA damage was evaluated. Hydrolytic deamination results in the conversion of cytosine to uracil in the DNA. It is one of the most abundant forms of damage in aged DNA molecules and manifests in the final DNA sequence data as conversion of cytosine to thymine (C-T) or, in the reverse complementary sequence, guanine to adenine (G-A), especially at the ends of molecules [19,20,21]. Finally, the rate of modern DNA contamination in the data was estimated by comparing the data to modern human mtDNA diversity. As mtDNA is a haploid marker, contaminating molecules can be more easily detected than is possible with diploid nuclear DNA.

2.3. Sample Preparation

A full description of the sample preparation is provided in the Supplementary Material.
At the FBI Laboratory, DNA was extracted from 105 mg of dentin powder using a silica column-based protocol [22]. In order to remove damaged bases, the extract and associated extraction reagent control (RB) were treated for one hour with the USERTM kit, which contains a mixture of Uracil DNA Glycosylase (UDG) and Endo VIII (New England BioLabs, a.k.a. NEB, Ipswich, MA, USA). At the library preparation stage, a third sample, containing only water, was introduced (negative control or NC). Illumina libraries were generated for each sample (Lib1, Lib1/RB, Lib1/NC) using the NEBNext Ultra II for Illumina kit and NEB looped adaptors. After amplification and indexing, all three libraries were subjected to hybridization capture to enrich for human mtDNA with a MyBait1 kit (Arbor Biosciences, Ann Arbor, MI, USA) and the post-capture products (cap-Lib1; cap-Lib1/RB; cap-Lib1/NC) were sequenced on an Illumina MiSeq FGX at the FBI Laboratory (Figure 2). Having established the presence of endogenous DNA at the FBI Laboratory, the tooth was sent to the HMS, where a second extract and library were prepared.
At the HMS, DNA was extracted from 68 mg of powder using the protocols published in [23] and [24]. The extract, and associated RB and NC, were treated with UDG according to [25]. By using this treatment, called “partial UDG removal”, almost all the uracil present in the ancient DNA molecules are removed except a few located at the ends of the molecule. Next, three Illumina libraries (Lib2, Lib2/RB, and Lib2/NC) were prepared according to [26]. Fifteen microliters of barcoded Lib2, and 10 µL of Lib2/RB and Lib2/NC were shipped to the FBI Laboratory where 5 μL of each library was dual-indexed with eight cycles of PCR. At the FBI Laboratory, mtDNA hybridization capture was performed on the three Harvard libraries using a MyBait1 kit. Cap1-Lib2/RB and cap1-Lib2/NC were sequenced on the MiSeq FGX while cap-Lib1 and cap1-Lib2 were sequenced together on a HiSeq 2500 at the National Bioforensic Analysis Center (NBFAC). A second mtDNA enrichment on Lib2, Lib2/RB and Lib2/NC was performed at the HMS using a pool of oligonucleotides, synthesized on a microarray (CustomArray Inc., Bothell, WA, USA). Cap2-Lib2, cap2-Lib2/RB, and cap2-Lib2/NC were sequenced at the HMS on a NextSeq 500.
Following confirmation of the mitochondrial DNA data, Lib1 and Lib2 were both shotgun-sequenced on a HiSeq 2500 at NBFAC.

2.4. Bioinformatics

Libraries sequenced on the HiSeq or the MiSeq FGX were analyzed at NBFAC. After conversion from BCL to FASTQ format, the HiSeq files were de-multiplexed and the adaptors removed. CutAdapt (v1.9.1; [27]) was used to remove the 7 bp barcodes in Lib2 and cap1-Lib2. For every run except the NextSeq run, paired reads were mapped to the human genome hg19 and the revised Cambridge reference sequence rCRS (NC012910; [28]) using the Burrows–Wheeler Aligner (BWA aln v.0.7.13-r1126; [29]) with the parameters recommended by [30] (“−l 16500 −n 0.01 −o 2”). SAMtools (v1.3.1; [31]) was used to extract mapped reads, merge the BAM files from all HiSeq lanes, and remove mapped reads with a quality score (Q) less than 30. The Picard program v.1.96 [32]) was used to remove duplicates. Data for cap2-Lib2, cap2-lib2/RB and cap2-lib2/NC were analyzed at HMS according to [33]. BAM files of all three alignments were sent to the FBI. At the FBI, all mappings were imported into the CLC Genomics Workbench program (v.10.0.1; Qiagen/CLC bio, Aarhus, Denmark) for visualization. Finally, variants were called using the CLC Fixed Ploidy tool.

2.4.1. Post-Capture Mitochondrial DNA Sequence Analysis

Sequence data produced with cap-Lib1 on two Illumina instruments (runs I and II) were merged into a single CAP-LIB1 file. Similarly, data produced with cap1-Lib2 and cap2-Lib2 (runs III and IV) were merged into a CAP-LIB2 file (Figure 2). For all data sets, duplicates were removed after merging. Reads shorter than or equal to 35 bp were removed from all data sets to avoid spurious microbial alignments (as suggested by [34]) and reads longer than 70 bp (that could potentially originate from modern contaminants) were removed before variant calling. Data observed in the reagent blanks and negative controls were mapped and analyzed using the same parameters.
To assess data authenticity, C-T and G-A errors due to deamination were quantified using MapDamage 2.0 [35]. Modern DNA contamination was estimated using ContamMix [36] and by evaluating the number of human sequences observed in the extraction (RB) and library preparation control (NC) samples.
The mtDNA haplogroup of the final profile was determined using Phylotree build 17 [37] and HaploGrep 2.0 [38]. Finally, the sequence was compared to other mtDNA sequences from modern and ancient populations stored in GenBank, the European Nucleotide Archives database of the EMBL (ENA) and EMPOP (EDNAP mtDNA population database).

2.4.2. Shotgun Sequence Analysis for Biological Sex Determination

Two established methods that have been successfully applied in a number of previous studies [39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56] were used to determine the sex of the individual from whom the tooth was obtained. Biological sex was first determined using RY [30], which is defined as the number of reads that mapped to chromosome Y (nY) divided by the number of reads mapped to both the X and Y chromosomes (nX + nY). The 95% confidence interval (CI) is defined as RY ± 1.96 standard error (SE), where SE = √[(RY x (1 − RY))/nX + nY]. The results are indicative of a female when the upper bound of the 95% CI is less than 0.016, whereas the results are indicative of a male when the lower bound of the CI is greater than 0.075.
The second method, RX [57], compares the number of sequences originating from the X chromosome to the number of sequences originating from the 22 autosomes. Molecular sex is assigned male if the upper bound of the 95% CI is less than 0.60, and female if the lower bound is greater than 0.80 (see RX CI definition in the Supplementary Material).
In order to mitigate against the possible impact of modern contaminants on the calculation, RX and RY were also determined using only reads showing signs of damage, as these reads are unlikely to have originated from modern contaminants [58]. Read filtering was performed using PMDtools (v0.50; [58]) with a threshold of three, as recommended by the author of the program.

3. Results

3.1. Data Authentication Based on Captured mtDNA Reads

Sequencing results for all four post-capture library runs are presented in Supplementary Table S1. The size distributions of the mapped mtDNA reads from the merged datasets (CAP1-LIB1 and CAP2-LIB2) are presented in Figure 3a,b. The vast majority of reads from the FBI library (CAP1-LIB1) were between 35 and 80 bp, with a mode around 47 bp. The majority of the HMS library reads, on the other hand, were between 25 and 70 bp, with a mode around 38 bp. The shorter average length of the captured mtDNA sequence in the Harvard library (Figure 3b) reflects the efficiency of the extraction protocol from [23] in retaining the smallest molecules.
The final numbers of unique reads (Q score >30; size ranging from 36 to 70 bp) were 14,408 for CAP-LIB1 and 26,502 for CAP-LIB2 (40,910 reads in total). The total number of mtDNA bases sequenced was ~1.8 million (703,905 for CAP-LIB1 and 1,096,128 for CAP-LIB2). When aligned to the mtGenome, coverage of any given position ranged from 5× to 247×, with an average coverage of 108× (Figure 4).

3.2. Quantification of Deamination

Despite the fact that DNA damage resulting from deamination slowly accumulates over time, deamination does not strictly correlate with sample age. Environmental conditions, in addition to age, can significantly impact DNA damage. High temperatures and humidity generally increase the speed of degradation, while cold and dry climates generally favor preservation. A 22,000-year-old bone preserved in permafrost can exhibit a C-T substitution frequency at the 5’ end of less than 0.20 [59], while a much younger, 6500-year-old human tooth found in Spain can exhibit rates greater than 0.25 [44]. Even if the C-T substitution frequency cannot be used as a strict dating tool, deamination rates are generally expected to be greater than 0.1 in unrepaired DNA extracted from samples greater than 500 years of age that have not been preserved in ice [60].
In the present study, the C-T substitution frequency at the 5’ end of the mtDNA reads in Cap-Lib1 was 0.27 (Figure S2A). This high value was somewhat unexpected since the extract had been treated with UDG (an enzyme that cleaves uracil) and only sporadic uracil bases should have remained. These results suggest that the reaction was partially or completely inhibited. Conversely, and as expected for the partial UDG treated cap1-Lib2, very few uracil bases remained in the mtDNA reads, resulting in a nucleotide misincorporation rate of 0.084 at the 5’ end (Figure S2C). Since it has been shown in [25] that a successful “partial UDG treatment” generally results in the misincorporation rate to be reduced by three-fold, it is likely that our USER reaction performed on Lib1 was completely inhibited.

3.3. Contamination

Sequencing statistics for the reagent blank and negative control libraries are presented in Table 1. No signs of measurable human contamination that could have impacted results from the tooth were observed. The vast majority of reads did not map to the human genome and most likely originate from environmental DNA that was present in the tooth, the laboratory, and/or the reagents.
Modern human contamination in the mtDNA data, as assessed using ContamMix, was estimated at 0.6% (95% CI 0.08–1.8%) in CAP-Lib1 and 11.6% in CAP-Lib2 (95% CI 9.7–13.6%). The value for Lib1 is consistent with the rates observed in other ancient DNA studies [44,46,48,49,50,51,52,61,62,63,64,65,66]. Though the estimated contamination rate in Lib2 was somewhat higher, it had no effect on consensus calling for the endogenous molecules, as the same mtGenome profile was recovered from both libraries.

3.4. Mitochondrial Haplotype

The mtGenome profile independently obtained from the tooth by the FBI and HMS laboratories were identical and can be found in Table S2. The haplotype (deposited in GenBank under accession number MG736653) belongs to mitochondrial DNA lineage U5b2b5, but the specific sequence has not been previously reported in the 35,942 mtGenomes stored in the NCBI GenBank database (as of October 2017). The sequence closest to the mummy’s belongs to a contemporary individual from Lebanon (KT779192 [67]); however, the two haplotypes still differ at five positions, three of them in the control region (CR). A comparison between the mummy CR and the 26,127 CR sequences from the EMPOP database produced no match.
To better understand the mtDNA lineage of the mummy in the context of known Egyptian mtDNA diversity, the mummy haplogroup was compared to the mtDNA haplogroup distribution of 668 Egyptians from various modern populations [68,69,70,71,72,73]. The dominant haplogroups among this dataset were haplogroup T (11.98%) and L3 (11.23%; Table S3). Out of the 64 individuals who belonged to haplogroup U, seven belonged to haplogroup U5 (1.05%), and three (0.5%) belonged to one of the U5b subgroups (U5b1c; U5b1d1a; U5b2a5).
The Djehutynakht sequence was also compared to available ancient human DNA sequences (Table S4). Not surprisingly, no direct matches to the Djehutynakht sequence have been reported. However, related U5b2b sequences have been observed in ancient human remains from Europe, and a haplogroup U5b2c1 haplotype was recently discovered in 2000-year-old remains from Phoenicia [67]. When only the mtDNA sequences recovered from ancient Egyptian human remains are considered, the Djehutynakht sequence most closely resembles a U5a lineage from sample JK2903, a 2000-year-old skeleton from Abusir el-Meleq [74].

3.5. Shotgun Sequencing

Shotgun sequencing statistics are presented in Table 2. Although aliquots of libraries Lib1 and Lib2 were mixed (50/50) and sequenced together on the HiSeq, the Lib2 data values were, for all measures, substantially higher than the Lib1 values. This is perhaps due to the smaller overall size of the HMS library fragments (Figure S3). As smaller fragments are known to be preferentially sequenced on Illumina platforms, the smaller average size of the HMS library likely explains the greater numbers of raw paired reads, reads that mapped to the human genome, and unique mapped reads.
It is also likely that the greater fraction of human DNA recovered from the HMS library (6.57% versus 2.24% in the FBI library) is a direct result of the increased recovery of smaller fragments.
Due to the extremely degraded state of the endogenous DNA, more of it was likely recovered in the HMS library.
Following the removal of duplicates and reads with low mapping quality scores, Lib1 yielded 1,595,239 reads (1,593,816 nuclear sequences and 1423 mtDNA sequences), while Lib2 yielded 7,691,326 reads (7,687,370 nuclear reads and 3956 mtDNA reads; Table S4). The percentage of mtDNA reads in each of the two libraries was <0.1%, in contrast to the mtDNA capture libraries, where it was 10.9% in Cap-Lib1 and 19.7% in Cap1-Lib2. The average coverage of the mtGenome using shotgun data was too low to produce a full profile. However, when variants could be called (minimum of five reads and a frequency ≥80%), they were in agreement with the profile produced by the hybridization capture data. The average coverage over the entire human genome hg19 was 0.02× for Lib1 and 0.09× for Lib2.

3.6. Biological Sex Determination

Both RY and RX were first calculated using all reads with Q-scores greater than 30 that aligned to the human genome (Table 3; for mapping details, see Table S5). Since the contamination rate in Lib2 was somewhat high (>11%), RY and RX were also calculated using only reads exhibiting signs of DNA damage (i.e., molecules assumed to be endogenous) (Table 4). For both shotgun libraries (Lib1 and Lib2), the calculated lower bound of the RY confidence interval was >0.075, regardless of whether all reads or only reads showing signs of deamination were used. Values greater than >0.075 point to a male. Similarly, and again regardless of whether or not undamaged molecules were included in the calculations, data from both Lib1 and Lib2 produced an upper bound for the RX confidence interval that was less than 0.6. For this calculation, values less than 0.6 are consistent with a male assignment.
The DNA length distribution of the shotgun data for Lib1 showed a spike at 35 bp, possibly indicative of spurious microbial alignments (Supplementary Figure S3A; [34]). Therefore, as an extra precaution, RX and RY were also calculated using only damaged reads >35 bp (467,923 reads). The results were as follows: RX = 0.4369 (95% CI 0.4134–0.4603) and RY = 0.0879 (95% CI 0.0827–0.093). The male sex assignment remained the same.

4. Discussion

4.1. Origins of the Mummy

For nearly a century, the contents of Tomb 10A have been studied by archaeologists, philologists, and art historians. Based on the styles of pottery, the many wooden sculptures and, of course, the decorations and texts on governor Djehutynakht’s coffins, the tomb has been firmly dated to within a generation between the end of the 11th and the beginning of the 12th Dynasties (1961–2010 BC; [3,75]). The head is clearly original to the tomb. While, for genetic testing purposes, carbon dating of the specimen would have been ideal to further establish its age, additional manipulations were unfortunately impossible, given the value of the head as an object of art. However, several features of both the head and the tomb point clearly to the head originating from one of the original tomb occupants. For one, there is no evidence that the tomb was later reused despite it having been looted in antiquity. No artifacts or tomb contents dating to a later period have been found; in fact, additional human bones found in one of the shafts of Tomb 10A during a more recent 2009 expedition at Deir el-Bersha still support the conclusion that the chamber was used as the resting place for only two individuals [6]. The preparation of the head itself is consistent with the dating of the artifacts in the tomb (1961–2010 BC). The head displays rather archaic mummification practices, namely, facial features that were modeled in linen and then painted. This type of treatment was common during the Old Kingdom, but progressively disappeared during the MK [76,77,78]. Indeed, the eyes and lips of Djehutynakht’s head were padded and his eyebrows painted (Figure 1). Similar features have been observed in another Deir-el Bersha mummy slightly predating Djehutynakht’s era. The mummy of Henu, recovered from a completely undisturbed late First Intermediate Period (~2100–2040 BC) tomb, exhibits the same characteristics regarding facial modeling [79]. As these particular mummification practices did not persist through the MK and later, this feature supports the position that the Tomb 10A head is contemporaneous with the rest of the tomb’s early MK contents [76]. All told, from art history, anthropological, and archaeological perspectives, the facts point to the mummy head belonging to one of the original Tomb 10A occupants: a MK nomarch of the late 11th or early 12th Dynasty or his wife.

4.2. Eurasian mtDNA Haplogroups in Ancient Egyptians

At the time DNA testing was performed on the tooth, and for reasons previously discussed, very little had been published on DNA recovery from ancient Egyptian human remains. Only one publication including HTS and quality control measures was available in early 2016, which described the mtGenome sequencing of an Egyptian mummy from the Greco-Roman period. The individual belonged to mtDNA haplogroup I2 [80]. Two other studies describing mtDNA recovery from ancient African samples were also available at the time, but centered on skeletons from more southern regions of the continent. One described the L0d2c1c lineage mtGenome of a 2330-year-old male skeleton from South Africa [81], while the other described the recovery of a L3x2a mtGenome from the remains of a 4500-year-old individual from Ethiopia [82].
Given limited available data and the fact that U5 is the dominant mitochondrial haplogroup found among hunter-gatherers in Europe [83,84], the recovery of a haplogroup U5b2b5 sequence from the mummy of Djehutynakht raises the question of data authenticity, despite the molecular metrics suggesting otherwise. When the mummy’s mtDNA sequence is viewed in the context of modern mtDNA diversity, however, the observed U5 lineage could potentially reflect interactions between Egypt and the Near East that date as far back as the Predynastic and Early Dynastic periods [85]. Trade between Egypt and the Near East is evidenced by, among other things, ceramic imports to Egypt [86]. In addition, dwellings similar to those found in Palestine suggest some immigration to Egypt from more arid Near Eastern areas from the late Predynastic to the Old Kingdom [85,87]. Both trade and immigration between Egypt and the Near East continued to increase over time. Demand in Egypt for cedar of Lebanon wood (a wood available and harvested in Lebanon and Syria during the MK) led to the further establishment of trade routes between Egypt and the Levant [85,86]. It is interesting, and perhaps not coincidental, that the individual with the mtDNA sequence most similar to Djehutynakht comes from a Lebanese individual.
On top of this historical information offering an explanation for the observed mtDNA data are now additional, recently published, mtGenomes from Africa, and Egypt in particular. MtDNA haplotypes recently obtained from ancient human remains from sub-Saharan Africa belong only to haplogroup L subgroups [65,88]. However, nearly all of the remains excavated in the Northern part of the continent belong to Eurasian mtDNA lineages [63,67,74,89,90]. In fact, of the 114 mtDNA genomes now available from northern African ancient human remains, only one belongs to an African lineage (L3 observed in a skeleton from Abusir el-Meleq [74]). The deep presence of Eurasian mtDNA lineages in Northern Africa has, therefore, been clearly established with these recent reports and offers further support for the authenticity of the Eurasian mtDNA sequence observed in the Djehutynakht mummy. In the present study, Near Eastern influence has been found in an individual of high social status who lived in Upper Egypt during the Middle Kingdom.

4.3. Perspectives for Forensic Laboratories

All currently employed capillary electrophoresis-based human DNA identification methods, as well as recently released commercial HTS assays developed for forensic applications, are based on targeted PCR approaches requiring endogenous DNA fragments, sometimes as small as 70 bp [91] but generally >100 bp. In forensic cases involving the most limited and degraded specimens (e.g., burnt bones, touch DNA, or single hair shafts), the DNA is often too damaged to yield amplicons of that size. Even in those instances for which DNA typing is successful with currently employed technologies, results are generally purposely limited to small regions of the mitochondrial DNA control region. Previous DNA analyses of the Djehutynakht tooth, prior to its receipt by the FBI, are consistent with these trends. Testing performed on the root end of the tooth (Figure 2) between 2009 and 2012 by two different laboratories that used traditional PCR approaches and amplicons as small as 80 bp failed.
The utility of capture-based and shotgun HTS approaches for these types of difficult samples has been recognized for some time in forensics. A number of studies have demonstrated the utility of HTS in recovering DNA from samples recalcitrant to standard PCR. DNA has been recovered from severely degraded human remains [92,93,94,95] as well as single shed hairs [96]. All work to date, however, has focused exclusively on the mitochondrial genome. The mtDNA results produced here further support the use of HTS techniques for the recovery of accurate and reliable complete mtGenome data from severely degraded specimens, and from capture-based data in particular. However, the data described here also expand on these mtDNA-only approaches by demonstrating the recovery of authentic data from the much more discriminating nuclear genome.
Nuclear DNA testing in forensics is currently based on autosomal short tandem repeat (STR) loci, and, thus, STRs would be the logical target marker for individual identification purposes. Unfortunately, due to their relatively large size, successful targeted amplification of STRs is rare with highly degraded samples. In addition, STRs are difficult to recover from shotgun data in these types of specimens. Not only is the full length of the repeat region rarely intact, but small recovered fragments of repetitive sequence are also unlikely to include enough genome-specific flanking regions to allow for proper mapping. Single nucleotide polymorphisms (SNPs), however, can be used for individual identification [97,98,99], and these could be typed using hybridization capture and HTS. SNPs are already being targeted by the ancient DNA community to determine Y haplogroups [74,88,100], assess kinship [89], and test models of human population history and evolution [41,61,62,64,65,66,100,101,102]. For these types of questions, thousands of nuclear SNPs are often targeted, and low coverage of even a subset is generally sufficient to answer the question at hand.
For forensics and other regulated disciplines (e.g., clinical genetics), where questions pertain specifically to individual genetic profiles, strict metrics for profile accuracy and reliability must be met. For example, adequate depth of coverage must be established to ensure SNP genotyping accuracy, and proper analytical thresholds must be set to capture stochasticity in allelic sampling, laboratory processes, and data analysis workflows [103,104]. These data quality metrics must furthermore be met with time-efficient and cost-effective workflows. Thus, while proof of principle studies in the ancient DNA and, now, forensics disciplines demonstrate recovery of endogenous nuclear DNA from severely compromised specimens, more developmental work is clearly required to enable standard forensic application of these techniques.

5. Conclusions

Over the past year, molecular techniques developed and routinely used by the ancient DNA community have finally permitted the recovery of endogenous DNA from ancient Egyptian remains. Here, those techniques were employed to recover the complete mtGenome of the 4000-year-old mummy, Djehutynakht, and to determine that the biological sex was male. Although these approaches have been recently adopted and implemented by the forensic community to develop probative mitochondrial DNA data from the most degraded specimens, further work is needed for the recovery of individually identifiable nuclear DNA markers. Assays will need to be optimized and cost-effective workflows developed to achieve accurate and reliable calls from limited quantities of damaged nuclear DNA. In the meantime, the work described here represents a first step towards further improving forensic DNA testing capabilities in such cases.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/9/3/135/s1. Figure S1: 3D views of the mummy’s skull; Figure S2: MapDamage analyses; Figure S3: DNA length distribution in shotgun sequencing data. Table S1: Sequencing statistics for post capture libraries; Table S2: Mitochondrial DNA profile of the Djehutynakht mummy; Table S3: Haplogroup distribution among 668 modern Egyptians; Table S4: U5b2b individuals in ancient populations; Table S5: Mapped read distributions used to determine RX and RY.

Acknowledgments

Above all, we would like to thank Rita Freed and the Museum of Fine Arts for entrusting us with the analysis of such a precious sample and for providing valuable historical and contextual information. We would also like to express our gratitude to the following people: Marleen De Meyer, for her constant assistance in educating us about Egyptology and for her input on the manuscript; Jeroen van der Perk, for assistance with bioinformatics analyses; Nicole Adamski, Nasreen Broomandkhoshbacht, Matthew Ferry, Ann-Marie Lawson, Megan Michel, Jonas Oppenheimer, Kristin Stewardson, Zhao Zhang, Kimberly Eskey, and Sana Enke, for their technical assistance; Richard Thomas, for studying the mummy CT scans and for his assistance with anthropology jargon; Brandon C. Letts, Rebecca Just, and Lilliana Moreno for their help with the manuscript or the figures; and finally, Raj Gupta, Paul Chapman, Fabio Nunes, and Marina Faerman for their support and for being so responsive to our multiple inquiries. This work was funded under Contract No. HSHQDC-15-C-00064 awarded by the Department of Homeland Security (DHS) Science and Technology Directorate (S&T) for the operation and management of the National Biodefense Analysis and Countermeasures Center (NBACC), a Federally Funded Research and Development Center.

Author Contributions

O.L. and J.A.I. conceived and designed the experiments; and N.H.B. and A.J.O. contributed reagents and materials. O.L. performed the experiments at the FBI. N.R., S.M., and D.R. supervised the laboratory and bioinformatics work at HMS. S.R., A.L.B., D.D.S., T.B.S., P.S., P.L.F.J., S.M., and O.L. processed and analyzed the data. J.A.I. and O.L. wrote the paper with input from A.L.B., P.L.F.J., and N.R. All authors read and approved the manuscript.

Conflicts of Interest

Names of commercial manufacturers are provided for identification purposes only and inclusion does not imply endorsement of the manufacturer, or its products or services by the FBI. The views expressed are those of the authors and do not necessarily reflect the official policy or position of the FBI or the U.S. Government. This is FBI Laboratory publication #18-13. The views and conclusions contained in this document are those of the author and should not be interpreted as necessarily representing the official policies, either expressed or implied, of the DHS or S&T. In no event shall DHS, NBACC, S&T, or Battelle National Biodefense Institute have any responsibility or liability for any use, misuse, inability to use, or reliance upon the information contained herein. DHS does not endorse any products or commercial services mentioned in this publication.

References

  1. Freed, R.E.; Berman, L.M.; Doxey, D.M.; Picardo, N.S. The Secrets of Tomb 10A: Egypt 2000 BC; MFA Publications: Boston, MA, USA, 2009. [Google Scholar]
  2. Willems, H.O. The nomarchs of the Hare nome and early Middle Kingdom history. In Proceedings of the Fourth International Congress of Egyptology, Munich, Germany, 26 August–1 September 1985. [Google Scholar]
  3. De Meyer, M.; Dils, P. Fowl for the Governor: the tomb of governor Djehutinakht IV or V at Dayr al-Barshā reinvestigated. Part I. J. Egyptian Archaeol. 2012, 98, 55–72. [Google Scholar] [CrossRef]
  4. Roth, A.M.; Roehrig, C.H. The Bersha procession: A new reconstruction. J. Mus. Fine Arts 1989, 1, 31–40. [Google Scholar]
  5. Museum of Fine Arts, Boston. Available online: http://www.mfa.org/collections/object/front-side-panel-of-outer-coffin-of-djehutynakht-142815 (accessed on 26 February 2018).
  6. De Meyer, M.; Linseele, V.; Vereecken, S.; Williams, L.J. Fowl for the governor. The tomb of governor Djehutinakht IV or V at Dayr al- Barsha reinvestigated. Part 2: Pottery, human remains, and faunal remains. J. Egyptian Archaeol. 2014, 100, 67–87. [Google Scholar] [CrossRef]
  7. Gupta, R.; Markowitz, Y.; Berman, L.; Chapman, P. High-resolution imaging of an ancient Egyptian mummified head: New insights into the mummification process. Am. J. Neuroradiol. 2008, 29, 705–713. [Google Scholar] [CrossRef] [PubMed]
  8. Thomas, R., Anthropologist at the Federal Bureau of Investigation. Personal communication, March 2017.
  9. Poinar, H.N.; Hoss, M.; Bada, J.L.; Pääbo, S. Amino acid racemization and the preservation of ancient DNA. Science 1996, 272, 864–866. [Google Scholar] [CrossRef] [PubMed]
  10. Pääbo, S. Molecular cloning of Ancient Egyptian mummy DNA. Nature 1985, 314, 644–645. [Google Scholar] [CrossRef] [PubMed]
  11. Van der Kuyl, A.C.; Dekker, J.; Attia, M.A.M.; Iskander, N.; Perizonius, W.R.K.; Goudsmit, J. DNA from ancient Egyptian monkey bones. Ancient DNA Newsletters 1994, 2, 19–20. [Google Scholar]
  12. Smith, C.I.; Chamberlain, A.T.; Riley, M.S.; Stringer, C.; Collins, M.J. The thermal history of human fossils and the likelihood of successful DNA amplification. J. Hum. Evol. 2003, 45, 203–217. [Google Scholar] [CrossRef]
  13. Lindahl, T.; Nyberg, B. Rate of depurination of native deoxyribonucleic acid. Biochemistry 1972, 11, 3610–3618. [Google Scholar] [CrossRef] [PubMed]
  14. Marota, I.; Basile, C.; Ubaldi, M.; Rollo, F. DNA decay rate in papyri and human remains from Egyptian archaeological sites. Am. J. Phys. Anthropol. 2002, 117, 310–318. [Google Scholar] [CrossRef] [PubMed]
  15. Gilbert, M.T.; Barnes, I.; Collins, M.J.; Smith, C.; Eklund, J.; Goudsmit, J.; Poinar, H.; Cooper, A. Long-term survival of ancient DNA in Egypt: Response to Zink and Nerlich. Am. J. Phys. Anthropol. 2005, 128, 110–114. [Google Scholar] [CrossRef] [PubMed]
  16. ForenSeqTM DNA Siganture Prep kit (TG-450–1001); Illumina, Inc.: San Diego, CA, USA, 2017.
  17. Precision ID Kits (Identity panel; ancestry panel; mtDNA panel); ThermoFisher: Waltham, MA, USA, 2017.
  18. Nunes, F., Assistant Professor of Clinical Pediatrics; Indiana University School of Medicine. Personal communication, January 2017.
  19. Hofreiter, M.; Jaenicke, V.; Serre, D.; von Haeseler, A.; Pääbo, S. DNA sequences from multiple amplifications reveal artifacts induced by cytosine deamination in ancient DNA. Nucleic Acids Res. 2001, 29, 4793–4799. [Google Scholar] [CrossRef] [PubMed]
  20. Stiller, M.; Green, R.E.; Ronan, M.; Simons, J.F.; Du, L.; He, W.; Egholm, M.; Rothberg, J.M.; Keates, S.G.; Ovodov, N.D.; et al. Patterns of nucleotide misincorporations during enzymatic amplification and direct large-scale sequencing of ancient DNA. Proc. Natl. Acad. Sci. USA 2006, 103, 13578–13584. [Google Scholar] [CrossRef] [PubMed]
  21. Briggs, A.W.; Stenzel, U.; Johnson, P.L.; Green, R.E.; Kelso, J.; Prufer, K.; Meyer, M.; Krause, J.; Ronan, M.T.; Lachmann, M.; et al. Patterns of damage in genomic DNA sequences from a Neandertal. Proc. Natl. Acad. Sci. USA 2007, 104, 14616–14621. [Google Scholar] [CrossRef] [PubMed]
  22. Loreille, O.M.; Parr, R.L.; McGregor, K.A.; Fitzpatrick, C.M.; Lyon, C.; Yang, D.Y.; Speller, C.F.; Grimm, M.R.; Grimm, M.J.; Irwin, J.A.; et al. Integrated DNA and fingerprint analyses in the identification of 60-year-old mummified human remains discovered in an Alaskan glacier. J. Forensic Sci. 2010, 55, 813–818. [Google Scholar] [CrossRef] [PubMed]
  23. Dabney, J.; Knapp, M.; Glocke, I.; Gansauge, M.T.; Weihmann, A.; Nickel, B.; Valdiosera, C.; Garcia, N.; Pääbo, S.; Arsuaga, J.L.; et al. Complete mitochondrial genome sequence of a Middle Pleistocene cave bear reconstructed from ultrashort DNA fragments. Proc. Natl. Acad. Sci. USA 2013, 110, 15758–15763. [Google Scholar] [CrossRef] [PubMed]
  24. Korlevic, P.; Gerber, T.; Gansauge, M.T.; Hajdinjak, M.; Nagel, S.; Aximu-Petri, A.; Meyer, M. Reducing microbial and human contamination in DNA extractions from ancient bones and teeth. Biotechniques 2015, 59, 87–93. [Google Scholar] [CrossRef] [PubMed]
  25. Rohland, N.; Harney, E.; Mallick, S.; Nordenfelt, S.; Reich, D. Partial uracil-DNA-glycosylase treatment for screening of ancient DNA. Phil. Trans. R. Soc. Lond. B Biol. Sci. 2015, 370, 20130624. [Google Scholar] [CrossRef] [PubMed]
  26. Rohland, N.; Reich, D. Cost-effective, high-throughput DNA sequencing libraries for multiplexed target capture. Genome Res. 2012, 22, 939–946. [Google Scholar] [CrossRef] [PubMed]
  27. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal 2011, 17, 10–12. [Google Scholar] [CrossRef]
  28. Andrews, R.M.; Kubacka, I.; Chinnery, P.F.; Lightowlers, R.N.; Turnbull, D.M.; Howell, N. Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat. Genet. 1999, 23, 147. [Google Scholar] [PubMed]
  29. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [PubMed]
  30. Skoglund, P.; Storå, J.; Götherström, A.; Jakobsson, M. Accurate sex identification of ancient human remains using DNA shotgun sequencing. J. Archaeol. Sci. 2013, 40, 4477–4482. [Google Scholar] [CrossRef]
  31. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. The sequence alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed]
  32. Picard. Available online: http://broadinstitute.github.io/picard (accessed on 26 February 2018).
  33. Meyer, M.; Kircher, M.; Gansauge, M.T.; Li, H.; Racimo, F.; Mallick, S.; Schraiber, J.G.; Jay, F.; Prufer, K.; de Filippo, C.; et al. A high-coverage genome sequence from an archaic Denisovan individual. Science 2012, 338, 222–226. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Meyer, M.; Arsuaga, J.L.; de Filippo, C.; Nagel, S.; Aximu-Petri, A.; Nickel, B.; Martínez, I.; Gracia, A.; Bermúdez de Castro, J.M.; Carbonell, E.; et al. Nuclear DNA sequences from the Middle Pleistocene Sima de los Huesos hominins. Nature 2016, 531, 504–507. [Google Scholar] [CrossRef] [PubMed]
  35. Jonsson, H.; Ginolhac, A.; Schubert, M.; Johnson, P.L.; Orlando, L. MapDamage2.0: Fast approximate Bayesian estimates of ancient DNA damage parameters. Bioinformatics 2013, 29, 1682–1684. [Google Scholar] [CrossRef] [PubMed]
  36. Fu, Q.; Mittnik, A.; Johnson, P.L.F.; Bos, K.; Lari, M.; Bollongino, R.; Sun, C.; Giemsch, L.; Schmitz, R.; Burger, J.; et al. A revised timescale for human evolution based on ancient mitochondrial genomes. Curr. Biol. 2013, 23, 553–559. [Google Scholar] [CrossRef] [PubMed]
  37. Van Oven, M. PhyloTree Build 17: Growing the human mitochondrial DNA tree. Forensic Sci. Int. Genet. Suppl. Ser. 2015, 5, e392–e394. [Google Scholar] [CrossRef]
  38. Weissensteiner, H.; Pacher, D.; Kloss-Brandstatter, A.; Forer, L.; Specht, G.; Bandelt, H.J.; Kronenberg, F.; Salas, A.; Schonherr, S. HaploGrep 2: Mitochondrial haplogroup classification in the era of high-throughput sequencing. Nucleic Acids Res. 2016, 44, W58–W63. [Google Scholar] [CrossRef] [PubMed]
  39. Gamba, C.; Jones, E.R.; Teasdale, M.D.; McLaughlin, R.L.; Gonzalez-Fortes, G.; Mattiangeli, V.; Domboroczki, L.; Kovari, I.; Pap, I.; Anders, A.; et al. Genome flux and stasis in a five millennium transect of European prehistory. Nat. Commun. 2014, 5, 5257. [Google Scholar] [CrossRef] [PubMed]
  40. Lazaridis, I.; Patterson, N.; Mittnik, A.; Renaud, G.; Mallick, S.; Kirsanow, K.; Sudmant, P.H.; Schraiber, J.G.; Castellano, S.; Lipson, M.; et al. Ancient human genomes suggest three ancestral populations for present-day Europeans. Nature 2014, 513, 409–413. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Haak, W.; Lazaridis, I.; Patterson, N.; Rohland, N.; Mallick, S.; Llamas, B.; Brandt, G.; Nordenfelt, S.; Harney, E.; Stewardson, K.; et al. Massive migration from the steppe was a source for Indo-European languages in Europe. Nature 2015, 522, 207–211. [Google Scholar] [CrossRef] [PubMed]
  42. Jones, E.R.; Gonzalez-Fortes, G.; Connell, S.; Siska, V.; Eriksson, A.; Martiniano, R.; McLaughlin, R.L.; Gallego Llorente, M.; Cassidy, L.M.; Gamba, C.; et al. Upper Palaeolithic genomes reveal deep roots of modern Eurasians. Nat. Commun. 2015, 6, 8912. [Google Scholar] [CrossRef] [PubMed]
  43. Mathieson, I.; Lazaridis, I.; Rohland, N.; Mallick, S.; Patterson, N.; Roodenberg, S.A.; Harney, E.; Stewardson, K.; Fernandes, D.; Novak, M.; et al. Genome-wide patterns of selection in 230 ancient Eurasians. Nature 2015, 528, 499–503. [Google Scholar] [CrossRef] [PubMed]
  44. Olalde, I.; Schroeder, H.; Sandoval-Velasco, M.; Vinner, L.; Lobon, I.; Ramirez, O.; Civit, S.; Garcia Borja, P.; Salazar-Garcia, D.C.; Talamo, S.; et al. A common genetic origin for early farmers from Mediterranean Cardial and Central European LBK Cultures. Mol. Biol. Evol. 2015, 32, 3132–3142. [Google Scholar] [CrossRef] [PubMed]
  45. Raghavan, M.; Steinrucken, M.; Harris, K.; Schiffels, S.; Rasmussen, S.; DeGiorgio, M.; Albrechtsen, A.; Valdiosera, C.; Avila-Arcos, M.C.; Malaspinas, A.S.; et al. Genomic evidence for the Pleistocene and recent population history of Native Americans. Science 2015, 349, aab3884. [Google Scholar] [CrossRef] [PubMed]
  46. Schroeder, H.; Avila-Arcos, M.C.; Malaspinas, A.S.; Poznik, G.D.; Sandoval-Velasco, M.; Carpenter, M.L.; Moreno-Mayar, J.V.; Sikora, M.; Johnson, P.L.; Allentoft, M.E.; et al. Genome-wide ancestry of 17th-century enslaved Africans from the Caribbean. Proc. Natl. Acad. Sci. USA 2015, 112, 3669–3673. [Google Scholar] [CrossRef] [PubMed]
  47. Kilinc, G.M.; Omrak, A.; Ozer, F.; Gunther, T.; Buyukkarakaya, A.M.; Bicakci, E.; Baird, D.; Donertas, H.M.; Ghalichi, A.; Yaka, R.; et al. The demographic development of the first farmers in Anatolia. Curr. Biol. 2016, 26, 2659–2666. [Google Scholar] [CrossRef] [PubMed]
  48. Martiniano, R.; Caffell, A.; Holst, M.; Hunter-Mann, K.; Montgomery, J.; Muldner, G.; McLaughlin, R.L.; Teasdale, M.D.; van Rheenen, W.; Veldink, J.H.; et al. Genomic signals of migration and continuity in Britain before the Anglo-Saxons. Nat. Commun. 2016, 7, 10326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Omrak, A.; Gunther, T.; Valdiosera, C.; Svensson, E.M.; Malmstrom, H.; Kiesewetter, H.; Aylward, W.; Stora, J.; Jakobsson, M.; Gotherstrom, A. Genomic evidence establishes Anatolia as the source of the European neolithic gene pool. Curr. Biol. 2016, 26, 270–275. [Google Scholar] [CrossRef] [PubMed]
  50. Skoglund, P.; Posth, C.; Sirak, K.; Spriggs, M.; Valentin, F.; Bedford, S.; Clark, G.R.; Reepmeyer, C.; Petchey, F.; Fernandes, D.; et al. Genomic insights into the peopling of the Southwest Pacific. Nature 2016, 538, 510–513. [Google Scholar] [CrossRef] [PubMed]
  51. Gonzalez-Fortes, G.; Jones, E.R.; Lightfoot, E.; Bonsall, C.; Lazar, C.; Grandal-d’Anglade, A.; Garralda, M.D.; Drak, L.; Siska, V.; Simalcsik, A.; et al. Paleogenomic evidence for multi-generational mixing between Neolithic farmers and Mesolithic hunter-gatherers in the Lower Danube basin. Curr. Biol. 2017, 27, 1801. [Google Scholar] [CrossRef] [PubMed]
  52. Hedenstierna-Jonson, C.; Kjellstrom, A.; Zachrisson, T.; Krzewinska, M.; Sobrado, V.; Price, N.; Gunther, T.; Jakobsson, M.; Gotherstrom, A.; Stora, J. A female Viking warrior confirmed by genomics. Am. J. Phys. Anthropol. 2017, 164, 853–860. [Google Scholar] [CrossRef] [PubMed]
  53. Juras, A.; Chylenski, M.; Krenz-Niedbala, M.; Malmstrom, H.; Ehler, E.; Pospieszny, L.; Lukasik, S.; Bednarczyk, J.; Piontek, J.; Jakobsson, M.; et al. Investigating kinship of Neolithic post-LBK human remains from Krusza Zamkowa, Poland using ancient DNA. Forensic Sci. Int. Genet. 2017, 26, 30–39. [Google Scholar] [CrossRef] [PubMed]
  54. Kennett, D.J.; Plog, S.; George, R.J.; Culleton, B.J.; Watson, A.S.; Skoglund, P.; Rohland, N.; Mallick, S.; Stewardson, K.; Kistler, L.; et al. Archaeogenomic evidence reveals prehistoric matrilineal dynasty. Nat. Commun. 2017, 8, 14115. [Google Scholar] [CrossRef] [PubMed]
  55. Lindo, J.; Achilli, A.; Perego, U.A.; Archer, D.; Valdiosera, C.; Petzelt, B.; Mitchell, J.; Worl, R.; Dixon, E.J.; Fifield, T.E.; et al. Ancient individuals from the North American Northwest Coast reveal 10,000 years of regional genetic continuity. Proc. Natl. Acad. Sci. USA 2017, 114, 4093–4098. [Google Scholar] [CrossRef] [PubMed]
  56. Sikora, M.; Seguin-Orlando, A.; Sousa, V.C.; Albrechtsen, A.; Korneliussen, T.; Ko, A.; Rasmussen, S.; Dupanloup, I.; Nigst, P.R.; Bosch, M.D.; et al. Ancient genomes show social and reproductive behavior of early Upper Paleolithic foragers. Science 2017, 358, 659–662. [Google Scholar] [CrossRef] [PubMed]
  57. Mittnik, A.; Wang, C.C.; Svoboda, J.; Krause, J. A Molecular Approach to the Sexing of the Triple Burial at the Upper Paleolithic Site of Dolni Vestonice. PLoS ONE 2016, 11, e0163019. [Google Scholar] [CrossRef] [PubMed]
  58. Skoglund, P.; Northoff, B.H.; Shunkov, M.V.; Derevianko, A.P.; Pääbo, S.; Krause, J.; Jakobsson, M. Separating endogenous ancient DNA from modern day contamination in a Siberian Neandertal. Proc. Natl. Acad. Sci. USA 2014, 111, 2229–2234. [Google Scholar] [CrossRef] [PubMed]
  59. Enk, J.; Rouillard, J.M.; Poinar, H. Quantitative PCR as a predictor of aligned ancient DNA read counts following targeted enrichment. BioTechniques 2013, 55, 300–309. [Google Scholar] [CrossRef] [PubMed]
  60. Sawyer, S.; Krause, J.; Guschanski, K.; Savolainen, V.; Pääbo, S. Temporal patterns of nucleotide misincorporations and DNA fragmentation in ancient DNA. PLoS ONE 2012, 7, e34131. [Google Scholar] [CrossRef] [PubMed]
  61. Cassidy, L.M.; Martiniano, R.; Murphy, E.M.; Teasdale, M.D.; Mallory, J.; Hartwell, B.; Bradley, D.G. Neolithic and Bronze Age migration to Ireland and establishment of the insular Atlantic genome. Proc. Natl. Acad. Sci. USA 2016, 113, 368–373. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Hofmanova, Z.; Kreutzer, S.; Hellenthal, G.; Sell, C.; Diekmann, Y.; Diez-Del-Molino, D.; van Dorp, L.; Lopez, S.; Kousathanas, A.; Link, V.; et al. Early farmers from across Europe directly descended from Neolithic Aegeans. Proc. Natl. Acad. Sci. USA 2016, 113, 6886–6891. [Google Scholar] [CrossRef] [PubMed]
  63. Molto, J.E.; Loreille, O.; Mallott, E.K.; Malhi, R.S.; Fast, S.; Daniels-Higginbotham, J.; Marshall, C.; Parr, R. Complete Mitochondrial Genome Sequencing of a Burial from a Romano-Christian Cemetery in the Dakhleh Oasis, Egypt: Preliminary Indications. Genes 2017, 8, 262. [Google Scholar] [CrossRef] [PubMed]
  64. Rodriguez-Varela, R.; Gunther, T.; Krzewinska, M.; Stora, J.; Gillingwater, T.H.; MacCallum, M.; Arsuaga, J.L.; Dobney, K.; Valdiosera, C.; Jakobsson, M.; et al. Genomic analyses of Pre-European conquest human remains from the Canary Islands reveal close affinity to modern North Africans. Curr. Biol. 2017, 27, 3396–3402. [Google Scholar] [CrossRef] [PubMed]
  65. Schlebusch, C.M.; Malmstrom, H.; Gunther, T.; Sjodin, P.; Coutinho, A.; Edlund, H.; Munters, A.R.; Vicente, M.; Steyn, M.; Soodyall, H.; et al. Southern African ancient genomes estimate modern human divergence to 350,000 to 260,000 years ago. Science 2017, 358, 652–655. [Google Scholar] [CrossRef] [PubMed]
  66. Siska, V.; Jones, E.R.; Jeon, S.; Bhak, Y.; Kim, H.M.; Cho, Y.S.; Kim, H.; Lee, K.; Veselovskaya, E.; Balueva, T.; et al. Genome-wide data from two early Neolithic East Asian individuals dating to 7700 years ago. Sci. Adv. 2017, 3, e1601877. [Google Scholar] [CrossRef] [PubMed]
  67. Matisoo-Smith, E.A.; Gosling, A.L.; Boocock, J.; Kardailsky, O.; Kurumilian, Y.; Roudesli-Chebbi, S.; Badre, L.; Morel, J.P.; Sebai, L.L.; Zalloua, P.A. A European mitochondrial haplotype identified in ancient Phoenician remains from Carthage, North Africa. PLoS ONE 2016, 11, e0155046. [Google Scholar] [CrossRef] [PubMed]
  68. Stevanovitch, A.; Gilles, A.; Bouzaid, E.; Kefi, R.; Paris, F.; Gayraud, R.P.; Spadoni, J.L.; El-Chenawi, F.; Beraud-Colomb, E. Mitochondrial DNA sequence diversity in a sedentary population from Egypt. Ann. Hum. Genet. 2004, 68, 23–39. [Google Scholar] [CrossRef] [PubMed]
  69. Saunier, J.L.; Irwin, J.A.; Strouss, K.M.; Ragab, H.; Sturk, K.A.; Parsons, T.J. Mitochondrial control region sequences from an Egyptian population sample. Forensic Sci. Int. Genet. 2009, 3, e97–e103. [Google Scholar] [CrossRef] [PubMed]
  70. Kujanova, M.; Pereira, L.; Fernandes, V.; Pereira, J.B.; Cerny, V. Near eastern neolithic genetic input in a small oasis of the Egyptian Western Desert. Am. J. Phys. Anthropol. 2009, 140, 336–346. [Google Scholar] [CrossRef] [PubMed]
  71. Pagani, L.; Schiffels, S.; Gurdasani, D.; Danecek, P.; Scally, A.; Chen, Y.; Xue, Y.; Haber, M.; Ekong, R.; Oljira, T.; et al. Tracing the route of modern humans out of Africa by using 225 human genome sequences from Ethiopians and Egyptians. Am. J. Hum. Genet. 2015, 96, 986–991. [Google Scholar] [CrossRef] [PubMed]
  72. Elmadawy, M.A.; Nagai, A.; Gomaa, G.M.; Hegazy, H.M.; Shaaban, F.E.; Bunai, Y. Investigation of mtDNA control region sequences in an Egyptian population sample. Leg. Med. 2013, 15, 338–341. [Google Scholar] [CrossRef] [PubMed]
  73. EMPOP (https://empop.online/). Population EMP00351; EMPOP: Ismailia, Egypt, 2017. [Google Scholar]
  74. Schuenemann, V.J.; Peltzer, A.; Welte, B.; van Pelt, W.P.; Molak, M.; Wang, C.C.; Furtwangler, A.; Urban, C.; Reiter, E.; Nieselt, K.; et al. Ancient Egyptian mummy genomes suggest an increase of Sub-Saharan African ancestry in post-Roman periods. Nat. Commun. 2017, 8, 15694. [Google Scholar] [CrossRef] [PubMed]
  75. Doxey, D.M. Funerary beliefs and practices in the Middle Kingdom. In The Secrets of Tomb 10A; MFA Publications: New York, NY, USA, 2009; pp. 39–64. [Google Scholar]
  76. Ikram, S.; Dodson, A. The Mummy in Ancient Egypt: Equipping the Dead for Eternity; Thames & Hudson Ltd.: London, UK, 1998. [Google Scholar]
  77. Aufderheide, A.C. The Scientific Study of Mummies; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  78. Bard, K.A. An Introduction to the Archaeology of Ancient Egypt; Wiley-Blackwell: Hoboken, NJ, USA, 2015. [Google Scholar]
  79. De Meyer, M. The tomb of Henu at Deir el-Bersha. Egyptian Archaeol. 2007, 31, 20–24. [Google Scholar]
  80. Khairat, R. Next generation sequencing of DNA extracted from mummified tissue. Ph.D. Thesis, Universität Tübingen, Tübingen, Germany, 2013. [Google Scholar]
  81. Morris, A.G.; Heinze, A.; Chan, E.K.; Smith, A.B.; Hayes, V.M. First ancient mitochondrial human genome from a prepastoralist southern African. Genome Biol. Evol. 2014, 6, 2647–2653. [Google Scholar] [CrossRef] [PubMed]
  82. Gallego Llorente, M.; Jones, E.R.; Eriksson, A.; Siska, V.; Arthur, K.W.; Arthur, J.W.; Curtis, M.C.; Stock, J.T.; Coltorti, M.; Pieruccini, P.; et al. Ancient Ethiopian genome reveals extensive Eurasian admixture throughout the African continent. Science 2015, 350, 820–822. [Google Scholar] [CrossRef] [PubMed]
  83. Bramanti, B.; Thomas, M.G.; Haak, W.; Unterlaender, M.; Jores, P.; Tambets, K.; Antanaitis-Jacobs, I.; Haidle, M.N.; Jankauskas, R.; Kind, C.J.; et al. Genetic discontinuity between local hunter-gatherers and central Europe’s first farmers. Science 2009, 326, 137–140. [Google Scholar] [CrossRef] [PubMed]
  84. Posth, C.; Renaud, G.; Mittnik, A.; Drucker, D.G.; Rougier, H.; Cupillard, C.; Valentin, F.; Thevenet, C.; Furtwangler, A.; Wissing, C.; et al. Pleistocene mitochondrial genomes suggest a single major dispersal of non-Africans and a late glacial population turnover in Europe. Curr. Biol. 2016, 26, 827–833. [Google Scholar] [CrossRef] [PubMed]
  85. Bietak, M. Egypt and the Levant in the Egyptian World; Wilkinson, T., Ed.; Routledge: New York, NY, USA, 2007; pp. 417–448. [Google Scholar]
  86. Picardo, N.S. Middle Kingdom History, Politics and Social Organization. In The Secrets of Tomb 10A; MFA Publications: New York, NY, USA, 2009; pp. 17–38. [Google Scholar]
  87. Wilkinson, T. Early Dynastic Egypt; Routledge: New York, NY, USA, 2001. [Google Scholar]
  88. Skoglund, P.; Thompson, J.C.; Prendergast, M.E.; Mittnik, A.; Sirak, K.; Hajdinjak, M.; Salie, T.; Rohland, N.; Mallick, S.; Peltzer, A.; et al. Reconstructing Prehistoric African Population Structure. Cell 2017, 171, 59–71. [Google Scholar] [CrossRef] [PubMed]
  89. Drosou, K.; Price, C.; Brown, T.A. The kinship of two 12th Dynasty mummies revealed by ancient DNA sequencing. J. Arch. Sci. 2018, 17, 793–797. [Google Scholar] [CrossRef]
  90. Fregel, R.; Mendez, F.L.; Bokbot, Y.; Martin-Soca, D.; Camalich-Massieu, M.D.; Ávila-Arcos, M.C.; Underhill, P.A.; Shapiro, B.; Wojcik, G.; Rodríguez-Santos, F.J.; et al. Neolithization of North Africa involved 1 the migration of people from both the Levant and Europe. bioRxiv 2018. [Google Scholar] [CrossRef]
  91. Fazi, A.; Gobeski, B.; Foran, D. Development of two highly sensitive forensic sex determination assays based on human DYZ1 and Alu repetitive DNA elements. Electrophoresis 2014, 35, 3028–3035. [Google Scholar] [CrossRef] [PubMed]
  92. Loreille, O.; Koshinsky, H.; Fofanov, V.Y.; Irwin, J.A. Application of next generation sequencing technologies to the identification of highly degraded unknown soldiers’ remains. Forensic Sci. Int. Genet. Suppl. Ser. 2011, 3, e540–e541. [Google Scholar] [CrossRef]
  93. Templeton, J.E.; Brotherton, P.M.; Llamas, B.; Soubrier, J.; Haak, W.; Cooper, A.; Austin, J.J. DNA capture and next-generation sequencing can recover whole mitochondrial genomes from highly degraded samples for human identification. Invest. Genet. 2013, 4, 26. [Google Scholar] [CrossRef] [PubMed]
  94. Eduardoff, M.; Xavier, C.; Strobl, C.; Casas-Vargas, A.; Parson, W. Optimized mtDNA control region primer extension capture analysis for forensically relevant samples and highly compromised mtDNA of different age and origin. Genes 2017, 8, 237. [Google Scholar] [CrossRef] [PubMed]
  95. Marshall, C.; Sturk-Andreaggi, K.; Daniels-Higginbotham, J.; Oliver, R.S.; Barritt-Ross, S.; McMahon, T.P. Performance evaluation of a mitogenome capture and Illumina sequencing protocol using non-probative, case-type skeletal samples: Implications for the use of a positive control in a next-generation sequencing procedure. Forensic Sci. Int. Genet. 2017, 31, 198–206. [Google Scholar] [CrossRef] [PubMed]
  96. Parson, W.; Huber, G.; Moreno, L.; Madel, M.B.; Brandhagen, M.D.; Nagl, S.; Xavier, C.; Eduardoff, M.; Callaghan, T.C.; Irwin, J.A. Massively parallel sequencing of complete mitochondrial genomes from hair shaft samples. Forensic Sci. Int. Genet. 2015, 15, 8–15. [Google Scholar] [CrossRef] [PubMed]
  97. Pakstis, A.J.; Speed, W.C.; Fang, R.; Hyland, F.C.; Furtado, M.R.; Kidd, J.R.; Kidd, K.K. SNPs for a universal individual identification panel. Hum. Genet. 2010, 127, 315–324. [Google Scholar] [CrossRef] [PubMed]
  98. Sanchez, J.J.; Phillips, C.; Borsting, C.; Balogh, K.; Bogus, M.; Fondevila, M.; Harrison, C.D.; Musgrave-Brown, E.; Salas, A.; Syndercombe-Court, D.; et al. A multiplex assay with 52 single nucleotide polymorphisms for human identification. Electrophoresis 2006, 27, 1713–1724. [Google Scholar] [CrossRef] [PubMed]
  99. Kidd, K.K.; Pakstis, A.J.; Speed, W.C.; Grigorenko, E.L.; Kajuna, S.L.; Karoma, N.J.; Kungulilo, S.; Kim, J.J.; Lu, R.B.; Odunsi, A.; et al. Developing a SNP panel for forensic identification of individuals. For. Sci. Int. 2006, 164, 20–32. [Google Scholar] [CrossRef] [PubMed]
  100. Lazaridis, I.; Nadel, D.; Rollefson, G.; Merrett, D.C.; Rohland, N.; Mallick, S.; Fernandes, D.; Novak, M.; Gamarra, B.; Sirak, K.; et al. Genomic insights into the origin of farming in the ancient Near East. Nature 2016, 536, 419–424. [Google Scholar] [CrossRef] [PubMed]
  101. Fu, Q.; Posth, C.; Hajdinjak, M.; Petr, M.; Mallick, S.; Fernandes, D.; Furtwangler, A.; Haak, W.; Meyer, M.; Mittnik, A.; et al. The genetic history of Ice Age Europe. Nature 2016, 534, 200–205. [Google Scholar] [CrossRef] [PubMed]
  102. Kanzawa-Kiriyama, H.; Kryukov, K.; Jinam, T.A.; Hosomichi, K.; Saso, A.; Suwa, G.; Ueda, S.; Yoneda, M.; Tajima, A.; Shinoda, K.I.; et al. A partial nuclear genome of the Jomons who lived 3000 years ago in Fukushima, Japan. J. Hum. Genet. 2017, 62, 213–221. [Google Scholar] [CrossRef] [PubMed]
  103. Gargis, A.S.; Kalman, L.; Lubin, I.M. Assuring the Quality of Next-Generation Sequencing in Clinical Microbiology and Public Health Laboratories. J. Clinic. Microbiol. 2016, 54, 2857–2865. [Google Scholar] [CrossRef] [PubMed]
  104. Richards, S.; Aziz, N.; Bale, S.; Bick, D.; Das, S.; Gastier-Foster, J.; Grody, W.W.; Hegde, M.; Lyon, E.; Spector, E.; et al. Standards and guidelines for the interpretation of sequence variants: A joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet. Med. 2015, 17, 405–424. [Google Scholar] [CrossRef] [PubMed]
Genes 09 00135 g001 550
Figure 1. Head of the Djehutynakht mummy (2010–1961 BC). Photograph © 2018 Museum of Fine Arts, Boston, USA.
Figure 1. Head of the Djehutynakht mummy (2010–1961 BC). Photograph © 2018 Museum of Fine Arts, Boston, USA.
Genes 09 00135 g001
Genes 09 00135 g002 550
Figure 2. Sequencing strategy in this study. Molecular work performed at the Federal Bureau of Investigation (FBI) Laboratory is in blue, while work performed at Harvard Medical School (HMS) is in pink. Cap1-Lib2 has both colors because the extraction and library preparation of Lib2 was done at the HMS, while the hybridization capture took place at the FBI. CAP-LIB1 is the result of merging data from runs I and II, while CAP-LIB2 is the result of merging runs III and IV. Lib1 and Lib2 were shotgun-sequenced together on the HiSeq at National Bioforensic Analysis Center (NBFAC, run V).
Figure 2. Sequencing strategy in this study. Molecular work performed at the Federal Bureau of Investigation (FBI) Laboratory is in blue, while work performed at Harvard Medical School (HMS) is in pink. Cap1-Lib2 has both colors because the extraction and library preparation of Lib2 was done at the HMS, while the hybridization capture took place at the FBI. CAP-LIB1 is the result of merging data from runs I and II, while CAP-LIB2 is the result of merging runs III and IV. Lib1 and Lib2 were shotgun-sequenced together on the HiSeq at National Bioforensic Analysis Center (NBFAC, run V).
Genes 09 00135 g002
Genes 09 00135 g003 550
Figure 3. (a) Size distribution of the reads that aligned to the mtGenome in CAP-LIB1. Reads ≤35 bp or >70 bp (shown in orange) were removed from the final data to avoid any impact from nonspecific [34] or contaminating reads. The length with the greatest number of reads was 47 bp. (b) Size distribution of the reads that aligned to the mtGenome in CAP-LIB2. Only reads >35bp and ≤70 bp (shown in blue) were retained in the final data. The length with the greatest number of reads was 38 bp.
Figure 3. (a) Size distribution of the reads that aligned to the mtGenome in CAP-LIB1. Reads ≤35 bp or >70 bp (shown in orange) were removed from the final data to avoid any impact from nonspecific [34] or contaminating reads. The length with the greatest number of reads was 47 bp. (b) Size distribution of the reads that aligned to the mtGenome in CAP-LIB2. Only reads >35bp and ≤70 bp (shown in blue) were retained in the final data. The length with the greatest number of reads was 38 bp.
Genes 09 00135 g003
Genes 09 00135 g004 550
Figure 4. Distribution of 40,910 reads over the entire mtGenome. Sequence coverage at each position ranged from 5× to 247× (average 108×).
Figure 4. Distribution of 40,910 reads over the entire mtGenome. Sequence coverage at each position ranged from 5× to 247× (average 108×).
Genes 09 00135 g004
Table
Table 1. Sequencing statistics for the extraction reagent control (RB) and high throughput sequencing (HTS) reagent negative controls (NC).
Table 1. Sequencing statistics for the extraction reagent control (RB) and high throughput sequencing (HTS) reagent negative controls (NC).
Total # Reads# Unique Reads Mapped to hg19# Unique Reads Mapped to the mtGenome
Cap-Lib1RB428,1925600
NC47,4141210
Cap2-Lib2RB744,22857010
NC250,2821287
Table
Table 2. Shotgun sequencing mapping statistics.
Table 2. Shotgun sequencing mapping statistics.
Sequence StatisticsFBI Shotgun-Lib1HMS Shotgun-Lib2
aNumber of raw paired reads164,451,485266,162,607
bNumber of reads mapped to the human genome hg19 and rCRS3,692,50419,485,309
cPercentage of endogenous human DNA 2.24%6.57%
dNumber of unique human reads with Q >30 1,595,2397,691,326
eAverage coverage hg19
Average coverage mtGenome		0.02×
4.21×		0.09×
8.93×
fNumber of unique mapped human reads with signs of damage
(PMDtools score >3)		518,381344,995
Table
Table 3. RX and RY determination using unique mapped reads with Q > 30.
Table 3. RX and RY determination using unique mapped reads with Q > 30.
SamplesLib1SexLib2Sex
Mapped reads1,595,239 7,691,326
Mapped to X37,605 176,181
Mapped to Y3732 16,469
RY0.090 0.0855
95% CI0.087–0.093♂0.084–0.087♂
RX0.45 0.433
95% CI0.429–0.471♂0.4–0.466♂
Bold symbolizes the CI used to determine sex.
Table
Table 4. RX and RY determination using unique mapped reads that showed signs of DNA damage.
Table 4. RX and RY determination using unique mapped reads that showed signs of DNA damage.
SamplesLib1SexLib2Sex
Mapped reads518,381 344,995
Mapped to X11,512 6688
Mapped to Y1133 637
RY0.0896 0.087
95% CI0.0846–0.0946♂0.0805–0.0934♂
RX0.4216 0.3638
95% CI0.3987–0.4447♂0.3261–0.4015♂
Bold symbolizes the CI used to determine sex.

Share and Cite

MDPI and ACS Style

Loreille, O.; Ratnayake, S.; Bazinet, A.L.; Stockwell, T.B.; Sommer, D.D.; Rohland, N.; Mallick, S.; Johnson, P.L.F.; Skoglund, P.; Onorato, A.J.; et al. 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. Genes 2018, 9, 135. https://doi.org/10.3390/genes9030135

AMA Style

Loreille O, Ratnayake S, Bazinet AL, Stockwell TB, Sommer DD, Rohland N, Mallick S, Johnson PLF, Skoglund P, Onorato AJ, et al. 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. Genes. 2018; 9(3):135. https://doi.org/10.3390/genes9030135

Chicago/Turabian Style

Loreille, Odile, Shashikala Ratnayake, Adam L. Bazinet, Timothy B. Stockwell, Daniel D. Sommer, Nadin Rohland, Swapan Mallick, Philip L.F. Johnson, Pontus Skoglund, Anthony J. Onorato, and et al. 2018. "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" Genes 9, no. 3: 135. https://doi.org/10.3390/genes9030135

APA Style

Loreille, O., Ratnayake, S., Bazinet, A. L., Stockwell, T. B., Sommer, D. D., Rohland, N., Mallick, S., Johnson, P. L. F., Skoglund, P., Onorato, A. J., Bergman, N. H., Reich, D., & Irwin, J. A. (2018). 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. Genes, 9(3), 135. https://doi.org/10.3390/genes9030135

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop