View metadata, citation and similar papers at core.ac.uk brought to you by CORE

provided by Elsevier - Publisher Connector

FEBS Letters 580 (2006) 6242–6246

Genetic analysis on remains excavated from Qilang Mountain Cemetery in Qahar Right Wing Middle Banner of Inner

Yu Changchuna,b, Xie Lib, Zhang Xiaoleia, Huia,c,*, Zhu Honga a Ancient DNA Laboratory, Research Center for Chinese Frontier Archaeology, University, Changchun 130012, P.R. b College of Life Science, Jilin Normal University, Siping 136000, P.R. China c College of Life Science, Jilin University, Changchun 130012, P.R. China

Received 17 July 2006; revised 5 October 2006; accepted 12 October 2006

Available online 20 October 2006

Edited by Takashi Gojobori

three southward migrations, and finally founded the Northern Abstract Sixteen sequences of the hypervariable segment I (HVS-I, 16039–16398) in mtDNA control region from ancient Wei Dynasty and controlled the northern region of China Tuoba Xianbei remains excavated from Qilang Mountain (386–534 A.D.). By the time of the they had Cemetery were analyzed. In which, 13 haplotypes were found largely merged with populace. by 25 polymorphic sites. The haplotype diversity and nucleotide At present, the origin of Tuoba Xianbei ethnic group and its diversity were 0.98 and 0.0189, respectively, and the mean of impacts on modern minorities of northern China are not yet nucleotide number differences was 6.25. analysis clear. It will be very helpful for understanding in the origin, indicates these remains mainly belong to haplogroup C formation and development process of northern minorities of (31.25%) and D (43.75%). According to the published data were China by analyzing the data of mtDNA in Tuoba Xianbeis. considered, we can suggest that the Tuoba Xianbei presented a In this study, we examined the sequences of HVS-I from 16 close genetic affinity to Oroqen, Outer Mongolian and Ewenki Tuoba Xianbei remains excavated in Qilang Mountain Ceme- populations, especially Oroqen. 2006 Federation of European Biochemical Societies. Published tery that was indicated in Fig. 1, and compared the data with by Elsevier B.V. All rights reserved. those of relevant populations.

Keywords: Qilang Mountain; Tuoba Xianbei; Mitochondrial DNA 2. Materials and methods

2.1. Sampling The thighbone samples from 18 ancient Xianbei remains, provided by Dr. Wei-Jian ( Institute of Archaeology), which 1. Introduction were excavated from Qilang Mountain Cemetery in Qahar Right Wing Middle Banner of Inner Mongolia from 1995 to 1996. During this The PCR and ancient DNA technique have made it possible period, 20 graves were excavated, and out of which 18 remains were to identify the trace amount of DNA which still survived in collected (numbered ZQ01–18). According to Wei et al. [10], these remains occupied between 4th and 5th centuries A.D., belong to archaeological remains [1–3]. Among the most informa- Tuoba tribe of Xianbei alliance. tive genetic systems, mitochondrial DNA is the preferred target for recovering DNA from archaeological remains and 2.2. Contamination controls characterizing the structure of a gene pool, which well-solved All stages of the work were carried out in term of aDNA-specific many important questions, such as human origin, migration requirements [11]. First, before beginning our work, 2 decontamina- and evolution, as well as population affinities and so on [4–9]. tion methods, such as irradiating the samples 2–3 h with UV-light irra- diation (254 nm) and removing the sample surface 3–4 mm depth with In ancient China, the Xianbei was another significant noma- dentistry drill, were employed to remove the surface DNA deriving dic population succeeded by , which actually con- from other workers who had handled the samples before they reached sisted of a federation of sizeable non-Han groups, and the our laboratory; secondly, the preparation of bone powders and the Tuoba Xianbei was a very important sub-tribe in Xianbei Alli- extraction and amplification of DNA were carried out in separate ance. The Tuoba Xianbeis are said to be the northern-most dedicated laboratories, routing sterilized by treatment with DNase away (Molecular Bioproduct, San Diego, CA), bleach and UV-light branch of the Xianbei nomads, originally dwelled to the north- irradiation (254 nm), and the gloves, face masks and full-body suits eastern-most of all Xianbei tribles, in the region of the Great must be used during all the handling processes; thirdly, the blank- Xianbei Mountains, which is now the Erguna River valley of extractions, zero-controls and negative-controls must be carried out the northern section of the Great Xing’ Range. The lan- during the extraction and amplification of DNA, and fourthly, all results must be identical in two independent extractions and two guage spoken by Tuoba Xianbei was almost the same with independent amplifications, and compared with the sequences from Mongolic. According to Chinese historical documents, from the handlers. the 1st to the 3rd centuries, the Tuoba tribe had undergone 2.3. DNA extraction and amplification The method of Kalmar et al. [12] was used for the bone powder preparation. The manual protocol of GENECLEAN Kit for ancient *Corresponding author. Fax: +86 431 8498031. DNA (BIO101, USA) was used for DNA extraction. As for PCR E-mail address: [email protected] (Z. Hui). amplification, the technique of Adcock et al. [13] was adopted, and

