Search

BackgroundThe Sahelian Fulani are the largest nomadic pastoral ethnic group. Their origins are still largely unknown and their Eurasian genetic component is usually explained by recent admixture events with northern African groups. However, it has also been proposed that Fulani may be the descendants of ancient groups settled in the Sahara during its last Green phase (12000-5000 BP), as also suggested by Y chromosome results.ResultsWe produced 23 high-coverage (30 ×) whole genomes from Fulani individuals from 8 Sahelian countries, plus 17 samples from other African groups and 3 Europeans as controls, for a total of 43 new whole genome sequences. These data have been compared with 814 published modern whole genomes and analyzed together with relevant published ancient individuals (for a total of > 1800 samples). These analyses showed that the non-sub-Saharan genetic ancestry component of Fulani cannot be only explained by recent admixture events, but it could be shaped at least in part by older events by events more ancient than previously reported, possibly tracing its origin to the last Green Sahara.ConclusionsAccording to our results, Fulani may be the descendants of Saharan cattle herders settled in that area during the last Green Sahara. The exact ancestry composition of such ghost Saharan population(s) cannot be completely unveiled from modern genomes only, but the joint analysis with the available African ancient samples suggested a similarity between ancient Saharans and Late Neolithic Moroccans.

The extent of the devastation of the Black Death pandemic (1346-53) on European populations is known from documentary sources and its bacterial source illuminated by studies of ancient pathogen DNA. What has remained less understood is the effect of the pandemic on human mobility and genetic diversity at local scale in the context of the social stratification of medieval communities. Here we study 275 newly reported ancient genomes from later medieval and post-medieval Cambridgeshire, from individuals buried before, during, and after the Black Death. The majority of individuals examined had local genetic ancestries. Consistent with the function of the institutions, we found a lack of close relatives among the friars and the inmates of the hospital in contrast to their abundance in general urban and rural parish communities. Accounting for the genetic component for height accentuates the disparities between social groups in stature estimated from long bones, as a proxy for health and the quality of life. While we detect long-term shifts in local genetic ancestry in Cambridgeshire that either pre- or postdate the Black Death, we find no evidence of major changes in genetic ancestry nor, in contrast to recent claims, higher differentiation of immune loci between cohorts living before and after the Black Death.

The Finnish population is a unique example of a genetic isolate affected by a recent founder event. Previous studies have suggested that the ancestors of Finnic-speaking Finns and Estonians reached the circum-Baltic region by the 1(st) millennium BC. However, high linguistic similarity points to a more recent split of their languages. To study genetic connectedness between Finns and Estonians directly, we first assessed the efficacy of imputation of low-coverage ancient genomes by sequencing a medieval Estonian genome to high depth (23×) and evaluated the performance of its down-sampled replicas. We find that ancient genomes imputed from >0.1× coverage can be reliably used in principal-component analyses without projection. By searching for long shared allele intervals (LSAIs; similar to identity-by-descent segments) in unphased data for >143,000 present-day Estonians, 99 Finns, and 14 imputed ancient genomes from Estonia, we find unexpectedly high levels of individual connectedness between Estonians and Finns for the last eight centuries in contrast to their clear differentiation by allele frequencies. High levels of sharing of these segments between Estonians and Finns predate the demographic expansion and late settlement process of Finland. One plausible source of this extensive sharing is the 8(th)–10(th) centuries AD migration event from North Estonia to Finland that has been proposed to explain uniquely shared linguistic features between the Finnish language and the northern dialect of Estonian and shared Christianity-related loanwords from Slavic. These results suggest that LSAI detection provides a computationally tractable way to detect fine-scale structure in large cohorts.

