Microsatellite diversity and the demographic history of modern humans

Mutation rates at Y chromosome specific microsatellites

A collaborative work was carried out by the Spanish and Portuguese ISFG Working Group (GEP-ISFG) to estimate Y-STR mutation rates. Seventeen Y chromosome STR loci (DYS19, DYS385, DYS389I and II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS460, DYS461, DYS635 [GATA C4], GATA H4, and GATA A10) were analyzed in a sample of 3,026 father/son pairs. Among 27,029 allele transfers, 54 mutations were observed, with an overall mutation rate across the 17 loci of 1.998 x 10(-3) (95% CI, 1.501 x 10(-3) to 2.606 x 10(-3)). With just one exception, all of the mutations were single-step, and they were observed only once per gametogenesis. Repeat gains were more frequent than losses, longer alleles were found to be more mutable, and the mutation rate seemed to increase with the father's age. Hum Mutat 26(6), 520-528, 2005. (c) 2005 Wiley-Liss, Inc.

Qatari DNA variation at a crossroad of human migrations

Genomic diversity of the Qatari population was investigated by screening 15 autosomal short tandem repeats (STRs). Significant departures from genetic equilibrium were detected at the D13S317, D19S433 and VWA loci, which persisted after applying Bonferroni-type corrections. Gene diversity (GD) values ranged from 0.6851 (TPOX) to 0.8813 (D2S1338), while observed heterozygosity (Ho) oscillated between 0.3388 (D19S433) and 0.8397 (D2S1338). Interestingly, Ho was lower than expected (He) for 14 of the loci analyzed. The information provided by these microsatellite markers was analyzed by means of genetic distances, multidimensional scaling, hierarchical analyses of the molecular variance (AMOVA) and admixture estimations to assess the genetic relationships of Qatar with European, Asian, African and other Middle Eastern populations. The main findings of the study were the genetic uniqueness of the Qatari population, its strong similarity to the United Arab Emirates (UAE) group, a slight genetic differentiation with respect to other Arab populations (Syria and Egypt) and Turkey, and a certain genetic affinity with sub-Saharan African populations. These results are discussed in light of two major issues: the high consanguinity rates characterizing the Qatari population and its strategic geographic position in the Arabian Peninsula close to major migratory routes, an important pivotal contact zone for bidirectional dispersals between Eurasia and Africa.

Deciphering host migrations and origins by means of their microbes

Mitochondrial DNA and microsatellite sequences are powerful genetic markers for inferring the genealogy and the population genetic structure of animals but they have only limited resolution for organisms that display low genetic variability due to recent strong bottlenecks. An alternative source of data for deciphering migrations and origins in genetically uniform hosts can be provided by some of their microbes, if their evolutionary history correlates closely with that of the host. In this review, we first discuss how a variety of viruses, and the bacterium Helicobacter pylori, can be used as genetic tracers for one of the most intensively studied species, Homo sapiens. Then, we review statistical problems and limitations that affect the calculation of particular population genetic parameters for these microbes, such as mutation rates, with particular emphasis on the effects of recombination, selection and mode of transmission. Finally, we extend the discussion to other host-parasite systems and advocate the adoption of an integrative approach to both sampling and analysis.

Y Chromosome Diversity in Brazilians: Switching Perspectives from Slow to Fast Evolving Markers

We have previously shown that the Y chromosomes of 'white' Brazilians have their immediate geographical origin in Europe, with low frequency of sub-Saharan African chromosomes and virtual absence of Amerindian contribution. The typing of slow evolving polymorphisms on the Y chromosome also revealed no differences between Brazilians and Portuguese, the bulk of European immigrants to Brazil, and even among Brazilians from distinct regions of Brazil, the latter being in sharp contrast with mtDNA data. In order to test if the lack of differentiation is a sex-biased and not a marker-biased phenomenon, we decided to study faster evolving Y chromosome markers in samples from Brazil and Portugal previously studied. The population structure revealed by this work confirmed that there were indeed no significant differences between Brazil and Portugal and no population differentiation within the four geographical regions of Brazil, suggesting that this phenomenon is unrelated to the nature of the markers typed. Nevertheless the fast evolving markers did uncover a higher within population diversity in Brazil than Portugal, which could be explained by the input of diverse European Y chromosomes carried by several migration waves to Brazil. Our present data highlight the significance of typing and combining Y markers that evolve according to distinct mutational paces to usefully assess the levels of diversity in a given population, and can be applied in the study of populations derived from distinct geographical origins such as the Brazilians.

Diversity and divergence among the tribal populations of India

Tribal populations of the Indian subcontinent have been of longstanding interest to anthropologists and human geneticists. To investigate the relationship of Indian tribes to Indian castes and continental populations, we analyzed 45 unlinked autosomal STR loci in 9 tribal groups, 8 castes, and 18 populations from Africa, Europe and East Asia. South Indian tribal populations demonstrate low within-population heterozygosity (range: 0.54 - 0.69), while tribal populations sampled further north and east have higher heterozygosity (range: 0.69 - 0.74). Genetic distance estimates show that tribal Indians are more closely related to caste Indians than to other major groups. Between-tribe differentiation is high and exceeds that for eight sub-Saharan African populations (4.8% vs. 3.7%). Telugu-speaking populations are less differentiated than non-Telugu speakers (F(ST): 0.029 vs. 0.079), but geographic distance was not predictive of genetic affinity between tribes. South Indian tribes show significant population structure, and individuals can be clustered statistically into groups that correspond with their tribal affiliation. These results are consistent with high levels of genetic drift and isolation in Indian tribal populations, particularly those of South India, and they imply that these populations may be potential candidates for linkage disequilibrium and association mapping.

Inference of human geographic origins using Alu insertion polymorphisms

The inference of an individual's geographic ancestry or origin can be critical in narrowing the field of potential suspects in a criminal investigation. Most current technologies rely on single nucleotide polymorphism (SNP) genotypes to accomplish this task. However, SNPs can introduce homoplasy into an analysis since they can be identical-by-state. We introduce the use of insertion polymorphisms based on short interspersed elements (SINEs) as a potential alternative to SNPs. SINE polymorphisms are identical-by-descent, essentially homoplasy-free, and inexpensive to genotype using a variety of approaches. Herein, we present results of a blind study using 100 Alu insertion polymorphisms to infer the geographic ancestry of 18 unknown individuals from a variety of geographic locations. Using a Structure analysis of the Alu insertion polymorphism-based genotypes, we were able to correctly infer the geographic affiliation of all 18 unknown human individuals with high levels of confidence. This technique to infer the geographic affiliation of unknown human DNA samples will be a useful tool in forensic genomics.

Microsatellite diversity suggests different histories for Mediterranean and Northern European cattle populations

