Evolutionary relationships of human populations on a global scale

Using gene frequency data for 29 polymorphic loci (121 alleles), we conducted a phylogenetic analysis of 26 representative populations from around the world by using the neighbor-joining (NJ) method. We also conducted a separate analysis of 15 populations by using data for 33 polymorphic loci. These analyses have shown that the first major split of the phylogenetic tree separates Africans from non-Africans and that this split occurs with a 100% bootstrap probability. The second split separates Caucasian populations from all other non-African populations, and this split is also supported by bootstrap tests. The third major split occurs between Native American populations and the Greater Asians that include East Asians (mongoloids), Pacific Islanders, and Australopapuans (native Australians and Papua New Guineans), but Australopapuans are genetically quite different from the rest of the Greater Asians. The second and third levels of population splitting are quite different from those of the phylogenetic tree obtained by Cavalli- Sforza et al. (1988), where Caucasians, Northeast Asians, and Ameridians from the Northeurasian supercluster and the rest of non- Africans form the Southeast Asian supercluster. One of the major factors that caused the difference between the two trees is that Cavalli-Sforza et al. used unweighted pair-group method with arithmetic mean (UPGMA) in phylogenetic inference, whereas we used the NJ method in which evolutionary rate is allowed to vary among different populations. Bootstrap tests have shown that the UPGMA tree receives poor statistical support whereas the NJ tree is well supported. Implications that the phylogenetic tree obtained has on the current controversy over the out-of-Africa and the multiregional theories of human origins are discussed.      

Tandem-repeat internal mapping (TRIM) of the involucrin gene: repeat number and repeat-pattern polymorphism within a coding region in human populations.

We have analyzed the human involucrin gene in 41 British African-Caribbeans and 37 British white Caucasians by tandem-repeat internal mapping and DNA sequencing. A point mutation (i.e., Bc) in the last B repeat unit was found in 98.6% of British white Caucasians and in 52.4% of British African-Caribbeans. The distribution of repeat patterns was also different between the two populations. Nine previously unreported repeat pattern alleles, 4 with and 5 without the Bc repeat, have been found, increasing the range of variation in humans to 15 reported repeat patterns, 6 with and 9 without the Caucasian mutation. Three further sequence variations, each occurring in a single individual, were found. The evolutionary significance of variation in the human involucrin gene is discussed.     

Polymorphism due to variable number of repeats in the human involucrin gene.

The coding region of the involucrin gene in higher primates contains a segment consisting of numerous tandem repeats of a 10-codon sequence. The process of repeat addition began in a common ancestor of all higher primates and subsequent repeats were added vectorially. As a result, the principal site of repeat addition has moved in the 3' to 5' direction and the most recently generated repeats (the late region) are close to the 5' end of the segment of repeats. In the human, most of the late region is made up of two different blocks, each consisting of nearly identical repeats. We describe here five polymorphic forms resulting from the addition of differing numbers of repeats to each block. As the variety and nature of the polymorphic alleles are different in different human populations, we postulate that the process of repeat addition is genetically determined.     

Minisatellite MS32 Alleles Show Population Specificity Among Thai,Chinese, and Japanese

Lineages of structurally related alleles at minisatellite MS32 in human populations show considerable differentiation at the continental level. However, the regional specificity of these lineages remains unknown. We now describe the comparison of allele structures in Thai, Han Chinese, and Japanese populations with lineages previously established for North Europeans and Africans. The great majority of alignable Asian alleles showed their closest structural relative in Asia, with few instances of preferential alignment of Asian with European alleles and only one isolated incident showing a best match with an African allele. Further, there was a strong tendency, most marked for Japanese, for Asian alleles to align preferentially with other alleles from the same population, indicating strong regional specificity of allele lineages. This rapidly evolving minisatellite can therefore serve as a lineage marker for exploring recent events in human population history and dissecting population structure at the fine-scale level, as well as being an extremely informative DNA marker for personal identification.     

The involucrin gene of the gibbon: the middle region shared by the hominoids