About one-quarter of the euchromatic portion of the male-specific region of the human Y chromosome consists of large duplicated sequences that are organized in eight palindromes (termed P1–P8), which undergo arm-to arm gene conversion, a proposed mechanism for maintaining their sequence integrity. Although the relevance of gene conversion in the evolution of palindromic sequences has been profoundly recognized, the dynamic of this mechanism is still nuanced. To shed light into the evolution of these genomic elements, we performed a high-depth (50×) targeted next-generation sequencing of the palindrome P6 in 157 subjects belonging to the most divergent evolutionary lineages of the Y chromosome. We found 118 new paralogous sequence variants, which were placed into the context of a robust Y chromosome phylogeny based on 7240 SNPs of the X-degenerate region. We mapped along the phylogeny 80 gene conversion events that shaped the diversity of P6 arms during recent human history. In contrast to previous studies, we demonstrated that arm-to-arm gene conversion, which occurs at a rate of 6.01 × 10 −6 conversions/base/year, is not biased toward the retention of the ancestral state of sequences. We also found a significantly lower mutation rate of the arms (6.18 × 10−10 mutations/base/year) compared with the spacer (9.16 × 10−10 mutations/base/year), a finding that may explain the observed higher inter-species conservation of arms, without invoking any bias of conversion. Finally, by formally testing the mutation/conversion balance in P6, we found that the arms of this palindrome reached a steady-state equilibrium between mutation and gene conversion.

The human MSY ampliconic region is mainly composed of large duplicated sequences that are organized in eight palindromes (termed P1–P8), and may undergo arm-to-arm gene conversion. Although the importance of these elements is widely recognized, their evolutionary dynamics are still nuanced. Here, we focused on the P8 palindrome, which shows a complex evolutionary history, being involved in intra- and inter-chromosomal gene conversion. To disclose its evolutionary complexity, we performed a high-depth (50×) targeted next-generation sequencing of this element in 157 subjects belonging to the most divergent lineages of the Y chromosome tree. We found a total of 72 polymorphic paralogous sequence variants that have been exploited to identify 41 Y-Y gene conversion events that occurred during recent human history. Through our analysis, we were able to categorize P8 arms into three portions, whose molecular diversity was modelled by different evolutionary forces. Notably, the outer region of the palindrome is not involved in any gene conversion event and evolves exclusively through the action of mutational pressure. The inner region is affected by Y-Y gene conversion occurring at a rate of 1.52 × 10−5 conversions/base/year, with no bias towards the retention of the ancestral state of the sequence. In this portion, GC-biased gene conversion is counterbalanced by a mutational bias towards AT bases. Finally, the middle region of the arms, in addition to intra-chromosomal gene conversion, is involved in X-to-Y gene conversion (at a rate of 6.013 × 10−8 conversions/base/year) thus being a major force in the evolution of the VCY/VCX gene family.

Y chromosome short tandem repeats (Y-STRs) typing is becoming increasingly popular in forensic casework mainly because it allows the recovery of male-specific genetic information from severely unbalanced male-female DNA mixtures. The relatively low discrimination power of conventional Y-STR multiplexes, due to linkage disequilibrium among polymorphic loci, has been partially overcome by the introduction of rapidly mutating Y microsatellites (RM Y-STRs) with mutation rates exceeding 1 × 10

Human herpes simplex virus 1 (HSV-1), a life-long infection spread by oral contact, today infects a majority of adults globally1, yet no ancient HSV-1 genomes have yet been published. Phylogeographic clustering of sampled diversity into European, pan-Eurasian, and African groups2, 3 has suggested that the virus co-diverged with anatomically modern humans migrating out of Africa4, although a much younger origin has also been proposed5. The lack of ancient HSV-1 genomes, high rates of recombination, and high mobility of humans in the modern era have impeded the understanding of HSV-1’s evolutionary history. Here we present three full ancient European HSV-1 genomes and one partial genome, dating to between the 3rd and 17th century CE, sequenced to up to 9.5× with paired human genomes up to 10.16×. These HSV-1 strains fall within modern Eurasian diversity. We estimate a mean mutation rate of 7.6 × 10-7 Introduction Results - Retrieved genomes are likely from typical infections - Demographic history of HSV-1 in a global context Discussion Material and Methods - Ethics statement - Sampling - Generation of aDNA libraries - Sequencing - aDNA authentication - Metagenomic screening - Targeted capture of HSV-1 - Alignment of viral data to the reference sequence - Genotyping - HSV-1 linkage disequilibrium and population genetic analysis - Compilation of comparative HSV data - Preparation of genome sequences - HSV-1 phylogenetic analysis and recombination filtering - Phylogenetic dating - Alignment of human data to the reference sequence and quality control - Genetic sex estimation, mtDNA, and Y haplotyping - Human variant calling and imputation of genotypes

