Reconstruction of patrilineages and matrilineages of Samaritans and other Israeli populations from Y-Chromosome and mitochondrial DNA sequence Variation

African mitochondrial haplogroup L7: a 100,000-year-old maternal human lineage discovered through reassessment and new sequencing

Archaeological and genomic evidence suggest that modern Homo sapiens have roamed the planet for some 300–500 thousand years. In contrast, global human mitochondrial (mtDNA) diversity coalesces to one African female ancestor (“Mitochondrial Eve”) some 145 thousand years ago, owing to the ¼ gene pool size of our matrilineally inherited haploid genome. Therefore, most of human prehistory was spent in Africa where early ancestors of Southern African Khoisan and Central African rainforest hunter-gatherers (RFHGs) segregated into smaller groups. Their subdivisions followed climatic oscillations, new modes of subsistence, local adaptations, and cultural-linguistic differences, all prior to their exodus out of Africa. Seven African mtDNA haplogroups (L0–L6) traditionally captured this ancient structure—these L haplogroups have formed the backbone of the mtDNA tree for nearly two decades. Here we describe L7, an eighth haplogroup that we estimate to be ~ 100 thousand years old and which has been previously misclassified in the literature. In addition, L7 has a phylogenetic sublineage L7a*, the oldest singleton branch in the human mtDNA tree (~ 80 thousand years). We found that L7 and its sister group L5 are both low-frequency relics centered around East Africa, but in different populations (L7: Sandawe; L5: Mbuti). Although three small subclades of African foragers hint at the population origins of L5'7, the majority of subclades are divided into Afro-Asiatic and eastern Bantu groups, indicative of more recent admixture. A regular re-estimation of the entire mtDNA haplotype tree is needed to ensure correct cladistic placement of new samples in the future.

National Biobanking in Qatar and Israel: Tracing how Global Scientific Institutions Mediate Local Ethnic Identities

Biobanks are a growing phenomenon in global biomedicine, as they are key tools of precision medicine initiatives. National biobanks, however, collect data and biological material from populations in specific regions, and the knowledge that national biobanks yield can impact understandings of identity, origins and belonging. Drawing on ethnographic work and documentary analysis examining the Israeli and Qatari national biobanks, I find that these two Middle Eastern biobanks aim to contribute to global biobanking trends, while at the same time, they reinforce local ethnic and national identities. The Israeli biobank reflects pre-existing ethnic identities in Israeli society, while the Qatari biobank predominantly emphasises the emergent national character of the Qatari population. Neither of the biobanks assert a high genetic homogeneity of the national population; rather, they both emphasise a genetically diverse national cohort that is a valuable resource for biomedical research. Through a comparative analysis of global biobanking and ethnic identities, this article demonstrates that biobanks are a rich site for tracking emergent national identities in the Middle East region.

Middle eastern genetic legacy in the paternal and maternal gene pools of Chuetas

Chuetas are a group of descendants of Majorcan Crypto-Jews (Balearic Islands, Spain) who were socially stigmatized and segregated by their Majorcan neighbours until recently; generating a community that, although after the seventeenth century no longer contained Judaic religious elements, maintained strong group cohesion, Jewishness consciousness, and endogamy. Collective memory fixed 15 surnames as a most important defining element of Chueta families. Previous studies demonstrated Chuetas were a differentiated population, with a considerable proportion of their original genetic make-up. Genetic data of Y-chromosome polymorphism and mtDNA control region showed, in Chuetas’ paternal lineages, high prevalence of haplogroups J2-M172 (33%) and J1-M267 (18%). In maternal lineages, the Chuetas hallmark is the presence of a new sub-branching of the rare haplogroup R0a2m as their modal haplogroup (21%). Genetic diversity in both Y-chromosome and mtDNA indicates the Chueta community has managed to avoid the expected heterogeneity decrease in their gene pool after centuries of isolation and inbreeding. Moreover, the composition of their uniparentally transmitted lineages demonstrates a remarkable signature of Middle Eastern ancestry—despite some degree of host admixture—confirming Chuetas have retained over the centuries a considerable degree of ancestral genetic signature along with the cultural memory of their Jewish origin.

Paternal portrait of populations of the middle Magdalena River region (Tolima and Huila, Colombia): New insights on the peopling of Central America and northernmost South America

The valley of the Magdalena River is one of the main population pathways in Colombia. The gene pool and spatial configuration of human groups in this territory have been outlined throughout three historical stages: the Native pre-Hispanic world, Spanish colonization, and XIX century migrations. This research was designed with the goal of characterizing the diversity and distribution pattern of Y-chromosome lineages that are currently present in the Tolima and Huila departments (middle Magdalena River region). Historic cartography was used to identify the main geographic sites where the paternal lineages belonging to this area have gathered. Twelve municipalities were chosen, and a survey that included genealogical information was administered. Samples collected from 83 male volunteers were analyzed for 48 Y-SNPs and 17 Y-STRs. The results showed a highly diverse region characterized by the presence of 16 sublineages within the major clades R, Q, J, G, T and E and revealed that 93% (n = 77) of haplotypes were different. Among these haplogroups, European-specific R1b-M269 lineages were the most representative (57.83%), with six different subhaplogroups and 43 unique haplotypes. Native American paternal ancestry was also detected based on the presence of the Q1a2-M3*(xM19, M194, M199) and Q1a2-M346*(xM3) lineages. Interestingly, all Q1a2-M346*(xM3) samples (n = 7, with five different haplotypes) carried allele six at the DYS391 locus. This allele has a worldwide frequency of 0.169% and was recently associated with a new Native subhaplogroup. An in-depth phylogenetic analysis of these samples suggests the Tolima and Huila region to be the principal area in all Central and South America where this particular Native lineage is found. This lineage has been present in the region for at least 1,809 (+/- 0,5345) years.

The peopling of the last Green Sahara revealed by high-coverage resequencing of trans-Saharan patrilineages

BACKGROUND

Little is known about the peopling of the Sahara during the Holocene climatic optimum, when the desert was replaced by a fertile environment.

RESULTS

In order to investigate the role of the last Green Sahara in the peopling of Africa, we deep-sequence the whole non-repetitive portion of the Y chromosome in 104 males selected as representative of haplogroups which are currently found to the north and to the south of the Sahara. We identify 5,966 mutations, from which we extract 142 informative markers then genotyped in about 8,000 subjects from 145 African, Eurasian and African American populations. We find that the coalescence age of the trans-Saharan haplogroups dates back to the last Green Sahara, while most northern African or sub-Saharan clades expanded locally in the subsequent arid phase.

