Current versus historical population sizes in vertebrate species with high gene flow: a comparison based on mitochondrial DNA lineages and inbreeding theory for neutral mutations

Using inbreeding theory as applied to neutral alleles inherited maternally, we generate expected probability distributions of times to identity by descent for random pairs of mitochondrial genotypes within a population or within an entire species characterized by high gene flow. For comparisons with these expectations, empirical distributions of times to most recent common ancestry were calculated (by conventional mtDNA clock calibrations) from mtDNA haplotype distances observed within each of three vertebrate species--American eels, hardhead catfish, and redwinged blackbirds. These species were chosen for analysis because census population size in each is currently large and because both genetic and life-history data are consistent with the postulate that historical gene flow within these species has been high. The observed molecular distances among mtDNA lineages were two to three orders of magnitude lower than predicted from census sizes of breeding females, suggesting that rate of mtDNA evolution is decelerated in these species and/or that long-term effective population size is vastly smaller than present-day population size. Several considerations point to the latter possibility as most likely. The genetic structure of any species is greatly influenced by historical demography; even for species that are currently abundant, mtDNA gene lineages appear to have been channeled through fairly small numbers of ancestors.      

Population structure in the American oyster as inferred by nuclear gene genealogies

Multiple haplotypes from each of three nuclear loci were isolated and sequenced from geographic populations of the American oyster, Crassostrea virginica. In tests of alternative phylogeographic hypotheses for this species, nuclear gene genealogies constructed for these haplotypes were compared to one another, to a mitochondrial gene tree, and to patterns of allele frequency variation in nuclear restriction site polymorphisms (RFLPs) and allozymes. Oyster populations from the Atlantic versus the Gulf of Mexico are not reciprocally monophyletic in any of the nuclear gene trees, despite considerable genetic variation and despite large allele frequency differences previously reported in several other genetic assays. If these populations were separated vicariantly in the past, either insufficient time has elapsed for neutral lineage sorting to have achieved monophyly at most nuclear loci, or balancing selection may have inhibited lineage extinction, or secondary gene flow may have moved haplotypes between regions. These and other possibilities are examined in light of available genetic evidence, and it is concluded that no simple explanation can account for the great variety of population genetic patterns across loci displayed by American oysters. Regardless of the source of this heterogeneity, this study provides an empirical demonstration that different sequences of DNA within the same organismal pedigree can have quite different phylogeographic histories.      

Genetic Variation of Sclerotinia sclerotiorum from Multiple Crops in the North Central United States

Sclerotinia sclerotiorum is an important pathogen of numerous crops in the North Central region of the United States. The objective of this study was to examine the genetic diversity of 145 isolates of the pathogen from multiple hosts in the region. Mycelial compatibility groups (MCG) and microsatellite haplotypes were determined and analyzed for standard estimates of population genetic diversity and the importance of host and distance for genetic variation was examined. MCG tests indicated there were 49 different MCGs in the population and 52 unique microsatellite haplotypes were identified. There was an association between MCG and haplotype such that isolates belonging to the same MCG either shared identical haplotypes or differed at no more than 2 of the 12 polymorphic loci. For the majority of isolates, there was a one-to-one correspondence between MCG and haplotype. Eleven MCGs shared haplotypes. A single haplotype was found to be prevalent throughout the region. The majority of genetic variation in the isolate collection was found within rather than among host crops, suggesting little genetic divergence of S. sclerotiorum among hosts. There was only weak evidence of isolation by distance. Pairwise population comparisons among isolates from canola, dry bean, soybean and sunflower suggested that gene flow between host-populations is more common for some crops than others. Analysis of linkage disequilibrium in the isolates from the four major crops indicated primarily clonal reproduction, but also evidence of genetic recombination for isolates from canola and sunflower. Accordingly, genetic diversity was highest for populations from canola and sunflower. Distribution of microsatellite haplotypes across the study region strongly suggest that specific haplotypes of S. sclerotiorum are often found on multiple crops, movement of individual haplotypes among crops is common and host identity is not a barrier to gene flow for S. sclerotiorum in the north central United States.     

Global population genetic structure and male-mediated gene flow in the green sea turtle (Chelonia mydas): analysis of microsatellite loci