Across Europe, the genetics of the Chalcolithic/Bronze Age transition is increasingly characterized in terms of an influx of Steppe-related ancestry. The effect of this major shift on the genetic structure of populations in the Italian Peninsula remains underexplored. Here, genome-wide shotgun data for 22 individuals from commingled cave and single burials in Northeastern and Central Italy dated between 3200 and 1500 BCE provide the first genomic characterization of Bronze Age individuals (n = 8; 0.001-1.2× coverage) from the central Italian Peninsula, filling a gap in the literature between 1950 and 1500 BCE. Our study confirms a diversity of ancestry components during the Chalcolithic and the arrival of Steppe-related ancestry in the central Italian Peninsula as early as 1600 BCE, with this ancestry component increasing through time. We detect close patrilineal kinship in the burial patterns of Chalcolithic commingled cave burials and a shift away from this in the Bronze Age (2200-900 BCE) along with lowered runs of homozygosity, which may reflect larger changes in population structure. Finally, we find no evidence that the arrival of Steppe-related ancestry in Central Italy directly led to changes in frequency of 115 phenotypes present in the dataset, rather that the post-Roman Imperial period had a stronger influence, particularly on the frequency of variants associated with protection against Hansen's disease (leprosy). Our study provides a closer look at local dynamics of demography and phenotypic shifts as they occurred as part of a broader phenomenon of widespread admixture during the Chalcolithic/Bronze Age transition. ispartof: CURRENT BIOLOGY vol:31 issue:12 pages:2576-+ ispartof: location:England status: published

aDNA reveals Eastern hunter-gatherers in NW Russia after Ice Age, Steppe and farmer ancestries in first forest pastoralists., The transition from Stone to Bronze Age in Central and Western Europe was a period of major population movements originating from the Ponto-Caspian Steppe. Here, we report new genome-wide sequence data from 30 individuals north of this area, from the understudied western part of present-day Russia, including 3 Stone Age hunter-gatherers (10,800 to 4250 cal BCE) and 26 Bronze Age farmers from the Corded Ware complex Fatyanovo Culture (2900 to 2050 cal BCE). We show that Eastern hunter-gatherer ancestry was present in northwestern Russia already from around 10,000 BCE. Furthermore, we see a change in ancestry with the arrival of farming—Fatyanovo Culture individuals were genetically similar to other Corded Ware cultures, carrying a mixture of Steppe and European early farmer ancestry. Thus, they likely originate from a fast migration toward the northeast from somewhere near modern-day Ukraine—the closest area where these ancestries coexisted from around 3000 BCE.

Funder: Sapienza Università di Roma, Although ancient DNA data have become increasingly more important in studies about past populations, it is often not feasible or practical to obtain high coverage genomes from poorly preserved samples. While methods of accurate genotype imputation from > 1 × coverage data have recently become a routine, a large proportion of ancient samples remain unusable for downstream analyses due to their low coverage. Here, we evaluate a two-step pipeline for the imputation of common variants in ancient genomes at 0.05–1 × coverage. We use the genotype likelihood input mode in Beagle and filter for confident genotypes as the input to impute missing genotypes. This procedure, when tested on ancient genomes, outperforms a single-step imputation from genotype likelihoods, suggesting that current genotype callers do not fully account for errors in ancient sequences and additional quality controls can be beneficial. We compared the effect of various genotype likelihood calling methods, post-calling, pre-imputation and post-imputation filters, different reference panels, as well as different imputation tools. In a Neolithic Hungarian genome, we obtain ~ 90% imputation accuracy for heterozygous common variants at coverage 0.05 × and > 97% accuracy at coverage 0.5 ×. We show that imputation can mitigate, though not eliminate reference bias in ultra-low coverage ancient genomes.