The evolution of the anthropoid involucrin gene has resulted largely from a process of vectorial addition of short tandem repeats. The coding region of the involucrin gene of the gibbon (Hylobates lar), including the segment of repeats, has been cloned and sequenced, and its repeat structure can now be compared with that of the other hominoids. In the gibbon, as in the others, repeat additions in the past can be assigned to early, middle, and late regions of the present-day segment of repeats. All 10 repeats of the gibbon early region were completed in a common anthropoid ancestor. All 17 repeats of the gibbon middle region were completed in a common hominoid ancestor. After divergence of the gibbon lineage, eight repeats were added to the middle region of the great ape-human lineages. Seven of these are shared by two to four species, according to the order of their divergences from each other. After its divergence, the gibbon lineage added a short species-specific late region. The gibbon also possesses an incomplete repeat just 3' of the early region, the only addition in this region in any hominoid. Comparison of the number of repeats added with the number of nucleotides substituted shows an inconstant relation between the two.     

Patterns of instability of expanded CAG repeats at the ERDA1 locus in general populations.

A highly polymorphic CAG repeat locus, ERDA1, was recently described on human chromosome 17q21.3, with alleles as large as 50-90 repeats and without any disease association in the general population. We have studied allelic distribution at this locus in five human populations and have characterized the mutational patterns by direct observation of 731 meioses. The data show that large alleles (>/=40 CAG repeats) are generally most common in Asian populations, less common in populations of European ancestry, and least common among Africans. We have observed a high intergenerational instability (46. 3%+/-5.1%) of the large alleles. Although the mutation rate is not dependent on parental sex, paternal transmissions have predominantly resulted in contractions, whereas maternal transmissions have yielded expansions. Within this class of large alleles, the mutation rate increases concomitantly with increasing allele size, but the magnitude of repeat size change does not depend on the size of the progenitor allele. Sequencing of specific alleles reveals that the intermediate-sized alleles (30-40 repeats) have CAT/CAC interruptions within the CAG-repeat array. These results indicate that expansion and instability of trinucleotide repeats are not exclusively disease-associated phenomena. The implications of the existence of massively expanded alleles in the general populations are not yet understood.     

Systematic repeat addition at a precise location in the coding region of the involucrin gene of wild mice reveals their phylogeny

The involucrin gene encodes a protein of terminally differentiated keratinocytes. Its segment of repeats, which represents up to 80% of the coding region, is highly polymorphic in mouse strains derived from wild progenitors. Polymorphism includes nucleotide substitutions, but is most strikingly due to the recent addition of a variable number of repeats at a precise location within the segment of repeats. Each mouse taxon examined showed consistent and distinctive patterns of evolution of its variable region: very rapid changes in most M. m. domesticus alleles, slow changes in M. m. musculus, and complete arrest in M. spretus. We conclude that changes in the variable region are controlled by the genetic background. One of the M. m. domesticus alleles (DIK-L), which is of M. m. musculus origin, has undergone a recent repeat duplication typical of M. m. domesticus. This suggests that the genetic background controls repeat duplications through trans-acting factors. Because the repeat pattern differs in closely related murine taxa, involucrin reveals with greater sensitivity than random nucleotide substitutions the evolutionary relations of the mouse and probably of all murids.     

Confirmation of the potential usefulness of two human beta globin pseudogene markers to estimate gene flows to and from sub-Saharan Africans

Two polymorphic sites, -107 and -100 with respect to the "cap" site of the human beta globin pseudogene, recently discovered in our laboratory, turned out to have an ethnically complementary distribution. The first site is polymorphic in Europeans, North Africans, Indians (Hindu), and Oriental Asians, and monomorphic in sub-Saharan Africans. Conversely, the second site is polymorphic in sub-Saharan African populations and monomorphic in the aforementioned populations. Here we report the gene frequencies of these two polymorphic sites in nine additional populations (Egyptians, Spaniards, Japanese, Chinese, Filipinos, Vietnamese, Africans from Togo and from Benin, and Pygmies), confirming their ethnospecificity and, through the analysis of these two markers in Oromo and Amhara of Ethiopia (two mixed populations), their usefulness in genetic admixture studies. Moreover, we studied another marker polymorphic in sub-Saharan African populations only, a TaqI restriction fragment length polymorphism located in the same region as the present markers, demonstrating the absence of linkage disequilibrium between it and the -100 site, so that we can exclude that the information they provide is redundant.     

