Spatial and Temporal Distribution of the Neutral Polymorphisms in the Last ZFX Intron: Analysis of the Haplotype Structure and Genealogy

Corresponding author: Damian Labuda, Centre de Recherche, Hôpital Sainte-Justine, 3175 Côte-Sainte Catherine, Montreal, Quebec H3T 1C5 Canada., labuda{at}ere.umontreal.ca (E-mail)

Communicating editor: A. G. CLARK

	ABSTRACT

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

With 10 segregating sites (simple nucleotide polymorphisms) in the last intron (1089 bp) of the ZFX gene we have observed 11 haplotypes in 336 chromosomes representing a worldwide array of 15 human populations. Two haplotypes representing 77% of all chromosomes were distributed almost evenly among four continents. Five of the remaining haplotypes were detected in Africa and 4 others were restricted to Eurasia and the Americas. Using the information about the ancestral state of the segregating positions (inferred from human-great ape comparisons), we applied coalescent analysis to estimate the age of the polymorphisms and the resulting haplotypes. The oldest haplotype, with the ancestral alleles at all the sites, was observed at low frequency only in two groups of African origin. Its estimated age of 740 to 1100 kyr corresponded to the time to the most recent common ancestor. The two most frequent worldwide distributed haplotypes were estimated at 550 to 840 and 260 to 400 kyr, respectively, while the age of the continentally restricted polymorphisms was 120 to 180 kyr and smaller. Comparison of spatial and temporal distribution of the ZFX haplotypes suggests that modern humans diverged from the common ancestral stock in the Middle Paleolithic era. Subsequent range expansion prevented substantial gene flow among continents, separating African groups from populations that colonized Eurasia and the New World.

A widespread effort to document the amount and geographic distribution of genetic variation in our species is motivated by our curiosity about the origins and prehistory of human populations and our interest in the genetic basis of different diseases. These issues are closely related because the genetic bases of human conditions are a function of the present-day structure of human populations, which by itself cannot be understood without knowing the mechanisms shaping the present genetic variation. Indeed, a starting point for studies that aim to explain the role of genetic variation in disease risk is a description of the quality, quantity, and organization of genetic variation within and between human populations. With the advent of new techniques of DNA analysis it is possible to investigate this variability directly. The greatest progress has been achieved in the analysis of mitochondrial DNA (see JORDE et al. 1998 for a recent review). Nuclear DNA studies have been focused on microsatellite markers (BOWCOCK et al. 1994 ; DEKA et al. 1995 ; JORDE et al. 1995 , JORDE et al. 1997 ; DI RIENZO et al. 1998 ; KIMMEL et al. 1998 ; REICH and GOLDSTEIN 1998 ) and on insertion polymorphisms resulting from Alu retropositions (BATZER et al. 1994 , BATZER et al. 1996 ; STONEKING et al. 1997 ). Studies of nonrecombinogenic portions of the Y chromosome are just now gaining momentum (see JOBLING and TYLER-SMITH 1995 ; HAMMER and ZEGURA 1996 for the reviews; UNDERHILL et al. 1997 ; POLONI et al. 1997 ; SCOZZARI et al. 1997 ; HAMMER et al. 1998 ; MALASPINA et al. 1998 ).

Genetic systems differ in their capacity to reveal the information pertinent to current structure and likely population history. Due to differences in effective population size, the time depth of the autosomal diversity is expected to be four times greater and the X-chromosomal diversity three times greater than that of the Y-chromosome or mitochondrial genome. How far one can look back in time depends upon mutation rates. The fast mutation rate of microsatellites makes them suitable for tracing recent evolutionary events, while the presence of recurrent parallel mutations makes homoplasies frequent at a longer time range. Slower mutation rate and virtually negligible probability of a recurrent mutation in human and great ape lineages characterize classical protein markers, simple nucleotide changes underlying restriction fragment length polymorphisms (RFLPs) and Alu insertion polymorphisms. Studies of these variants provided abundant data on polymorphisms dispersed throughout the genome. However, they often suffered from ascertainment bias. Polymorphisms were initially characterized in a small number of samples primarily of European origin, which led to genotyping of extant world populations for polymorphisms known from one human group. Neglecting polymorphisms that were endemic in the compared populations resulted in an inadvertent loss of pertinent phylogenetic information. Although the importance of ascertaining the polymorphisms in an unbiased manner in all groups of populations under study is now widely recognized (BOWCOCK et al. 1991 ; MOUNTAIN and CAVALLI-SFORZA 1994 ; ROGERS and JORDE 1996 ), only a few non-Y-chromosome nuclear DNA segments have been systematically surveyed for the presence of sequence polymorphisms in a number of distantly related populations. Variable sites were assessed by direct sequencing in an ~3-kb segment of the ß-globin locus (FULLERTON et al. 1994 ; HARDING et al. 1997 ) and a 10-kb portion of the lipoprotein lipase gene (CLARK et al. 1998 ) in 349 and 142 chromosomes, respectively. In another survey of 13 populations from four continents, we screened ~20 independent chromosomes from each of these populations (total of 250) by single-strand conformational polymorphism (SSCP) analysis for polymorphisms in an 8-kb segment of the dystrophin locus (ZIETKIEWICZ et al. 1997 ). Most of these studies indicate the effective population size to be in the range of ~10,000 individuals, which is an order of magnitude lower than that concluded from diversity data at the major histocompatibility complex (TAKAHATA 1993 ). Mitochondrial and microsatellite data suggest recent demographic expansion; this, however, is not the case of simple nuclear polymorphisms. The highest molecular diversity is consistently seen in sub-Saharan Africans, but it is variably interpreted as being due to the older age of this population, the higher population size in Africa, and/or the greater gene flow on this continent. Because the history of a single locus is not sufficient to make conclusions about the history of populations, more variability data from independent nuclear loci are needed.

Studying DNA diversity on the X chromosome offers an advantage of straightforward and unequivocal determination of haplotypes in hemizygous (male) samples. The X-chromosome-specific zinc finger protein (ZFX) locus became of particular interest because of the virtual absence of DNA diversity in its Y chromosome homologue, the ZFY gene (DORIT et al. 1995 ; HUANG et al. 1998 ). Therefore, we have characterized DNA variation in the last introns of ZFX and ZFY in a worldwide sample of human chromosomes (J. JARUZELSKA, E. ZIETKIEWICZ and D. LABUDA, unpublished results; see also HUANG et al. 1998 ). In the present study, the ZFX polymorphisms were combined into haplotypes. The frequencies of ZFX haplotypes and their geographical stratification in 15 globally distributed populations were analyzed. The age of the underlying mutations and the time to the most recent common ancestor (TMRCA) were estimated from the coalescent haplotype tree. The common distribution of the most frequent haplotypes indicates a diversity that characterized the ancestral population of modern humans. Differential patterns of haplotype variability in Africans and non-Africans reflect more recent events. The data are consistent with the divergence of modern human groups from a common ancestral population ~100 kyr ago, followed by range expansion leading to isolation by distance.

	MATERIALS AND METHODS

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