Based on archaeological evidence, the spread of agropastoralism across Europe followed two main paths: the Danubian route, along which Neolithic farmers expanded north across the central European plains; and the Mediterranean route, where migration occurred along the coast of the Mediterranean sea. Here we examine 20 cattle breeds from the continent and assess the genetic diversity levels and relationships among the breeds using 19 microsatellite markers. Additionally, we show evidence that concords with two distinct cattle migrations from the Near East, and also demonstrate that Mediterranean cattle breeds may have had more recent input from both the Near East and Africa.

Microsatellite data support subpopulation structuring among Basques

Genomic diversity based on 13 short tandem repeat (STR) loci (D3S1358, vWA, FGA, D8S1179, D21S11, D18S51, D5S818, D13S317, D7S820, D16S539, TH01, TPOX, and CSF1PO) is reported for the first time in Basques from the provinces of Guipúzcoa and Navarre (Spain). STR data from previous studies on Basques from Alava and Vizcaya provinces were also examined using hierarchal analysis of molecular variance (AMOVA) and genetic admixture estimations to ascertain whether the Basques are genetically heterogeneous. To assess the genetic position of Basques in a broader geographic context, we conducted phylogenetic analyses based on F(ST) genetic distances [neighbor-joining trees and multidimensional scaling (MDS)] using data compiled in previous publications. The genetic profile of the Basque groups revealed distinctive regional partitioning of short tandem repeat (STR) diversity. Consistent with the above, native Basques clearly segregated from other populations from Europe (including Spain), North Africa, and the Middle East. The main line of genetic discontinuity inferred from the spatial variability of the microsatellite diversity in Basques significantly overlapped the geographic distribution of the Basque language. The genetic heterogeneity among native Basque groups correlates with the peculiar geography of peopling and marital structure in rural Basque zones and with language boundaries resulting from the uneven impact of Romance languages in the different Basque territories.

Influence of language and ancestry on genetic structure of contiguous populations: a microsatellite based study on populations of Orissa

Background

We have examined genetic diversity at fifteen autosomal microsatellite loci in seven predominant populations of Orissa to decipher whether populations inhabiting the same geographic region can be differentiated on the basis of language or ancestry. The studied populations have diverse historical accounts of their origin, belong to two major ethnic groups and different linguistic families. Caucasoid caste populations are speakers of Indo-European language and comprise Brahmins, Khandayat, Karan and Gope, while the three Australoid tribal populations include two Austric speakers: Juang and Saora and a Dravidian speaking population, Paroja. These divergent groups provide a varied substratum for understanding variation of genetic patterns in a geographical area resulting from differential admixture between migrants groups and aboriginals, and the influence of this admixture on population stratification.

Results

The allele distribution pattern showed uniformity in the studied groups with approximately 81% genetic variability within populations. The coefficient of gene differentiation was found to be significantly higher in tribes (0.014) than caste groups (0.004). Genetic variance between the groups was 0.34% in both ethnic and linguistic clusters and statistically significant only in the ethnic apportionment. Although the populations were genetically close (F_ST = 0.010), the contemporary caste and tribal groups formed distinct clusters in both Principal-Component plot and Neighbor-Joining tree. In the phylogenetic tree, the Orissa Brahmins showed close affinity to populations of North India, while Khandayat and Gope clustered with the tribal groups, suggesting a possibility of their origin from indigenous people.

Conclusions

The extent of genetic differentiation in the contemporary caste and tribal groups of Orissa is highly significant and constitutes two distinct genetic clusters. Based on our observations, we suggest that since genetic distances and coefficient of gene differentiation were fairly small, the studied populations are indeed genetically similar and that the genetic structure of populations in a geographical region is primarily influenced by their ancestry and not by socio-cultural hierarchy or language. The scenario of genetic structure, however, might be different for other regions of the subcontinent where populations have more similar ethnic and linguistic backgrounds and there might be variations in the patterns of genomic and socio-cultural affinities in different geographical regions.

Structure, diversity, and evolution of the 45-bp VNTR in intron 5 of the USH1C gene

Usher syndrome type IC is a rare, autosomal recessive sensorineural disorder caused by mutations in the USH1C gene, which encodes a PDZ-domain protein named harmonin. The Acadian-specific 216G-->A mutation in exon 3 and a variant 9-repeat VNTR allele (designated VNTR(t,t)) in intron 5 are in complete linkage disequilibrium. (The usual form of the allele is referred to as VNTR(t).) To gain insight into the structure, diversity, and evolution of the VNTR, we analyzed individuals from seven different populations, as well as nonhuman primates and rodents. The 2-, 3-, and 6-repeat VNTR alleles were the most common. Four novel alleles containing 1, 5, 7, and 10 repeats were detected with frequencies of 0.002, 0.029, 0.005, and 0.001, respectively. The USH1C VNTR region is highly conserved among primates, but not between primates and rodents. Five unrelated individuals had a 3-repeat VNTR(t,t) allele. Haplotype analysis indicates that the 9-repeat VNTR(t,t) and the 3-repeat VNTR(t,t) alleles arose independently. However, the 9-repeat VNTR(t,t) and 6-repeat VNTR(t) alleles shared the same haplotype, suggesting an expansion from 6(t) to 9(t,t).

Genetic variation, classification and 'race'

New genetic data has enabled scientists to re-examine the relationship between human genetic variation and 'race'. We review the results of genetic analyses that show that human genetic variation is geographically structured, in accord with historical patterns of gene flow and genetic drift. Analysis of many loci now yields reasonably accurate estimates of genetic similarity among individuals, rather than populations. Clustering of individuals is correlated with geographic origin or ancestry. These clusters are also correlated with some traditional concepts of race, but the correlations are imperfect because genetic variation tends to be distributed in a continuous, overlapping fashion among populations. Therefore, ancestry, or even race, may in some cases prove useful in the biomedical setting, but direct assessment of disease-related genetic variation will ultimately yield more accurate and beneficial information.

Evidence for gradients of human genetic diversity within and among continents

Genetic variation in humans is sometimes described as being discontinuous among continents or among groups of individuals, and by some this has been interpreted as genetic support for "races." A recent study in which >350 microsatellites were studied in a global sample of humans showed that they could be grouped according to their continental origin, and this was widely interpreted as evidence for a discrete distribution of human genetic diversity. Here, we investigate how study design can influence such conclusions. Our results show that when individuals are sampled homogeneously from around the globe, the pattern seen is one of gradients of allele frequencies that extend over the entire world, rather than discrete clusters. Therefore, there is no reason to assume that major genetic discontinuities exist between different continents or "races."

Polymorphic Alu insertions in six Brazilian African-derived populations

