Multiple origins of Ashkenazi Levites: Y chromosome evidence for both Near Eastern and European ancestries

Previous Y chromosome studies have shown that the Cohanim, a paternally inherited Jewish priestly caste, predominantly share a recent common ancestry irrespective of the geographically defined post-Diaspora community to which they belong, a finding consistent with common Jewish origins in the Near East. In contrast, the Levites, another paternally inherited Jewish caste, display evidence for multiple recent origins, with Ashkenazi Levites having a high frequency of a distinctive, non-Near Eastern haplogroup. Here, we show that the Ashkenazi Levite microsatellite haplotypes within this haplogroup are extremely tightly clustered, with an inferred common ancestor within the past 2,000 years. Comparisons with other Jewish and non-Jewish groups suggest that a founding event, probably involving one or very few European men occurring at a time close to the initial formation and settlement of the Ashkenazi community, is the most likely explanation for the presence of this distinctive haplogroup found today in >50% of Ashkenazi Levites.     

Association of Y chromosome haplogroup I with HIV progression, and HAART outcome

The host genetic basis of differential outcomes in HIV infection, progression, viral load set point and highly active retroviral therapy (HAART) responses was examined for the common Y haplogroups in European Americans and African Americans. Accelerated progression to acquired immune deficiency syndrome (AIDS) and related death in European Americans among Y chromosome haplogroup I (Y-I) subjects was discovered. Additionally, Y-I haplogroup subjects on HAART took a longer time to HIV-1 viral suppression and were more likely to fail HAART. Both the accelerated progression and longer time to viral suppression results observed in haplogroup Y-I were significant after false-discovery-rate corrections. A higher frequency of AIDS-defining illnesses was also observed in haplogroup Y-I. These effects were independent of the previously identified autosomal AIDS restriction genes. When the Y-I haplogroup subjects were further subdivided into six I subhaplogroups, no one subhaplogroup accounted for the effects on HIV progression, viral load or HAART response. Adjustment of the analyses for population stratification found significant and concordant haplogroup Y-I results. The Y chromosome haplogroup analyses of HIV infection and progression in African Americans were not significant. Our results suggest that one or more loci on the Y chromosome found on haplogroup Y-I have an effect on AIDS progression and treatment responses in European Americans. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Exploring the legacy of African and Indigenous Caribbean admixture in Puerto Rico

Objectives

From an anthropological genetic perspective, little is known about the ethnogenesis of African descendants in Puerto Rico. Furthermore, historical interactions between Indigenous Caribbean and African descendant peoples that may be reflected in the ancestry of contemporary populations are understudied. Given this dearth of genetic research and the precedence for Afro-Indigenous interactions documented by historical, archeological, and other lines of evidence, we sought to assess the biogeographic origins of African descendant Puerto Ricans and to query the potential for Indigenous ancestry within this community.

Materials and Methods

Saliva samples were collected from 58 self-identified African descendant Puerto Ricans residing in Puerto Rico. We sequenced whole mitochondrial genomes and genotyped Y chromosome haplogroups for each male individual (n = 25). Summary statistics, comparative analyses, and network analysis were used to assess diversity and variation in haplogroup distribution between the sample and comparative populations.

Results

As indicated by mitochondrial haplogroups, 66% had African, 5% had European, and 29% had Indigenous American matrilines. Along the Y chromosome, 52% had African, 28% had Western European, 16% had Eurasian, and, notably, 4% had Indigenous American patrilines. Both mitochondrial and Y chromosome haplogroup frequencies were significantly different from several comparative populations.

Discussion

Biogeographic origins are consistent with historical accounts of African, Indigenous American, and European ancestry. However, this first report of Indigenous American paternal ancestry in Puerto Rico suggests distinctive features within African descendant communities on the island. Future studies expanding sampling and incorporating higher resolution genetic markers are necessary to more fully understand African descendant history in Puerto Rico.     

Machine-learning approaches for classifying haplogroup from Y chromosome STR data

Genetic variation on the non-recombining portion of the Y chromosome contains information about the ancestry of male lineages. Because of their low rate of mutation, single nucleotide polymorphisms (SNPs) are the markers of choice for unambiguously classifying Y chromosomes into related sets of lineages known as haplogroups, which tend to show geographic structure in many parts of the world. However, performing the large number of SNP genotyping tests needed to properly infer haplogroup status is expensive and time consuming. A novel alternative for assigning a sampled Y chromosome to a haplogroup is presented here. We show that by applying modern machine-learning algorithms we can infer with high accuracy the proper Y chromosome haplogroup of a sample by scoring a relatively small number of Y-linked short tandem repeats (STRs). Learning is based on a diverse ground-truth data set comprising pairs of SNP test results (haplogroup) and corresponding STR scores. We apply several independent machine-learning methods in tandem to learn formal classification functions. The result is an integrated high-throughput analysis system that automatically classifies large numbers of samples into haplogroups in a cost-effective and accurate manner.     

