Abstract In November and December 2013, unidentified human skeletal remains buried in a mokgwakmyo (a traditional wooden coffin) were unearthed while conducting an archaeological investigation near Gyeongju, which was the capital of the Silla Kingdom (57 BCE– 660 CE) of ancient Korea. The human skeletal remains were preserved in relatively intact condition. In an attempt to obtain biological information on the skeleton, physical anthropological, mitochondrial DNA, stable isotope and craniofacial analyses were carried out. The results indicated that the individual was a female from the Silla period, of 155 ± 5 cm height, who died in her late thirties. The maternal lineage belonged to the haplogroup F1b1a, typical for East Asia, and the diet had been more C 3 - (wheat, rice and potatoes) than C 4 -based (maize, millet and other tropical grains). Finally, the face of the individual was reconstructed utilizing the skull (restored from osseous fragments) and three-dimensional computerized modeling system. This study, applying multi-dimensional approaches within an overall bio-anthropological analysis, was the first attempt to collect holistic biological information on human skeletal remains dating to the Silla Kingdom period of ancient Korea.

Citation: Lee W-J, Woo EJ, Oh CS, Yoo JA, Kim Y-S, Hong JH, et al. (2016) Bio-Anthropological Studies on Human Skeletons from the 6th Century Tomb of Ancient Silla Kingdom in South Korea. PLoS ONE 11(6): e0156632. https://doi.org/10.1371/journal.pone.0156632 Editor: Mark Spigelman, Hebrew University, ISRAEL Received: November 6, 2015; Accepted: May 17, 2016; Published: June 1, 2016 Copyright: © 2016 Lee et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability: All relevant data are within the paper and its Supporting Information files. Funding: This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (NRF-2015R1C1A1A02036740). WJL received the funding. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing interests: The authors have declared that no competing interests exist.

Introduction The Silla Kingdom, beginning its history as a small city-state, prospered during the first millennium CE (57 BCE– 935 CE), becoming a well-organized state that ruled most of the Korean peninsula. Understanding the history of the Silla Kingdom from an archaeological perspective is very important, as it remains highly relevant. In fact, much of present-day Korea’s cultural heritage is rooted in its traditions. In Silla funeral culture, various types of tombs corresponding to the different social classes were constructed. Among them, wood coffins (mokgwanmyo) and stone-mounted wooden chamber tombs (jeokseok mokgwangmyo) are representative. In spite of the importance of Silla tombs from the academic perspective, very few bio-anthropological studies have been carried out, due to the typically degenerated condition of the bones contained therein. This is largely the result of the acidic soil and the alternate hot/wet, cold/dry weather patterns, which tend to accelerate the decomposition of a human cadaver [1]. Archaeologists/anthropologists in South Korea therefore have not had enough opportunities to explore intact ancient human bones from the Silla period, except for a few anthropometrical studies [2], cremation and scanning electron microscope (SEM) studies [3], and stable isotope analysis for dietary pattern studies [4]. These previous studies however are limited, as, having applied only a single method rather than utilizing a holistic approach with a variety of analysis tools, they can provide only partial biological information. The present investigation represents a rare chance to examine intact human remains dating to the Silla period. Researchers from varied fields collaborated to acquire biological data based on physical anthropological, mitochondrial DNA (mtDNA), stable isotope and craniofacial analysis. This study can provide significant information on the ancient Silla people, which is very meaningful not only to academics but also to the general Korean public.

