A test of genomic modularity among life‐history adaptations promoting speciation with gene flow

Speciation with gene flow may require adaptive divergence of multiple traits to generate strong ecologically based reproductive isolation. Extensive negative pleiotropy or physical linkage of genes in the wrong phase affecting these diverging traits may therefore hinder speciation, while genetic independence or “modularity” among phenotypic traits may reduce constraints and facilitate divergence. Here, we test whether the genetics underlying two components of diapause life history, initial diapause intensity and diapause termination timing, constrain differentiation between sympatric hawthorn and apple‐infesting host races of the fly Rhagoletis pomonella through analysis of 10,256 SNPs measured via genotyping‐by‐sequencing (GBS). Loci genetically associated with diapause termination timing were mainly observed for SNPs mapping to chromosomes 1–3 in the genome, most notably for SNPs displaying higher levels of linkage disequilibrium (LD), likely due to inversions. In contrast, selection on initial diapause intensity affected loci on all five major chromosomes of the genome, specifically those showing low levels of LD. This lack of overlap in genetically associated loci suggests that the two diapause phenotypes are largely modular. On chromosome 2, however, intermediate level LD loci and a subgroup of high LD loci displayed significant negative relationships between initial diapause intensity and diapause termination time. These gene regions on chromosome 2 therefore affected both traits, while most regions were largely independent. Moreover, loci associated with both measured traits also tended to exhibit highly divergent allele frequencies between the host races. Thus, the presence of nonoverlapping genetic modules likely facilitates simultaneous, adaptive divergence for the measured life‐history components.     

物种形成过程中的分化基因组岛及其形成机制

理解物种形成机制是生态和进化领域的重要任务。得益于测序技术的快速发展, 越来越多研究发现分化种群(亚种、物种)间的基因组常呈现异质性分化景观, 存在分化基因组岛, 这被认为是基因流存在下的歧化选择引起的, 支持基因流存在下的成种假说。然而, 基因渐渗、祖先多态性的差异分选、连锁选择等其他进化过程也可导致分化基因组岛的形成。现有实证研究在解析分化基因组岛的形成机制时, 往往忽略了上述其他进化过程的作用。为此, 本文在辨析分化基因组岛相关概念的基础上, 总结了利用种群基因组数据鉴定分化基因组岛的方法, 对比了不同进化过程形成分化基因组岛的特征, 指出在区分不同机制时联用基因渐渗程度、绝对分化指数(d_XY)、相对节点深度(RND)、重组率等多个指标的必要性, 归纳了物种形成过程中分化基因组岛形成机制解析的研究思路, 并对未来在生殖隔离机制上的深入探索以及实证研究的整合分析等方面进行了展望。     

Interpreting the genomic landscape of speciation: a road map for finding barriers to gene flow

Speciation, the evolution of reproductive isolation among populations, is continuous, complex, and involves multiple, interacting barriers. Until it is complete, the effects of this process vary along the genome and can lead to a heterogeneous genomic landscape with peaks and troughs of differentiation and divergence. When gene flow occurs during speciation, barriers restricting gene flow locally in the genome lead to patterns of heterogeneity. However, genomic heterogeneity can also be produced or modified by variation in factors such as background selection and selective sweeps, recombination and mutation rate variation, and heterogeneous gene density. Extracting the effects of gene flow, divergent selection and reproductive isolation from such modifying factors presents a major challenge to speciation genomics. We argue one of the principal aims of the field is to identify the barrier loci involved in limiting gene flow. We first summarize the expected signatures of selection at barrier loci, at the genomic regions linked to them and across the entire genome. We then discuss the modifying factors that complicate the interpretation of the observed genomic landscape. Finally, we end with a road map for future speciation research: a proposal for how to account for these modifying factors and to progress towards understanding the nature of barrier loci. Despite the difficulties of interpreting empirical data, we argue that the availability of promising technical and analytical methods will shed further light on the important roles that gene flow and divergent selection have in shaping the genomic landscape of speciation.     

Low‐coverage genomic data resolve the population divergence and gene flow history of an Australian rain forest fig wasp

