The role of recombination,niche‐specific gene pools and flexible genomes in the ecological speciation of bacteria

Bacteria diversify into genetic clusters analogous to those observed in sexual eukaryotes, but the definition of bacterial species is an ongoing problem. Recent work has focused on adaptation to distinct ecological niches as the main driver of clustering, but there remains debate about the role of recombination in that process. One view is that homologous recombination occurs too rarely for gene flow to constrain divergent selection. Another view is that homologous recombination is frequent enough in many bacterial populations that barriers to gene flow are needed to permit divergence. Niche‐specific gene pools have been proposed as a general mechanism to limit gene flow. We use theoretical models to evaluate additional hypotheses that evolving genetic architecture, specifically the effect sizes of genes and gene gain and loss, can limit gene flow between diverging populations. Our model predicts that (a) in the presence of gene flow and recombination, ecological divergence is concentrated in few loci of large effect and (b) high rates of gene flow plus recombination promote gene loss and favor the evolution of niche‐specific genes. The results show that changing genetic architecture and gene loss can facilitate ecological divergence, even without niche‐specific gene pools. We discuss these results in the context of recent studies of sympatric divergence in microbes.     

Localizing FST outliers on a QTL map reveals evidence for large genomic regions of reduced gene exchange during speciation‐with‐gene‐flow

Populations that maintain phenotypic divergence in sympatry typically show a mosaic pattern of genomic divergence, requiring a corresponding mosaic of genomic isolation (reduced gene flow). However, mechanisms that could produce the genomic isolation required for divergence‐with‐gene‐flow have barely been explored, apart from the traditional localized effects of selection and reduced recombination near centromeres or inversions. By localizing F_ST outliers from a genome scan of wild pea aphid host races on a Quantitative Trait Locus (QTL) map of key traits, we test the hypothesis that between‐population recombination and gene exchange are reduced over large ‘divergence hitchhiking’ (DH) regions. As expected under divergence hitchhiking, our map confirms that QTL and divergent markers cluster together in multiple large genomic regions. Under divergence hitchhiking, the nonoutlier markers within these regions should show signs of reduced gene exchange relative to nonoutlier markers in genomic regions where ongoing gene flow is expected. We use this predicted difference among nonoutliers to perform a critical test of divergence hitchhiking. Results show that nonoutlier markers within clusters of F_ST outliers and QTL resolve the genetic population structure of the two host races nearly as well as the outliers themselves, while nonoutliers outside DH regions reveal no population structure, as expected if they experience more gene flow. These results provide clear evidence for divergence hitchhiking, a mechanism that may dramatically facilitate the process of speciation‐with‐gene‐flow. They also show the power of integrating genome scans with genetic analyses of the phenotypic traits involved in local adaptation and population divergence.     

Reanalysis suggests that genomic islands of speciation are due to reduced diversity,not reduced gene flow

The metaphor of ‘genomic islands of speciation’ was first used to describe heterogeneous differentiation among loci between the genomes of closely related species. The biological model proposed to explain these differences was that the regions showing high levels of differentiation were resistant to gene flow between species, while the remainder of the genome was being homogenized by gene flow and consequently showed lower levels of differentiation. However, the conditions under which such differentiation can occur at multiple unlinked loci are restrictive; additionally, essentially, all previous analyses have been carried out using relative measures of divergence, which can be misleading when regions with different levels of recombination are compared. Here, we test the model of differential gene flow by asking whether absolute divergence is also higher in the previously identified ‘islands’. Using five species pairs for which full sequence data are available, we find that absolute measures of divergence are not higher in genomic islands. Instead, in all cases examined, we find reduced diversity in these regions, a consequence of which is that relative measures of divergence are abnormally high. These data therefore do not support a model of differential gene flow among loci, although islands of relative divergence may represent loci involved in local adaptation. Simulations using the program IMa2 further suggest that inferences of any gene flow may be incorrect in many comparisons. We instead present an alternative explanation for heterogeneous patterns of differentiation, one in which postspeciation selection generates patterns consistent with multiple aspects of the data.     

