Using nested clade analysis to assess the history of colonization and the persistence of populations of an Iberian Lizard

The distribution of the lizard Lacerta schreiberi is likely to have been severely affected by the climatic cycles that have influenced the Iberian Peninsula. Information about the species ecology and Iberian physiogeography was used to generate specific hypotheses about episodes of colonization and subsequent population persistence. These hypotheses generated predictions about the distribution of genetic variation, which were tested using nested clade analysis (NCA) supplemented by analysis of molecular variance (amova). Two predictions were confirmed by NCA; that is those that specified multiple and allopatric refugia. However, the remaining three predictions were not corroborated by the analyses. Firstly, a simple analysis of the distribution of genetic variability failed to detect an expected difference in the pattern of colonization between the inland mountain system and the coastal region. Moreover, while NCA did detect the expected genetic pattern in southern coastal populations, it was explained in terms of long-distance migration, which seems implausible because of the extent of unsuitable habitat. A more likely cause of the pattern is population fragmentation and a reduction in population size caused during the Holocene. Finally, NCA also failed to detect a northwestern population expansion, which is supported by other evidence. We conclude that NCA has a limited ability to detect range expansion led by individuals with more ancestral (interior) haplotypes.     

Phylogeography of Cavernularia hultenii: evidence of slow genetic drift in a widely disjunct lichen

Population structure and history is poorly known in most lichenized ascomycetes. Many species display large-scale infraspecific disjunctions, which have been explained alternately by range fragmentation in species of high age and widespread long-distance dispersal. Using the lichen Cavernularia hultenii, which is widely disjunct across North America and Europe, Pleistocene and Holocene population history was inferred. The internal transcribed spacer (ITS) and part of the the intergenic spacer (IGS) region of the nuclear ribosomal DNA were sequenced in 300 individuals representing 62 populations across the range of the species. While four ancestral haplotypes are found in all areas, none of the observed tip haplotypes is present in more than one of the three part ranges. Although this is evidence for a past fragmentation event, nested clade analysis (NCA) remains equivocal in the choice between allopatric fragmentation and long-distance dispersal. Mismatch distributions indicate exponential population growth, probably during postglacial invasion of C. hultenii into formerly glaciated areas of western North America. The presence of one southern and at least one northern glacial refugium in South Central Alaska is inferred. Evidence for another refugium in the Queen Charlotte Islands or Alexander Archipelago is inconclusive because of sparse sampling. However, a range expansion was not confirmed unambiguously by NCA. The limited power of NCA to infer past range fragmentations and expansions is due apparently to the shallow haplotype network and widespread ancestral haplotypes. This can be explained by slow genetic drift causing incomplete removal of ancestral haplotypes from the postfragmentation and postexpansion areas.     

Phylogeographic structure suggests multiple glacial refugia in northern Victoria Land for the endemic Antarctic springtail Desoria klovstadi (Collembola, Isotomidae)

We carried out a phylogeographic study using mtDNA (COII) for the endemic springtail Desoria klovstadi (formerly Isotoma klovstadi ) from northern Victoria Land, Antarctica. Low levels of sequence divergence (≤ 1.6%) across 26 unique haplotypes (from 69 individuals) were distributed according to geographic location. Cape Hallett and Daniell Peninsula contained the highest nucleotide (both > 0.004) and haplotype (both > 0.9) diversity with 10 (of 16) and 8 (of 12) unique haplotypes, respectively. All other populations (Football Saddle, Crater Cirque, Cape Jones) had lower diversity with 2–4 unique haplotypes. Across the 69 individuals from five populations there was only a single haplotype shared between two populations (Daniell Peninsula and Football Saddle). Furthermore, nested clade analyses revealed that some of the Daniell Peninsula haplotypes were more closely related to Football Saddle haplotypes than to any other population. Such discrete haplotype groupings suggest historical (rare) dispersal across the Pleistocene (1.8 mya−11 kya) and Holocene (11 kya–present), coupled with repeated extinction, range contraction and expansion events, and/or incomplete sampling across the species range. The nested clade analyses reveal that a common pattern of climatic and geological history over long-term glacial habitat fragmentation has determined the geographic and haplotype distributions found for D. klovstadi .     

Evolutionary history of the European whitefish Coregonus lavaretus (L.) species complex as inferred from mtDNA phylogeography and gill-raker numbers

