Analysis of Y chromosome STR haplotypes in the European part of Russia reveals high diversities but non-significant genetic distances between populations.

A total of 17 Y-specific STR loci were studied in 12 districts of the European part of Russia aiming to ascertain the amount of substructure required for the construction of a representative regional database. All groups exhibited high haplotype diversities but low inter-population variance as measured by an analysis of molecular variance. However, when Western Russia is taken as a whole, the genetic distances to the neighbouring populations were significant. Whereas gradual change in the Y chromosome pool exists between Russia and the Slavic-speaking populations to the West, remarkable discontinuities were observed with neighbouring populations in the East, North and South.

Introduction

Y chromosome short tandem repeat (STR) haplotypes show a substantial inter-population differentiation both on a worldwide and continental scale [1, 2]. Because Y-STR haplotype frequencies are required to provide statistical estimates of the significance of a match between forensic samples, local databases must be developed taking the demographic history of the investigated population into account. Large metadatabases as the Y Chromosome Haplotype Reference Database (YHRD) pool different local databases based on intrinsic (genetic) and contextual (geographic, linguistic) information [3]. However, assembling different regional pools is only valid if there is no population substructure, i.e. no statistically significant difference between the Y-STR haplotype distributions in different regions [4]. Here, we present haplotype analysis and genetic differentiation tests of a large Russian population consisting of 12 subgroups leading to a sensible decision on the assignment of the analysed population to a framework of worldwide metapopulations.

Population

A total of 545 unrelated males from 12 Western Russian populations previously typed for Y-chromosome single nucleotide polymorphism markers [5] were analysed for 17 Y-STR loci evaluated in forensic routine diagnostics [6]. The sampling was carried out in the administrative centers of the following districts (oblasts): Smolenskaja (Smo, n = 43), Brianskaja (Bri, 43), Ivanovskaja (Iva, 40), Lipezkaja (Lip, 47), Penzenskaja (Pen, 81), Ryazanskaja (Rya, 36), Orlovskaja (Orl, 42), Tverskaja (Tve, 43), Vologodskaja (Vol, 40), Tambovskaja (Tam, 48), Archangelskaja (Arch, 42) and Nowgorodskaja (Now, 40; see Fig. 1). Informed consent and information about the birthplace of the donor and his parents and grandparents were obtained.

Fig. 1

Map of Europe depicting the analysed populations from the European part of Russia and neighbouring territories

Materials and methods

Deoxyribonucleic acid (DNA) was extracted from whole blood samples using the QIAmp DNA Blood Mini Kit (Qiagen, Hilden, Germany) following the manufacturer’s recommendations. All samples were genotyped for 17 Y chromosomal STR loci (DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS385a, DYS385b, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635, GATA H4) using the AmpFlSTR®YFiler® PCR Amplification kit (Applied Biosystems, Darmstadt, Germany) according to the manufacturer’s instructions. STR products were analysed on an ABI Prism 3100 AVANT automated sequencer with GeneScan and GenoTyper v. 3.7 (Applied Biosystems). The updated recommendations of the DNA Commission of the International Society of Forensic Genetics for analysis of Y-STR systems were followed [4].

Quality control

Proficiency testing of the German DNA Profiling group (www.gednap.de) and the YHRD (www.yhrd.org) trials was carried out.

Analysis of data