Transition from the Stone to the Bronze Age in Central and Western Europe was a period of major population movements originating from the Ponto-Caspian Steppe. Here, we report new genome-wide sequence data from 28 individuals from the territory north of this source area – from the under-studied Western part of present-day Russia, including Stone Age hunter-gatherers (10,800–4,250 cal BC) and Bronze Age farmers from the Corded Ware complex called Fatyanovo Culture (2,900–2,050 cal BC). We show that Eastern hunter-gatherer ancestry was present in Northwestern Russia already from around 10,000 BC. Furthermore, we see a clear change in ancestry with the arrival of farming – the Fatyanovo Culture individuals were genetically similar to other Corded Ware cultures, carrying a mixture of Steppe and European early farmer ancestry and thus likely originating from a fast migration towards the northeast from somewhere in the vicinity of modern-day Ukraine, which is the closest area where these ancestries coexisted from around 3,000 BC.

Y chromosome short tandem repeats (Y-STRs) are commonly used to identify male lineages for investigative and judicial purposes and could represent the only source of male-specific genetic information from unbalanced female-male mixtures. The Yfiler Plus multiplex, which includes twenty conventional and seven rapidly-mutating Y-STRs, represents the most discriminating patrilineal system commercially available to date. Over the past five years, this multiplex has been used to analyze several Eurasian populations, with a reported discrimination capacity (DC) approaching or corresponding to the highest possible value. However, despite the inclusion of rapidly mutating Y-STRs, extensive haplotype sharing was still reported for some African populations due to a number of different factors affecting the effective population size. In the present study, we analyzed 27 Y-STRs included in the Yfiler Plus multiplex and 82 Y-SNPs in central Sahel (northern Cameroon and western Chad), an African region characterized by a strong ethnic fragmentation and linguistic diversity. We evaluated the effects of population sub-structuring on genetic diversity by stratifying a sample composed of 431 males according to their ethnicity (44 different ethnic groups) and urbanization degree (four villages and four towns). Overall, we observed a low discrimination capacity (DC = 0.90), with 71 subjects (16.5 %) sharing 27 Y-STR haplotypes. Haplotype sharing was essentially limited to subjects with the same binary haplogroup, coming from the same location and belonging to the same ethnic group. Haplotype sharing was much higher in rural areas (average DC = 0.83) than urban settlements (average DC = 0.96) with a significant correlation between DC and census size (r = 0.89; p = 0.003). Notably, we found that genetic differentiation between villages from the same country ( Φ ST = 0.14) largely exceeded that found among countries ( Φ ST = 0.02). These findings have important implications for the choice of the appropriate reference population database to evaluate the statistical relevance of forensic Y-haplotype matches.

The male-specific northern African genetic pool is characterised by a high frequency of the E-M81 haplogroup, which expanded in very recent times (2-3 kiloyears ago). As a consequence of their recent coalescence, E-M81 chromosomes often cannot be completely distinguished on the basis of their Y-STR profiles, unless rapidly-mutating Y-STRs (RM Y-STRs) are analysed. In this study, we used the Yfiler® Plus kit, which includes 7 RM Y-STRs and 20 standard Y-STR, to analyse 477 unrelated males coming from 11 northern African populations sampled from Morocco, Algeria, Libya and Egypt. The Y chromosomes were assigned to monophyletic lineages after the analysis of 72 stable biallelic polymorphisms and, as expected, we found a high proportion of E-M81 subjects (about 46%), with frequencies decreasing from west to east. We found low intra-population diversity indexes, in particular in the populations that experienced long-term isolation. The AMOVA analysis showed significant differences between the countries and between most of the 11 populations, with a rough differentiation between northwestern Africa and northeastern Africa, where the Egyptians Berbers from Siwa represented an outlier population. The comparison between the Yfiler® and the Yfiler® Plus network of the E-M81 Y chromosomes confirmed the high power of discrimination of the latter kit, thanks to higher variability of the RM Y-STRs: indeed, the number of chromosomes sharing the same haplotype was drastically reduced from 201 to 81 and limited, in the latter case, to subjects from the same population.

