Multilocus analysis of nucleotide variation and speciation in three closely related Populus (Salicaceae) species

Historical tectonism and climate oscillations can isolate and contract the geographical distributions of many plant species, and they are even known to trigger species divergence and ultimately speciation. Here, we estimated the nucleotide variation and speciation in three closely related Populus species, Populus tremuloides, P. tremula and P. davidiana, distributed in North America and Eurasia. We analysed the sequence variation in six single‐copy nuclear loci and three chloroplast (cpDNA) fragments in 497 individuals sampled from 33 populations of these three species across their geographic distributions. These three Populus species harboured relatively high levels of nucleotide diversity and showed high levels of nucleotide differentiation. Phylogenetic analysis revealed that P. tremuloides diverged earlier than the other two species. The cpDNA haplotype network result clearly illustrated the dispersal route from North America to eastern Asia and then into Europe. Molecular dating results confirmed that the divergence of these three species coincided with the sundering of the Bering land bridge in the late Miocene and a rapid uplift of the Qinghai‐Tibetan Plateau around the Miocene/Pliocene boundary. Vicariance‐driven successful allopatric speciation resulting from historical tectonism and climate oscillations most likely played roles in the formation of the disjunct distributions and divergence of these three Populus species.     

Extreme population structure and high interspecific divergence of the Silene Y chromosome

Previous studies have demonstrated that the diversity of Y-linked genes is substantially lower than that of their X-linked homologs in the plant Silene latifolia. This difference has been attributed to selective sweeps, Muller's ratchet, and background selection, processes that are predicted to severely affect the evolution of the nonrecombining Y chromosome. We studied the DNA diversity of a noncoding region of the homologous genes DD44Y and DD44X, sampling S. latifolia populations from a wide geographical area and also including the closely related species S. dioica, S. diclinis, and S. heuffelii. On the Y chromosome of S. latifolia, we found substantial DNA diversity. Geographical population structure was far higher than on the X chromosome and differentiation between the species was also higher for the Y than for the X chromosome. Our findings indicate that the loss of genetic diversity on the Y chromosome in Silene occurs within local populations rather than within entire species. These results are compatible with background selection, Muller's ratchet, and local selective sweeps, but not with species-wide selective sweeps. The higher interspecific divergence of DD44Y, compared to DD44X, supports the hypothesis that Y chromosome differentiation between incipient species precedes reproductive isolation of the entire genome, forming an early stage in the process of speciation.     

Role of ancient lakes in genetic and phenotypic diversification of freshwater snails

Endemic organisms of ancient lakes have been studied as models to understand processes of speciation and adaptive radiation. However, it remains unclear how ancient lakes play roles in genetic and phenotypic diversity of freshwater mollusks. In the present study, we focus on viviparid freshwater snails in the ancient lakes of East and Southeast Asia (Japan and China) to address this question. Using molecular phylogenetic analyses based on mitochondrial (COI, 16S) and nuclear genes (18S, 28S, H3), we show that patterns of species diversification in viviparid lineages. Colonization to ancient lakes occurred independently in China and Japan at least four times, with subsequent diversification into more than two species within each lake group. Morphological analyses of fossil related viviparids suggest parallel phenotypic evolution occurred in the different lakes and ages. Each lake contained a single lineage which was phenotypically diversified relative to those from other sites. Using genome‐wide SNPs obtained by MIG‐seq, we also examined the genetic structure of three Japanese viviparids, including two endemic species of ancient Lake Biwa. The results suggest that these two species diversified from the population of the third species living in wetlands surrounding the lake. These findings suggest that rapid diversification of lineages and phenotypic divergence can occur in ancient lakes compared to other habitats. Formation of large lakes probably promotes speciation and phenotypic divergence as a result of adaptation into different microhabitats. High numbers of ancient lakes could be a driver of species diversity in Asian viviparid snails.     

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.     

Differential divergence in autosomes and sex chromosomes is associated with intra‐island diversification at a very small spatial scale in a songbird lineage