We assessed the degree of population subdivision among global populations of green sea turtles, Chelonia mydas, using four microsatellite loci. Previously, a single-copy nuclear DNA study indicated significant male-mediated gene flow among populations alternately fixed for different mitochondrial DNA haplotypes and that genetic divergence between populations in the Atlantic and Pacific Oceans was more common than subdivisions among populations within ocean basins. Even so, overall levels of variation at single-copy loci were low and inferences were limited. Here, the markedly more variable microsatellite loci confirm the presence of male-mediated gene flow among populations within ocean basins. This analysis generally confirms the genetic divergence between the Atlantic and Pacific. As with the previous study, phylogenetic analyses of genetic distances based on the microsatellite loci indicate a close genetic relationship among eastern Atlantic and Indian Ocean populations. Unlike the single-copy study, however, the results here cannot be attributed to an artifact of general low variability and likely represent recent or ongoing migration between ocean basins. Sequence analyses of regions flanking the microsatellite repeat reveal considerable amounts of cryptic variation and homoplasy and significantly aid in our understanding of population connectivity. Assessment of the allele frequency distributions indicates that at least some of the loci may not be evolving by the stepwise mutation model.     

中国甜菜夜蛾地理种群的遗传分化与基因流

【目的】为了明确我国甜菜夜蛾Spodoptera exigua地理种群间的遗传分化及基因流,阐明该种害虫在我国的种群历史动态。【方法】本研究对采自我国20个地理种群的529头甜菜夜蛾样品进行线粒体COI基因序列测定与分析,利用DnaSP 5.0和Arlequin 3.11软件分析种群间遗传多样性、遗传分化、基因流水平及分子变异,构建了单倍型系统发育树与单倍型网络图。【结果】在所分析的所有529个序列样本中,共检测出10个单倍型,其中Hap_1为所有种群所共享。总群体遗传多样性较低（Hd=0.257±0.025,Pi=0.0007±0.0001,Kxy=0.323）,群体间遗传分化较小（F_ST=0.211）,基因流水平较高（N_m=1.870）。AMOVA分析表明,甜菜夜蛾遗传变异主要来自种群内,种群间变异水平较低。各种群间遗传分化程度与地理距离无显著相关性（R²=0.005,P>0.05）。各单倍型相互散布在不同种群中,未形成明显系统地理结构。中性检验（Tajima’s D=-2.177, P<0.05; Fu’s F_S=-8.629, P<0.05)与错配分布分析表明,我国甜菜夜蛾种群曾经历种群近期扩张。【结论】研究结果揭示,甜菜夜蛾各种群间的基因交流并未受到地理距离的影响,验证了甜菜夜蛾具有高度的迁飞能力。     

Population analysis of vitamin D receptor polymorphisms and the role of genetic ancestry in an admixed population

The vitamin D receptor (VDR) is an essential protein related to bone metabolism. Some VDR alleles are differentially distributed among ethnic populations and display variable patterns of linkage disequilibrium (LD). In this study, 200 unrelated Brazilians were genotyped using 21 VDR single nucleotide polymorphisms (SNPs) and 28 ancestry informative markers. The patterns of LD and haplotype distribution were compared among Brazilian and the HapMap populations of African (YRI), European (CEU) and Asian (JPT+CHB) origins. Conditional regression and haplotype-specific analysis were performed using estimates of individual genetic ancestry in Brazilians as a quantitative trait. Similar patterns of LD were observed in the 5' and 3' gene regions. However, the frequency distribution of haplotype blocks varied among populations. Conditional regression analysis identified haplotypes associated with European and Amerindian ancestry, but not with the proportion of African ancestry. Individual ancestry estimates were associated with VDR haplotypes. These findings reinforce the need to correct for population stratification when performing genetic association studies in admixed populations.     

AFLPs and Mitochondrial Haplotypes Reveal Local Adaptation to Extreme Thermal Environments in a Freshwater Gastropod

The way environmental variation shapes neutral and adaptive genetic variation in natural populations is a key issue in evolutionary biology. Genome scans allow the identification of the genetic basis of local adaptation without previous knowledge of genetic variation or traits under selection. Candidate loci for divergent adaptation are expected to show higher F_ST than neutral loci influenced solely by random genetic drift, migration and mutation. The comparison of spatial patterns of neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection among populations living in contrasting environments. Using the gastropod Radix balthica as a system, we analyzed 376 AFLP markers and 25 mtDNA COI haplotypes for candidate loci and associations with local adaptation among contrasting thermal environments in Lake Mývatn, a volcanic lake in northern Iceland. We found that 2% of the analysed AFLP markers were under directional selection and 12% of the mitochondrial haplotypes correlated with differing thermal habitats. The genetic networks were concordant for AFLP markers and mitochondrial haplotypes, depicting distinct topologies at neutral and candidate loci. Neutral topologies were characterized by intense gene flow revealed by dense nets with edges connecting contrasting thermal habitats, whereas the connections at candidate loci were mostly restricted to populations within each thermal habitat and the number of edges decreased with temperature. Our results suggest microgeographic adaptation within Lake Mývatn and highlight the utility of genome scans in detecting adaptive divergence.     

