Origin of Amerindian Y-chromosomes as inferred by the analysis of six polymorphic markers

Uniparental origins of the admixed Argentine Patagonia

OBJECTIVES

We aimed to contribute to the understanding of the ancient geographic origins of the uniparentally inherited markers in modern admixed Argentinian populations from central Patagonia with new information provided for the city of Trelew. We attempted to highlight the importance of combining different genetic markers when studying population history.

METHODS

The mtDNA control region sequence was typified in 89 individuals and 12 Y-STR and 15 Y-SNP loci were analyzed in 66 males. With these data, analysis of molecular variance and Network analyses were carried out. We exhaustively compared the modern data with ancient mtDNA information. Finally, we tested the differences in continental origins estimated by uniparental and previously published biparental markers.

RESULTS

Native American mtDNAs (53.9%) increased when maternal ancestors were born in the northern (81.8%) and southern (58.5%) regions of Argentina or in Chile (77.8%). Population substructure was only observed for Y-chromosome haplotypes. Some mtDNA haplogroups have been present in the area for at least ca. 2762-2430 and ca. 500 (D1g and D1g4 haplogroups) and ca. 6736 and ca. 6620 (C1b and C1c haplogroups) years, respectively. In contrast, haplogroups B2i2 and C1b13, frequent in modern Patagonia populations, had not been found in previous ancient DNA studies of the region.

CONCLUSIONS

The results suggest that Native American ancestry is well preserved in the region. Trelew samples had characteristic native mtDNA haplogroups previously described in Chilean and Argentine Patagonian populations, but not observed in ancient samples until now. These findings support the idea that these lineages have a recent regional origin. Finally, the estimated proportions of continental ancestry depend on the genetic marker analyzed.

Identifying new lineages in the Y chromosome of Colombian Amazon indigenous populations

OBJECTIVES

The Y chromosome has highly informative markers, such as single-nucleotide polymorphisms (SNPs), that are useful for making historical inferences about the settlement of the Americas. However, the scarcity of these markers has limited their use. This study aims to identify new SNPs and increase the phylogenetic resolution of haplogroup Q for the Americas, mainly focusing on the lineages of the Amazon region.

MATERIALS AND METHODS

Next-generation sequencing was performed on two Y chromosomes belonging to haplogroup Q-M3 using samples with divergent short tandem repeat haplotypes from the Colombian Amazon, and 14 of the new variants identified were selected for characterization in 207 samples of indigenous Colombians belonging to haplogroup Q-M3.

RESULTS

This methodology allowed us to establish nine new lineages within Q-M3, including its paragroups. The most basal lineages were predominant in communities of Andean origin, such as the Embera-Katio, the Nasas, and the Pastos. In contrast, the most distal lineages were restricted to inhabitants of the Amazon region of Vaupés.

DISCUSSION

The SNPs reported here advance the development of subhaplogroups of Q-M3 with a higher level of phylogenetic resolution than has been previously reported, which allowed the differentiation between populations that inhabit two regions of Vaupes area: the Pirá-Paraná region and the upper and middle sections of the Vaupés River, and the region encompassing the Papurí River and the lower Vaupés. They are very useful for the microevolutionary analysis of the Amerindian populations of Colombia and of the Americas.

Greater prevalence of Y chromosome Q1a3a haplogroup in Y-microdeleted Chilean men: a case-control study

PURPOSE

To determine the prevalence of South Amerindian Y chromosome in Chilean patients with spermatogenic failure and their association with classical and/or AZFc-partial Y chromosome deletions.

METHODS

We studied 400 men, 218 with secretory azo/oligozoospermia (cases) and 182 controls (116 fertile and/or normozoospermic, and 66 azoospermic with normal spermatogenesis). After a complete testicular characterization (physical evaluation, hormonal and/or biopsy) peripheral blood was drawn to obtain DNA for Y chromosome microdeletions, AZFc-partial deletions and biallelic analysis by allele specific polymerase chain reaction (PCR) of the M3 (rs3894) single nucleotide polymorphism (SNP).

RESULTS

Classical AZF microdeletions were found in 23 cases (Y-microdeleted). AZFc-partial deletions were observed in 10 cases (6 "gr/gr", 3 "b2/b3" and 1 "b1/b3") and 4 controls (4 "gr/gr"). The AZFc-partial deletions were mainly associated with the absence of DAZ1/DAZ2 (64 %). No significant differences in the prevalence of AZFc-partial deletions were observed between cases and controls. We observed a significant higher proportion of the Q1a3a haplogroup in Y-microdeleted men compared to patients with spermatogenic failure without deletions and control men (P<0.01 and P<0.05, respectively by Bonferroni test). Among them, patients with AZFb deletions had an increased prevalence of the Q1a3a haplogroup compared to controls, cases without deletions and to those with complete or partial-AZFc deletions (P<0.01, Bonferroni test).

CONCLUSIONS

The Q1a3a South Amerindian lineage seems to increase the susceptibility to non AZFc microdeletions. On the other hand, in Chilean population the AZFc-partial deletions ("gr/gr", "b1/b3" and/or "b2/b3") does not seem to predispose to severe spermatogenic impairment.

Y-chromosome analysis reveals genetic divergence and new founding native lineages in Athapaskan- and Eskimoan-speaking populations

