Search

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.

International audience; It has long been believed that the male-specific region of the human Y chromosome (MSY) is genetically independent from the X chromosome. This idea has been recently dismissed due to the discovery that X-Y gametologous gene conversion may occur. However, the pervasiveness of this molecular process in the evolution of sex chromosomes has yet to be exhaustively analyzed. In this study, we explored how pervasive X-Y gene conversion has been during the evolution of the youngest stratum of the human sex chromosomes. By comparing about 0.5 Mb of human-chimpanzee gametologous sequences, we identified 19 regions in which extensive gene conversion has occurred. From our analysis, two major features of these emerged: 1) Several of them are evolutionarily conserved between the two species and 2) almost all of the 19 hotspots overlap with regions where X-Y crossing-over has been previously reported to be involved in sex reversal. Furthermore, in order to explore the dynamics of X-Y gametologous conversion in recent human evolution, we resequenced these 19 hotspots in 68 widely divergent Y haplogroups and used publicly available single nucleotide polymorphism data for the X chromosome. We found that at least ten hotspots are still active in humans. Hence, the results of the interspecific analysis are consistent with the hypothesis of widespread reticulate evolution within gametologous sequences in the differentiation of hominini sex chromosomes. In turn, intraspecific analysis demonstrates that X-Y gene conversion may modulate human sex-chromosome-sequence evolution to a greater extent than previously thought.

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.

Forty-seven mtDNAs collected in the Dominican Republic and belonging to the African-specific haplogroup L2 were studied by high-resolution RFLP and control-region sequence analyses. Four sets of diagnostic markers that subdivide L2 into four clades (L2a–L2d) were identified, and a survey of published African data sets appears to indicate that these clades encompass all L2 mtDNAs and harbor very different geographic/ethnic distributions. One mtDNA from each of the four clades was completely sequenced by means of a new sequencing protocol that minimizes time and expense. The phylogeny of the L2 complete sequences showed that the two mtDNAs from L2b and L2d seem disproportionately derived, compared with those from L2a and L2c. This result is not consistent with a simple model of neutral evolution with a uniform molecular clock. The pattern of nonsynonymous versus synonymous substitutions hints at a role for selection in the evolution of human mtDNA. Regardless of whether selection is shaping the evolution of modern human mtDNAs, the population screening of L2 mtDNAs for the mutations identified by our complete sequence study should allow the identification of marker motifs of younger age with more restricted geographic distributions, thus providing new clues about African prehistory and the origin and relationships of African ethnic groups.

Founder analysis is a method for analysis of nonrecombining DNA sequence data, with the aim of identification and dating of migrations into new territory. The method picks out founder sequence types in potential source populations and dates lineage clusters deriving from them in the settlement zone of interest. Here, using mtDNA, we apply the approach to the colonization of Europe, to estimate the proportion of modern lineages whose ancestors arrived during each major phase of settlement. To estimate the Palaeolithic and Neolithic contributions to European mtDNA diversity more accurately than was previously achievable, we have now extended the Near Eastern, European, and northern-Caucasus databases to 1,234, 2, 804, and 208 samples, respectively. Both back-migration into the source population and recurrent mutation in the source and derived populations represent major obstacles to this approach. We have developed phylogenetic criteria to take account of both these factors, and we suggest a way to account for multiple dispersals of common sequence types. We conclude that (i) there has been substantial back-migration into the Near East, (ii) the majority of extant mtDNA lineages entered Europe in several waves during the Upper Palaeolithic, (iii) there was a founder effect or bottleneck associated with the Last Glacial Maximum, 20,000 years ago, from which derives the largest fraction of surviving lineages, and (iv) the immigrant Neolithic component is likely to comprise less than one-quarter of the mtDNA pool of modern Europeans.

Four X-linked loci showing homology with a previously described Y-linked polymorphic locus (DYS413) were identified and characterized. By fluorescent in situ hybridization (FISH), somatic cell hybrids, and YAC screening, the X-linked members of this small family of sequences (CAIII) all map in Xp22, while the Y members map in Yq11. These loci contribute to the overall similarity of the two genomic regions. All of the CAIII loci contain an internal microsatellite of the (CA)n type. The microsatellites display extensive length polymorphism in two of the X-linked members as well as in the Y members. In addition, common sequence variants are found in the portions flanking the microsatellites in two of the X-linked members. Our results indicate that, during the evolution of this family, length variation on the Y chromosome was accumulated at a rate not slower than that on the X chromosome. Finally, these sequences represent a model system with which to analyze human populations for similar X- and Y-linked polymorphisms.