0014-5793/$32.00 2006 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved. doi:10.1016/j.febslet.2006.10.030 Y. Changchun et al. / FEBS Letters 580 (2006) 6242–6246 6243

Fig. 1. The location of Qahar Right Wing Middle Banner Qilang Mountain Cemetery.

modified slightly. Four sets of primers for PCR were L16035–16055/ puted by Arlequin 2000. The net genetic distances and phylogenetic 16162–H16142, L16112–16131/H16237–16218, L16166–16185/16307– tree were carried out using MEGA 3.0. To visualize the Fst distances H16286 and L16261–16281/H16398–16379. The PCR conditions were among populations, the multidimensional scaling (MDS) analysis was described as follow: initial denaturizing at 94 C for 3 min, followed by performed by SPSS 12.0. For comparison, the published data from 40 cycles of denaturizing at 94 C for 1 min, annealing at 55 C for populations of Evenk [16], Inner Mongolian, Daur, Oroqen and Ewen- 45 s, and extension at 72 C for 1 min, followed by the final extension ki [17], Outer Mongolian [15], Buryat [18], Kazak [19], Uzbek [20], 10 min. The quality of amplification was estimated on a 2% agarose Turkey [21], (Zhangjiang Han and Wuhan Han [22]), gel, and the PCR products were purified using ‘QIAEX II’ gel extrac- Northern Han (Qingdao Han, Fengcheng Han [22] and Xi’an Han tion kit (QIAGen, Germany). [23]) were used.

2.4. Cloning and sequencing The PCR products were cloned using TOPO TA Cloning kit 3. Results and discussion (Invitrogen) according to the manufacture’s instruction, which can be used for assessing for damage, contaminations and jumping PCR [11]. Screening of white recombinant colonies was accomplished by 3.1. Sequence diversity and the mean number of pairwise PCR. The colonies were transferred into a 30 ll reactions mix with differences 67 mM Tris–HCl (pH 8.8), 2 mM MgCl2,1lM of each primer (M13 Sixteen sequences of HVS-I (position 16039–16398) were forward and reverse universal primers), 0.125 mM of each dNTP, obtained from 18 Tuoba Xianbei remains. Of these 18 subjects, and 0.75 U of Taq polymerase. After 5 min at 92 C, 30 cycles of PCR (30 s at 90 C, 1 min at 50 C, 1 min at 72 C) were carried out, we could not obtained PCR products from ZQ07. The the clones with an insert of the expected size were identified by aga- sequence from ZQ01 was excluded, because it was ambiguous rose-gel electrophoresis. After purification of these PCR products with by the clone test. Among the 16 sequences, 13 different hyplo- Microcon PCR diveces (Amicon), a volume of 1.5 ll was cycle-se- types were identified by 25 polymorphic sites relative to the ref- quenced using the BigDye Terminater cycle sequencing kit (Applied erence sequence [24], and thereinto, ZQ02 and ZQ10, ZQ08 Biosystems). Six clones of the same PCR products were sequenced, and finally, the sequences were analyzed on an Applied Biosystems and ZQ12, as well as ZQ16 and ZQ17 shared the same haplo- DNA 377 sequencer. types. The haplotype diversity and the mean number of pair- wise differences are 0.98 and 6.25. According to Pakendorf 2.5. RFLP analysis et al. [9], the gene diversity of the Tuoba Xoanbai (0.98) fell In order to confirm the haplogroup status of Tuoba Xianbei sam- into those of the steppe populations (around 0.99) and the ples, three restriction sites (AluI 5176, DdeI 10394 and AluI 10397) Siberian populations (0.93–0.95). in coding region were further detected. The primers for AluI 5176 were cited from Batista et al. [14], and for DdeI 10394 and AluI 10397 came from Kolman et al. [15]. 3.2. Nucleotide diversity Of these 25 mutation sites, 24 were transitions, only 1 was 2.6. Data analysis the transversion (16183 A fi C). The ratio of transition:trans- The editing of sequences was performed by Clustal X1.83. The final information used for data analysis about each individual was the HVS- version is 24:1, which agrees with other data obtained from I of 360 bp, nucleotide positions 16039–16398. Fst genetic distances, humans [25]. Among the transitions, there were a high propor- Genetic diversity, the mean number of pairwise difference were com- tion of C M T substitutions compared to A M G, and the 6244 Y. Changchun et al. / FEBS Letters 580 (2006) 6242–6246