CONCLUSIONS

Our findings suggest that the Green Sahara promoted human movements and demographic expansions, possibly linked to the adoption of pastoralism. Comparing our results with previously reported genome-wide data, we also find evidence for a sex-biased sub-Saharan contribution to northern Africans, suggesting that historical events such as the trans-Saharan slave trade mainly contributed to the mtDNA and autosomal gene pool, whereas the northern African paternal gene pool was mainly shaped by more ancient events.

Reconstructing Druze population history

The Druze are an aggregate of communities in the Levant and Near East living almost exclusively in the mountains of Syria, Lebanon and Israel whose ~1000 year old religion formally opposes mixed marriages and conversions. Despite increasing interest in genetics of the population structure of the Druze, their population history remains unknown. We investigated the genetic relationships between Israeli Druze and both modern and ancient populations. We evaluated our findings in light of three hypotheses purporting to explain Druze history that posit Arabian, Persian or mixed Near Eastern-Levantine roots. The biogeographical analysis localised proto-Druze to the mountainous regions of southeastern Turkey, northern Iraq and southeast Syria and their descendants clustered along a trajectory between these two regions. The mixed Near Eastern–Middle Eastern localisation of the Druze, shown using both modern and ancient DNA data, is distinct from that of neighbouring Syrians, Palestinians and most of the Lebanese, who exhibit a high affinity to the Levant. Druze biogeographic affinity, migration patterns, time of emergence and genetic similarity to Near Eastern populations are highly suggestive of Armenian-Turkish ancestries for the proto-Druze.

New native South American Y chromosome lineages

Many single-nucleotide polymorphisms (SNPs) in the non-recombining region of the human Y chromosome have been described in the last decade. High-coverage sequencing has helped to characterize new SNPs, which has in turn increased the level of detail in paternal phylogenies. However, these paternal lineages still provide insufficient information on population history and demography, especially for Native Americans. The present study aimed to identify informative paternal sublineages derived from the main founder lineage of the Americas-haplogroup Q-L54-in a sample of 1841 native South Americans. For this purpose, we used a Y-chromosomal genotyping multiplex platform and conventional genotyping methods to validate 34 new SNPs that were identified in the present study by sequencing, together with many Y-SNPs previously described in the literature. We updated the haplogroup Q phylogeny and identified two new Q-M3 and three new Q-L54*(xM3) sublineages defined by five informative SNPs, designated SA04, SA05, SA02, SA03 and SA29. Within the Q-M3, sublineage Q-SA04 was mostly found in individuals from ethnic groups belonging to the Tukanoan linguistic family in the northwest Amazon, whereas sublineage Q-SA05 was found in Peruvian and Bolivian Amazon ethnic groups. Within Q-L54*, the derived sublineages Q-SA03 and Q-SA02 were exclusively found among Coyaima individuals (Cariban linguistic family) from Colombia, while Q-SA29 was found only in Maxacali individuals (Jean linguistic family) from southeast Brazil. Furthermore, we validated the usefulness of several published SNPs among indigenous South Americans. This new Y chromosome haplogroup Q phylogeny offers an informative paternal genealogy to investigate the pre-Columbian history of South America.Journal of Human Genetics advance online publication, 31 March 2016; doi:10.1038/jhg.2016.26.

Y-chromosome phylogeographic analysis of the Greek-Cypriot population reveals elements consistent with Neolithic and Bronze Age settlements

Background

The archeological record indicates that the permanent settlement of Cyprus began with pioneering agriculturalists circa 11,000 years before present, (ca. 11,000 y BP). Subsequent colonization events followed, some recognized regionally. Here, we assess the Y-chromosome structure of Cyprus in context to regional populations and correlate it to phases of prehistoric colonization.

Results

Analysis of haplotypes from 574 samples showed that island-wide substructure was barely significant in a spatial analysis of molecular variance (SAMOVA). However, analyses of molecular variance (AMOVA) of haplogroups using 92 binary markers genotyped in 629 Cypriots revealed that the proportion of variance among the districts was irregularly distributed. Principal component analysis (PCA) revealed potential genetic associations of Greek-Cypriots with neighbor populations. Contrasting haplogroups in the PCA were used as surrogates of parental populations. Admixture analyses suggested that the majority of G2a-P15 and R1b-M269 components were contributed by Anatolia and Levant sources, respectively, while Greece Balkans supplied the majority of E-V13 and J2a-M67. Haplotype-based expansion times were at historical levels suggestive of recent demography.

Conclusions

Analyses of Cypriot haplogroup data are consistent with two stages of prehistoric settlement. E-V13 and E-M34 are widespread, and PCA suggests sourcing them to the Balkans and Levant/Anatolia, respectively. The persistent pre-Greek component is represented by elements of G2-U5(xL30) haplogroups: U5*, PF3147, and L293. J2b-M205 may contribute also to the pre-Greek strata. The majority of R1b-Z2105 lineages occur in both the westernmost and easternmost districts. Distinctively, sub-haplogroup R1b- M589 occurs only in the east. The absence of R1b- M589 lineages in Crete and the Balkans and the presence in Asia Minor are compatible with Late Bronze Age influences from Anatolia rather than from Mycenaean Greeks.

Electronic supplementary material

The online version of this article (doi:10.1186/s13323-016-0032-8) contains supplementary material, which is available to authorized users.

Portuguese crypto-Jews: the genetic heritage of a complex history

The first documents mentioning Jewish people in Iberia are from the Visigothic period. It was also in this period that the first documented anti-Judaic persecution took place. Other episodes of persecution would happen again and again during the long troubled history of the Jewish people in Iberia and culminated with the Decrees of Expulsion and the establishment of the Inquisition: some Jews converted to Catholicism while others resisted and were forcedly baptized, becoming the first Iberian Crypto-Jews. In the 18th century the official discrimination and persecution carried out by the Inquisition ended and several Jewish communities emerged in Portugal. From a populational genetics point of view, the worldwide Diaspora of contemporary Jewish communities has been intensely studied. Nevertheless, very little information is available concerning Sephardic and Iberian Crypto-Jewish descendants. Data from the Iberian Peninsula, the original geographic source of Sephardic Jews, is limited to two populations in Portugal, Belmonte, and Bragança district, and the Chueta community from Mallorca. Belmonte was the first Jewish community studied for uniparental markers. The construction of a reference model for the history of the Portuguese Jewish communities, in which the genetic and classical historical data interplay dynamically, is still ongoing. Recently an enlarged sample covering a wide region in the Northeast Portugal was undertaken, allowing the genetic profiling of male and female lineages. A Jewish specific shared female lineage (HV0b) was detected between the community of Belmonte and Bragança. In contrast to what was previously described as a hallmark of the Portuguese Jews, an unexpectedly high polymorphism of lineages was found in Bragança, showing a surprising resistance to the erosion of genetic diversity typical of small-sized isolate populations, as well as signs of admixture with the Portuguese host population.