Human DNA samples (nonnominative, characterized only by their origin) represented 15 human populations from four continents. Europeans included Polish (30 unrelated chromosomes), French-Canadians from the Province of Quebec (23), Italians (19), and French (21); Asians were represented by Siberian Nentsi (25), Japanese (21), and Chinese from mainland China (25); Amerindians by Ojibwa (19), Maya (23), and Karitiana from Brazil (16); Africans by Biaka Pygmies from Central Africa Republic (22), M'Buti Pygmies from Congo (22), Rimaibe (19), and Mossi from Burkina Faso (27), and African-Americans from Michigan (24).

DNA variation within 1089 bp of the last intron of the ZFX gene on Xp21.3 (SHIMMIN et al. 1993 ) in the above sample of 336 worldwide chromosomes was investigated by SSCP/heteroduplex approach, as described earlier (ZIETKIEWICZ et al. 1992A , ZIETKIEWICZ et al. 1997 ). Briefly, the analyzed segment was amplified in three overlapping fragments using the following primer pairs: (i) 5'-CGGCAGACTGGCTAAACAA and 5'-ATGCTTATAACATATTTGAGGG (349 bp, annealing temperature 55°); (ii) 5'-AGGACATGGCTGAAACAT and 5'-GTGACAAAAATTTCCACTG (306 bp, 55°); (iii) 5'-GCTGTAAGTTAACGTAAGT and 5'-CTGTTCCAGTTTCTTTGCG (666 bp, 50°); amplification of the first two fragments was done in the presence of 2% formamide. Prior to electrophoresis in a 6% polyacrylamide gel (acrylamide to bisacrylamide ratio 50:1 and 30:1), the 666-bp PCR product was digested with EcoRI into fragments of 319 and 347 bp. Mutations underlying gel mobility variants were identified by dideoxy sequencing.

The mutation rate in the 1089-bp ZFX segment, estimated from the human-chimpanzee and human-orangutan comparisons, was 2.55 (±0.44) x 10^-5 per generation assuming a generation time of 20 yr. The ancestral state of the segregating sites that were due to nucleotide substitutions was inferred by comparison with the orthologous positions in chimpanzee and orangutan DNA. Given the mutation rate in the order of 10^-9 per nucleotide position per year, there is a very small chance that a recurrent event has taken place since the divergence of these species from human lineage 5 and 12 mya. In other words, the probability is negligible that a human allele identical with the corresponding site in chimpanzee and orangutan results from a back mutation or that three identical mutations occurred independently after the separation of these lineages. The identity by state of a human allele and the corresponding positions in great apes is in practice tantamount to their identity by descent. Thus, the human allele identical by state with the chimpanzee and orangutan orthologues was considered ancestral, whereas the other allele of the same polymorphic site was considered new.

The derivation of haplotypes from the genotypes was straightforward in hemizygous males (76 chromosomes), in homozygous female samples (112 chromosomes), and those heterozygous at a single position (92 chromosomes). In multiple heterozygotes (female samples having distinct alleles at more than one polymorphic site), the genotype data were easily resolved into haplotypes when two polymorphisms were present within a single SSCP-analyzed fragment (21 individuals). The remaining haplotypes (7 chromosome pairs) were inferred assuming most likely combinations of the unequivocally resolved haplotypes and taking into account their distribution among populations.