Linked vs unlinked markers: multilocus microsatellite haplotype-sharing as a tool to estimate gene flow and introgression

We have explored the use of multilocus microsatellite haplotypes to study introgression from cultivated (Malus domestica) into wild apple (Malus sylvestris), and to study gene flow among remnant populations of M. sylvestris. A haplotype consisted of alleles at microsatellite loci along one chromosome. As destruction of haplotypes through recombination occurs much faster than loss of alleles due to genetic drift, the lifespan of a multilocus haplotype is much shorter than that of the underlying alleles. When different populations share the same haplotype, this may indicate recent gene flow between populations. Similarly, haplotypes shared between two species would be a strong signal for introgression. As the expected lifespan of a haplotype depends on the strength of the linkage, the length [in centiMorgans (cM)] of the haplotype shared contains information on the number of generations passed. This application of shared haplotypes is distinct from using haplotype-sharing to detect association between markers and a certain trait. We inferred haplotypes for four to eight microsatellite loci on Linkage Group 10 of apple from genotype data using the program phase, and then identified those haplotypes shared between populations and species. Compared with a Bayesian analysis of unlinked microsatellite loci using the program structure, haplotype-sharing detected a partially different set of putative hybrids. Cultivated haplotypes present in M. sylvestris were short (< 1.5 cM), indicating that introgression had taken place many generations ago, except for two Belgian plants that contained a haplotype of 47.1 cM, indicating recent introgression. In the estimation of gene flow, F(ST) based on unlinked loci indicated small (0.032-0.058) but statistically significant differentiation between some populations only. However, various M. sylvestris haplotypes were shared in nearly all pairwise comparisons of populations, and their length indicated recent gene flow. Hence, all Dutch populations should be considered as one conservation unit. The added value of using sharing of multilocus microsatellite haplotypes as a source of population genetic information is discussed.     

Can microsatellites be used to infer phylogenies? Evidence from population affinities of the Western Canary Island lizard (Gallotia galloti)

Population phylogeographic studies are generally based solely on mtDNA without corroboration, from an independent segregating unit (i.e., nuclear genes), that the mtDNA gene tree represents the organismal phylogeny. This paper attempts to evaluate the utility of microsatellites for this process by use of the Western Canary Island lacertid (Gallotia galloti) as a model. The geological times of island eruptions are known, and well-supported mtDNA phylogenies exist (corroborated as the organismal phylogeny rather than just a gene tree by nuclear random amplified polymorphic DNAs (RAPDs)). The allelic variation in 12 populations from four islands (representing five haplotype lineages) was investigated in five unlinked microsatellite loci. Analysis of molecular variance showed this data to be highly structured. A series of genetic distances among populations was computed based on both the variance in allele frequency (i.e., F(st) related) and the variance in repeat numbers (i.e., R(st) related). The genetic distances based on the former were more highly correlated with the mtDNA genetic distances than those based on the latter. All trees based on both models supported the primary division shown by mtDNA and RAPDs, which is dated at ca. 2.8 to 5.6 mybp (depending on calibration of the mtDNA clock) and which could, under the evolutionary species concept, be regarded separate species. This was achieved despite theoretical problems posed by the use of few loci, suspected bottlenecks, and large population sizes. The finer details were less consistently represented. Nevertheless, this study demonstrates that even a small number of microsatellites can be useful in corroborating the deeper divisions of a population phylogeny.     

<Emphasis Type="Italic">FMR1</Emphasis> haplotype analyses among Indians: a weak founder effect and other findings

