The genetic variation in the R1a clade among the Ashkenazi Levites' Y chromosome

The Qatari population's genetic structure and gene flow as revealed by the Y chromosome

The Y-chromosome has been widely used in forensic genetic applications and human population genetic studies due to its uniparental origins. A large database on the Qatari population was created for comparison with other databases from the Arabian Peninsula, the Middle East, and Africa. We provide a study of 23 Y-STR loci included in PowerPlex Y23 (Promega, USA) that were genotyped to produce haplotypes in 379 unrelated males from Qatar, a country at the crossroads of migration patterns. Overall, the most polymorphic locus provided by the Promega kit was DYS458, with a genetic diversity value of 0.85 and a haplotype diversity of 0.998924. Athey's Haplogroup Predictor tool was used to predict haplogroups from Y-STR haplotypes in the Qatari population. In a median-joining network, the haplogroup J1 predominance (49%) in Qatar generated a star-like expansion cluster. The graph of population Q-matrix was developed using Y-STR data from 38 Middle Eastern and 97 African populations (11,305 individuals), and it demonstrated a stronger sub-grouping of countries within each ethnic group and showed the effect of Arabs on the indigenous Berbers of North Africa. The estimated migration rate between the Qatari and other Arabian populations was inferred using Bayesian coalescence theory in the Migrate-n program. According to the Gene Flow study, the main migration route was from Yemen to Kuwait through Qatar. Our research, using the PowerPlex Y23 database, shows the importance of gene diversity, as well as regional and social structuring, in determining the utility of demographic and forensic databases.

Breast cancer polygenic risk scores derived in White European populations are not calibrated for women of Ashkenazi Jewish descent

PURPOSE

Polygenic risk scores (PRSs) are a major component of accurate breast cancer (BC) risk prediction but require ethnicity-specific calibration. Ashkenazi Jewish (AJ) population is assumed to be of White European (WE) origin in some commercially available PRSs despite differing effect allele frequencies (EAFs). We conducted a case-control study of WE and AJ women from the Predicting Risk of Cancer at Screening Study. The Breast Cancer in Northern Israel Study provided a separate AJ population-based case-control validation series.

METHODS

All women underwent Illumina OncoArray single-nucleotide variation (SNV; formerly single-nucleotide polymorphism [SNP]) analysis. Two PRSs were assessed, SNV142 and SNV78. A total of 221 of 2243 WE women (discovery: cases = 111; controls = 110; validation: cases = 651; controls = 1772) and 221 AJ women (cases = 121; controls = 110) were included from the UK study; the Israeli series consisted of 2045 AJ women (cases = 1331; controls = 714). EAFs were obtained from the Genome Aggregation Database.

RESULTS

In the UK study, the mean SNV142 PRS demonstrated good calibration and discrimination in WE population, with mean PRS of 1.33 (95% CI 1.18-1.48) in cases and 1.01 (95% CI 0.89-1.13) in controls. In AJ women from Manchester, the mean PRS of 1.54 (1.38-1.70) in cases and 1.20 (1.08-1.32) in controls demonstrated good discrimination but overestimation of BC relative risk. After adjusting for EAFs for the AJ population, mean risk was corrected (mean SNV142 PRS cases = 1.30 [95% CI 1.16-1.44] and controls = 1.02 [95% CI 0.92-1.12]). This was recapitulated in the larger Israeli data set with good discrimination (area under the curve = 0.632 [95% CI 0.607-0.657] for SNV142).

CONCLUSION

AJ women should not be given BC relative risk predictions based on PRSs calibrated to EAFs from the WE population. PRSs need to be recalibrated using AJ-derived EAFs. A simple recalibration using the mean PRS adjustment ratio likely performs well.

Genome-wide data from medieval German Jews show that the Ashkenazi founder event pre-dated the 14th century

We report genome-wide data from 33 Ashkenazi Jews (AJ), dated to the 14^th century, obtained following a salvage excavation at the medieval Jewish cemetery of Erfurt, Germany. The Erfurt individuals are genetically similar to modern AJ, but they show more variability in Eastern European-related ancestry than modern AJ. A third of the Erfurt individuals carried a mitochondrial lineage common in modern AJ and eight carried pathogenic variants known to affect AJ today. These observations, together with high levels of runs of homozygosity, suggest that the Erfurt community had already experienced the major reduction in size that affected modern AJ. The Erfurt bottleneck was more severe, implying substructure in medieval AJ. Overall, our results suggest that the AJ founder event and the acquisition of the main sources of ancestry pre-dated the 14^th century and highlight late medieval genetic heterogeneity no longer present in modern AJ.

