Isolates of Cryptococcus neoformans from infected animals reveal genetic exchange in unisexual, alpha mating type populations

Sexual reproduction and genetic exchange are important for the evolution of fungal pathogens and for producing potentially infective spores. Studies to determine whether sex occurs in the pathogenic yeast Cryptococcus neoformans var. grubii have produced enigmatic results, however: basidiospores are the most likely infective propagules, and clinical isolates are fertile and genetically diverse, consistent with a sexual species, but almost all populations examined consist of a single mating type and have little evidence for genetic recombination. The choice of population is critical when looking for recombination, particularly when significant asexual propagation is likely and when latency may complicate assessing the origin of an isolate. We therefore selected isolates from infected animals living in the region of Sydney, Australia, with the assumption that the relatively short life spans and limited travels of the animal hosts would provide a very defined population. All isolates were mating type α and were of molecular genotype VNI or VNII. A lack of linkage disequilibrium among loci suggested that genetic exchange occurred within both genotype groups. Four diploid VNII isolates that produced filaments and basidium-like structures when cultured in proximity to an a mating type strain were found. Recent studies suggest that compatible α-α unions can occur in C. neoformans var. neoformans populations and in populations of the sibling species Cryptococcus gattii. As a mating type strains of C. neoformans var. grubii have never been found in Australia, or in the VNII molecular type globally, the potential for α-α unions is evidence that α-α unisexual mating maintains sexual recombination and diversity in this pathogen and may produce infectious propagules.     

Analysis of genetic structure in soil populations of Rhizobium leguminosarum recovered from the USA and the UK

Although many studies have shown that animal-associated bacterial species exhibit linkage disequilibrium at chromosomal loci, recent studies indicate that both animal-associated and soil-borne bacterial species can display a nonclonal genetic structure in which alleles at chromosomal loci are in linkage equilibrium. To examine the situation in soil-borne species further, we compared genetic structure in two soil populations of Rhizobium leguminosarum bv. trifolii and two populations of R. leguminosarum bv. viciae from two sites in Oregon, with genetic structure in R. leguminosarum bv. viciae populations recovered from peas grown at a site in Washington, USA, and at a site in Norfolk, UK. A total of 234 chromosomal types (ET) were identified among 682 strains analysed for allelic variation at 13 enzyme-encoding chromosomal loci by multilocus enzyme electrophoresis (MLEE). Chi-square tests for heterogeneity of allele frequencies showed that the populations were not genetically uniform. A comparison of the genetic diversity within combined and individual populations confirmed that the Washington population was the primary cause of genetic differentiation between the populations. Each individual population exhibited linkage disequilibrium, with the magnitude of the disequilibrium being greatest in the Washington population and least in the UK population of R. leguminosarum bv. viciae. Linkage disequilibrium in the UK population was created between two clusters of 9 and 23 ETs, which, individually, were in linkage equilibrium. Strong linkage disequilibrium between the two major clusters of 8 and 12 ETs in the Washington population was caused by the low genetic diversity of the ETs within each cluster relative to the inter-cluster genetic distance. Because neither the magnitude of genetic diversity nor of linkage disequilibrium increased as hierarchical combinations of the six local populations were analysed, we conclude that the populations have not been isolated from each other for sufficient time, nor have they been exposed to enough selective pressure to develop unique multilocus genetic structure.     

Genetic variability and distribution of mating type alleles in field populations of Leptosphaeria maculans from France

Leptosphaeria maculans is the most ubiquitous fungal pathogen of Brassica crops and causes the devastating stem canker disease of oilseed rape worldwide. We used minisatellite markers to determine the genetic structure of L. maculans in four field populations from France. Isolates were collected at three different spatial scales (leaf, 2-m2 field plot, and field) enabling the evaluation of spatial distribution of the mating type alleles and of genetic variability within and among field populations. Within each field population, no gametic disequilibrium between the minisatellite loci was detected and the mating type alleles were present at equal frequencies. Both sexual and asexual reproduction occur in the field, but the genetic structure of these populations is consistent with annual cycles of randomly mating sexual reproduction. All L. maculans field populations had a high level of gene diversity (H = 0.68 to 0.75) and genotypic diversity. Within each field population, the number of genotypes often was very close to the number of isolates. Analysis of molecular variance indicated that >99.5% of the total genetic variability was distributed at a small spatial scale, i.e., within 2-m2 field plots. Population differentiation among the four field populations was low (GST < 0.02), suggesting a high degree of gene exchange between these populations. The high gene flow evidenced here in French populations of L. maculans suggests a rapid countrywide diffusion of novel virulence alleles whenever novel resistance sources are used.     