Hitching a lift on the road to speciation

Understanding how speciation can take place in the presence of homogenizing gene flow remains a major challenge in evolutionary biology. In the early stages of ecological speciation, reproductive isolation between populations occupying different habitats is expected to be concentrated around genes for local adaptation. These genomic regions will show high divergence while gene exchange in other regions of the genome should continue relatively unimpaired, resulting in low levels of differentiation. The problem is to explain how speciation progresses from this point towards complete reproductive isolation, allowing genome‐wide divergence. A new study by Via and West (2008) on speciation between host races of the pea aphid, Acyrthosiphon pisum, introduces the mechanism of ‘divergence hitchhiking’ which can generate large ‘islands of differentiation’ and facilitate the build‐up of linkage disequilibrium, favouring increased reproductive isolation. This idea potentially removes a major stumbling block to speciation under continuous gene flow.     

THE GENETIC ARCHITECTURE OF REPRODUCTIVE ISOLATION DURING SPECIATION‐WITH‐GENE‐FLOW IN LAKE WHITEFISH SPECIES PAIRS ASSESSED BY RAD SEQUENCING

During speciation‐with‐gene‐flow, effective migration varies across the genome as a function of several factors, including proximity of selected loci, recombination rate, strength of selection, and number of selected loci. Genome scans may provide better empirical understanding of the genome‐wide patterns of genetic differentiation, especially if the variance due to the previously mentioned factors is partitioned. In North American lake whitefish (Coregonus clupeaformis), glacial lineages that diverged in allopatry about 60,000 years ago and came into contact 12,000 years ago have independently evolved in several lakes into two sympatric species pairs (a normal benthic and a dwarf limnetic). Variable degrees of reproductive isolation between species pairs across lakes offer a continuum of genetic and phenotypic divergence associated with adaptation to distinct ecological niches. To disentangle the complex array of genetically based barriers that locally reduce the effective migration rate between whitefish species pairs, we compared genome‐wide patterns of divergence across five lakes distributed along this divergence continuum. Using restriction site associated DNA (RAD) sequencing, we combined genetic mapping and population genetics approaches to identify genomic regions resistant to introgression and derive empirical measures of the barrier strength as a function of recombination distance. We found that the size of the genomic islands of differentiation was influenced by the joint effects of linkage disequilibrium maintained by selection on many loci, the strength of ecological niche divergence, as well as demographic characteristics unique to each lake. Partial parallelism in divergent genomic regions likely reflected the combined effects of polygenic adaptation from standing variation and independent changes in the genetic architecture of postzygotic isolation. This study illustrates how integrating genetic mapping and population genomics of multiple sympatric species pairs provide a window on the speciation‐with‐gene‐flow mechanism.     

Spatially explicit models of divergence and genome hitchhiking

Strong barriers to genetic exchange can exist at divergently selected loci, whereas alleles at neutral loci flow more readily between populations, thus impeding divergence and speciation in the face of gene flow. However, ‘divergence hitchhiking’ theory posits that divergent selection can generate large regions of differentiation around selected loci. ‘Genome hitchhiking’ theory suggests that selection can also cause reductions in average genome‐wide rates of gene flow, resulting in widespread genomic divergence (rather than divergence only around specific selected loci). Spatial heterogeneity is ubiquitous in nature, yet previous models of genetic barriers to gene flow have explored limited combinations of spatial and selective scenarios. Using simulations of secondary contact of populations, we explore barriers to gene flow in various selective and spatial contexts in continuous, two‐dimensional, spatially explicit environments. In general, the effects of hitchhiking are strongest in environments with regular spatial patterning of starkly divergent habitat types. When divergent selection is very strong, the absence of intermediate habitat types increases the effects of hitchhiking. However, when selection is moderate or weak, regular (vs. random) spatial arrangement of habitat types becomes more important than the presence of intermediate habitats per se. We also document counterintuitive processes arising from the stochastic interplay between selection, gene flow and drift. Our results indicate that generalization of results from two‐deme models requires caution and increase understanding of the genomic and geographic basis of population divergence.     