Genetics and the history of the Samaritans: Y-chromosomal microsatellites and genetic affinity between Samaritans and Cohanim

The Samaritans are a group of some 750 indigenous Middle Eastern people, about half of whom live in Holon, a suburb of Tel Aviv, and the other half near Nablus. The Samaritan population is believed to have numbered more than a million in late Roman times but less than 150 in 1917. The ancestry of the Samaritans has been subject to controversy from late Biblical times to the present. In this study, liquid chromatography/electrospray ionization/quadrupole ion trap mass spectrometry was used to allelotype 13 Y-chromosomal and 15 autosomal microsatellites in a sample of 12 Samaritans chosen to have as low a level of relationship as possible, and 461 Jews and non-Jews. Estimation of genetic distances between the Samaritans and seven Jewish and three non-Jewish populations from Israel, as well as populations from Africa, Pakistan, Turkey, and Europe, revealed that the Samaritans were closely related to Cohanim. This result supports the position of the Samaritans that they are descendants from the tribes of Israel dating to before the Assyrian exile in 722-720 BCE. In concordance with previously published single-nucleotide polymorphism haplotypes, each Samaritan family, with the exception of the Samaritan Cohen lineage, was observed to carry a distinctive Y-chromosome short tandem repeat haplotype that was not more than one mutation removed from the six-marker Cohen modal haplotype.

Afghanistan from a Y-chromosome perspective

Central Asia has served as a corridor for human migrations providing trading routes since ancient times. It has functioned as a conduit connecting Europe and the Middle East with South Asia and far Eastern civilizations. Therefore, the study of populations in this region is essential for a comprehensive understanding of early human dispersal on the Eurasian continent. Although Y- chromosome distributions in Central Asia have been widely surveyed, present-day Afghanistan remains poorly characterized genetically. The present study addresses this lacuna by analyzing 190 Pathan males from Afghanistan using high-resolution Y-chromosome binary markers. In addition, haplotype diversity for its most common lineages (haplogroups R1a1a*-M198 and L3-M357) was estimated using a set of 15 Y-specific STR loci. The observed haplogroup distribution suggests some degree of genetic isolation of the northern population, likely due to the Hindu Kush mountain range separating it from the southern Afghans who have had greater contact with neighboring Pathans from Pakistan and migrations from the Indian subcontinent. Our study demonstrates genetic similarities between Pathans from Afghanistan and Pakistan, both of which are characterized by the predominance of haplogroup R1a1a*-M198 (>50%) and the sharing of the same modal haplotype. Furthermore, the high frequencies of R1a1a-M198 and the presence of G2c-M377 chromosomes in Pathans might represent phylogenetic signals from Khazars, a common link between Pathans and Ashkenazi groups, whereas the absence of E1b1b1a2-V13 lineage does not support their professed Greek ancestry.

Paternal genetic history of the Basque population of Spain

This study examines the genetic variation in Basque Y chromosome lineages using data on 12 Y-short tandem repeat (STR) loci in a sample of 158 males from four Basque provinces of Spain (Alava, Vizcaya, Guipuzcoa, and Navarre). As reported in previous studies, the Basques are characterized by high frequencies of haplogroup R1b (83%). AMOVA analysis demonstrates genetic homogeneity, with a small but significant amount of genetic structure between provinces (Y-short tandem repeat loci STRs: 1.71%, p = 0.0369). Gene and haplotype diversity levels in the Basque population are on the low end of the European distribution (gene diversity: 0.4268; haplotype diversity: 0.9421). Post-Neolithic contribution to the paternal Basque gene pool was estimated by measuring the proportion of those haplogroups with a Time to Most Recent Common Ancestor (TMRCA) previously dated either prior (R1b, I2a2) or subsequent to (E1b1b, G2a, J2a) the Neolithic. Based on these estimates, the Basque provinces show varying degrees of post-Neolithic contribution in the paternal lineages (10.9% in the combined sample).

Mitochondrial DNA reveals distinct evolutionary histories for Jewish populations in Yemen and Ethiopia

Southern Arabia and the Horn of Africa are important geographic centers for the study of human population history because a great deal of migration has characterized these regions since the first emergence of humans out of Africa. Analysis of Jewish groups provides a unique opportunity to investigate more recent population histories in this area. Mitochondrial DNA is used to investigate the maternal evolutionary history and can be combined with historical and linguistic data to test various population histories. In this study, we assay mitochondrial control region DNA sequence and diagnostic coding variants in Yemenite (n = 45) and Ethiopian (n = 41) Jewish populations, as well as in neighboring non-Jewish Yemeni (n = 50) and Ethiopian (previously published Semitic speakers) populations. We investigate their population histories through a comparison of haplogroup distributions and phylogenetic networks. A high frequency of sub-Saharan African L haplogroups was found in both Jewish populations, indicating a significant African maternal contribution unlike other Jewish Diaspora populations. However, no identical haplotypes were shared between the Yemenite and Ethiopian Jewish populations, suggesting very little gene flow between the populations and potentially distinct maternal population histories. These new data are also used to investigate alternate population histories in the context of historical and linguistic data. Specifically, Yemenite Jewish mitochondrial diversity reflects potential descent from ancient Israeli exiles and shared African and Middle Eastern ancestry with little evidence for large-scale conversion of local Yemeni. In contrast, the Ethiopian Jewish population appears to be a subset of the larger Ethiopian population suggesting descent primarily through conversion of local women.

Alcohol-metabolizing enzyme gene polymorphisms in the Basque Country, Morocco, and Ecuador

BACKGROUND

Genes ADH1B and ADH1C have certain functional SNPs that are related to alcoholism. The frequencies of these polymorphisms vary between populations, so studying them in populations made up of groups with different phylogeographic origins requires an individualized analysis of each group. In the Basque Country, various recently arrived foreign groups live side by side with the original Southern European population, particularly North Africans from Morocco and Hispanics from Ecuador. This study sets out to examine the distribution of the frequencies of alleles that encode alcohol dehydrogenase with different metabolization rates, as higher rates make for greater susceptibility to alcoholism.