Genetic diversity and linkage disequilibrium in the Polynesian population of Niue Island

Isolated populations that recently have been derived from small homogeneous groups of founders should have low genetic diversity and high levels of linkage disequilibrium and should be ideal for mapping ancestral polymorphisms that influence complex genetic disease susceptibility. Populations that fulfill these criteria have been difficult to identify. We have been looking for Polynesian populations with these characteristics, because Polynesians have high rates of complex genetic diseases. In Niue Islanders all ancestral female (mitochondrial HSVI sequence) and 90.4% of ancestral male (Y-chromosome haplogroup) lineages are of Southeast Asian origin. The frequency of European Y-chromosome haplogroups is 7.2%. The diversities of mitochondrial HSV1 sequences (h = 0.18 +/- 0.05) and Y-chromosome haplo-groups (h = 0.18 +/- 0.05) are lower than values published for any other population. Ten autosomal microsatellites spaced over 5.8 cM show low allele numbers in Niue Islanders relative to Europeans (55 vs. 88 total alleles, respectively) and a modest reduction in heterozygous loci (0.71 +/- 0.02 vs. 0.78 +/- 0.02, p = 0.04). The higher linkage disequilibrium (d2) between these loci in Niue Islanders relative to Europeans (p = 0.001) is negatively correlated (r = -0.47, p = 0.01) with genetic distance. In summary, Niue Islanders are genetically isolated and have a homogeneous Southeast Asian ancestry. They have reduced autosomal genetic diversity and high levels of linkage disequilibrium that are consistent with the influence of genetic drift mechanisms, such as a founder effect or bottlenecks. High-powered linkage disequilibrium studies designed to map ancestral polymorphisms that influence complex genetic disease susceptibility may be feasible in this population.     

A global perspective on genetic variation at the ADH genes reveals unusual patterns of linkage disequilibrium and diversity

Variants of different Class I alcohol dehydrogenase (ADH) genes have been shown to be associated with an effect that is protective against alcoholism. Previous work from our laboratory has shown that the two sites showing the association are in linkage disequilibrium and has identified the ADH1B Arg47His site as causative, with the ADH1C Ile349Val site showing association only because of the disequilibrium. Here, we describe an initial study of the nature of linkage disequilibrium and genetic variation, in population samples from different regions of the world, in a larger segment of the ADH cluster (including the three Class I ADH genes and ADH7). Linkage disequilibrium across approximately 40 kb of the Class I ADH cluster is moderate to strong in all population samples that we studied. We observed nominally significant pairwise linkage disequilibrium, in some populations, between the ADH7 site and some Class I ADH sites, at moderate values and at a molecular distance as great as 100 kb. Our data indicate (1) that most ADH-alcoholism association studies have failed to consider many sites in the ADH cluster that may harbor etiologically significant alleles and (2) that the relevance of the various ADH sites will be population dependent. Some individual sites in the Class I ADH cluster show Fst values that are among the highest seen among several dozen unlinked sites that were studied in the same subset of populations. The high Fst values can be attributed to the discrepant frequencies of specific alleles in eastern Asia relative to those in other regions of the world. These alleles are part of a single haplotype that exists at high (>65%) frequency only in the eastern-Asian samples. It seems unlikely that this haplotype, which is rare or unobserved in other populations, reached such high frequency because of random genetic drift alone.     

Genetic variability of HLA in the Dariusleut Hutterites. A comparative genetic analysis of the Hutterities, the Amish, and other selected Caucasian populations.