Migration-induced phenotypic divergence: the migration-selection balance of correlated traits

Genetically correlated traits are known to respond to indirect selection pressures caused by directional selection on other traits. It is however unclear how local adaptation in populations diverging along some phenotypic traits but not others is affected by the joint action of gene flow and genetic correlations among traits. This simulation study shows that although gene flow is a potent constraining mechanism of population adaptive divergence, it may induce phenotypic divergence in traits under homogeneous selection among habitats if they are genetically correlated with traits under divergent selection. This correlated phenotypic divergence is a nonmonotonous function of migration and increases with mutational correlation among traits. It also increases with the number of divergently selected traits provided their genetic autonomy relative to the uniformly selected trait is reduced by specific patterns of genetic covariances: populations with lower effective trait dimensionality are more likely to generate very large correlated divergence. The correlated divergence is likely to be picked up by Q(ST)-F(ST) analysis of population genetic differentiation and be erroneously ascribed to adaptive divergence under divergent selection. This study emphasizes the necessity to understand the interaction between selection and the genetic basis of adaptation in a multivariate rather than univariate context.     

A genomic island linked to ecotype divergence in Atlantic cod

The genomic architecture underlying ecological divergence and ecological speciation with gene flow is still largely unknown for most organisms. One central question is whether divergence is genome‐wide or localized in ‘genomic mosaics’ during early stages when gene flow is still pronounced. Empirical work has so far been limited, and the relative impacts of gene flow and natural selection on genomic patterns have not been fully explored. Here, we use ecotypes of Atlantic cod to investigate genomic patterns of diversity and population differentiation in a natural system characterized by high gene flow and large effective population sizes, properties which theoretically could restrict divergence in local genomic regions. We identify a genomic region of strong population differentiation, extending over approximately 20 cM, between pairs of migratory and stationary ecotypes examined at two different localities. Furthermore, the region is characterized by markedly reduced levels of genetic diversity in migratory ecotype samples. The results highlight the genomic region, or ‘genomic island’, as potentially associated with ecological divergence and suggest the involvement of a selective sweep. Finally, we also confirm earlier findings of localized genomic differentiation in three other linkage groups associated with divergence among eastern Atlantic populations. Thus, although the underlying mechanisms are still unknown, the results suggest that ‘genomic mosaics’ of differentiation may even be found under high levels of gene flow and that marine fishes may provide insightful model systems for studying and identifying initial targets of selection during ecological divergence.     

Genomic islands of divergence linked to ecotypic variation in sockeye salmon

Regions of the genome displaying elevated differentiation (genomic islands of divergence) are thought to play an important role in local adaptation, especially in populations experiencing high gene flow. However, the characteristics of these islands as well as the functional significance of genes located within them remain largely unknown. Here, we used data from thousands of SNPs aligned to a linkage map to investigate genomic islands of divergence in three ecotypes of sockeye salmon (Oncorhynchus nerka) from a single drainage in southwestern Alaska. We found ten islands displaying high differentiation among ecotypes. Conversely, neutral structure observed throughout the rest of the genome was low and not partitioned by ecotype. One island on linkage group So13 was particularly large and contained six SNPs with F_ST > 0.14 (average F_ST of neutral SNPs = 0.01). Functional annotation revealed that the peak of this island contained a nonsynonymous mutation in a gene involved in growth in other species (TULP4). The islands that we discovered were relatively small (80–402 Kb), loci found in islands did not show reduced levels of diversity, and loci in islands displayed slightly elevated linkage disequilibrium. These attributes suggest that the islands discovered here were likely generated by divergence hitchhiking; however, we cannot rule out the possibility that other mechanisms may have produced them. Our results suggest that islands of divergence serve an important role in local adaptation with gene flow and represent a significant advance towards understanding the genetic basis of ecotypic differentiation.     

Divergent selection and heterogeneous genomic divergence

