Mitochondrial DNA control region analysis of three ethnic groups in the Republic of Macedonia (original) (raw)

Abstract

A total of 444 individuals representing three ethnic groups (Albanians, Turks and Romanies) in the Republic of Macedonia were sequenced in the mitochondrial control region. The mtDNA haplogroup composition differed between the three groups. Our results showed relatively high frequencies of haplogroup H12 in Albanians (8.8%) and less in Turks (3.3%), while haplogroups M5a1 and H7a1a were dominant in Romanies (13.7% and 10.3%, respectively) but rare in the former two. This highlights the importance of regional sampling for forensic mtDNA databasing purposes. These population data will be available on EMPOP under accession numbers EMP00644 (Albanians), EMP00645 (Romanies) and EMP00646 (Turks).

Keywords: mtDNA, EMPOP, Haplogroups, Albanians, Turks, Romanies

1. Introduction

According to demographic data from the register in 2002 the population of the Republic of Macedonia is 2,022,547 with 64.2% ethnical Macedonian, 25.2% ethnical Albanians, 3.9% ethnical Turks, 2.7% ethnical Romanies and a small percentage of other ethnic groups. People from different ethnic communities rarely have marriages between each other due to their national and religious determination. Whether or not this has an effect on the distribution of mitochondrial lineages has yet not been studied for Macedonia. An earlier study described mitochondrial (mt)DNA control region variation for ethnical Macedonians, which brought a similar haplogroup distribution to other West-Eurasian populations [1]. Here, we describe mtDNA control region variation in carefully selected samples of the three other major ethnic groups (148 Albanians, 150 Turks and 146 Romanies) and thus add a total of 444 high quality mtDNA lineages to the body of world-wide mtDNA database. The data will also be made available for forensic searches via EMPOP [2] under accession numbers EMP00644 (Albanians), EMP00645 (Romanies), and EMP00646 (Turks).

2. Materials and methods

2.1. Samples and DNA extraction

This study was reviewed and approved by the ethics commission of the University “St.Cyril and Methodius” (study classification number 03-5904/2 from 01.02.13, session number XXVI). All participants (N = 444) gave their written consent before a buccal swab was taken. Study participants were sampled from different geographic locations in the Republic of Macedonia (Fig. S1), Albanians derived mostly from the western part, Turks originated from the eastern, western and southern parts and Romanies derived from the central and northern parts of the country. DNA was extracted using the QIAamp mini kit (Qiagen, Hilden, Germany) following the manufacturer's recommendations.

2.2. Amplification, sequencing and data analysis

PCR amplification and mtDNA control region sequencing were carried out following the EMPOP protocol [3] updated in [4]. Nucleotide sequences were analysed and interpreted using Sequencher (Version 5.1, Gene Codes Corporation) and aligned relative to the rCRS [5] following phylogenetic alignment rules defined in [6]. The mtDNA haplotypes were determined using EMMA [7] based on Phylotree (www.phylotree.org; Build 16; [8]). The random match probability was calculated as sum of squares of the haplotype frequencies [9]. Genetic diversity indices were calculated using the ARLEQUIN software (Version 3.5) [10]. C-Stretch length variants in HVS-I (around 16,193), HVS-II (around 309) and HVS-III (around 573) were ignored for calculating random match probabilities and genetic diversity indices.

3. Results and discussion

The mtDNA control region sequence analysis in three Macedonian ethnic groups consisting of 444 individuals (148 Albanians, 150 Turks and 146 Romanies) showed 108 different haplotypes (73%) in Albanians, 100 (66.7%) in Turks and 64 (43.8%) in Romanies, respectively (Tables 1 and S1). Thereof, 87 (80.6%), 74 (74%) and 42 (65.6%) were unique and haplotype diversity was 0.983, 0.986 and 0.966 respectively (Table 1). AMOVA was performed taking into consideration the following published datasets: Macedonia [1], Greece [11], Cyprus [11], Hungarian Ashkenazi [12], Hungarian Baranya Romany [13], Hungarians from Budapest [13], Romanian Csango [14] and Romanian Szekely [14]. Fst comparison, pairwise differences and shared haplotypes are given in ESM 1.

Table 1.

Diversity measures for control region haplotypes (16024-576) from three Macedonian ethnic groups.

The distribution of observed lineages differed between the three investigated populations (Table 2). Albanians showed a relatively high abundance of hg H12 lineages (8.8%) that were generally rare elsewhere, 1.3% in northern Greeks [11] and 3% in Orthodox Macedonians [1]. Romanies showed high frequencies of hgs H7a1a (10.3%) and M5a1 (13.7%) that is common in the South Asian phylogeny [15]. This emphasizes the requirement of regional databases when assessing haplotype frequencies in a forensic context.

Table 2.

Distribution of selected haplogroups in three Macedonian ethnic groups.

Acknowledgements

The authors would like to thank all volunteers that participated in this study. This work leading to these results has received funding from the European Union Seventh Framework Programme (FP7/2007-2013) under grant agreement n° 285487 (EUROFORGEN-NoE) and was in part supported by the Austrian Science Fund (FWF) [P22880-B12]. Also, we would like to thank colleagues from Macedonia, especially to d-r Agim Ramadani and Sefedin Biljali for their help during samples collection.