At least 25 African-derived populations (quilombo remnants) are believed to exist in the Ribeira River Valley, located in the southern part of São Paulo State, Brazil. We studied four Alu polymorphic loci (APO, ACE, TPA25, and FXIIIB) in individuals belonging to six quilombo remnants in addition to individuals sampled from the city of São Paulo. The allelic frequencies observed in the quilombo remnants were similar to those previously observed in African-derived populations from Central and North America. Genetic variability indexes (Fst and Gst values) in our quilombos were higher than the reported values for the majority of other populations analyzed for the same kind of markers, but lower than the variability usually observed in Amerindian groups. The observed high degree of genetic differentiation may be due to genetic drift, especially the founder effect. Our results suggest that these populations behave genetically as semi-isolates. The degree of genetic variability within populations was larger than among them, a finding described in other studies. In the neighbor-joining tree, some of the Brazilian quilombos clustered with the African and African-derived populations (São Pedro and Galvão), others with the Europeans (Pilões, Maria Rosa, and Abobral). Pedro Cubas was placed in an isolated branch. Principal component analysis was also performed and confirmed the trends observed in the neighbor-joining tree. Overall, the quilombos showed a higher degree of gene flow than average when compared to other worldwide populations, but similar to other African-derived populations.

Integrating ethics and science in the International HapMap Project

Genomics resources that use samples from identified populations raise scientific, social and ethical issues that are, in many ways, inextricably linked. Scientific decisions about which populations to sample to produce the HapMap, an international genetic variation resource, have raised questions about the relationships between the social identities used to recruit participants and the biological findings of studies that will use the HapMap. The sometimes problematic implications of those complex relationships have led to questions about how to conduct genetic variation research that uses identified populations in an ethical way, including how to involve members of a population in evaluating the risks and benefits posed for everyone who shares that identity. The ways in which these issues are linked is increasingly drawing the scientific and ethical spheres of genomics research closer together.

Haplotype diversity across 100 candidate genes for inflammation, lipid metabolism, and blood pressure regulation in two populations

Recent studies have suggested that a significant fraction of the human genome is contained in blocks of strong linkage disequilibrium, ranging from ~5 to >100 kb in length, and that within these blocks a few common haplotypes may account for >90% of the observed haplotypes. Furthermore, previous studies have suggested that common haplotypes in candidate genes are generally shared across populations and represent the majority of chromosomes in each population. The conclusions drawn from these preliminary studies, however, are based on an incomplete knowledge of the variation in the regions examined. To bridge this gap in knowledge, we have completely resequenced 100 candidate genes in a population of African descent and one of European descent. Although these genes have been well studied because of their medical importance, we demonstrate that a large amount of sequence variation has not yet been described. We also report that the average number of inferred haplotypes per gene, when complete data is used, is higher than in previous reports and that the number and proportion of all haplotypes represented by common haplotypes per gene is variable. Furthermore, we demonstrate that haplotypes shared between the two populations constitute only a fraction of the total number of haplotypes observed and that these shared haplotypes represent fewer of the African-descent chromosomes than was expected from previous studies. Finally, we show that restricting variation discovery to coding regions does not adequately describe all common haplotypes or the true haplotype block structure observed when all common variation is used to infer haplotypes. These data, derived from complete knowledge of genetic variation in these genes, suggest that the haplotype architecture of candidate genes across the human genome is more complex than previously suggested, with important implications for candidate gene and genomewide association studies.

Distinguishing human ethnic groups by means of sequences from Helicobacter pylori: lessons from Ladakh

The history of mankind remains one of the most challenging fields of study. However, the emergence of anatomically modern humans has been so recent that only a few genetically informative polymorphisms have accumulated. Here, we show that DNA sequences from Helicobacter pylori, a bacterium that colonizes the stomachs of most humans and is usually transmitted within families, can distinguish between closely related human populations and are superior in this respect to classical human genetic markers. H. pylori from Buddhists and Muslims, the two major ethnic communities in Ladakh (India), differ in their population-genetic structure. Moreover, the prokaryotic diversity is consistent with the Buddhists having arisen from an introgression of Tibetan speakers into an ancient Ladakhi population. H. pylori from Muslims contain a much stronger ancestral Ladakhi component, except for several isolates with an Indo-European signature, probably reflecting genetic flux from the Near East. These signatures in H. pylori sequences are congruent with the recent history of population movements in Ladakh, whereas similar signatures in human microsatellites or mtDNA were only marginally significant. H. pylori sequence analysis has the potential to become an important tool for unraveling short-term genetic changes in human populations.

Recent developments in genetic data analysis: what can they tell us about human demographic history?

Over the last decade, a number of new methods of population genetic analysis based on likelihood have been introduced. This review describes and explains the general statistical techniques that have recently been used, and discusses the underlying population genetic models. Experimental papers that use these methods to infer human demographic and phylogeographic history are reviewed. It appears that the use of likelihood has hitherto had little impact in the field of human population genetics, which is still primarily driven by more traditional approaches. However, with the current uncertainty about the effects of natural selection, population structure and ascertainment of single-nucleotide polymorphism markers, it is suggested that likelihood-based methods may have a greater impact in the future.

The International HapMap Project

The goal of the International HapMap Project is to determine the common patterns of DNA sequence variation in the human genome and to make this information freely available in the public domain. An international consortium is developing a map of these patterns across the genome by determining the genotypes of one million or more sequence variants, their frequencies and the degree of association between them, in DNA samples from populations with ancestry from parts of Africa, Asia and Europe. The HapMap will allow the discovery of sequence variants that affect common disease, will facilitate development of diagnostic tools, and will enhance our ability to choose targets for therapeutic intervention.

Ethnicity-independent genetic basis of functional psychoses: a genotype-to-phenotype approach

The functional psychoses schizophrenia, schizoaffective disorder, and bipolar illness represent complex clinical syndromes that are characterized by phenotypic heterogeneity. Yet evidence from numerous studies suggests that (1) the prevalence of schizophrenia and bipolar illness is with 1% very similar across ethnicities, and (2) a strong genetic component is involved in the disorders' pathogenesis. Using data from different US-American ethnicities (77 families with a total of 17 unaffected and 170 affected sib pairs; 276 marker loci), we searched for ethnicity-independent oligogenic susceptibility loci for which the between-sib genetic similarity in affected sib pairs deviated from the expected values. Specifically, we addressed the question of the extent to which genetic risk factors and their interactions constitute multigenic inheritance of functional psychoses across populations and might constitute universal targets for treatment. Our novel multivariate genotype-to-phenotype search strategy was based on a genetic similarity function that allowed us to quantify the inter-individual genetic distances d(x(i), x(j)) between the allelic genotype patterns x(i), x(j) of any two subjects i, j with respect to n loci l(1), l(2), em leader l(n). Thus, we were able to assess the between-ethnicity, the within-ethnicity, and the within-family genetic similarities. The problem of ethnicity-independent vulnerability was addressed by treating the Afro-American families as "training" samples, while the non-Afro-American families served as independent "test" samples. We evaluated the between-sib similarities, which were expected to deviate from "0.5" in affected sib pairs if the region of interest contained markers close to vulnerability genes. The reference value "0.5" was derived from the parent-offspring similarities that are always 0.5, irrespective of the affection status of parents and offspring. We found 12 vulnerability loci on chromosomes 1, 4, 5, 6, 13, 14, 18, and 20, that were reproducible across the two samples under comparison and therefore, likely to constitute an ethnicity-independent, oligogenic vulnerability model of functional psychoses. The elevated vulnerability appeared to be unspecific and to act in such a way that exogenous factors become more likely to trigger the onset of psychiatric illnesses.