Levels of genetic differentiation between populations can be highly variable across the genome, with divergent selection contributing to such heterogeneous genomic divergence. For example, loci under divergent selection and those tightly physically linked to them may exhibit stronger differentiation than neutral regions with weak or no linkage to such loci. Divergent selection can also increase genome‐wide neutral differentiation by reducing gene flow (e.g. by causing ecological speciation), thus promoting divergence via the stochastic effects of genetic drift. These consequences of divergent selection are being reported in recently accumulating studies that identify: (i) ‘outlier loci’ with higher levels of divergence than expected under neutrality, and (ii) a positive association between the degree of adaptive phenotypic divergence and levels of molecular genetic differentiation across population pairs [‘isolation by adaptation’ (IBA)]. The latter pattern arises because as adaptive divergence increases, gene flow is reduced (thereby promoting drift) and genetic hitchhiking increased. Here, we review and integrate these previously disconnected concepts and literatures. We find that studies generally report 5–10% of loci to be outliers. These selected regions were often dispersed across the genome, commonly exhibited replicated divergence across different population pairs, and could sometimes be associated with specific ecological variables. IBA was not infrequently observed, even at neutral loci putatively unlinked to those under divergent selection. Overall, we conclude that divergent selection makes diverse contributions to heterogeneous genomic divergence. Nonetheless, the number, size, and distribution of genomic regions affected by selection varied substantially among studies, leading us to discuss the potential role of divergent selection in the growth of regions of differentiation (i.e. genomic islands of divergence), a topic in need of future investigation.     

Extensive linkage disequilibrium and parallel adaptive divergence across threespine stickleback genomes

Population genomic studies are beginning to provide a more comprehensive view of dynamic genome-scale processes in evolution. Patterns of genomic architecture, such as genomic islands of increased divergence, may be important for adaptive population differentiation and speciation. We used next-generation sequencing data to examine the patterns of local and long-distance linkage disequilibrium (LD) across oceanic and freshwater populations of threespine stickleback, a useful model for studies of evolution and speciation. We looked for associations between LD and signatures of divergent selection, and assessed the role of recombination rate variation in generating LD patterns. As predicted under the traditional biogeographic model of unidirectional gene flow from ancestral oceanic to derived freshwater stickleback populations, we found extensive local and long-distance LD in fresh water. Surprisingly, oceanic populations showed similar patterns of elevated LD, notably between large genomic regions previously implicated in adaptation to fresh water. These results support an alternative biogeographic model for the stickleback radiation, one of a metapopulation with appreciable bi-directional gene flow combined with strong divergent selection between oceanic and freshwater populations. As predicted by theory, these processes can maintain LD within and among genomic islands of divergence. These findings suggest that the genomic architecture in oceanic stickleback populations may provide a mechanism for the rapid re-assembly and evolution of multi-locus genotypes in newly colonized freshwater habitats, and may help explain genetic mapping of parallel phenotypic variation to similar loci across independent freshwater populations.     

Gene flow and selection interact to promote adaptive divergence in regions of low recombination

Adaptation to new environments often occurs in the face of gene flow. Under these conditions, gene flow and recombination can impede adaptation by breaking down linkage disequilibrium between locally adapted alleles. Theory predicts that this decay can be halted or slowed if adaptive alleles are tightly linked in regions of low recombination, potentially favouring divergence and adaptive evolution in these regions over others. Here, we compiled a global genomic data set of over 1,300 individual threespine stickleback from 52 populations and compared the tendency for adaptive alleles to occur in regions of low recombination between populations that diverged with or without gene flow. In support of theory, we found that putatively adaptive alleles (F_ST and d_XY outliers) tend to occur more often in regions of low recombination in populations where divergent selection and gene flow have jointly occurred. This result remained significant when we employed different genomic window sizes, controlled for the effects of mutation rate and gene density, controlled for overall genetic differentiation, varied the genetic map used to estimate recombination and used a continuous (rather than discrete) measure of geographic distance as proxy for gene flow/shared ancestry. We argue that our study provides the first statistical evidence that the interaction of gene flow and selection biases divergence toward regions of low recombination.     