For decades, the peopling of the Americas has been explored through the analysis of uniparentally inherited genetic systems in Native American populations and the comparison of these genetic data with current linguistic groupings. In northern North America, two language families predominate: Eskimo-Aleut and Na-Dene. Although the genetic evidence from nuclear and mtDNA loci suggest that speakers of these language families share a distinct biological origin, this model has not been examined using data from paternally inherited Y chromosomes. To test this hypothesis and elucidate the migration histories of Eskimoan- and Athapaskan-speaking populations, we analyzed Y-chromosomal data from Inuvialuit, Gwich'in, and Tłįch populations living in the Northwest Territories of Canada. Over 100 biallelic markers and 19 chromosome short tandem repeats (STRs) were genotyped to produce a high-resolution dataset of Y chromosomes from these groups. Among these markers is an SNP discovered in the Inuvialuit that differentiates them from other Aboriginal and Native American populations. The data suggest that Canadian Eskimoan- and Athapaskan-speaking populations are genetically distinct from one another and that the formation of these groups was the result of two population expansions that occurred after the initial movement of people into the Americas. In addition, the population history of Athapaskan speakers is complex, with the Tłįch being distinct from other Athapaskan groups. The high-resolution biallelic data also make clear that Y-chromosomal diversity among the first Native Americans was greater than previously recognized.

YAP insertion signature in South Asia

A total of 2169 samples from 21 tribal populations from different regions of India were scanned for the Y-chromosome Alu polymorphism. This study reports, for the first time, high frequencies (8-65%) of Y Alu polymorphic (YAP) insertion in northeast Indian tribes. All seven Jarawa samples from the Andaman and Nicobar islands had the YAP insertion, in conformity with an earlier study of Andaman Islanders. One isolated case with haplotype E* was found in Dungri Bhill, a western Indian population, while YAP insertion in northeast India and Andaman tribes was found in association with haplotype D* (M168, M174). YAP insertion frequencies reported in the mainland Indian populations are negligible, according to previous studies. Genetic drift may be the causative factor for the variable frequency of the YAP insertion in the mainland populations, while the founder effect may have resulted in the highest incidence of haplotype D among the Andaman Islanders. The results of YAP insertion and the evidence of previous mtDNA studies indicate an early out of Africa migration to the Andaman and Nicobar Islands. The findings of YAP insertion in northeast Indian tribes are very significant for understanding the evolutionary history of the region.

Genetic admixture, relatedness, and structure patterns among Mexican populations revealed by the Y-chromosome

Y-linked markers are suitable loci to analyze genetic diversity of human populations, offering knowledge of medical, forensic, and anthropological interest. In a population sample of 206 Mestizo males from western Mexico, we analyzed two binary loci (M3 and YAP) and six Y-STRs, adding to the analysis data of Mexican Mestizos and Amerindians, and relevant worldwide populations. The paternal ancestry estimated in western Mexican-Mestizos was mainly European (60-64%), followed by Amerindian (25-21%), and African ( approximately 15%). Significant genetic heterogeneity was established between Mestizos from western (Jalisco State) and northern Mexico (Chihuahua State) compared with Mexicans from the center of the Mexican Republic (Mexico City), this attributable to higher European ancestry in western and northern than in central and southeast populations, where higher Amerindian ancestry was inferred. This genetic structure has important implications for medical and forensic purposes. Two different Pre-Hispanic evolutionary processes were evident. In Mesoamerican region, populations presented higher migration rate (N(m) = 24.76), promoting genetic homogeneity. Conversely, isolated groups from the mountains and canyons of the Western and Northern Sierra Madre (Huichols and Tarahumaras, respectively) presented a lower migration rate (N(m) = 10.27) and stronger genetic differentiation processes (founder effect and/or genetic drift), constituting a Pre-Hispanic population substructure. Additionally, Tarahumaras presented a higher frequency of Y-chromosomes without Q3 that was explained by paternal European admixture (15%) and, more interestingly, by a distinctive Native-American ancestry. In Purepechas, a special admixture process involving preferential integration of non-Purepecha women in their communities could explain contrary genetic evidences (autosomal vs. Y-chromosome) for this tribe.

Polymorphic Alu insertions among Mayan populations

The Mayan homeland within Mesoamerica spans five countries: Belize, El Salvador, Guatemala, Honduras and Mexico. There are indications that the people we call the Maya migrated from the north to the highlands of Guatemala as early as 4000 B.C. Their existence was village-based and agricultural. The culture of these Preclassic Mayans owes much to the earlier Olmec civilization, which flourished in the southern portion of North America. In this study, four different Mayan groups were examined to assess their genetic variability. Ten polymorphic Alu insertion (PAI) loci were employed to ascertain the genetic affinities among these Mayan groups. North American, African, European and Asian populations were also examined as reference populations. Our results suggest that the Mayan groups examined in this study are not genetically homogeneous.

South Amerindian craniofacial morphology: Diversity and implications for Amerindian evolution

The most compelling models concerning the peopling of the Americas consider that modern Amerindians share a common biological pattern, showing affinities with populations of the Asian Northeast. The aim of the present study was to assess the degree of variation of craniofacial morphology of South American Amerindians in a worldwide context. Forty-three linear variables were analyzed on crania derived from American, Asian, Australo-Melanesian, European, South-Saharan African, and Polynesian regions. South America was represented by seven Amerindian samples. In order to understand morphologic diversity among Amerindians of South America, variation was estimated using regions and local populations as units of analysis. Variances and F(ST) values were calculated for each unit, respectively. Both analyses indicated that morphologic variation in Southern Amerindians is extremely high: an F(ST) of 0.01531 was obtained for Southern Amerindians, and values from 0.0371-0.1205 for other world regions. Some aspects linked to the time and mode of the peopling of the Americas and various microevolutionary processes undergone by Amerindians are discussed. Some of the alternatives proposed to explain this high variation include: a greater antiquity of the peopling than what is mostly accepted, a peopling by several highly differentiated waves, an important effect of genetic drift, and gene flow with Paleoamericans. A combination of some of these alternatives explains at least some of the variation.

Phylogenetic information in polymorphic L1 andAlu insertions from East Asians and Native American populations