Population divergence and gene flow are key processes in evolution and ecology. Model‐based analysis of genome‐wide data sets allows discrimination between alternative scenarios for these processes even in nonmodel taxa. We used two complementary approaches (one based on the blockwise site frequency spectrum [bSFS], the second on the pairwise sequentially Markovian coalescent [PSMC]) to infer the divergence history of a fig wasp, Pleistodontes nigriventris. Pleistodontes nigriventris and its fig tree mutualist Ficus watkinsiana are restricted to rain forest patches along the eastern coast of Australia and are separated into The Northern population is to the north of the Southern populations by two dry forest corridors (the Burdekin and St. Lawrence Gaps). We generated whole genome sequence data for two haploid males per population and used the bSFS approach to infer the timing of divergence between northern and southern populations of P. nigriventris, and to discriminate between alternative isolation with migration (IM) and instantaneous admixture (ADM) models of postdivergence gene flow. Pleistodontes nigriventris has low genetic diversity (π = 0.0008), to our knowledge one of the lowest estimates reported for a sexually reproducing arthropod. We find strongest support for an ADM model in which the two populations diverged ca. 196 kya in the late Pleistocene, with almost 25% of northern lineages introduced from the south during an admixture event ca. 57 kya. This divergence history is highly concordant with individual population demographies inferred from each pair of haploid males using PSMC. Our analysis illustrates the inferences possible with genome‐level data for small population samples of tiny, nonmodel organisms and adds to a growing body of knowledge on the population structure of Australian rain forest taxa.     

Population history of Berthelot's pipit: colonization, gene flow and morphological divergence in Macaronesia

The fauna of oceanic islands provide exceptional models with which to examine patterns of dispersal, isolation and diversification, from incipient speciation to species level radiations. Here, we investigate recent differentiation and microevolutionary change in Berthelot's pipit (Anthus berthelotii), an endemic bird species inhabiting three Atlantic archipelagos. Mitochondrial DNA sequence data and microsatellite markers were used to deduce probable colonization pathway, genetic differentiation, and gene flow among the 12 island populations. Phenotypic differentiation was investigated based on eight biologically important morphological traits. We found little mitochondrial DNA variability, with only one and four haplotypes for the control region and cytochrome b, respectively. However, microsatellite data indicated moderate population differentiation (FST=0.069) between the three archipelagos that were identified as genetically distinct units with limited gene flow. Both results, combined with the estimated time of divergence (2.5 millions years ago) from the Anthus campestris (the sister species), suggest that this species has only recently dispersed throughout these islands. The genetic relationships, patterns of allelic richness and exclusive alleles among populations suggest the species originally colonized the Canary Islands and only later spread from there to the Madeiran archipelago and Selvagen Islands. Differentiation has also occurred within archipelagos, although to a lesser degree. Gene flow was observed more among the eastern and central islands of the Canaries than between these and the western islands or the Madeiran Islands. Morphological differences were also more important between than within archipelagos. Concordance between morphological and genetic differentiation provided ambiguous results suggesting that genetic drift alone was not sufficient to explain phenotypic differentiation. The observed genetic and morphological differences may therefore be the result of differing patterns of selection pressures between populations, with Berthelot's pipit undergoing a process of incipient differentiation.     

Transcription in space – environmental vs. genetic effects on differential immune gene expression

Understanding how organisms adapt to their local environment is one of the key goals in molecular ecology. Adaptation can be achieved through qualitative changes in the coding sequence and/or quantitative changes in gene expression, where the optimal dosage of a gene's product in a given environment is being selected for. Differences in gene expression among populations inhabiting distinct environments can be suggestive of locally adapted gene regulation and have thus been studied in different species (Whitehead & Crawford 2006 ; Hodgins‐Davis & Townsend 2009 ). However, in contrast to a gene's coding sequence, its expression level at a given point in time may depend on various factors, including the current environment. Although critical for understanding the extent of local adaptation, it is usually difficult to disentangle the heritable differences in gene regulation from environmental effects. In this issue of Molecular Ecology, Stutz et al. ( 2015 ) describe an experiment in which they reciprocally transplanted three‐spined sticklebacks (Gasterosteus aculeatus) between independent pairs of small and large lakes. Their experimental design allows them to attribute differences in gene expression among sticklebacks either to lake of origin or destination lake. Interestingly, they find that translocated sticklebacks show a pattern of gene expression more similar to individuals from the destination lake than to individuals from the lake of origin, suggesting that expression of the targeted genes is more strongly regulated by environmental effects than by genetics. The environmental effect by itself is not entirely surprising; however, the relative extent of it is. Especially when put in the context of local adaptation and population differentiation, as done here, these findings cast a new light onto the heritability of differential gene expression and specifically its relative importance during population divergence and ultimately ecological speciation.     