Genetic variation in gorillas

This review summarizes what is currently known concerning genetic variation in gorillas, on both inter- and intraspecific levels. Compared to the human species, gorillas, along with the other great apes, possess greater genetic variation as a consequence of a demographic history of rather constant population size. Data and hence conclusions from analysis of mitochondrial DNA (mtDNA), the usual means of describing intraspecific patterns of genetic diversity, are limited at this time. An important task for future studies is to determine the degree of confidence with which gorilla mtDNA can be analyzed, in view of the risk that one will inadvertently analyze artifactual rather than genuine sequences. The limited information available from sequences of nuclear genomic segments does not distinguish western from eastern gorillas, and, in comparison with results from the two chimpanzee species, suggests a relatively recent common ancestry for all gorillas. In the near future, the greatest insights are likely to come from studies aimed at genetic characterization of all individual members of social groups. Such studies, addressing topics such as behavior of individuals with kin and non-kin, and the actual success of male reproductive strategies, will provide a link between behavioral and genetic studies of gorillas.

The effective mutation rate at Y chromosome short tandem repeats, with application to human population-divergence time

We estimate an effective mutation rate at an average Y chromosome short-tandem repeat locus as 6.9x10-4 per 25 years, with a standard deviation across loci of 5.7x10-4, using data on microsatellite variation within Y chromosome haplogroups defined by unique-event polymorphisms in populations with documented short-term histories, as well as comparative data on worldwide populations at both the Y chromosome and various autosomal loci. This value is used to estimate the times of the African Bantu expansion, the divergence of Polynesian populations (the Maoris, Cook Islanders, and Samoans), and the origin of Gypsy populations from Bulgaria.

Race, ancestry, and genes: implications for defining disease risk

Geneticists are interested in finding genes associated with disease. Because of widespread health disparities, race is a variable that is often said to be relevant in this context. The idea is that members of a preconceived "race" share common ancestry that may include genetic risk factors. Human variation has been shaped by the long-term processes of population history, and population samples that reflect that history carry statistical information about shared genetic variation or "ancestry." But race is an elusive concept and a term difficult even to define rigorously. Unfortunately, these problems are neither new nor related to recent genetic knowledge. Race is also one of the most politically charged subjects in American life because its associated sociocultural component has notoriously led to categorical treatment that has been misleading and politically misused. There are ways in which the concept of race (whether or not the term is used) can be a legitimate tool in the search for disease-associated genes. But in that context race reflects deeply confounded cultural as well as biological factors, and a careful distinction must be made between race as a statistical risk factor and causal genetic variables.

Haplotypes in the dystrophin DNA segment point to a mosaic origin of modern human diversity

Although Africa has played a central role in human evolutionary history, certain studies have suggested that not all contemporary human genetic diversity is of recent African origin. We investigated 35 simple polymorphic sites and one T(n) microsatellite in an 8-kb segment of the dystrophin gene. We found 86 haplotypes in 1,343 chromosomes from around the world. Although a classical out-of-Africa topology was observed in trees based on the variant frequencies, the tree of haplotype sequences reveals three lineages accounting for present-day diversity. The proportion of new recombinants and the diversity of the T(n) microsatellite were used to estimate the age of haplotype lineages and the time of colonization events. The lineage that underwent the great expansion originated in Africa prior to the Upper Paleolithic (27,000-56,000 years ago). A second group, of structurally distinct haplotypes that occupy a central position on the tree, has never left Africa. The third lineage is represented by the haplotype that lies closest to the root, is virtually absent in Africa, and appears older than the recent out-of-Africa expansion. We propose that this lineage could have left Africa before the expansion (as early as 160,000 years ago) and admixed, outside of Africa, with the expanding lineage. Contemporary human diversity, although dominated by the recently expanded African lineage, thus represents a mosaic of different contributions.

Reconstruction of human evolutionary tree using polymorphic autosomal microsatellites

Allelic frequencies of 182 tri- and tetra-autosomal microsatellites were used to examine phylogenetic relationships among 19 extant human populations. In particular, because the languages of the Basques and Hunza Burusho have been suggested to have an ancient relationship, this study sought to explore the genetic relationship between these two major language isolate populations and to compare them with other human populations. The work presented here shows that the microsatellite allelic diversity and the number of unique alleles were highest in sub-Saharan Africans. Neighbor-joining trees based on genetic distances and principal component analyses separated populations from different continents, and are consistent with an African origin for modern humans. For the first time, with biparentally transmitted markers, the microsatellite tree also shows that the San are the first branch of the human tree before the branch leading to all other Africans. In contrast to an earlier study, these results provided no evidence of a genetic relationship among the two language isolate groups. Genetic relationships, as ascertained by these microsatellites, are dictated primarily by geographic proximity rather than by remote linguistic origin, Mantel test, R(0) = 0.484, g = 3.802 (critical g value = 1.645; P = 0.05).

Allelic variation at alcohol metabolism genes ( ADH1B, ADH1C, ALDH2) and alcohol dependence in an American Indian population

Enzymes encoded by two gene families, alcohol dehydrogenase ( ADH) and aldehyde dehydrogenase ( ALDH), mediate alcohol metabolism in humans. Allelic variants have been identified that alter metabolic rates and influence risk for alcoholism. Specifically, ADH1B*47His (previously ADH2-2) and ALDH2-2 have been shown to confer protection against alcoholism, presumably through accumulation of acetaldehyde in the blood and a resultant 'flushing response' to alcohol consumption. In the current study, variants at ADH1B (previously ADH2), ADH1C (previously ADH3), and ALDH2 were assayed in DNA extracts from participants belonging to a Southwest American Indian tribe ( n=490) with a high prevalence of alcoholism. Each subject underwent a clinical interview for diagnosis of alcohol dependence, as well as evaluation of intermediate phenotypes such as binge drinking and flushing response to alcohol consumption. Detailed haplotypes were constructed and tested against alcohol dependence and related intermediate phenotypes using both association and linkage analysis. ADH and ALDH variants were also assayed in three Asian and one African population (no clinical data) in order to provide an evolutionary context for the haplotype data. Both linkage and association analysis identified several ADH1C alleles and a neighboring microsatellite marker that affected risk of alcohol dependence and were also related to binge drinking. These data strengthen the support for ADH as a candidate locus for alcohol dependence and suggest further productive study.

The global distribution of length polymorphisms of the promoters of the glucuronosyltransferase 1 gene (UGT1A1): hematologic and evolutionary implications