A number of summary statistics (haplotype diversity G, mean number of pairwise differences , nucleotide diversity , and different estimates of , the scaled mutation parameter) were compared using the Arlequin package v.1.1. (SCHNEIDER et al. 1997 ). Arlequin was also used to perform neutrality tests according to EWENS 1972 , WATTERSON 1978 , TAJIMA 1989 , and CHAKRABORTY 1990 , as well as to assess population genetic structure by analysis of molecular variance (AMOVA; EXCOFFIER et al. 1992 ). For detailed description of the parameters and tests please refer to the Arlequin help manual (http://anthropologie.unige.ch/arlequin/).

The maximum-likelihood estimate of using the full information in the sequence data set (a haplotype tree, including the information about the ancestral state) was found using a computational method proposed by GRIFFITHS and TAVARE 1994 , implemented in the program genetree and available from the mathematical genetics group web page at http://stats.ox.ac.uk/. genetree was also used to find the distribution of the TMRCA and the distribution of the ages of the underlying mutations, conditional on the gene tree, using the methods described in GRIFFITHS and TAVARE 1999 ; alternatively, the average age of the mutations was computed from the frequencies of their new (nonancestral) alleles (KIMURA and OHTA 1973 ; GRIFFITHS and TAVARE 1998 ; ZIETKIEWICZ et al. 1998 ).

Maximum-likelihood estimates of obtained with a constant population size model were compared with those using a model of exponential growth. In the exponential growth model there is a decline of the population size backward in time from a current size N₀ = /3µ (22,000 for the maximum-likelihood estimate of 1.66), such that the population size at time t ago is N_t = N₀e^-ßt. Using the genetree program with the expansion rate parameter ß of 0.1, 0.2, and 0.3 and the generating value of 1.66, we obtained no improvement in the log-likelihood compared with the constant population size model (i.e., no evidence for exponential population size expansion was found). The arbitrary ß values selected were lower than these suggested in the genetree manual (0.5 and above) to keep them close to more realistic range.

	RESULTS

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

Ten single-nucleotide polymorphisms within the 1089-bp sequence of the last ZFX intron were previously ascertained in a sample of 336 worldwide-distributed chromosomes (J. JARUZELSKA, E. ZIETKIEWICZ and D. LABUDA, unpublished results). The ancestral and the new allele at each site were inferred by comparison with the orthologous great ape sequences (see MATERIALS AND METHODS). The new allele frequencies in the world population varied from 0.3 to 99.1% (Table 1). The average heterozygosity per polymorphic site was 8.9%, corresponding to the nucleotide diversity (or the mean number of pairwise differences per nucleotide position) of 0.082%. The latter values ranged from 0.06% in Europeans to 0.10% in Asians, while African and Amerindian groups were characterized by of 0.08% (Table 2). The number of polymorphisms in particular populations ranged from two to four (Figure 1B). One polymorphic site (700) was shared by all the groups; at another (1093), the new allele was fixed in all non-African populations. Three sites (514, 632, and 757) had the new allele present only in non-Africans, three others (203, 400, and 491) in Africans, while at two sites (502 and 716) the new allele had a patchy distribution, presumably resulting from admixture (see DISCUSSION). Thus, the distribution of the polymorphisms distinguished Africans from other continental groups.

View larger version (20K): In this window In a new window Download PPT slide   Figure 1. Maximum parsimony network of ZFX haplotypes. Circle sizes in the tree are proportional to the haplotype frequencies in the world (A), individual populations (B), or continents (C). In A, small numbers next to the lines connecting individual haplotypes indicate the underlying mutation sites. Note that most of the haplotypes can be derived from each other by a single change; the order of two mutations separating haplotypes H1 and H8 is unknown; H5/1 is assumed to result from the recombination between H1 and H5 (broken line).

Nucleotide diversity parameters reported here could be underestimated because of limitations of the SSCP method. On the other hand, use of various gel conditions, analyzing partially overlapping fragments, as well as screening a large sample of chromosomes largely improve the detection rate of DNA polymorphisms. We estimate the efficiency of the SSCP/heteroduplex approach at well above 80%, especially given our long experience in using this technique (e.g., ZIETKIEWICZ et al. 1992A , ZIETKIEWICZ et al. 1992B , ZIETKIEWICZ et al. 1997 ; AKALIN et al. 1994 ). This is a conservative estimate in light of the success rate of up to 98% reported by other authors (HAYASHI 1991 ; RAVNIK-GLAVAC et al. 1994 ; WELSH et al. 1997 ). Moreover, HUANG et al. 1998 sequenced the same ZFX segment in a sample of 29 worldwide chromosomes revealing only the ubiquitous site 700. In this context the probability that a frequent polymorphism was missed in our analysis is very low.

The 10 polymorphic positions described above cosegregated as 11 haplotype variants (Table 1). The phylogenetic relations between these variants are represented by a tree (Figure 1A) that summarizes the underlying history of mutations. The tree is rooted at the ancestral haplotype (H0), in which the alleles at all the polymorphic positions are identical by state with the orthologous positions in nonhuman primates. Two most common haplotypes, H1 and H2 are separated from the ancestral one by one and two mutational steps, respectively, and from each other by a single mutation. The only exception is the single copy of H5/1, which could be regarded as a recombinant of H5 and H1. Given the mutation rate in the order of 10^-9 per nucleotide per year, the alternative routes of generating H5/1 through reverse/recurrent mutation at position 700 or at position 502 are much less likely.

The worldwide distribution of haplotype frequencies (Table 1), with H1 and H2 representing 77% of chromosomes, fits very well the infinite-allele model according to Ewens' sampling theory (EWENS 1972 ). The contribution of these two frequent haplotypes was high in all local populations with the exception of the Chinese, in which H1 was not observed. Except for the rare ancestral haplotype H0 (found at low frequency in Mossi and African-Americans), all others were derived from H1 and H2 (Figure 1A) and represented young variants based on coalescent analysis of both the full data and the segregating sites separately (see below). The geographical distribution of these young haplotypes was not uniform (Figure 1B), resembling that of the underlying polymorphisms that differentiated between Sub-Saharan Africans and populations from other continents. Haplotype H9 was restricted to the Chinese population and haplotypes H7 and H8 to Pygmies. The only haplotype resulting from a recombination, H5/1, was found in a single copy in Mossi. A relatively rare haplotype H6 was present in African-Americans in addition to populations in Europe and Siberia, while a single copy of the otherwise African haplotype H5 was found in Maya. The patchy distribution of these two haplotypes presumably represents recent admixture events (see DISCUSSION).

In the virtual absence of recombinations, the haplotype diversity within the analyzed ZFX segment results solely from mutation. Leaving aside the H5/1 recombinant, the remaining haplotypes represent sequences that can be treated formally as mitochondrial DNA and/or haplotypes in the nonrecombining portion on the Y chromosome. The lower overall rate of evolution and lower diversity of the ZFX intron, as compared to other genomic segments, suggests that selection may have acted on this sequence (J. JARUZELSKA, E. ZIETKIEWICZ and D. LABUDA, unpublished results). However, this effect appears too small to influence the population variability of ZFX, which, according to Ewens-Watterson and Tajima tests (not shown and Table 2, respectively), does not differ from that of a neutral locus.

We used a coalescent model to infer the time scale of the origin and evolution of polymorphic variation within the ZFX segment. The distribution of coalescence times and the ages of mutations depend on the mutation rate : = 2Nµ_g, where N is the number of chromosomes in the population (for the X-linked locus, N corresponds to 1.5N_e, assuming equal contribution of both sexes to N_e), and µ_g is the mutation rate per DNA segment per generation. The N_e estimated from , using µ_g = 2.55 x 10^-5 per DNA segment per generation (J. JARUZELSKA, E. ZIETKIEWICZ and D. LABUDA, unpublished results) is ~12,000 (i.e., ~18,000 chromosomes). The corresponding N_e estimate (of ~20,000 individuals) based upon _S is almost twice as large (Table 2). The maximum-likelihood estimate of conditional on the ZFX gene tree (_ML) is 1.66 (±0.43), and the corresponding N_e of ~22,000 is close to that based on _S. In what follows, we compare the distribution of mutational ages using different estimates of (and hence N_e). In principle it is possible to allow for the effect of variability in when inferring the ages of mutations and coalescence times (cf., TAVARE et al. 1997 ), but here we treat the values of as known.

The mutation tree shown in Figure 2 summarizes the results obtained by applying the approach of GRIFFITHS and TAVARE 1999 to infer mutational ages. This tree is rooted because the ancestral state at each site is known. The distribution of haplotypes among the continental groups is indicated below the tree. Using the maximum-likelihood estimate of (_ML), assuming neutrality, random mating, and a constant population size, the expected height of the tree (TMRCA) is estimated at 2.5, and the standard deviation of the height of the tree is 0.9, time being measured in units of N_e generations (Table 2). The corresponding expected TMRCA value using the estimate of is 3.1 N_e generations. Of note, we found no significant evidence of exponential population expansion (see MATERIALS AND METHODS). Conversion of the age estimates above into years (Figure 2) depends on the effective population size N_e.

View larger version (38K): In this window In a new window Download PPT slide   Figure 2. Gene tree of the ZFX polymorphisms in the sample of 335 worldwide chromosomes. Mutations are indicated on the branches at the distances proportional to their estimated average age (2 x 107 simulation runs), using ML (maximum-likelihood estimate of ) or . The order of mutations at positions 400 and 491 cannot be resolved. The values at the right correspond to the coalescent times in Ne generations or, between parentheses, in years (with a generation time of 20 yr and Ne conditional on ML or , i.e., Ne = 22,000 or 12,000, respectively). The corresponding haplotypes are indicated below the tree; the presence of a specific haplotype in the continental population is indicated by a filled box whereas a hatched box indicates the presumed admixture (see text).

In Figure 2 we report the expected ages of the mutations obtained using different values of . The values shown alongside the tree indicate the expected ages in N_e generation units, obtained using the _ML estimate and the estimate (the corresponding values in kyr, using 20 yr per generation, are between parentheses). The expected TMRCA values of 2.5 and 3.1 N_e generation units correspond to 1100 and 740 kyr, respectively. TMRCA also represents the expected age of the ancestral haplotype H0. Two globally distributed haplotypes H1 and H2, which arose as the result of mutation at sites 1093 and 700, are 550 to 840 kyr and 200 to 400 kyr old, respectively. The estimated ages of the remaining haplotypes, i.e., those with the restricted or patchy geographical distribution, vary from 5–9 to 120–180 kyr (Figure 2). The average age of each polymorphism estimated from the full data is very similar to the average age computed from just the frequencies of its new (nonancestral) alleles (KIMURA and OHTA 1973 ; GRIFFITHS and TAVARE 1998 ; ZIETKIEWICZ et al. 1998 ). Figure 3 compares the average ages of polymorphisms (in N_e generation units) and associated 5 and 95% percentiles obtained from the new allele frequencies with those obtained from the full data using the estimates _ML and of . The expected age of the recombinant H5/1 is the same as that of haplotype H9.

View larger version (22K): In this window In a new window Download PPT slide   Figure 3. Comparison of the different estimates of the average ages of polymorphisms and of the time to the most recent common ancestor (TMRCA). Polymorphisms are ordered according to increasing frequency of the new allele. Age estimates obtained from the new allele frequencies (KIMURA and OHTA 1973 ) are represented by the filled bars and coalescent estimates from the full data conditional on ML and by the empty and hatched bars, respectively. The range between the 5th and 95th percentiles is indicated on the bars. Percentiles for estimates based on the full data were obtained from the genetree program; those for the ages estimated from the new allele frequencies were calculated as in GRIFFITHS and TAVARE 1998 .

Table 2 summarizes the population parameters estimated for the world, continents, and all non-Africans together. The estimates of population parameters based on are generally (except for Asians) lower than these based on S, but this difference is not statistically significant. This holds, whether we consider the whole sample or its subsets as representative of the world population. Thus, these population groups look quantitatively similar in spite of the qualitative differences in distribution of young haplotypes. The respective parameters calculated for each of the 15 populations analyzed (not shown) are within the range indicated by the standard deviation values for the larger, continental groups. It is noteworthy that the coalescent time (Table 2) for the world and for the continental samples is almost identical, suggesting again that these groups spread recently from a common ancestral stock and share most of their earlier genetic history (see also DI RIENZO et al. 1998 ).

To assess the level of population structure, we estimated the F_ST parameter, describing the contribution of variance among populations to total variance. The F_ST values were computed using an AMOVA test (EXCOFFIER et al. 1992 ) based on the pairwise distances between haplotypes. The significance of the components of variance was assessed by comparing observed levels with the distribution of 16,000 values obtained by randomization. F_ST calculated for the world population, considered to be composed of 15 subpopulations, was 0.081, indicating that only 8.1% of the total variance resulted from the differences among populations and ~92% was due to the variance within populations. When the world was divided into four continental populations or only into Africans and non-Africans, F_ST values were 0.067 and 0.048, respectively. The F_ST value for Africans was 0.064 and for the cluster combining the remaining groups, 0.062. The results for Europe, Asia, and the Americas were nonsignificant, indicating the absence of population structure in these groups, which is consistent with the conclusions based on other summary statistics above.

	DISCUSSION

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

The geographical distribution of the DNA diversity in the 1.1-kb ZFX segment is not uniform. Two out of 11 haplotypes (H1 and H2) are ubiquitous and represent 77% of the worldwide sample of 336 chromosomes. Similar prevalence of common haplotypes was also found at other genomic loci. At the 2.67-kb ß-globin segment (HARDING et al. 1997 ), 5 common haplotypes contributed ~77% to the worldwide diversity in a sample of 349 chromosomes, and at the 8-kb-long dys44 segment, 6 frequent haplotypes represented ~76% of the worldwide sample of 973 chromosomes (our unpublished data). The remaining 9 haplotypes at the ZFX locus have patchy or continentally restricted distribution. There is good reason to believe that the occurrence of haplotype H6 in African-Americans is due to European admixture in this population (CHAKRABORTY et al. 1992 ), although this will have to be confirmed by studies of a greater number of local populations and/or chromosomes. Similarly, haplotype H5 in Maya may be the result of an African-American admixture. The observation of SEIELSTAD et al. 1994 , who found evidence of African admixture in the same Mayan samples, reinforces our conclusion; African-American admixture was also observed among other Amerindian populations by A. RUIZ-LINARES (personal communication). We can therefore tentatively assign 5 of the haplotypes (H0, H5, H7, H8, and H5/1) to populations from Sub-Saharan Africa. Haplotypes H3, H4, H6, and H9 appear specific to the group of non-African populations. Quantitatively, the haplotype diversity (G values in Table 2) is similar in Africa and the remaining continents. Summary statistics based on the mean number of pairwise differences () indicate the greatest diversity in Asia. However, similar to what was observed at the ß-globin locus (HARDING et al. 1997 ), the estimates based on the number of segregating sites as well as maximum-likelihood estimates from the whole data set show that the diversity is greatest in Africa (_S and _MLE).

All continentally restricted haplotypes except for the ancestral H0 are estimated to have arisen between 120 and 180 kya and the present, on the basis of a coalescent analysis of the full data, and between 160 kya and the present using the mutation ages estimated from the new allele frequency (Figure 2 and Figure 3). In contrast, variants H1 and H2 estimated at 550 to 840 and 200 to 400 kyr, respectively, are found at high proportions in all populations except from Chinese (Figure 2). The absence of H1 from Chinese, where it is replaced by the much younger haplotype H3, may reflect founder effect associated with peopling of Asia and its subregions (e.g., WANG et al. 1991 ). The presence of H1 in two other Asian populations, the paucity of H3 in Europeans, as well as the tree structure in Figure 1A (indicating that H1 "missing" in Chinese is an intermediary form between H0 and both H2 and H3) preclude the interpretation that the two most common haplotypes in the Chinese reflect ancestral haplotype frequencies. Rather, it is the presence of H1 and H2 and their frequency profile that are similar across the continents (Figure 1C), together with the oldest variant H0 that is still detectable in Africa, which can be regarded as the ancestral configuration.

Considering the continents separately, one obtains the same coalescent time (Table 2) as for the pooled sample representing the world. A similar observation was made with microsatellite data (DI RIENZO et al. 1998 ), which are consistent with the shared roots of the contributing populations throughout the substantial period of their earlier genetic history. Worldwide similarity of the frequency profiles of the most common haplotypes could reflect extensive gene flow between the continents. An AMOVA test performed to provide estimates of variance components and of the F-statistics representing the correlation of haplotypes at different levels of hierarchical groupings revealed low variation among constituting populations or groups of populations as compared to the variation within populations. This is similar to the observation at the dystrophin locus (ZIETKIEWICZ et al. 1997 ) and could again indicate a high degree of gene flow connecting the populations as well as continents. However, if such intercontinental genetic exchange existed during the last 100 to 200 kyr, why would a number of distinct young haplotypes persist in isolation, either in Africa or outside, in Eurasia, and the Americas? Likewise, why would the oldest H0 occur at detectable frequency only among Sub-Saharan Africans? We argue that if substantial genetic exchange occurred, it was at the time of ancestral population divergence, preceding the expansion of modern groups. Thereafter, gene flow was restricted by geographic isolation and physical distance between populations colonizing new continents.

With reduced genetic drift, population growth could be conducive to maintaining the frequency profile of polymorphisms characteristic of the ancestral stock from which expanding populations diverged. The excess of the number of segregating sites over the nucleotide diversity (i.e., the negative value of Tajima's D parameter) and the excess of the observed over expected number of haplotypes, E(k), according to EWENS 1972 sampling (see World, Table 2) may suggest that a weak trace of the recent demographic growth is reflected in our sample. However, selection could also be evoked to explain these observations (J. JARUZELSKA, E. ZIETKIEWICZ and D. LABUDA, unpublished results). Moreover, other non-Y-chromosome nuclear loci studied in a similar way provide no evidence for population expansion. HEY 1997 indicated that most nuclear genes show slightly positive skew in frequency spectrum, and studies of ß-globin, dystrophin, and LPL loci (CLARK et al. 1998 ; HARDING et al. 1997 ; ZIETKIEWICZ et al. 1997 , ZIETKIEWICZ et al. 1998 ) reported positive D associated with the non-Y-chromosome nuclear diversity. On the other hand, these observations do not allow firm conclusions. The positive D values reported for ß-globin, dystrophin, and LPL, as well as the negative ones observed here for ZFX were not significantly different from zero. This is further illustrated by the fact that the D values calculated here for individual populations were both positive and negative (not shown). The coalescent estimates from ß-globin data (HARDING et al. 1997 ) did not perform better under a model of expanding population than under a constant population size model, which was also the case for ZFX (see MATERIALS AND METHODS). Thus, if there was any population size expansion, it was not detected by the exponential model or it was probably too recent to be detectable in the observed patterns of ZFX and/or ß-globin diversity. On the other hand, a model that assumes expansion from the Middle/Upper Paleolithic up to the present is suggested by paleontological studies (LAHR and FOLEY 1998 ). A recent population growth is also consistent with genetic studies of mitochondrial DNA or nuclear microsatellites (DI RIENZO and WILSON 1991 ; SLATKIN and HUDSON 1991 ; ROGERS and HARPENDING 1992 ; JORDE et al. 1997 ; DI RIENZO et al. 1998 ; KIMMEL et al. 1998 ). These genetic systems, because of the higher mutation rates, are more suitable for detecting expansion at relatively short evolutionary time scales in the range of 100 kyr. In conclusion, most nuclear loci characterized for simple nucleotide polymorphisms imply a different past in comparison to mitochondrial or microsatellite studies (JORDE et al. 1995 ; HEY 1997 ).

Interestingly, the answer to these discrepancies may lie in the analysis of extended haplotypes, where nuclear diversity generated by rare mutational events would be enhanced by recombinational events, giving insight into more recent population histories, and corresponding to the time frame of mitochondrial and microsatellite mutations. It should be emphasized that the virtual absence of recombinants in ZFX does not necessarily indicate lack of recombination in this region. Assuming the genomic average of 44 crossovers per female meiosis per generation (BROMAN et al. 1998 ) and taking into account that only two-thirds of the X chromosome population participates in female meioses, an average genomic crossover rate is 10^-8 per base pair per generation, i.e., ~10^-5 per 1089 bp in the ZFX segment analyzed. Yet, informative chromosomal pairs that would lead to detectable recombinants represent only 3% of all ZFX haplotype combinations (in other words, 97% of recombinations between "noninformative" ZFX segments remain undetected). Thus, an apparent recombination rate in ZFX is 3 x 10^-7, a hundred times less than the mutation rate of 2.6 x 10^-5 in the same segment.

The ZFX data truncation (a single chromosome representing the recombined haplotype 5/1 was excluded) introduced for the sake of the coalescent analysis was minimal. In some loci, however, the contribution of crossovers to haplotype diversity may be non-negligible. Such an effect is very pronounced in the recently studied LPL locus (CLARK et al. 1998 ; NICKERSON et al. 1998 ), as well as in the dys44 locus (our unpublished data). In this context the good fit between the mutational age estimates based on the full ZFX data and those using just the frequencies of the new allele (Figure 3) is noteworthy. Similar concordance was shown for the ß-globin data of HARDING et al. 1997 (see Figure 4 in ZIETKIEWICZ et al. 1998 ). Although the latter age estimate has a larger variance, it can be calculated for each mutation in the haplotype without knowing the underlying tree.

The genetic diversity measures, pointing to the population size of ~10,000 to 20,000, presumably reflect the size of the ancestral population or speciation bottleneck rather than the harmonic mean of the effective sizes of the world population from Middle Pleistocene to the present. Speciation bottleneck and/or divergence and range expansion from a relatively small ancestral population would have a greater impact on the number of segregating sites than on the nucleotide diversity. If such an event had taken place relatively recently, it would have left no time (depending also upon the size of the expanding populations) for the system to fully recover and reach the equilibrium. As a result, we may observe either the shortage of segregating sites, as in non-African populations at the dys44 locus (ZIETKIEWICZ et al. 1998 ) or a number of new sites that contribute to _S but very little to nucleotide diversity, as in Africans at dys44 (ZIETKIEWICZ et al. 1998 ) and in this study. The net effect on the world diversity reflected by and S may be that the data fit well and can be thus analyzed in terms of a constant population size model.

Our interpretation of the ZFX data is supported by similar results obtained for the ß-globin (HARDING et al. 1997 ) and dys44 (ZIETKIEWICZ et al. 1998 ) loci. Taken together with mitochondrial, microsatellite, and paleontological evidence, these data suggest the divergence of modern human groups from a common ancestral population ~100 kyr ago (taking into account that population trees are usually shorter than these of the underlying genes) followed by range expansion leading to the isolation by distance. What is the geographic origin of the ancestral population? Greater genetic diversity in Africans has frequently been taken as evidence of this population's antiquity (NEI et al. 1975 ; BATZER et al. 1994 ; HORAI et al. 1995 ; TISHKOFF et al. 1996 ). The fact that genetic diversity outside of Africa can be interpreted as being a subset of the African set has been seen as even more compelling evidence for the out-of-Africa origin (STONEKING et al. 1997 ; JORDE et al. 1998 ). However, in the dys44 study, the excess of the African diversity was shown to have accumulated only following the divergence of modern humans (ZIETKIEWICZ et al. 1998 ). Greater diversity in one group may be due to its better representation in the worldwide sample analyzed (BOWCOCK et al. 1991 ; MOUNTAIN and CAVALLI-SFORZA 1994 ; ROGERS and JORDE 1996 ) or to the way we group or define the populations. It may also reflect greater long-term effective population size following expansion (RELETHFORD 1995 ; RELETHFORD and HARPENDING 1995 ; ROGERS and JORDE 1995 ) or particular history of founder effects, extinctions, and expansions of local populations examined (TAKAHATA 1994 ).

Assuming that the ancestral population diverged and underwent range expansion in the Middle/Upper Paleolithic, old haplotypes could have been carried away from the place of origin and preferentially preserved in populations that colonized new areas conducive to relatively uninterrupted growth and demographic prosperity. Another model, which can be described as a single-species model (see FOLEY 1998 ), favors extensive gene flow among distant hominid populations over the whole Paleolithic and thus the contribution of dispersed ancient populations to modern human origins (reviewed in TEMPLETON 1997 ). The argument becomes circular if assuming gene flow we subsequently argue that the presence of a haplotype in a geographically limited sample points to its ancestral place of origin there. Thus, a demonstration of an archaic lineage on one or another continent does not prove anything under any of the above scenarios. With migration and range expansion playing an important role in all the underlying models, the geographic location of the modern human origin cannot be inferred from the genetic evidence alone. Paleontological findings and paleoclimatic data (TATTERSALL 1995 ; LAHR and FOLEY 1998 ) place the origin of the ancestral population in Africa or the Middle East, between the Klasies River and Qafzeh. Our data neither contradict nor prove this scenario. We believe, however, that integrating data from a variety of genetic systems and interpreting them in the context of the outside evidence will be more illuminating in this respect. The need for more extensive sampling of human populations from different geographic areas and of different loci from the nuclear genome appears more than evident, and this study is one of the contributions to fulfill this task.

	FOOTNOTES

¹ These authors have contributed equally to this work.

	ACKNOWLEDGMENTS

Thanks to Dominik Gehl for his assistance in computing and to Raffaela Ballarano for typing the manuscript. This work was supported by National Science Foundation grants DMS 95-04393 to S.T. and SBR-9610147 to M.B., by the Canadian Genome Analysis and Technology Program, and by the Medical Research Council of Canada grant to D.L.

Manuscript received November 19, 1998; Accepted for publication March 19, 1999.

	LITERATURE CITED

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

AKALIN, N., E. ZIETKIEWICZ, W. MAKALOWSKI, and D. LABUDA, 1994  Are CpG sites mutation hot spots in the dystrophin gene? Hum. Mol. Genet. 3:1425-1426[Free Full Text].

BATZER, M. A., M. STONEKING, M. ALEGRIA-HARTMAN, H. BAZAN, and D. H. KASS et al., 1994  African origin of human-specific polymorphic Alu insertions. Proc. Natl. Acad. Sci. USA 91:12288-12292[Abstract/Free Full Text].

BATZER, M. A., S. S. ARCOT, J. W. PHINNEY, M. ALEGRIA-HARTMAN, and D. H. KASS et al., 1996  Genetic variation of recent Alu insertions in human populations. J. Mol. Evol. 42:22-29[Medline].

BOWCOCK, A. M., J. KIDD, J. L. MOUNTAIN, J. M. HEBERT, and L. CAROTENUTO et al., 1991  Drift, admixture, and selection in human evolution: a study with DNA polymorphisms. Proc. Natl. Acad. Sci. USA 88:839-843[Abstract/Free Full Text].

BOWCOCK, A. M., A. RUIZ-LINARES, J. TOMFOHRDE, E. MINCH, and J. R. KIDD et al., 1994  High resolution of human evolutionary trees with polymorphic microsatellites. Nature 368:455-457[Medline].

BROMAN, K. W., J. C. MURRAY, V. C. SHEFFIELD, R. L. WHITE, and J. L. WEBER, 1998  Comprehensive human genetic maps: individual and sex-specific variation in recombination. Am. J. Hum. Genet. 63:861-869[Medline].

CHAKRABORTY, R., 1990  Mitochondrial DNA polymorphism reveals hidden heterogeneity within some Asian populations. Am. J. Hum. Genet. 47:87-94[Medline].

CHAKRABORTY, R., M. I. KAMBOH, M. NWANKWO, and R. E. FERRELL, 1992  Caucasian genes in American Blacks: new data. Am. J. Hum. Genet. 50:145-155[Medline].

CLARK, A., K. M. WEISS, D. A. NICKERSON, S. L. TAYLOR, and A. BUCHANAN et al., 1998  Haplotype structure and population genetic inferences from nucleotide-sequence variation in human lipoprotein lipase. Am. J. Hum. Genet. 63:595-612[Medline].

DEKA, R., L. JIN, M. D. SHRIVER, L. M. YU, and S. DECROO et al., 1995  Population genetics of dinucleotide (dC-dA)_n ·(dG-dT)_n polymorphisms in world populations. Am. J. Hum. Genet. 56:461-474[Medline].

DI RIENZO, A. and A. C. WILSON, 1991  Branching pattern in the evolutionary tree for human mitochondrial DNA. Proc. Natl. Acad. Sci. USA 88:1597-1601[Abstract/Free Full Text].

DI RIENZO, A., P. DONNELLY, C. TOOMAJIAN, B. SISK, and A. HILL et al., 1998  Heterogeneity of microsatellite mutations within and between loci, and implications for human demographic histories. Genetics 148:1269-1284[Abstract/Free Full Text].

DORIT, R. L., H. AKASHI, and W. GILBERT, 1995  Absence of polymorphism at the ZFY locus on the human Y chromosome. Science 268:1183-1185[Abstract/Free Full Text].

EWENS, W. J., 1972  The sampling theory of selectively neutral alleles. Theor. Popul. Biol. 3:87-112[Medline].

EXCOFFIER, L., P. SMOUSE, and J. QUATTRO, 1992  Analysis of molecular variance inferred from metric distance among DNA haplotypes: application to human mitochondria DNA restriction data. Genetics 131:479-491[Abstract].

FOLEY, R., 1998  Genes, evolution and diversity: yet another look at the problem of modern human origins. Evol. Anthropol. 6:191-193.

FULLERTON, S. M., R. M. HARDING, A. J. BOYCE, and J. B. CLEGG, 1994  Molecular and population genetic analysis of allelic sequence diversity at the human ß-globin locus. Proc. Natl. Acad. Sci. USA 91:1805-1809[Abstract/Free Full Text].

GRIFFITHS, R. C. and S. TAVARÉ, 1994  Ancestral inference in population genetics. Stat. Sci. 9:307-319.

GRIFFITHS, R. C. and S. TAVARÉ, 1998  The age of a mutation in a general coalescent tree. Stochastic Models 14:273-295.

GRIFFITHS, R. C. and S. TAVARÉ, 1999  The ages of mutations in gene trees. Ann. Appl. Prob. in press.

HAMMER, M. F. and S. L. ZEGURA, 1996  The role of the Y chromosome in human evolutionary studies. Evol. Anthropol. 5:116-134.

HAMMER, M. F., T. KARAFET, A. RASANAYAGAM, E. T. WOOD, and T. K. ALTHEIDE et al., 1998  Out of Africa and back again: nested cladistic analysis of human Y chromosome variation. Mol. Biol. Evol. 15(4):427-441[Abstract].

HARDING, R. M., S. M. FULLERTON, R. C. GRIFFITHS, J. BOND, and M. J. COX et al., 1997  Archaic African and Asian lineages in the genetic ancestry of modern humans. Am. J. Hum. Genet. 60:772-789[Medline].

HAYASHI, K., 1991  PCR-SSCP: A simple and sensitive method for detection of mutations in the genomic DNA. PCR Methods Appl. 1:34-38[Medline].

HEY, J., 1997  Mitochondrial and nuclear genes present conflicting portraits of human origins. Mol. Biol. Evol. 14:166-172[Abstract].

HORAI, S., K. HAYASAKA, R. KONDO, K. TSUGANE, and N. TAKAHATA, 1995  Recent African origin of modern humans revealed by complete sequences of hominoid mitochondrial DNAs. Proc. Natl. Acad. Sci. USA 92:532-536[Abstract/Free Full Text].

HUANG, W., Y.-X. FU, B. H. CHANG, X. GU, and L. B. JORDE et al., 1998  Sequence variation in ZFX introns in human populations. Mol. Biol. Evol. 15(2):138-142[Abstract].

JOBLING, M. A. and C. TYLER-SMITH, 1995  Fathers and sons: the Y chromosome and human evolution. Trends. Genet. 11:449-456[Medline].

JORDE, L. B., M. J. BAMSHAD, W. S. WATKINS, R. ZENGER, and A. E. FRALEY et al., 1995  Origins and affinities of modern humans: a comparison of mitochondrial and nuclear genetic data. Am. J. Hum. Genet. 57:523-538[Medline].

JORDE, L. B., A. R. ROGERS, M. BAMSHAD, W. S. WATKINS, and P. KRAKOWIAK et al., 1997  Microsatellite diversity and the demographic history of modern humans. Proc. Natl. Acad. Sci. USA 94:3100-3103[Abstract/Free Full Text].

JORDE, L. B., M. BAMSHAD, and A. R. ROGERS, 1998  Using mitochondrial and nuclear DNA markers to reconstruct human evolution. Bioessays 20:126-136[Medline].

KIMMEL, M., R. CHAKRABORTY, J. P. KING, M. BAMSHAD, and W. S. WATKINS et al., 1998  Signatures of population expansion in microsatellite repeat data. Genetics 148:1921-1930[Abstract/Free Full Text].

KIMURA, M. and T. OHTA, 1973  The age of a neutral mutant persisting in a finite population. Genetics 75:199-212[Abstract/Free Full Text].

LAHR, M. M. and R. FOLEY, 1998  Towards a theory of modern human origins: geography, demography and diversity in recent human evolution. Am. J. Phys. Anthropol. 27(Suppl.):137-176.

MALASPINA, P., F. CRUCIANI, B. M. CIMINELLI, L. TERRANATO, and P. SANTOLAMAZZA et al., 1998  Network analyses of Y-chromosomal types in Europe, Northern Africa, and Western Asia reveal specific patterns of geographic distribution. Am. J. Hum. Genet. 63:847-860[Medline].

MOUNTAIN, J. L. and L. L. CAVALLI-SFORZA, 1994  Inference of human evolution through cladistic analysis of nuclear DNA restriction polymorphisms. Proc. Natl. Acad. Sci. USA 91:6515-6519[Abstract/Free Full Text].

NEI, M., T. MARUYAMA, and R. CHAKRABORTY, 1975  The bottleneck effect and genetic variability in populations. Evolution 29:1-10.

NICKERSON, D. A., S. L. TAYLOR, K. M. WEISS, A. G. CLARK, and T. G. HUTCHINSON et al., 1998  Genome resequencing and variation analysis in a 9.7 kb region of the human lipoprotein lipase gene. Nat. Genet. 19:233-240[Medline].

POLONI, E. S., O. SEMINO, G. PASSARINO, A. S. SANTACHIARA-BENERECETTI, and I. DUPANLOUP et al., 1997  Human genetic affinities for Y-chromosome P49a,f/TaqI haplotypes show strong correspondence with linguistics. Am. J. Hum. Genet. 61:1015-1035[Medline].

RAVNIK-GLAVAC, M., D. GLAVAC, and M. DEAN, 1994  Sensitivity of single-strand conformation polymorphism and heteroduplex method for mutation detection in the cystic fibrosis gene. Hum. Mol. Genet. 3:801-807[Abstract/Free Full Text].

REICH, D. E. and D. B. GOLDSTEIN, 1998  Genetic evidence for a Paleolithic human population expansion in Africa. Proc. Natl. Acad. Sci. USA 95:8119-8123[Abstract/Free Full Text].

RELETHFORD, J. H., 1995  Genetics and modern human origins. Evol. Anthropol. 4:53-63.

RELETHFORD, J. H. and H. C. HARPENDING, 1995  Ancient differences in population size can mimic a recent African origin of modern humans. Curr. Anthropol. 36:667-673.

ROGERS, A. R. and H. HARPENDING, 1992  Population growth makes waves in the distribution of pairwise genetic differences. Mol. Biol. Evol. 9:552-569[Abstract].

ROGERS, A. R. and L. B. JORDE, 1995  Genetic evidence on modern human origins. Hum. Biol. 67:1-36[Medline].

ROGERS, A. R. and L. B. JORDE, 1996  Ascertainment bias in estimates of average heterozygosity. Am. J. Hum. Genet. 58:1033-1041[Medline].

SCHNEIDER, S., J.-M. KUEFFER, D. ROESSLI and L. EXCOFFIER, 1997 Arlequin ver. 1.1: a software for population genetic data analysis. Genetics and Biometry Laboratory, University of Geneva, Switzerland.

SCOZZARI, R., F. CRUCIANI, P. MALASPINA, P. SANTOLAMAZZA, and B. M. CIMINELLI et al., 1997  Differential structuring of human populations for homologous X and Y microsatellite loci. Am. J. Hum. Genet. 61:719-733[Medline].

SEIELSTAD, M. T., J. M. HEBERT, A. A. LIN, P. A. UNDERHILL, and M. IBRAHIM et al., 1994  Construction of human Y-chromosomal haplotypes using a new polymorphic A to G transition. Hum. Mol. Genet. 3(12):2159-2161[Abstract/Free Full Text].

SHIMMIN, L. C., B. H. CHANG, and W. H. LI, 1993  Male-driven evolution of DNA sequences. Nature 362:745-777[Medline].

SLATKIN, M. and R. R. HUDSON, 1991  Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics 129:555-562[Abstract].

STONEKING, M., J. J. FONTIUS, S. L. CLIFFORD, H. SOODYALL, and S. S. ARCOT et al., 1997  Alu insertion polymorphisms and human evolution: evidence for a larger population size in Africa. Genome Res. 7:1061-1071[Abstract/Free Full Text].

TAJIMA, F., 1989  Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123:585-595[Abstract/Free Full Text].

TAKAHATA, N., 1993  Allelic genealogy and human evolution. Mol. Biol. Evol. 10:2-22[Abstract].

TAKAHATA, N., 1994  Repeated failures that led to the eventual success in human evolution. Mol. Biol. Evol. 11:803-805[Medline].

TATTERSALL, I., 1995 The Fossil Trail: How We Know What We Think We Know About Human Evolution. American Museum of Natural History. Oxford University Press, New York.

TAVARÉ, S., D. J. BALDING, R. C. GRIFFITHS, and P. DONNELLY, 1997  Inferring coalescence times for molecular sequence data. Genetics 145:505-518[Abstract].

TEMPLETON, A. R., 1997  Out of Africa? What do genes tell us? Curr. Opin. Genet. Dev. 7:841-847[Medline].

TISHKOFF, S. A., E. DIETZSCH, W. SPEED, A. J. PAKSTIS, and J. R. KIDD et al., 1996  Global patterns of linkage disequilibrium at the CD4 locus and modern human origins. Science 271:1380-1387[Abstract].

UNDERHILL, P. A., L. JIN, A. A. LIN, S. Q. MEHDI, and T. JENKINS et al., 1997  Detection of numerous Y chromosome biallelic polymorphisms by denaturing high-performance liquid chromatography. Genome Res. 7:996-1005[Abstract/Free Full Text].

WANG, T., Y. OKANO, R. C. EISENSMITH, M. L. HARVEY, and W. H. Y. LO et al., 1991  Founder effect of a prevalent phenylketonuria mutation in the Oriental population. Proc. Natl. Acad. Sci. USA 88:2146-2150[Abstract/Free Full Text].

WATTERSON, G. A., 1978  The homozygosity test of neutrality. Genetics 88:405-417[Abstract/Free Full Text].

WELSH, J. A., K. CASTREN, and K. H. VAHAKANGAS, 1997  Single-strand conformation polymorphism analysis to detect p53 mutations: characterisation and development of controls. Clin. Chem. 43:2251-2255[Abstract/Free Full Text].

ZIETKIEWICZ, E., D. SINNETT, C. RICHER, G. MITCHELL, and M. VANASSE et al., 1992a  Single-strand conformational polymorphisms (SSCP): detection of useful polymorphisms at the dystrophin locus. Hum. Genet. 89:453-456[Medline].

ZIETKIEWICZ, E., L. R. SIMARD, S. B. MELANÇON, M. VANASSE, and D. LABUDA, 1992b  Carrier status diagnosis in Duchenne muscular dystrophy with "conformational" DNA polymorphism. Lancet 339:134[Medline].

ZIETKIEWICZ, E., V. YOTOVA, M. JARNIK, M. KORAB-LASKOWSKA, and K. K. KIDD et al., 1997  Nuclear DNA diversity in worldwide distributed human populations. Gene 205:161-171[Medline].

ZIETKIEWICZ, E., V. YOTOVA, M. JARNIK, M. KORAB-LASKOWSKA, and K. K. KIDD et al., 1998  Genetic structure of the ancestral population of modern humans. J. Mol. Evol. 47:146-155[Medline].

This article has been cited by other articles:

				K.M. Kusz, L. Tomczyk, M. Sajek, A. Spik, A. Latos-Bielenska, P. Jedrzejczak, L. Pawelczyk, and J. Jaruzelska The highly conserved NANOS2 protein: testis-specific expression and significance for the human male reproduction Mol. Hum. Reprod., March 1, 2009; 15(3): 165 - 171. [Abstract] [Full Text] [PDF]

				A. R. Rogers, S. Wooding, C. D. Huff, M. A. Batzer, and L. B. Jorde Ancestral Alleles and Population Origins: Inferences Depend on Mutation Rate Mol. Biol. Evol., April 1, 2007; 24(4): 990 - 997. [Abstract] [Full Text] [PDF]

				M. K. Shimada, K. Panchapakesan, S. A. Tishkoff, A. Q. Nato Jr, and J. Hey Divergent Haplotypes and Human History as Revealed in a Worldwide Survey of X-Linked DNA Sequence Variation Mol. Biol. Evol., March 1, 2007; 24(3): 687 - 698. [Abstract] [Full Text] [PDF]

				M. F. Hammer, D. Garrigan, E. Wood, J. A. Wilder, Z. Mobasher, A. Bigham, J. G. Krenz, and M. W. Nachman Heterogeneous Patterns of Variation Among Multiple Human X-Linked Loci: The Possible Role of Diversity-Reducing Selection in Non-Africans Genetics, August 1, 2004; 167(4): 1841 - 1853. [Abstract] [Full Text] [PDF]

				M. Le Jossec, T. Wambach, D. Labuda, D. Sinnett, and E. Levy Genetic Diversity Patterns in the SR-BI/II Locus Can Be Explained by a Recent Selective Sweep Mol. Biol. Evol., April 1, 2004; 21(4): 760 - 769. [Abstract] [Full Text] [PDF]

				G. Zhou, Y. Zhai, X. Dong, X. Zhang, F. He, K. Zhou, Y. Zhu, H. Wei, Z. Yao, S. Zhong, et al. Haplotype Structure and Evidence for Positive Selection at the Human IL13 Locus Mol. Biol. Evol., January 1, 2004; 21(1): 29 - 36. [Abstract] [Full Text] [PDF]

				M. A. Saunders, M. F. Hammer, and M. W. Nachman Nucleotide Variability at G6pd and the Signature of Malarial Selection in Humans Genetics, December 1, 2002; 162(4): 1849 - 1861. [Abstract] [Full Text] [PDF]

				N. Yu, F.-C. Chen, S. Ota, L. B. Jorde, P. Pamilo, L. Patthy, M. Ramsay, T. Jenkins, S.-K. Shyue, and W.-H. Li Larger Genetic Differences Within Africans Than Between Africans and Eurasians Genetics, May 1, 2002; 161(1): 269 - 274. [Abstract] [Full Text] [PDF]

				B. A. Payseur and M. W. Nachman Gene Density and Human Nucleotide Polymorphism Mol. Biol. Evol., March 1, 2002; 19(3): 336 - 340. [Abstract] [Full Text] [PDF]

				L.B. Jorde, W.S. Watkins, and M.J. Bamshad Population genomics: a bridge from evolutionary history to genetic medicine Hum. Mol. Genet., October 1, 2001; 10(20): 2199 - 2207. [Abstract] [Full Text] [PDF]

				S. Alonso and J. A. L. Armour A highly variable segment of human subterminal 16p reveals a history of population growth for modern humans outside Africa PNAS, December 14, 2000; (2000) 11244998. [Abstract] [Full Text]

				D. A. Nickerson, S. L. Taylor, S. M. Fullerton, K. M. Weiss, A. G. Clark, J. H. Stengard, V. Salomaa, E. Boerwinkle, and C. F. Sing Sequence Diversity and Large-Scale Typing of SNPs in the Human Apolipoprotein E Gene Genome Res., October 1, 2000; 10(10): 1532 - 1545. [Abstract] [Full Text]

				D. Labuda, E. Zietkiewicz, and V. Yotova Archaic Lineages in the History of Modern Humans Genetics, October 1, 2000; 156(2): 799 - 808. [Abstract] [Full Text]

				M. W. Nachman and S. L. Crowell Contrasting Evolutionary Histories of Two Introns of the Duchenne Muscular Dystrophy Gene, Dmd, in Humans Genetics, August 1, 2000; 155(4): 1855 - 1864. [Abstract] [Full Text]

				J. D. Wall and M. Przeworski When Did the Human Population Size Start Increasing? Genetics, August 1, 2000; 155(4): 1865 - 1874. [Abstract] [Full Text]

				R. Thomson, J. K. Pritchard, P. Shen, P. J. Oefner, and M. W. Feldman Recent common ancestry of human Y chromosomes: Evidence from DNA sequence data PNAS, June 20, 2000; 97(13): 7360 - 7365. [Abstract] [Full Text] [PDF]

				S. Alonso and J. A. L. Armour A highly variable segment of human subterminal 16p reveals a history of population growth for modern humans outside Africa PNAS, January 30, 2001; 98(3): 864 - 869. [Abstract] [Full Text] [PDF]

				R. Martinez-Arias, F. Calafell, E. Mateu, D. Comas, A. Andres, and J. Bertranpetit Sequence Variability of a Human Pseudogene Genome Res., June 1, 2001; 11(6): 1071 - 1085. [Abstract] [Full Text] [PDF]