We compared mitochondrial DNA and gill-raker number variation in populations of the European whitefish Coregonus lavaretus (L.) species complex to illuminate their evolutionary history, and discuss mechanisms behind diversification. Using single-strand conformation polymorphism (SSCP) and sequencing 528 bp of combined parts of the cytochrome oxidase b (cyt b) and NADH dehydrogenase subunit 3 (ND3) mithochondrial DNA (mtDNA) regions, we documented phylogeographic relationships among populations and phylogeny of mtDNA haplotypes. Demographic events behind geographical distribution of haplotypes were inferred using nested clade analysis (NCA) and mismatch distribution. Concordance between operational taxonomical groups, based on gill-raker numbers, and mtDNA patterns was tested. Three major mtDNA clades were resolved in Europe: a North European clade from northwest Russia to Denmark, a Siberian clade from the Arctic Sea to southwest Norway, and a South European clade from Denmark to the European Alps, reflecting occupation in different glacial refugia. Demographic events inferred from NCA were isolation by distance, range expansion, and fragmentation. Mismatch analysis suggested that clades which colonized Fennoscandia and the Alps expanded in population size 24 500-5800 years before present, with minute female effective population sizes, implying small founder populations during colonization. Gill-raker counts did not commensurate with hierarchical mtDNA clades, and poorly with haplotypes, suggesting recent origin of gill-raker variation. Whitefish designations based on gill-raker numbers were not associated with ancient clades. Lack of congruence in morphology and evolutionary lineages implies that the taxonomy of this species complex should be reconsidered.     

Dispersal, habitat differences, and comparative phylogeography of Southeast Asian seahorses (Syngnathidae: Hippocampus)