Human DNA sequence variation in a 6.6-kb region containing the melanocortin 1 receptor promoter.

An approximately 6.6-kb region located upstream from the melanocortin 1 receptor (MC1R) gene and containing its promoter was sequenced in 54 humans (18 Africans, 18 Asians, and 18 Europeans) and in one chimpanzee, gorilla, and orangutan. Seventy-six polymorphic sites were found among the human sequences and the average nucleotide diversity (pi) was 0.141%, one of the highest among all studies of nuclear sequence variation in humans. Opposite to the pattern observed in the MC1R coding region, in the present region pi is highest in Africans (0.136%) compared to Asians (0.116%) and Europeans (0.122%). The distributions of pi, theta, and Fu and Li's F-statistic are nonuniform along the sequence and among continents. The pattern of genetic variation is consistent with a population expansion in Africans. We also suggest a possible phase of population size reduction in non-Africans and purifying selection acting in the middle subregion and parts of the 5' subregion in Africans. We hypothesize diversifying selection acting on some sites in the 5' and 3' subregions or in the MC1R coding region in Asians and Europeans, though we cannot reject the possibility of relaxation of functional constraints in the MC1R gene in Asians and Europeans. The mutation rate in the sequenced region is 1.65 x 10(-9) per site per year. The age of the most recent common ancestor for this region is similar to that for the other long noncoding regions studied to date, providing evidence for ancient gene genealogies. Our population screening and phylogenetic footprinting suggest potentially important sites for the MC1R promoter function.     

Geographic differences in the allele frequencies of the human Y-linked tetranucleotide polymorphism DYS19

We have studied the allele frequency distribution of the microsatellite locus DYS 19 in several populations with different geographical origins worldwide. Three new alleles were found. In addition, remarkable geographic and ethnic differences were observed in the allele frequency profiles and DNA marker (gene) diversity among populations and major ethnic groups. Amerindians showed an overwhelming predominance of the A allele, while in Caucasians the B allele was modal, and in Greater Asians and Africans allele C became predominant. Even within these geographic regions there were significant gradients, as exemplified by the decreasing frequency profile of the B allele from Great Britain over Germany to Slovakia. Thus, DYS 19 emerges as a useful tool for studying the structure and dynamics of human populations.     

The involucrin genes of the white-fronted capuchin and cottontop tamarin: the platyrrhine middle region.

In all anthropoid species, the coding region of the involucrin gene contains a segment of short tandem repeats that were added sequentially, beginning in a common anthropoid ancestor. The involucrin coding region of each of two platyrrhine species, the white-fronted capuchin (Cebus albifrons) and the cottontop tamarin (Saguinus oedipus), has now been cloned and sequenced. These genes share with the genes of the catarrhines the repeats added in the common anthropoid lineage (the early region). After their divergence, the platyrrhines, like the catarrhines, continued to add repeats vectorially 5' of the early region, to form a middle region. The mechanism that was established in the common anthropoid lineage for the addition of repeats at a definite site in the coding region was transmitted to both platyrrhines and catarrhines, enabling each to generate its middle region independently. The process of vectorial repeat addition continued in two platyrrhine sublineages after their divergence from each other.     

A novel tetrameric short tandem repeat located in the 3' flanking region of the human ABO-secretor gene (FUT2) and association between FUT2 and FUT2/01 loci

We found a novel polymorphic short tandem repeat (FUT2/01), 3.8 kb downstream of the coding region of FUT2. Seventeen length and 33 sequence variants were identified in 300 individuals representing three major human populations. Africans (Xhosa) and Europeans were characterized by high microvariation, and Japanese were characterized by a simple repeat structure. All exhibited high haplotype diversity.     

Y-chromosome specific alleles and haplotypes in European and Asian populations: Linkage disequilibrium and geographic diversity

