Improving structural variant clustering to reduce the negative effect of the breakpoint uncertainty problem.

BACKGROUND: Structural variants (SVs) represent an important source of genetic variation. One of the most critical problems in their detection is breakpoint uncertainty associated with the inability to determine their exact genomic position. Breakpoint uncertainty is a characteristic issue of structural variants detected via short-read sequencing methods and complicates subsequent population analyses. The commonly used heuristic strategy reduces this issue by clustering/merging nearby structural variants of the same type before the data from individual samples are merged. RESULTS: We compared the two most used dissimilarity measures for SV clustering in terms of Mendelian inheritance errors (MIE), kinship prediction, and deviation from Hardy-Weinberg equilibrium. We analyzed the occurrence of Mendelian-inconsistent SV clusters that can be collapsed into one Mendelian-consistent SV as a new measure of dataset consistency. We also developed a new method based on constrained clustering that explicitly identifies these types of clusters. CONCLUSIONS: We found that the dissimilarity measure based on the distance between SVs breakpoints produces slightly better results than the measure based on SVs overlap. This difference is evident in trivial and corrected clustering strategy, but not in constrained clustering strategy. However, constrained clustering strategy provided the best results in all aspects, regardless of the dissimilarity measure used.

Subsistence strategy was the main factor driving population differentiation in the bidirectional corridor of the African Sahel.

OBJECTIVES: The Sahel belt is occupied by populations who use two types of subsistence strategy, nomadic pastoralism and sedentary farming, and who belong to three linguistic families, Niger-Congo, Nilo-Saharan, and Afro-Asiatic. Little is known, however, about the origins of these two populations and their mutual genetic relationships. MATERIALS AND METHODS: We have built a large dataset of mitochondrial DNA sequences and Y chromosomal STR haplotypes of pastoralists and farmers belonging to all three linguistic phyla in the western, central, and eastern parts of the Sahel. We calculated pairwise genetic, geographic, and linguistic distances between populations and analyzed the effects of geography, language, and subsistence on population genetic structure. RESULTS: We found that subsistence mode significantly contributed to the generally low population structure in the Sahel and that language affiliation plays a more important role for pastoralists than for farmers. We also demonstrated that geographic isolation significantly influenced the population structure of sedentary farmers but not of nomadic pastoralists. Finally, we found haplotypes shared between the Fulani and Arabic-speaking Baggara, supporting the theory of Baggarization, which explains the recent adaptation of Arabic-speaking nomads in the Sahel region through contact with autochthonous sub-Saharan populations. CONCLUSIONS: Based on various genetic and archaeological evidence pertaining to the Sahel, we suggest that the idea of a bidirectional Sahelian corridor is valid, but that pastoralists made a more important contribution to its population structure. It is also possible that agropastoralists diverged into farmers and pastoralists in the early stages of formation of the Sahelian gene pool.

Dispersals of the Siberian Y-chromosome haplogroup Q in Eurasia.

The human Y-chromosome has proven to be a powerful tool for tracing the paternal history of human populations and genealogical ancestors. The human Y-chromosome haplogroup Q is the most frequent haplogroup in the Americas. Previous studies have traced the origin of haplogroup Q to the region around Central Asia and Southern Siberia. Although the diversity of haplogroup Q in the Americas has been studied in detail, investigations on the diffusion of haplogroup Q in Eurasia and Africa are still limited. In this study, we collected 39 samples from China and Russia, investigated 432 samples from previous studies of haplogroup Q, and analyzed the single nucleotide polymorphism (SNP) subclades Q1a1a1-M120, Q1a2a1-L54, Q1a1b-M25, Q1a2-M346, Q1a2a1a2-L804, Q1a2b2-F1161, Q1b1a-M378, and Q1b1a1-L245. Through NETWORK and BATWING analyses, we found that the subclades of haplogroup Q continued to disperse from Central Asia and Southern Siberia during the past 10,000 years. Apart from its migration through the Beringia to the Americas, haplogroup Q also moved from Asia to the south and to the west during the Neolithic period, and subsequently to the whole of Eurasia and part of Africa.

Internal diversification of non-Sub-Saharan haplogroups in Sahelian populations and the spread of pastoralism beyond the Sahara.

BACKGROUND: Today, African pastoralists are found mainly in the Sahel/Savannah belt spanning 6,000 km from west to east, flanked by the Sahara to the north and tropical rainforests to the south. The most significant group among them are the Fulani who not only keep cattle breeds of possible West Eurasian ancestry, but form themselves a gene pool containing some paternally and maternally-transmitted West Eurasian haplogroups. MATERIALS AND METHODS: We generated complete sequences for 33 mitogenomes belonging to haplogroups H1 and U5 (23 and 10, respectively), and genotyped 16 STRs in 65 Y chromosomes belonging to haplogroup R1b-V88. RESULTS: We show that age estimates of the maternal lineage H1cb1, occurring almost exclusively in the Fulani, point to the time when the first cattle herders settled the Sahel/Savannah belt. Similar age estimates were obtained for paternal lineage R1b-V88, which occurs today in the Fulani but also in other, mostly pastoral populations. Maternal clade U5b1b1b, reported earlier in the Berbers, shows a shallower age, suggesting another possibly independent input into the Sahelian pastoralist gene pool. CONCLUSIONS: Despite the fact that animal domestication originated in the Near East ∼ 10 ka, and that it was from there that animals such as sheep, goats as well as cattle were introduced into Northeast Africa soon thereafter, contemporary cattle keepers in the Sahel/Savannah belt show uniparental genetic affinities that suggest the possibility of an ancient contact with an additional ancestral population of western Mediterranean ancestry.