Four distinct phylogeographical patterns across Southeast Asia were observed for four species of seahorse (genus Hippocampus) with differing ecologies. For all species, genetic differentiation (based on cytochrome b sequence comparisons) was significantly associated with sample site (Phi(ST) = 0.190-0.810, P < 0.0001) and with geographical distance (Mantel's r = 0.37-0.59, P < 0.019). Geographic locations of genetic breaks were inconsistent across species in 7/10 comparisons, although some similarities across species were also observed. The two shallow-water species (Hippocampus barbouri and Hippocampus kuda) have colonized the Sunda Shelf to a lesser degree than the two deeper-water species (Hippocampus spinosissimus and Hippocampus trimaculatus). In all species the presence of geographically restricted haplotypes in the Philippines could indicate past population fragmentation and/or long-distance colonization. A nested clade analysis (NCA) revealed that long-distance colonization and/or fragmentation were likely the dominant forces that structure populations of the two shallow-water species, whereas range expansion and restricted dispersal with isolation by distance were proportionally more important in the history of the two deeper-water species. H. trimaculatus has the most widespread haplotypes [average clade distance (D(c)) of nonsingleton haplotypes = 1169 km], indicating potentially high dispersal capabilities, whereas H. barbouri has the least widespread haplotypes (average D(c) = 67 km) indicating potentially lower dispersal capabilities. Pleistocene separation of marine basins and postglacial flooding of the Sunda Shelf are extrinsic factors likely to have contributed to the phylogeographical structure observed, whereas differences among the species appear to reflect their individual ecologies.     

Phylogeography of the leaf beetle Chrysolina virgata in wetlands of Japan inferred from the distribution of mitochondrial haplotypes

The genetic differentiation among populations of the leaf beetle Chrysolina virgata living in wetlands of Japan was studied based on the sequence data of the mitochondrial cytochrome oxidase subunit I gene region (750 bp). Two distinct lineages of mitochondrial haplotypes were found: one (clade A) consisted of 26 haplotypes distributed over the distribution range of C. virgata between north‐east Honshu and Kyushu, whereas the other (clade B) was monotypic and confined to a small region in north‐east Honshu where it coexisted with clade A. Nested clade analysis for these haplotypes suggested that range expansion and following differentiation due to isolation by distance might have resulted in the present distribution pattern of the haplotypes in clade A. We discuss the evolutionary process leading to the occurrence of two distinct haplotype clades in Japan in terms of repeated colonization from the continent and range expansion and contraction during climatic changes.     

Population Genetics and Phylogeography of Bufo gargarizans in China

The Chinese toad Bufo gargarizans has an extensive range, covering nearly all of China. This study addresses the population genetic structure of this toad and the relation to geographic location. Partial sequences of the mtDNA control region were obtained from 14 populations, 29 haplotypes were defined, and 73 variable sites were found, including 66 parsimony informative sites. All population genetic analyses indicated no clear geographic pattern in the distribution of haplotypes. The genetic divergence between the populations was significant. Phylogeographic analyses suggested that past fragmentation and/or long-distance colonization seemed to have shaped the present-day distribution of the haplotypes of B. gargarizans.     

Phylogeographic study reveals the origin and evolutionary history of a Rhododendron species complex in Taiwan

This study infers a single origin and a once-widespread distribution of the Rhododendron pseudochrysanthum species complex in Taiwan based on chloroplast DNA (cpDNA) variation. In total, 124 individuals from five endemic Rhododendron species were used for amplifications of two chloroplast intergenic spacers: trnL-trnF and atpB-rbcL. The haplotype and nucleotide diversities were much lower for the R. pseudochrysanthum complex, comprised of the species R. pseudochrysanthum, R. morii, Rhododendron rubropunctatum, and Rhododendron hyperythrum, than for Rhododendron formosanum. Two measures of pairwise population differentiation, N(ST) and F(ST), consistently revealed mostly non-significant levels of genetic divergence between populations of the R. pseudochrysanthum complex. No genetic difference was found among the four species of the R. pseudochrysanthum complex by analysis of molecular variance (AMOVA), which is concordant with the parsimonious topology of cpDNA haplotypes for the complex. Nested clade analysis (NCA) of the cpDNA haplotypes indicated that restricted gene flow with isolation-by-distance characterized the recolonization pattern of the R. pseudochrysanthum complex. In contrast, the NCA analysis indicated a contiguous range expansion for cpDNA haplotypes of R. formosanum. This research suggests a once-widespread distribution of the R. pseudochrysanthum complex probably via north-to-south colonization of mid-elevations during low-temperature periods of the Pleistocene. Population fragmentation followed the warmer climate which began in the Holocene and resulted in the present-day range contraction into high elevations.     

Complex patterns of colonization and refugia revealed for European grayling Thymallus thymallus,based on complete sequencing of the mitochondrial DNA control region

The complete mitochondrial DNA (mtDNA) control region (1043 base pairs) and 162-bp of flanking transfer RNA genes were sequenced in 316 European grayling, Thymallus thymallus, from 44 populations throughout the Western European range of the species. A total of 58 haplotypes were revealed with pairwise divergence ranging from 0.001 to 0.038. An inferred intraspecific phylogenetic tree revealed two well-supported clades within the Danube basin, one highly divergent clade in the Adriatic basin, and one large, diverse group representing most other populations. A deeply divergent haplotype fixed in the Loire basin in central France, more groups of haplotypes from distinct Danubian tributaries, and a relatively ancestral haplotype fixed in former tributaries of the Elbe in Denmark all suggest a complex pattern of interglacial and postglacial expansions originating from disjunct refugia throughout central Europe. Despite some evidence of human-mediated stock transfers, parsimony-network-based nested-clade analysis (NCA) supported specific inferences relating to corridors of postglacial expansion such as the lower Rhine (Moselle) and Elbe systems (Danish populations) serving as sources for expansion into the Baltic to the north as well as the upper Rhine and Danube to the south; and specific Rhine populations (Doller, Orbe and Reuss) serving as sources for colonization of the Rhone. The multiple divergent clades representing populations in the upper Danube, as well as the deeply divergent haplotypes found in the Adriatic and Loire basins (> 5% divergence from Asian outgroups) support the theory that European grayling have had a long history in Western Europe, pre-dating Pleistocene glacial cycles. The patterns of mtDNA divergence shown here support a perspective of rich inter- and intrabasin genetic diversity that should be protected from current trends to translocate brood stocks for rearing and release in response to declining populations, especially in southern European basins.     

Phylogeography of the rice frog, Fejervarya multistriata (Anura: Ranidae), from China based on mtDNA D-loop sequences

The rice frog, Fejervarya multistriata, is an amphibian widely distributed in China. In this study, we sampled the species across its distributional area in China and sequenced the mtDNA D-loop to investigate the genetic diversity and geographical pattern of the frog population. The results revealed 38 haplotypes in the population, with K2P values varying from 0.19% to 4.22%. Both a phylogenetic analysis and a nested clade analysis (NCA) detected two geographically isolated lineages respectively distributed around the Yangtze drainage (Yangtze lineage) and the south of China (southern lineage). NCA inferred a contiguous range expansion within the Yangtze lineage and allopatric fragmentation within the southern lineage, which might be partly due to the limited samples from this lineage. Accordingly, Fu's Fs test also indicated a population expansion after glacial movement. Therefore, we assumed that the species history responding to glacial events shaped the present population pattern of F. multistriata on the Chinese mainland.     

Phylogeographic analysis of the endemic species Sibiraea angustata reveals a marginal refugium in the Qinghai–Tibet Plateau

The phylogeography of Sibiraea angustata, an endemic shrub species, was studied in the Qinghai–Tibet plateau (QTP). We investigated 466 individuals of S. angustata from 39 populations basically covering its total distribution area. Eight haplotypes (A–H) were detected by sequencing the intergenic chloroplast spacer trnS–trnG (600 bp), and one ancestral haplotype (A) was found to be widely distributed. The level of differentiation among populations was very high (G_ST=0.768; N_ST=0.850) and a significant phylogeographical structure was revealed (N_ST>G_ST). Analysis of molecular variance (AMOVA) similarily revealed a high level of differentiation among populations (84%, F_ST=0.842), indicating that little gene flow has occurred among populations mutually isolated by high mountains and rivers in the QTP. On the QTP platform there was only one widespread haplotype (A) in most populations, while populations along the eastern and southeastern edges had high diversity and unique haplotypes. Our results suggest that a glacial refugium may have been located on the eastern or southeastern edges of QTP during the last glaciation, and that interglacial and postglacial range expansion occurred from that refugium. Nested clade analysis (NCA) also suggests this scenario, which indicates that the current spatial distribution of cpDNA haplotypes and populations mainly resulted from long distance colonization, possibly coupled with subsequent or past fragmentation followed by range expansion and allopatric fragmentation.     

Genetic structure of the star sea squirt, Botryllus schlosseri, introduced in southern European harbours

The introduction of new genetic variants or species is often caused by maritime transport between harbours. Botryllus schlosseri is a cosmopolitan ascidian species that is found in both harbours and open shore habitats. In order to determine the influence of ship traffic on the genetic structure and phylogeography of B. schlosseri in southern Europe, we analyzed the variability of a fragment of the mitochondrial gene cytochrome c oxidase subunit I (COI). We sampled seven Atlanto-Mediterranean harbour populations and three open-shore populations. In addition, we sequenced some colonies from the US-Atlantic coast and from other Mediterranean localities to perform phylogenetic analyses. Although the number of polymorphic sites recorded (25.8%) was within the range observed in other population studies based on ascidian COI sequences, the haplotypic diversity (16 haplotypes out of 181 sequences) was much lower. Moreover, a lack of intermediate haplotypes was observed. This pattern of high nucleotide diversity and low haplotype diversity was consistent with introduction events of a few divergent haplotypes. We found a strong genetic structure in the study populations. Gene flow was only appreciable between some harbour populations. Harbour- and open-shore populations were well differentiated, although there was no evidence for isolation by distance. A nested clade analysis pointed to long-distance colonization, possibly coupled with subsequent fragmentation, as the underlying process. Our results suggest that B. schlosseri entered the study area via harbour-hopping, possibly through recurrent introduction events. The haplotypes from North America and most of the European ones were grouped in the same phylogenetic clade. This suggests occasional gene flow between both continents, probably through ship transport.     

Cryptic lineages and Pleistocene population expansion in a Brazilian Cerrado frog

Diversification of South American species endemic to open habitats has been attributed to both Tertiary events and Pleistocene climatic fluctuations. Nonetheless, phylogeographical studies of taxa in these regions are few, precluding generalizations about the timing and processes leading to differentiation and speciation. We inferred population structure of Hypsiboas albopunctatus, a frog widely distributed in the Brazilian Cerrado. Three geographically distinct lineages were recovered in our phylogeny. The Chapada dos Guimarães (CG) clade was the first to diverge from other populations and contains multiple haplotypes from a single population in western Cerrado, probably representing a cryptic species. The southeast clade (SE) includes populations along the southeastern limit of the range within the historical distribution of the Brazilian Atlantic forest. Finally, the Central Cerrado (CC) group includes haplotypes from the interior of Brazil that are paraphyletic relative to the SE clade. Analyses of historical demography indicate significant population expansion in the CC and SE populations, likely associated with colonization of newly formed open habitats. The divergence of populations in the CG clade occurred in the late Miocene, concordant with the uplift of the central Brazilian plateau. Divergence of the SE clade from the CC occurred during the mid‐Pleistocene. Thus, both Tertiary geological events and Pleistocene climatic fluctuations promoted divergences among lineages. Our study reveals a complex history of diversification in the Cerrado, a morphoclimatic domain highly threatened because of anthropogenic habitat alteration. We identified surprisingly deep divergences in a widely distributed frog, indicating that the Cerrado is not a barrier‐free habitat and that its diversity is likely underestimated.     

Genetic variability of Aedes aegypti in the Americas using a mitochondrial gene: evidence of multiple introductions

To analyze the genetic relatedness and phylogeographic structure of Aedes aegypti, we collected samples from 36 localities throughout the Americas (Brazil, Peru, Venezuela, Guatemala, US), three from Africa (Guinea, Senegal, Uganda), and three from Asia (Singapore, Cambodia, Tahiti). Amplification and sequencing of a fragment of the mitochondrial NADH dehydrogenase subunit 4 gene identified 20 distinct haplotypes, of which 14 are exclusive to the Americas, four to African/Asian countries, one is common to the Americas and Africa, and one to the Americas and Asia. Nested clade analysis (NCA), pairwise distribution, statistical parsimony, and maximum parsimony analyses were used to infer evolutionary and historic processes, and to estimate phylogenetic relationships among haplotypes. Two clusters were found in all the analyses. Haplotypes clustered in the two clades were separated by eight mutational steps. Phylogeographic structure detected by the NCA was consistent with distant colonization within one clade and fragmentation followed by range expansion via long distance dispersal in the other. Three percent of nucleotide divergence between these two clades is suggestive of a gene pool division that may support the hypothesis of occurrence of two subspecies of Ae. aegypti in the Americas.     

Phylogeography of Japanese water crowfoot based on chloroplast DNA haplotypes

Water crowfoot, Ranunculus subgenus Batrachium, are submerged macrophytes in the Ranunculaceae. We aimed to infer the Quaternary history of these aquatic macrophytes in the Japanese archipelago. We studied 212 individuals of three perennial and one annual species from 46 populations covering the entire geographic range and found eight haplotypes based on approximately 1800 bp of four spacers in chloroplast DNA. The relationships among haplotypes were resolved using maximum parsimony and parsimony network analyses. To identify the zones of clear genetic boundaries, Monmonier's algorithm was used. The eight haplotypes were distinguished from adjacent haplotypes by one substitution or indel. Each of the 46 populations was fixed for a single haplotype. Inconsistency between cpDNA haplotypes and Batrachium taxa was found, except in one annual species. The distribution of the haplotypes in perennial species was highly geographically structured. An abrupt genetic change was detected between the Tohoku region and more southerly regions. Since the perennial Batrachium are cold-adapted, this genetic differentiation may be due to historical changes in their distributions caused by Quaternary climatic oscillations: interglacial retreats in colder refugia and glacial range expansions. The single haplotype composition of each population may have been shaped by founder effects during colonization and/or by genetic drift. The interglacial refugial populations must have been small enough to deplete haplotype diversity. Inconsistency between cpDNA haplotypes and Batrachium taxa may be due to incomplete lineage sorting of ancestral polymorphic haplotypes.     

Genetic structure of desert ground squirrels over a 20-degree-latitude transect from Oregon through the Baja California peninsula

The genetic structure of populations over a wide geographical area should reflect the demographic and evolutionary processes that have shaped a species across its range. We examined the population genetic structure of antelope ground squirrels (Ammospermophilus leucurus) across the complex of North American deserts from the Great Basin of Oregon to the cape region of the Baja California peninsula. We sampled 73 individuals from 13 major localities over this 2500-km transect, from 43 to 22 degrees north. Our molecular phylogeographical analysis of 555 bp of the mitochondrial cytochrome b gene and 510 bp of the control region revealed great genetic uniformity in a single clade that extends from Oregon to central Baja California. A second distinct clade occupies the southern half of the peninsula. The minimal geographical structure of the northern clade, its low haplotype diversity and the distribution of pairwise differences between haplotypes suggest a rapid northward expansion of the population that must have followed a northward desert habitat shift associated with the most recent Quaternary climate warming and glacial retreat. The higher haplotype diversity within the southern clade and distribution of pairwise differences between haplotypes suggest that the southern clade has a longer, more stable history associated with a southern peninsular refugium. This system, as observed, reflects both historical and contemporary ecological and evolutionary responses to physical environmental gradients within genetically homogeneous populations.     

Cryptic Neogene vicariance and Quaternary dispersal of the red-spotted toad (Bufo punctatus): insights on the evolution of North American warm desert biotas

We define the geographical distributions of mitochondrial DNA (mtDNA) lineages embedded within a broadly distributed, arid-dwelling toad, Bufo punctatus. These patterns were evaluated as they relate to hypothesized vicariant events leading to the formation of desert biotas within western North America. We assessed mtDNA sequence variation among 191 samples from 82 sites located throughout much of the species' range. Parsimony-based haplotype networks of major identified lineages were used in nested clade analysis (NCA) to further elucidate and evaluate shallow phylogeographic patterns potentially associated with Quaternary (Pleistocene-Holocene) vicariance and dispersal. Phylogenetic analyses provided strong support for three monophyletic lineages (clades) within B. punctatus. The geographical distributions of the clades showed little overlap and corresponded to the general boundaries of the Peninsular Desert, and two continental desert regions, Eastern (Chihuahuan Desert-Colorado Plateau) and Western (Mojave-Sonoran deserts), geographically separated along the Rocky Mountains and Sierra Madre Occidental. The observed divergence levels and congruence with postulated events in earth history implicate a late Neogene (latest Miocene-early Pliocene) time frame for separation of the major mtDNA lineages. Evaluation of nucleotide and haplotype diversity and interpretations from NCA reveal that populations on the Colorado Plateau resulted from a recent, likely post-Pleistocene, range expansion from the Chihuahuan Desert. Dispersal across historical barriers separating major continental clades appear to be recent, resulting in secondary contacts in at least two areas. Given the observed contact between major clades, we speculated as to why the observed deep phylogeographic structure has not been eroded during the multiple previous interglacials of the Pleistocene.     

Phylogeography of the Vosges mountains populations of Gonioctena pallida (Coleoptera: Chrysomelidae): a nested clade analysis of mitochondrial DNA haplotypes

The pattern of genetic variation in the leaf beetle Gonioctena pallida was investigated inside the Vosges mountains using a highly variable 363 bp DNA fragment of the mitochondrial control region. Sequencing of 242 individuals, sampled in a geographical area of 100 x 40 km, identified 61 haplotypes whose genealogy was inferred. The resulting haplotype network exhibits four star-like phylogenies, two of which may be indicative of a population having recently expanded in size from a small number of founders. Nested clade analysis suggested multiple past expansion events, but also isolation by distance and possibly past fragmentation events, as the causes of the detected geographical associations of haplotypes. These results indicate the existence of effective barriers to gene flow inside the investigated area. Because the oldest demographic events inferred in the nested clade analysis were identified as expansion events, we hypothesize that a small population of founders have expanded not only in size, but also in geographical range from the south towards the north and east of the Vosges.     

Hierarchical patterns of global human Y-chromosome diversity.

We examined 43 biallelic polymorphisms on the nonrecombining portion of the Y chromosome (NRY) in 50 human populations encompassing a total of 2,858 males to study the geographic structure of Y-chromosome variation. Patterns of NRY diversity varied according to geographic region and method/level of comparison. For example, populations from Central Asia had the highest levels of heterozygosity, while African populations exhibited a higher level of mean pairwise differences among haplotypes. At the global level, 36% of the total variance of NRY haplotypes was attributable to differences among populations (i.e., Phi(ST) = 0.36). When a series of AMOVA analyses was performed on different groupings of the 50 populations, high levels of among-groups variance (Phi(CT)) were found between Africans, Native Americans, and a single group containing all 36 remaining populations. The same three population groupings formed distinct clusters in multidimensional scaling plots. A nested cladistic analysis (NCA) demonstrated that both population structure processes (recurrent gene flow restricted by isolation by distance and long-distance dispersals) and population history events (contiguous range expansions and long-distance colonizations) were instrumental in explaining this tripartite division of global NRY diversity. As in our previous analyses of smaller NRY data sets, the NCA detected a global contiguous range expansion out of Africa at the level of the total cladogram. Our new results support a general scenario in which, after an early out-of-Africa range expansion, global-scale patterns of NRY variation were mainly influenced by migrations out of Asia. Two other notable findings of the NCA were (1) Europe as a "receiver" of intercontinental signals primarily from Asia, and (2) the large number of intracontinental signals within Africa. Our AMOVA analyses also supported the hypothesis that patrilocality effects are evident at local and regional scales, rather than at intercontinental and global levels. Finally, our results underscore the importance of subdivision of the human paternal gene pool and imply that caution should be exercised when using models and experimental strategies based on the assumption of panmixia.