This study on allelic/haplotypic fragile X associations evaluated using STR (DXS548, FRAXAC1, FRAXAC2) and SNP (ATL1) markers flanking the (CGG)(n) locus of FMR1is the first report from the large ethnically complex Indian population. Results have been compared with allele/haplotype distributions reported for other major ethnic groups, including White Caucasians, Africans, and Pacific Asians. Though overall allele frequency distributions at the individual loci are more similar to Western Caucasians compared with others, significant differences are observed in haplotypic associations with the mutated X. The striking findings are: (1) high diversity and heterozygosity of haplotypes among fragile X chromosomes ( n=40) and controls ( n=262), including four haplotypes found exclusively in this study sample; (2) weak association of DXS548-FRAXAC1-FRAXAC2 haplotypes, 2-1-3, 6-3-3+ and 7-4-6+ with the disorder, and absence of White Caucasian fragile X haplotypes 6-4-4 and 6-4-5; (3) weak founder effect for the fragile X expansion mutation in the Indians; (4) lack of a continuum of haplotype-based FMR1 alleles between intermediate (CGG)(n) size ranges and expanded alleles; (5) exclusion of ATL1 as a candidate genetic indicator of FMR1 instability. The high STR-based haplotype diversity observed among fragile X lineages, irrespective of ethnic alliances, strongly suggests the inappropriateness of using STR haplotypes to infer predisposition to instability among ethnically separated fragile X pedigrees and may reiterate the need for identifying newer SNPs from this region to not only determine true founder effects for the fragile X mutation, but also decipher possible mechanisms leading to CGG instability.     

Localization of the familial Mediterranean fever gene (FMF) to a 250-kb interval in non-Ashkenazi Jewish founder haplotypes. The French FMF Consortium.

Chromosome 16p13.3 harbors a gene (MEF) associated with familial Mediterranean fever (FMF), a recessive disease very common in populations of Mediterranean ancestry. In the course of positional cloning of MEF, we genotyped 26 non-Ashkenazi Jewish FMF pedigrees (310 meioses) with 15 microsatellite markers, most of which were recently developed by Généthon. Identification of recombination events in the haplotypes allowed narrowing of the MEF interval to a region between D16S3124 (telomeric) and D16S475 (centromeric). Two markers, D16S3070 and D16S3275, a microsatellite marker isolated from a YAC that also contains D16S3070, showed no recombination with the disease. Linkage disequilibrium and haplotype analysis highlighted the existence of a founder haplotype in our population. The core ancestral alleles were present in 71% of MEF-bearing chromosomes at loci D16S3070 and D16S3275. Furthermore, identification of historical crossing-over events in these pedigrees indicated that MEF is located between these two loci, which are both contained in a 250-kb genomic fragment.     

Contribution of mutation, recombination, and gene conversion to chicken MHC-B haplotype diversity

The Mhc is a highly conserved gene region especially interesting to geneticists because of the rapid evolution of gene families found within it. High levels of Mhc genetic diversity often exist within populations. The chicken Mhc is the focus of considerable interest because of the strong, reproducible infectious disease associations found with particular Mhc-B haplotypes. Sequence data for Mhc-B haplotypes have been lacking thereby hampering efforts to systematically resolve which genes within the Mhc-B region contribute to well-defined Mhc-B-associated disease responses. To better understand the genetic factors that generate and maintain genomic diversity in the Mhc-B region, we determined the complete genomic sequence for 14 Mhc-B haplotypes across a region of 59 kb that encompasses 14 gene loci ranging from BG1 to BF2. We compared the sequences using alignment, phylogenetic, and genome profiling methods. We identified gene structural changes, synonymous and non-synonymous polymorphisms, insertions and deletions, and allelic gene rearrangements or exchanges that contribute to haplotype diversity. Mhc-B haplotype diversity appears to be generated by a number of mutational events. We found evidence that some Mhc-B haplotypes are derived by whole- and partial-allelic gene conversion and homologous reciprocal recombination, in addition to nucleotide mutations. These data provide a framework for further analyses of disease associations found among these 14 haplotypes and additional haplotypes segregating and evolving in wild and domesticated populations of chickens.     

Can Microsatellites Be Used to Infer Phylogenies? Evidence from Population Affinities of the Western Canary Island Lizard (Gallotia galloti)