A fi G only occurred one time at position 16175 in ZQ03 (Table 1). The nucleotide diversity is 0.0189, located between those of eastern Asian populations [19], which agreed with the conclusion of gene diversity.

3.3. Haplogroup analysis n.d. B B A According to Wallace et al. [26], in Asia, the in- clude A, B, C, D, E, F and G, etc. Haplogroups A and B lack both the site DdeI 10394 and AluI 10397, whereas C and D I10397 contain these sites. Comparing haplogroup G and D, D lack Alu AluI 5176 site, and G has it. For haplogroup A, its character- istic mutation sites are at positions 16362 (T to C), 16319 (G to A), 16290 (C to T) and 16223 (C to T), for B, at 16217 (T to C) I10394 and 16189 (T to C), for C, at 16327 (C to T), 16298 (T to C) ÀÀ ++D ++D ++D ++D ++D ++D ++D Dde and 16223 (C to T), and for D, at 16362 (T to C) and 16223 (C to T). Together with Kivisild et al. [27], Kong et al. [17] and Yao et al. [22], we can conclude that the samples ZQ02, I5176 05, 09, 10, 15, 16 and 17 belong to D, ZQ03, 04, 08, 11 and ÀÀ À ÀÀ À À À À À À À Restriction sitesAlu À Haplo-groups 12 belong to C, ZQ13 and 14 belong to B, and ZQ06 belong to A, but the ZQ18 was not determined (Table 1). The haplo- group frequencies of the Tuoba Xianbei were 43.75% D, 31.25% C, 12.5% B, 6.25% A and 6.25% other. By comparing with the relevant populations, we found that the Tuoba Xian- bei ethnic group presented a closer similarity to the Oroqen population (Table 2).

3.4. Population comparison To compare with the populations published, the Fst genetic distances and the matrix of significant Fst P value (significance level = 0.05) were calculated and listed in Table 3. The Tuoba –––––––C+ + + –––––––+ –––C–––C –––––––C ––––––AC –––C–––C –––C–––C Xianbei presented the shortest Fst genetic distance to the Oro- qen (0.017), Outer Mongolian (0.019) and Ewenki (0.020), the longest to the Turkeys (0.146), implying that there was a closer genetic affinity among the four populations. According to his- tory documents, the forefathers of the Oroqen, Ewenki and Mongolian can be traced to the ancient ‘‘Shiwei’’ alliance, and the Shiwei and the Tuoba Xianbei all originated from

–––––––C––––+ ‘‘Eastern ’’ population [28,29], which indicated that our re- sults agreed with history records. The genetic affinity among populations can be reflected well in the phylogenetic tree. An unrooted Neighbor-Joining tree of the 14 populations, based on the net average distance (the data not shown), was built using MEGA 3.0. The Xianbei population was located between the Oroqen and Outer Mongolian populations in the tree (Fig. 2). The MDS analysis can really reflects mtDNA similarities and dissimilarities among populations, and the Fst distances and net distances among the Xianbei and the populations com- pared can be interpreted well in MDS analysis plot (Fig. 3). In MDS map, the Xianbai population is situated very close to the Oroqen, Ewenki and Outer Mongolian populations, which supported the conclusion of the population tree.

4. Conclusion