Increased resolution of Y chromosome haplogroup T defines relationships among populations of the Near East, Europe, and Africa

Increasing phylogenetic resolution of the Y chromosome haplogroup tree has led to finer temporal and spatial resolution for studies of human migration. Haplogroup T, initially known as K2 and defined by mutation M70, is found at variable frequencies across West Asia, Africa, and Europe. While several SNPs were recently discovered that extended the length of the branch leading to haplogroup T, only two SNPs are known to mark internal branches of haplogroup T. This low level of phylogenetic resolution has hindered studies of the origin and dispersal of this interesting haplogroup, which is found in Near Eastern non-Jewish populations, Jewish populations from several communities, and in the patrilineage of President Thomas Jefferson. Here we map 10 new SNPs that, together with the previously known SNPs, mark 11 lineages and two large subclades (T1a and T1b) of haplogroup T. We also report a new SNP that links haplogroups T and L within the major framework of Y chromosome evolution. Estimates of the timing of the branching events within haplogroup T, along with a comprehensive geographic survey of the major T subclades, suggest that this haplogroup began to diversify in the Near East -25 kya. Our survey also points to a complex history of dispersal of this rare and informative haplogroup within the Near East and from the Near East to Europe and sub-Saharan Africa. The presence of T1a2 chromosomes in Near Eastern Jewish and non-Jewish populations may reflect early exiles between the ancient lands of Israel and Babylon. The presence of different subclades of T chromosomes in Europe may be explained by both the spread of Neolithic farmers and the later dispersal of Jews from the Near East. Finally, the moderately high frequency (-18%) of T1b* chromosomes in the Lemba of southern Africa supports the hypothesis of a Near Eastern, but not necessarily a Jewish, origin for their paternal line.     

Y-SNP haplogroups related to the Yqh+ heteromorphism in the Mexican northwestern population

Morphological variation of the Y chromosome has been observed in different populations. This variation is mostly related to the heteromorphic Yq12 band, which is composed of a variable block of constitutive heterochromatin. The Yqh+ heteromorphism has a worldwide frequency of 2.85% and is considered clinically innocuous. The aim of this study was to identify the ancestry of the Yqh+ heteromorphism present in individuals from western Mexico. For this purpose, 17 Y-chromosome single nucleotide polymorphisms were analysed by multiplex polymerase chain reaction and SNaPshot assays. In 28 Yqh+ males, only five haplogroups were observed; with a haplogroup diversity of 0.4841 ± 0.1094, which was less than that observed in a study of unselected Mexican mestizo population. Differences were specifically conferred by the high frequencies of haplogroups R1b1 and P*(xQ,R), and by the absence of the Amerindian haplogroup Q (Q*(xQ1a3a) plus Q1a3a) from the Yqh+ group. This study suggests a post-1492 incorporation for Yqh+ chromosomes into the Mexican northwestern population.     

A genome-wide genetic signature of Jewish ancestry perfectly separates individuals with and without full Jewish ancestry in a large random sample of European Americans

Background  

It was recently shown that the genetic distinction between self-identified Ashkenazi Jewish and non-Jewish individuals is a prominent component of genome-wide patterns of genetic variation in European Americans. No study however has yet assessed how accurately self-identified (Ashkenazi) Jewish ancestry can be inferred from genomic information, nor whether the degree of Jewish ancestry can be inferred among individuals with fewer than four Jewish grandparents.     

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).     

Contrasting patterns of Y chromosome variation in Ashkenazi Jewish and host non-Jewish European populations