To shed light on the structure of the basal backbone of the human Y chromosome phylogeny, we sequenced about 200 kb of the male-specific region of the human Y chromosome (MSY) from each of seven Y chromosomes belonging to clades A1, A2, A3, and BT. We detected 146 biallelic variant sites through this analysis. We used these variants to construct a patrilineal tree, without taking into account any previously reported information regarding the phylogenetic relationships among the seven Y chromosomes here analyzed. There are several key changes at the basal nodes as compared with the most recent reference Y chromosome tree. A different position of the root was determined, with important implications for the origin of human Y chromosome diversity. An estimate of 142 KY was obtained for the coalescence time of the revised MSY tree, which is earlier than that obtained in previous studies and easier to reconcile with plausible scenarios of modern human origin. The number of deep branchings leading to African-specific clades has doubled, further strengthening the MSY-based evidence for a modern human origin in the African continent. An analysis of 2204 African DNA samples showed that the deepest clades of the revised MSY phylogeny are currently found in central and northwest Africa, opening new perspectives on early human presence in the continent.

Haplogroup E1b1, defined by the marker P2, is the most represented human Y chromosome haplogroup in Africa. A phylogenetic tree showing the internal structure of this haplogroup was published in 2008. A high degree of internal diversity characterizes this haplogroup, as well as the presence of a set of chromosomes undefined on the basis of a derived character. Here we make an effort to update the phylogeny of this highly diverse haplogroup by including seven mutations which have been newly discovered by direct resequencing. We also try to incorporate five previously-described markers which were not, however, reported in the 2008 tree. Additionally, during the process of mapping, we found that two previously reported SNPs required a new position on the tree. There are three key changes compared to the 2008 phylogeny. Firstly, haplogroup E-M2 (former E1b1a) and haplogroup E-M329 (former E1b1c) are now united by the mutations V38 and V100, reducing the number of E1b1 basal branches to two. The new topology of the tree has important implications concerning the origin of haplogroup E1b1. Secondly, within E1b1b1 (E-M35), two haplogroups (E-V68 and E-V257) show similar phylogenetic and geographic structure, pointing to a genetic bridge between southern European and northern African Y chromosomes. Thirdly, most of the E1b1b1* (E-M35*) paragroup chromosomes are now marked by defining mutations, thus increasing the discriminative power of the haplogroup for use in human evolution and forensics.

More than 2700 unrelated individuals from Europe, northern Africa and western Asia were analyzed for the marker M269, which defines the Y chromosome haplogroup R1b1b2. A total of 593 subjects belonging to this haplogroup were identified and further analyzed for two SNPs, U106 and U152, which define haplogroups R1b1b2g and R1b1b2h, respectively. These haplogroups showed quite different frequency distribution patterns within Europe, with frequency peaks in northern Europe (R1b1b2g) and northern Italy/France (R1b1b2h).

Although human Y chromosomes belonging to haplogroup R1b are quite rare in Africa, being found mainly in Asia and Europe, a group of chromosomes within the paragroup R-P25(*) are found concentrated in the central-western part of the African continent, where they can be detected at frequencies as high as 95%. Phylogenetic evidence and coalescence time estimates suggest that R-P25(*) chromosomes (or their phylogenetic ancestor) may have been carried to Africa by an Asia-to-Africa back migration in prehistoric times. Here, we describe six new mutations that define the relationships among the African R-P25(*) Y chromosomes and between these African chromosomes and earlier reported R-P25 Eurasian sub-lineages. The incorporation of these new mutations into a phylogeny of the R1b haplogroup led to the identification of a new clade (R1b1a or R-V88) encompassing all the African R-P25(*) and about half of the few European/west Asian R-P25(*) chromosomes. A worldwide phylogeographic analysis of the R1b haplogroup provided strong support to the Asia-to-Africa back-migration hypothesis. The analysis of the distribution of the R-V88 haplogroup in1800 males from 69 African populations revealed a striking genetic contiguity between the Chadic-speaking peoples from the central Sahel and several other Afroasiatic-speaking groups from North Africa. The R-V88 coalescence time was estimated at 9.2-5.6 [corrected] kya, in the early mid Holocene. We suggest that R-V88 is a paternal genetic record of the proposed mid-Holocene migration of proto-Chadic Afroasiatic speakers through the Central Sahara into the Lake Chad Basin, and geomorphological evidence is consistent with this view.