Population phylogeographic studies are generally based solely on mtDNA without corroboration, from an independent segregating unit (i.e., nuclear genes), that the mtDNA gene tree represents the organismal phylogeny. This paper attempts to evaluate the utility of microsatellites for this process by use of the Western Canary Island lacertid (Gallotia galloti) as a model. The geological times of island eruptions are known, and well-supported mtDNA phylogenies exist (corroborated as the organismal phylogeny rather than just a gene tree by nuclear random amplified polymorphic DNAs (RAPDs)). The allelic variation in 12 populations from four islands (representing five haplotype lineages) was investigated in five unlinked microsatellite loci. Analysis of molecular variance showed this data to be highly structured. A series of genetic distances among populations was computed based on both the variance in allele frequency (i.e., F_st related) and the variance in repeat numbers (i.e., R_st related). The genetic distances based on the former were more highly correlated with the mtDNA genetic distances than those based on the latter. All trees based on both models supported the primary division shown by mtDNA and RAPDs, which is dated at ca. 2.8 to 5.6 mybp (depending on calibration of the mtDNA clock) and which could, under the evolutionary species concept, be regarded separate species. This was achieved despite theoretical problems posed by the use of few loci, suspected bottlenecks, and large population sizes. The finer details were less consistently represented. Nevertheless, this study demonstrates that even a small number of microsatellites can be useful in corroborating the deeper divisions of a population phylogeny.     

Beta-globin haplotype analysis suggests that a major source of Malagasy ancestry is derived from Bantu-speaking Negroids.

The origins of the inhabitants of Madagascar have not been fully resolved. Anthropological studies and preliminary genetic data point to two main sources of ancestry of the Malagasy, namely, Indonesian and African, with additional contributions from India and Arabia. The sickle-cell (beta s) mutation is found in populations of African and Indian origin. The frequency of the beta s-globin gene, derived from 1,425 Malagasy individuals, varies from 0 in some highland populations to .25 in some coastal populations. The beta s mutation is thought to have arisen at least five times, on the basis of the presence of five distinct beta s-associated haplotypes, each found in a separate geographic area. Twenty-five of the 35 Malagasy beta s haplotypes were of the typical "Bantu" type, 1 "Senegal" haplotype was found, and 2 rare or atypical haplotypes were observed; the remaining 7 haplotypes were consistent with the Bantu haplotype. The Bantu beta s mutation is thought to have been introduced into Madagascar by Bantu-speaking immigrants (colonists or slaves) from central or east Africa. The Senegal beta s mutation may have been introduced to the island via Portuguese naval explorers. This study provides the first definitive biological evidence that a major component of Malagasy ancestry is derived from African populations, in particular, Bantu-speaking Negroids. beta A haplotypes are also consistent with the claim for a significant African contribution to Malagasy ancestry but are also suggestive of Asian/Oceanic and Caucasoid admixture within the Malagasy population.     

High-resolution Y chromosome haplotypes of Israeli and Palestinian Arabs reveal geographic substructure and substantial overlap with haplotypes of Jews

High-resolution Y chromosome haplotype analysis was performed in 143 paternally unrelated Israeli and Palestinian Moslem Arabs (I&P Arabs) by screening for 11 binary polymorphisms and six microsatellite loci. Two frequent haplotypes were found among the 83 detected: the modal haplotype of the I&P Arabs (approximately 14%) was spread throughout the region, while its one-step microsatellite neighbor, the modal haplotype of the Galilee sample (approximately 8%), was mainly restricted to the north. Geographic substructuring within the Arabs was observed in the highlands of Samaria and Judea. Y chromosome variation in the I&P Arabs was compared to that of Ashkenazi and Sephardic Jews, and to that of North Welsh individuals. At the haplogroup level, defined by the binary polymorphisms only, the Y chromosome distribution in Arabs and Jews was similar but not identical. At the haplotype level, determined by both binary and microsatellite markers, a more detailed pattern was observed. Single-step microsatellite networks of Arab and Jewish haplotypes revealed a common pool for a large portion of Y chromosomes, suggesting a relatively recent common ancestry. The two modal haplotypes in the I&P Arabs were closely related to the most frequent haplotype of Jews (the Cohen modal haplotype). However, the I&P Arab clade that includes the two Arab modal haplotypes (and makes up 32% of Arab chromosomes) is found at only very low frequency among Jews, reflecting divergence and/or admixture from other populations.     

Haplotype identity between individuals who share a CFTR mutation allele “identical by descent”: demonstration of the usefulness of the haplotype-sharing concept for gene mapping in real populations