METHODS

Four SNPs: rs1229984, rs2066702, rs698, and rs1693482 using Taqman technology with a Rt-PCR were studied in a sample of 114 European individuals originating from the Basque Country, 100 North Africans from Morocco, and 109 Hispanics from Ecuador. The allele and genotype frequencies were calculated using Genepop v4.0. The most frequent haplotypes were estimated using the ELB algorithm with Arlequin v3.01. A breakdown of the complete disequilibrium commonly observed between the 2 missense polymorphisms that distinguish the common ADH1C alleles rs698 and rs1693482 was observed and confirmed by sequencing in 2 individuals from the Basque Country.

RESULTS

A higher frequency of protective allele ADH1C*1 was found in the North African population group. Haplotype combinations are also studied, and the rare association of alleles ADH1B*2-ADH1C*2 was observed in the Southern European group in the Basque Country, along with an allele not hitherto described in the ADH1C locus.

CONCLUSIONS

This study provides the first data published on the allele and genotype frequencies of the ADH1C locus in the Moroccan population and on the ADH1B and ADH1C loci in the Ecuadorian population. The study shows differences in the distribution of the frequency of allele ADH1C*1 between the Basque Country and Moroccan populations, and a new allele not described to date.

The origin of Eastern European Jews revealed by autosomal, sex chromosomal and mtDNA polymorphisms

Background

This study aims to establish the likely origin of EEJ (Eastern European Jews) by genetic distance analysis of autosomal markers and haplogroups on the X and Y chromosomes and mtDNA.

Results

According to the autosomal polymorphisms the investigated Jewish populations do not share a common origin, and EEJ are closer to Italians in particular and to Europeans in general than to the other Jewish populations. The similarity of EEJ to Italians and Europeans is also supported by the X chromosomal haplogroups. In contrast according to the Y-chromosomal haplogroups EEJ are closest to the non-Jewish populations of the Eastern Mediterranean. MtDNA shows a mixed pattern, but overall EEJ are more distant from most populations and hold a marginal rather than a central position. The autosomal genetic distance matrix has a very high correlation (0.789) with geography, whereas the X-chromosomal, Y-chromosomal and mtDNA matrices have a lower correlation (0.540, 0.395 and 0.641 respectively).

Conclusions

The close genetic resemblance to Italians accords with the historical presumption that Ashkenazi Jews started their migrations across Europe in Italy and with historical evidence that conversion to Judaism was common in ancient Rome. The reasons for the discrepancy between the biparental markers and the uniparental markers are discussed.

Reviewers

This article was reviewed by Damian Labuda (nominated by Jerzy Jurka), Kateryna Makova and Qasim Ayub (nominated by Dan Graur).

Identification of population substructure among Jews using STR markers and dependence on reference populations included

BACKGROUND

Detecting population substructure is a critical issue for association studies of health behaviors and other traits. Whether inherent in the population or an artifact of marker choice, determining aspects of a population's genetic history as potential sources of substructure can aid in design of future genetic studies. Jewish populations, among which association studies are often conducted, have a known history of migrations. As a necessary step in understanding population structure to conduct valid association studies of health behaviors among Israeli Jews, we investigated genetic signatures of this history and quantified substructure to facilitate future investigations of these phenotypes in this population.

RESULTS

Using 32 autosomal STR markers and the program STRUCTURE, we differentiated between Ashkenazi (AJ, N = 135) and non-Ashkenazi (NAJ, N = 226) Jewish populations in the form of Northern and Southern geographic genetic components (AJ north 73%, south 23%, NAJ north 33%, south 60%). The ability to detect substructure within these closely related populations using a small STR panel was contingent on including additional samples representing major continental populations in the analyses.

CONCLUSIONS

Although clustering programs such as STRUCTURE are designed to assign proportions of ancestry to individuals without reference population information, when Jewish samples were analyzed in the absence of proxy parental populations, substructure within Jews was not detected. Generally, for samples with a given grandparental country of birth, STRUCTURE assignment values to Northern, Southern, African and Asian clusters agreed with mitochondrial DNA and Y-chromosomal data from previous studies as well as historical records of migration and intermarriage.

Extended Y chromosome haplotypes resolve multiple and unique lineages of the Jewish priesthood

It has been known for over a decade that a majority of men who self report as members of the Jewish priesthood (Cohanim) carry a characteristic Y chromosome haplotype termed the Cohen Modal Haplotype (CMH). The CMH has since been used to trace putative Jewish ancestral origins of various populations. However, the limited number of binary and STR Y chromosome markers used previously did not provide the phylogenetic resolution needed to infer the number of independent paternal lineages that are encompassed within the Cohanim or their coalescence times. Accordingly, we have genotyped 75 binary markers and 12 Y-STRs in a sample of 215 Cohanim from diverse Jewish communities, 1,575 Jewish men from across the range of the Jewish Diaspora, and 2,099 non-Jewish men from the Near East, Europe, Central Asia, and India. While Cohanim from diverse backgrounds carry a total of 21 Y chromosome haplogroups, 5 haplogroups account for 79.5% of Cohanim Y chromosomes. The most frequent Cohanim lineage (46.1%) is marked by the recently reported P58 T->C mutation, which is prevalent in the Near East. Based on genotypes at 12 Y-STRs, we identify an extended CMH on the J-P58* background that predominates in both Ashkenazi and non-Ashkenazi Cohanim and is remarkably absent in non-Jews. The estimated divergence time of this lineage based on 17 STRs is 3,190 ± 1,090 years. Notably, the second most frequent Cohanim lineage (J-M410*, 14.4%) contains an extended modal haplotype that is also limited to Ashkenazi and non-Ashkenazi Cohanim and is estimated to be 4.2 ± 1.3 ky old. These results support the hypothesis of a common origin of the CMH in the Near East well before the dispersion of the Jewish people into separate communities, and indicate that the majority of contemporary Jewish priests descend from a limited number of paternal lineages.

The complex and diversified mitochondrial gene pool of Berber populations