Genetic and phylogeographic evidence for Jewish Holocaust victims at the Sobibór death camp

Six million Jews were killed by Nazi Germany and its collaborators during World War II. Archaeological excavations in the area of the death camp in Sobibór, Poland, revealed ten sets of human skeletal remains presumptively assigned to Polish victims of the totalitarian regimes. However, their genetic analyses indicate that the remains are of Ashkenazi Jews murdered as part of the mass extermination of European Jews by the Nazi regime and not of otherwise hypothesised non-Jewish partisan combatants. In accordance with traditional Jewish rite, the remains were reburied in the presence of a Rabbi at the place of their discovery.

Origin and diffusion of human Y chromosome haplogroup J1-M267

Human Y chromosome haplogroup J1-M267 is a common male lineage in West Asia. One high-frequency region-encompassing the Arabian Peninsula, southern Mesopotamia, and the southern Levant-resides ~ 2000 km away from the other one found in the Caucasus. The region between them, although has a lower frequency, nevertheless demonstrates high genetic diversity. Studies associate this haplogroup with the spread of farming from the Fertile Crescent to Europe, the spread of mobile pastoralism in the desert regions of the Arabian Peninsula, the history of the Jews, and the spread of Islam. Here, we study past human male demography in West Asia with 172 high-coverage whole Y chromosome sequences and 889 genotyped samples of haplogroup J1-M267. We show that this haplogroup evolved ~ 20,000 years ago somewhere in northwestern Iran, the Caucasus, the Armenian Highland, and northern Mesopotamia. The major branch-J1a1a1-P58-evolved during the early Holocene ~ 9500 years ago somewhere in the Arabian Peninsula, the Levant, and southern Mesopotamia. Haplogroup J1-M267 expanded during the Chalcolithic, the Bronze Age, and the Iron Age. Most probably, the spread of Afro-Asiatic languages, the spread of mobile pastoralism in the arid zones, or both of these events together explain the distribution of haplogroup J1-M267 we see today in the southern regions of West Asia.

Subdividing Y-chromosome haplogroup R1a1 reveals Norse Viking dispersal lineages in Britain

The influence of Viking-Age migrants to the British Isles is obvious in archaeological and place-names evidence, but their demographic impact has been unclear. Autosomal genetic analyses support Norse Viking contributions to parts of Britain, but show no signal corresponding to the Danelaw, the region under Scandinavian administrative control from the ninth to eleventh centuries. Y-chromosome haplogroup R1a1 has been considered as a possible marker for Viking migrations because of its high frequency in peninsular Scandinavia (Norway and Sweden). Here we select ten Y-SNPs to discriminate informatively among hg R1a1 sub-haplogroups in Europe, analyse these in 619 hg R1a1 Y chromosomes including 163 from the British Isles, and also type 23 short-tandem repeats (Y-STRs) to assess internal diversity. We find three specifically Western-European sub-haplogroups, two of which predominate in Norway and Sweden, and are also found in Britain; star-like features in the STR networks of these lineages indicate histories of expansion. We ask whether geographical distributions of hg R1a1 overall, and of the two sub-lineages in particular, correlate with regions of Scandinavian influence within Britain. Neither shows any frequency difference between regions that have higher (≥10%) or lower autosomal contributions from Norway and Sweden, but both are significantly overrepresented in the region corresponding to the Danelaw. These differences between autosomal and Y-chromosomal histories suggest either male-specific contribution, or the influence of patrilocality. Comparison of modern DNA with recently available ancient DNA data supports the interpretation that two sub-lineages of hg R1a1 spread with the Vikings from peninsular Scandinavia.

Determination of the phylogenetic origins of the Árpád Dynasty based on Y chromosome sequencing of Béla the Third