According to our results and the published data [9,19], 9022377788811125789911246 3469624534937835480819762 0111111111122222222233333 we found that the Tuoba Xianbei bore the typical genetic characteristics from Northeastern Asia, and present some closer genetic affinities to the populations of Oroqen, Outer Mongolian and Ewenki, especially Oroqen. Since genetic Only the last three digits are given, e.g. site 093 means 16093. ZQ13ZQ06ZQ18 – – – – T – – – C C – – – – – – – – C – C – – – – – C – – – – C – – – – – – – – T T – – – – – T – – A – – C + ZQ08ZQ12ZQ04 –ZQ11 –ZQ03 – CZQ14 – – C – – – – A – – A – – – – – – – – – – – – – – – – – – – – – T C – – – – – – – – – – – – – G – – – – – – – – – C – – – – – – – – – – – C – – T – – – T – – – C T – – – T T – – T – – – – – – – – – – – – C – C – C – – – C – – A C C A T C – – T – – – T – – – T – – T + – – + – – + – + + + + + + + + + + + + C C C C C ZQ05 – – – – – – – – – T – A – – T A A ZQ09 – – – A – – – – – – – – – – T – – ZQ10ZQ15 C C – – – – – – – – – – – – – – – C – – – C – – – – – – T T – – – – – T – – – – – C ZQ02 C – – – – – – – – – – – – – T – – – T – – – – – C ZQ16 C – – – – – – – – – – – – – T – – Table 1 Variable sites of HSV-I in 16Mitotypes Tuoba Xianbei individuals (position Polymorphic 16039–16398) nucleotide and positions restriction sites R TCTGTTCAACTGTCCGGTCTTGCGTCRS ZQ17 C – – – – – – – – – – – –informations – T – – in Tuoba Xianbei remains hold the key to study Y. Changchun et al. / FEBS Letters 580 (2006) 6242–6246 6245

Table 2 The haplogroup distribution frequencies (%) between Tuoba Xianbei and 11 populations compared Haplogroups Xib, Dau, Oro, Ewn, Kor, IM, OM, Bur, Yak, SH, NH, Uzb, N =16 N =45a N =44a N =47a N =48a N =48a N = 100b N = 126c N = 117c N =72d N = 101d N =58e D 43.75 24.4 43.1 31.9 33.3 39.6 21 19.9 28.9 11.1 29.7 22.4 C 31.25 6.7 29.5 19.1 6.3 14 15.9 38.5 2.8 1.0 1.7 A 6.25 2.2 4.5 4.3 14.6 8.3 4 7.1 1.7 9.7 5.0 B 12.5 15.6 2.3 10.6 10.5 10.5 10 3.6 19.4 13.9 other 6.25 54.1 20.6 34.1 41.6 35.3 51 54.1 30.9 57 50.4 85.9 Population Tuoba Xianbei, Daur, Oroqen, Ewenki, Koran, Inner Mongolian, Outer Mongolian, Buryat, Yakut, Southern Han, Northern Han and Uzbek are abbreviated as Xib, Dau, Oro, Ewn, Kor, IM, OM, Bur, Yak, SH, NH and Uzb, respectively. aData from Ref. [17]. bData from Ref. [15]. cData from Ref. [18]. dData from Ref. [22]. eData from Ref. [20].

