Gene Flow at the Crossroads of Humanity: Mtdna Sequence Diversity and Alu Insertion Polymorphism Frequencies in Uzbekistan
Total Page:16
File Type:pdf, Size:1020Kb
The Open Genomics Journal, 2009, 2, 1-11 1 Open Access Gene Flow at the Crossroads of Humanity: mtDNA Sequence Diversity and Alu Insertion Polymorphism Frequencies in Uzbekistan Eric J. Devor*,1, Ibrokhim Abdurakhmonov2, Mark Zlojutro3, Meredith P. Millis1, Jessica J. Galbraith1, Michael H. Crawford4, Shukhrat Shermatov2, Zabardast Buriev2 and Abdusattor Abdukarimov2 1Molecular Genetics and Bioinformatics, Integrated DNA Technologies, Coralville, Iowa 52241 USA 2Laboratory of Genetic Engineering and Biotechnology, Institute of Genetics and Plant Experimental Biology, Academy of Sciences of Uzbekistan, Yuqori Yuz, Tashkent-702151, Uzbekistan 3Department of Genetics, Southwest Foundation for Biomedical Research, San Antonio, Texas 78227, USA 4Department of Anthropology, University of Kansas, Lawrence, Kansas 66045 USA Abstract: Mitochondrial DNA (mtDNA) HVS-I region sequences were obtained from 47 unrelated individuals represent- ing 10 of 12 viloyats of Uzbekistan. In addition, frequencies for five Alu insertion polymorphisms were determined for the same 47 individuals. These data were used to assess the genetic position of Uzbekistan relative to other Central Asians and, more generally, to Eurasian groups. Results show that the Uzbek sample has an approximate balance of west Eura- sian (53.2%) and Asian (46.8%) mtDNA lineages, which is reflected by their intermediate position relative to other Eura- sian groups in MDS plots based on genetic distance matrices. The Uzbeks also exhibit high sequence diversity, a frag- mented median-joining network, and a low regional ST score, all of which suggests a high degree of gene flow from neighboring gene pools. This is consistent with Central Asia’s history of repeated incursions by various nomadic peoples from the Asian steppe and the location of the well-traveled Silk Road within the present-day borders of Uzbekistan. In ad- dition, a star-like cluster within haplogroup H was identified in the Uzbek network, potentially representing a west Eura- sian genetic substratum of an earlier occupation by Caucasoid peoples. Keywords: Mitochondrial DNA, Alu polymorphism, genetic history, Uzbekistan. INTRODUCTION distance matrices that position the Uzbeks intermediate to western and eastern Eurasian populations, strengthening the Central Asia, which encompasses the territories between view of Uzbekistan as a crossroad of humanity. the Caspian Sea and the western boundaries of present-day China, has been a conduit for cultural, economic, and bio- SUBJECTS AND METHODS logical exchanges between Europe and Asia for thousands of years [1]. While much of Central Asia is composed of deso- Population Sample late high altitude regions and broad cold deserts devoid of A panel of 50 DNA samples was prepared from whole significant human habitation, its river basins and steppes blood draws at the Institute of Genetics and Plant Experi- have been occupied since the Lower Paleolithic [2]. How- mental Biology of the Academy of Sciences in Tashkent. ever, in spite of its historical importance during the rise of The origin of each sample was assigned by the birthplace of several far-flung empires and its long record of human habi- the donor’s mother. Three samples were eliminated because tation, the area is poorly characterized genetically. Several maternal birthplace was not located within Uzbekistan. The reports have begun to address this deficiency [3-6], revealing remaining 47 individuals are from ten of the twelve viloyats that the region has been at the center of widespread demo- or political districts of Uzbekistan, with the majority (33 of graphic processes responsible for molding the present-day 47 samples, 70.2%) representing the viloyats of Bukhara, genetic variation of local populations, as well as neighboring Kashkadarya, and Tashkent (Fig. 1). Thus, the panel pro- Eurasian groups. In this study, we contribute to the genetic vides a broader geographic coverage of Uzbekistan than understanding of Central Asia by reporting mitochondrial other samples that are in the literature. DNA sequence diversity and Alu element insertion frequen- mtDNA Sequencing and RFLP Analyses cies in Uzbekistan. Our results show comparable frequencies of European and East Asian mtDNA lineages which is re- A 447 bp mtDNA amplicon containing HVS-I segment flected by ordination analyses of the pairwise genetic of the mitochondrial control region was obtained for each member of the Uzbek panel using the primer set mtFOR: 5’- TCCACCATTAGCACCCAAAGCTA-3’ (pos. 15,976 to *Address correspondence to this author at the Molecular Genetics and Bio- o 15,998; TM= 58.4 C) and mtREV: 5’- ATTGATTTCACG- informatics, Integrated DNA Technologies, 1710 Commercial Park, GAGGATGG-3’ (pos. 16,422 to 16,402; T = 52.6oC). Am- Coralville, Iowa 52241, USA; Tel: 319-626-8450; M E-mail: [email protected] plicons were then sequenced in both directions on an Ap- 1875-693X/09 2009 Bentham Open 2 The Open Genomics Journal, 2009, Volume 2 Devor et al. Fig. (1). Map of Uzbekistan. Historic silk routes are shown as dashed lines along with principal regional sites. mtDNA samples were as- signed by mother’s birthplace to Kashkadarya (13), Tashkent (12), Bukhara (8), Surkhandarya (4), Namagan (3), Samarkand (2), Khorezm (2), Djizakh (1), Andijan (1), and Sirdarya (1). plied Biosystems Model 310 automated fluorescence se- Alu Insertion Polymorphisms quencer. Alu elements are a class of retroposons unique to primate Additional characterization of mtDNA haplogroup as- genomes [7]. There are more than one million Alu elements signments was obtained, as needed, by RFLP typing. The in the human genome [8] and a few are young enough that 10,397 Alu I and 10,394 Dde I cut sites were typed by re- they have yet to become fixed. Polymorphic Alu insertions striction enzyme digestion of a 269 bp amplicon produced by are two allele loci easily typed by PCR amplification. We the primers 10,290 FOR: 5’-GCCCTACAAACAACTAAC chose five loci for typing the Uzbek panel; FXIIIB, TPA25, o CTGCCA-3’ (TM= 59.2 C) and 10,558 REV: 5’- AGGGAG APO, PV92, and ACE. Summary information on each locus, o GATATGAGGTGTGAGCGAATA-3’ (TM= 60.1 C). The including PCR primer sequences and PCR reaction condi- 12,406 Hinc II/Hpa I cut site was typed by restriction en- tions, are presented in Table 1. zyme digestion of a 215 bp amplicon produced by primers Comparison Population Data 12,361 FOR: 5’-ACCGAGAAAGCTCACAAGAACTGC-3’ o (TM= 59.2 C) and 12,575 REV: 5’-AGCTTAAGGGAGAGC Comparative information for mtDNA sequences and Alu o TGGGTTGTTT-3’ (TM= 60.0 C). These primer sets were insertion polymorphism frequencies was assembled for a designed and analyzed for melting temperature and secon- number of populations both within and outside Central Asia. dary structure using the on-line SciTools software at Inte- MtDNA sequences were obtained for the Central Asian grated DNA Technologies (www.idtdna.com). PCR condi- populations Kazakhs, Kirghiz, and Uighurs [3]. Southern o 5:00 o 0:30 o 0:30 o 0:30 o 7:00 tions were 94 C [94 C ; 53 C ; 72 C ]3572 C for Siberian groups are represented by the Altai [9], Tuva and the HVR I amplifications and 94oC5:00[94oC0:30; 58oC0:30; Buryats [10], while Eastern Asia is represented by the Mon- o 0:30 o 7:00 72 C ]3572 C for the RFLP amplifications. Restriction gols [11] and Han Chinese [12]. European groups consist of enzyme digests of the RFLP amplicons were carried out as Russian [13], Greek [14], and Anatolian Turk [15] samples. per supplier recommendations. Uzbekistan Population Genetics The Open Genomics Journal, 2009, Volume 2 3 Table 1. Primer Sequences and PCR Conditions for the Alu Insertion Polymorphisms Alu ID Chr. 5’ Primer Sequence 3’ Tm(oC)a Ta(oC)b Allelesc TCAACTCCATGAGATTTTCAGAAGT 56.2 FXIIIB 1 54.0 500;200 CTGGAAAAAATGTATTCAGGTGAGT 55.3 GTAAGAGTTCCGTAACAGGACAGCT 59.8 TPA25 8 58.0 424;104 CCCCACCCTAGGAGAACTTCTCTTT 61.3 AAGTGCTGTAGGCCATTTAGATTAG 56.5 APO 11 55.0 409;83 AGTCTTCGATGACAGCGTATACAGA 58.9 AACTGGGAAAATTTGAAGAGAAAGT 55.3 PV92 16 55.0 442;119 TGAGTTCTCAACTCCTGTGTGTTAG 58.0 CTGGAGACCACTCCCATCCTTTCT 61.6 ACE 17 58.0 470;173 GATGTGGCCATCACATTCGTCAGAT 61.0 aSequence melting temperature estimated by OligoAnalyzer 3.0 (see SciTools at www.idtdna.com). b o 5:00 o 0:30 o 0:30 o 0:30 o 7:00 Temperature set for primer annealing step of the PCR assay. Constant conditions are: 94 C [94 C ; Ta( C) ; 72 C ]3572 C . cAmplicon length with and without the Alu insertion. FXIIIB, TPA25, APO, PV92, and ACE Alu insertion sional scaling (MDS) using the NTSYS statistical package frequencies used for population comparison were obtained [23]. from the ALFRED database at Yale University (http://alfred. For the Alu polymorphism frequencies, the coefficient of med.yale.edu/alfred/index.asp), including Uighur, Tajik, gene differentiation, G = (H – H )/H , was computed for Buryat, Turk and Greek populations. In addition, several ST T S T each locus using the executable DISPAN [24], where HT is Asian population samples were typed in the present study. the gene diversity for the total population (i.e., average allele These include the Yakuts (n = 50) and Mongols (n = 50), as frequencies for the entire data set) and H is the average of well as the Chuvash (n = 50) and Koryaks (n = 19) that have S the gene diversities computed for the individual not been previously typed for these particular markers. (sub)populations [25]. Nei’s modified Cavalli-Sforza genetic Polymorphism frequencies from these additional samples distance [26] was calculated between the populations and have been deposited in ALFRED. visualized using the MDS methodology. Data Analysis Results Using the software package Arlequin ver. 2.0 [16], three mtDNA Sequence Diversity different sequence diversity measures were employed: (i) 2 mtDNA sequences obtained from the Uzbek panel are haplotypic diversity [(n/n-1)(1-p )], where p is the fre- i i presented in Fig.