Speciational history of Australian grass finches (Poephila) inferred from thirty gene trees

Multilocus genealogical approaches are still uncommon in phylogeography and historical demography, fields which have been dominated by microsatellite markers and mitochondrial DNA, particularly for vertebrates. Using 30 newly developed anonymous nuclear loci, we estimated population divergence times and ancestral population sizes of three closely related species of Australian grass finches (Poephila) distributed across two barriers in northern Australia. We verified that substitution rates were generally constant both among lineages and among loci, and that intralocus recombination was uncommon in our dataset, thereby satisfying two assumptions of our multilocus analysis. The reconstructed gene trees exhibited all three possible tree topologies and displayed considerable variation in coalescent times, yet this information provided the raw data for maximum likelihood and Bayesian estimation of population divergence times and ancestral population sizes. Estimates of these parameters were in close agreement with each other regardless of statistical approach and our Bayesian estimates were robust to prior assumptions. Our results suggest that black-throated finches (Poephila cincta) diverged from long-tailed finches (P. acuticauda and P. hecki) across the Carpentarian Barrier in northeastern Australia around 0.6 million years ago (mya), and that P. acuticauda diverged from P. hecki across the Kimberley Plateau-Arnhem Land Barrier in northwestern Australia approximately 0.3 mya. Bayesian 95% credibility intervals around these estimates strongly support Pleistocene timing for both speciation events, despite the fact that many gene divergences across the Carpentarian region clearly predated the Pleistocene. Estimates of ancestral effective population sizes for the basal ancestor and long-tailed finch ancestor were large (about 521,000 and about 384,000, respectively). Although the errors around the population size parameter estimates are considerable, they are the first for birds taking into account multiple sources of variance.     

Flickers of speciation: Sympatric colour morphs of the arc‐eye hawkfish,Paracirrhites arcatus,reveal key elements of divergence with gene flow

One of the primary challenges of evolutionary research is to identify ecological factors that favour reproductive isolation. Therefore, studying partially isolated taxa has the potential to provide novel insight into the mechanisms of evolutionary divergence. Our study utilizes an adaptive colour polymorphism in the arc‐eye hawkfish (Paracirrhites arcatus) to explore the evolution of reproductive barriers in the absence of geographic isolation. Dark and light morphs are ecologically partitioned into basaltic and coral microhabitats a few metres apart. To test whether ecological barriers have reduced gene flow among dark and light phenotypes, we evaluated genetic variation at 30 microsatellite loci and a nuclear exon (Mc1r) associated with melanistic coloration. We report low, but significant microsatellite differentiation among colour morphs and stronger divergence in the coding region of Mc1r indicating signatures of selection. Critically, we observed greater genetic divergence between colour morphs on the same reefs than that between the same morphs in different geographic locations. We hypothesize that adaptation to the contrasting microhabitats is overriding gene flow and is responsible for the partial reproductive isolation observed between sympatric colour morphs. Combined with complementary studies of hawkfish ecology and behaviour, these genetic results indicate an ecological barrier to gene flow initiated by habitat selection and enhanced by assortative mating. Hence, the arc‐eye hawkfish fulfil theoretical expectations for the earliest phase of speciation with gene flow.     

Genomic differentiation and patterns of gene flow between two long‐tailed tit species (Aegithalos)

Patterns of heterogeneous genomic differentiation have been well documented between closely related species, with some highly differentiated genomic regions (“genomic differentiation islands”) spread throughout the genome. Differential levels of gene flow are proposed to account for this pattern, as genomic differentiation islands are suggested to be resistant to gene flow. Recent studies have also suggested that genomic differentiation islands could be explained by linked selection acting on genomic regions with low recombination rates. Here, we investigate genomic differentiation and gene‐flow patterns for autosomes using RAD‐seq data between two closely related species of long‐tailed tits (Aegithalos bonvaloti and A. fuliginosus) in both allopatric and contact zone populations. The results confirm recent or ongoing gene flow between these two species. However, there is little evidence that the genomic regions that were found to be highly differentiated between the contact zone populations are resistant to gene flow, suggesting that differential levels of gene flow is not the cause of the heterogeneous genomic differentiation. Linked selection may be the cause of genomic differentiation islands between the allopatric populations with no or very limited gene flow, but this could not account for the heterogeneous genomic differentiation between the contact zone populations, which show evidence of recent or ongoing gene flow.     

