pmc logo imageJournal ListSearchpmc logo image
Logo of ajhgJournal URL: redirect3.cgi?&&auth=0Vm6-HkPH92GquwsVeLyAAv2bt4MBeOH5-pBQx5bd&reftype=publisher&artid=1180365&article-id=1180365&iid=120389&issue-id=120389&jid=203&journal-id=203&FROM=Article|Banner&TO=Publisher|Other|N%2FA&rendering-type=normal&&http://www.journals.uchicago.edu/AJHG/index.html
Am J Hum Genet. 2003 August; 73(2): 247–260.
Published online 2003 July 10.
PMCID: PMC1180365
Nuclear and Mitochondrial DNA Analysis of a 2,000-Year-Old Necropolis in the Egyin Gol Valley of Mongolia
Christine Keyser-Tracqui,1 Eric Crubézy,2 and Bertrand Ludes1,2
1Laboratoire d’Anthropologie Moléculaire, Institut de Médecine Légale, Strasbourg, France, and 2Anthropobiologie, Université Paul Sabatier, CNRS, UMR 8555, Toulouse, France
Address for correspondence and reprints: Dr. Christine Keyser-Tracqui, Laboratoire d’Anthropologie Moléculaire, Institut de Médecine Légale, 11, rue Humann, 67085 Strasbourg Cedex, France. E-mail: ckeyser/at/mageos.com
Received February 26, 2003; Accepted May 7, 2003.
Abstract
DNA was extracted from the skeletal remains of 62 specimens excavated from the Egyin Gol necropolis, in northern Mongolia. This burial site is linked to the Xiongnu period and was used from the 3rd century b.c. to the 2nd century a.d. Three types of genetic markers were used to determine the genetic relationships between individuals buried in the Egyin Gol necropolis. Results from analyses of autosomal and Y chromosome short tandem repeats, as well as mitochondrial DNA, showed close relationships between several specimens and provided additional background information on the social organization within the necropolis as well as the funeral practices of the Xiongnu people. To the best of our knowledge, this is the first study using biparental, paternal, and maternal genetic systems to reconstruct partial genealogies in a protohistoric necropolis.
Introduction

In recent years, molecular studies have become widely employed to investigate parentage relationships within burial groups (Fily et al. 1998; Stone and Stoneking 1999; Schultes et al. 2000; Clisson et al. 2002), because morphological indicators of kinship are much less precise than the genetic data potentially available by analysis of ancient DNA. Understanding genetic relationships within and between burial sites helps us to understand the organization of sepulchral places and the origin of human remains recovered (e.g., unrelated individuals or members of a single or a limited number of family groups). This should be the first step of any work devoted to the history of settlement based on the investigation of remains from a cemetery, because every external inclusion in a group of subjects sharing a common parentage may introduce a bias (Crubézy et al. 2000).

In the present study, we examined biological kinship in a necropolis from the Xiongnu period, a culture known mainly through the graves discovered in 1943 by a joint Mongolian-Russian expedition in the Noin-Ula Mountains in northern Mongolia (Rudenko 1970) but also through other funerary sites of the Selenge Basin (Konovalov 1976). The Xiongnu were an ancient nomadic Turkomongolian tribe who were first described in Chinese manuscripts as early as the 4th century b.c. (Minajev 1996). In the 3rd century b.c., Xiongnu tribes rose to great power and created the first empire governed by Central Asian nomads. They ruled over a territory that extended from Lake Baikal in the north to the Gobi desert in the south and from western Manchuria in the east to the Pamirs in the west. During the newly established Han dynasty (206 b.c. to a.d. 220), China expanded its borders, and the Xiongnu empire lost ground (Marx 2000).

According to radiocarbon dating, the Egyin Gol site was used from the 3rd century b.c. to the 2nd century a.d. (i.e., over the whole Xiongnu period). It is located in northern Mongolia, in a cold environment favorable to a good preservation of the DNA (Burger et al. 1999; Leonard et al. 2000). We studied genetic diversity at the Egyin Gol site, first by use of autosomal STRs. Autosomal STRs consist of tandemly organized repeats of short nucleotide patterns (2–6 bp), which are transmitted according to a Mendelian mode of inheritance. These genetic markers took precedence in our study, owing to their excellent power of discrimination for the study of close parentage relationships. They also represent propitious markers for ancient DNA studies because of their small size and because they allow detection of sample contamination (Hummel et al. 2000). Moreover, they can be simultaneously amplified, reducing to an absolute minimum the amount of sample material necessary for kinship analysis. Although both maternal and paternal genetic contributions can be assessed with autosomal markers, such as STRs, we also studied the genetic diversity by typing the nonrecombining part of the Y chromosome, as well as the hypervariable region I (HVI) of the mtDNA. We studied paternal and maternal transmitted polymorphisms to complete the data obtained by autosomal STR analyses and, above all, to confirm the authenticity of the molecular data obtained from the ancient Egyin Gol specimens. These polymorphisms also provided additional information on the genetic history of the Xiongnu tribes.

Material and Methods