The mitochondrial DNA variation of 295 Berber-speakers from Morocco (Asni, Bouhria and Figuig) and the Egyptian oasis of Siwa was evaluated by sequencing a portion of the control region (including HVS-I and part of HVS-II) and surveying haplogroup-specific coding region markers. Our findings show that the Berber mitochondrial pool is characterized by an overall high frequency of Western Eurasian haplogroups, a somehow lower frequency of sub-Saharan L lineages, and a significant (but differential) presence of North African haplogroups U6 and M1, thus occupying an intermediate position between European and sub-Saharan populations in PCA analysis. A clear and significant genetic differentiation between the Berbers from Maghreb and Egyptian Berbers was also observed. The first are related to European populations as shown by haplogroup H1 and V frequencies, whereas the latter share more affinities with East African and Nile Valley populations as indicated by the high frequency of M1 and the presence of L0a1, L3i, L4*, and L4b2 lineages. Moreover, haplogroup U6 was not observed in Siwa. We conclude that the origins and maternal diversity of Berber populations are old and complex, and these communities bear genetic characteristics resulting from various events of gene flow with surrounding and migrating populations.

Mitochondrial DNA variation in Jordanians and their genetic relationship to other Middle East populations

BACKGROUND

The Levant is a crucial region in understanding human migrations between Africa and Eurasia. Although some mitochondrial DNA (mtDNA) studies have been carried out in this region, they have not included the Jordan area. This paper deals with the mtDNA composition of two Jordan populations.

AIM

The main objectives of this article are: first, to report mtDNA sequences of an urban and an isolate sample from Jordan and, second, to compare them with each other and with other nearby populations.

SUBJECTS AND METHODS

The analyses are based on HVSI and HVSII mtDNA sequences and diagnostic RFLPs to unequivocally classify into haplogroups 101 Amman and 44 Dead Sea unrelated individuals from Jordan.

RESULTS

Statistical analysis revealed that, whereas the sample from Amman did not significantly differ from their Levantine neighbours, the Dead Sea sample clearly behaved as a genetic outlier in the region. Its outstanding Eurasian haplogroup U3 frequency (39%) and its south-Saharan Africa lineages (19%) are the highest in the Middle East. On the contrary, the lack ((preHV)1) or comparatively low frequency (J and T) of Neolithic lineages is also striking. Although strong drift by geographic isolation could explain the anomalous mtDNA pool of the Dead Sea sample, the fact that its mtDNA lineage composition mirrors, in geographic origin and haplogroup frequencies, its Y-chromosome pool, points to founder effect as the main cause. Ancestral M1 lineages detected in Jordan that have affinities with those recently found in Northwest but not East Africa question the African origin of the M1 haplogroup.

CONCLUSION

Results are in agreement with an old human settlement in the Jordan region. However, in spite of the attested migratory spreads, genetically divergent populations, such as that of the Dead Sea, still exist in the area.

Counting the founders: the matrilineal genetic ancestry of the Jewish Diaspora

The history of the Jewish Diaspora dates back to the Assyrian and Babylonian conquests in the Levant, followed by complex demographic and migratory trajectories over the ensuing millennia which pose a serious challenge to unraveling population genetic patterns. Here we ask whether phylogenetic analysis, based on highly resolved mitochondrial DNA (mtDNA) phylogenies can discern among maternal ancestries of the Diaspora. Accordingly, 1,142 samples from 14 different non-Ashkenazi Jewish communities were analyzed. A list of complete mtDNA sequences was established for all variants present at high frequency in the communities studied, along with high-resolution genotyping of all samples. Unlike the previously reported pattern observed among Ashkenazi Jews, the numerically major portion of the non-Ashkenazi Jews, currently estimated at 5 million people and comprised of the Moroccan, Iraqi, Iranian and Iberian Exile Jewish communities showed no evidence for a narrow founder effect, which did however characterize the smaller and more remote Belmonte, Indian and the two Caucasus communities. The Indian and Ethiopian Jewish sample sets suggested local female introgression, while mtDNAs in all other communities studied belong to a well-characterized West Eurasian pool of maternal lineages. Absence of sub-Saharan African mtDNA lineages among the North African Jewish communities suggests negligible or low level of admixture with females of the host populations among whom the African haplogroup (Hg) L0-L3 sub-clades variants are common. In contrast, the North African and Iberian Exile Jewish communities show influence of putative Iberian admixture as documented by mtDNA Hg HV0 variants. These findings highlight striking differences in the demographic history of the widespread Jewish Diaspora.

Differential Y-chromosome Anatolian influences on the Greek and Cretan Neolithic

The earliest Neolithic sites of Europe are located in Crete and mainland Greece. A debate persists concerning whether these farmers originated in neighboring Anatolia and the role of maritime colonization. To address these issues 171 samples were collected from areas near three known early Neolithic settlements in Greece together with 193 samples from Crete. An analysis of Y-chromosome haplogroups determined that the samples from the Greek Neolithic sites showed strong affinity to Balkan data, while Crete shows affinity with central/Mediterranean Anatolia. Haplogroup J2b-M12 was frequent in Thessaly and Greek Macedonia while haplogroup J2a-M410 was scarce. Alternatively, Crete, like Anatolia showed a high frequency of J2a-M410 and a low frequency of J2b-M12. This dichotomy parallels archaeobotanical evidence, specifically that while bread wheat (Triticum aestivum) is known from Neolithic Anatolia, Crete and southern Italy; it is absent from earliest Neolithic Greece. The expansion time of YSTR variation for haplogroup E3b1a2-V13, in the Peloponnese was consistent with an indigenous Mesolithic presence. In turn, two distinctive haplogroups, J2a1h-M319 and J2a1b1-M92, have demographic properties consistent with Bronze Age expansions in Crete, arguably from NW/W Anatolia and Syro-Palestine, while a later mainland (Mycenaean) contribution to Crete is indicated by relative frequencies of V13.

Mitochondrial DNA structure in the Arabian Peninsula

BACKGROUND

Two potential migratory routes followed by modern humans to colonize Eurasia from Africa have been proposed. These are the two natural passageways that connect both continents: the northern route through the Sinai Peninsula and the southern route across the Bab al Mandab strait. Recent archaeological and genetic evidence have favored a unique southern coastal route. Under this scenario, the study of the population genetic structure of the Arabian Peninsula, the first step out of Africa, to search for primary genetic links between Africa and Eurasia, is crucial. The haploid and maternally inherited mitochondrial DNA (mtDNA) molecule has been the most used genetic marker to identify and to relate lineages with clear geographic origins, as the African Ls and the Eurasian M and N that have a common root with the Africans L3.

RESULTS