There are three endogamous subdivisions of the Hutterite population, a North American religious isolate. These individuals live on communal farms, and residence is strictly patrilocal. We report on the distributions of HLA-A and B alleles and haplotypes in 203 married women from one subdivision--the Dariusleut--in Alberta, Canada. We demonstrate that there is significant linkage disequilbrium among a large fraction of the distinct haplotypes in the Dariusleut Hutterite data; there is a restriction in the number of distinct haplotypes present in the Dariusleut; the Hutterites and the Old Order Amish (Lancaster County, Pennsylvania) are the most genetically distant pair of populations in an ensemble of 11 Caucasian populations; and, finally, the Old Order Amish and the Hutterites are approximately as distant from the Indiana Amish as they are from the eight other Caucasian populations, which are tightly clustered in the space of gene frequencies. These results are consistent with the fact that the Amish and the Hutterites are genetic isolates with small numbers of founders. Certain haplotypes show significant linkage disequilibrium in these as well as in other Caucasian samples. Thus, some of the linkage disequilibrium antedates the formation of these Anabaptist sects.     

Genetic diversity and bottleneck analysis of Indian bellary sheep by microsatellite markers

Bellary sheep population variability and structure was investigated genetically utilizing FAO recommended microsatellite markers. Genetic variation at 20 microsatellite loci, population structure, and genetic bottleneck hypothesis were examined. Estimates of genetic variability such as effective number of alleles and gene diversities revealed substantial genetic variation frequently displayed by microsatellite markers. A total of 133 alleles were detected. Average polymorphism across the studied loci and expected gene diversity in the population were 1.419 +/- 0.405 and 0.684 +/- 0.140, respectively. No significant genotypic linkage disequilibrium was detected across population, suggesting no evidence of linkage between loci. The population was observed to be significantly differentiated into different groups, showed fairly high level of inbreeding (f = 0.253 +/- 0.050) and global heterozygote deficit. Population structure analysis indicated the intermixing/introduction of unique/rare alleles in these migrating flocks. A normal L-shaped distribution of mode-shift test, non-significant heterozygosity excess on the basis of different models, as revealed from Sign, Standardized differences and Wilcoxon sign rank tests suggested that there was no recent bottleneck. The study revealed that even breed with increasing population trend needs genetic management for the conservation and improvement.     

Mating Patterns of Plasmodium falciparum

Recent empirical data have enabled a more informed debate over the extent of clonality in Plasmodium falciparum populations. Oocyst heterozygosity data reveal that the mating structure of malaria populations varies according to the transmission intensity. This finding provides a more detailed picture of the malaria mating structure than previous conclusions, which were based on indirect measures of population mating structure, ie. linkage disequilibrium analyses. In this article, Ric Paul and Karen Day discuss aspects of the genetic structure of malaria populations as evidenced by oocyst heterozygosity and linkage disequilibrium data. They address the difficulties of performing genetic analyses of malaria parasite population structure inherent in parasite sampling, why two identical parasites are rarely observed in the field and how features of the epidemiology determine parasite population structure.     

Genetic diversity and structure of natural fragmented Chamaecyparis obtusa populations as revealed by microsatellite markers

The genetic diversity and population structure of hinoki (Chamaecyparis obtusa) in Japan were investigated by examining the distribution of alleles at 13 microsatellite loci in 25 natural populations from Iwaki in northern Japan to Yakushima Island in southern Japan. On average, 26.9 alleles per locus were identified across all populations and 4.0% of the genetic variation was retained among populations (G _ST = 0.040). According to linkage disequilibrium analysis, estimates of effective population size and detected evidence of bottleneck events, the genetic diversity of some populations may have declined as a result of fragmentation and/or over-exploitation. The central populations located in the Chubu district appear to have relatively large effective population sizes, while marginal populations, such as the Yakushima, Kobayashi and Iwaki populations, have smaller effective population sizes and are isolated from the other populations. Microsatellite analysis revealed the genetic uniqueness of the Yakushima population. Although genetic differentiation between populations was low, we detected a gradual cline in the genetic structure and found that locus Cos2619 may be non-neutral with respect to natural selection.     

Population genetic structure of clinical and environmental isolates of Blastomyces dermatitidis, based on 27 polymorphic microsatellite markers

