Search tips
Search criteria 


Logo of croatmedjFree full text at www.cmj.hrAboutSubscribeSubmitInfo for AuthorsFree full text at
Croat Med J. 2010 June; 51(3): 202–208.
PMCID: PMC2897085

Population Genetic Analysis of Haplotypes Based on 17 Short Tandem Repeat Loci on Y Chromosome in Population Sample from Eastern Croatia



To investigate the population genetics of 17 short tandem repeat (STR) loci on the Y chromosome in the population of eastern Croatia.


We carried out a statistical analysis of the data from previously performed genetic analysis collected during routine forensic work by the Forensic Science Centre “Ivan Vučetić.” A total of 220 unrelated healthy men from eastern Croatia were selected for the purpose of this study. Genomic DNA was extracted by Chelex from FTA® cards. Y-chromosomal STRs were determined using the AmpFISTR Yfiler PCR amplification kit. The haplotype frequencies were determined by direct counting and analyzed using Arlequin 3.1 and analysis of molecular variance calculated with the Y-chromosome haplotype reference database online analysis tool.


A total of 207 haplotypes were recorded, 197 of which were unique (90%). Haplotype diversity was 0.9993, with the most frequent haplotype found in 4 of 220 men (1.8%). Average locus diversity was 0.600, and it ranged from 0.256 for DYS392 to 0.780 for DYS458. Our results were compared with the pattern of Y-chromosome variability in publicly available population samples based on a minimal European haplotype set of 9 STRs and the greatest resemblance was found with samples from the Croatian capital of Zagreb, from Bosnia and Herzegovina, and from Serbia.


This is the first description of Y chromosome haplotyping of the population of eastern Croatia, which may serve as a basis for genetic epidemiology and forensic studies. Further studies are needed for characterization of the genetic structure of the Y-chromosome in the modern Croatian population.

The Y chromosome is much smaller than the X chromosome and is largely composed of repeating sequences (1). It makes up only about 2% of the total haploid genome and spans approximately 60 Mb. The non-recombining region (NRY) of the human Y chromosome comprises approximately 95% of the chromosome (1). Most of the genes present on the Y chromosome have their homologues on the X chromosome and these genes on the X chromosome are not subjected to X inactivation. The NRY is flanked on both sides by pseudoautosomal regions, where X-Y crossing-over is a normal and frequent event in male meiosis.