This study attempts to ascertain genetic affinities between Native American and East Asian populations by analyzing four polymorphic Alu insertions (PAIs) and three L1 polymorphic loci. These two genetic systems demonstrated strong congruence when levels of diversity and genetic distances were considered. Overall, genetic relatedness within Native American groups does not correlate with geographical and linguistic structure, although strong grouping for Native Americans with East Asians was demonstrated, with clear discrimination from African and European groups. Most of the variation was assigned to differences occurring within groups, but the interpopulation variation found for South Amerindians was recognizably higher in comparison to the other sampled groups of populations. Our data suggest that bottleneck events followed by strong influence of genetic drift in the process of the peopling of the Americas may have been determinant factors in delineating the genetic background of present-day South Amerindians. Since no clear subgroups were detected within Native Americans and East Asians, there is no indication of multiple waves in the early colonization of the New World.

DYS19 and DYS199 loci in a Chilean population of mixed ancestry

The current Chilean population originated from admixture between aboriginal populations (Amerindians) and Spanish conquerors of European origin. Consequently, the unions that gave rise to the Chilean population were chiefly between Spanish males and aboriginal females, and not the converse. To test the hypothesis that the Y chromosome of the Chilean population is mainly of Spanish origin, while the other chromosomes are from mixed (European and aboriginal) origin, we studied the DYS19 and DYS199 loci in two samples. One sample was obtained from a high socioeconomic stratum, while a second sample was from a low stratum. We studied male blood donors (N = 187) from Santiago, the capital of the country. Subjects were typed for the autosomal ABO and Rh (locus D) blood groups, and for the Y-linked DYS19 and the DYS199 loci, reported as Y-chromosome haplotypes. The aboriginal admixture was estimated for each genetic marker. The percentage of aboriginal admixture was 38.17% for the ABO system and 31.28% for the Rh system in the low socioeconomic stratum and 19.22% and 22.5%, respectively, in the high stratum. Y-chromosome haplotype frequencies constructed from the DYS19 and DYS199 loci demonstrated that the main haplotypes were DYS19*14/DYS199 C, as is often the case with many European populations, and DYS19*13/DYS199 C. The aboriginal admixture from Y-haplotype frequencies was estimated to be 15.83% in the low socioeconomic stratum and 6.91% in the high stratum. These values are lower than the values found using autosomal genetic markers, and are consistent with the historical background of the population studied. This study highlights the population genetic consequences of the asymmetric pattern of genome admixture between two ancestral populations (European and Amerindian).

Mitochondrial DNA and Y chromosome diversity and the peopling of the Americas: Evolutionary and demographic evidence

A number of important insights into the peopling of the New World have been gained through molecular genetic studies of Siberian and Native American populations. While there is no complete agreement on the interpretation of the mitochondrial DNA (mtDNA) and Y chromosome (NRY) data from these groups, several generalizations can be made. To begin with, the primary migration of ancestral Asians expanded from south-central Siberia into the New World and gave rise to ancestral Amerindians. The initial migration seems to have occurred between 20,000-15,000 calendar years before present (cal BP), i.e., before the emergence of Clovis lithic sites (13,350-12,895 cal BP) in North America. Because an interior route through northern North America was unavailable for human passage until 12,550 cal BP, after the last glacial maximum (LGM), these ancestral groups must have used a coastal route to reach South America by 14,675 cal BP, the date of the Monte Verde site in southern Chile. The initial migration appears to have brought mtDNA haplogroups A-D and NRY haplogroups P-M45a and Q-242/Q-M3 to the New World, with these genetic lineages becoming widespread in the Americas. A second expansion that perhaps coincided with the opening of the ice-free corridor probably brought mtDNA haplogroup X and NRY haplogroups P-M45b, C-M130, and R1a1-M17 to North and Central America. Finally, populations that formerly inhabited Beringia expanded into northern North America after the LGM, and gave rise to Eskimo-Aleuts and Na-Dené Indians.

Y-chromosome evidence for differing ancient demographic histories in the Americas

To scrutinize the male ancestry of extant Native American populations, we examined eight biallelic and six microsatellite polymorphisms from the nonrecombining portion of the Y chromosome, in 438 individuals from 24 Native American populations (1 Na Dené and 23 South Amerinds) and in 404 Mongolians. One of the biallelic markers typed is a recently identified mutation (M242) characterizing a novel founder Native American haplogroup. The distribution, relatedness, and diversity of Y lineages in Native Americans indicate a differentiated male ancestry for populations from North and South America, strongly supporting a diverse demographic history for populations from these areas. These data are consistent with the occurrence of two major male migrations from southern/central Siberia to the Americas (with the second migration being restricted to North America) and a shared ancestry in central Asia for some of the initial migrants to Europe and the Americas. The microsatellite diversity and distribution of a Y lineage specific to South America (Q-M19) indicates that certain Amerind populations have been isolated since the initial colonization of the region, suggesting an early onset for tribalization of Native Americans. Age estimates based on Y-chromosome microsatellite diversity place the initial settlement of the American continent at approximately 14,000 years ago, in relative agreement with the age of well-established archaeological evidence.

Early population differentiation in extinct aborigines from Tierra del Fuego-Patagonia: Ancient mtDNA sequences and Y-Chromosome STR characterization

Ancient mtDNA was successfully recovered from 24 skeletal samples of a total of 60 ancient individuals from Patagonia-Tierra del Fuego, dated to 100-400 years BP, for which consistent amplifications and two-strand sequences were obtained. Y-chromosome STRs (DYS434, DYS437, DYS439, DYS393, DYS391, DYS390, DYS19, DYS389I, DYS389II, and DYS388) and the biallelic system DYS199 were also amplified, Y-STR alleles could be characterized in nine cases, with an average of 4.1 loci per sample correctly typed. In two samples of the same ethnic group (Aonikenk), an identical and complete eight-loci haplotype was recovered. The DYS199 biallelic system was used as a control of contamination by modern DNA and, along with DYS19, as a marker of American origin. The analysis of both mtDNA and Y-STRs revealed DNA from Amerindian ancestry. The observed polymorphisms are consistent with the hypothesis that the ancient Fuegians are close to populations from south-central Chile and Argentina, but their high nucleotide diversity and the frequency of single lineages strongly support early genetic differentiation of the Fuegians through combined processes of population bottleneck, isolation, and/or migration, followed by strong genetic drift. This suggests an early genetic diversification of the Fuegians right after their arrival at the southernmost extreme of South America.