Site
The necropolis is located in the Egyin Gol valley near the Egyin Gol river, ~10 km from its confluence with the Selenge, a main tributary of Lake Baikal (fig. 1. The valley's position is −49° 27′ N, 103° 30′ E, and it has a continental climate, with an average annual temperature of −1°C. The winter (October to April) is cold (with temperatures often dropping to −30°C in January and February), whereas Summer (July to September) is pleasant (with temperatures sometimes as high as 22°C). Precipitation is light (300–400 mm per year). Because of its relatively high altitude (885 m), the valley floor is covered with snow from mid-November to April, and ice thickness on the Selenge reaches 1.8 m during this period. Permafrost was found in some areas by the geologists who were present on the site.
Figure  1Figure 1
Location of the Egyin Gol site

From 1997 to 1999, the burial site was wholly excavated by a French-Mongolian expedition, under the sponsorship of UNESCO, after preliminary boring revealed the excellent preservation of the graves (Crubézy et al. 1996). The necropolis comprised a total of 103 graves, among which 84 were excavated by the archaeological mission. The 19 remaining graves had been explored before the arrival of the mission in Mongolia, and no data were available on these spots. Graves were organized on both sides of a small depression on the river valley, in four sectors that were designated “A,” “B,” “C,” and “D” (figs. 2, 3, and 4). The southern sector (A) was composed of four double graves (32/32A, 33/33A, 37/37A, and 38/38A), each of which contained two sets of remains that were probably buried within the same period (Murail et al. 2000). Grave 27 was surmounted by a standing stone and was found to conceal exceptional furniture. In eight graves (18, 47, 49, 54, 59, 69, 83, and 85), secondary deposits (bones of very young children) were found beside the deceased.

Figure  2Figure 2
Map of the necropolis showing the autosomal STR data. Graves are represented by circles. Letters A, B, C, and D refer to the four sectors distinguished in the present study. Dotted lines define the boundary of these sectors.
Figure  3Figure 3
Map of the necropolis showing the Y chromosome STR data. Graves containing specimens of the same patrilineage are represented by an identical geometric figure.
Figure  4Figure 4
Map of the necropolis showing the mtDNA sequences data. Graves containing specimens of the same matrilineage are represented by an identical geometric figure.

The associated funeral material was of great interest and allowed us to link the necropolis to the Xiongnu culture (Crubézy et al. 1996). Bone samples from 31 specimens scattered across the necropolis were dated by carbon 14 (14C) determinations. The projection of the 31 mean values corresponding to each radiocarbon datum were linearly extrapolated, by use of UNIRAS software, to establish clines, which are represented by shades of gray on the necropolis map (fig. 5). This diagram suggests a topographical development of the burial ground, with a progressive expansion from south to north. Indeed, grave 28, slightly remote in the southern sector, was found to be the oldest of the necropolis, followed by grave 27 and the double burials. Therefore, sector A is probably the oldest, even though some graves located around it appear more recent. Sectors B and C seem more recent, although some graves situated near the center of sector B might have been implanted earlier.

Figure  5Figure 5
Radiocarbon dating map

The graves were at a depth of 2–5 m and were delimited by stones set in circles with diameters of several meters. They were protected by several layers of stones included in a loessial sediment. Chests and coffins were perfectly visible and relatively well preserved, as were most of the artifacts made of perishable matter (e.g., horn and bone) that were found in the graves.

Samples
Excavation of the 84 unexplored graves resulted in the recovery of 99 human skeletons (including double graves and secondary deposits). In most instances, complete and articulated skeletons were recovered, but, in some cases (e.g., secondary deposits or looted graves), numerous bones were missing or severely damaged. Most of the skeletal material was in an excellent state of preservation, as was confirmed by the mineral/organic composition of the bones, which did not differ significantly from that of contemporary bones. For instance, the mean ± SD crystallinity index was 0.07 ± 0.02, close to that of ice-preserved ancient bones (Person et al. 1996). Mean ± SD percentages of carbon and nitrogen were 13.8% ± 0.8 and 4.2% ± 0.2, respectively, quite similar to those of reference material obtained from surgical samples.

Sex was established according to the methodology developed by Murail et al. (1999). Age at death was estimated using dental calcification for the children and epiphyseal fusion for the adolescents (Crubézy et al. 2000). The age distribution of the skeletons did not correspond to expected human mortality patterns (for a 30-year life expectancy), since the 0–9-year-old group was underrepresented. Moreover, the total number of subjects was surprisingly low for such a long period of use (at least 400 years). These findings suggest that only specific members of the Xiongnu community were buried in this necropolis.

Samples for DNA analysis were collected from such skeletal elements as astragalus, calcaneus, rib, vertebrae, and teeth during the first year of the excavation; samples from more substantial long cortical bones, such as femur, tibia, and humerus, were collected during the next 2 years. After authorization from the Mongolian authorities, bone samples from 80 skeletal remains (taken, when possible, in duplicate) were transferred to Strasbourg, France, under appropriate storage conditions. On arrival in the laboratory, highly damaged bones (showing extreme fragility and porosity) and severely deteriorated teeth were excluded from the genetic analysis.

DNA Extraction and Purification
DNA was extracted from 79 bone samples corresponding to 62 individuals (some individuals were typed from two independent samples).

To eliminate surface contamination, the outer surface of the bones was removed to almost 2–3 mm of depth with a sanding machine (Dremel). Powdered bone was generated by grinding bone fragments under liquid nitrogen in a 6800 Freezer Mill (Fischer Bioblock) or with a drill fitted with a surgical trepan to avoid overheating.

DNA was carefully extracted according to a published protocol (Fily et al. 1998). In brief, 2 g of the pulverized material was incubated at 50°C overnight in 5 ml of a solution containing 5 mmol EDTA, 2% SDS, 10 mmol Tris HCl (pH 8.0), 0.3 mol sodium acetate, and 1 ml proteinase K/ml. A phenol/chloroform/isoamyl alcohol (25/24/1, v/v) extraction was performed on the supernatant. The aqueous phase was then purified with the Cleanmix kit (Talent), which relies on the strong affinity of DNA to silica in the presence of guanidium thiocyanate. After the elution step with 400 μl sterile water, the DNA was concentrated by passing through a Microcon YM30 filter (Millipore).

To ensure the accuracy and reliability of the results, all samples were amplified (for each marker) at least six times (more when apparent homozygotes were found by autosomal STR analysis) from three independent DNA extracts and, when possible, from two different bones of the same individual.

Autosomal STR Analysis
Autosomal STRs were amplified using the AmpFlSTR profiler Plus kit (Applied Biosystems). Nine STRs (D3S1358, vWA, FGA, D8S1179, D21S11, D18S51, D5S818, D13S317, and D7S820) and the sex determination marker amelogenin were simultaneously amplified.

PCRs were performed according to the manufacturer’s protocol (Applied Biosystems), except that 37 cycles were used instead of 28 in a reaction volume of 10 μl, thereby reducing the volume of the DNA samples and improving the amplification yield.

For three samples (57, 58, and 59), further analyses were performed using the AmpFlSTR SGM Plus kit (Applied Biosystems), which allows the simultaneous amplification of 10 STR loci (4 more than with the previous kit). The genetic relationships between individuals were tested by pairwise comparison of the profiles.

Y Chromosome STR analysis
The DNA of male individuals was analyzed at eight Y chromosome STR loci. Six of them (DYS19, DYS389-II, DYS390, DYS391, DYS393, and DYS385) were coamplified in a multiplex reaction, using the Y-Plex6 kit, according to the manufacturers’ recommendations (ReliaGene Technologies). The two others (YCAII and DYS392) were amplified by singleplex PCR. Primer sequences were those described by de Knijff et al. (1997). For PCR amplification (using a Biometra thermocycler), the following conditions were used: predenaturation at 94°C for 3 min; 30 annealing cycles at 94°C for 30 s, 56°C for 30 s, and 72°C for 90 s; and a final extension at 72°C for 7 min. The allele nomenclature was the one recommended by the International Society of Forensic Genetics (Gill et al. 2001).

mtDNA Analysis
The HVI of the mitochondrial control region was amplified and sequenced from nucleotide positions 16009 to 16390 (Anderson et al. 1981), using primers L15989 and H16410 (Gabriel et al. 2001). When no amplification was obtained with these primers, presumably because of DNA degradation, the additional primers H16239 (Ivanov et al. 1996) and L16190 (Gabriel et al. 2001) were used to amplify the HVI fragment in two steps. PCR was performed with AmpliTaq Gold polymerase, as follows: predenaturation at 94°C for 10 min; 38 annealing cycles at 94°C for 30 s, 48°C or 51°C for 30 s, and 72°C for 45 s; and final extension at 72°C for 10 min. Amplification products were checked on a 1% agarose gel and purified with Microcon-PCR filters (Millipore). The sequence reaction was performed with the same primers on each strand with the ABI Prism BigDye Terminator Cycle Sequencing kit (Applied Biosystems).

Amplification Product Analysis
PCR products were analyzed on an ABI Prism 3100 (Applied Biosystems) automated DNA sequencer. Fragment sizes were determined automatically by use of GeneMapper software and were typed by comparison with allelic ladders (provided in the kits or obtained by the mixture of previously sequenced samples for the most common alleles). mtDNA sequences were analyzed using the Sequencing Analysis and Sequence Navigator software packages.

Measures Taken to Avoid Contamination
Because the possibility of performing genetic analyses had been considered before beginning the archaeological work, precautions were taken to reduce contamination during excavation and curation, for example, samples were handled with gloves by a reduced number of anthropologists wearing face masks. To check for possible modern contamination, the DNA extracted from saliva samples of all people handling the material or working in the laboratory was genetically typed and then compared with the profiling results of all ancient samples.

The entire process of DNA extraction and PCR amplification was performed in an isolated laboratory dedicated to work with ancient DNA, where all staff wore lab coats, face masks, and gloves and where strict cleaning procedures were respected (frequent treatment with DNAse Away and UV light and frequent change of gloves). Autoclaved disposable plasticware, dedicated reagents, and pipettes with aerosol-resistant tips were used; extraction and template blanks were included in every PCR assay; and positive PCRs were never performed. Multiple extractions from the same samples were undertaken at different times, and PCR products were never brought into the ancient DNA laboratory.

Results

Autosomal STR Analysis
Of the 62 individual remains analyzed by multiplex amplification, 8 DNA samples (from graves 32, 34, 51, 60, 78, 83bis, 84bis, and 85 [fig. 2]) appeared severely degraded, since no amplifiable product could be obtained (from at least three independent extracts). One sample (from grave 18) was excluded from further analyses, because it was considered a likely case of contamination (the multiallelic profile matched that of one of the staff, despite multiple independent extractions of this vertebral sample). Four other DNA samples (from graves 26, 27, 67, and 81) were found to contain too few template DNA molecules to provide reproducible results (data not shown). The remaining extracted samples gave 49 more or less complete allelic profiles. Consensus data are reported in table 1. In most cases, these 49 DNA profiles were obtained from diaphyses, but vertebrae provided the genetic profiles in 4 cases, calcaneus in 1 case, and clavicle in 1 case. Long cortical bones (such as femur, tibia, and humerus) thus appeared to be good sources of ancient DNA, whereas rib samples and other thin bones did not. When apparent homozygotes were obtained, amplifications were repeated as many as eight times to avoid the possibility that one allele of an heterozygote was not detected.
Table 1Table 1
Consensus Allelic Profiles of 49 Specimens Recovered from the Egyin Gol Necropolis[Note]

Morphological and molecular typing results for sex determination were in accordance with each other, which is indicative of authentic ancient DNA extracts. We used the amelogenin locus to deduce the sex of six juvenile skeletons for which morphological indicators of sex were absent: three were male (graves 18A, 36, and 84.1), and three others were female (graves 41, 83, and 91). For three other adolescent remains (graves 74, 75, and 76), even molecular determination of sex was ambiguous, despite nearly complete STR profiles.

Comparison of the profiles in pairs allowed us to identify a family composed of the father (grave 57), the mother (grave 59), and one child (grave 58). For these three specimens, all nine STR loci were amplified; and, at each locus, alleles of the child profile could have been assigned to the genotypes from either grave 57 or grave 59 (table 1). To confirm this result, further analysis was performed on these three samples, using the AmpFlSTR SGM Plus Kit (PE Biosystems). This additional analysis confirmed and completed previous results, with one allele contributed from each parent (data not shown), and proved the parental relationships.

It was also possible to determine other familial relationships; for instance, the child from grave 36 is probably the son of the female individual buried near him (grave 35), since the genotypes of these two subjects shared a common allele at each of the nine loci tested. No putative father was found among the profiles of table 1. In the same manner, the genetic profile of the male skeleton retrieved from grave 50 shared one allele at each locus with individuals from graves 46, 52, and 54 and is probably the father of these three individuals. It also shared eight alleles with individuals from graves 48, 65, and 66. Individual profiles from graves 46 and 48 and from 47 and 48 also indicated a parent/child relationship, with one common allele at each locus, as did profiles from the following pairs of graves: 63 and 65, 32A and 33, 70 and 72, 72 and 94, 88 and 94, and 93 and 95. Other individuals could have been closely related parents: the two adolescents from graves 74 and 76 shared one allele at seven or eight of the nine STR markers, as did those from the following sets of graves: 46, 52, and 54; 53 and 54; 50 and 52; 82 and 83; 53 and 69; 65 and 66; and, finally, 94 and 95. The incomplete genotyping of some samples probably hampered the search for other familial relationships.

Y Chromosome STR Analysis
To identify male lineages, an analysis of polymorphic STR systems located on the male-specific part of the Y chromosome was performed. Eight Y-specific STRs were typed and used to construct haplotypes. Of the 35 individuals who were male or whose sex could not be determined, 27 could be typed at more than three loci (table 2). Among them, 18 different haplotypes could be identified (even when incomplete, most haplotypes could be differentiated). The loci DYS385 and DYS389II often could not be amplified, probably because they are expressed in the higher molecular weight range. Such an inverse dependence of the amplification efficiency on the size of the segment to be amplified is typical of DNA retrieved from ancient remains and results from damage and degradation of the DNA.
Table 2Table 2
Y Chromosome STR Haplotypes Determined for 27 of the Ancient Male Specimens[Note]

The most common haplotype was observed in six male specimens buried in the C sector (graves 46, 47, 50, 52, 53, and 54 [fig. 3; table 2]), suggesting a grouping of individuals belonging to the same paternal lineage. Three of these individuals (graves 46, 52, and 54), shared, in addition, the same mtDNA sequence (see below); they were therefore considered to be brothers. Since their autosomal allelic profile showed one allele in common, at each locus, with that of the male from grave 50, the latter was considered to be their father. The others (from graves 47 and 53) were probably more distant paternal relatives (half-brothers, nephew and uncle, or grandfather and grandson).

The study of these uniparentally inherited STR markers also showed that the individual in grave 58 had the same Y haplotype as his putative father in grave 57: all seven regions of the Y chromosome tested matched, confirming the autosomal typing results (paternity). Three other adult individuals, buried in the northern part of the necropolis (graves 70, 72, and 73) were found to share an identical six-locus haplotype (fig. 3; table 2). Two of them (from graves 70 and 72), who shared at least one common allele at each locus (see “Autosomal STR Analysis” section), may be considered to be a father and son. Close to them, two other specimens (graves 88 and 94) sharing an identical six-locus haplotype (table 2) and half of their autosomal alleles (table 1) were supposed to be genetically linked by a father/son relationship.

Some DNA samples failed to yield any amplification results. Among them were the DNA extracted from the adolescent remains from graves 74 and 75. Since the sex of both specimens could not be clearly established (either morphologically or genetically), one can suppose that these adolescents were female individuals. Conversely, the individual from grave 76 gave incomplete but consistent results, suggesting that this specimen was a male.

mtDNA Analysis
Sequence variation in the HVI was investigated in 56 of the ancient specimens. Reproducible HVI sequences were obtained for 46 of them. Among these 46 individuals, a total of 28 different sequences, defined by 44 variable positions, were identified (table 3). The most frequent mtDNA type was scored in four individuals: three of them (from graves 46, 52, and 54) belong to the C sector and were considered to be brothers, since they also shared an identical Y haplotype (table 2); the fourth (from grave 57) was the father of the little family identified in the middle of the necropolis. Twelve other mtDNA types were shared by at least two individuals, the remaining 15 mtDNA types being represented by just one individual. Two of these unique mtDNA types (from graves 48 and 61) differed from the most frequent one (from graves 46, 52, 54, and 57) by a single mutation and may be considered to arise from the same maternal lineage.
Table 3Table 3
mtDNA HVI Sequences of the Egyin Gol Specimens[Note]

Differences between sequences could be mostly attributed to transitional substitutions (90%) and concerned mainly the pyrimidines; however, at position 183, transversions occurred in several different sequences. The C→T transition at position 16223 (nucleotide position in the reference sequence of Anderson et al. [1981]) was shared by most of the ancient specimens (35/46 individuals), as was the T→C transition at position 16362 (27/46 individuals). Two instances of a transition and transversion at the same site were also observed: position 16129 showed both G→A and G→C mutations, and position 16232 showed both C→T and C→A mutations, as previously reported by Kolman et al. (1996). Insertion of a C residue was found once, between positions 16193 and 16194.

Polymorphic sites shared by two individuals allowed us to confirm or to reconsider close genetic affinities between some specimens. For instance the individuals from graves 59 and 58 showed an identical HVI sequence, confirming the maternal relationship deduced from the autosomal STR typing. The child of grave 36 had the same mtDNA sequence as her presumed mother (grave 35) and the female specimen from grave 37A (double grave) (fig. 4). Other complete matches were noted between individuals from graves 83bis and 91; 65 and 77; 32A and 72; 28, 73, and 74; 70, 88, and 94; 53 and 69; 83 and 82; and 76 and 86 (table 3), even though autosomal STR data did not always clearly show any parental relationship. The male individuals from graves 88 and 94, who were thought to be a father-son pair, since they share at least one allele at every locus (see the “Autosomal STR Analysis” section) and an identical five-locus haplotype, may in fact be brothers, since their mitochondrial haplotype is identical. Similarly, the two adolescents from graves 74 and 76 who were thought to be siblings are obviously not, since they do not share an identical HVI sequence.

Heteroplasmies were found within the mtDNA sequences of individuals from graves 39, 41, and 18A. Specimens from graves 39 and 18A were grouped, respectively, with individuals from graves 49 and 27, since the remaining nucleotides perfectly matched each other; a comparison with nuclear data to decide whether or not graves 18A and 27 contained maternal relatives was not possible, because of the incompleteness of the autosomal DNA profiles. Individuals from graves 39 and 49 were considered to be maternal relatives on the basis of the genotyping results.

Although the mtDNA sequences obtained could not be assigned with certainty to mtDNA haplogroups (since they encompassed only the HVI of the control region), three (A, C, and D) of the four major haplogroups observed in Native American (Torroni et al. 1993) and Siberian (Starikovskaya et al. 1998; Schurr et al. 1999) populations were detected in the ancient samples tested (with haplogroup D being the most prevalent). A few sequences belonging to subclusters B4b, D5 or D5a, F1b, J1, G2a, U2 or U5a1a—and some that probably belonged to cluster M—were also observed (Richards et al. 2000; Yao et al. 2002). No member of the major European cluster H, which occurs in >40% of most European populations (Richards et al. 1996) was found. Interestingly, some of the haplotypes reported here are similar to those found in previous studies of the area (Kolman et al. 1996; Comas et al. 1998).

Discussion

In the present article, partial genealogical reconstruction was obtained using biparental, paternal, and maternal genetic systems in a sample of 62 human skeletal remains exhumed from a cemetery dating from >2,000 years ago. To the best of our knowledge, no equivalent molecular analysis has been undertaken so far. Such a study was possible because the Egyin Gol necropolis was mainly composed of relatively well-preserved skeletons. Indeed, the climatic conditions (cold and dry) and the archaeological context (architectural structure of the graves) encountered at this site had undoubtedly protected the recovered specimens against DNA degradation. Regarding DNA retrieval, PCR amplification results showed three kinds of samples, as described elsewhere (Burger et al. 1999): (i) samples in which sufficient DNA molecules were preserved and for which definite and reproducible genotypes or haplotypes could be determined; (ii) samples in which DNA could be detected sporadically but without reproducible results; and (iii) samples in which no DNA (or almost none) was preserved.

The choice of autosomal STR markers as a first approach for analyzing individual remains of the Egyin Gol necropolis was based on (i) their high discriminatory power, in comparison with mtDNA analysis, as a means to investigate close familial relationship; (ii) their small size, which facilitates amplifications in old or degraded DNA; (iii) the possibility of amplifying several of them simultaneously from minute amounts of DNA; and (iv) their ability to indicate a result’s authenticity (notably by comparing amplified products to the profiles of all persons involved in the investigation).

The multiallelic DNA profiles obtained for 49 of the ancient specimens were compared with each other. A direct parental link was considered plausible if a pair shared an allele at each of the nine loci tested; in this manner, a total of nine pairs were identified as representing possible parent-child relationships. Moreover, we could identify three children thought to have a common parent (individuals from graves 46, 52, and 54, fathered by the individual from grave 50). The traditional parentage trio, with both parents, was encountered only once. This is not surprising, considering that not all profiles were complete and that the number of inhumations that occurred throughout more than four centuries (the duration of the necropolis’s use) is relatively small.

To verify the accuracy of the biological relationships deduced from autosomal STR data and to gain a higher power of discrimination, the study of nonrecombining marker systems, such as Y chromosome STRs and mtDNA was undertaken. These analyses confirmed the child or sibling status of some individuals (graves 35 and 36; 57, 58, and 59; 50, 46, 52, and 54; 70 and 72; and 88 and 94) and the close genetic relationship between some others (graves 53 and 54, 82 and 83, and 53 and 69). In one case, however, discordant results between the biparental and the two uniparental systems were observed. Indeed, the multiallelic profiles obtained from the skeletal remains recovered from graves 94 and 95 supported a biological relationship between them. Nevertheless, because these two “relatives” had neither the same Y haplotype nor the same mtDNA sequence, we had to consider the possibility that the two specimens were genealogically unrelated (unless each of the two sets of parents were siblings). For other pairs of individuals (from graves 47 and 48, 63 and 65, 65 and 66, and 93 and 95), the validation of a close parental link was not possible without knowing which was the parent and which was the child.

Among other results, the Y chromosomal STR analysis revealed that one of the defined topographical sectors was exclusively composed of males of the same patrilineage (the individuals from graves 48, 49, and 51 were female, and no bone samples were available for specimens from graves 43–46A). The other males found to share the same Y haplotype (graves 57 and 58; 88 and 94; and 70, 72, and 73) were also buried close to each other. Such a grouping of male relatives has never been demonstrated before for ancient specimens and provides an insight into the funeral practices of ancient Eurasian tribes.

The maternal genetic inheritance, which was tested through sequencing of the mtDNA HVI, revealed some biological links. For example, children from graves 83bis and 91 might be considered to be relatives, since they share an identical HVI sequence (table 3). Other maternal links were revealed, such as those between individuals from graves 65 and 77; 32A and 72; 28, 73, and 74; 70 and both 88 and 94; and 18A and 27. In some cases, heteroplasmies were observed. The possibility that these heteroplasmies resulted from contamination was invalidated by the fact that identical results were obtained from DNA samples extracted and amplified in triplicate at different time intervals.

Nevertheless, the pitfall of contamination from extraneous human DNA is a major concern for researchers working with human remains (Handt et al. 1994; Kolman and Tuross 2000) and should not be underrated. Since some erroneous ancient DNA results have been published, a number of “criteria of authenticity” need to be fulfilled before results from ancient DNA analyses can be taken to be genuine (Handt et al. 1994; Cooper and Poinar 2000). In the present study, extensive precautions (described in the “Material and Methods” section) were taken to avoid the amplification of contaminating contemporary DNA molecules. Despite the fact that not all reported criteria of authenticity could be met, the possibility that our data arose from contaminating DNA was considered highly unlikely for the following reasons: (i) reproducible PCR results were obtained from multiple extractions and amplifications of the same samples made at different times; (ii) multiallelic profiles were not mixtures of different individuals’ DNA and were not found to correspond to someone involved in the present work (except once); (iii) the results of both sex typing methods (morphologic and genetic) were in accordance with each other; (iv) an inverse relationship between amplification efficiency and length of the amplification products was observed, especially with STR markers; (v) the crystallinity index and the carbon/nitrogen ratio determinations indicated no significant alteration of the bones (Nielsen-Marsh et al. 2000); (vi) a concordance was observed between data obtained with the markers inherited biparentally, paternally, and maternally; (vii) mtDNA analysis of the ancient sample revealed that most of the haplogroups present were of Asian origin and that European maternal lineages identical to those of the excavators or laboratory personnel (all of whom were of European origin) were absent; (viii) the 16223 thymine-cytosine transition was found in 76% of the ancient Egyin Gol samples, a result close to that of Kolman et al. (1996), who found it in 65% of Mongolian samples (compared with 7% of European samples), and (ix) the level of genetic diversity detected in the protohistoric population, as well as some of the haplotypes reported, are similar to those obtained in modern Mongolian populations (authors' unpublished data; Kolman et al. 1996).

Nevertheless, the test that confers the greatest level of robustness is duplicate analysis by two independent laboratories. This was not feasible in the present study, because of the large number of subjects tested. Cloning of the PCR products was not conceivable for the same reason and is not really adapted to the study of STR. Methods such as amino acid racemization and DNA quantitation were not applied, mainly because they do not allow the distinction between contaminated and uncontaminated samples (Kolman and Tuross 2000). Moreover, the proposed use of amino acid racemization to estimate DNA survival in archaelogical bones is challenged by some authors (Collins et al. 1999). The fact that the multiallelic profiles were repeatable from the same—and different—DNA extracts of a specimen allowed us to consider that the number of starting templates was high enough to obtain reliable results and to analyze mtDNA (for which the contamination problem is worse). On the other hand, we subscribed to other important criteria, such as the reiteration of the extraction and amplification steps, the sexing of the samples tested, and, foremost, the use of a molecular combined approach. Moreover, the bone samples studied are not fossil remains and, consequently, are not prone to high rates of DNA alteration.

A majority (89%) of the Xiongnu sequences can be classified as belonging to an Asian haplogroup (A, B4b, C, D4, D5 or D5a, or F1b), and nearly 11% belong to European haplogroups (U2, U5a1a, and J1). This finding indicates that the contacts between European and Asian populations were anterior to the Xiongnu culture, and it confirms results reported for two samples from an early 3rd century b.c. Scytho-Siberian population (Clisson et al. 2002).

The genetic data obtained in the present study, in addition to the topographical and radiocarbon data, suggested hypotheses concerning the social history of the necropolis. Around the 3rd century b.c., the grave of an adult male (grave 28) had been dug on the southern part of the Egyin Gol valley (A sector). At a short distance from him, a privileged man was also buried (grave 27), as were other individuals, including those found in double graves (graves 32/32A, 33/33A, 37/37A, and 38/38A). Some of these surrounding graves could be sacrificial burials, as has been reported elsewhere for one of them (Murail et al. 2000). This tradition of having double graves near an opulent one in cemeteries containing individuals of high social class is well documented, notably in the Sakka (another group of nomadic people of the Eurasian steppes) and the Pazyryk cultures (Francfort et al. 2000). This ritual, at the first developmental step of the cemetery, suggests that the cultural influence of the “old Scythian spirit” was already present in some nomadic families at the beginning of the Xiongnu empire. Although close genetic relationships could not be clearly established between ancient specimens of the A sector (because of the lack of amplification results), a parent/child link was nevertheless shown between individuals in graves 32A and 33, suggesting the possibility of burials based on familial relationships.

Some years later, a new sector of interment seems to have been created (sector B). In this second sector, some individuals shared mtDNA sequences with individuals from the most ancient graves. (Thus, individuals from graves 73 and 74, 18A, and 72 were found to share mtDNA sequences with those from graves 28, 27, and 32A, respectively.) This result suggests that maternal relatives of the individuals first buried might participate in the extension of a new cemetery area. It should be noted that, except for individuals from grave 18A, these maternal relatives were buried close to each other. Two of them (from graves 73 and 72) were from the same paternal lineage. We can imagine that the creation of this new burial area could be the result of tensions between members of the ruling family. The fact that no double graves were built could reflect the cultural rupture with the “old Scythian spirit.”

From the 2nd century b.c. to the 1st century a.d., the social organization of the necropolis cannot be clearly deduced from the genetic data. Male and female individuals from different paternal and maternal lineages were buried from south to north in the A sector and from north to south in the B sector. Genetic analyses revealed several familial groups buried close to each other, notably, one consisting of a father, mother, and son. Since these three individuals were probably not buried at the same time (the son appeared to be as old as his parents), the family grouping was organized at least one generation earlier. From these results, it seems reasonable to speculate that the Xiongnu buried relatives together although the practice was not systematic. During this period the site might have been the cemetery of a social group with a significant genetic diversity, of which only certain members were buried (the small number of inhumations that occurred throughout more than three centuries suggests that this Xiongnu tribe probably had other sepulchral places).

After the fusion of the A and B sectors, new graves were dug in the west. These graves correspond to a group of genetically linked individuals, since they belong to a single paternal lineage. Interestingly, this paternal lineage has been, at least in part (6 of 7 STRs), found in a present-day Turkish individual (Henke et al. 2001). Moreover, the mtDNA sequence shared by four of these paternal relatives (from graves 46, 52, 54, and 57) were also found in a Turkish individuals (Comas et al. 1996), suggesting a possible Turkish origin of these ancient specimens. Two other individuals buried in the B sector (graves 61 and 90) were characterized by mtDNA sequences found in Turkish people (Calafell 1996; Richards et al. 2000). These data might reflect the emergence at the end of the necropolis of a Turkish component in the Xiongnu tribe.

In conclusion, our study shows how the use of genetic markers of different mutability might provide an insight into the history of past necropolises. It also provides genetic data on ancient Eurasian specimens that could help to confirm or disprove models developed from modern genetic data to explain population history. Finally, it provides an excellent tool to select samples of interest for interpopulation analyses.

Acknowledgments

This research was supported by a grant (Aide à Projet Nouveau) of the Centre National de la Recherche Scientifique (to E.C.). Additional support was provided by the Institut of Legal Medicine of Strasbourg (postdoctoral contract) to C.K.-T. We are indebted to P. Blandin, F.-X. Ricaut, and E. Tissier for their help in this work. We would also like to thank P.-H. Giscard, for the management of the excavations, P. Murail, for the classic anthropological study (ages and sex), and J. P. Verdier, for the topographic study. The staff of J. Jaubert is also gratefully acknowledged for help with the geological study of the valley and the discovery of permafrost. We thank V. Balter and A. Person for the crystallographic studies. The fieldwork was made possible thanks to the Ministère des Affaires Etrangères (France) and UNESCO. We also thank S. Erdenebaatar and D. Turbat from the Mongolian Academy for the codirection, with P.-H. Giscard and E. Crubézy, of the team in the field.

References
Anderson S, Bankier AT, Barrell BG, de Bruijn MH, Coulson AR, Drouin J, Eperon IC, Nierlich DP, Roe BA, Sanger F, Schreier PH, Smith AJ, Staden R, Young IG (1981) Sequence and organisation of the human mitochondrial genome. Nature 290:457–465 [PubMed].
Burger J, Hummel S, Herrmann B, Henke W (1999) DNA preservation: a microsatellite-DNA study on ancient skeletal remains. Electrophoresis 20:1722–1728 [PubMed].
Calafell F, Underhill P, Tolun A, Angelicheva D, Kalaydjieva L (1996) From Asia to Europe: mitochondrial DNA sequence variability in Bulgarians and Turks. Ann Hum Genet 60:35–49 [PubMed].
Clisson I, Keyser C, Francfort HP, Crubézy E, Samashev Z, Ludes B (2002) Genetic analysis of human remains from a double inhumation in a frozen kurgan in Kazakhstan (Berel site, early 3rd century BC). Int J Legal Med 116:304–308 [PubMed].
Collins MJ, Waite ER, van Duin AC (1999) Predicting protein decomposition: the case of aspartic-acid racemization kinetics. Philos Trans R Soc Lond B Biol Sci 354:51–64 [PubMed].
Comas D, Calafell F, Mateu E, Pérez-Lezaun A, Bertranpetit J (1996) Geographic variation in human mitochondrial DNA control region sequence: the population history of Turkey and its relationship to the European populations. Mol Biol Evol 13:1067–1077 [PubMed].
Comas D, Calafell F, Mateu E, Pérez-Lezaun A, Bosch E, Martinez-Arias R (1998) Trading genes along the Silk Road: mtDNA sequences and the origin of Central Asian populations. Am J Hum Genet 63:1824–1838 [PubMed].
Cooper A, Poinar HN (2000) Ancient DNA: do it right or not at all. Science 289:1139 [PubMed].
Crubézy E, Martin H, Giscard PH, Batsaikhan S, Erdenebaatar D, Verdier JP, Maureille B (1996) Funeral practices and animal sacrifices in Mongolia at the Uigur period: archeological and ethno-historical study of a kurgan in the Egyin Gol valley. Antiquity 70:891–899.
Crubézy E, Masset C, Lorans E, Perrin F, Tranoy L (2000) Archéologie funéraire. Editions Errance, Paris.
de Knijff P, Kayser M, Caglià A, Corach D, Fretwell N, Gehrig C, Graziosi G et al (1997) Chromosome Y microsatellites: population genetic and evolutionary aspects. Int J Legal Med 110:134–149 [PubMed].
Fily ML, Crubézy E, Courtaud P, Keyser C, Ebrard D, Ludes B (1998) Analyse paléogénétique des sujets de la grotte sépulcrale d’Elzarreko Karbia (Bronze ancien, Pays Basque). CR Acad Sci III 321:79–85 [PubMed].
Francfort HP, Ligabue G, Samashev Z (2000) La fouille d’un kourgane Scythe gelé du IVe siècle av. notre ère à Berel dans l’Altaï (Kazakhstan). Académie des Inscriptions et Belles Lettres, Paris, pp 775–806.
Gabriel MN, Huffine EF, Ryan JH, Holland MM, Parsons TJ (2001) Improved mtDNA sequence analysis of forensic remains using “mini-primer set” amplification strategy. J Forensic Sci 46:247–253 [PubMed].
Gill P, Brenner C, Brinkmann B, Budowle B, Carracedo A, Jobling MA, de Knijff P, Kayser M, Krawczak M, Mayr WR, Morling N, Olaisen B, Pascali V, Prinz M, Roewer L, Schneider PM, Sajantila A, Tyler-Smith C (2001) DNA Commission of the International Society of Forensic Genetics: recommendations on forensic analysis using Y-chromosome STRs. Int J Legal Med 114:305–309 [PubMed].
Handt O, Höss M, Krings M, Pääbo S (1994) Ancient DNA: methodological challenges. Experientia 50:524–529 [PubMed].
Henke J, Henke L, Chatthopadhyay P, Kayser M, Dülmer M, Cleef S, Pöche H, Felske-Zech H (2001) Application of Y-chromosomal STR haplotypes to forensic genetics. Croat Med J 42:292–297 [PubMed].
Hummel S, Bramanti B, Schultes T, Kahle M, Haffner S, Herrmann B (2000) Megaplex DNA typing can provide a strong indication of the authenticity of ancient DNA amplifications by clearly recognizing any possible type of modern contamination. Anthropol Anz 58:15–21 [PubMed].
Ivanov PL, Wadhams MJ, Holland MM, Weedn VW, Parsons TJ (1996) Mitochondrial DNA sequence heteroplasmy in the Grand Duke of Russia Georgij Romanov establishes the authenticity of the remains of Tsar Nicholas II. Nat Genet 12:417–420 [PubMed].
Kolman CJ, Sambuughin N, Bermingham E (1996) Mitochondrial DNA analysis of Mongolian populations and implications for the origin of New World founders. Genetics 142:1321–1334 [PubMed].
Kolman CJ, Tuross N (2000) Ancient DNA analysis of human populations. Am J Phys Anthropol 111:5–23 [PubMed].
Konovalov NB (1976) Khuny V Zabayikale Akademia nauk SSSR, Sibirskoe Otoelenie, Buriatskue Knijnoe izdatelstvo, Ulan-Ude.
Leonard JA, Wayne RK, Cooper A (2000) Population genetics of ice age brown bears. Proc Natl Acad Sci USA 97:1651–1654 [PubMed].
Marx I (2000) The Xiongnu culture—third century BCE. http://www.silk-road.com/toc/index.html (accessed June 12, 2003).
Minajev SS (1996) Les Xiongnu. Dossiers d’Archéologie 212:74–85.
Murail P, Bruzek J, Braga J (1999) A new approach to sexual diagnosis in past populations: practical adjustments from Van Vark’s procedure. Int J Osteoarchaeol 9:39–53.
Murail P, Crubézy E, Martin H, Haye L, Bruzek J, Giscard PH, Turbat T, Erdenebaatar D (2000) The man, the woman and the hyoid bone: from archaeology to the burial practices of the Xiongnu people (Egiin Gol valley, Mongolia). Antiquity 74:531–536.
Nielsen-Marsh C, Gernaey A, Turner-Walker G, Hedges R, Pike A, Collins M (2000) The chemical degradation of bone. In: Cox M, Mays S (eds) Human osteology in archeology and forensic science. London, pp 439–454.
Person A, Bocherens H, Mariotti A, Renard M (1996) Diagenetic evolution and experimental heating of bone phosphate. Palaeogeography, Palaeoclimatology, Palaeoecology 126:135–150.
Richards M, Côrte-Real H, Forster P, Macauley V, Wilkinson-Herbots H, Demaine A, Papiha S, Hedges R, Bandelt HJ, Sykes B (1996) Paleolithic and Neolithic lineages in the European mitochondrial gene pool. Am J Hum Genet 59:185–203 [PubMed].
Richards M, Macaulay V, Hickey E, Vega E, Sykes B, Guida V, Rengo C, et al (2000) Tracing European founder lineages in the Near Eastern mtDNA pool. Am J Hum Genet 67:1251–1276 [PubMed].
Rudenko SI (1970) Frozen tombs of Siberia: the Pazyryk burials of Iron Age horsemen. London, Dent and Son.
Schultes T, Hummel S, Herrmann B (2000) Ancient DNA-typing approaches for the determination of kinship in a disturbed collective burial site. Anthropol Anz 58:37–44 [PubMed].
Schurr TG, Sukernik RI, Starikovskaya YB, Wallace DC (1999) Mitochondrial DNA variation in Koryaks and Itel’men: population replacement in the Okhotsk Sea-Bering Sea region during the Neolithic. Am J Phys Anthropol 108:1–39 [PubMed].
Starikovskaya YB, Sukernik RI, Schurr TG, Kogelnik AM, Wallace DC (1998) mtDNA diversity in Chukchi and Siberian Eskimos: implications for the genetic history of ancient Beringia and the peopling of the New World. Am J Hum Genet 63:1473–1491 [PubMed].
Stone AC, Stoneking M (1999) Analysis of ancient DNA from a prehistoric Amerindian cemetery. Philos Trans R Soc Lond B Biol Sci 354:153–159 [PubMed].
Torroni A, Schurr TG, Cabell MF, Brown MD, Neel JV, Larsen M, Smith DG, Vullo CM, Wallace DC (1993) Asian affinities and continental radiation of the four founding Native American mtDNAs. Am J Hum Genet 53:563–590 [PubMed].
Yao YG, Kong QP, Bandelt HJ, Kivisild T, Zhang YP (2002) Phylogeographic differentiation of mitochondrial DNA in Han Chinese. Am J Hum Genet 70:635–651 [PubMed].