Blastomyces dermatitidis, a thermally dimorphic fungus, is the etiologic agent of North American blastomycosis. Clinical presentation is varied, ranging from silent infections to fulminant respiratory disease and dissemination to skin and other sites. Exploration of the population genetic structure of B. dermatitidis would improve our knowledge regarding variation in virulence phenotypes, geographic distribution, and difference in host specificity. The objective of this study was to develop and test a panel of microsatellite markers to delineate the population genetic structure within a group of clinical and environmental isolates of B. dermatitidis. We developed 27 microsatellite markers and genotyped B. dermatitidis isolates from various hosts and environmental sources (n=112). Assembly of a neighbor-joining tree of allele-sharing distance revealed two genetically distinct groups, separated by a deep node. Bayesian admixture analysis showed that two populations were statistically supported. Principal coordinate analysis also reinforced support for two genetic groups, with the primary axis explaining 61.41% of the genetic variability. Group 1 isolates average 1.8 alleles/locus, whereas group 2 isolates are highly polymorphic, averaging 8.2 alleles/locus. In this data set, alleles at three loci are unshared between the two groups and appear diagnostic. The mating type of individual isolates was determined by PCR. Both mating type-specific genes, the HMG and α-box domains, were represented in each of the genetic groups, with slightly more isolates having the HMG allele. One interpretation of this study is that the species currently designated B. dermatitidis includes a cryptic subspecies or perhaps a separate species.     

Genetic structure of a natural population of Dictyostelium discoideum, a cellular slime mould

Dictyostelium discoideum is a eukaryotic microbe feeding on soil bacteria. A first step towards describing the genetic structure of populations of this species was made by examining multiple isolates from a single locale. The isolates were grown clonally and their RFLP patterns compared, using a probe specific for a family of tRNA genes. Thirty-nine types were distinguished in 54 isolates. To determine if genetic exchange occurs among members of the population, an analysis of linkage disequilibrium was performed on the RFLP data. Little disequilibrium was found, implying gene flow in the population. In conflict with this result is the finding that no recombinant progeny were recovered from many attempted crosses between pairs of isolates. The tentative conclusion is that genetic exchange does not in fact occur, and that the observed shuffling of RFLP bands is caused by insertion and excision of transposons known to be associated with the tRNA genes of Dictyostelium.     

Genetic diversity and bottleneck analysis of Indian bellary sheep by microsatellite markers

Bellary sheep population variability and structure was investigated genetically utilizing FAO recommended microsatellite markers. Genetic variation at 20 microsatellite loci, population structure, and genetic bottleneck hypothesis were examined. Estimates of genetic variability such as effective number of alleles and gene diversities revealed substantial genetic variation frequently displayed by microsatellite markers. A total of 133 alleles were detected. Average polymorphism across the studied loci and expected gene diversity in the population were 1.419 ± 0.405 and 0.684 ± 0.140, respectively. No significant genotypic linkage disequilibrium was detected across population, suggesting no evidence of linkage between loci. The population was observed to be significantly differentiated into different groups, showed fairly high level of inbreeding (f = 0.253 ± 0.050) and global heterozygote deficit. Population structure analysis indicated the intermixing/introduction of unique/rare alleles in these migrating flocks. A normal L-shaped distribution of mode-shift test, non-significant heterozygosity excess on the basis of different models, as revealed from sign, standardized differences and Wilcoxon sign rank tests suggested that there was no recent bottleneck. The study revealed that even a breed with increasing population trend needs genetic management for the conservation and improvement. The text was submitted by the authors in English.     

Marker selection for the transmission/disequilibrium test, in recently admixed populations.

Recent admixture between genetically differentiated populations can result in high levels of association between alleles at loci that are <=10 cM apart. The transmission/disequilibrium test (TDT) proposed by Spielman et al. (1993) can be a powerful test of linkage between disease and marker loci in the presence of association and therefore could be a useful test of linkage in admixed populations. The degree of association between alleles at two loci depends on the differences in allele frequencies, at the two loci, in the founding populations; therefore, the choice of marker is important. For a multiallelic marker, one strategy that may improve the power of the TDT is to group marker alleles within a locus, on the basis of information about the founding populations and the admixed population, thereby collapsing the marker into one with fewer alleles. We have examined the consequences of collapsing a microsatellite into a two-allele marker, when two founding populations are assumed for the admixed population, and have found that if there is random mating in the admixed population, then typically there is a collapsing for which the power of the TDT is greater than that for the original microsatellite marker. A method is presented for finding the optimal collapsing that has minimal dependence on the disease and that uses estimates either of marker allele frequencies in the two founding populations or of marker allele frequencies in the current, admixed population and in one of the founding populations. Furthermore, this optimal collapsing is not always the collapsing with the largest difference in allele frequencies in the founding populations. To demonstrate this strategy, we considered a recent data set, published previously, that provides frequency estimates for 30 microsatellites in 13 populations.     