Table 3 Population pairwise Fsts (below the diagonal) and the matrix of significant Fst P values (above the diagonal) Xib Eve Bur OM IM Dau Oro Ewn Kor NH SH Kaz Uzb Tur Xib 0.00000 + + + + + ÀÀ++++++ Eve 0.05937 0.00000 + + ++++++++++ Bur 0.05942 0.08885 0.00000 + À +++++++++ OM 0.01888 0.04989 0.02488 0.00000 + À +++++À ++ IM 0.03000 0.07758 0.00582 0.01025 0.00000 À +++À ++++ Dau 0.04690 0.07484 0.02606 0.00992 0.01044 0.00000 + À + À ++++ Oro 0.01744 0.03679 0.03291 0.02706 0.03138 0.03728 0.00000 + + + + + + + Ewn 0.01994 0.05896 0.03055 0.01204 0.02235 0.01210 0.03192 0.00000 + + + + + + Kor 0.04151 0.10005 0.03233 0.02282 0.01027 0.01622 0.03740 0.03187 0.00000 + + + + + NH 0.02968 0.08101 0.01983 0.01053 0.00288 0.00531 0.02950 0.01795 0.00755 0.00000 + + + + SH 0.05461 0.10773 0.06481 0.02893 0.02288 0.02066 0.05847 0.03565 0.02487 0.01131 0.00000 + + + Kaz 0.05139 0.09065 0.04140 0.00678 0.02098 0.01704 0.05159 0.01970 0.02715 0.01299 0.01770 0.00000 À + Uzb 0.08808 0.11408 0.04992 0.01618 0.03277 0.03178 0.08515 0.03265 0.05662 0.02554 0.03748 0.00244 0.00000 À Tur 0.14590 0.16614 0.11987 0.05694 0.09215 0.07753 0.14585 0.06970 0.11704 0.07197 0.07151 0.02870 0.01195 0.00000 Populations Tuoba Xianbei, Evenk, Buryat, Outer Mongolian, Inner Mongolian, Daur, Oroqen, Ewenki, Korean, Northern Han, Southern Han, Kazakh, Uzbek and Turkey are abbreviated as Xib, Eve, Bur, OM, IM, Dau, Oro, Ewn, Kor, NH, SH, Kaz, Uzb and Tur, respectively.

Turkey

Uzbek ZQ-Xianbei

Oroqen Kazakh Ewenki Outer-Mongolian Southern-Han

Northern-Han Inner-Mongolian

Daur Korean 0.0005

Fig. 2. Unrooted neighbor-joining tree of Tuoba Xianbei group and 13 populations compared. The ZQ-Xianbei indicates the Tuoba Xianbei.

the origin of the nomadic populations from northern China, Acknowledgements: This study was supported by the National Science further mtDNA analyses on remains from other Xianbei cem- Fund for Fostering Talents in Basic Research (Grant Number, eteries are necessary in order to have a complete understanding J0530184), Ministry of Education, China (the code, 00–07). We are grateful to Wei-jian who provided us the samples of Tuoba Xianbei on Tuoba Xianbei genetic structure. remains. 6246 Y. Changchun et al. / FEBS Letters 580 (2006) 6242–6246

2 Evn

Tur 1

Uzb Ewn Xib Kaz OM 0 Oro Dau Dimension 2 NH Bur IM -1 SH Kor

-4 -2 0 2 4 Dimension 1

Fig. 3. MDS plot of Tuoba Xianbei group and 13 populations compared. Populations Tuoba Xianbei, Evenk, Buryat, Outer Mongolian, Inner Mongolian, Daur, Oroqen, Ewenki, Korean, Northern Han, Southern Han, Kazakh, Uzbek and Turkey are abbreviated as Xib, Eve, Bur, OM, IM, Dau, Oro, Ewn, Kor, NH, SH, Kaz, Uzb and Tur, respectively.