We set out to identify the origins of the Árpád Dynasty based on genome sequencing of DNA derived from the skeletal remains of Hungarian King Béla III (1172–1196) and eight additional individuals (six males, two females) originally interred at the Royal Basilica of Székesfehérvár. Y-chromosome analysis established that two individuals, Béla III and HU52 assign to haplogroups R-Z2125 whose distribution centres near South Central Asia with subsidiary expansions in the regions of modern Iran, the Volga Ural region and the Caucasus. Out of a cohort of 4340 individuals from these geographic areas, we acquired whole-genome data from 208 individuals derived for the R-Z2123 haplogroup. From these data we have established that the closest living kin of the Árpád Dynasty are R-SUR51 derived modern day Bashkirs predominantly from the Burzyansky and Abzelilovsky districts of Bashkortostan in the Russian Federation. Our analysis also reveals the existence of SNPs defining a novel Árpád Dynasty specific haplogroup R-ARP. Framed within the context of a high resolution R-Z2123 phylogeny, the ancestry of the first Hungarian royal dynasty traces to the region centering near Northern Afghanistan about 4500 years ago and identifies the Bashkirs as their closest kin, with a separation date between the two populations at the beginning of the first millennium CE.

The spectrum of BRCA1 and BRCA2 pathogenic sequence variants in Middle Eastern, North African, and South European countries

BRCA1 BRCA2 mutational spectrum in the Middle East, North Africa, and Southern Europe is not well characterized. The unique history and cultural practices characterizing these regions, often involving consanguinity and inbreeding, plausibly led to the accumulation of population-specific founder pathogenic sequence variants (PSVs). To determine recurring BRCA PSVs in these locales, a search in PUBMED, EMBASE, BIC, and CIMBA was carried out combined with outreach to researchers from the relevant countries for unpublished data. We identified 232 PSVs in BRCA1 and 239 in BRCA2 in 25 of 33 countries surveyed. Common PSVs that were detected in four or more countries were c.5266dup (p.Gln1756Profs), c.181T>G (p.Cys61Gly), c.68_69del (p.Glu23Valfs), c.5030_5033del (p.Thr1677Ilefs), c.4327C>T (p.Arg1443Ter), c.5251C>T (p.Arg1751Ter), c.1016dup (p.Val340Glyfs), c.3700_3704del (p.Val1234Glnfs), c.4065_4068del (p.Asn1355Lysfs), c.1504_1508del (p.Leu502Alafs), c.843_846del (p.Ser282Tyrfs), c.798_799del (p.Ser267Lysfs), and c.3607C>T (p.Arg1203Ter) in BRCA1 and c.2808_2811del (p.Ala938Profs), c.5722_5723del (p.Leu1908Argfs), c.9097dup (p.Thr3033Asnfs), c.1310_1313del (p. p.Lys437Ilefs), and c.5946del (p.Ser1982Argfs) for BRCA2. Notably, some mutations (e.g., p.Asn257Lysfs (c.771_775del)) were observed in unrelated populations. Thus, seemingly genotyping recurring BRCA PSVs in specific populations may provide first pass BRCA genotyping platform.

Estimating Y-Str Mutation Rates and Tmrca Through Deep-Rooting Italian Pedigrees

In the population genomics era, the study of Y-chromosome variability is still of the greatest interest for several fields ranging from molecular anthropology to forensics and genetic genealogy. In particular, mutation rates of Y-chromosomal Short Tandem Repeats markers (Y-STRs) are key parameters for different interdisciplinary applications. Among them, testing the patrilineal relatedness between individuals and calculating their Time of Most Recent Common Ancestors (TMRCAs) are of the utmost importance. To provide new valuable estimates and to address these issues, we typed 47 Y-STRs (comprising Yfiler, PowerPlex23 and YfilerPlus loci, the recently defined Rapidly Mutating [RM] panel and 11 additional markers often used in genetic genealogical applications) in 135 individuals belonging to 66 deep-rooting paternal genealogies from Northern Italy. Our results confirmed that the genealogy approach is an effective way to obtain reliable Y-STR mutation rate estimates even with a limited number of samples. Moreover, they showed that the impact of multi-step mutations and backmutations is negligible within the temporal scale usually adopted by forensic and genetic genealogy analyses. We then detected a significant association between the number of mutations within genealogies and observed TMRCAs. Therefore, we compared observed and expected TMRCAs by implementing a Bayesian procedure originally designed by Walsh (2001) and showed that the method yields a good performance (up to 96.72%), especially when using the Infinite Alleles Model (IAM).