Recently, thirteen rapidly mutating Y-STRs (RM Y-STRs), characterized by a mutation rate higher than 10−2/STR/generation, have been proven to be extremely useful in distinguishing among close male relatives. Six of these RM-YSTRs have been included in the Yfiler™ Plus PCR Amplification Kit, which shows the highest discriminatory power among Y-STR commercial multiplex currently available. In this study, we used Yfiler Plus to analyze 1437 males from eastern, northern and sub-Saharan Africa (462, 477 and 498 subjects, respectively). 242 out of 1437 subjects were found to share 102 Y-STR haplotypes, resulting in a low discrimination capacity (DC = 0.90) as compared to other continental regions. With few exceptions, Y-STR haplotype sharing was limited to subjects coming from the same country and ethnic group, and belonging to the same binary Y-SNP haplogroup, suggesting possible familial relationships along the male lineage. In order to evaluate the presence of hidden familial relationships, all the 242 subjects sharing a Y-STR haplotype were further genotyped for 16 autosomal STRs using the AmpFlSTR® NGM SElect™ PCR Amplification Kit. Blind search and simulation test analyses for kinship using the Familias software revealed the presence of a low proportion of close relatives (second degree or closer) in our sample set. These findings show that close relatedness explains only a relatively small proportion of the observed Y-STR haplotype sharing, suggesting that a higher number of RM Y-STRs should be included in commercial kits to improve the discrimination power of male-specific markers in regions characterized by high levels of endogamy.

In order to improve the phylogeography of the male-specific genetic traces of Greek and Phoenician colonizations on the Northern coasts of the Mediterranean, we performed a geographically structured sampling of seven subclades of haplogroup J in Turkey, Greece and Italy. We resequenced 4.4 Mb of Y-chromosome in 58 subjects, obtaining 1079 high quality variants. We did not find a preferential coalescence of Turkish samples to ancestral nodes, contradicting the simplistic idea of a dispersal and radiation of Hg J as a whole from the Middle East. Upon calibration with an ancient Hg J chromosome, we confirmed that signs of Holocenic Hg J radiations are subtle and date mainly to the Bronze Age. We pinpointed seven variants which could potentially unveil star clusters of sequences, indicative of local expansions. By directly genotyping these variants in Hg J carriers and complementing with published resequenced chromosomes (893 subjects), we provide strong temporal and distributional evidence for markers of the Greek settlement of Magna Graecia (J2a-L397) and Phoenician migrations (rs760148062). Our work generated a minimal but robust list of evolutionarily stable markers to elucidate the demographic dynamics and spatial domains of male-mediated movements across and around the Mediterranean, in the last 6,000 years.