The promoter region of the UDP glucuronosyltransferase 1 gene (UGT1A1) contains a run of thymine-adenine (TA) repeats, usually six (TA)(6). As well as its relationship to Gilbert's syndrome, homozygosity for the extended sequence, (TA)(7) (TA)(7), has been found to be an important risk factor for hyperbilirubinemia and gallstones in patients with hemoglobin E-beta-thalassemia and other intermediate forms of beta thalassemia. To assess the importance of this polymorphism in these common disorders a wide-scale population study of the relative frequency of the size alleles of the UGT1A1 promoter has been carried out. Homozygosity for the (TA)(7) allele occurs in 10-25% of the populations of Africa and the Indian subcontinent, with a variable frequency in Europe. It occurs at a much lower frequency in Southeast Asia, Melanesia, and the Pacific Islands, ranging from 0 to 5%. African populations show a much greater diversity of length alleles than other populations. These findings define those populations with a high frequency of hemoglobin E-beta-thalassemia and related disorders that are at increased risk for hyperbilirubinemia and gall bladder disease and provide evolutionary insights into how these polymorphisms have arisen and are so unequally distributed among human populations.

Genetic variation among world populations: inferences from 100 Alu insertion polymorphisms

We examine the distribution and structure of human genetic diversity for 710 individuals representing 31 populations from Africa, East Asia, Europe, and India using 100 Alu insertion polymorphisms from all 22 autosomes. Alu diversity is highest in Africans (0.349) and lowest in Europeans (0.297). Alu insertion frequency is lowest in Africans (0.463) and higher in Indians (0.544), E. Asians (0.557), and Europeans (0.559). Large genetic distances are observed among African populations and between African and non-African populations. The root of a neighbor-joining network is located closest to the African populations. These findings are consistent with an African origin of modern humans and with a bottleneck effect in the human populations that left Africa to colonize the rest of the world. Genetic distances among all pairs of populations show a significant product-moment correlation with geographic distances (r = 0.69, P < 0.00001). F(ST), the proportion of genetic diversity attributable to population subdivision is 0.141 for Africans/E. Asians/Europeans, 0.047 for E. Asians/Indians/Europeans, and 0.090 for all 31 populations. Resampling analyses show that approximately 50 Alu polymorphisms are sufficient to obtain accurate and reliable genetic distance estimates. These analyses also demonstrate that markers with higher F(ST) values have greater resolving power and produce more consistent genetic distance estimates.

Inferring recent human phylogenies using forensic STR technology

STR loci are characterized by extremely high mutation rates and thus, high levels of length polymorphism both within and among populations. In addition, much of the observed variation is believed to be nearly selectively neutral. Because of these features, STRs are ideal markers for genetic mapping, intra-species phylogenetic reconstructions and forensic analysis. In the present study, we investigate the application of five STR loci (CS1PO, TH01, TPOX, FGA and vWA) routinely used in forensic analysis for delineating the phylogenetic relationships of 10 human populations representing the three major racial groups (African-Caribbean, Croatian from the island of Hvar, East Asian, Han Chinese, Italian, Japanese, Portuguese, UK Caucasian, US Caucasian and Zimbabwe). The resulting tree topology exhibited strong geographic and racial partitioning consistent with that obtained with mtDNA haplotypes, Y-chromosome markers, SNPs, PAIs (polymorphic Alu insertions) as well as classic genetic polymorphisms. These findings suggest that forensic STR loci may be particularly powerful tools and provide the necessary fine resolution for the reconstruction of recent human evolutionary history.

Differentiation of Caucasians and Chinese at bone mass candidate genes: implication for ethnic difference of bone mass

Bone mineral density (BMD) is an important risk factor for osteoporosis and has strong genetic determination. While average BMD differs among major ethnic groups, several important candidate genes have been shown to underlie BMD variation within populations of the same ethnicity. To investigate whether important candidate genes may contribute to ethnic differences in BMD, we studied the degree of genetic differentiation among several important candidate genes between two major ethnic groups: Caucasians and Chinese. The genetic variability of these two populations (1131 randomly selected individuals) was studied at six restriction sites exhibiting polymorphisms of five important candidate genes for BMD: the BsaHI polymorphism of the calcium-sensing receptor (CASR) gene, the SacI polymorphism of the alpha2HS-glycoprotein (AHSG) gene, the PvuII and XbaI polymorphisms of the estrogen receptor alpha (ESR1) gene, the ApaI polymorphism of the vitamin D receptor (VDR) gene, and the BstBI polymorphism of the parathyroid hormone (PTH) gene. The two ethnic groups showed significant allelic and genotypic differentiation of all the polymorphisms studied. The mean FST was 0.103, which significantly differed from zero (P < 0.01). The Chinese population had lower mean heterozygosity (0.331) than the Caucasian one (0.444); the CASR-BsaHI and PTH-BstBI polymorphisms contributed most significantly to this difference. Analysis of the intra- and inter-population variability suggests that various types of natural selection may affect the observed patterns of variation at some loci. If some of the candidate genes we studied indeed underlie variation in BMD, their population differentiation revealed here between ethnic groups may contribute to understanding ethnic difference in BMD.

Features of evolution and expansion of modern humans, inferred from genomewide microsatellite markers

We study data on variation in 52 worldwide populations at 377 autosomal short tandem repeat loci, to infer a demographic history of human populations. Variation at di-, tri-, and tetranucleotide repeat loci is distributed differently, although each class of markers exhibits a decrease of within-population genetic variation in the following order: sub-Saharan Africa, Eurasia, East Asia, Oceania, and America. There is a similar decrease in the frequency of private alleles. With multidimensional scaling, populations belonging to the same major geographic region cluster together, and some regions permit a finer resolution of populations. When a stepwise mutation model is used, a population tree based on TD estimates of divergence time suggests that the branches leading to the present sub-Saharan African populations of hunter-gatherers were the first to diverge from a common ancestral population (approximately 71-142 thousand years ago). The branches corresponding to sub-Saharan African farming populations and those that left Africa diverge next, with subsequent splits of branches for Eurasia, Oceania, East Asia, and America. African hunter-gatherer populations and populations of Oceania and America exhibit no statistically significant signature of growth. The features of population subdivision and growth are discussed in the context of the ancient expansion of modern humans.

Variation at 16 STR loci in Rwandans (Hutu) and implications on profile frequency estimation in Bantu-speakers

A data set of 16 autosomal STRs (the 13 CODIS loci plus HumCD4, HumFES, HumF13A1) was obtained in a sample of 52 unrelated Hutus from Rwanda. Genotypes at all loci met Hardy-Weinberg expectations with the exception of HumCSF1PO. No significant evidence of association across alleles at independent loci was obtained. Statistical parameters demonstrated the forensic usefulness of the analysed systems (combined PE=0.9999996, combined PD=1:2.27 x 10(18)). Pairwise comparisons showed that the Hutu gene pool differs substantially from that of other Bantu-speaking populations suggesting the use of ethnic-specific population databases in forensic casework analysis. The introduction of a non-negligible bias was confirmed by calculating the differences between multiple-locus profile frequencies of western and eastern Bantoids using local and non-local reference databases.

Extensive linkage disequilibrium, a common 16.7-kilobase deletion, and evidence of balancing selection in the human protocadherin alpha cluster