Characteristic parameters for each population, consisting of the number of different haplotypes, the discrimination capacity (D) and the haplotype diversity (h), were calculated (Table 1). Pairwise values of Φ, an analogue of Wrights Fst that takes the evolutionary distance between individual haplotypes into account [7, 8], were calculated to measure genetic distances between 17 locus haplotypes of 12 Western Russian populations. Subsequently, these populations were treated as one metapopulation and the minimal haplotypes (which includes the nine loci: DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS385ab) were compared to published data from nine neighbouring regions, namely 271 minimal haplotypes from Belarus [9], 502 from the Caucasus [10], 133 from Estonia, 145 from Latvia and 157 from Lithuania [11], 399 samples from Finland [12], 3,021 from Poland [13-15], 370 from Siberia [16, 17] and 368 from Ukraine [18], with the statistical significance determined by a permutation test (10,000 replicates; Table 2). We used our own implementation of analysis of molecular variance (AMOVA; available at http://rprojekt.org/amova/). The DYS389I allele length was obtained by subtracting the shorter allele from the longer allele at DYS389I/II. The statistics and analytics software package (STATISTICA package; StatSoft Inc.) was used for multi-dimensional scaling (MDS) analysis [19] based on pairwise Φ values (Fig. 2).

Table 1

Forensic parameters

Population	Longitude	Latitude	Number (n)	diff Ht	D	h
Archangelskaja	40.36	64.35	42	36	0.857	0.9930
Brianskaja	34.22	53.18	43	43	1,000	1.0000
Ivanovskaja	41.01	57.02	40	40	1,000	1.0000
Lipezkaja	39.37	52.38	47	46	0.978	0.9991
Novgorodskaja	31.10	58.34	40	39	0.975	0.9987
Orlovskaja	36.04	52.58	42	42	1,000	1.0000
Penzenskaja	44.28	53.11	81	81	1,000	1.0000
Ryasanskaja	39.38	54.36	36	35	0.972	0.9984
Smolenskaja	32.03	54.47	43	43	1,000	1.0000
Tambovskaja	41.26	52.45	48	47	0.979	0.9991
Tverskaja	35.48	56.53	43	40	0.930	0.9967
Wologodskaja	39.53	59.15	40	39	0.975	0.9987
Pooled Sample			545	494	0.906	0.9994

diff Ht number of different haplotypes, D discrimination capacity, h haplotype diversity

Table 2

AMOVA pairwise distance based on Φ values between western Russia (n = 545) and nine neighbouring populations (n = 5,366)

	Finland	Siberia	Caucasus	Ukraine	Belarus	Poland	Latvia	Lithuania	Estonia	Western Russia
Finland	–	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
Siberia	0.1530	–	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
Caucasus	0.3305	0.4515	–	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
Ukraine	0.2972	0.4276	0.1414	–	0.3346*	0.0000	0.0000	0.0000	0.0000	0.0001
Belarus	0.2680	0.4003	0.1459	0.0003*	–	0.0007	0.0000	0.0000	0.0000	0.0952*
Poland	0.2598	0.4145	0.2200	0.0210	0.0104	–	0.0000	0.0000	0.0000	0.0000
Latvia	0.1674	0.3026	0.2164	0.0637	0.0394	0.0299	–	0.1907*	0.0022	0.0005
Lithuania	0.1907	0.3034	0.2667	0.0943	0.0654	0.0436	0.0032*	–	0.0000	0.0000
Estonia	0.0670	0.2452	0.2071	0.1159	0.0884	0.0985	0.0261	0.0516	–	0.0000
Western Russia	0.2152	0.3688	0.1672	0.0123	0.0023*	0.0089	0.0197	0.0413	0.0586	–

p values are shown above, Φ values below diagonal

*p > 0.05

Fig. 2

Multi-dimensional scaling (MDS) plot of the Western Russian metapopulation and nine reference populations, from pairwise Φ values. Stress value = 0.002. Acronyms are as follows: Rus–Russians, Bel–Belarus, Ukr–Ukrainians, Lat–Latvians, Lit–Lithuanians, Est–Estonians, Cau-Caucasus, Pol–Polish, Fin-Finns, Sib–Siberia

Results and discussion

A total of 545 samples from 12 Western Russian populations were investigated in this study and 494 different 17 locus haplotypes (D = 0.906, h = 0.9994) were detected. Two haplotypes occurred six times, one occurred five times, two occurred four times, three occurred three times, 25 twice and 461 were unique (Supplementary Table S1). The eight most frequent haplotypes occurring 34 times in this sample are closely related and belong (with the exception of one haplotype belonging to hg N3-Tat) to haplogroup R1a1-M17 (printed in bold in Table S1). This typical Eastern European haplogroup is the most frequent in all analysed populations with frequencies ranging between 0.31 and 0.56, followed by hg I-M170 (0.09–0.31) and N3-Tat (0.06–0.29) [5]. We observed high haplotype diversities (h) in all 12 populations ranging from 0.993 to 1.000 (Table 1). Three duplications (12 and 13 at locus GATA H4, 15 and 16 at DYS19 and 20 and 21 at DYS448), one “null” allele at DYS19 and several intermediate-sized alleles were observed.

No genetic substructure was found among the Russian populations. All pairwise comparisons were non-significant (p > 0.05). In contrast, significant variation between populations was observed in the comparison of Western Russia, treated as one homogeneous metapopulation, with neighbouring groups (Table 2). Because only reduced haplotype formats were available for such reference populations, we performed AMOVA based on 545 minimal 9-locus haplotypes from Western Russia with the previously published 5,366 haplotypes from 11 neighbouring regions (in clockwise direction: Ukraine, Belarus, Lithuania, Latvia, Estonia, Finland, Siberia and the Caucasus region) [9-18]. All pairwise Φ comparisons between Russia and these neighbours (with the exception of Russia vs. Belarus) were significant with values ranging between 0.0089 (Russia vs. Poland) and 0.3688 (Russia vs. eight Altaic- and Uralic-speaking groups from Siberia). MDS plot based on pairwise Φ values shows a closely related core group consisting of Slavic-speaking populations (Russia, Ukraine, Belarus and Poland) with an elevated distance to Baltic populations and a large span to linguistically different groups from Estonia, Finland, Siberia and the Caucasus (Fig. 2). From these analyses, we conclude that autochthonous Russian-speaking populations residing for centuries in the European part of Russia can be pooled to form a representative regional reference database for assessment of Y chromosomal matches in forensic analyses. However, on a level of quite low but significant Φ values, Western Russian populations can be grouped together with a much larger metapopulation defined within the forensic YHRD. This genetically distinct “Eastern European” metapopulation (n = 5,993, YHRD release 22 from 2007-08-10) comprises 56 Balto-Slavic-speaking populations from Eastern Europe [2]. Because the population frequency of a given haplotype is positively correlated with the combined frequency of closely related (surrounding) haplotypes in the population of origin [20], forensic databases collecting haploid markers have to be tested for substructure. This and other studies demonstrate that AMOVA is a sensitive method to detect the extent of inter-population differences accrued from the genetic and demographic history of the populations. The assignment of each population sample to a set of populations sharing a common linguistic, demographic, geographic and genetic background (metapopulations) facilitates the statistical evaluation of haplotype matches due to a significant enlargement of sample sizes.

The data obtained in this study from Russian populations were submitted to the Y Chromosome Haplotype Reference Database (www.yhrd.org) and they were assigned to the Eurasian–European–Eastern European metapopulation [3].

Below is the link to the electronic supplementary material.

17 locus haplotypes and associated haplogroups [5] in 545 Russian males from 12 populations. The eight most frequent haplotypes all belonging to hg R1a1-M17 are marked in bold (XLS 189 kb).