Dissimilarities in the process of formation of Curiaú, a semi-isolated Afro-Brazilian population of the Amazon region

The genetic consequences of the social policy of the past in relation to the formation of Afro-Brazilian societies are interesting and have been studied at various biological levels (classical polymorphisms and the mitochondrial and nuclear levels. These allow the estimation of the contribution of African genes and the participation of other ethnic groups in the formation of these communities. With this objective, uniparental systems of exclusively maternal (mtDNA) or paternal (Y-DNA) inheritance in the Curiaú community were analyzed. The results demonstrate a differential contribution of the maternal and paternal genetic systems. Thirty-three sequences were identified by mtDNA analysis; 53% showing an African and 47% an Amerindian origin. For the paternal system, 57% were of African, 37% of European, and 6% of Amerindian origin.

Mosaic AZF deletions and susceptibility to testicular tumors

We tested for azoospermia factor (AZF) deletions 17 loci corresponding to AZF subintervals a-d in 17 cases of testicular tumors occurring in Finns. While DNA samples from 48 CEPH and 32 Finnish males showed no deletions, patients with testicular cancer displayed AZF deletion mosaicisms in various non-tumor tissues (13 cases) and specific deletion haplotypes in tumor tissues (10 cases). Two of the cases with AZF deletions were testicular non-Hodgkin lymphomas indicating that Y-microdeletions appear also in malignancies other than seminoma and non-seminoma tumors. In good agreement with this assumption, we detected one AZF deletion in normal cells from 1 of 5 HNPCC cases, heterozygous for an MLH1 mutation. We propose that AZF deletions occur in early embryogenesis due to mutations of TSPY, mismatch repair (MMR), or X-specific genes. Since fathers of testicular, tumor cases did not exhibit AZF deletions, we assumed they were not carriers of the mutation inducing AZF deletion-mosaicisms. Therefore, tumor cases should have received the MMR gene or X mutations via the maternal lineage, or for the case of TSPY and MMR genes via a sperm carrying a mutation occurred in the paternal germ-cell line. We consider AZF microdeletions in non-tumor cells to be part of a broader pattern of chromosome instability producing susceptibility to testicular tumors. Clonal transformation and expansion of one of these tumor-susceptible cell lineages give rise to testicular tumors showing genome anomalies characteristic of testicular cancers (i12p, LOH and genetic imbalance for various autosomal regions, Y- and autosomal MSI, specific AZF deletion haplotypes).

Molecular variability in Amerindians: widespread but uneven information

A review was made in relation to the molecular variability present in North, Central, and South American Indian populations. It involved results from ancient DNA, mitochondrial DNA in extant populations, HLA and other autosomal markers, X and Y chromosome variation, as well as data from parasitic viruses which could show coevolutionary changes. The questions considered were their origin, ways in which the early colonization of the continent took place, types and levels of the variability which developed, peculiarities of the Amerindian evolutionary processes, and eventual genetic heterogeneity which evolved in different geographical areas. Although much information is already available, it is highly heterogeneous in relation to populations and types of genetic systems investigated. Unfortunately, the present trend of favoring essentially applied research suggest that the situation will not basically improve in the future.

Y-chromosomal DNA variation in Pakistan

Eighteen binary polymorphisms and 16 multiallelic, short-tandem-repeat (STR) loci from the nonrecombining portion of the human Y chromosome were typed in 718 male subjects belonging to 12 ethnic groups of Pakistan. These identified 11 stable haplogroups and 503 combination binary marker/STR haplotypes. Haplogroup frequencies were generally similar to those in neighboring geographical areas, and the Pakistani populations speaking a language isolate (the Burushos), a Dravidian language (the Brahui), or a Sino-Tibetan language (the Balti) resembled the Indo-European-speaking majority. Nevertheless, median-joining networks of haplotypes revealed considerable substructuring of Y variation within Pakistan, with many populations showing distinct clusters of haplotypes. These patterns can be accounted for by a common pool of Y lineages, with substantial isolation between populations and drift in the smaller ones. Few comparative genetic or historical data are available for most populations, but the results can be compared with oral traditions about origins. The Y data support the well-established origin of the Parsis in Iran, the suggested descent of the Hazaras from Genghis Khan's army, and the origin of the Negroid Makrani in Africa, but do not support traditions of Tibetan, Syrian, Greek, or Jewish origins for other populations.

The dual origin and Siberian affinities of Native American Y chromosomes

The Y chromosomes of 549 individuals from Siberia and the Americas were analyzed for 12 biallelic markers, which defined 15 haplogroups. The addition of four microsatellite markers increased the number of haplotypes to 111. The major Native American founding lineage, haplogroup M3, accounted for 66% of male Y chromosomes and was defined by the biallelic markers M89, M9, M45, and M3. The founder haplotype also harbored the microsatellite alleles DYS19 (10 repeats), DYS388 (11 repeats), DYS390 (11 repeats), and DYS391 (10 repeats). In Siberia, the M3 haplogroup was confined to the Chukotka peninsula, adjacent to Alaska. The second major group of Native American Y chromosomes, haplogroup M45, accounted for about one-quarter of male lineages. M45 was subdivided by the biallelic marker M173 and by the four microsatellite loci alleles into two major subdivisions: M45a, which is found throughout the Americas, and M45b, which incorporates the M173 variant and is concentrated in North and Central America. In Siberia, M45a haplotypes, including the direct ancestor of haplogroup M3, are concentrated in Middle Siberia, whereas M45b haplotypes are found in the Lower Amur River and Sea of Okhotsk regions of eastern Siberia. Among the remaining 5% of Native American Y chromosomes is haplogroup RPS4Y-T, found in North America. In Siberia, this haplogroup, along with haplogroup M45b, is concentrated in the Lower Amur River/Sea of Okhotsk region. These data suggest that Native American male lineages were derived from two major Siberian migrations. The first migration originated in southern Middle Siberia with the founding haplotype M45a (10-11-11-10). In Beringia, this gave rise to the predominant Native American lineage, M3 (10-11-11-10), which crossed into the New World. A later migration came from the Lower Amur/Sea of Okhkotsk region, bringing haplogroup RPS4Y-T and subhaplogroup M45b, with its associated M173 variant. This migration event contributed to the modern genetic pool of the Na-Dene and Amerinds of North and Central America.

Craniometric variation and the settlement of the Americas: testing hypotheses by means of R-matrix and matrix correlation analyses

New archaeological findings and the incorporation of new South American skull samples have raised fundamental questions for the classical theories of the Americas' settlement. The aim of this study was to estimate craniometric variability among several Asian and Native American populations in order to test goodness of fit of the data to different models of ancient population entries and dispersions into the New World. Our data set includes Howells' variables recorded on East Asian, North American, and South American natives (except for Na-Dene speakers). Five Fuego-Patagonian samples and one Paleoamerican sample were also included. A multivariate extension of the R-matrix method for quantitative traits was used to obtain Fst values, which were considered estimations of intergroup variation. Three main models for the peopling of the New World were represented in hypothetical design matrices. Matrix permutation tests were performed to quantify the fit of the observed data with 1) geographical separation of the samples and 2) three ways of settlement, which were the Three Migration Model (TMM), the Single Wave Migration model (SWM), and the Two Components Settlement Model (TCS). R-matrix results showed high levels of heterogeneity among Native Americans. Matrix permutation analyses suggested that the model involving high Amerindian heterogeneity and two different morphological patterns or components (derived "Mongoloid" vs. generalized "non-Mongoloid") explains better the variation observed, even when the effects of geographical separation are removed. Whether these patterns arose as a result of two separate migration events or by local evolution from Paleoamericans to Amerindians remains unresolved.