Recently diverged taxa showing marked phenotypic and ecological diversity provide optimal systems to understand the genetic processes underlying speciation. We used genome‐wide markers to investigate the diversification of the Reunion grey white‐eye (Zosterops borbonicus) on the small volcanic island of Reunion (Mascarene archipelago), where this species complex exhibits four geographical forms that are parapatrically distributed across the island and differ strikingly in plumage colour. One form restricted to the highlands is separated by a steep ecological gradient from three distinct lowland forms which meet at narrow hybrid zones that are not associated with environmental variables. Analyses of genomic variation based on single nucleotide polymorphism data from genotyping‐by‐sequencing and pooled RAD‐seq approaches show that signatures of selection associated with elevation can be found at multiple regions across the genome, whereas most loci associated with the lowland forms are located on the Z sex chromosome. We identified TYRP1, a Z‐linked colour gene, as a likely candidate locus underlying colour variation among lowland forms. Tests of demographic models revealed that highland and lowland forms diverged in the presence of gene flow, and divergence has progressed as gene flow was restricted by selection at loci across the genome. This system holds promise for investigating how adaptation and reproductive isolation shape the genomic landscape of divergence at multiple stages of the speciation process.     

种间杂种染色体配对所揭示的披碱草属植物StY基因组分化及其进化意义

在同倍体(homoploid)植物杂交-分化的物种形成(speciation)过程中,杂交后代与亲本之间的有效生殖隔离是新物种形成的关键。杂种与亲本在时间、空间、生态环境和基因水平上的隔离,保证了杂种后代的分化和稳定,并逐渐形成新物种。为了研究披碱草属(Elymus)含StY基因组四倍体物种的系统演化关系,本文对来自亚洲不同地理分布区的26种披碱草属植物进行了大规模的种间杂交和杂种F1减数分裂染色体配对行为的分析。结果表明各物种之间有不同程度的杂交亲合力,杂交结实率在各杂交组合之间有较大的变异(在4．8％-100％之间);但各物种之间的杂种F1完全不育。证明各物种之间形成了明显的生殖隔离。种间杂种F1减数分裂中期-I染色体配对分析的结果进一步表明,StY基因组随各披碱草属物种地理分布的不同而有不同程度的分化。来自同一分布区(如东亚或西亚分布区之内)物种的StY基因组分化程度较低,但来自不同分布区(如东亚和西亚)物种之间相同的StY基因组具有显著的分化。表明地理隔离对含StY基因组物种的分化起到了十分重要的作用。通过对种间和种内杂种F1的减数分裂异常现象和染色体配对频率变化规律的分析,作者认为细胞学水平的变化,如基因组同源性的分化和染色体结构的变异等都在杂种后代与亲本之间产生生殖隔离并逐渐形成新物种的进化过程中起到了积极的作用。     

Speciation genomics and a role for the Z chromosome in the early stages of divergence between Mexican ducks and mallards

Speciation is a continuous and dynamic process, and studying organisms during the early stages of this process can aid in identifying speciation mechanisms. The mallard (Anas platyrhynchos) and Mexican duck (A. [p.] diazi) are two recently diverged taxa with a history of hybridization and controversial taxonomy. To understand their evolutionary history, we conducted genomic scans to characterize patterns of genetic diversity and divergence across the mitochondrial DNA (mtDNA) control region, 3523 autosomal loci and 172 Z‐linked sex chromosome loci. Between the two taxa, Z‐linked loci (Φ_ST = 0.088) were 5.2 times more differentiated than autosomal DNA (Φ_ST = 0.017) but comparable to mtDNA (Φ_ST = 0.092). This elevated Z differentiation deviated from neutral expectations inferred from simulated data that incorporated demographic history and differences in effective population sizes between marker types. Furthermore, 3% of Z‐linked loci, compared to <0.1% of autosomal loci, were detected as outlier loci under divergent selection with elevated relative (Φ_ST) and absolute (d_XY) estimates of divergence. In contrast, the ratio of Z‐linked and autosomal differentiation among the seven Mexican duck sampling locations was close to 1:1 (Φ_ST = 0.018 for both markers). We conclude that between mallards and Mexican ducks, divergence at autosomal markers is largely neutral, whereas greater divergence on the Z chromosome (or some portions thereof) is likely the product of selection that has been important in speciation. Our results contribute to a growing body of literature indicating elevated divergence on the Z chromosome and its likely importance in avian speciation.     

Population structure and demographic history of a tropical lowland rainforest tree species Shorea parvifolia (Dipterocarpaceae) from Southeastern Asia