Regions of extensive linkage disequilibrium (LD) appear to be a common feature of the human genome. However, the mechanisms that maintain these regions are unknown. In an effort to understand whether gene density contributes to LD, we determined the degree of promoter sequence variation in a large tandem-arrayed gene family, the human protocadherin alpha cluster, on chromosome 5. These genes are expressed at synaptic junctions in the developing brain and the adult brain and may be involved in the determination of synaptic complexity. We sequenced the promoters of all 13 alpha protocadherin genes in 96 European Americans and identified polymorphisms in the promoters alpha 1, alpha 3, alpha 4, alpha 5, alpha 7, alpha 9, alpha 11, and alpha 13. In these promoters, 11 common SNPs are in extensive LD, forming two 48-kb haplotypes of equal frequency, in this population, that extend from the alpha1 through alpha 7 genes. We sequenced these promoters in East Asians and African Americans, and we estimated haplotype frequencies and calculated LD statistics for all three populations. Our results indicate that, although extensive LD is an ancient feature of the alpha cluster, it has eroded over time. SNPs 3' of alpha 7 are involved in ancestral recombination events in all populations, and overall alpha-cluster LD is reduced in African Americans. We obtained significant positive values for Tajima's D test for all alpha promoter SNPs in Europeans (D=3.03) and East Asians (D=2.64), indicating an excess of intermediate-frequency variants, which is a signature of balancing selection. We also discovered a 16.7-kb deletion that truncates the alpha 8 gene and completely removes the alpha 9 and alpha 10 genes. This deletion appears in unaffected individuals from multiple populations, suggesting that a reduction in protocadherin gene number is not obviously deleterious.

Structural analysis of insulin minisatellite alleles reveals unusually large differences in diversity between Africans and non-Africans

The insulin minisatellite (INS VNTR) associates with susceptibility to a variety of diseases. We have developed a high-resolution system for analyzing variant repeat distributions applicable to all known minisatellite alleles, irrespective of size, which allows lineages of related alleles to be identified. This system has previously revealed extremely low structural diversity in the minisatellite among northern Europeans from the United Kingdom, with all alleles belonging to one of only three highly diverged lineages called "I," "IIIA," and "IIIB." To explore the origins of this remarkably limited lineage diversity, we have characterized an additional 780 alleles from three non-African and three African populations. In total, 22 highly diverged lineages were identified, with structural intermediates absent from extant populations, suggesting a bottleneck within the ancestry of all humans. The difference between levels of diversity in Africans and non-Africans is unusually large, with all 22 lineages identified in Africa compared with only three lineages seen not only in the United Kingdom but also in the other non-African populations. We also find evidence for overrepresentation of lineage I chromosomes in non-Africans. These data are consistent with a common out-of-Africa origin and an unusually tight bottleneck within the ancestry of all non-African populations, possibly combined with differential and positive selection for lineage I alleles in non-Africans. The important implications of these data for future disease-association studies are discussed.

Evidence for balancing selection from nucleotide sequence analyses of human G6PD

Glucose-6-phosphate dehydrogenase (G6PD) mutations that result in reduced enzyme activity have been implicated in malarial resistance and constitute one of the best examples of selection in the human genome. In the present study, we characterize the nucleotide diversity across a 5.2-kb region of G6PD in a sample of 160 Africans and 56 non-Africans, to determine how selection has shaped patterns of DNA variation at this gene. Our global sample of enzymatically normal B alleles and A, A-, and Med alleles with reduced enzyme activities reveals many previously uncharacterized silent-site polymorphisms. In comparison with the absence of amino acid divergence between human and chimpanzee G6PD sequences, we find that the number of G6PD amino acid polymorphisms in human populations is significantly high. Unlike many other G6PD-activity alleles with reduced activity, we find that the age of the A variant, which is common in Africa, may not be consistent with the recent emergence of severe malaria and therefore may have originally had a historically different adaptive function. Overall, our observations strongly support previous genotype-phenotype association studies that proposed that balancing selection maintains G6PD deficiencies within human populations. The present study demonstrates that nucleotide sequence analyses can reveal signatures of both historical and recent selection in the genome and may elucidate the impact that infectious disease has had during human evolution.

DNA sequence variation in a 3.7-kb noncoding sequence 5' of the CYP1A2 gene: implications for human population history and natural selection

CYP1A2 is a cytochrome P450 gene that is involved in human physiological responses to a variety of drugs and toxins. To investigate the role of population history and natural selection in shaping genetic diversity in CYP1A2, we sequenced a 3.7-kb region 5' from CYP1A2 in a diverse collection of 113 individuals from three major continental regions of the Old World (Africa, Asia, and Europe). We also examined sequences in the 90-member National Institutes of Health DNA Polymorphism Discovery Resource (PDR). Eighteen single-nucleotide polymorphisms (SNPs) were found. Most of the high-frequency SNPs found in the Old World sample were also found in the PDR sample. However, six SNPs were detected in the Old World sample but not in the PDR sample, and two SNPs found in the PDR sample were not found in the Old World sample. Most pairs of SNPs were in complete linkage disequilibrium with one another, and there was no indication of a decline of disequilibrium with physical distance in this region. The average +/- SD nucleotide diversity in the Old World sample was 0.00043+/-0.00026. The African population had the highest level of nucleotide diversity and the lowest level of linkage disequilibrium. Two distinct haplotype clusters with broadly overlapping geographical distributions were present. Of the 17 haplotypes found in the Old World sample, 12 were found in the African sample, 8 were found in Indians, 5 were found in non-Indian Asians, and 5 were found in Europeans. Haplotypes found outside Africa were mostly a subset of those found within Africa. These patterns are all consistent with an African origin of modern humans. Seven SNPs were singletons, and the site-frequency spectrum showed a significant departure from neutral expectations, suggesting population expansion and/or natural selection. Comparison with outgroup species showed that four derived SNPs have achieved high (>0.90) frequencies in human populations, a trend consistent with the action of positive natural selection. These patterns have a number of implications for disease-association studies in CYP1A2 and other genes.

Is a Multivariate Consensus Representation of Genetic Relationships Among Populations Always Meaningful?

To determine the relationships among closely related populations or species, two methods are commonly used in the literature: phylogenetic reconstruction or multivariate analysis. The aim of this article is to assess the reliability of multivariate analysis. We describe a method that is based on principal component analysis and Mantel correlations, using a two-step process: The first step consists of a single-marker analysis and the second step tests if each marker reveals the same typology concerning population differentiation. We conclude that if single markers are not congruent, the compromise structure is not meaningful. Our model is not based on any particular mutation process and it can be applied to most of the commonly used genetic markers. This method is also useful to determine the contribution of each marker to the typology of populations. We test whether our method is efficient with two real data sets based on microsatellite markers. Our analysis suggests that for closely related populations, it is not always possible to accept the hypothesis that an increase in the number of markers will increase the reliability of the typology analysis.