Heterogeneity of Y chromosome markers among Brazilian Amerindians

The allele frequency distribution of DYS19 and DYS199 loci were analyzed in 59 Brazilian Amerindians from five tribes from the Amazon region (Zoé, Awá-Guajá, Urubú-Kaapór, Katuena, and Kayapó, Xikrin of Bacajá village). Three different alleles of the DYS19 microsatellite (182-bp, 186-bp, and 190-bp) were found at average frequencies of 0.08, 0.85, and 0.07, respectively. The DYS199-T allele was identified in 78% of the Amerindians studied (43/55), the frequencies varying from 0.46-0.93. Four different haplotypes were found, the combination DYS19-186/DYS199-T being the most common (average frequency of 0.65), followed by DYS19-186/DYS199-C with an average frequency of 0.22. These four haplotypes have been found in five other Brazilian tribes, and most of them were also identified in Native populations from South, Central and North America. The observed variability at the DYS19 microsatellite is probably due to forward or back mutations from the putative ancestral 186-bp allele, since the mutation rate of this locus is high and the post-Columbian admixture of the Brazilian tribes studied is very low or undetectable to explain these data. On the other hand, the DYS19/DYS199 haplotype distribution may suggest that the two most common haplotypes (186-bp/T and 186-bp/C) were present among the population(s) that peopled the New World. Am. J. Hum. Biol. 11:481-487, 1999. Copyright 1999 Wiley-Liss, Inc.

Correlation between molecular and conventional genealogies in Aicuña: a rural population from Northwestern Argentina

Aicuña is a village in the northwest of Argentina, located about 300 km south of La Rioja city, in the province of La Rioja. The population of Aicuña derives from a founder couple established in the uninhabited Aicuña valley in the early years of the 17th century. Due to land ownership litigation, the descendants maintained a well-documented genealogy that extends for 12 generations, comprising more than 8,000 individuals. From the historical pedigree of Aicuña, we selected 14 males with direct patrilineal descent from the 2 most ancient male founders, and 23 donors (9 females and 14 males) with direct matrilineal descent from the most ancient female founder. All 3 founders lived in the 17th century. We collected DNA from buccal swabs and characterized the mitochondrial DNA (mtDNA) and Y haplotypes using 14 Y-specific markers, 11 mtDNA polymorphic markers and sequencing of the mt hypervariable regions 1 and 2. We found four different Y haplotypes: Y1 and Y2 haplotypes of European origin corresponding to the founder ancestors Francisco Páez de Espinoza and Apolinario Ormeño, which were shared by 6 and 3 donors, respectively. Three males selected as Ormeño patrilineal descendants showed a different Y haplotype (Y3), probably originated by erroneous paternity registration due to illegitimacy. The remaining case (haplotype Y4), also assumed to belong to the Ormeño lineage, was probably also due to an erroneously registered paternity. Twenty-two donors showed an association of mtDNA markers corresponding to the Amerindian haplotype A2. The founder of this matrilineage could be traced back for more than 14 generations. The haplotype B of one remaining female did not correspond with the historical pedigree and could be due to an error in the genealogy registration. Our results showed an 85% agreement between conventional and molecular genealogies, with mtDNA markers being Amerindian, and Y markers being European. The methodology used in this report is a tool which could potentially be employed as a precedent for land ownership by Aicuña villagers and Amerindian populations.

Y-chromosomal SNPs in Finno-Ugric-speaking populations analyzed by minisequencing on microarrays