Variation on the Y chromosome may permit our understanding the evolution of the human paternal lineage and male gene flow. This study reports upon the distribution and non random association of alleles at four Y-chromosome specific loci in four populations, three Caucasoid (Italian, Greek and Slav) and one Asian. The markers include insertion/deletion (p12f), point mutation (92R7 and pYαI), and repeat sequence (p21A1) polymorphisms. Our data confirm that the p12f/TaqI 8 kb allele is a Caucasoid marker and that Asians are monomorphic at three of the loci (p12f, 92R7, and pYαI). The alleles at 92R7 and pYαI were found to be in complete disequilibrium in Europeans. Y-haplotype diversity was highly significant between Asians and all three European groups (P < 0.001), but the Greeks and Italians were also significantly different with respect to some alleles and haplotypes (P < 0.02). We find strong evidence that the p12f/TaqI 8 kb allele may have arisen only once, as a deletion event, and, additionally, that the present-day frequency distribution of Y chromosomes carrying the p12f/8 kb allele suggests that it may have been spread by colonising sea-faring peoples from the Near East, possibly the Phoenicians, rather than by expansion of Neolithic farmers into continental Europe. The p12f deletion is the key marker of a unique Y chromosome, found only in Caucasians to date, labelled ‘Mediterranean’ and this further increases the level of Y-chromosome diversity seen among Caucasoids when compared to the other major population groups. Am J Phys Anthropol 104:167–176, 1997. © 1997 Wiley-Liss, Inc.     

Genome-wide SNP typing reveals signatures of population history

Single-nucleotide polymorphism (SNP) arrays have become a popular technology for disease-association studies, but they also have potential for studying the genetic differentiation of human populations. Application of the Affymetrix GeneChip Human Mapping 500K Array Set to a population of 102 individuals representing the major ethnic groups in the United States (African, Asian, European, and Hispanic) revealed patterns of gene diversity and genetic distance that reflected population history. We analyzed allelic frequencies at 388,654 autosomal SNP sites that showed some variation in our study population and 10% or fewer missing values. Despite the small size (23-31 individuals) of each subpopulation, there were no fixed differences at any site between any two subpopulations. As expected from the African origin of modern humans, greater gene diversity was seen in Africans than in either Asians or Europeans, and the genetic distance between the Asian and the European populations was significantly lower than that between either of these two populations and Africans. Principal components analysis applied to a correlation matrix among individuals was able to separate completely the major continental groups of humans (Africans, Asians, and Europeans), while Hispanics overlapped all three of these groups. Genes containing two or more markers with extraordinarily high genetic distance between subpopulations were identified as candidate genes for health differences between subpopulations. The results show that, even with modest sample sizes, genome-wide SNP genotyping technologies have great promise for capturing signatures of gene frequency difference between human subpopulations, with applications in areas as diverse as forensics and the study of ethnic health disparities.     

The involucrin gene of the owl monkey: origin of the early region

A large part of the coding region of the hominoid involucrin gene is of recent origin. This part of the gene, which we have called the modern segment, contains numerous repeats of a sequence of 10 codons, created by multiple duplications some of which consist of 3-12 repeats. We have sequenced two alleles of the involucrin gene in the owl monkey and found that the involucrin gene of this species also possesses a modern segment. By comparing the modern segment of the owl monkey with that of the hominoids, we find that only a part of this segment is shared by the two species. We call this part the early region because it must have originated in a common ancestor of the anthropoids. The rest of the hominoid modern segment does not correspond to any groups of repeats in the owl monkey and was therefore created after divergence of the two lineages. As in the hominoids, the latest additions to the modern segment of the owl monkey have been in its 5' half, which possesses different duplication patterns in the two alleles. Lineage divergences within the anthropoids can be detected at different sites within the modern segment.      

Ethnic differences in the frequency of prostate cancer susceptibility alleles at SRD5A2 and CYP3A4