Isolation and gene flow: inferring the speciation history of European house mice

Inferring the history of isolation and gene flow during species differentiation can inform us on the processes underlying their formation. Following their recent expansion in Europe, two subspecies of the house mouse (Mus musculus domesticus and Mus musculus musculus) have formed a hybrid zone maintained by hybrid incompatibilities and possibly behavioural reinforcement, offering a good model of incipient speciation. We reconstruct the history of their divergence using an approximate Bayesian computation framework and sequence variation at 57 autosomal loci. We find support for a long isolation period preceding the advent of gene flow around 200,000 generations ago, much before the formation of the European hybrid zone a few thousand years ago. The duration of the allopatric episode appears long enough (74% of divergence time) to explain the accumulation of many post-zygotic incompatibilities expressed in the present hybrid zone. The ancient contact inferred could have played a role in mating behaviour divergence and laid the ground for further reinforcement. We suggest that both subspecies originally colonized the Middle East from the northern Indian subcontinent, domesticus settling on the shores of the Persian Gulf and musculus on those of the Caspian Sea. Range expansions during interglacials would have induced secondary contacts, presumably in Iran, where they must have also interacted with Mus musculus castaneus. Future studies should incorporate this possibility, and we point to Iran and its surroundings as a hot spot for house mouse diversity and speciation studies.     

Genomic evidence for convergent evolution of a key trait underlying divergence in island birds

Reproductive isolation can be initiated by changes in one or a few key traits that prevent random mating among individuals in a population. During the early stages of speciation, when isolation is often incomplete, there will be a heterogeneous pattern of differentiation across regions of the genome between diverging populations, with loci controlling these key traits appearing the most distinct as a result of strong diversifying selection. In this study, we used Illumina‐sequenced ddRAD tags to identify genomewide patterns of differentiation in three recently diverged island populations of the Monarcha castaneiventris flycatcher of the Solomon Islands. Populations of this species have diverged in plumage colour, and these differences in plumage colour, in turn, are used in conspecific recognition and likely important in reproductive isolation. Previous candidate gene sequencing identified point mutations in MC1R and ASIP, both known pigmentation genes, to be associated with the difference in plumage colour between islands. Here, we show that background levels of genomic differentiation based on over 70,000 SNPs are extremely low between populations of distinct plumage colour, with no loci reaching the level of differentiation found in either candidate gene. Further, we found that a phylogenetic analysis based on these SNPs produced a taxonomy wherein the two melanic populations appear to have evolved convergently, rather than from a single common ancestor, in contrast to their original classification as a single subspecies. Finally, we found evidence that the pattern of low genomic differentiation is the result of both incomplete lineage sorting and gene flow between populations.     

A parapatric propensity for breeding precludes the completion of speciation in common teal (Anas crecca, sensu lato)

Speciation is a process in which genetic drift and selection cause divergence over time. However, there is no rule dictating the time required for speciation, and even low levels of gene flow hinder divergence, so that taxa may be poised at the threshold of speciation for long periods of evolutionary time. We sequenced mitochondrial DNA (mtDNA) and eight nuclear introns (nuDNA) to estimate genomic levels of differentiation and gene flow between the Eurasian common teal (Anas crecca crecca) and the North American green‐winged teal (Anas crecca carolinensis). These ducks come into contact in Beringia (north‐eastern Asia and north‐western North America) and have probably done so, perhaps cyclically, since the Pliocene–Pleistocene transition, ～2.6 Ma, when they apparently began diverging. They have diagnosable differences in male plumage and are 6.9% divergent in the mtDNA control region, with only 1 of 58 crecca and 2 of 86 carolinensis having haplotypes grouping with the other. Two nuclear loci were likewise strongly structured between these teal (Φ_st ≥ 0.35), but six loci were undifferentiated or only weakly structured (Φ_st = 0.0–0.06). Gene flow between crecca and carolinensis was ～1 individual per generation in both directions in mtDNA, but was asymmetrical in nuDNA, with ～1 and ～20 individuals per generation immigrating into crecca and carolinensis, respectively. This study illustrates that species delimitation using a single marker oversimplifies the complexity of the speciation process, and it suggests that even with divergent selection, moderate levels of gene flow may stall the speciation process short of completion.     