Different X-homologous regions of the male-specific portion of the human Y chromosome (MSY) are characterized by a different content of putative single nucleotide polymorphisms (SNPs), as reported in public databases. The possible role of X-to-Y nonallelic gene conversion in contributing to these differences remains poorly understood. We explored this issue by analyzing sequence variation in three regions of the MSY characterized by a different degree of X―Y similarity and a different density of putative SNPs: the PCDH11Y gene in the X-transposed (X―Y identity 99%, high putative SNP content); the TBL1Y gene in the X-degenerate (X―Y identity 86―88%, low putative SNP content); and VCY genes―containing region in the P8 palindrome (X―Y identity 95%, low putative SNP content). Present findings do not provide any evidence for gene conversion in the PCDH11Y and TBL1Y genes; they also strongly suggest that most putative SNPs of the PCDH11Y gene (and possibly the entire X-transposed region) are most likely X―Y paralogous sequence variants, which have been entered in the databases as SNPs. On the other hand, clear evidence for the VCY genes in the P8 palindrome having acted as an acceptor of X-to-Y gene conversion was obtained. A rate of 1.8 x 10 ―7 X-to-Y conversions/ bp/year was estimated for these genes. These findings indicate that in the VCY region of the MSY, X-to-Y gene conversion can be highly effective to increase the level of diversity among human Y chromosomes and suggest an additional explanation for the ability of the Y chromosome to retard degradation during evolution. Present data are expected to pave the way for future investigations on the role of nonallelic gene conversion in double-strand break repair and the maintenance of Y chromosome integrity.

Thirty-five mitochondrial (mt) DNAs from Spain that harbor the mutation A3243G in association with either MELAS (mitochondrial myopathy, encephalopathy, lactic acidosis, and strokelike episodes) syndrome or a wide array of disease phenotypes (ranging from diabetes and deafness to a mixture of chronic progressive external ophthalmoplegic symptoms and strokelike episodes) were studied by use of high-resolution restriction fragment length polymorphism analysis and control-region sequencing. A total of 34 different haplotypes were found, indicating that all instances of the A3243G mutation are probably due to independent mutational events. Haplotypes were distributed into 13 haplogroups whose frequencies were close to those of the general Spanish population. Moreover, there was no statistically significant difference in haplogroup distribution between patients with MELAS and those with disease phenotypes other than MELAS. Overall, these data indicate that the A3243G mutation harbors all the evolutionary features expected from a severely deleterious mtDNA mutation under strong negative selection, and they reveal that European mtDNA backgrounds do not play a substantial role in modulating the mutation's phenotypic expression.

Mitochondrial HVS-I sequences from 10, 365 subjects belonging to 56 populations/geographical regions of western Eurasia and northern Africa were first surveyed for the presence of the T-*C transition at nucleotide position 16298, a mutation which has previously been shown to characterize haplogroup V mtDNAs. All mtDNAs with this mutation were then screened for a number of diagnostic RFLP sites, revealing two major subsets of mtDNAs. One is haplogroup V proper, and the other has been termed "pre~V, " since it predates V phylogenetically. The rather uncommon pre*V tends to be scattered throughout Europe (and northwestern Africa), whereas V attains two peaks of frequency: one situated in southwestern Europe and one in the Saami of northern Scandinavia. Geographical distributions and ages support the scenario that pre*V originated in Europe before the Last Glacial Maximum (LGM), whereas the more recently derived haplogroup V arose in a southwestern European refugium soon after the LGM. The arrival of V in eastern/central Europe, however, occurred much later, possibly with (post-)Neohthic contacts. The distribution of haplogroup V mtDNAs in modern European populations would thus, at least in part, reflect the pattern of postglacial human recolonization from that refugium, affecting even the Saami. Overall, the present study shows that the dissection of mtDNA variation into small and well-defined evolutionary units is an essential step in the identification of spatial frequency patterns. Mass screening of a few markers identified using complete mtDNA sequences promises to be an efficient strategy for inferring features of human prehistory.

In January 2010, we published in this journal a report1 on the frequency distribution of the Y chromosome haplogroup R1b1a (R-V88) in Africa, where it can be found at frequencies as high as about 90%. This haplogroup (or its ancestor) most likely traces its origins back to Eurasia, but is presently found very rarely outside Africa.