Materials and Methods Archaeological Information In November and December 2013, rescue archaeological excavation was carried out at 94 Gyo-dong, Gyeongju-si, Gyeongsangbuk-do province of South Korea before the construction of new driveway by the members of the Silla Cultural Heritage Research Institute (Fig 1A and 1B). The archaeological investigation was conducted under the permission of the Cultural Heritage Administration of Korea (Approval number: 2013–0285). At the district A in the archaeological investigation site, a complete wooden coffin (mokgwanmyo) was discovered in a tomb, and human skeletons were found within the coffin (Fig 1C). The size of the coffin was 230 × 80 × 20 cm in length, width and depth respectively. The significance of this tomb was that the human remains were preserved in almost complete form. The surroundings included a marshy area which might contribute to the isolation of the dead body from external destructive environments for a long time. The tomb was also occupied by diverse relics besides the human remains. The analyses of the relics indicated that the tomb was constructed in the Silla period, which also suggest that the human skeletons were from the same period. PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Fig 1. Location of the excavation site and the uncovered coffin. A: the red dot and arrow indicate the location of archaeological excavation site in Gyeongju-si in South Korea; B: the red dot and arrow indicate the location of the wooden coffin (mokgwanmyo) at the excavation site; C: the human skeletons in the coffin were preserved in a relatively intact condition. https://doi.org/10.1371/journal.pone.0156632.g001 Prior to the excavation of the human remains, a public announcement was made based on the Funeral Services Related Act of Korea (Section 2) in order to find possible descendants of the decedent. As the skeletal remains had not been claimed 60 days after the announcement by any possible descendants, the excavation was allowed to proceed with the permission of the Cultural Heritage Administration of Korea (Approval number: 2013–0285). The skeletal remains were moved to the Bio-anthropology Laboratory of Seoul National University Medical College for detailed examination. The specimen and repository numbers were given as GJGD-1 and #321. The repository was involved in the Joseon Dynasty Human Sample Collection of the Department of Anatomy at the Seoul National University College of Medicine (103 Daehak-ro, Chongno-gu, Seoul, South Korea, 03080). The human sample has not been publicly deposited but been allowed to access by other relevant researchers. This study was performed by the permission of the Institutional Review Board of Seoul National University (IRB number: H-0909-049-295). Since the osseous fragments were commingled with mud in the archaeological site (Fig 2), all human remains were recovered with the surrounding soil for the transfer to the laboratory. To minimize any contamination, the skeletons in the archaeological site were not allowed to contact with anyone until the removal and transfer from the site to the laboratory had been completed by the excavator (CSO). Every researcher who took part in this study was identified in order to verify cross-contamination during the DNA analysis. PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Fig 2. Fragmentary skull portion. The upper body part of the human skeletal remains was magnified from Fig 1C, showing the comingling with other archaeological relics and mud. https://doi.org/10.1371/journal.pone.0156632.g002 Physical Anthropological Analysis The analysis procedure was conducted in accordance with the Vermillion Accord on Human Remains adopted by World Archaeological Congress [5]. Skeletal preparation (cleaning) followed the suggestions described in the text written by White and Folkens [6]. After pre-treatment of the bones, physical anthropological analysis of the specimen was conducted. Stature was determined by measuring the maximum length of the femur following Fujii’s [7] and Pearson’s [8] methods. Sex was determined by morphological analysis using the pelvis and the restored skull [6]. Primary indicators for the sex determination were the greater sciatic notch and pre-auricular sulcus on the hip bone. The auxiliary indicators were morphological features of the skull. Age was estimated by the degenerative changes of the auricular and pubic symphyseal surface on the hip bone [6]. As all the teeth were remained in good condition except for the mandibular left central and lateral incisors, right canine and both 2nd molars, the age at death was estimated by the examination of the length ratios and degenerative change of the pulp cavity of the remaining teeth [9]. Restoration of the Fragmentary Skull Before craniofacial analysis and reconstruction, each osseous piece of the skull was cleaned utilizing lukewarm water with brushes and wooden probes. The washed osseous pieces were dried in room temperature for 48 hours, and prepared for the skull restoration (Fig 3). The restoration procedure was performed in two steps: at first, the virtual restoration with computerized three-dimensional (3D) modelling system (Geomagic® FreeForm® from 3D Systems®, USA) to simulate how the fragmentary osseous pieces could be assembled accurately; at second, the real restoration of the skull with actual fragmentary bones referring to the completed virtual restoration. PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Fig 3. Osseous fragments of the skull after pre-treatment. https://doi.org/10.1371/journal.pone.0156632.g003 Each skull fragment were scanned using a cone beam CT (Dentri™ from HDXILL®, Seoul, Korea) with field of view (FOV) of 16 x 8 cm, a voxel size of 0.20mm and +150 to +200 of Housefield units (Hus). The scanned 3D images of the fragments formatted as Digital Imaging and Communications in Medicine (DICOM) data were converted into stereolithography (STL) formats, and imported into the computerized 3D modelling system for the virtual restoration. Later, the actual fragments of the skull were assembled manually referring to the virtually restored 3D skull model. Metric and Paleopathological Analyses Metric information for the skeleton was measured by digital vernier & spreading caliper (Mitutoyo®; Tokyo, Japan) and osteometric board. Available measurements were examined as far as the relevant skeletal portions were observable. In addition, the restored skull and the postcranial skeleton were examined macroscopically to investigate any non-metric traits and pathological changes. Craniofacial Reconstruction The facial morphology and appearance of the individual was depicted by the technique of craniofacial reconstruction (CFR). The CFR using the restored skull was conducted by employing the computerized 3D modeling system. The sex, age and ancestry for the CFR followed the outcomes of physical anthropological analysis. An average facial tissue depths dataset for Korean adults was applied to locate tissue depth gauges at corresponding anatomical landmarks onto the surface of the skull [10]. Each major facial muscle was rebuilt as accurately as possible following anatomical structure of the skull. There are a number of guidelines to predict facial components such as eyes, nose, mouth and ears. The description in detail was stated in the previous studies into Korean CFRs [11,12]. For the final stage, a skin layer was added over the muscle and skull structure referring to the tissue depth gauges by utilizing haptic tools in the computerized modelling system. Stable Isotope Analysis Collagen is extracted from archaeologically obtained bone. Sample was prepared for isotope analysis, following a modified Longin method, with the addition of an ultrafiltration step [13]. Briefly, 200mg of powered bone were demineralized in 0.5M HCl solution at 4°C for 16 hour; and then in a pH 3 water, collagen was gelatinized at 75°C for 48 hour. After removing insoluble residues with a 5–8 Ezee® mesh filter, the remaining solution was filtered to remove impurities using 30kDa filters once again [14]. Isotope content and ratio of C and N were determined using a continuous-flow stable isotope ratio mass spectrometer (IsoPrime-EA, Micromass, UK) linked with a CN analyzer (NA Series 2, CE Instruments, Italy). Carbon and nitrogen isotope compositions (δ13C and δ15N) were calculated as: δ(‰) = [(R sample /R standard ) - 1] × 1000 where R is the ratio of 13C/12C or 15N/14N; and the standards were the Pee Dee Belemnite (PDB) for carbon and atmospheric (AIR) for nitrogen. Two replicate analyses were carried out for the sample and the average was used for the statistical analysis. Multiple replicate analyses indicated that standard deviations for the δ13C and δ15N measurements were <0.1‰ and <0.2‰, respectively. Mitochondrial DNA Analysis The femur fragment (sample No. 321) was collected with a sterilized knife. This were then exposed to UV irradiation for 20 min, and subsequently immersed in 5.4% (w/v) sodium hypochlorite. When the samples were washed with distilled water and absolute ethanol, they were then air-dried and pulverized to a fine powder using a SPEX 6750 Freezer / Mill (SPEX SamplePrep, Metuchen, NJ) [15]. Bone powder (0.5 g) was incubated in 1 mL of lysis buffer (EDTA 50 mM, pH 8.0; 1 mg /mL of proteinase K; SDS 1%; 0.1 M DTT) at 56°C for 24 h. Total DNA was extracted with an equal volume of phenol / chloroform / isoamyl alcohol (25:24:1), and then was treated with chloroform / isoamyl alcohol (24:1). DNA isolation and purification was performed using a QIAmp PCR purification kit (Qiagen, Hilden, Germany). The purified DNA was eluted in 40 μl of EB buffer (Qiagen) [15]. After quantification was performed by NanoDrop™ ND-1000 Spectrophotometer (Thermo Fisher Scientific, MA, USA), 40 ng of the extracted DNA was mixed with the reagent premix containing 1X AmpliTaq Gold® 360 Master Mix (Life Technologies, USA) and 10 pmol of each primer (Integrated DNA Technology, USA). PCR conditions used in this study were as follows: pre-denaturation at 94°C for 10 min; 45 cycles of pre-denaturation at 94°C for 10 min; 50 cycles of annealing at 45–56°C for 30 sec; extension at 72°C for 30 sec; final extension at 72°C for 10 min. PCR amplification was performed using a PTC-200 DNA Engine (Bio-Rad Laboratories, Hercules, CA). Primer sets used in this study are summarized in supporting information 1 (S1 Table). The PCR products were separated on 2.5% agarose gel and then isolated using a Qiagen gel extraction kit (Qiagen, Germany). The sequencing of each amplicon was performed by ABI Prism® 3100 Genetic Analyzer (Applied Biosystems, USA), using ABI Prism® BigDye™ Terminator Cycle Sequencing Ready Reaction Kit (Applied Biosystems, USA). Multiple replications of DNA extraction, PCR amplication, cloning and sequencing were performed with an independent laboratory. The obtained DNA sequences were compared with the revised Cambridge Reference Sequence (rCRS; accession number: NC_012920), to identify the sequence differences. The resultant control region mutation motifs were imported into the online-based program MitoTool [16] for mtDNA haplogrouping. In order to minimize any modern DNA contamination of ancient samples, the mtDNA profiles of the researcher (CSO) was determined with the permission of the Institutional Review Board of Seoul National University (IRB number: H-0909-049-295). The researcher was the only person who sampled and carried out mtDNA analysis in this study. The researcher’s mtDNA haplotype was compared with the mtDNA profile from the Silla bone to rule out the possibility of modern DNA contamination. During the mtDNA work, the researcher always wore protection gloves, masks, gowns and head caps. Our mtDNA lab facilities were set up in accordance with the protocol of Hofreiter et al. [17]. The rooms for mtDNA extraction or PCR preparation were physically separated from our main PCR lab. The DNA extraction/PCR preparation rooms were equipped with night UV irradiation, isolated ventilation, and a laminated flow hood.