The influence of gene flow and drift on genetic and phenotypic divergence in two species of Zosterops in Vanuatu

Colonization of an archipelago sets the stage for adaptive radiation. However, some archipelagos are home to spectacular radiations, while others have much lower levels of diversification. The amount of gene flow among allopatric populations is one factor proposed to contribute to this variation. In island colonizing birds, selection for reduced dispersal ability is predicted to produce changing patterns of regional population genetic structure as gene flow-dominated systems give way to drift-mediated divergence. If this transition is important in facilitating phenotypic divergence, levels of genetic and phenotypic divergence should be associated. We consider population genetic structure and phenotypic divergence among two co-distributed, congeneric (Genus: Zosterops) bird species inhabiting the Vanuatu archipelago. The more recent colonist, Z. lateralis, exhibits genetic patterns consistent with a strong influence of distance-mediated gene flow. However, complex patterns of asymmetrical gene flow indicate variation in dispersal ability or inclination among populations. The endemic species, Z. flavifrons, shows only a partial transition towards a drift-mediated system, despite a long evolutionary history on the archipelago. We find no strong evidence that gene flow constrains phenotypic divergence in either species, suggesting that levels of inter-island gene flow do not explain the absence of a radiation across this archipelago.     

Divergence,gene flow,and speciation in eight lineages of trans‐Beringian birds

Determining how genetic diversity is structured between populations that span the divergence continuum from populations to biological species is key to understanding the generation and maintenance of biodiversity. We investigated genetic divergence and gene flow in eight lineages of birds with a trans‐Beringian distribution, where Asian and North American populations have likely been split and reunited through multiple Pleistocene glacial cycles. Our study transects the speciation process, including eight pairwise comparisons in three orders (ducks, shorebirds and passerines) at population, subspecies and species levels. Using ultraconserved elements (UCEs), we found that these lineages represent conditions from slightly differentiated populations to full biological species. Although allopatric speciation is considered the predominant mode of divergence in birds, all of our best divergence models included gene flow, supporting speciation with gene flow as the predominant mode in Beringia. In our eight lineages, three were best described by a split‐migration model (divergence with gene flow), three best fit a secondary contact scenario (isolation followed by gene flow), and two showed support for both models. The lineages were not evenly distributed across a divergence space defined by gene flow (M) and differentiation (F_ST), instead forming two discontinuous groups: one with relatively shallow divergence, no fixed single nucleotide polymorphisms (SNPs), and high rates of gene flow between populations; and the second with relatively deeply divergent lineages, multiple fixed SNPs, and low gene flow. Our results highlight the important role that gene flow plays in avian divergence in Beringia.     

Fine‐grained adaptive divergence in an amphibian: genetic basis of phenotypic divergence and the role of nonrandom gene flow in restricting effective migration among wetlands

Adaptive ecological differentiation among sympatric populations is promoted by environmental heterogeneity, strong local selection and restricted gene flow. High gene flow, on the other hand, is expected to homogenize genetic variation among populations and therefore prevent local adaptation. Understanding how local adaptation can persist at the spatial scale at which gene flow occurs has remained an elusive goal, especially for wild vertebrate populations. Here, we explore the roles of natural selection and nonrandom gene flow (isolation by breeding time and habitat choice) in restricting effective migration among local populations and promoting generalized genetic barriers to neutral gene flow. We examined these processes in a network of 17 breeding ponds of the moor frog Rana arvalis, by combining environmental field data, a common garden experiment and data on variation in neutral microsatellite loci and in a thyroid hormone receptor (TRβ) gene putatively under selection. We illustrate the connection between genotype, phenotype and habitat variation and demonstrate that the strong differences in larval life history traits observed in the common garden experiment can result from adaptation to local pond characteristics. Remarkably, we found that haplotype variation in the TRβ gene contributes to variation in larval development time and growth rate, indicating that polymorphism in the TRβ gene is linked with the phenotypic variation among the environments. Genetic distance in neutral markers was correlated with differences in breeding time and environmental differences among the ponds, but not with geographical distance. These results demonstrate that while our study area did not exceed the scale of gene flow, ecological barriers constrained gene flow among contrasting habitats. Our results highlight the roles of strong selection and nonrandom gene flow created by phenological variation and, possibly, habitat preferences, which together maintain genetic and phenotypic divergence at a fine‐grained spatial scale.     