Genomics of adaptive divergence with chromosome‐scale heterogeneity in crossover rate

Genetic differentiation between divergent populations is often greater in chromosome centres than peripheries. Commonly overlooked, this broadscale differentiation pattern is sometimes ascribed to heterogeneity in crossover rate and hence linked selection within chromosomes, but the underlying mechanisms remain incompletely understood. A literature survey across 46 organisms reveals that most eukaryotes indeed exhibit a reduced crossover rate in chromosome centres relative to the peripheries. Using simulations of populations diverging into ecologically different habitats through sorting of standing genetic variation, we demonstrate that such chromosome‐scale heterogeneity in crossover rate, combined with polygenic divergent selection, causes stronger hitchhiking and especially barriers to gene flow across chromosome centres. Without requiring selection on new mutations, this rapidly leads to elevated population differentiation in the low‐crossover centres relative to the high‐crossover peripheries of chromosomes (“Chromosome Centre‐Biased Differentiation”, CCBD). Using simulated and empirical data, we then show that strong CCBD between populations can provide evidence of polygenic adaptive divergence with a phase of gene flow. We further demonstrate that chromosome‐scale heterogeneity in crossover rate impacts analyses beyond that of population differentiation, including the inference of phylogenies and parallel adaptive evolution among populations, the detection of genetic loci under selection, and the interpretation of the strength of selection on genomic regions. Overall, our results call for a greater appreciation of chromosome‐scale heterogeneity in crossover rate in evolutionary genomics.     

Characterization of a genomic divergence island between black-and-yellow and gopher Sebastes rockfishes

Islands of high genomic divergence that contain genes of evolutionary significance may form between diverging species. The gopher rockfish, Sebastes carnatus, and black-and-yellow rockfish, S. chrysomelas, are sympatrically distributed temperate marine species inhabiting rocky reefs and kelp forests on the west coast of the United States. Prior studies documented low levels of genetic divergence between the two species, except at a single microsatellite locus that displayed high divergence, Sra.7-2. To better characterize genome wide divergence, we scored 25 additional microsatellite loci. Mean neutral divergence between species (F(ST) = 0.01) changed little from prior estimates. Sra.7-2 continued to be an extreme divergence outlier. Five novel microsatellites within ± 15 kb of Sra.7-2 were characterized. High divergence, consistently low diversity in S. chrysomelas, and linkage disequilibrium were detected at these loci, suggesting the influence of recent selection. However, coalescent modelling of divergence at neutral and Sra.7-2 regions showed that initial divergence at Sra.7-2 was ancient, likely predating divergence at neutral regions. It is therefore unlikely that Sra.7-2 divergence represents solely recent ecological divergence within one species and may represent the action of recurrent selection. Introgressive gene flow (2N(E) m) was much higher (>1) at neutral than Sra.7-2 regions (<1) despite evidence that two S. carnatus individuals have recently mixed ancestry at the Sra.7-2 region. The Sra.7-2 genomic region is likely one of several regions containing genes involved in initiating and maintaining species integrity. Completion of the final stages of speciation appears to be a slow and ongoing process for these species.     

Heterogeneity and concordance in locus‐specific differentiation and introgression between species of towhees