Distribution of CGG repeats and FRAXAC1/DXS548 alleles in South American populations

In order to assess the molecular variability related to fragile X (FMR1 locus), we investigated the distribution of CGG repeats and DXS548/FRAXAC1 haplotypes in normal South American populations of different ethnic backgrounds. Special attention was given to Amerindian Wai-Wai (Northern Brazil) and Ache (Paraguay), as well as to Brazilian isolated communities of African ancestry, the remnants of quilombos. Comparison of samples from quilombos, Amerindians, and the ethnically mixed, but mainly European-derived population of São Paulo revealed that the 30-copy allele of the fragile X gene is the most frequent in all groups. A second peak at 20 repeats was present in the population of São Paulo only, confirming this as a European peculiarity. The distribution of DXS548 and FRAXAC1 alleles led to a high expected heterozygosity in African Brazilians, followed by that observed in the population of São Paulo. Amerindians showed the lowest diversity in CGG repeats and DXS548/FRAXAC1 haplotypes. Some rare alleles, for example, the 148-bp (FRAXAC1) or 200-bp (DXS548) variants, which seem to be almost absent in Europe, occurred in higher frequencies among African Brazilians. This suggests a general trend for higher genetic diversity among Africans; these rarer alleles could be African in origin and would have been lost or possibly were not present in the groups that gave rise to the Europeans.

A strong signature of balancing selection in the 5' cis-regulatory region of CCR5

CCR5 encodes a cell surface chemokine receptor molecule that serves as the principal coreceptor, with CD4, for HIV-type 1 (HIV-1). Varied HIV-1 susceptibility and time to progression to AIDS have been associated with polymorphisms in CCR5. Many of these polymorphisms are located in the 5' cis-regulatory region of CCR5, suggesting that it may have been a target of natural selection. We characterized CCR5 sequence variation in this region in 400 chromosomes from worldwide populations and compared it to a genome-wide analysis of 100 Alu polymorphisms typed in the same populations. Variation was substantially higher than expected and characterized by an excess of intermediate-frequency alleles. A genealogy of CCR5 haplotypes had deep branch lengths despite markedly little differentiation among populations. This finding suggested a deviation from neutrality not accounted for by population structure, which was confirmed by tests for natural selection. These results are strong evidence that balancing selection has shaped the pattern of variation in CCR5 and suggest that HIV-1 resistance afforded by CCR5 5' cis-regulatory region haplotypes may be the consequence of adaptive changes to older pathogens.

Genetic perspectives on human origins and differentiation

This is a review of genetic evidence about the ancient demography of the ancestors of our species and about the genesis of worldwide human diversity. The issue of whether or not a population size bottleneck occurred among our ancestors is under debate among geneticists as well as among anthropologists. The bottleneck, if it occurred, would confirm the Garden of Eden (GOE) model of the origin of modern humans. The competing model, multiregional evolution (MRE), posits that the number of human ancestors has been large, occupying much of the temperate Old World for the last two million years. While several classes of genetic marker seem to contain a strong signal of demographic recovery from a small number of ancestors, other nuclear loci show no such signal. The pattern at these loci is compatible with the existence of widespread balancing selection in humans. The study of human diversity at (putatively) neutral genetic marker loci has been hampered since the beginning by ascertainment bias since they were discovered in Europeans. The high levels of polymorphism at microsatellite loci means that they are free of this bias. Microsatellites exhibit a clear almost linear diversity gradient away from Africa, so that New World populations are approximately 15% less diverse than African populations. This pattern is not compatible with a model of a single large population expansion and colonization of most of the Earth by our ancestors but suggests, instead, gradual loss of diversity in successive colonization bottlenecks as our species grew and spread.

Categorization of humans in biomedical research: genes, race and disease

A debate has arisen regarding the validity of racial/ethnic categories for biomedical and genetic research. An epidemiologic perspective on the issue of human categorization in biomedical and genetic research strongly supports the continued use of self-identified race and ethnicity.

Genetic sampling error of distance (delta(mu))2 and variation in mutation rate among microsatellite loci

An expression is obtained for the time-dependent variance of the microsatellite genetic distance (delta(mu))2 when the mutation rate is allowed to vary randomly among loci. An estimator is presented for the coefficient of variation, C(w), in the mutation rate. Estimated values of C(w) from genetic distances between African and non-African populations were less than 100%. Caveats to this conclusion are discussed.

Demographic history of Diadema antillarum, a keystone herbivore on Caribbean reefs

The sea urchin Diadema antillarum was the most important herbivore on Caribbean reefs until 1983, when mass mortality reduced its populations by more than 97%. Knowledge of its past demography is essential to reconstruct reef ecology as it was before human impact, which has been implicated as having caused high pre-mortality Diadema abundance. To determine the history of its population size, we sequenced the ATPase 6 and 8 region of mitochondrial DNA from populations in the Caribbean and in the eastern Atlantic (which was not affected by the mass mortality), as well as from the eastern Pacific D. mexicanum. The Caribbean population harbours an order of magnitude more molecular diversity than those of the eastern Pacific or the eastern Atlantic and, despite the recent mass mortality, its DNA sequences bear the genetic signature of a previous population expansion. By estimating mutation rates from divergence between D. antillarum and D. mexicanum, that were separated at a known time by the Isthmus of Panama, and by using estimates of effective population size derived from mismatch distributions and a maximum likelihood coalescence algorithm, we date the expansion as having occurred no more recently than 100 000 years before the present. Thus, Diadema was abundant in the Caribbean long before humans could have affected ecological processes; the genetic data contain no evidence of a recent, anthropogenically caused, population increase.

The forensic DNA implications of genetic differentiation between endogamous communities

In many indigenous minority populations, and among migrants from Asian and African populations now resident in western Europe, North America and Australia, there is a strong tradition of endogamy and a preference for consanguineous unions. These marriage practices can result in F(ST) values greatly in excess of the maximum value (0.01) currently recommended for forensic DNA purposes under guidelines established by the National Research Council (NRC) of the USA. To examine the possible extent of deviation from this accepted norm, three co-resident Pakistani communities were studied using 10 autosomal dinucleotide markers and six tetranucleotide markers on the Y-chromosome. The mean population subdivision coefficient (FST) value was 0.13 for the autosomal loci, and Y-chromosome loci exhibited even stronger differentiation with unique alleles identified in all three communities. The data indicate that even when sub-populations are virtually indistinguishable in terms of anthropology, geography, ethnicity or culture, they may still exhibit major genetic differentiation. Where significant population stratification is known to exist, more detailed genetic databases should be developed for forensic DNA purposes, based on reference data from each of the appropriate sub-populations and not on random or combined samples.

The HUMFIBRA (FGA) polymorphism in an italian population and a world-wide frequency distribution analysis

OBJECTIVE

The study investigated the highly polymorphic HumFGA short tandem repeat in a sample of 219 unrelated and native individuals from Bologna, and analysed a complete database of FGA allele frequency distributions in 57 world-wide populations collected from the literature.