The molecular basis of more than 25 genetic diseases has been described in Ashkenazi Jewish populations. Most of these diseases are characterized by one or two major founder mutations that are present in the Ashkenazi population at elevated frequencies. One explanation for this preponderance of recessive diseases is accentuated genetic drift resulting from a series of dispersals to and within Europe, endogamy, and/or recent rapid population growth. However, a clear picture of the manner in which neutral genetic variation has been affected by such a demographic history has not yet emerged. We have examined a set of 32 binary markers (single nucleotide polymorphisms; SNPs) and 10 microsatellites on the non-recombining portion of the Y chromosome (NRY) to investigate the ways in which patterns of variation differ between Ashkenazi Jewish and their non-Jewish host populations in Europe. This set of SNPs defines a total of 20 NRY haplogroups in these populations, at least four of which are likely to have been part of the ancestral Ashkenazi gene pool in the Near East, and at least three of which may have introgressed to some degree into Ashkenazi populations after their dispersal to Europe. It is striking that whereas Ashkenazi populations are genetically more diverse at both the SNP and STR level compared with their European non-Jewish counterparts, they have greatly reduced within-haplogroup STR variability, especially in those founder haplogroups that migrated from the Near East. This contrasting pattern of diversity in Ashkenazi populations is evidence for a reduction in male effective population size, possibly resulting from a series of founder events and high rates of endogamy within Europe. This reduced effective population size may explain the high incidence of founder disease mutations despite overall high levels of NRY diversity.Electronic Supplementary Material Supplementary material is available in the online version of this article at D.M. Behar and D. Garrigan contributed equally to this workElectronic database information: URLs for the data in this article are as follows:ARLEQUIN,     

Breast cancer risk and 6q22.33: combined results from Breast Cancer Association Consortium and Consortium of Investigators on Modifiers of BRCA1/2

Recently, a locus on chromosome 6q22.33 (rs2180341) was reported to be associated with increased breast cancer risk in the Ashkenazi Jewish (AJ) population, and this association was also observed in populations of non-AJ European ancestry. In the present study, we performed a large replication analysis of rs2180341 using data from 31,428 invasive breast cancer cases and 34,700 controls collected from 25 studies in the Breast Cancer Association Consortium (BCAC). In addition, we evaluated whether rs2180341 modifies breast cancer risk in 3,361 BRCA1 and 2,020 BRCA2 carriers from 11 centers in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). Based on the BCAC data from women of European ancestry, we found evidence for a weak association with breast cancer risk for rs2180341 (per-allele odds ratio (OR)?=?1.03, 95% CI 1.00-1.06, p?=?0.023). There was evidence for heterogeneity in the ORs among studies (I(2)?=?49.3%; p?=?<0.004). In CIMBA, we observed an inverse association with the minor allele of rs2180341 and breast cancer risk in BRCA1 mutation carriers (per-allele OR?=?0.89, 95%CI 0.80-1.00, p?=?0.048), indicating a potential protective effect of this allele. These data suggest that that 6q22.33 confers a weak effect on breast cancer risk.     

Mitochondrial Haplogroups and Control Region Polymorphisms Are Not Associated with Prostate Cancer in Middle European Caucasians

Background

Besides being responsible for energy production in the cell, mitochondria are central players in apoptosis as well as the main source of harmful reactive oxygen species. Therefore, it can be hypothesised that sequence variation in the mitochondrial genome is a contributing factor to the etiology of diseases related to these different cellular events, including cancer. The aim of the present study was to assess the frequency of haplogroups and polymorphisms in the control region (CR) of mitochondrial DNA of peripheral blood mononuclear cells from patients with prostate carcinoma (n = 304) versus patients screened for prostate disease but found to be negative for cancer on biopsy (n = 278) in a Middle European population.

Methodology/Principal Findings

The nine major European haplogroups and the CR polymorphisms were identified by means of primer extension analysis and DNA sequencing, respectively. We found that mitochondrial haplogroup frequencies and CR polymorphisms do not differ significantly between patients with or without prostate cancer, implying no impact of inherited mitochondrial DNA variation on predisposition to prostate carcinoma in a Middle European population.

Conclusions/Significance

Our results contrast with a recent report claiming an association between mtDNA haplogroup U and prostate cancer in a North American population of caucasian descent.     

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.     

Characterization of mitochondrial haplogroups in a large population-based sample from the United States

Mitochondrial DNA (mtDNA) haplogroups are valuable for investigations in forensic science, molecular anthropology, and human genetics. In this study, we developed a custom panel of 61 mtDNA markers for high-throughput classification of European, African, and Native American/Asian mitochondrial haplogroup lineages. Using these mtDNA markers, we constructed a mitochondrial haplogroup classification tree and classified 18,832 participants from the National Health and Nutrition Examination Surveys (NHANES). To our knowledge, this is the largest study to date characterizing mitochondrial haplogroups in a population-based sample from the United States, and the first study characterizing mitochondrial haplogroup distributions in self-identified Mexican Americans separately from Hispanic Americans of other descent. We observed clear differences in the distribution of maternal genetic ancestry consistent with proposed admixture models for these subpopulations, underscoring the genetic heterogeneity of the United States Hispanic population. The mitochondrial haplogroup distributions in the other self-identified racial/ethnic groups within NHANES were largely comparable to previous studies. Mitochondrial haplogroup classification was highly concordant with self-identified race/ethnicity (SIRE) in non-Hispanic whites (94.8 %), but was considerably lower in admixed populations including non-Hispanic blacks (88.3 %), Mexican Americans (81.8 %), and other Hispanics (61.6 %), suggesting SIRE does not accurately reflect maternal genetic ancestry, particularly in populations with greater proportions of admixture. Thus, it is important to consider inconsistencies between SIRE and genetic ancestry when performing genetic association studies. The mitochondrial haplogroup data that we have generated, coupled with the epidemiologic variables in NHANES, is a valuable resource for future studies investigating the contribution of mtDNA variation to human health and disease.     