To assess the role of the Arabian Peninsula in the southern route, we genetically analyzed 553 Saudi Arabs using partial (546) and complete mtDNA (7) sequencing, and compared the lineages obtained with those present in Africa, the Near East, central, east and southeast Asia and Australasia. The results showed that the Arabian Peninsula has received substantial gene flow from Africa (20%), detected by the presence of L, M1 and U6 lineages; that an 18% of the Arabian Peninsula lineages have a clear eastern provenance, mainly represented by U lineages; but also by Indian M lineages and rare M links with Central Asia, Indonesia and even Australia. However, the bulk (62%) of the Arabian lineages has a Northern source.

CONCLUSION

Although there is evidence of Neolithic and more recent expansions in the Arabian Peninsula, mainly detected by (preHV)1 and J1b lineages, the lack of primitive autochthonous M and N sequences, suggests that this area has been more a receptor of human migrations, including historic ones, from Africa, India, Indonesia and even Australia, than a demographic expansion center along the proposed southern coastal route.

Y-chromosome diversity characterizes the Gulf of Oman

Arabia has served as a strategic crossroads for human disseminations, providing a natural connection between the distant populations of China and India in the east to the western civilizations along the Mediterranean. To explore this region's critical role in the migratory episodes leaving Africa to Eurasia and back, high-resolution Y-chromosome analysis of males from the United Arab Emirates (164), Qatar (72) and Yemen (62) was performed. The role of the Levant in the Neolithic dispersal of the E3b1-M35 sublineages is supported by the data, and the distribution and STR-based analyses of J1-M267 representatives points to their spread from the north, most likely during the Neolithic. With the exception of Yemen, southern Arabia, South Iran and South Pakistan display high diversity in their Y-haplogroup substructure possibly a result of gene flow along the coastal crescent-shaped corridor of the Gulf of Oman facilitating human dispersals. Elevated rates of consanguinity may have had an impact in Yemen and Qatar, which experience significant heterozygote deficiencies at various hypervariable autosomal STR loci.

The phylogeography of African Brazilians

BACKGROUND/AIMS

Approximately four million Africans were taken as slaves to Brazil, where they interbred extensively with Amerindians and Europeans. We have previously shown that while most White Brazilians carry Y chromosomes of European origin, they display high proportions of African and Amerindian mtDNA lineages, because of sex-biased genetic admixture.

METHODS

We studied the Y chromosome and mtDNA haplogroup structure of 120 Black males from Sao Paulo, Brazil.

RESULTS

Only 48% of the Y chromosomes, but 85% of the mtDNA haplogroups were characteristic of sub-Saharan Africa, confirming our previous observation of sexually biased mating. We mined literature data for mtDNA and Y chromosome haplogroup frequencies for African native populations from regions involved in Atlantic Slave Trade. Principal Components Analysis and Bayesian analysis of population structure revealed no genetic differentiation of Y chromosome marker frequencies between the African regions. However, mtDNA examination unraveled considerable genetic structure, with three clusters at Central-West Africa, West Africa and Southeast Africa. A hypothesis is proposed to explain this structure.

CONCLUSION

Using these mtDNA data we could obtain for the first time an estimate of the relative ancestral contribution of Central-West (0.445), West (0.431) and Southeast Africa (0.123) to African Brazilians from Sao Paulo. These estimates are consistent with historical information.

Mitochondrial lineage M1 traces an early human backflow to Africa

Background

The out of Africa hypothesis has gained generalized consensus. However, many specific questions remain unsettled. To know whether the two M and N macrohaplogroups that colonized Eurasia were already present in Africa before the exit is puzzling. It has been proposed that the east African clade M1 supports a single origin of haplogroup M in Africa. To test the validity of that hypothesis, the phylogeographic analysis of 13 complete mitochondrial DNA (mtDNA) sequences and 261 partial sequences belonging to haplogroup M1 was carried out.

Results

The coalescence age of the African haplogroup M1 is younger than those for other M Asiatic clades. In contradiction to the hypothesis of an eastern Africa origin for modern human expansions out of Africa, the most ancestral M1 lineages have been found in Northwest Africa and in the Near East, instead of in East Africa. The M1 geographic distribution and the relative ages of its different subclades clearly correlate with those of haplogroup U6, for which an Eurasian ancestor has been demonstrated.

Conclusion

This study provides evidence that M1, or its ancestor, had an Asiatic origin. The earliest M1 expansion into Africa occurred in northwestern instead of eastern areas; this early spread reached the Iberian Peninsula even affecting the Basques. The majority of the M1a lineages found outside and inside Africa had a more recent eastern Africa origin. Both western and eastern M1 lineages participated in the Neolithic colonization of the Sahara. The striking parallelism between subclade ages and geographic distribution of M1 and its North African U6 counterpart strongly reinforces this scenario. Finally, a relevant fraction of M1a lineages present today in the European Continent and nearby islands possibly had a Jewish instead of the commonly proposed Arab/Berber maternal ascendance.

Molecular characterization of a polymorphic 3-Mb deletion at chromosome Yp11.2 containing the AMELY locus in Singapore and Malaysia populations

Amelogenin paralogs on Chromosome X (AMELX) and Y (AMELY) are commonly used sexing markers. Interstitial deletion of Yp involving the AMELY locus has previously been reported. The combined frequency of the AMELY null allele in Singapore and Malaysia populations is 2.7%, 0.6% in Indian and Malay ethnic groups respectively. It is absent among 541 Chinese screened. The null allele in this study belongs to 3 Y haplogroups; J2e1 (85.7%), F* (9.5%) and D* (4.8%). Low and high-resolution STS mapping, followed by sequence analysis of breakpoint junction confirmed a large deletion of 3 to 3.7-Mb located at the Yp11.2 region. Both breakpoints were located in TSPY repeat arrays, suggesting a non-allelic homologous recombination (NAHR) mechanism of deletion. All regional null samples shared identical breakpoint sequences according to their haplogroup affiliation, providing molecular evidence of a common ancestry origin for each haplogroup, and at least 3 independent deletion events recurred in history. The estimated ages based on Y-SNP and STR analysis were approximately 13.5 +/- 3.1 kyears and approximately 0.9 +/- 0.9 kyears for the J2e1 and F* mutations, respectively. A novel polymorphism G > A at Y-GATA-H4 locus in complete linkage disequilibrium with J2e1 null mutations is a more recent event. This work re-emphasizes the need to include other sexing markers for gender determination in certain regional populations. The frequency difference among global populations suggests it constitutes another structural variation locus of human chromosome Y. The breakpoint sequences provide further information to a better understanding of the NAHR mechanism and DNA rearrangements due to higher order genomic architecture.