An increasing number of single nucleotide polymorphisms (SNPs) on the Y chromosome are being identified. To utilize the full potential of the SNP markers in population genetic studies, new genotyping methods with high throughput are required. We describe a microarray system based on the minisequencing single nucleotide primer extension principle for multiplex genotyping of Y-chromosomal SNP markers. The system was applied for screening a panel of 25 Y-chromosomal SNPs in a unique collection of samples representing five Finno--Ugric populations. The specific minisequencing reaction provides 5-fold to infinite discrimination between the Y-chromosomal genotypes, and the microarray format of the system allows parallel and simultaneous analysis of large numbers of SNPs and samples. In addition to the SNP markers, five Y-chromosomal microsatellite loci were typed. Altogether 10,000 genotypes were generated to assess the genetic diversity in these population samples. Six of the 25 SNP markers (M9, Tat, SRY10831, M17, M12, 92R7) were polymorphic in the analyzed populations, yielding six distinct SNP haplotypes. The microsatellite data were used to study the genetic structure of two major SNP haplotypes in the Finns and the Saami in more detail. We found that the most common haplotypes are shared between the Finns and the Saami, and that the SNP haplotypes show regional differences within the Finns and the Saami, which supports the hypothesis of two separate settlement waves to Finland.

Y-chromosomal diversity in Europe is clinal and influenced primarily by geography, rather than by language

Clinal patterns of autosomal genetic diversity within Europe have been interpreted in previous studies in terms of a Neolithic demic diffusion model for the spread of agriculture; in contrast, studies using mtDNA have traced many founding lineages to the Paleolithic and have not shown strongly clinal variation. We have used 11 human Y-chromosomal biallelic polymorphisms, defining 10 haplogroups, to analyze a sample of 3,616 Y chromosomes belonging to 47 European and circum-European populations. Patterns of geographic differentiation are highly nonrandom, and, when they are assessed using spatial autocorrelation analysis, they show significant clines for five of six haplogroups analyzed. Clines for two haplogroups, representing 45% of the chromosomes, are continentwide and consistent with the demic diffusion hypothesis. Clines for three other haplogroups each have different foci and are more regionally restricted and are likely to reflect distinct population movements, including one from north of the Black Sea. Principal-components analysis suggests that populations are related primarily on the basis of geography, rather than on the basis of linguistic affinity. This is confirmed in Mantel tests, which show a strong and highly significant partial correlation between genetics and geography but a low, nonsignificant partial correlation between genetics and language. Genetic-barrier analysis also indicates the primacy of geography in the shaping of patterns of variation. These patterns retain a strong signal of expansion from the Near East but also suggest that the demographic history of Europe has been complex and influenced by other major population movements, as well as by linguistic and geographic heterogeneities and the effects of drift.

Y-chromosomal DNA haplotypes in infertile European males carrying Y-microdeletions

We have determined Y-chromosomal DNA haplotypes in 73 infertile European males carrying Y microdeletions and compared them with the haplotypes of 299 infertile males lacking microdeletions. Chromosomes were typed with a set of 11 binary Y markers, which identified 8 haplogroups in the sample. Haplogroup frequencies were compared between 3 microdeletion classes and the non-deleted infertile males. Deletions arise on many different haplotypic backgrounds. No statistically significant differences in frequency were seen, although the small number of AZFa deletions lay predominantly on one branch of the Y haplotype tree.

Y chromosome sequence variation and the history of human populations

Binary polymorphisms associated with the non-recombining region of the human Y chromosome (NRY) preserve the paternal genetic legacy of our species that has persisted to the present, permitting inference of human evolution, population affinity and demographic history. We used denaturing high-performance liquid chromatography (DHPLC; ref. 2) to identify 160 of the 166 bi-allelic and 1 tri-allelic site that formed a parsimonious genealogy of 116 haplotypes, several of which display distinct population affinities based on the analysis of 1062 globally representative individuals. A minority of contemporary East Africans and Khoisan represent the descendants of the most ancestral patrilineages of anatomically modern humans that left Africa between 35,000 and 89,000 years ago.

Demography, genetic diversity, and population relationships among Argentinean Mapuche Indians

Fertility, mortality and migration data from four Mapuche Indian communities located along a 215-km NE-SW linear area in the Province of Río Negro, Argentina, were collated with genetic information furnished by nine blood group systems and by mtDNA haplogroups. The demographic and genetic data indicated a clear dichotomy, which split the four populations into two groups of two. Differing degrees of non-Indian exchanges was probably the main determining factor for this separation. Total genetic variability was very similar in all groups, and the interpopulational variability accounted for only 10% of the total variability. A low prevalence of the Diego(a) antigen among the Mapuche was confirmed. The fact that significant genetic heterogeneity and population clusters were found in such a small territorial region attests to the sensitivity of demographic and genetic approaches in unraveling human history.

Genetic structure of Quechua-speakers of the Central Andes and geographic patterns of gene frequencies in South Amerindian populations

A sample of 141 Quechua-speaking individuals of the population of Tayacaja, in the Peruvian Central Andes, was typed for the following 16 genetic systems: ABO, Rh, MNSs, P, Duffy, AcP1, EsD, GLOI, PGM1, AK, 6-PGD, Hp, Gc, Pi, C3, and Bf. The genetic structure of the population was analyzed in relation to the allele frequencies available for other South Amerindian populations, using a combination of multivariate and multivariable techniques. Spatial autocorrelation analysis was performed independently for 13 alleles to identify patterns of gene flow in South America as a whole and in more specific geographic regions. We found a longitudinal cline for the AcP1*a and EsD*1 alleles which we interpreted as the result of an ancient longitudinal expansion of a putative ancestral population of modern Amerindians. Monmonnier's algorithm, used to identify areas of sharp genetic discontinuity, suggested a clear east-west differentiation of native South American populations, which was confirmed by analysis of the distribution of genetic distances. We suggest that this pattern of genetic structures is the consequence of the independent peopling of western and eastern South America or to low levels of gene flow between these regions, related to different environmental and demographic histories.

Beta-globin gene cluster haplotypes in two North American indigenous populations