High-resolution Y chromosome haplotypes of Israeli and Palestinian Arabs reveal geographic substructure and substantial overlap with haplotypes of Jews

High-resolution Y chromosome haplotype analysis was performed in 143 paternally unrelated Israeli and Palestinian Moslem Arabs (I&P Arabs) by screening for 11 binary polymorphisms and six microsatellite loci. Two frequent haplotypes were found among the 83 detected: the modal haplotype of the I&P Arabs (approximately 14%) was spread throughout the region, while its one-step microsatellite neighbor, the modal haplotype of the Galilee sample (approximately 8%), was mainly restricted to the north. Geographic substructuring within the Arabs was observed in the highlands of Samaria and Judea. Y chromosome variation in the I&P Arabs was compared to that of Ashkenazi and Sephardic Jews, and to that of North Welsh individuals. At the haplogroup level, defined by the binary polymorphisms only, the Y chromosome distribution in Arabs and Jews was similar but not identical. At the haplotype level, determined by both binary and microsatellite markers, a more detailed pattern was observed. Single-step microsatellite networks of Arab and Jewish haplotypes revealed a common pool for a large portion of Y chromosomes, suggesting a relatively recent common ancestry. The two modal haplotypes in the I&P Arabs were closely related to the most frequent haplotype of Jews (the Cohen modal haplotype). However, the I&P Arab clade that includes the two Arab modal haplotypes (and makes up 32% of Arab chromosomes) is found at only very low frequency among Jews, reflecting divergence and/or admixture from other populations.     

Abraham's Children in the Genome Era: Major Jewish Diaspora Populations Comprise Distinct Genetic Clusters with Shared Middle Eastern Ancestry

For more than a century, Jews and non-Jews alike have tried to define the relatedness of contemporary Jewish people. Previous genetic studies of blood group and serum markers suggested that Jewish groups had Middle Eastern origin with greater genetic similarity between paired Jewish populations. However, these and successor studies of monoallelic Y chromosomal and mitochondrial genetic markers did not resolve the issues of within and between-group Jewish genetic identity. Here, genome-wide analysis of seven Jewish groups (Iranian, Iraqi, Syrian, Italian, Turkish, Greek, and Ashkenazi) and comparison with non-Jewish groups demonstrated distinctive Jewish population clusters, each with shared Middle Eastern ancestry, proximity to contemporary Middle Eastern populations, and variable degrees of European and North African admixture. Two major groups were identified by principal component, phylogenetic, and identity by descent (IBD) analysis: Middle Eastern Jews and European/Syrian Jews. The IBD segment sharing and the proximity of European Jews to each other and to southern European populations suggested similar origins for European Jewry and refuted large-scale genetic contributions of Central and Eastern European and Slavic populations to the formation of Ashkenazi Jewry. Rapid decay of IBD in Ashkenazi Jewish genomes was consistent with a severe bottleneck followed by large expansion, such as occurred with the so-called demographic miracle of population expansion from 50,000 people at the beginning of the 15^th century to 5,000,000 people at the beginning of the 19^th century. Thus, this study demonstrates that European/Syrian and Middle Eastern Jews represent a series of geographical isolates or clusters woven together by shared IBD genetic threads.     

Two sources of the Russian patrilineal heritage in their Eurasian context