Footnotes

Appendix A. Supplementary data

The following are supplementary data to this article:

References

• 1.Zimmermann B., Brandstätter A., Duftner N., Niederwieser D., Spiroski M., Arsov T., Parson W. Mitochondrial DNA control region population data from Macedonia. Forensic Sci. Int. Genet. 2007;1:e4–e9. doi: 10.1016/j.fsigen.2007.03.002. [DOI] [PubMed] [Google Scholar]
• 2.Parson W., Dür A. EMPOP—a forensic mtDNA database. Forensic Sci. Int. Genet. 2007;1:88–92. doi: 10.1016/j.fsigen.2007.01.018. [DOI] [PubMed] [Google Scholar]
• 3.Brandstätter A., Salas A., Niederstätter H., Gassner C., Carracedo A., Parson W. Dissection of mitochondrial superhaplogroup H using coding region SNPs. Electrophoresis. 2006;27:2541–2550. doi: 10.1002/elps.200500772. [DOI] [PubMed] [Google Scholar]
• 4.Parson W., Bandelt H. Extended guidelines for mtDNA typing population data in forensic science. Forensic Sci. Int. Genet. 2007;1:13–19. doi: 10.1016/j.fsigen.2006.11.003. [DOI] [PubMed] [Google Scholar]
• 5.Andrews R.M., Kubacka I., Chinnery P.F., Lightowlers R.N., Turnbull D.M., Howell N. Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat. Genet. 1999;23:147. doi: 10.1038/13779. [DOI] [PubMed] [Google Scholar]
• 6.Bandelt H., Parson W. Consistent treatment of length variants in the human mtDNA control region: a reappraisal. Int. J. Legal Med. 2008;122:11–21. doi: 10.1007/s00414-006-0151-5. [DOI] [PubMed] [Google Scholar]
• 7.Röck A., Dür A., vanOven M., Parson W. Concept for estimating mitochondrial DNA haplogroups using a maximum likelihood approach (EMMA) Forensic Sci. Int. Genet. 2013;7:601–609. doi: 10.1016/j.fsigen.2013.07.005. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 8.van Oven M., Kayser M. Updated comprehensive phylogenetic tree of global human mitochondrial DNA variation. Hum. Mut. 2008;30:e386–e394. doi: 10.1002/humu.20921. [DOI] [PubMed] [Google Scholar]
• 9.Stoneking M., Hedgecock D., Higuchi R.G., Vigilant L., Erlich H.A. Population variation of human mtDNA control region sequences detected by enzymatic amplification and sequence-specific oligonucleotide probes. Am. J. Hum. Genet. 1991;48:370–382. [PMC free article] [PubMed] [Google Scholar]
• 10.Schneider S., Roessli D., Excoffier L. Genetics and Biometry Laboratory, University of Geneva; Switzerland: 2000. Arlequin Ver. 2.0: A Software for Population Genetic Data Analysis. [Google Scholar]
• 11.Irwin J., Saunier J., Strouss K., Paintner C., Diegoli T., Sturk K., Kovatsi L., Brandstätter A., Cariolou M.A., Parson W., Parsons T.J. Mitochondrial control region sequences from northern Greece and Greek Cypriots. Int. J. Legal Med. 2008;122:87–89. doi: 10.1007/s00414-007-0173-7. [DOI] [PubMed] [Google Scholar]
• 12.Brandstätter A., Egyed B., Zimmermann B., Tordai A., Padar Z., Parson W. Mitochondrial DNA control region variation in Ashkenazi Jews from Hungary. Forensic Sci. Int. Genet. 2008;2:e4–e6. doi: 10.1016/j.fsigen.2007.07.006. [DOI] [PubMed] [Google Scholar]
• 13.Irwin J., Egyed B., Saunier J., Szamosi G., O’Callaghan J., Padar Z., Parsons T.J. Hungarian mtDNA population databases from Budapest and the Baranya country Roma. Int. J. Legal Med. 2007;121:377–383. doi: 10.1007/s00414-006-0128-4. [DOI] [PubMed] [Google Scholar]
• 14.Brandstätter A., Egyed B., Zimmermann B., Duftner N., Padar Z., Parson W. Migration rates and genetic structure of two Hungarian Ethnic groups in Transylvania, Romania. Ann. Hum. Genet. 2007;71:1–13. doi: 10.1111/j.1469-1809.2007.00371.x. [DOI] [PubMed] [Google Scholar]
• 15.Mendizabal I., Valente C., Gusmao A., Alves C., Gomes V., Goios A., Parson W., Calafell F., Alvarez L., Amorium A., Gusmao L., Comas D., Prata M. Reconstructing the Indian Origin and Dispersal of the European Roma: a maternal genetic perspective. PLoS ONE. 2011;6:1–10. doi: 10.1371/journal.pone.0015988. [DOI] [PMC free article] [PubMed] [Google Scholar]

Associated Data

This section collects any data citations, data availability statements, or supplementary materials included in this article.

Supplementary Materials