Stark sexual display divergence among jumping spider populations in the face of gene flow

Gene flow can inhibit evolutionary divergence by eroding genetic differences between populations. A current aim in speciation research is to identify conditions in which selection overcomes this process. We focused on a state of limited differentiation, asking whether selection enables divergence with gene flow in a set of Habronattus americanus jumping spider populations that exhibit three distinct male sexual display morphs. We found that each population is at high frequency or fixed for a single morph. These strong phenotypic differences contrast with low divergence at 210 AFLP markers, suggesting selection has driven or maintains morph divergence. Coinciding patterns of isolation by distance and ‘isolation by phenotype’ (i.e. increased genetic divergence among phenotypically contrasting populations) across the study area support several alternative demographic hypotheses for display divergence, each of which entails gene flow. Display‐associated structure appears broadly distributed across the genome and the markers producing this pattern do not stand out from background levels of differentiation. Overall, the results suggest selection can promote stark sexual display divergence in the face of gene flow among closely related populations.     

Pleistocene speciation with and without gene flow in Euphaea damselflies of subtropical and tropical East Asian islands

Climatic oscillations during the Pleistocene period could have had a profound impact on the origin of tropical species by the alternation of allopatric isolation and interpopulation gene flow cycles. However, whether tropical speciation involves strictly allopatric isolation, or proceeds in the face of homogenizing gene flow, is relatively unclear. Here, we investigated geographical modes of speciation in four closely related Euphaea damselfly species endemic to the subtropical and tropical East Asian islands using coalescent analyses of a multilocus data set. The reconstructed phylogenies demonstrated distinct species status for each of the four species and the existence of two sister species pairs, Euphaea formosa/E. yayeyamana and E. decorata/E. ornata. The species divergence time of the sibling Euphaea damselflies dates back to within the last one Mya of the Middle to Lower Pleistocene. The speciation between the populous E. formosa of Taiwan and the less numerous E. yayeyamana of the Yaeyama islands occurred despite significant bidirectional, asymmetric gene flow, which is strongly inconsistent with a strictly allopatric model. In contrast, speciation of the approximately equal-sized populations of E. decorata of the southeast Asian mainland and E. ornata of Hainan is inferred to have involved allopatric divergence without gene flow. Our findings suggest that differential selection of natural or sexual environments is a prominent driver of species divergence in subtropical E. formosa and E. yayeyamana; whereas for tropical E. decorata and E. ornata at lower latitudes, allopatric isolation may well be a pivotal promoter of species formation.     

Estimating divergence time with the use of microsatellite genetic distances: impacts of population growth and gene flow

Genetic distances play an important role in estimating divergence time of bifurcated populations. However, they can be greatly affected by demographic processes, such as migration and population dynamics, which complicate their interpretation. For example, the widely used distance for microsatellite loci, (deltamu)2, assumes constant population size, no gene flow, and mutation-drift equilibrium. It is shown here that (deltamu)2 strongly underestimates divergence time if populations are growing and/or connected by gene flow. In recent publications, the average estimate of divergence time between African and non-African populations obtained by using (deltamu)2 is about 34,000 years, although archaeological data show a much earlier presence of modern humans out of Africa. I introduce a different estimator of population separation time based on microsatellite statistics, T(D), that does not assume mutation-drift equilibrium, is independent of population dynamics in the absence of gene flow, and is robust to weak migration flow for growing populations. However, it requires a knowledge of the variance in the number of repeats at the beginning of population separation, V(0). One way to overcome this problem is to find minimal and maximal bounds for the variance and thus obtain the earliest and latest bounds for divergence time (this is not a confidence interval, and it simply reflects an uncertainty about the value of V(0) in an ancestral population). Another way to avoid the uncertainty is to choose from among present populations a reference whose variation is presumably close to what it might have been in an ancestral population. A different approach for using T(D) is to estimate the time difference between adjacent nodes on a phylogenetic population tree. Using data on variation at autosomal short tandem repeat loci with di-, tri-, and tetranucleotide repeats in worldwide populations, T(D) gives an estimate of 57,000 years for the separation of the out-of-Africa branch of modern humans from Africans based on the value of V(0) in the Southern American Indian populations; the earliest bound for this event has been estimated to be about 135,000 years. The data also suggest that the Asian and European populations diverged from each other about 20,000 years, after the occurrence of the out-of-Africa branch.