METHOD

The HumFGA polymorphism was screened by automated fluorescence analysis of PCR-amplified labelled sample fragments performed with an ABI PRISM 310 Genetic Analyser. Genetic distances (Dsw, delta mu2 and Fst) between populations were computed with the MSAT.2 program. Non-metric multidimensional scaling (nmMDS) and neighbour-joining trees (NJTs) were used to investigate patterns of population affinities. Correspondence analysis of the genetic relationships among populations was also performed.

MAIN RESULTS AND CONCLUSIONS

The FGA microsatellite locus is a population marker with a high degree of polymorphism throughout the world. Fourteen HumFGA alleles, ranging in size from 18 to 29 repeats, were identified and sequenced in the Bologna population. The sample was in Hardy-Weinberg equilibrium and had a heterozygosity value of 0.86. Results obtained from the multivariate analyses were consistent in showing great similarity among Europeans. The few African populations investigated were characterized by an even higher level of polymorphism, probably related to the ancient peopling of that continent.

Comparison of DNA and protein polymorphisms between humans and chimpanzees

To examine the nucleotide diversity at silent (synonymous + intron + untranslated) and non-silent (nonsynonymous) sites in chimpanzees and humans, genes at six nuclear loci from two chimpanzees were sequenced. The average silent diversity was 0.19%, which was significantly higher than that in humans (0.05%). This observation suggests a significantly larger effective population size and a higher extent of neutral polymorphism in chimpanzees than in humans. On the other hand, the non-silent nucleotide diversity is similar in both species, resulting in a larger fraction of neutral mutations at non-silent sites in humans than in chimpanzees. Other types of polymorphism data were collected from the literature or databases to examine whether or not they are consistent with the nuclear DNA sequence polymorphism observed here. The nucleotide diversity at both silent and non-silent sites in mitochondrial (mt) DNA genes was compatible with that of the nuclear genes. Microsatellite loci showed a similar high extent of heterozygosity in both species, perhaps due to the combined effect of a high mutation rate and a recent population expansion in humans. At protein loci, humans are more heterozygous than chimpanzees, and the estimated fraction of neutral alleles in humans (0.84) is much larger than that in chimpanzees (0.26). These data show that the neutral fraction in non-silent changes is relatively large in the human population. This difference may be due to a relaxation of the functional constraint against proteins in the human lineage. To evaluate this possibility, it will be necessary to examine nucleotide sequences in relation to the physiological or biochemical properties of proteins.

Genetic evidence on the origins of Indian caste populations

The origins and affinities of the approximately 1 billion people living on the subcontinent of India have long been contested. This is owing, in part, to the many different waves of immigrants that have influenced the genetic structure of India. In the most recent of these waves, Indo-European-speaking people from West Eurasia entered India from the Northwest and diffused throughout the subcontinent. They purportedly admixed with or displaced indigenous Dravidic-speaking populations. Subsequently they may have established the Hindu caste system and placed themselves primarily in castes of higher rank. To explore the impact of West Eurasians on contemporary Indian caste populations, we compared mtDNA (400 bp of hypervariable region 1 and 14 restriction site polymorphisms) and Y-chromosome (20 biallelic polymorphisms and 5 short tandem repeats) variation in approximately 265 males from eight castes of different rank to approximately 750 Africans, Asians, Europeans, and other Indians. For maternally inherited mtDNA, each caste is most similar to Asians. However, 20%-30% of Indian mtDNA haplotypes belong to West Eurasian haplogroups, and the frequency of these haplotypes is proportional to caste rank, the highest frequency of West Eurasian haplotypes being found in the upper castes. In contrast, for paternally inherited Y-chromosome variation each caste is more similar to Europeans than to Asians. Moreover, the affinity to Europeans is proportionate to caste rank, the upper castes being most similar to Europeans, particularly East Europeans. These findings are consistent with greater West Eurasian male admixture with castes of higher rank. Nevertheless, the mitochondrial genome and the Y chromosome each represents only a single haploid locus and is more susceptible to large stochastic variation, bottlenecks, and selective sweeps. Thus, to increase the power of our analysis, we assayed 40 independent, biparentally inherited autosomal loci (1 LINE-1 and 39 Alu elements) in all of the caste and continental populations (approximately 600 individuals). Analysis of these data demonstrated that the upper castes have a higher affinity to Europeans than to Asians, and the upper castes are significantly more similar to Europeans than are the lower castes. Collectively, all five datasets show a trend toward upper castes being more similar to Europeans, whereas lower castes are more similar to Asians. We conclude that Indian castes are most likely to be of proto-Asian origin with West Eurasian admixture resulting in rank-related and sex-specific differences in the genetic affinities of castes to Asians and Europeans.

Genetic Evidence on the Origins of Indian Caste Populations

The origins and affinities of the ∼1 billion people living on the subcontinent of India have long been contested. This is owing, in part, to the many different waves of immigrants that have influenced the genetic structure of India. In the most recent of these waves, Indo-European-speaking people from West Eurasia entered India from the Northwest and diffused throughout the subcontinent. They purportedly admixed with or displaced indigenous Dravidic-speaking populations. Subsequently they may have established the Hindu caste system and placed themselves primarily in castes of higher rank. To explore the impact of West Eurasians on contemporary Indian caste populations, we compared mtDNA (400 bp of hypervariable region 1 and 14 restriction site polymorphisms) and Y-chromosome (20 biallelic polymorphisms and 5 short tandem repeats) variation in ∼265 males from eight castes of different rank to ∼750 Africans, Asians, Europeans, and other Indians. For maternally inherited mtDNA, each caste is most similar to Asians. However, 20%–30% of Indian mtDNA haplotypes belong to West Eurasian haplogroups, and the frequency of these haplotypes is proportional to caste rank, the highest frequency of West Eurasian haplotypes being found in the upper castes. In contrast, for paternally inherited Y-chromosome variation each caste is more similar to Europeans than to Asians. Moreover, the affinity to Europeans is proportionate to caste rank, the upper castes being most similar to Europeans, particularly East Europeans. These findings are consistent with greater West Eurasian male admixture with castes of higher rank. Nevertheless, the mitochondrial genome and the Y chromosome each represents only a single haploid locus and is more susceptible to large stochastic variation, bottlenecks, and selective sweeps. Thus, to increase the power of our analysis, we assayed 40 independent, biparentally inherited autosomal loci (1 LINE-1 and 39 Aluelements) in all of the caste and continental populations (∼600 individuals). Analysis of these data demonstrated that the upper castes have a higher affinity to Europeans than to Asians, and the upper castes are significantly more similar to Europeans than are the lower castes. Collectively, all five datasets show a trend toward upper castes being more similar to Europeans, whereas lower castes are more similar to Asians. We conclude that Indian castes are most likely to be of proto-Asian origin with West Eurasian admixture resulting in rank-related and sex-specific differences in the genetic affinities of castes to Asians and Europeans.