Haplotypes derived from five polymorphic restriction sites were determined in 50 Carrier-Sekani and 70 Mvskoke chromosomes, and the results were integrated with those previously obtained for 11 South American Indian populations. Eleven haplotypes were identified in the Mvskokes, while five were observed in the Carrier-Sekani. As in South American natives, haplotype 2 (+----) and 6 (-++ -+) were the most prevalent among the Mvskoke (46% and 30%, respectively). In the Carrier-Sekani, haplotype 2 was also the most common, but haplotype 5 (-+ -++) was somewhat more frequent (18%) than 6 (12%). High heterozygosities, as well as genetic differentiation, were observed among these two North American and two other South American groups (Mapuche and Xavante). They could be due to non-Indian admixture in the Mvskoke and Mapuche, but the findings in the other two populations require some other type of explanation.

Origin of YAP+ lineages of the human Y-chromosome

We screened a total of 841 Y-chromosomes representing 36 human populations of wide geographical distribution for the presence of a Y-specific Alu insert (YAP+ chromosomes). The Alu element was found in 77 cases. We tested 5 biallelic and 8 polyallelic markers in 70 out of the 77 YAP+ chromosomes. We could identify the existence of a hierarchical and chronological structuring of ancestral and derived YAP+ lineages, giving rise to 4 haplogroups, 14 subhaplogroups and 60 haplotypes. Moreover, we propose a monophyletic origin for each one of the YAP+ lineages. Out-of-Africa and out-of-Asia models have been suggested to explain the origin and evolution of ancestral and derived YAP+ elements. We analyze the evidence supporting these two hypotheses, and we conclude that the information available does not allow one to decide between the out-of-Asia or out-of-Africa models.

Variation in short tandem repeats is deeply structured by genetic background on the human Y chromosome

Eleven biallelic polymorphisms and seven short-tandem-repeat (STR) loci mapping on the nonrecombining portion of the human Y chromosome have been typed in men from northwestern Africa. Analysis of the biallelic markers, which represent probable unique events in human evolution, allowed us to characterize the stable backgrounds or haplogroups of Y chromosomes that prevail in this geographic region. Variation in the more rapidly mutating genetic markers (STRs) has been used both to estimate the time to the most recent common ancestor for STR variability within these stable backgrounds and to explore whether STR differentiation among haplogroups still retains information about their phylogeny. When analysis of molecular variance was used to study the apportionment of STR variation among both genetic backgrounds (i.e., those defined by haplogroups) and population backgrounds, we found STR variability to be clearly structured by haplogroups. More than 80% of the genetic variance was found among haplogroups, whereas only 3.72% of the genetic variation could be attributed to differences among populations-that is, genetic variability appears to be much more structured by lineage than by population. This was confirmed when two population samples from the Iberian Peninsula were added to the analysis. The deep structure of the genetic variation in old genealogical units (haplogroups) challenges a population-based perspective in the comprehension of human genome diversity. A population may be better understood as an association of lineages from a deep and population-independent gene genealogy, rather than as a complete evolutionary unit.

Recent male-mediated gene flow over a linguistic barrier in Iberia, suggested by analysis of a Y-chromosomal DNA polymorphism

We have examined the worldwide distribution of a Y-chromosomal base-substitution polymorphism, the T/C transition at SRY-2627, where the T allele defines haplogroup 22; sequencing of primate homologues shows that the ancestral state cannot be determined unambiguously but is probably the C allele. Of 1,191 human Y chromosomes analyzed, 33 belong to haplogroup 22. Twenty-nine come from Iberia, and the highest frequencies are in Basques (11%; n=117) and Catalans (22%; n=32). Microsatellite and minisatellite (MSY1) diversity analysis shows that non-Iberian haplogroup-22 chromosomes are not significantly different from Iberian ones. The simplest interpretation of these data is that haplogroup 22 arose in Iberia and that non-Iberian cases reflect Iberian emigrants. Several different methods were used to date the origin of the polymorphism: microsatellite data gave ages of 1,650, 2,700, 3,100, or 3,450 years, and MSY1 gave ages of 1,000, 2,300, or 2,650 years, although 95% confidence intervals on all of these figures are wide. The age of the split between Basque and Catalan haplogroup-22 chromosomes was calculated as only 20% of the age of the lineage as a whole. This study thus provides evidence for direct or indirect gene flow over the substantial linguistic barrier between the Indo-European and non-Indo-European-speaking populations of the Catalans and the Basques, during the past few thousand years.

Microsatellites provide evidence for Y chromosome diversity among the founders of the New World

Recently, Y chromosome markers have begun to be used to study Native American origins. Available data have been interpreted as indicating that the colonizers of the New World carried a single founder haplotype. However, these early studies have been based on a few, mostly complex polymorphisms of insufficient resolution to determine whether observed diversity stems from admixture or diversity among the colonizers. Because the interpretation of Y chromosomal variation in the New World depends on founding diversity, it is important to develop marker systems with finer resolution. Here we evaluate the hypothesis of a single-founder Y haplotype for Amerinds by using 11 Y-specific markers in five Colombian Amerind populations. Two of these markers (DYS271, DYS287) are reliable indicators of admixture and detected three non-Amerind chromosomes in our sample. Two other markers (DYS199, M19) are single-nucleotide polymorphisms mostly restricted to Native Americans. The relatedness of chromosomes defined by these two markers was evaluated by constructing haplotypes with seven microsatellite loci (DYS388 to 394). The microsatellite backgrounds found on the two haplogroups defined by marker DYS199 demonstrate the existence of at least two Amerind founder haplotypes, one of them (carrying allele DYS199 T) largely restricted to Native Americans. The estimated age and distribution of these haplogroups places them among the founders of the New World.

Ancestral Asian source(s) of new world Y-chromosome founder haplotypes