By using the recently introduced 6-dye Yfiler® Plus multiplex, we analyzed 462 males belonging to 20 ethnic groups from four eastern African countries (Eritrea, Ethiopia, Djibouti and Kenya). Through a Y-STR sequence analysis, combined with 62 SNP-based haplogroup information, we were able to classify observed microvariant alleles at four Y-STR loci as either monophyletic (DYF387S1 and DYS458) or recurrent (DYS449 and DYS627). We found evidence of non-allelic gene conversion among paralogous STRs of the two-copy locus DYF387S1. Twenty-two diallelic and triallelic patterns observed at 13 different loci were found to be significantly over-represented (p

Long Terminal Repeats (LTRs) are nearly identical DNA sequences found at either end of Human Endogenous Retroviruses (HERVs). The high sequence similarity that exists among different LTRs suggests they could be substrate of ectopic gene conversion events. To understand the extent to which gene conversion occurs and to gain new insights into the evolutionary history of these elements in humans, we performed an intra-species phylogenetic study of 52 LTRs on different unrelated Y chromosomes. From this analysis, we obtained direct evidence that demonstrates the occurrence of ectopic gene conversion in several LTRs, with donor sequences located on both sex chromosomes and autosomes. We also found that some of these elements are characterized by an extremely high density of polymorphisms, showing one of the highest nucleotide diversities in the human genome, as well as a complex patchwork of sequences derived from different LTRs. Finally, we highlighted the limits of current short-read NGS studies in the analysis of genetic diversity of the LTRs in the human genome. In conclusion, our comparative re-sequencing analysis revealed that ectopic gene conversion is a common event in the evolution of LTR elements, suggesting complex genetic links among LTRs from different chromosomes.

Factors affecting the rate and pattern of the mutational process are being identified for human autosomes, but the same relationships for the male specific portion of the Y chromosome (MSY) are not established. We considered 3,390 mutations occurring in 19 sequence bins identified by sequencing 1.5 Mb of the MSY from each of 104 present-day chromosomes. The occurrence of mutations was not proportional to the amount of sequenced bases in each bin, with a 2-fold variation. The regression of the number of mutations per unit sequence against a number of indicators of the genomic features of each bin, revealed the same fundamental patterns as in the autosomes. By considering the sequences of the same region from two precisely dated ancient specimens, we obtained a calibrated region-specific substitution rate of 0.716 × 10-9/site/year. Despite its lack of recombination and other peculiar features, the MSY then resembles the autosomes in displaying a marked regional heterogeneity of the mutation rate. An immediate implication is that a given figure for the substitution rate only makes sense if bound to a specific DNA region. By strictly applying this principle we obtained an unbiased estimate of the antiquity of lineages relevant to the genetic history of the human Y chromosome. In particular, the two deepest nodes of the tree highlight the survival, in Central-Western Africa, of lineages whose coalescence (291 ky, 95% C.I. 253–343) predates the emergence of anatomically modern features in the fossil record.

One hundred and forty-six previously detected mutations were more precisely positioned in the human Y chromosome phylogeny by the analysis of 51 representative Y chromosome haplogroups and the use of 59 mutations from literature. Twenty-two new mutations were also described and incorporated in the revised phylogeny. This analysis made it possible to identify new haplogroups and to resolve a deep trifurcation within haplogroup B2. Our data provide a highly resolved branching in the African-specific portion of the Y tree and support the hypothesis of an origin in the north-western quadrant of the African continent for the human MSY diversity.

Haplogroup E, defined by mutation M40, is the most common human Y chromosome clade within Africa. To increase the level of resolution of haplogroup E, we disclosed the phylogenetic relationships among 729 mutations found in 33 haplogroup DE Y-chromosomes sequenced at high coverage in previous studies. Additionally, we dissected the E-M35 subclade by genotyping 62 informative markers in 5,222 samples from 118 worldwide populations. The phylogeny of haplogroup E showed novel features compared with the previous topology, including a new basal dichotomy. Within haplogroup E-M35, we resolved all the previously known polytomies and assigned all the E-M35* chromosomes to five new different clades, all belonging to a newly identified subhaplogroup (E-V1515), which accounts for almost half of the E-M35 chromosomes from the Horn of Africa. Moreover, using a Bayesian phylogeographic analysis and a single nucleotide polymorphism-based approach we localized and dated the origin of this new lineage in the northern part of the Horn, about 12ka. Time frames, phylogenetic structuring, and sociogeographic distribution of E-V1515 and its subclades are consistent with a multistep demic spread of pastoralism within north-eastern Africa and its subsequent diffusion to subequatorial areas. In addition, our results increase the discriminative power of the E-M35 haplogroup for use in forensic genetics through the identification of new ancestry-informative markers., The authors are grateful to all the donors for providing DNA samples and to the people that contributed to the sample collection. In particular, they thank João Lavinha (for the Portuguese samples); Farha El Chennawi, Anne Cambon-Thomsen, M.S. Issad, Eric Crubézy, Abdellatif Baali, Mohammed Cherkaoui, and Mohammed Melhaoui for their help in the collection of the Moroccan, Algerian, and Egyptian Berbers samples; and the National Laboratory for the Genetics of Israeli Populations. This work was supported by the Italian Ministry of Education, Progetti PRIN (grant numbers 2012JA4BTY_004 to F.C. and 2012JA4BTY_003 to A.N.) and by Sapienza University of Rome (grant number C26A13S9AR to F.C.). The sampling of the Berbers was made within the framework of the Inserm "Réseau Nord/Sud" No. 490NS1 (Mozabite Berbers), "The Origin of Man, Language and Languages," EUROCORES Programme and benefited from funding by the Région Midi-Pyrénées (Toulouse, France), the CNRS, and the E.C. Sixth Framework Programme under Contract ERASCT-2003-980409.