Discussion Not a few discoveries have been made in relation to the tombs originally constructed for the royal families and nobles during the era of the ancient Korean kingdoms. In those archaeological cases, the human remains were also exhumed from the tombs for the purpose of bio-anthropological studies [14,23]. Although the previous pioneering studies have revealed invaluable information on the biological aspects of the ancient Korean people, our knowledge of the societies of the first millennium is still very incomplete. In this regard, the current study, having employed various investigative tools, provides vivid biological information on one individual from the Silla period. The restoration of the fragmentary skull was the first important step in this study to acquire physical anthropological information. A number of virtual techniques have been applied for examination of poorly-preserved human specimens. Representatives are reverse engineering, computer-aided design (CAD), and rapid prototyping (RP) techniques [30,31]. However, those tools are rather unnecessarily sophisticated and time-consuming instead, the present study employed a method of virtual skull restoration utilizing a computerized 3D modeling system. By this method, the fragmentary human remains could be assembled very accurately and rapidly. By this technique moreover, missing portions could be repaired efficiently, particularly where an opposite portion of the structure is intact. Of course, the application of what is known as the “mirroring technique” must be carried out with caution where significant asymmetry is observed in the remaining structure. In the current case, a more reliable craniometric analysis was possible due to the virtually restored skull and the obtained primary craniometric measurements. The cephalic index (CI), at 73.5, suggested a dolichocephalic trait differed from the average or modern Korean adults, who typically demonstrate the brachycephalic or mesochephalic type [19,20]. Reportedly, cranial indices among Koreans changed, from mesochephalic (CI = 77.6 ± 4.9) during the ancient period (4th-7th century) to brachychephalic (CI = 82.3 ± 5.4) during the Joseon period (1392–1910 CE) [22,32]. When considering these and the other relevant previous studies [19,20,22,32], the long and narrow cranium and narrower facial shape of the individual in this study should be regarded as an idiosyncratic characteristic of the skull and not a typical feature of East Asians in general or of Koreans in particular. Two hypotheses are possible. The first is that the ancient Koreans living in the region during the Silla period had a narrower cranium than did people living contemporaneously in other regions, the second is that the head of this particular had been artificially narrowed in a cultural rite practiced at that time. The latter however was discounted because, other than the inclined forehead, the cranium did not manifest any typical changes caused by artificial cranial deformation. The former hypothesis seems the more persuasive. One individual however cannot account for any possible variability of the Silla population relative to others of the same era. Thus, further investigation with significantly more numerous skull samples from the period will be necessary. CFR is a technique used to rebuild a living facial appearance onto a skull in order to recognize or identify an individual. The accuracy of CFR has been of primary importance to maintaining the reliability of its applications to forensic investigation and archaeological research. The most recent study employing quantitative geometric surface comparison found that 75% of CFR surfaces showed less than 2.0 mm of error relative to the corresponding actual faces [11,12,33]. In stable isotope analysis, the carbon isotope ratio δ13C shows the proportions of C 3 - (wheat, rice and potatoes) and C 4 -based (maize, millet and other tropical grains) foods in the diets of historical individuals. C 4 -based foods show a higher δ13C value than do C 3 -based foods [34]. A previous study determined that the δ13C value of C 3 foods is -25.4 ± 1.6‰ while that of C 4 grains is -10.0 ± 0.5‰ [35]. The δ13C value for the Silla individual examined in the present study was -19.72 ± 0.12‰, indicating that she had consumed a more C 3 -based diet. In fact, this δ13C data is similar to the results of studies on Joseon skeletons [25–27], and slightly higher than those for ancient skeletons from Yeanri [14], Imdang [23], and Neukdo [24]. On the other hand, the present δ 15N value (7.89 ± 0.38‰) is markedly lower than the data on any of the ancient or medieval skeletal remains discovered in South Korea. Judging from a previous study on human bones from a Neolithic settlement in China and modern human hair from England and Germany, it is likely that the present δ 15N value represents a vegetarian diet [36]. Considering the prevalence of Buddhism among the Silla people of the first millennium, the individual examined in the current study probably practiced strict vegetarianism. Regarding mtDNA analysis, relevant previous studies already have reported ancient DNA findings on Korean skeletons. Lee et al. [37] presented a genetic characterization of 11 ancient Korean skeletal remains dating to times ranging from the Paleolithic era to the Goryeo Dynasty (918–1392 CE): G3a (Paleolithic), B4f1, D4c, B4b1a1 (Neolithic), D4, D4c1b (Bronze Age), D4e1a (Baekje, 18 BCE-660 CE), and F1a1, D6, A5c, N9a1 (Goryeo). In the current study, however, the Silla individual’s maternal lineage was that of the haplogroup F1b1a. Whereas this is typical of East Asian populations, we must note that it is not the dominant group among modern Koreans [38]. Actually, according to Hong et al. [38], the most common mtDNA sub-haplogroups found in Northeast Asian populations are D4a and D4, followed by A5a and Y1. We admit that our mtDNA results will be truly significant only once a sufficient number of additional cases have been examined in future research. The forthcoming ancient DNA studies on archaeological remains will broaden the information on the genetic traits of ancient Koreans, thereby expanding the knowledge base concerning the migrations of East Asian peoples as well.