The maintenance or breakdown of reproductive isolation is an observable outcome of secondary contact between species. In cases where hybrids beyond the F1 are formed, the representation of each species' ancestry can vary dramatically among genomic regions. This genomic heterogeneity in ancestry and introgression can offer insight into evolutionary processes, particularly if introgression is compared in multiple hybrid zones. Similarly, considerable heterogeneity exists across the genome in the extent to which populations and species have diverged, reflecting the combined effects of different evolutionary processes on genetic variation. We studied hybridization across two hybrid zones of two phenotypically well‐differentiated bird species in Mexico (Pipilo maculatus and P. ocai), to investigate genomic heterogeneity in differentiation and introgression. Using genotyping‐by‐sequencing (GBS) and hierarchical Bayesian models, we genotyped 460 birds at over 41 000 single nucleotide polymorphism (SNP) loci. We identified loci exhibiting extreme introgression relative to the genome‐wide expectation using a Bayesian genomic cline model. We also estimated locus‐specific F_ST and identified loci with exceptionally high genetic divergence between the parental species. We found some concordance of locus‐specific introgression in the two independent hybrid zones (6–20% of extreme loci shared across zones), reflecting areas of the genome that experience similar gene flow when the species interact. Additionally, heterogeneity in introgression and divergence across the genome revealed another subset of loci under the influence of locally specific factors. These results are consistent with a history in which reproductive isolation has been influenced by a common set of loci in both hybrid zones, but where local environmental and stochastic factors also lead to genomic differentiation.     

Disentangling structural genomic and behavioural barriers in a sea of connectivity

Genetic divergence among populations arises through natural selection or drift and is counteracted by connectivity and gene flow. In sympatric populations, isolating mechanisms are thus needed to limit the homogenizing effects of gene flow to allow for adaptation and speciation. Chromosomal inversions act as an important mechanism maintaining isolating barriers, yet their role in sympatric populations and divergence with gene flow is not entirely understood. Here, we revisit the question of whether inversions play a role in the divergence of connected populations of the marine fish Atlantic cod (Gadus morhua), by exploring a unique data set combining whole‐genome sequencing data and behavioural data obtained with acoustic telemetry. Within a confined fjord environment, we find three genetically differentiated Atlantic cod types belonging to the oceanic North Sea population, the western Baltic population and a local fjord‐type cod. Continuous behavioural tracking over 4 year revealed temporally stable sympatry of these types within the fjord. Despite overall weak genetic differentiation consistent with high levels of gene flow, we detected significant frequency shifts of three previously identified inversions, indicating an adaptive barrier to gene flow. In addition, behavioural data indicated that North Sea cod and individuals homozygous for the LG12 inversion had lower fitness in the fjord environment. However, North Sea and fjord‐type cod also occupy different depths, possibly contributing to prezygotic reproductive isolation and representing a behavioural barrier to gene flow. Our results provide the first insights into a complex interplay of genomic and behavioural isolating barriers in Atlantic cod and establish a new model system towards an understanding of the role of genomic structural variants in adaptation and diversification.     

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.     

Genomic divergence across ecological gradients in the Central African rainforest songbird (Andropadus virens)

The little greenbul, a common rainforest passerine from sub‐Saharan Africa, has been the subject of long‐term evolutionary studies to understand the mechanisms leading to rainforest speciation. Previous research found morphological and behavioural divergence across rainforest–savannah transition zones (ecotones), and a pattern of divergence with gene flow suggesting divergent natural selection has contributed to adaptive divergence and ecotones could be important areas for rainforests speciation. Recent advances in genomics and environmental modelling make it possible to examine patterns of genetic divergence in a more comprehensive fashion. To assess the extent to which natural selection may drive patterns of differentiation, here we investigate patterns of genomic differentiation among populations across environmental gradients and regions. We find compelling evidence that individuals form discrete genetic clusters corresponding to distinctive environmental characteristics and habitat types. Pairwise F_ST between populations in different habitats is significantly higher than within habitats, and this differentiation is greater than what is expected from geographic distance alone. Moreover, we identified 140 SNPs that showed extreme differentiation among populations through a genomewide selection scan. These outliers were significantly enriched in exonic and coding regions, suggesting their functional importance. Environmental association analysis of SNP variation indicates that several environmental variables, including temperature and elevation, play important roles in driving the pattern of genomic diversification. Results lend important new genomic evidence for environmental gradients being important in population differentiation.     

Genomewide patterns of variation in genetic diversity are shared among populations,species and higher‐order taxa

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (N_e) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage‐specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in N_e. Utilizing population‐level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in N_e would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population‐scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (F_ST, PBS, d_xy) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.     

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.