Distribution of tropical rainforests in Southeastern Asia has changed over geo-logical time scale, due to movement of tectonic plates and/or global climatic changes. Shorea parvifolia is one of the most common tropical lowland rainforest tree species in Southeastern Asia. To infer population structure and demographic history of S. parvifolia, as indicators of temporal changes in the distribution and extent of tropical rainforest in this region, we studied levels and patterns of nucleotide polymorphism in the following five nuclear gene regions: GapC, GBSSI, PgiC, SBE2, and SODH. Seven populations from peninsular Malaysia, Sumatra, and eastern Borneo were included in the analyses. STRUCTURE analysis revealed that the investigated populations are divided into two groups: Sumatra-Malay and Borneo. Furthermore, each group contained one admixed population. Under isolation with migration model, divergence of the two groups was estimated to occur between late Pliocene (2.6 MYA) and middle Pleistocene (0.7 MYA). The log-likelihood ratio tests of several demographic models strongly supported model with population expansion and low level of migration after divergence of the Sumatra-Malay and Borneo groups. The inferred demographic history of S. parvifolia suggested the presence of a scarcely forested land bridge on the Sunda Shelf during glacial periods in the Pleistocene and predominance of tropical lowland rainforest at least in Sumatra and eastern Borneo.     

FEMALE HETEROGAMETY AND SPECIATION: REDUCED INTROGRESSION OF THE Z CHROMOSOME BETWEEN TWO SPECIES OF NIGHTINGALES

Several lines of evidence suggest that the X chromosome plays a large role in intrinsic postzygotic isolation. The role of the Z chromosome in speciation is much less understood. To explore the role of the Z chromosome in reproductive isolation, we studied nucleotide variation in two closely related bird species, the Thrush Nightingale ( Luscinia luscinia ) and the Common Nightingale ( L. megarhynchos ). These species are isolated by incomplete prezygotic isolation and female hybrid sterility. We sequenced introns of four Z-linked and eight autosomal loci and analyzed patterns of polymorphism and divergence using a divergence-with-gene flow framework. Our results suggest that the nightingale species diverged approximately 1.8 Mya. We found strong evidence of gene flow after divergence in both directions, although more introgression occurred from L. megarhynchos into L. luscinia . Gene flow was significantly higher on the autosomes than on the Z chromosome. Our results support the idea that the Z chromosome plays an important role in intrinsic postzygotic isolation in birds, although it may also contribute to the evolution of prezygotic isolation through sexual selection. This highlights the similarities in the genetic basis of reproductive isolation between organisms with heterogametic males and organisms with heterogametic females during the early stages of speciation.     

Impacts of climate warming on hybrid zone movement: Geographically diffuse and biologically porous “species borders”