Conclusion This study, having applied multi-dimensional approaches within an overall bio-anthropological analysis, demonstrates, for the first time, holistic features in the form of biological information on an individual from the ancient Silla period in Korea. The human skeleton examined was estimated to have been that of a female of 155 ± 5 cm height who died in her late 30s. The fragmentary skull was successfully reassembled by reconstructing the damaged or missing portions utilizing computerized 3D modeling software According to a cephalometric analysis, the cephalic index (CI) of the individual was 73.5, which suggests the dolichocephalic trait (long-headedness). This, as compared with similar demographic groups of the era, is unusual. Further, the face of the individual was obtained by craniofacial reconstruction (CFR) technique employing the reconstructed skull and the above-noted computerized 3D modeling. A stable isotope analysis demonstrated that the Silla individual had consumed a more C 3 (rather than C 4 )-based diet, which, in this particular case, probably was vegetarian. An additional mtDNA analysis showed that the Silla individual’s maternal lineage was that of the haplogroup F1b1a typical of East Asia.

Acknowledgments The authors would like to thank Dr. Ji-Woong Chung (Dacapo Oral & Maxillofacial Clinic in Seoul for the CT scan of the skull), Dr. Seong-Nam Kim (Chicago Children’s Dental Clinic in Seoul for the radiographs and cephalometic analysis of the skull).

Author Contributions Conceived and designed the experiments: WJL EJW CSO JAY YSK JHH DHS. Performed the experiments: WJL EJW CSO JAY YSK JHH AYY SDL DHS. Analyzed the data: WJL EJW CSO JAY YSK JHH AYY DHS. Contributed reagents/materials/analysis tools: WJL EJW CSO JAY YSK JHH AYY CMW JOJ SJC SDL DHS. Wrote the paper: WJL EJW CSO JAY YSK JHH CMW DHS.