The influence of recombination on the population structure and evolution of the human pathogen Neisseria meningitidis.

The extent to which recombination disrupts the bifurcating treelike phylogeny and clonal structure imposed by binary fission on bacterial populations remains contentious. Here, we address this question with a study of nucleotide sequence data from 107 isolates of the human pathogen Neisseria meningitidis. Gene fragments from 12 house-keeping loci distributed around the meningococcal chromosome were analyzed, showing that (1) identical alleles are disseminated among genetically diverse isolates, with no evidence for linkage disequilibrium; (2) different loci give distinct and incongruent phylogenetic trees; and (3) allele sequences are incompatible with a bifurcating treelike phylogeny at all loci. These observations are consistent with the hypothesis that meningococcal populations comprise organisms assembled from a common gene pool, with alleles and allele fragments spreading independently, together with the occasional importation of genetic material from other species. Further, they support the view that recombination is an important genetic mechanism in the generation new meningococcal clones and alleles. Consequently, for anything other than the short-term evolution of this species, a bifurcating treelike phylogeny is not an appropriate model.     

Length Variation of Cag/Caa Trinucleotide Repeats in Natural Populations of Drosophila Melanogaster and Its Relation to the Recombination Rate

Eleven genes distributed along the Drosophila melanogaster chromosome 2 and showing exonic tandem repeats of glutamine codons (CAG or CAA) were surveyed for length variation in a sample of four European and African populations. Only one gene was monomorphic. Eight genes were polymorphic in all populations, with a total number of alleles varying between five and 12 for 120 chromosomes. The average heterozygozity per locus and population was 0.41. Selective neutrality in length variation could not be rejected under the assumptions of the infinite allele model. Significant population subdivision was found though no geographical pattern emerged, all populations being equally different. Significant linkage disequilibrium was found in four out of seven cases where the genetic distance between loci was <1 cM and was negligible when the distance was larger. There is evidence that these associations were established after the populations separated. An unexpected result was that variation at each locus was independent of the coefficient of exchange, although the latter ranged from zero to the relatively high value of 6.7%. This would indicate that background selection and selective hitchhiking, which are thought to affect levels of nucleotide substitution polymorphism, have no effect on trinucleotide repeat variation.     

The IS1167 Insertion Sequence Is a Phylogenetically Informative Marker Among Isolates of Serotype 6B Streptococcus pneumoniae

The phylogenetic utility of the IS1167 insertion sequence was examined with restriction fragment length polymorphism (RFLP) analyses of a sample of 50, predominantly invasive, capsular serotype 6B Streptococcus pneumoniae isolates previously characterized by multilocus enzyme electrophoresis (MLEE). The strains represented a genetically diverse assemblage of 34 distinct clonotypes composed of 26 restriction fragment types and 23 multilocus enzyme types. All isolates carried the IS1167 insertion sequence, with an average of 9.5 copies. The cross-classification of isolates based on RFLP and MLEE typing schemes was 81% concordant. Phylogenetic analyses demonstrated a significant (P < 0.0001) association between strains of a given RFLP lineage with those of a given MLEE lineage. A significant correlation (P < 0.00004) was also found between the proportion of restriction fragments shared by any given pair of isolates and their genetic distances estimated from the MLEE data. Parity between the two genetic markers implied that the sampled isolates were in linkage disequilibrium. The existence of nonrandom associations among genetic loci was confirmed by Monte Carlo analyses of the MLEE data. These studies, thus, demonstrated that invasive pneumococcal isolates of a single capsule type recovered on a regional scale can retain a largely clonal population structure over a period of 8 years. The ability to detect linkage disequilibrium and generate relatively congruent dendrograms based on distance and parsimony methods suggested that the restriction fragment data were robust to phylogenetic analysis. Received: 20 May 1997 / Accepted: 20 November 1997     

Genetic structure of South Australian Pyrenophora teres populations as revealed by microsatellite analyses