Cystic fibrosis (CF) patients with the A455E mutation, in both the French Canadian and the Dutch population, share a common haplotype over distances of up to 25 cM. French Canadian patients with the 621+1G→T mutation share a common haplotype of more than 14 cM. In contrast, haplotypes containing the ΔF508 mutation show haplotype identity over a much shorter genomic distance within and between populations, probably because of the multiple introduction of this most common mutation. Haplotype analysis for specific mutations in CF or in other recessive diseases can be used as a model for studying the occurrence of genetic drift conditional on gene frequencies. Moreover, from our results, it can be inferred that analysis of shared haplotypes is a suitable method for genetic mapping in general. Received: 30 November 1995 / Revised: 11 April 1996     

Study of genetic relationships between wild and domesticated grapevine distributed from Middle East Regions to European countries

Archaeobotanical-archaeological, cultural and historical data indicate that grapevine domestication can be dated back from 6000 to 7000 years ago and that it took place in the Caucasian and Middle East Regions. However, events leading to the domestication of this crop species are still an open issue. In this paper, 6 chloroplast microsatellites have been used to assess genetic similarities among, and within, domesticated and wild grapevine accessions representative of 7 distinct geographical regions from the Middle-East to Western Europe. Results show that 2 out of the 6 analyzed chloroplast loci are polymorphic within the 193 domesticated individuals and the 387 samples of 69 wild populations. Allele variants of the Cp-SSR loci combine in a total of 6 different haplotypes. The data show that the haplotype distribution is not homogeneous: the 6 haplotypes are present in the domesticated varieties, but only 5 (haplotype VI is absent) are observed in wild populations. The analysis of haplotype distribution allows discussion of the relationships between the two grape subspecies. The contribution of the wild grape germplasm to the domesticated gene pool still growing in different geographical regions can be, in cases, made evident, suggesting that beside domestication, gene introgression has also played a role in shaping the current varietal landscape of the European viticulture.

---

     

A test for concordance between the multilocus genealogies of genes and microsatellites in the pathogenic fungus Coccidioides immitis

Uncovering the correct phylogeny of closely related species requires analysis of multiple gene genealogies or, alternatively, genealogies inferred from the multiple alleles found at highly polymorphic loci, such as microsatellites. However, a concern in using microsatellites is that constraints on allele sizes may occur, resulting in homoplasious distributions of alleles, leading to incorrect phylogenies. Seven microsatellites from the pathogenic fungus Coccidioides immitis were sequenced for 20 clinical isolates chosen to represent the known genetic diversity of the pathogen. An organismal phylogeny for C. immitis was inferred from microsatellite-flanking sequence polymorphisms and other restriction fragment length polymorphism-containing loci. Two microsatellite genetic distances were then used to determine phylogenies for C. immitis, and the trees found by these three methods were compared. Congruence between the organismal and microsatellite phylogenies occurred when microsatellite distances were based on simple allele frequency data. However, complex mutation events at some loci made distances based on stepwise mutation models unreliable. Estimates of times of divergence for the two species of C. immitis based on microsatellites were significantly lower than those calculated from flanking sequence, most likely due to constraints on microsatellite allele sizes. Flanking-sequence insertions/deletions significantly decreased the accuracy of genealogical information inferred from microsatellite loci and caused interspecific length homoplasies at one of the seven loci. Our analysis shows that microsatellites are useful phylogenetic markers, although care should be taken to choose loci with appropriate flanking sequences when they are intended for use in evolutionary studies.     

Genomewide distribution of high-frequency, completely mismatching SNP haplotype pairs observed to be common across human populations

Knowledge of human haplotype structure has important implications for strategies of disease-gene mapping and for understanding human evolutionary history. Many attributes of SNPs and haplotypes appear to exhibit highly nonrandom behavior, suggesting past operation of selection or other nonneutral forces. We report the exceptional abundance of a particular haplotype pattern in which two high-frequency haplotypes have different alleles at every SNP site (hence the name "yin yang haplotypes"). Analysis of common haplotypes in 62 random genomic loci and 85 gene coding regions in humans shows that the proportion of the genome spanned by yin yang haplotypes is 75%-85%. Population data of 28 genomic loci in Drosophila melanogaster reveal a similar pattern. The high recurrence (>/=85%) of these haplotype patterns in four distinct human populations suggests that the yin yang haplotypes are likely to predate the African diaspora. The pattern initially appeared to suggest deep population splitting or maintenance of ancient lineages by selection; however, coalescent simulation reveals that the yin yang phenomenon can be explained by strictly neutral evolution in a well-mixed population.