References implications for the origin of new world founders. Genetics 142, 1321–1334. [1] Higuchi, R., Bowman, B. and Freiberger, M., et al. (1984) DNA [16] Kaessmann, H., Zollner, S. and Gustafsson, A.C., et al. (2002) sequences from the quagga, an extinct member of the horse Extensive linkage disequilibrium in small human populations in family. Nature 312, 282–284. . Am. J. Hum. Genet. 70, 673–685. [2] Paabo, S., Higuchi, R.G. and Wilson, A.C. (1989) Ancient DNA [17] Kong, Q.P., Yao, Y.G. and Liu, M., et al. (2003) Mitochondrial and the polymerase chain reaction: The emerging field of DNA sequence polymorphisms of five ethnic populations from molecular archaeology. J. Biol. Chem. 264, 9709–9712. northern China. Hum. Genet. 113, 391–405. [3] Cipollaro, M., Bernardo, G. and Galano, G., et al. (1998) [18] Pakendorf, B., Wiebe, V. and Tarskaia, L.A., et al. (2003) Ancient DNA in human bone remains from Pompeii archaeo- Mitochondrial DNA evidences for admixed origin of Central logical site. Biochem. Biophys. Res. Commun. 247, 901–904. Siberian populations. Am. J. Phys. Anthropol. 120, 211–224. [4] Colson, I.B., Richards, M.B. and Bailey, J.F., et al. (1997) DNA [19] Comas, D., Calafell, F. and Mateu, E., et al. (1998) Trading analysis of seven human skeletons excavated from the Terp of genes along the Silk Road: mtDNA sequences and the origin Wijnaldum. J. Archaeolog. Sci. 24, 911–917. of Central Asian populations. Am. J. Hum. Genet. 63, 1824– [5] Vigilant, L., Stoneking, M. and Harpending, H., et al. (1991) 1838. African populations and the evolution of mitochondrial DNA. [20] Yao, Y.G., Kong, Q.P. and , C.Y., et al. (2004) Different Science 253, 1503–1507. matrilineal contributions to genetic structure of ethnic groups in [6] Bertranpetit, J., Sala, J. and Calafell, F., et al. (1995) Human the silk road region in China. Mol. Biol. Evol. 21, 2265–2280. mitochondrial DNA variation and the origin of the Basques. Ann. [21] Comas, D., Calafell, F. and Mateu, E., et al. (1996) Geographic Hum. Genet. 59, 63–81. variation in human mitochondrial DNA control region sequence: [7] Seielstad, M.T., Minch, E. and Cavalli-Sforza, L.L. (1998) the population history of Turkey and its relationship to the Genetic evidence for a higher female migration rate in human. European populations. Mol. Biol. Evol. 13, 1067–1077. Nat Genet 20, 219–220. [22] Yao, Y.G., Kong, Q.P., Bandelt, H.J. and Kivisild, T. (2002) [8] Yao, Y.G., , Z.G. and Zhou, Z.D., et al. (2002) Genetic Phylogeographic differentiation of mitochondrial DNA in Han relationship of Chinese ethnic populations revealed by mtDNA Chinese. Am. J. Hum. Genet. 70, 635–651. sequence diversity. Am. J. Phys. Anthropol. 118, 63–76. [23] Oota, H., Kitano, T. and Jin, F., et al. (2002) Extreme mtDNA [9] Pakendorf, B., Wiebe, V. and Larskaia, A.T., et al. (2003) homogeneity in continental Asian populations. Am. J. Phys. Mitochondrial DNA evidences for admixed origin of Central Anthropol. 118, 146–153. Siberian populations. Am. J. Phys. Anthropol. 120, 211–224. [24] Anderson, S., Bankier, A.T. and Barrell, B.G., et al. (1981) [10] Wei, J. (2004) The Finding and Researching on the Xianbei Sequence and organization of the human mitochondrial genome. Cemeteries in Inner Mongolian Region, Science Press, . Nature 290, 457–465. [11] Gilbert, M.T., Bandelt, H.J., Hofreiter, M. and Barnes, I. (2005) [25] Ward, R.H., Frazier, B.L. and Dew-Jager, K., et al. (1991) Assessing ancient DNA studies. Trends Ecol. Evol. 20, 541–544. Extensive mitochondrial diversity within a single Amerindian [12] Kalmar, T., Bachrati, C.Z., Marcsik, A. and Rasko, I. (2000) A tribe. Proc. Natl. Acad. Sci. USA 88, 8720–8724. simple and efficient method for PCR amplifiable DNA extraction [26] Wallace, D.C., Brown, M.D. and Lott, M.T. (1999) Mitochon- from ancient bones. Nucleic Acids Res. 28, e67. drial DNA variation in human evolution and disease. Gene 238, [13] Adcock, G.J., Dennis, E.S. and Easteal, S., et al. (2001) 211–230. Mitochondrial DNA sequences in ancient Australians: Implica- [27] Kivisild, T., Tolk, H.V. and Parik, J., et al. (2002) The emerging tions for modern human origins. Proc. Natl. Acad. Sci. USA 98, limbs and twigs of the East Asian mtDNA tree. Mol. Biol. Evol. 537–542. 19, 1737–1751. [14] Batista, O., Kolman, C.J. and Bermingham, E. (1995) Mitochon- [28] Fu, L.Y. and Yang, Y. (1983) The Brief History on Northeast drial DNA diversity in the Kuna Amerinds of Panama. Hum. Ethnic Groups, Jilin People Press, Changchun. Mol. Genet. 4, 921–929. [29] , R. and Yip, V.F. (1993) Ethnic groups in China, Science [15] Kolman, C.J., Sanbunghin, N. and Bermingham, E. (1996) Press, Beijing. Mitochondrial DNA analysis of Mongolian populations and