Progress in the mapping of population genetic substructure provides a core source of data for the reconstruction of the demographic history of our species and for the discovery of common signals relevant to disease research: These two aspects of enquiry overlap in their empirical data content and are especially informative at continental and subcontinental levels. In the present study of the variation of the Y chromosome pool of ethnic Russians, we show that the patrilineages within the pre-Ivan the Terrible historic borders of Russia have two main distinct sources. One of these antedates the linguistic split between West and East Slavonic-speaking people and is common for the two groups; the other is genetically highlighted by the pre-eminence of haplogroup (hg) N3 and is most parsimoniously explained by extensive assimilation of (or language change in) northeastern indigenous Finno-Ugric tribes. Although hg N3 is common for both East European and Siberian Y chromosomes, other typically Siberian or Mongolian hgs (Q and C) have negligible influence within the studied Russian Y chromosome pool. The distribution of all frequent Y chromosome haplogroups (which account for 95% of the Y chromosomal spectrum in Russians) follows a similar north-south clinal pattern among autosomal markers, apparent from synthetic maps. Multidimensional scaling (MDS) plots comparing intra ethnic and interethnic variation of Y chromosome in Europe show that although well detectable, intraethnic variation signals do not cross interethnic borders, except between Poles, Ukrainians, and central-southern Russians, thereby revealing their overwhelmingly shared patrilineal ancestry.     

A new topology of the human Y chromosome haplogroup E1b1 (E-P2) revealed through the use of newly characterized binary polymorphisms

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.     

The relationship between Y chromosome DNA haplotypes and Y chromosome deletions leading to male infertility

Microdeletions on the short arm of the Y chromosome have defined three non-overlapping regions (AZFa, b, c) recurrently deleted among infertile males. These regions contain several genes or gene families involved in male germ-cell development and maintenance. Even though a meiotic origin for these microdeletions is assumed, the mechanisms and causes leading to microdeletion formation are largely unknown. In order to assess whether some Y chromosome groups (or haplogroups) are predisposed to, or protected against, deletion formation during male meiosis, we have defined and compared Y chromosome haplogroup distribution in a group of infertile/subfertile males harbouring Yq deletions and in a relevant Northwestern European control population. Our analyses suggest that Y chromosome deletion formation is, at least in the study populations, a stochastic event independent of the Y chromosome background on which they arise and may be caused by other genetic and/or environmental factors.     

Y chromosome--a tool in infertility studies of Latvian population

Human Y chromosome is used as a tool in male infertility and population genetic studies. The aims of this research were to analyse the prevalence of Y chromosome microdeletions among infertile Latvian men, and to identify possible lineages of Y chromosome that may be at increased risk of developing infertility. A study encompassed 105 infertile men with different spermatogenic disturbances. Deletions on Y chromosome were detected in 5 out of 105 (approximately 5%) cases analysed in this study. Three of them carried deletion in AZFc region and two individuals had AZFa + b + c deletion. Study of Y chromosome haplogroups showed that N3a1 and R1a1 lineages were found less frequently in the infertile male group compared to ethnic Latvian group, however K* cluster was predominantly found in infertile male Y chromosomes. Conclusions: 1) Our study advocates running Y chromosome microdeletion analyses only in cases of severe form of infertility; 2) Y chromosome haplogroup analysis showed statistically significant tendencies that some haplogroups are more common in ethnic male group, but others are more common in infertile males.     

Excavating Y-chromosome haplotype strata in Anatolia

Analysis of 89 biallelic polymorphisms in 523 Turkish Y chromosomes revealed 52 distinct haplotypes with considerable haplogroup substructure, as exemplified by their respective levels of accumulated diversity at ten short tandem repeat (STR) loci. The major components (haplogroups E3b, G, J, I, L, N, K2, and R1; 94.1%) are shared with European and neighboring Near Eastern populations and contrast with only a minor share of haplogroups related to Central Asian (C, Q and O; 3.4%), Indian (H, R2; 1.5%) and African (A, E3*, E3a; 1%) affinity. The expansion times for 20 haplogroup assemblages was estimated from associated STR diversity. This comprehensive characterization of Y-chromosome heritage addresses many multifaceted aspects of Anatolian prehistory, including: (1) the most frequent haplogroup, J, splits into two sub-clades, one of which (J2) shows decreasing variances with increasing latitude, compatible with a northward expansion; (2) haplogroups G1 and L show affinities with south Caucasus populations in their geographic distribution as well as STR motifs; (3) frequency of haplogroup I, which originated in Europe, declines with increasing longitude, indicating gene flow arriving from Europe; (4) conversely, haplogroup G2 radiates towards Europe; (5) haplogroup E3b3 displays a latitudinal correlation with decreasing frequency northward; (6) haplogroup R1b3 emanates from Turkey towards Southeast Europe and Caucasia and; (7) high resolution SNP analysis provides evidence of a detectable yet weak signal (<9%) of recent paternal gene flow from Central Asia. The variety of Turkish haplotypes is witness to Turkey being both an important source and recipient of gene flow.