Human Y-chromosome short tandem repeats (Y-STR) are tandemly repeated arrays of 2-7 base pair units on the NRY. They are present only in men and are transmitted from father to son unchanged, as a haplotype, except for occasional mutations. Therefore, Y-STR haplotyping is of great importance for forensic applications, such as identification of unknown persons; sexual assault cases, where Y-STRs provide a male-specific DNA profile; paternity testing, even in cases when the alleged father is deceased; verification of amelogenin Y-deficient men; human migration and evolutionary studies; and historical and genealogical research (2). The frequency of individual haplotype is important in a relevant population, especially in forensic cases when the evidence and suspect match and the frequency of the Y-STR haplotype is needed to calculate a match probability (3). The increased use of human Y-STR markers in forensic analysis and in anthropological and archeological research has prompted a collaborative effort to collect haplotype data from different populations and to create the Y Chromosome Haplotype Reference Database (YHRD, (4).

For Croatia, the YHRD previously contained haplotypes only from the capital city of Zagreb. We have expanded the YHRD to include Y-STR data covering all Croatian regions, so we have recently added a total of 1100 haplotypes from eastern, western, southern, northern, and central Croatia (220 haplotypes for each region). As a first step in analyzing these new data, the present study investigated the population genetics of the 17 STR loci on the Y chromosome in the eastern Croatian population.

Participants and methods

Study sample

We carried out a statistical analysis of the data from previously performed genetic analysis collected during routine forensic work by the Forensic Science Centre “Ivan Vučetić.” A total of 220 samples from the following 5 counties in eastern Croatia were evaluated: Virovitičko-podravska, Požeško-slavonska, Brodsko-posavska, Osječko-baranjska, and Vukovarsko-srijemska County. The participants were not related and the samples were of sufficient quality and quantity to be included in the statistical analysis. The study was approved by the Ethics Committee of the Institute for Medical Research and Occupational Health, Zagreb, Croatia.

DNA analysis

Genomic DNA from all samples of the materials expertise was extracted from FTA cards (Whatman, Maidstone, Kent, UK) using Chelex (5). After isolation, the genomic DNA content of each sample was determined by quantitative real-time polymerase chain reaction (PCR) using the Quantifiler™ Human Male DNA Quantification kit (Applied Biosystems, Foster City, CA, USA), which includes an internal positive control to test for the presence of PCR inhibitors in the DNA extracts. Quantitative real-time PCR was performed on a 7500 Real-Time PCR System (Applied Biosystems). Genomic DNA (1 ng) was amplified using the AmpFISTR Yfiler PCR amplification kit (Applied Biosystems), which coamplifies 17 Y-STRs: DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS385, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635, and GATAH4. These amplification reactions were performed using a GeneAmp PCR System 9700 (Applied Biosystems). Y-STR amplification products were analyzed on 3130xl Genetic Analyzer (Applied Biosystems). Analysis of the data was performed using Genemapper® software (version 3.2, Applied Biosystems). Amplicon sizing was performed using an internal size standard (GeneScan-500 LIZ, Applied Biosystems), and the amplicons were compared with the AmpFlSTR Yfiler allelic ladder for unambiguous allele designation.

Statistical analysis

Allelic and haplotype frequencies were estimated by direct counting. Locus and haplotype diversities were calculated using Arlequin 3.1 (6), according to the formula:

graphic position="anchor" xlink:href="CroatMedJ_51_0202-M1" xmlns:xlink=""/>

where n is the population size and pi is the frequency of the i-th haplotype. Locus diversity was calculated according to the same formula, using allele frequencies instead of haplotype frequencies.

Furthermore, we used the YHRD ( to characterize the Y-STRs in the samples and compare them with the countries in the region. Due to the limited number of markers reported in other publicly available samples, the comparative analysis was performed on a minimal European Y-STR haplotype comprising 9 loci: DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, and DYS385ab (7). A total of 10 population samples were included in this analysis, to give a total of 1469 samples. The samples were from the following areas (the YHRD designation is given in parentheses): Budapest, Hungary (“Hungarian,” n = 200 haplotypes); Doboj-Banja Luka-Bijeljina, Bosnia and Herzegovina (“Bosnian,” n = 31); Ljubljana, Slovenia (“Slovenian,” n = 180); Mostar, Bosnia and Herzegovina (“Bosnian,” n = 34); Novi Sad, Serbia (“Serbian,” n = 215); Romania (“Romanian,” n = 104); Sarajevo, Bosnia and Herzegovina (“Bosnian,” n = 35); Szeged, Hungary (“Hungarian,” n = 100); Zagreb, Croatia (“Croatian,” n = 150); and the study sample from eastern Croatia, Croatia (“Croatian,” n = 220). Since analyses conducted in the YHRD are limited to 10 populations per analysis, we also performed a second-stage analysis, with a number of other more geographically and historically distant populations. These comparisons were also made using the minimal haplotype set of 9 loci and the publicly available population sets in the YHRD. This analysis yielded significant differences from the eastern Croatian sample: Albania (“Albanian,” n = 111), Republic of Macedonia (“Macedonian,” n = 250); Graz, Austria (“Austrian,” n = 65); Bonn, Germany (“German,” n = 90); and the Czech Republic (“Czech,” n = 69).

We calculated population pairwise genetic distances (Rst). This is an extension of the commonly used Fst measure (8), defined as Rst = (Sb - Sw)/Sb, where Sw is the sum over all loci of twice the weighted mean of the within-population variances V(A) and V(B), and Sb is the sum over all loci of twice the variance V(A+B) of the combined population (8). P values were calculated using analysis of molecular variance (AMOVA), with 10 000 permutations, using an online tool of the YHRD. Significance was set at P < 0.05.


The initial analysis indicated a wide extent of variance of the number of alleles in different loci, with the most recorded alleles in DYS458 (Table 1). Locus diversity for the entire sample was 0.600, ranging from 0.256 for DYS392 to 0.780 for DYS458. A total of 207 haplotypes were recorded, 197 of which were unique (90%). Total haplotype diversity was 0.9993, with the commonest recorded haplotype found in 4 of 220 men (1.8%) (YHRD accession code YA003594). Discrimination capacity was 94.1%.

Table 1
Allele frequencies and locus diversity values at Y chromosome short tandem repeat (Y-STR) loci in the eastern Croatian population*

AMOVA indicated that our eastern Croatian sample was closest to populations from the Croatian capital of Zagreb, Bosnia and Herzegovina, and Serbia. The remaining pair-wise comparisons showed significant differences (Table 2). A second-stage analysis yielded significant differences between our eastern Croatian sample and the following samples (P < 0.001 in all cases): Albania (“Albanian,” Rst = 0.052), Republic of Macedonia (“Macedonian,” Rst = 0.034), Graz, Austria (”Austrian,“ Rst=0.102), Bonn, Germany (“German,” Rst = 0.101), and the Czech Republic (“Czech,” Rst = 0.144).

Table 2
Analysis of molecular variance pair-wise distances based on Rst values between the eastern Croatian population in this study and selected comparison populations


This study is the first attempt to characterize the genetic structure of the eastern Croatian population in terms of Y-chromosome STRs. Such information is of crucial importance for any forensic work, for genetic epidemiology and population genetic purposes.

Nowadays, a great effort is made in order to improve the discrimination power of Y chromosome haplotypes as an increased number of Y-STRs is used and introduced in commercially available typing kits. The Y-STR analysis has been facilitated when commercial kits, such as the AmpFlSTR Yfiler, has become commercially available (9-11,20).

We have recently expanded the YHRD to include Y-STR data from all 5 regions of Croatia (western, northern, southern, central, and eastern). Our results from the eastern region of Croatia can be directly compared with those published for Splitsko-dalmatinska County in the southern region (12), because this study examined the same Y-STR polymorphic loci using the same AmpFlSTR Yfiler amplification kit. These two studies showed the same pattern of most frequent alleles across all loci, except DYS391 and DYS385 (12). The different results at locus DYS385 are consistent with the fact that DYS385 is one of the most polymorphic Y-STR markers.

A study in western Croatia (Primorsko-goranska County) analyzed fewer loci than our study (13). Nevertheless, they found the same pattern of the most frequent alleles for DYS389I, DYS390, DYS19, DYS393, DYS391, and DYS392 loci, but not for DYS389II and DYS385. Together with the results from eastern and southern Croatia, this study shows that the DYS385 locus diversity values differ between the three regions of Croatia.

DYS458 locus is a heterogeneous locus composed of a polymorphic number of tetranucleotide repeats (GAAA) (14). A large deletion in the Yp11.2 band of at least 1.13 Mb was reported to include amelogenin, the NRY1 minisatellite, and the DYS458 locus (15). Intermediate alleles, which contain incomplete repeat units, have been described at the DYS458 locus (16) and some of them are denoted with a “0.2.” Such intermediate alleles arise by insertion/deletion events most probably caused by slipped-strand mis-pairing within the locus during spermatogenesis. These intermediate alleles were shown to be most frequent in north and east Africa, and in the Caucasus, while their occurrence in Europe is less widespread (16). This study found for the first time the variant alleles 12.2, 16.2, 17.2, and 19.2 at the DYS458 locus in the Croatian population. The 12.2 allelic variant at the DYS458 locus has not been previously reported in the YHRD (4). The 16.2 allelic variant at the DYS458 locus has been previously described in Germany, Poland, Portugal, and Russia (4,17,18). The 19.2 allelic variant at DYS458 locus has been found in the Republic of Macedonia, Poland, Russia, Greece, Italy, and Austria (4,17-19,21,22), while the 17.2 variant has been reported in Greece, Italy, Poland, Portugal, Germany, Albania, Republic of Macedonia, and Russia (4,17-19,21,23).

Our study showed a great similarity between our and the samples from the Croatian capital of Zagreb, Bosnia and Herzegovina, and Serbia, while there were difference from other neighboring countries (24-30). These results are in line with expectations, suggesting that neighboring populations shared similar Y-STR origins. Populations that were geographically distant from the eastern Croatia sample were also more genetically distant, while the second-stage analysis showed that even more geographically distant populations were considerably genetically different (31-33).

The limitations of this study include the possibility that the sample was not ethnically representative of the population of eastern Croatia. Furthermore, the comparison with the other populations was based on the samples available in the YHRD, which were occasionally rather small and thus perhaps not representative of their corresponding populations. In light of this, further studies are needed to characterize the genetic structure of the entire Croatian population in order to provide the basis for future research in forensics, genetic epidemiology, and understanding of some major demographic and historical events in this region (34).


This study was funded by the Ministry of Interior of Croatia. The authors thank Sunčica Papić and Marija Valić for technical support.


1. Charlesworth B, Charlesworth D. The degeneration of Y chromosomes. Philos Trans R Soc Lond B Biol Sci. 2000;355:1563–72. doi: 10.1098/rstb.2000.0717. [PMC free article] [PubMed] [Cross Ref]
2. Butler JM. Recent developments in Y-short tandem repeat and Y-single nucleotide polymorphism analysis. Forensic Science Review. 2003;15:91–111. [PubMed]
3. Gusmao L, Butler JM, Carracedo A, Gill P, Kayser M, Mayr WR, et al. DNA Commission of the International Society of Forensic Genetics (ISFG): an update of the recommendations on the use of Y-STRs in forensic analysis. Forensic Sci Int. 2006;157:187–97. doi: 10.1016/j.forsciint.2005.04.002. [PubMed] [Cross Ref]
4. Roewer L, Krawczak M, Willuweit S, Nagy M, Alves C, Amorim A, et al. Online reference database of European Y-chromosomal short tandem repeat (STR) haplotypes. Forensic Sci Int. 2001;118:106–13. doi: 10.1016/S0379-0738(00)00478-3. [PubMed] [Cross Ref]
5. Walsh PS, Metzger DA, Higuchi R. Chelex 100 as a medium for simple extraction of DNA for PCR-based typing from forensic material. Biotechniques. 1991;10:506–13. [PubMed]
6. Excoffier L, Laval G, Schneider S. Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinform Online. 2005;1:47–50. [PMC free article] [PubMed]
7. Kayser M, Caglia A, Corach D, Fretwell N, Gehrig C, Graziosi G, et al. Evaluation of Y-chromosomal STRs: a multicenter study Int J Legal Med 1997. 11012533,141910.1007/s004140050051 [PubMed] [Cross Ref]
8. Slatkin M. A measure of population subdivision based on microsatellite allele frequencies. Genetics. 1995;139:457–62. [PubMed]
9. Mulero JJ, Chang CW, Calandro LM, Green RL, Li Y, Johnson CL, et al. Development and validation of the AmpFlSTR Yfiler PCR amplification kit: a male specific, single amplification 17 Y-STR multiplex system. J Forensic Sci. 2006;51:64–75. doi: 10.1111/j.1556-4029.2005.00016.x. [PubMed] [Cross Ref]
10. Kayser M, Roewer L, Hedman M, Henke L, Henke J, Brauer S, et al. Characteristics and frequency of germline mutations at microsatellite loci from the human Y chromosome, as revealed by direct observation in father/son pairs. Am J Hum Genet. 2000;66:1580–8. doi: 10.1086/302905. [PubMed] [Cross Ref]
11. Kayser M, Sajantila A. Mutations at Y-STR loci: implications for paternity testing and forensic analysis. Forensic Sci Int. 2001;118:116–21. doi: 10.1016/S0379-0738(00)00480-1. [PubMed] [Cross Ref]
12. Ljubkovic J, Stipisic A, Sutlovic D, Definis-Gojanovic M, Bucan K, Andelinovic S. Y-chromosomal short tandem repeat haplotypes in southern Croatian male population defined by 17 loci. Croat Med J. 2008;49:201–6. doi: 10.3325/cmj.2008.2.201. [PMC free article] [PubMed] [Cross Ref]
13. Lovrecic L, Ristic S, Brajenovic B, Kapovic M, Peterlin B. Human Y-specific STR haplotypes in the Western Croatian population sample. Forensic Sci Int. 2005;149:257–61. doi: 10.1016/j.forsciint.2004.06.026. [PubMed] [Cross Ref]
14. Myres NM, Ekins JE, Lin AA, Cavalli-Sforza LL, Woodward SR, Underhill PA. Y-chromosome short tandem repeat DYS458.2 non-consensus alleles occur independently in both binary haplogroups J1-M267 and R1b3-M405. Croat Med J. 2007;48:450–9. doi: 10.1016/S0379-0738(02)00347-X. [PMC free article] [PubMed] [Cross Ref]
15. Chang YM, Perumal R, Keat PY, Yong RY, Kuehn DL, Burgoyne L. A distinct Y-STR haplotype for Amelogenin negative males characterized by a large Y(p)11.2 (DYS458-MSY1-AMEL-Y) deletion. Forensic Sci Int. 2007;166:115–20. doi: 10.1016/j.forsciint.2006.04.013. [PubMed] [Cross Ref]
16. Ferri G, Robino C, Alu M, Luiselli D, Tofanelli S, Caciagli L, et al. Molecular characterisation and population genetics of the DYS458. 2 allelic variant. Forensic Science International: Genetics Supplement Series. 2008;1:203–5. doi: 10.1016/j.fsigss.2007.10.217. [Cross Ref]
17. Wolanska-Nowak P, Branicki W, Parys-Proszek A, Kupiec T. A population data for 17 Y-chromosome STR loci in South Poland population sample–some DYS458.2 variants uncovered and sequenced. Forensic Sci Int; Genet. 2009;4:e43–4. doi: 10.1016/j.fsigen.2009.04.009. [PubMed] [Cross Ref]
18. Roewer L, Willuweit S, Krüger C, Nagy M, Rychkov S, Morozowa I, et al. Analysis of Y chromosome STR haplotypes in the European part of Russia reveals high diversities but non-significant genetic distances between populations. Int J Legal Med. 2008;122:219–23. doi: 10.1007/s00414-007-0222-2. [PMC free article] [PubMed] [Cross Ref]
19. Kovatsi L, Saunier JL, Irwin JA. Population genetics of Y-chromosome STRs in a population of Northern Greeks. Forensic Sci Int; Genet. 2009;4:e21–2. doi: 10.1016/j.fsigen.2009.01.001. [PubMed] [Cross Ref]
20. Sanchez-Diz P, Alves C, Carvalho E, Carvalho M, Espinheira R, Garcia O, et al. Population and segregation data on 17 Y-STRs: results of a GEP-ISFG collaborative study. Int J Legal Med. 2008;122:529–33. doi: 10.1007/s00414-008-0265-z. [PubMed] [Cross Ref]
21. Onofri V, Alessandrini F, Turchi C, Fraternale B, Buscemi L, Pesaresi M, et al. Y-chromosome genetic structure in sub-Apennine populations of Central Italy by SNP and STR analysis. Int J Legal Med. 2007;121:234–7. doi: 10.1007/s00414-007-0153-y. [PubMed] [Cross Ref]
22. Berger B, Lindinger A, Niederstatter H, Grubwieser P, Parson W. Y-STR typing of an Austrian population sample using a 17-loci multiplex PCR assay. Int J Legal Med. 2005;119:241–6. doi: 10.1007/s00414-005-0546-8. [PubMed] [Cross Ref]
23. Ferri G, Ceccardi S, Lugaresi F, Bini C, Ingravallo F, Cicognani A, et al. Male haplotypes and haplogroups differences between urban (Rimini) and rural area (Valmarecchia) in Romagna region (North Italy). Forensic Sci Int. 2008;175:250–5. doi: 10.1016/j.forsciint.2007.06.007. [PubMed] [Cross Ref]
24. Marjanovic D, Fornarino S, Montagna S, Primorac D, Hadziselimovic R, Vidovic S, et al. The peopling of modern Bosnia-Herzegovina: Y-chromosome haplogroups in the three main ethnic groups. Ann Hum Genet. 2005;69:757–63. doi: 10.1111/j.1529-8817.2005.00190.x. [PubMed] [Cross Ref]
25. Veselinovic IS, Zgonjanin DM, Maletin MP, Stojkovic O, Djurendic-Brenesel M, Vukovic RM, et al. Allele frequencies and population data for 17 Y-chromosome STR loci in a Serbian population sample from Vojvodina province. Forensic Sci Int. 2008;176:e23–8. doi: 10.1016/j.forsciint.2007.04.003. [PubMed] [Cross Ref]
26. Roewer L, Croucher PJ, Willuweit S, Lu TT, Kayser M, Lessig R, et al. Signature of recent historical events in the European Y-chromosomal STR haplotype distribution. Hum Genet. 2005;116:279–91. doi: 10.1007/s00439-004-1201-z. [PubMed] [Cross Ref]
27. Furedi S, Woller J, Padar Z, Angyal M. Y-STR haplotyping in two Hungarian populations. Int J Legal Med. 1999;113:38–42. doi: 10.1007/s004140050276. [PubMed] [Cross Ref]
28. Barbarii LE, Rolf B, Dermengiu D. Y-chromosomal STR haplotypes in a Romanian population sample. Int J Legal Med. 2003;117:312–5. doi: 10.1007/s00414-003-0397-0. [PubMed] [Cross Ref]
29. Beer Z, Csete K, Varga T. Y-chromosome STR haplotype in Szekely population. Forensic Sci Int. 2004;139:155–8. doi: 10.1016/j.forsciint.2003.10.010. [PubMed] [Cross Ref]
30. Sterlinko H, Pajnic IZ, Balazic J, Komel R. Human Y-specific STR haplotypes in a Slovenian population sample. Forensic Sci Int. 2001;120:226–8. doi: 10.1016/S0379-0738(01)00390-5. [PubMed] [Cross Ref]
31. Spiroski M, Arsov T, Krüger C, Willuweit S, Roewer L. Y-chromosomal STR haplotypes in Macedonian population samples. Forensic Sci Int. 2005;148:69–73. doi: 10.1016/j.forsciint.2004.04.067. [PubMed] [Cross Ref]
32. Junge A, Madea B. Population studies of the Y-chromosome specific polymorphisms DYS19, DYS389 I + II, DYS390 and DYS393 in a western German population (Bonn area). Forensic Sci Int. 1999;101:195–201. doi: 10.1016/S0379-0738(99)00029-8. [PubMed] [Cross Ref]
33. Robino C, Gino S, Ricci U, Grignani P, Previdere C, Torre C. Y-chromosomal STR haplotypes in an Albanian population sample. Forensic Sci Int. 2002;129:128–30. doi: 10.1016/S0379-0738(02)00224-4. [PubMed] [Cross Ref]
34. Gršković B, Mršić G. Y chromosome: from evolution to forensics- an overwiev. Acta Med Croatica. 2010;64:33–40. [PubMed]

Articles from Croatian Medical Journal are provided here courtesy of Medicinska Naklada