The aim of this study was to determine the genetic structure of South Australian field populations of the barley net blotch pathogens, Pyrenophora teres f. sp. teres (PTT) and P. teres f. sp. maculata (PTM), using microsatellite DNA markers. Three PTT populations (76 isolates total) and two PTM populations (43 isolates total) were sampled from separate fields during a single growing season. The results showed that of the 20 microsatellite loci examined, 17 (85 %) were polymorphic within the PTT and PTM populations. In total, 120 distinct alleles were identified of which only 11 (9 %) were shared between the two population types. Nei’s measure of gene diversity across the PTT and PTM populations was similar at 0.38 and 0.40, respectively, and also much higher than previously reported from studies in which other types of molecular markers were used. The coefficient of genetic differentiation among both populations was the same (G_ST = 0.03) and the low and insignificant estimates of F_ST, as indicated by θ, between populations of the same type (PTT: θ < 0.008, PTM: θ = 0.014) indicated that isolates sampled from different areas within the same field were genetically similar. In contrast, high and significant genetic differentiation was observed among and between populations of different type (G_ST = 0.42, θ > 0.567). The high number of unique multilocus haplotypes observed within the PTT (84 %) and PTM (100 %) populations, combined with a 1:1 distribution of both mating types, suggested that sexual reproduction was predominant among these populations. However, tests for multilocus associations showed that both PTT and PTM populations were in significant linkage disequilibrium. Although the levels of disequilibrium were low, an asexual reproductive component could not be excluded.     

Genetic Variability and Distribution of Mating Type Alleles in Field Populations of Leptosphaeria maculans from France

Leptosphaeria maculans is the most ubiquitous fungal pathogen of Brassica crops and causes the devastating stem canker disease of oilseed rape worldwide. We used minisatellite markers to determine the genetic structure of L. maculans in four field populations from France. Isolates were collected at three different spatial scales (leaf, 2-m² field plot, and field) enabling the evaluation of spatial distribution of the mating type alleles and of genetic variability within and among field populations. Within each field population, no gametic disequilibrium between the minisatellite loci was detected and the mating type alleles were present at equal frequencies. Both sexual and asexual reproduction occur in the field, but the genetic structure of these populations is consistent with annual cycles of randomly mating sexual reproduction. All L. maculans field populations had a high level of gene diversity (H = 0.68 to 0.75) and genotypic diversity. Within each field population, the number of genotypes often was very close to the number of isolates. Analysis of molecular variance indicated that >99.5% of the total genetic variability was distributed at a small spatial scale, i.e., within 2-m² field plots. Population differentiation among the four field populations was low (G_ST < 0.02), suggesting a high degree of gene exchange between these populations. The high gene flow evidenced here in French populations of L. maculans suggests a rapid countrywide diffusion of novel virulence alleles whenever novel resistance sources are used.     

Microsatellite DNA markers for the snow vole (Chionomys nivalis)

A total of 14 dinucleotide microsatellite loci were characterized in the snow vole (Chionomys nivalis). Allelic polymorphism across all loci and 28 individuals representing a single population in the Swiss Alps was high (mean = 10.1 alleles). No significant linkage disequilibrium between pairs of loci and no departure from Hardy–Weinberg equilibrium were found. These loci will be useful for describing mating systems and population structure and to investigate the genetic consequences of a species living in a highly fragmented habitat.     

Nucleotide variability at G6pd and the signature of malarial selection in humans

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common enzymopathy in humans. Deficiency alleles for this X-linked disorder are geographically correlated with historical patterns of malaria, and the most common deficiency allele in Africa (G6PD A-) has been shown to confer some resistance to malaria in both hemizygous males and heterozygous females. We studied DNA sequence variation in 5.1 kb of G6pd from 47 individuals representing a worldwide sample to examine the impact of selection on patterns of human nucleotide diversity and to infer the evolutionary history of the G6PD A- allele. We also sequenced 3.7 kb of a neighboring locus, L1cam, from the same set of individuals to study the effect of selection on patterns of linkage disequilibrium. Despite strong clinical evidence for malarial selection maintaining G6PD deficiency alleles in human populations, the overall level of nucleotide heterozygosity at G6pd is typical of other genes on the X chromosome. However, the signature of selection is evident in the absence of genetic variation among A- alleles from different parts of Africa and in the unusually high levels of linkage disequilibrium over a considerable distance of the X chromosome. In spite of a long-term association between Plasmodium falciparum and the ancestors of modern humans, patterns of nucleotide variability and linkage disequilibrium suggest that the A- allele arose in Africa only within the last 10,000 years and spread due to selection.