Haplotypes constructed from Y-chromosome markers were used to trace the origins of Native Americans. Our sample consisted of 2,198 males from 60 global populations, including 19 Native American and 15 indigenous North Asian groups. A set of 12 biallelic polymorphisms gave rise to 14 unique Y-chromosome haplotypes that were unevenly distributed among the populations. Combining multiallelic variation at two Y-linked microsatellites (DYS19 and DXYS156Y) with the unique haplotypes results in a total of 95 combination haplotypes. Contra previous findings based on Y- chromosome data, our new results suggest the possibility of more than one Native American paternal founder haplotype. We postulate that, of the nine unique haplotypes found in Native Americans, haplotypes 1C and 1F are the best candidates for major New World founder haplotypes, whereas haplotypes 1B, 1I, and 1U may either be founder haplotypes and/or have arrived in the New World via recent admixture. Two of the other four haplotypes (YAP+ haplotypes 4 and 5) are probably present because of post-Columbian admixture, whereas haplotype 1G may have originated in the New World, and the Old World source of the final New World haplotype (1D) remains unresolved. The contrasting distribution patterns of the two major candidate founder haplotypes in Asia and the New World, as well as the results of a nested cladistic analysis, suggest the possibility of more than one paternal migration from the general region of Lake Baikal to the Americas.

Characterization of ancestral and derived Y-chromosome haplotypes of New World native populations

We analyze the allelic polymorphisms in seven Y-specific microsatellite loci and a Y-specific alphoid system with 27 variants (alphah I-XXVII), in a total of 89 Y chromosomes carrying the DYS199T allele and belonging to populations representing Amerindian and Na-Dene linguistic groups. Since there are no indications of recurrence for the DYS199C-->T transition, it is assumed that all DYS199T haplotypes derive from a single individual in whom the C-->T mutation occurred for the first time. We identified both the ancestral founder haplotype, 0A, of the DYS199T lineage and seven derived haplogroups diverging from the ancestral one by one to seven mutational steps. The 0A haplotype (5.7% of Native American chromosomes) had the following constitution: DYS199T, alphah II, DYS19/13, DYS389a/10, DYS389b/27, DYS390/24, DYS391/10, DYS392/14, and DYS393/13 (microsatellite alleles are indicated as number of repeats). We analyzed the Y-specific microsatellite mutation rate in 1,743 father-son transmissions, and we pooled our data with data in the literature, to obtain an average mutation rate of.0012. We estimated that the 0A haplotype has an average age of 22,770 years (minimum 13,500 years, maximum 58,700 years). Since the DYS199T allele is found with high frequency in Native American chromosomes, we propose that 0A is one of the most prevalent founder paternal lineages of New World aborigines.

European Y-chromosomal lineages in Polynesians: a contrast to the population structure revealed by mtDNA

We have used Y-chromosomal polymorphisms to trace paternal lineages in Polynesians by use of samples previously typed for mtDNA variants. A genealogical approach utilizing hierarchical analysis of eight rare-event biallelic polymorphisms, seven microsatellite loci, and internal structural analysis of the hypervariable minisatellite, MSY1, has been used to define three major paternal-lineage clusters in Polynesians. Two of these clusters, both defined by novel MSY1 modular structures and representing 55% of the Polynesians studied, are also found in coastal Papua New Guinea. Reduced Polynesian diversity, relative to that in Melanesians, is illustrated by the presence of several examples of identical MSY1 codes and microsatellite haplotypes within these lineage clusters in Polynesians. The complete lack of Y chromosomes having the M4 base substitution in Polynesians, despite their prevalence (64%) in Melanesians, may also be a result of the multiple bottleneck events during the colonization of this region of the world. The origin of the M4 mutation has been dated by use of two independent methods based on microsatellite-haplotype and minisatellite-code diversity. Because of the wide confidence limits on the mutation rates of these loci, the M4 mutation cannot be conclusively dated relative to the colonization of Polynesia, 3,000 years ago. The other major lineage cluster found in Polynesians, defined by a base substitution at the 92R7 locus, represents 27% of the Polynesians studied and, most probably, originates in Europe. This is the first Y-chromosomal evidence of major European admixture with indigenous Polynesian populations and contrasts sharply with the picture given by mtDNA evidence.

Beta-globin gene cluster haplotypes in the Mapuche Indians of Argentina

Haplotypes derived from five polymorphic restriction sites in the beta-globin gene cluster were investigated in 86 chromosomes from the Argentinian Mapuche. These results were integrated with those previously obtained for ten Brazilian Indian tribes. Eight haplotypes were identified, the most frequent being 2 (57%) and 6 (27%). The presence of haplotype 3 in 2% of the Mapuche chromosomes is probably an evidence of admixture with individuals of African ancestry. Due to the high number of haplotypes observed, heterozygosity as measured by the Gini-Simpson index was higher in the Mapuche than in Brazilian Indians. The haplotypic distribution in the Mapuche was also significantly different from those of all Brazilian tribes investigated. This heterogeneity could be at least partially explained by admixture with non-Indian populations.

Network analyses of Y-chromosomal types in Europe, northern Africa, and western Asia reveal specific patterns of geographic distribution

In a study of 908 males from Europe, northern Africa, and western Asia, the variation of four Y-linked dinucleotide microsatellites was analyzed within three "frames" that are defined by mutations that are nonrecurrent, or nearly so. The rapid generation and extinction of new dinucleotide length variants causes the haplotypes within each lineage to diverge from one another. We constructed networks of "adjacent" haplotypes within each frame, by assuming changes of a single dinucleotide unit. Two small and six large networks were obtained, the latter including 94.9% of the sampled Y chromosomes. We show that the phenetic relationships among haplotypes, represented as a network, result largely from common descent and subsequent molecular radiation. The grouping of haplotypes of the same network thus fits an evolutionarily relevant criterion. Notably, this method allows the total diversity within a sample to be partitioned. Networks can be considered optimal markers for population studies, because reliable frequency estimates can be obtained in small samples. We present synthetic maps describing the incidence of different Y-chromosomal lineages in the extant human populations of the surveyed areas. Dinucleotide diversity also was used to infer time intervals for the coalescence of each network.