Analysis of 8 X-chromosomal markers in the population of central Croatia.

AIM: To analyze 8 X-linked short tandem repeat (STR) markers in the population of central Croatia and to evaluate their forensic efficiency. METHODS: We carried out a statistical analysis of the data from previously performed genetic analyses, collected during routine forensic work by the Forensic Science Centre ''Ivan Vucetic.'' Mentype® Argus X-8 PCR amplification kit was used for typing the data of 99 unrelated healthy women and 78 men from central Croatia. Haplotype frequencies were calculated only in male samples. Arlequin 3.5 software was used to assess Hardy-Weinberg equilibrium (HWE), linkage disequilibrium (LD), observed and expected heterozygosity. Power of discrimination (PD) for men and women, polymorphism information content (PIC), power of exclusion, and mean exclusion chance for deficiency cases, normal trios, and duos were determined using online database ChrX-STR.org. RESULTS: In female samples, deviations from HWE (P=0.006) for each locus were not found. LD test performed both on female and male samples revealed no significant association between markers (P=0.002). DXS10135 was the most polymorphic locus (PIC=0.931). PD varied from 0.692 to 0.935 in male and from 0.845 to 0.992 in female samples. Combined PD reached 99.999999% in men and 99.9999999999% in women. CONCLUSION: Performed analyses revealed that the studied marker set contained polymorphic markers with high power of discrimination. We can conclude that Mentype® Argus X-8 PCR is suitable for application in the population of central Croatia. Results of this study, together with collected allele and haplotype frequencies, are the first step in establishing a national reference X-STR database based on 8 X-STR loci.

Croatian national reference Y-STR haplotype database.

A reference Y-chromosome short tandem repeat (STR) haplotype database is needed for Y-STR match interpretation as well as for national and regional characterization of populations. The aim of this study was to create a comprehensive Y-STR haplotype database of the Croatian contemporary population and to analyze substructure between the five Croatian regions. We carried out a statistical analysis of the data from previously performed genetic analyses collected during routine forensic work by the Forensic Science Centre "Ivan Vucetic". A total of 1,100 unrelated men from eastern, western, northern, southern and central Croatia were selected for the purpose of this study. Y-STRs were typed using the AmpFISTR Yfiler PCR amplification kit. Analysis of molecular variance calculated with the Y chromosome haplotype reference database online analysis tool included 16 population samples with 20,247 haplotypes. A total of 947 haplotypes were recorded, 848 of which were unique (89.5%). Haplotype diversity was 0.998, with the most frequent haplotype found in 9 of 1,100 men (0.82%). Locus diversity varied from 0.266 for DYS392 to 0.868 for DYS385. Discrimination capacity was 86.1%. Our results suggested high level of similarity among regional subpopulations within Croatia, except for mildly different southern Croatia. Relative resemblance was found with Bosnia and Herzegovina and Serbia. Whit Atheys' Haplogroup Predictor was used to estimate the frequencies of Y-chromosome haplogroups. I2a, R1a, E1b1b and R1b haplogroups were most frequent in all Croatian regions. These results are important in forensics and contribute to the population genetics and genetic background of the contemporary Croatian population.

Y-chromosomal diversity of the Valachs from the Czech Republic: model for isolated population in Central Europe.

AIM: To evaluate Y-chromosomal diversity of the Moravian Valachs of the Czech Republic and compare them with a Czech population sample and other samples from Central and South-Eastern Europe, and to evaluate the effects of genetic isolation and sampling. METHODS: The first sample set of the Valachs consisted of 94 unrelated male donors from the Valach region in northeastern Czech Republic border-area. The second sample set of the Valachs consisted of 79 men who originated from 7 paternal lineages defined by surname. No close relatives were sampled. The third sample set consisted of 273 unrelated men from the whole of the Czech Republic and was used for comparison, as well as published data for other 27 populations. The total number of samples was 3244. Y-short tandem repeat (STR) markers were typed by standard methods using PowerPlex® Y System (Promega) and Yfiler® Amplification Kit (Applied Biosystems) kits. Y-chromosomal haplogroups were estimated from the haplotype information. Haplotype diversity and other intra- and inter-population statistics were computed. RESULTS: The Moravian Valachs showed a lower genetic variability of Y-STR markers than other Central European populations, resembling more to the isolated Balkan populations (Aromuns, Csango, Bulgarian, and Macedonian Roma) than the surrounding populations (Czechs, Slovaks, Poles, Saxons). We illustrated the effect of sampling on Valach paternal lineages, which includes reduction of discrimination capacity and variability inside Y-chromosomal haplogroups. Valach modal haplotype belongs to R1a haplogroup and it was not detected in the Czech population. CONCLUSION: The Moravian Valachs display strong substructure and isolation in their Y chromosomal markers. They represent a unique Central European population model for population genetics.

Allele frequency data for 17 short tandem repeats in a Czech population sample.

Allele frequencies for 17 short tandem repeats (STRs) autosomal loci (D2S1338, D3S1358, D5S818, D7S820, D8S1179, D13S317, D16S539, D18S51, D19S433, D21S11, CSF1PO, FGA, PentaD, PentaE, TH01, TPOX, vWA) were studied in an extensive sample (max. N=1411) of unrelated individuals originating from the Czech Republic. Population and forensic parameters were estimated. Except for FGA and Penta E loci, no deviations from the Hardy-Weinberg equilibrium were detected. A comparative analysis with published data revealed significant differences in allele frequencies for some loci from the Polish population and three Hungarian populations (Ashkenazim population and Romany populations from Debrecen and Baranya County, respectively). A combination of these 17 STR loci provides a powerful tool for forensic identification in the native Czech population.