OBJECTIVES: Ethnic differences in prostate cancer incidence are well documented, with African-Americans having among the highest rates in the world. Ethnic differences in genotypes for genes associated with androgen metabolism including SRD5A2 and CYP3A4 also may exist. The aim of this study was to evaluate differences in these genotypes by ethnicity. METHODS: We studied cancer-free controls representative of four groups: 147 African Americans, 410 Caucasian-Americans, 129 Ghanaians, and 178 Senegalese. PCR-based genotype analysis was undertaken to identify two alleles (V89L, A49T) at SRD5A2 and *1B allele at CYP3A4. RESULTS: Differences were observed for V89L (variant frequency of 30% in Caucasians, 27% in African Americans, 19% in Ghanaians, and 18% in Senegalese, p = 0.002) and were observed for CYP3A4*1B (variant frequencies of 8% in Caucasians, 59% in African Americans, 81% in Ghanaians, and 78% in Senegalese, p = 0.0001). Pooled data combining the present data and previously published data from from Asian, Hispanic, and Arab cancer-free controls showed significant ethnic differences for SRD5A2 and CYP3A4 polymorphisms. Overall, Asians were least likely to have alleles associated with increased prostate cancer risk, while Africans were most likely to have those alleles. CONCLUSIONS: These results suggest that ethnicity-specific differences in genotype frequencies exist for SRD5A2 and CYP3A4. Africans and African-Americans have the highest frequency of those alleles that have previously been associated with increased prostate cancer risk. Future studies should address whether allele frequency differences in part explain differences in prostate cancer incidence in these populations.     

Modulating action of the new polymorphism L436F detected in the<Emphasis Type="Italic"> GLB1</Emphasis> gene of a type-II GM1 gangliosidosis patient

Mutations in the GJB2 gene encoding connexin 26 (Cx26) are a major cause of autosomal recessive and sporadic cases of congenital deafness in most populations. The 235delC mutation of GJB2 is the most frequent known mutation in some east Asian populations, with a carrier frequency of approximately 1%. In order to study the origin of 235delC among east Asians, we analyzed single-nucleotide polymorphisms (SNPs) within the coding region of GJB2 and flanking the 235delC mutation. We observed significant linkage disequilibrium between 235delC and five linked polymorphic markers, suggesting that 235delC arose from a common founder. The detection of 235delC only in east Asians, but not in Caucasians, and the small chromosomal interval of the shared haplotype suggest that 235delC is an ancient mutation that arose after the divergence of Mongoloids and Caucasians. Similarly, the finding that this mutation appears on a single haplotype provides no support for the possibility that recurrent mutation is the explanation for the high frequency of the allele.D. Yan, H.-J. Park and X.M. Ouyang contributed equally to this work     

Population differences in the human functional olfactory repertoire

Olfactory receptors (OR) constitute the molecular basis for the sense of smell. They are encoded by a large multigene family that in humans includes approximately 400 functional genes and approximately 600 putative pseudogenes, distributed on all but two chromosomes. To examine the ethnogeographic variability in the functional chemosensory repertoire, we resequenced 32 OR loci reported to contain a single coding region disruption in seven Caucasians and seven Pygmies. Thirteen of the 32 OR loci were found to have an interrupted coding region in all 28 alleles sampled, seven had an intact form in all the individuals examined, and 12 were polymorphic, segregating both the intact and the null variants. Among the latter loci, the frequency of the null allele was higher in Caucasians than in Pygmies, suggesting that African populations may have a larger repertoire of functional OR genes. Interestingly, when analyzing the entire OR coding regions, we find an excess of high-frequency derived alleles at many loci in the Caucasian sample but less so in the Pygmy sample. Our observations are unlikely to be accounted for by simple demographic models but may be explained by positive selection acting on OR loci in Caucasians.     

The STR polymorphism (AAAAT)n within the intron 1 of the tumor protein 53 (TP53) locus in 17 populations of different ethnic groups of Africa, America, Asia and Europe

The STR (AAAAT)n within intron 1 of the TP53 locus was screened in 17 populations from 3 main ethnic groups: Europeans, Asiatics, and Africans, and from the hybrid population of Costa Rica (1968 samples). Three alleles, 126/7 (bp/copies of the repeat), 131/8 and 136/9 were the most prevalent in all populations. Other alleles rarely reached frequencies of 10% or higher. Observed heterozygosities ranged between 0.351 and 0.829. Patterns of diversity fit well with both the geographic origin of the samples and the history of the populations screened. A statistical test suggests that single-step mutational events have been the main mechanism producing new alleles at this locus. Fixation indexes (R(ST)) for this marker showed an effect of population subdivision on divergence only within the Asiatic group; they were insensitive at the level of major ethnic groups as well as within Africans and within Europeans.