Abstract The ecology and evolutionary biology of insect–plant associations has realized extensive attention, especially during the past 60 years. The classifications (categorical designations) of continuous variation in biodiversity, ranging from global patterns (e.g., latitudinal gradients in species richness/diversity and degree of herbivore feeding specialization) to localized insect–plant associations that span the biospectrum from polyphenisms, polymorphisms, biotypes, demes, host races, to cryptic species, remain academically contentious. Semantic and biosystematic (taxonomical) disagreements sometimes detract from more important ecological and evolutionary processes that drive diversification, the dynamics of gene flow and local extinctions. This review addresses several aspects of insect specialization, host‐associated divergence and ecological (including “hybrid”) speciation, with special reference to the climate warming impacts on species borders of hybridizing swallowtail butterflies (Papilionidae). Interspecific hybrid introgression may result in collapse of multi‐species communities or increase species numbers via homoploid hybrid speciation. We may see diverging, merging, or emerging genotypes across hybrid zones, all part of the ongoing processes of evolution. Molecular analyses of genetic mosaics and genomic dynamics with “divergence hitchhiking”, combined with ecological, ethological and physiological studies of “species porosity”, have already begun to unveil some answers for some important ecological/evolutionary questions. (i) How rapidly can host‐associated divergence lead to new species (and why doesn't it always do so, e.g., resulting in “incomplete” speciation)? (ii) How might “speciation genes” function, and how/where would we find them? (iii) Can oscillations from specialists to generalists and back to specialists help explain global diversity in herbivorous insects? (iv) How could recombinant interspecific hybridization lead to divergence and speciation? From ancient phytochemically defined angiosperm affiliations to recent and very local geographical mosaics, the Papilionidae (swallowtail butterflies) have provided a model for enhanced understanding of ecological patterns and evolutionary processes, including host‐associated genetic divergence, genomic mosaics, genetic hitchhiking and sex‐linked speciation genes. Apparent homoploid hybrid speciation in Papilio appears to have been catalyzed by climate warming‐induced interspecific introgression of some, but not all, species diagnostic traits, reflecting strong divergent selection (discordant), especially on the Z (= X) chromosome. Reproductive isolation of these novel recombinant hybrid genotypes appears to be accomplished via a delayed post‐diapause emergence or temporal isolation, and is perhaps aided by the thermal landscape. Changing thermal landscapes appear to have created (and may destroy) novel recombinant hybrid genotypes and hybrid species.     

Ecological divergence in the presence of gene flow in two closely related Oryza species (Oryza rufipogon and O. nivara)

Ecological divergence plays a prominent role in the process of speciation, but how divergence occurs in the face of gene flow is still less clear, and remains controversial among evolutionists. Here we investigated the nucleotide diversity, divergence and gene flow between Oryza nivara and O. rufipogon using sequences of seven chloroplast and nuclear loci. By analysing samples from 26 wild populations across the geographic ranges of the two species, we showed that both species were highly structured and O. rufipogon maintained a higher level of species‐wide diversity than O. nivara. Notably, phylogenetic, amova and F_ST analyses were unable to detect significant nucleotide differentiation between the two species. We estimated that the two species began to diverge at c. 0.16 million years ago. Our coalescent‐based simulations strongly rejected the simple isolation model of zero migration between species, but rather provided unambiguous evidence of bidirectional gene flow between species, particularly from O. rufipogon to O. nivara. Our simulations also indicated that gene flow was recurrent during the divergence process rather than arising from secondary contact after allopatric divergence. In conjunction with different morphological and life‐history traits and habitat preference in the two species, this study supports the hypothesis that these Oryza species are better treated as ecotypes that diverged quite recently and are still under the process of divergence. Importantly, we demonstrate the ecological divergence between O. rufipogon and O. nivara in the presence of significant gene flow, implying that natural selection plays a primary role in driving the divergence of the two Oryza species.     

Evolutionary history of the greater white-toothed shrew (Crocidura russula) inferred from analysis of mtDNA, Y, and X chromosome markers

We investigate the evolutionary history of the greater white-toothed shrew across its distribution in northern Africa and mainland Europe using sex-specific (mtDNA and Y chromosome) and biparental (X chromosome) markers. All three loci confirm a large divergence between eastern (Tunisia and Sardinia) and western (Morocco and mainland Europe) lineages, and application of a molecular clock to mtDNA divergence estimates indicates a more ancient separation (2.25 M yr ago) than described by some previous studies, supporting claims for taxonomic revision. Moroccan ancestry for the mainland European population is inconclusive from phylogenetic trees, but is supported by greater nucleotide diversity and a more ancient population expansion in Morocco than in Europe. Signatures of rapid population expansion in mtDNA, combined with low X and Y chromosome diversity, suggest a single colonization of mainland Europe by a small number of Moroccan shrews >38 K yr ago. This study illustrates that multilocus genetic analyses can facilitate the interpretation of species' evolutionary history but that phylogeographic inference using X and Y chromosomes is restricted by low levels of observed polymorphism.     

Glacial cycles as an allopatric speciation pump in north‐eastern American freshwater fishes

Allopatric speciation may be the principal mechanism generating new species. Yet, it remains difficult to judge the generality of this process because few studies have provided evidence that geographic isolation has triggered the development of reproductive isolation over multiple species of a regional fauna. Here, we first combine results from new empirical data sets (7 taxa) and published literature (9 taxa) to show that the eastern Great Lakes drainage represents a multispecies suture zone for glacial lineages of freshwater fishes with variable levels of genetic divergence. Second, we performed amplified fragment length polymorphism analyses among four pairs of lineages. Results indicate that lineages with relatively deep levels of mtDNA 5′ COI (barcode) sequence divergence (>2%) developed strong reproductive barriers, while lineages with lower levels of divergence show weaker reproductive isolation when found in sympatry. This suggests that a threshold of 2% sequence divergence at mtDNA could be used as a first step to flag cryptic species in North American freshwater fishes. By describing different levels of divergence and reproductive isolation in different co‐occurring fishes, we offer strong evidence that allopatric speciation has contributed significantly to the diversification of north‐eastern American freshwater fishes and confirm that Pleistocene glacial cycles can be viewed as a ‘speciation pump’ that played a predominant role in generating biodiversity.     

Vascular plant diversity in eastern Asia and North America: historical and ecological explanations

Taxonomic diversity of vascular plants (ferns, gymnosperms and angiosperms) was compared between eastern Asia and North America. Eastern Asia has significantly higher species richness in all three classes but the difference was greatest in ferns and least in angiosperms. Differences in taxonomic treatments between the two continents are not likely contributors to these patterns. The relationship of regional to global species richness across the three plant classes suggested that diversity patterns were relatively homogeneous at three taxonomic levels. Thus, differences in species richness are established at the family level and are therefore relatively old. The previously noted fact that eastern Asia has a higher proportion of primitive taxa was shown by analyses both among and within plant classes. Diversity patterns across three taxonomic levels (i.e. family, genus and species) of the three classes may reflect the relative historical positions of the two continents (following continental drift) to the centre(s) of their origin, neighbouring land masses, differential speciation/extinction rates, and switches in dominance levels associated with climate change (including glaciation), as well as reproductive/dispersal mechanisms of the three plant classes.     

Genetic association of physically unlinked islands of genomic divergence in incipient species of Anopheles gambiae

Previous efforts to uncover the genetic underpinnings of ongoing ecological speciation of the M and S forms of the African malaria vector Anopheles gambiae revealed two centromere‐proximal islands of genetic divergence on X and chromosome 2. Under the assumption of considerable ongoing gene flow between M and S, these persistently divergent genomic islands were widely considered to be ‘speciation islands’. In the course of microarray‐based divergence mapping, we discovered a third centromere‐associated island of divergence on chromosome 3, which was validated by targeted re‐sequencing. To test for genetic association between the divergence islands on all three chromosomes, SNP‐based assays were applied in four natural populations of M and S spanning West, Central and East Africa. Genotyping of 517 female M and S mosquitoes revealed nearly complete linkage disequilibrium between the centromeres of the three independently assorting chromosomes. These results suggest that despite the potential for inter‐form gene flow through hybridization, actual (realized) gene flow between M and S may be substantially less than commonly assumed and may not explain most shared variation. Moreover, the possibility of very low gene flow calls into question whether diverged pericentromeric regions—characterized by reduced levels of variation and recombination—are in fact instrumental rather than merely incidental to the speciation process.     

Influence of gene flow on divergence dating – implications for the speciation history of Takydromus grass lizards

Dating the time of divergence and understanding speciation processes are central to the study of the evolutionary history of organisms but are notoriously difficult. The difficulty is largely rooted in variations in the ancestral population size or in the genealogy variation across loci. To depict the speciation processes and divergence histories of three monophyletic Takydromus species endemic to Taiwan, we sequenced 20 nuclear loci and combined with one mitochondrial locus published in GenBank. They were analysed by a multispecies coalescent approach within a Bayesian framework. Divergence dating based on the gene tree approach showed high variation among loci, and the divergence was estimated at an earlier date than when derived by the species‐tree approach. To test whether variations in the ancestral population size accounted for the majority of this variation, we conducted computer inferences using isolation‐with‐migration (IM) and approximate Bayesian computation (ABC) frameworks. The results revealed that gene flow during the early stage of speciation was strongly favoured over the isolation model, and the initiation of the speciation process was far earlier than the dates estimated by gene‐ and species‐based divergence dating. Due to their limited dispersal ability, it is suggested that geographical isolation may have played a major role in the divergence of these Takydromus species. Nevertheless, this study reveals a more complex situation and demonstrates that gene flow during the speciation process cannot be overlooked and may have a great impact on divergence dating. By using multilocus data and incorporating Bayesian coalescence approaches, we provide a more biologically realistic framework for delineating the divergence history of Takydromus.     

Speciation as a sieve for ancestral polymorphism

Because they are considered rare, balanced polymorphisms are often discounted as crucial constituents of genome‐wide variation in sequence diversity. Despite its perceived rarity, however, long‐term balancing selection can elevate genetic diversity and significantly affect observed divergence between species. Here, we discuss how ancestral balanced polymorphisms can be “sieved” by the speciation process, which sorts them unequally across descendant lineages. After speciation, ancestral balancing selection is revealed by genomic regions of high divergence between species. This signature, which resembles that of other evolutionary processes, can potentially confound genomic studies of population divergence and inferences of “islands of speciation.”     

Robbery in progress: Historical museum collections bring to light a mitochondrial capture within a bird species widespread across southern Australia,the Copperback Quail‐thrush Cinclosoma clarum

We surveyed mitochondrial, autosomal, and Z chromosome diversity within and between the Copperback Quail‐thrush Cinclosoma clarum and Chestnut Quail‐thrush C. castanotum, which together span the arid and semi‐arid zones of southern Australia, and primarily from specimens held in museum collections. We affirm the recent taxonomic separation of the two species and then focus on diversity within the more widespread of the two species, C. clarum. To guide further study of the system and what it offers to understanding the genomics of the differentiation and speciation processes, we develop and present a hypothesis to explain mitonuclear discordance that emerged in ourdata. Following a period of historical allopatry, secondary contact has resulted in an eastern mitochondrial genome replacing the western mitochondrial genome in western populations. This is predicted under a population‐level invasion in the opposite direction, that of the western population invading the range of the eastern one. Mitochondrial captures can be driven by neutral, demographic processes, or adaptive mechanisms, and we favor the hypothesized capture being driven by neutral means. We cannot fully reject the adaptive process but suggest how these alternatives may be further tested. We acknowledge an alternative hypothesis, which finds some support in phenotypic data published elsewhere, namely that outcomes of secondary contact have been more complex than our current genomic data suggest. Discriminating and reconciling these two alternative hypotheses, which may not be mutually exclusive, could be tested with closer sampling at levels of population, individual, and nucleotide than has so far been possible. This would be further aided by knowledge of the genetic basis to phenotypic variation described elsewhere.     

Inter‐specific gene flow dynamics during the Pleistocene‐dated speciation of forest‐dependent mosquitoes in Southeast Asia

Tropical forests have undergone repeated fragmentation and expansion during Pleistocene glacial and interglacial periods, respectively. The effects of this repeated forest fragmentation in driving vicariance in tropical taxa have been well studied. However, relatively little is known about how often this process results in allopatric speciation, since it may be inhibited by recurrent gene flow during repeated secondary contact, or to what extent Pleistocene‐dated speciation results from ecological specialization in the face of gene flow. Here, divergence times and gene flow between three closely‐related mosquito species of the Anopheles dirus species complex endemic to the forests of Southeast Asia, are inferred using coalescent based Bayesian analysis. An Isolation with Migration model is applied to sequences of two mitochondrial and three nuclear genes, and 11 microsatellites. The divergence of An. scanloni has occurred despite unidirectional nuclear gene flow from this species into An. dirus. The inferred asymmetric gene flow may result from the unique evolutionary adaptation of An. scanloni to limestone karst habitat, and therefore the fitness advantage of this species over An. dirus in regions of sympatry. Mitochondrial introgression has led to the complete replacement of An. dirus haplotypes with those of An. baimaii through a recent (～62 kya) selective sweep. Speciation of An. baimaii and An. dirus is inferred to have involved allopatric divergence throughout much of the Pleistocene. Secondary contact and bidirectional gene flow has occurred only within the last 100 000 years, by which time the process of allopatric speciation seems to have been largely completed.     

Genetic variability,local selection and demographic history: genomic evidence of evolving towards allopatric speciation in Asian seabass

Genomewide analysis of genetic divergence is critically important in understanding the genetic processes of allopatric speciation. We sequenced RAD tags of 131 Asian seabass individuals of six populations from South‐East Asia and Australia/Papua New Guinea. Using 32 433 SNPs, we examined the genetic diversity and patterns of population differentiation across all the populations. We found significant evidence of genetic heterogeneity between South‐East Asian and Australian/Papua New Guinean populations. The Australian/Papua New Guinean populations showed a rather lower level of genetic diversity. F_ST and principal components analysis revealed striking divergence between South‐East Asian and Australian/Papua New Guinean populations. Interestingly, no evidence of contemporary gene flow was observed. The demographic history was further tested based on the folded joint site frequency spectrum. The scenario of ancient migration with historical population size changes was suggested to be the best fit model to explain the genetic divergence of Asian seabass between South‐East Asia and Australia/Papua New Guinea. This scenario also revealed that Australian/Papua New Guinean populations were founded by ancestors from South‐East Asia during mid‐Pleistocene and were completely isolated from the ancestral population after the last glacial retreat. We also detected footprints of local selection, which might be related to differential ecological adaptation. The ancient gene flow was examined and deemed likely insufficient to counteract the genetic differentiation caused by genetic drift. The observed genomic pattern of divergence conflicted with the ‘genomic islands’ scenario. Altogether, Asian seabass have likely been evolving towards allopatric speciation since the split from the ancestral population during mid‐Pleistocene.