Thomas Jefferson's Y chromosome belongs to a rare European lineage

We have characterized the Y chromosome carried by President Thomas Jefferson, the general rarity of which supported the idea that he, or a patrilineal relative, fathered the last son of his slave Sally Hemings. It belongs to haplogroup K2, a lineage representing only approximately 1% of chromosomes worldwide, and most common in East Africa and the Middle East. Phylogenetic network analysis of its Y-STR (short tandem repeat) haplotype shows that it is most closely related to an Egyptian K2 haplotype, but the presence of scattered and diverse European haplotypes within the network is nonetheless consistent with Jefferson's patrilineage belonging to an ancient and rare indigenous European type. This is supported by the observation that two of 85 unrelated British men sharing the surname Jefferson also share the President's Y-STR haplotype within haplogroup K2. Our findings represent a cautionary tale in showing the difficulty of assigning individual ancestry based on a Y-chromosome haplotype, particularly for rare lineages where population data are scarce.

Paleolithic Y-haplogroup heritage predominates in a Cretan highland plateau

The island of Crete, credited by some historical scholars as a central crucible of western civilization, has been under continuous archeological investigation since the second half of the nineteenth century. In the present work, the geographic stratification of the contemporary Cretan Y-chromosome gene pool was assessed by high-resolution haplotyping to investigate the potential imprints of past colonization episodes and the population substructure. In addition to analyzing the possible geographic origins of Y-chromosome lineages in relatively accessible areas of the island, this study includes samples from the isolated interior of the Lasithi Plateau--a mountain plain located in eastern Crete. The potential significance of the results from the latter region is underscored by the possibility that this region was used as a Minoan refugium. Comparisons of Y-haplogroup frequencies among three Cretan populations as well as with published data from additional Mediterranean locations revealed significant differences in the frequency distributions of Y-chromosome haplogroups within the island. The most outstanding differences were observed in haplogroups J2 and R1, with the predominance of haplogroup R lineages in the Lasithi Plateau and of haplogroup J lineages in the more accessible regions of the island. Y-STR-based analyses demonstrated the close affinity that R1a1 chromosomes from the Lasithi Plateau shared with those from the Balkans, but not with those from lowland eastern Crete. In contrast, Cretan R1b microsatellite-defined haplotypes displayed more resemblance to those from Northeast Italy than to those from Turkey and the Balkans.

Niger-Congo speaking populations and the formation of the Brazilian gene pool: mtDNA and Y-chromosome data

We analyzed sequence variation in the mitochondrial DNA (mtDNA) hypervariable segment I (HVS-I) from 201 Black individuals from two Brazilian cities (Rio de Janeiro and Porto Alegre), and compared these data with published information from 21 African populations. A subset of 187 males of the sample was also characterized for 30 Y-chromosome biallelic polymorphisms, and the data were compared with those from 48 African populations. The mtDNA data indicated that respectively 69% and 82% of the matrilineages found in Rio de Janeiro and Porto Alegre originated from West-Central/Southeast Africa. These estimates are in close agreement with historical records which indicated that most of the Brazilian slaves who arrived in Rio de Janeiro were from West-Central Africa. In contrast to mtDNA, Y-chromosome haplogroup analysis did not allow discrimination between places of origin in West or West-Central Africa. Thus, when comparing these two major African regions, there seems to be higher genetic structure with mtDNA than with Y-chromosome data.

Whole-mtDNA Genome Sequence Analysis of Ancient African Lineages

Studies of human mitochondrial (mt) DNA genomes demonstrate that the root of the human phylogenetic tree occurs in Africa. Although 2 mtDNA lineages with an African origin (haplogroups M and N) were the progenitors of all non-African haplogroups, macrohaplogroup L (including haplogroups L0-L6) is limited to sub-Saharan Africa. Several L haplogroup lineages occur most frequently in eastern Africa (e.g., L0a, L0f, L5, and L3g), but some are specific to certain ethnic groups, such as haplogroup lineages L0d and L0k that previously have been found nearly exclusively among southern African "click" speakers. Few studies have included multiple mtDNA genome samples belonging to haplogroups that occur in eastern and southern Africa but are rare or absent elsewhere. This lack of sampling in eastern Africa makes it difficult to infer relationships among mtDNA haplogroups or to examine events that occurred early in human history. We sequenced 62 complete mtDNA genomes of ethnically diverse Tanzanians, southern African Khoisan speakers, and Bakola Pygmies and compared them with a global pool of 226 mtDNA genomes. From these, we infer phylogenetic relationships amongst mtDNA haplogroups and estimate the time to most recent common ancestor (TMRCA) for haplogroup lineages. These data suggest that Tanzanians have high genetic diversity and possess ancient mtDNA haplogroups, some of which are either rare (L0d and L5) or absent (L0f) in other regions of Africa. We propose that a large and diverse human population has persisted in eastern Africa and that eastern Africa may have been an ancient source of dispersion of modern humans both within and outside of Africa.

Iran: tricontinental nexus for Y-chromosome driven migration

Due to its pivotal geographic position, present day Iran likely served as a gateway of reciprocal human movements. However, the extent to which the deserts within the Iranian plateau and the mountain ranges surrounding Persia inhibited gene flow via this corridor remains uncertain. In order to assess the magnitude of this region's role as a nexus for Africa, Asia and Europe in human migrations, high-resolution Y-chromosome analyses were performed on 150 Iranian males. Haplogroup data were subsequently compared to regional populations characterized at similar phylogenetic levels. The Iranians display considerable haplogroup diversity consistent with patterns observed in populations of the Middle East overall, reinforcing the notion of Persia as a venue for human disseminations. Admixture analyses of geographically targeted, regional populations along the latitudinal corridor spanning from Anatolia to the Indus Valley demonstrated contributions to Persia from both the east and west. However, significant differences were uncovered upon stratification of the gene donors, including higher proportions from central east and southeast Turkey as compared to Pakistan. In addition to the modulating effects of geographic obstacles, culturally mediated amalgamations consistent with the diverse spectrum of a variety of historical empires may account for the distribution of haplogroups and lineages observed. Our study of high-resolution Y-chromosome genotyping allowed for an in-depth analysis unattained in previous studies of the area, revealing important migratory and demographic events that shaped the contemporary genetic landscape.