SummaryTo define Y-chromosome haplotypes, we studied seven biallelic polymorphic sites. We combined data with those from four dinucleotide-repeat polymorphisms, to establish Y-chromosome compound superhaplotypes. Eight biallelic haplotypes that matched the dendrogram proposed by other investigators were identified in 762 Y chromosomes from 25 African populations. For each biallelic site, coalescence time of lineages carrying the derived allele was estimated and compared with previous estimates. The “ancestral” haplotype (haplotype 1A) was observed among Ethiopians, “Khoisan” (!Kung and Khwe), and populations from northern Cameroon. Microsatellite distributions within this haplotype showed that the Khoisan haplotypes 1A are widely divergent from those of the other two groups. Populations from northern Africa and northern Cameroon share a haplotype (i.e., 1C), which is not observed in other African populations but represents a major Eurasian cluster. Haplotypes 1C of northern Cameroon are clearly distinct from those of Europe, whereas haplotypes 1C of northern African are well intermingled with those of the other two groups. Apportionment of diversity for the Y-chromosomal biallelic haplotypes was calculated after populations were clustered into different configurations. Despite some correspondence between language affiliation and genetic similarity, geographic proximity seems to be a better predictor of genetic affinity.

mtDNA sequence variation was studied in 419 individuals from nine Eurasian populations, by high-resolution RFLP analysis, and it was followed by sequencing of the control region of a subset of these mtDNAs and a detailed survey of previously published data from numerous other European populations. This analysis revealed that a major Paleolithic population expansion from the "Atlantic zone" (southwestern Europe) occurred 10,000-15,000 years ago, after the Last Glacial Maximum. As an mtDNA marker for this expansion we identified haplogroup V, an autochthonous European haplogroup, which most likely originated in the northern Iberian peninsula or southwestern France at about the time of the Younger Dryas. Its sister haplogroup, H, which is distributed throughout the entire range of Caucasoid populations and which originated in the Near East approximately 25,000-30,000 years ago, also took part in this expansion, thus rendering it by far the most frequent (40%-60%) haplogroup in western Europe. Subsequent migrations after the Younger Dryas eventually carried those "Atlantic" mtDNAs into central and northern Europe. This scenario, already implied by archaeological records, is given overwhelming support from both the distribution of the autochthonous European Y chromosome type 15, as detected by the probes 49a/f, and the synthetic maps of nuclear data.

A polymorphic MspI site was located in the X-specific region proximal to the pseudoautosomal boundary. Two alleles defined by the absence or the presence of the MspI site were detected by the polymerase chain reaction (PCR) in a sample of Bantu-speaking individuals from Cameroon. No variation for this site was observed in a population sample from Italy.

Correction to: Eur J Hum Genet advance online publication, 6 January 2010; doi:10.1038/ejhg.2009.231 Since the publication of the above paper, the authors have identified two errors: (1) In the Abstract, on line 12, ‘9200–5600 kya’ should be ‘9.2–5.6 kya’. (2) On page 7, first column, line 37, ‘7000kya’ should be ‘7.

Summary The mtDNA variation of 50 Spanish and 4 Cuban families affected by nonsyndromic sensorineural deafness due to the A1555G mutation in the 12S rRNA gene was studied by high-resolution RFLP analysis and sequencing of the control region. Phylogenetic analyses of haplotypes and detailed survey of population controls revealed that the A1555G mutation can be attributed to ⩾30 independent mutational events among the 50 Spanish families and that it occurs on mtDNA haplogroups that are common in all European populations. This indicates that the relatively high detection rate of this mutation in Spain is not due to sampling biases or to a single major founder event. Moreover, the distribution of these mutational events on different haplogroups is compatible with a random occurrence of the A1555G mutation and tends to support the conclusion that mtDNA backgrounds do not play a significant role in the expression of the mutation. Overall, these findings appear to indicate that the rare detection of this mutation in other populations is most likely due to inadequacy in patient ascertainment and molecular screening. This probable lack of identification of the A1555G mutation in subjects affected by sensorineural hearing loss implies that their maternally related relatives are not benefiting from presymptomatic detection and information concerning their increased risk of ototoxicity due to aminoglycoside treatments.

The electrophoretic and thermostability polymorphisms of the PGM1 locus were examined in about 700 Czechoslovakians (Prague) and 3000 Italians. The Italian sample consisted of individuals from Pavia (Northern Italy), Viareggio and Rome (Central Italy) and Naples (Southern Italy). The eight PGM1 alleles, PGM1Str1, PGM1Sts1, PGM1Ftr1, PGM1Fts1, PGM2Str1, PGM2Sts1, PGM2Ftr1, PGM2Fts1, have been considered as combinations of mutations at three different sites, 1/2, S/F and tr/ts, within the PGM1 gene and their frequencies discussed in terms of linkage disequilibrium between these sites. All pairwise differences between the samples were significant except for Pavia-Viareggio and Viareggio-Rome. The frequencies of the PGMts1 alleles have been found to range from 0.0981 (Prague) to 0.0546 (Naples) and can be ordered according to a North-South cline.

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.