A distinct Y-STR haplotype for Amelogenin negative males characterized by a large Y(p)11.2 (DYS458-MSY1-AMEL-Y) deletion

The use of STR multiplexes with the incorporated gender marker Amelogenin is common practice in forensic DNA analysis. However, when a known male sample shows a dropout of the Amelogenin Y-allele, the STR system falsely genotypes it as a female. To date, our laboratory has observed 18 such cases: 12 from our Y-STR database and six from casework. A study on 980 male individuals in the Malaysian population using the AmpFlSTR Y-filer has revealed a distinct Y-chromosome haplotype associated with the Amelogenin nulls. Our results showed that whilst the Amelogenin nulls were noticeably absent among the Chinese, both the Indians and Malays exhibited such mutations at 3.2 and 0.6%, respectively. It was also found that the Amelogenin negative individuals predominantly belonged to the J2e lineage, suggesting the possibility of a common ancestor for at least some of these chromosomes. The null frequencies showed concordance with the data published in Chang et al. [Higher failures of Amelogenin sex test in an Indian population group, J. Forensic Sci. 48 (2003) 1309-1313] on a smaller Malaysian population of 338 males which used a Y-STR triplex. In the current study, apart from the absence of the Amelogenin Y-locus, a complete absence of the DYS458 locus in all the nulls was also observed. This study together with the 2003 study has indicated a similar deletion region exists on the Y(p)11.2 band in all the 18 Y-chromosomes. Using bioinformatics, this deletion has been mapped to a region of at least 1.13 Mb on the Y(p)11.2 encompassing the Amelogenin, MSY1 minisatellite and DYS458 locus. Further, the Y-filer haplotypes revealed an additional null at Y-GATA H4 in two of the Indian males presented here.

The role of selection in the evolution of human mitochondrial genomes

High mutation rate in mammalian mitochondrial DNA generates a highly divergent pool of alleles even within species that have dispersed and expanded in size recently. Phylogenetic analysis of 277 human mitochondrial genomes revealed a significant (P < 0.01) excess of rRNA and nonsynonymous base substitutions among hotspots of recurrent mutation. Most hotspots involved transitions from guanine to adenine that, with thymine-to-cytosine transitions, illustrate the asymmetric bias in codon usage at synonymous sites on the heavy-strand DNA. The mitochondrion-encoded tRNAThr varied significantly more than any other tRNA gene. Threonine and valine codons were involved in 259 of the 414 amino acid replacements observed. The ratio of nonsynonymous changes from and to threonine and valine differed significantly (P = 0.003) between populations with neutral (22/58) and populations with significantly negative Tajima's D values (70/76), independent of their geographic location. In contrast to a recent suggestion that the excess of nonsilent mutations is characteristic of Arctic populations, implying their role in cold adaptation, we demonstrate that the surplus of nonsynonymous mutations is a general feature of the young branches of the phylogenetic tree, affecting also those that are found only in Africa. We introduce a new calibration method of the mutation rate of synonymous transitions to estimate the coalescent times of mtDNA haplogroups.

Isolates in a corridor of migrations: a high-resolution analysis of Y-chromosome variation in Jordan

A high-resolution, Y-chromosome analysis using 46 binary markers has been carried out in two Jordan populations, one from the metropolitan area of Amman and the other from the Dead Sea, an area geographically isolated. Comparisons with neighboring populations showed that whereas the sample from Amman did not significantly differ from their Levantine neighbors, the Dead Sea sample clearly behaved as a genetic outlier in the region. Its high R1*-M173 frequency (40%) has until now only been found in northern Cameroonian samples. This contrasts with the comparatively low presence of J representatives (9%), which is the modal clade in Middle Eastern populations, including Amman. The Dead Sea sample also showed a high presence of E3b3a-M34 lineages (31%), which is only comparable to that found in Ethiopians. Although ancient and recent ties with sub-Saharan and eastern Africans cannot be discarded, it seems that isolation, strong drift, and/or founder effects are responsible for the anomalous Y-chromosome pool of this population. These results demonstrate that, at a fine scale, the smooth, continental clines detected for several Y-chromosome markers are often disrupted by genetically divergent populations.

MITOCHONDRIAL DNA AND HUMAN EVOLUTION

Several unique properties of human mitochondrial DNA (mtDNA), including its high copy number, maternal inheritance, lack of recombination, and high mutation rate, have made it the molecule of choice for studies of human population history and evolution. Here we review the current state of knowledge concerning these properties, how mtDNA variation is studied, what we have learned, and what the future likely holds. We conclude that increasingly, mtDNA studies are (and should be) supplemented with analyses of the Y-chromosome and other nuclear DNA variation. Some serious issues need to be addressed concerning nuclear inserts, database quality, and the possible influence of selection on mtDNA variation. Nonetheless, mtDNA studies will continue to play an important role in such areas as examining socio-cultural influences on human genetic variation, ancient DNA, certain forensic DNA applications, and in tracing personal genetic history.

Ethiopian mitochondrial DNA heritage: tracking gene flow across and around the gate of tears

Approximately 10 miles separate the Horn of Africa from the Arabian Peninsula at Bab-el-Mandeb (the Gate of Tears). Both historic and archaeological evidence indicate tight cultural connections, over millennia, between these two regions. High-resolution phylogenetic analysis of 270 Ethiopian and 115 Yemeni mitochondrial DNAs was performed in a worldwide context, to explore gene flow across the Red and Arabian Seas. Nine distinct subclades, including three newly defined ones, were found to characterize entirely the variation of Ethiopian and Yemeni L3 lineages. Both Ethiopians and Yemenis contain an almost-equal proportion of Eurasian-specific M and N and African-specific lineages and therefore cluster together in a multidimensional scaling plot between Near Eastern and sub-Saharan African populations. Phylogeographic identification of potential founder haplotypes revealed that approximately one-half of haplogroup L0-L5 lineages in Yemenis have close or matching counterparts in southeastern Africans, compared with a minor share in Ethiopians. Newly defined clade L6, the most frequent haplogroup in Yemenis, showed no close matches among 3,000 African samples. These results highlight the complexity of Ethiopian and Yemeni genetic heritage and are consistent with the introduction of maternal lineages into the South Arabian gene pool from different source populations of East Africa. A high proportion of Ethiopian lineages, significantly more abundant in the northeast of that country, trace their western Eurasian origin in haplogroup N through assorted gene flow at different times and involving different source populations.