Population genomics and speciation

How Does the Application of Beauveria bassiana and Compost on Corn Crops Affect the Survival and Genetic Diversity of Phyllophaga obsoleta (Coleoptera: Melolonthinae)?

The knowledge of the diversity and genetic structure of pest insects under management contributes to the improvement of control strategies. An experiment was run to investigate whether the addition of the fungus Beauveria bassiana (Balsamo) Vuillemin (Hypocreales: Cordycipitaceae) (BB) and compost to soil affects the presence and genetic diversity of adults and larvae of Phyllophaga obsoleta Blanch (Coleoptera: Melolonthinae) larvae in maize crops. We collected adults in and used mating pairs under four treatments (BB, compost, soil, blank). Genetic diversity and structure were determined through five allo/iso-enzymatic loci. Beauveria bassiana affected the presence and mortality of P. obsoleta in the laboratory but not under field conditions. The genetic diversity of P. obsoleta ranged from moderate to high (Ho = 0.26-0.31), with a low genetic differentiation among localities or treatments (Phi < 0.05), indicating high levels of gene flow. Our results showed a weak effect of B. bassiana on P. obsoleta in the field. Still, our laboratory observations suggest that the fungus may be a suitable alternative for biological control.

Comparative studies on speciation: 30 years since Coyne and Orr

Understanding the processes of population divergence and speciation remains a core question in evolutionary biology. For nearly a hundred years evolutionary geneticists have characterized reproductive isolation (RI) mechanisms and specific barriers to gene flow required for species formation. The seminal work of Coyne and Orr provided the first comprehensive comparative analysis of speciation. By combining phylogenetic hypotheses and species range data with estimates of genetic divergence and multiple mechanisms of RI across Drosophila, Coyne and Orr's influential meta-analyses answered fundamental questions and motivated new analyses that continue to push the field forward today. Now 30 years later, we revisit the five questions addressed by Coyne and Orr, identifying results that remain well supported and others that seem less robust with new data. We then consider the future of speciation research, with emphasis on areas where novel methods and data motivate potential progress. While the literature remains biased towards Drosophila and other model systems, we are enthusiastic about the future of the field.

The role of copy-number variation in the reinforcement of sexual isolation between the two European subspecies of the house mouse

Reinforcement has the potential to generate strong reproductive isolation through the evolution of barrier traits as a response to selection against maladaptive hybridization, but the genetic changes associated with this process remain largely unexplored. Building upon the increasing evidence for a role of structural variants in adaptation and speciation, we addressed the role of copy-number variation in the reinforcement of sexual isolation evidenced between the two European subspecies of the house mouse. We characterized copy-number divergence between populations of Mus musculus musculus that display assortative mate choice, and those that do not, using whole-genome resequencing data. Updating methods to detect deletions and tandem duplications (collectively: copy-number variants, CNVs) in Pool-Seq data, we developed an analytical pipeline dedicated to identifying genomic regions showing the expected pattern of copy-number displacement under a reinforcement scenario. This strategy allowed us to detect 1824 deletions and seven tandem duplications that showed extreme differences in frequency between behavioural classes across replicate comparisons. A subset of 480 deletions and four tandem duplications were specifically associated with the derived trait of assortative mate choice. These 'Choosiness-associated' CNVs occur in hundreds of genes. Consistent with our hypothesis, such genes included olfactory receptors potentially involved in the olfactory-based assortative mate choice in this system as well as one gene, Sp110, that is known to show patterns of differential expression between behavioural classes in an organ used in mate choice-the vomeronasal organ. These results demonstrate that fine-scale structural changes are common and highly variable within species, despite being under-studied, and may be important targets of reinforcing selection in this system and others. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Genomic Resources for Darters (Percidae: Etheostominae) Provide Insight into Postzygotic Barriers Implicated in Speciation

Comparative genomic approaches are increasingly being used to study the evolution of reproductive barriers in nonmodel species. Although numerous studies have examined prezygotic isolation in darters (Percidae), investigations into postzygotic barriers have remained rare due to long generation times and a lack of genomic resources. Orangethroat and rainbow darters naturally hybridize and provide a remarkable example of male-driven speciation via character displacement. Backcross hybrids suffer from high mortality, which appears to promote behavioral isolation in sympatry. To investigate the genomic architecture of postzygotic isolation, we used Illumina and PacBio sequencing to generate a chromosome-level, annotated assembly of the orangethroat darter genome and high-density linkage maps for orangethroat and rainbow darters. We also analyzed genome-wide RADseq data from wild-caught adults of both species and laboratory-generated backcrosses to identify genomic regions associated with hybrid incompatibles. Several putative chromosomal translocations and inversions were observed between orangethroat and rainbow darters, suggesting structural rearrangements may underlie postzygotic isolation. We also found evidence of selection against recombinant haplotypes and transmission ratio distortion in backcross hybrid genomes, providing further insight into the genomic architecture of genetic incompatibilities. Notably, regions with high levels of genetic divergence between species were enriched for genes associated with developmental and meiotic processes, providing strong candidates for postzygotic isolating barriers. These findings mark significant contributions to our understanding of the genetic basis of reproductive isolation between species undergoing character displacement. Furthermore, the genomic resources presented here will be instrumental for studying speciation in darters, the most diverse vertebrate group in North America.

Evidence for adaptive introgression of exons across a hybrid swarm in deer

Background

Secondary contact between closely related lineages can result in a variety of outcomes, including hybridization, depending upon the strength of reproductive barriers. By examining the extent to which different parts of the genome introgress, it is possible to infer the strength of selection and gain insight into the evolutionary trajectory of lineages. Following secondary contact approximately 8000 years ago in the Pacific Northwest, mule deer (Odocoileus hemionus hemionus) and black-tailed deer (O. h. columbianus) formed a hybrid swarm along the Cascade mountain range despite substantial differences in body size (up to two times) and habitat preference. In this study, we examined genetic population structure, extent of introgression, and selection pressures in freely interbreeding populations of mule deer and black-tailed deer using mitochondrial DNA sequences, 9 microsatellite loci, and 95 SNPs from protein-coding genes.

Results

We observed bi-directional hybridization and classified approximately one third of the 172 individuals as hybrids, almost all of which were beyond the F1 generation. High genetic differentiation between black-tailed deer and mule deer at protein-coding genes suggests that there is positive divergent selection, though selection on these loci is relatively weak. Contrary to predictions, there was not greater selection on protein-coding genes thought to be associated with immune function and mate choice. Geographic cline analyses were consistent across genetic markers, suggesting long-term stability (over hundreds of generations), and indicated that the center of the hybrid swarm is 20-30 km to the east of the Cascades ridgeline, where there is a steep ecological transition from wet, forested habitat to dry, scrub habitat.

Conclusions

Our data are consistent with a genetic boundary between mule deer and black-tailed deer that is porous but maintained by many loci under weak selection having a substantial cumulative effect. The absence of clear reproductive barriers and the consistent centering of geographic clines at a sharp ecotone suggests that ecology is a driver of hybrid swarm dynamics. Adaptive introgression in this study (and others) promotes gene flow and provides valuable insight into selection strength on specific genes and the evolutionary trajectory of hybridizing taxa.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1497-x) contains supplementary material, which is available to authorized users.

Genome-Wide Assessment of Diversity and Divergence Among Extant Galapagos Giant Tortoise Species

Genome-wide assessments allow for fuller characterization of genetic diversity, finer-scale population delineation, and better detection of demographically significant units to guide conservation compared with those based on "traditional" markers. Galapagos giant tortoises (Chelonoidis spp.) have long provided a case study for how evolutionary genetics may be applied to advance species conservation. Ongoing efforts to bolster tortoise populations, which have declined by 90%, have been informed by analyses of mitochondrial DNA sequence and microsatellite genotypic data, but could benefit from genome-wide markers. Taking this next step, we used double-digest restriction-site associated DNA sequencing to collect genotypic data at >26000 single nucleotide polymorphisms (SNPs) for 117 individuals representing all recognized extant Galapagos giant tortoise species. We then quantified genetic diversity, population structure, and compared results to estimates from mitochondrial DNA and microsatellite loci. Our analyses detected 12 genetic lineages concordant with the 11 named species as well as previously described structure within one species, C. becki. Furthermore, the SNPs provided increased resolution, detecting admixture in 4 individuals. SNP-based estimates of diversity and differentiation were significantly correlated with those derived from nuclear microsatellite loci and mitochondrial DNA sequences. The SNP toolkit presented here will serve as a resource for advancing efforts to understand tortoise evolution, species radiations, and aid conservation of the Galapagos tortoise species complex.

Chemical cues and genetic divergence in insects on plants: conceptual cross pollination between mutualistic and antagonistic systems

Cascading or reciprocal genetic diversification of herbivores, parasitoids, and pollinators can track chemotypic variation in host resources, and can lead to non-overlapping communities. Because plants simultaneously interact with both pollinators and herbivores, models investigating the genetic divergence of antagonistic herbivores and mutualistic pollinators should be merged in order to study how both processes interact using a common conceptual and methodological approach. We expect insects to mediate divergence in many systems, with outcomes depending on the level of pollinator or herbivore specialisation, and the relative selective pressures they impose. Applying approaches widely used to study insect pollinators, for example genomic tools and integration of behavioural, genetic and chemical data, to both pollinators and herbivores in the same system will facilitate our understanding of patterns of genetic divergence across multiple interacting species.

Genomic signatures of local adaptation to the degree of environmental predictability in rotifers

Environmental fluctuations are ubiquitous and thus essential for the study of adaptation. Despite this, genome evolution in response to environmental fluctuations —and more specifically to the degree of environmental predictability– is still unknown. Saline lakes in the Mediterranean region are remarkably diverse in their ecological conditions, which can lead to divergent local adaptation patterns in the inhabiting aquatic organisms. The facultatively sexual rotifer Brachionus plicatilis shows diverging local adaptation in its life-history traits in relation to estimated environmental predictability in its habitats. Here, we used an integrative approach —combining environmental, phenotypic and genomic data for the same populations– to understand the genomic basis of this diverging adaptation. Firstly, a novel draft genome for B. plicatilis was assembled. Then, genome-wide polymorphisms were studied using genotyping by sequencing on 270 clones from nine populations in eastern Spain. As a result, 4,543 high-quality SNPs were identified and genotyped. More than 90 SNPs were found to be putatively under selection with signatures of diversifying and balancing selection. Over 140 SNPs were correlated with environmental or phenotypic variables revealing signatures of local adaptation, including environmental predictability. Putative functions were associated to most of these SNPs, since they were located within annotated genes. Our results reveal associations between genomic variation and the degree of environmental predictability, providing genomic evidence of adaptation to local conditions in natural rotifer populations.

Speciation with gene flow via cycles of isolation and migration: insights from multiple mangrove taxa

Allopatric speciation requiring an unbroken period of geographical isolation has been the standard model of neo-Darwinism. While doubts have been repeatedly raised, strict allopatry without any gene flow remains a plausible mechanism in most cases. To rigorously reject strict allopatry, genomic sequences superimposed on the geological records of a well-delineated geographical barrier are necessary. The Strait of Malacca, narrowly connecting the Pacific and Indian Ocean coasts, serves at different times either as a geographical barrier or a conduit of gene flow for coastal/marine species. We surveyed 1700 plants from 29 populations of 5 common mangrove species by large-scale DNA sequencing and added several whole-genome assemblies. Speciation between the two oceans is driven by cycles of isolation and gene flow due to the fluctuations in sea level leading to the opening/closing of the Strait to ocean currents. Because the time required for speciation in mangroves is longer than the isolation phases, speciation in these mangroves has proceeded through many cycles of mixing-isolation-mixing, or MIM, cycles. The MIM mechanism, by relaxing the condition of no gene flow, can promote speciation in many more geographical features than strict allopatry can. Finally, the MIM mechanism of speciation is also efficient, potentially yielding mⁿ (m > 1) species after n cycles.

Can genomic data alone tell us whether speciation happened with gene flow?

The allopatric model, which requires a period of geographical isolation for speciation to complete, has been the standard model in the modern era. Recently, "speciation with gene flow" has been widely discussed in relation to the model of "strict allopatry" and the level of DNA divergence across genomic regions. We wish to caution that genomic data by themselves may only permit the rejection of the simplest form of allopatry. Even a slightly more complex and realistic model that starts with subdivided populations would be impossible to reject by the genomic data alone. To resolve this central issue of speciation, other forms of observations such as the sequencing of reproductive isolation genes or the identification of geographical barrier(s) will be necessary.

Highly localized divergence within supergenes in Atlantic cod (Gadus morhua) within the Gulf of Maine

BACKGROUND

Atlantic cod (Gadus morhua), is known to vary genetically across the North Atlantic, Greenland, and Newfoundland. This genetic variation occurs both spatially and temporally through decades of heavy fishing, and is concentrated in three linkage disequilibrium blocks, previously defined by pedigreed linkage mapping analysis. Variation within these genomic regions is correlated with both seawater temperature and behavioral ecotype. The full extent and nature of these linkage groups is important information for interpreting cod genetic structure as a tool for future fisheries management.

RESULTS

We conducted whole genome sequencing for 31 individual cod from three sub-populations in the Gulf of Maine. Across the genome, we found 3,390,654 intermediate to high frequency Single Nucleotide Polymorphisms (SNPs). We show that pairwise linkage analysis among these SNPs is a powerful tool to detect linkage disequilibrium clusters by recovering the three previously detected linkage groups and identifying the 1031 genes contained therein. Across these genes, we found significant population differentiation among spawning groups in the Gulf of Maine and between Georges Bank and Gulf of Maine. Coordinated divergence among these genes and their differentiation at both short and long spatial scales suggests that they are acting as linked supergenes in local adaptation of cod populations.

CONCLUSIONS

Differentiation between SNPs in linkage disequilibrium blocks is the major signal of genetic differentiation between all groups tested within the Gulf of Maine. Our data provide a map of genes contained in these blocks, allowing an enhanced search for neutral genetic structure for demographic inference and fisheries modeling. Patterns of selection and the history of populations may be possible to identify in cod using this description of linkage disequilibrium blocks and future data sets to robustly separate neutral and selected genetic markers.

Targeted re-sequencing confirms the importance of chemosensory genes in aphid host race differentiation

Host‐associated races of phytophagous insects provide a model for understanding how adaptation to a new environment can lead to reproductive isolation and speciation, ultimately enabling us to connect barriers to gene flow to adaptive causes of divergence. The pea aphid (Acyrthosiphon pisum) comprises host races specializing on legume species and provides a unique system for examining the early stages of diversification along a gradient of genetic and associated adaptive divergence. As host choice produces assortative mating, understanding the underlying mechanisms of choice will contribute directly to understanding of speciation. As host choice in the pea aphid is likely mediated by smell and taste, we use capture sequencing and SNP genotyping to test for the role of chemosensory genes in the divergence between eight host plant species across the continuum of differentiation and sampled at multiple locations across western Europe. We show high differentiation of chemosensory loci relative to control loci in a broad set of pea aphid races and localities, using a model‐free approach based on principal component analysis. Olfactory and gustatory receptors form the majority of highly differentiated genes and include loci that were already identified as outliers in a previous study focusing on the three most closely related host races. Consistent indications that chemosensory genes may be good candidates for local adaptation and barriers to gene flow in the pea aphid open the way to further investigations aiming to understand their impact on gene flow and to determine their precise functions in response to host plant metabolites.

Targeted resequencing reveals geographical patterns of differentiation for loci implicated in parallel evolution

Parallel divergence and speciation provide evidence for the role of divergent selection in generating biological diversity. Recent studies indicate that parallel phenotypic divergence may not have the same genetic basis in different geographical locations - 'outlier loci' (loci potentially affected by divergent selection) are often not shared among parallel instances of phenotypic divergence. However, limited sharing may be due, in part, to technical issues if false-positive outliers occur. Here, we test this idea in the marine snail Littorina saxatilis, which has evolved two partly isolated ecotypes (adapted to crab predation vs. wave action) in multiple locations independently. We argue that if the low extent of sharing observed in earlier studies in this system is due to sampling effects, we expect outliers not to show elevated FST when sequenced in new samples from the original locations and also not to follow predictable geographical patterns of elevated FST . Following a hierarchical sampling design (within vs. between country), we applied capture sequencing, targeting outliers from earlier studies and control loci. We found that outliers again showed elevated levels of FST in their original location, suggesting they were not generated by sampling effects. Outliers were also likely to show increased FST in geographically close locations, which may be explained by higher levels of gene flow or shared ancestral genetic variation compared with more distant locations. However, in contrast to earlier findings, we also found some outlier types to show elevated FST in geographically distant locations. We discuss possible explanations for this unexpected result.

Unravelling the effects of gene flow and selection in highly connected populations of the silver-lip pearl oyster (Pinctada maxima)

Many marine organisms often display weak levels of population genetic structuring as a result of both environmental characteristics (e.g., ocean currents) and life history traits (e.g., widely dispersed planktonic larval stages) maintaining high levels of gene flow. This can lead to the assumption that these organisms can be managed as a single stock based on high levels of population connectivity. However, this neglects to account for other micro-evolutionary forces such as selection, which also shape these populations. This study utilizes 1130 genome-wide SNP loci to unravel the effects of gene flow and selection shaping three highly connected populations of the silver-lip pearl oyster (Pinctada maxima) in the ecologically and economically important Indo-Pacific region (Aru, Bali, and West Papua). Twenty-two loci under directional selection were identified amongst the populations, providing further supporting evidence of strong local adaptation (i.e., G×E effects) among populations in this region. Global Fst values for directional outliers (0.348) were up to eight times greater than for neutral markers (0.043). Pairwise Fst comparisons between Aru and Bali revealed the largest directional differences (0.488), while Bali and West Papua had the least (0.062). Unrooted neighbour-joining (NJ) distance trees and genetic diversity indices of directional outliers revealed that individuals from Bali and West Papua had reduced allelic variation (MAFavg=0.144, Ho=0.238 and MAFavg=0.232, Ho=0.369, respectively) compared to Aru (MAFavg=0.292, Ho=0.412). This indicates that directional selection is most likely acting upon genetic variation within the Bali and West Papua populations. NJ distance trees, discriminant analysis of principal components, and Fst analyses of directional outliers revealed two divergent groups ("Bali/West Papua"; "Aru") that had previously gone unrecognized. This study not only illustrates that relatively strong local adaptive forces are occurring despite high gene flow, but identifies the populations that are most likely experiencing selection. Additionally, this study highlights the need to understand all micro-evolutionary forces acting on populations when resolving stock structure.

A population genomic scan in Chorthippus grasshoppers unveils previously unknown phenotypic divergence

Understanding the genetics of speciation and the processes that drive it is a central goal of evolutionary biology. Grasshoppers of the Chorthippus species group differ strongly in calling song (and corresponding female preferences) but are exceedingly similar in other characteristics such as morphology. Here, we performed a population genomic scan on three Chorthippus species (Chorthippus biguttulus, C. mollis and C. brunneus) to gain insight into the genes and processes involved in divergence and speciation in this group. Using an RNA-seq approach, we examined functional variation between the species by calling SNPs for each of the three species pairs and using FST -based approaches to identify outliers. We found approximately 1% of SNPs in each comparison to be outliers. Between 37% and 40% of these outliers were nonsynonymous SNPs (as opposed to a global level of 17%) indicating that we recovered loci under selection. Among the outliers were several genes that may be involved in song production and hearing as well as genes involved in other traits such as food preferences and metabolism. Differences in food preferences between species were confirmed with a behavioural experiment. This indicates that multiple phenotypic differences implicating multiple evolutionary processes (sexual selection and natural selection) are present between the species.

Genomic patterns of species diversity and divergence in Eucalyptus

We examined genome-wide patterns of DNA sequence diversity and divergence among six species of the important tree genus Eucalyptus and investigated their relationship with genomic architecture. Using c. 90 range-wide individuals of each Eucalyptus species (E. grandis, E. urophylla, E. globulus, E. nitens, E. dunnii and E. camaldulensis), genetic diversity and divergence were estimated from 2840 polymorphic diversity arrays technology markers covering the 11 chromosomes. Species differentiating markers (SDMs) identified in each of 15 pairwise species comparisons, along with species diversity (HHW ) and divergence (FST ), were projected onto the E. grandis reference genome. Across all species comparisons, SDMs totalled 1.1-5.3% of markers and were widely distributed throughout the genome. Marker divergence (FST and SDMs) and diversity differed among and within chromosomes. Patterns of diversity and divergence were broadly conserved across species and significantly associated with genomic features, including the proximity of markers to genes, the relative number of clusters of tandem duplications, and gene density within or among chromosomes. These results suggest that genomic architecture influences patterns of species diversity and divergence in the genus. This influence is evident across the six species, encompassing diverse phylogenetic lineages, geography and ecology.

Towards the identification of the loci of adaptive evolution

1. Establishing the genetic and molecular basis underlying adaptive traits is one of the major goals of evolutionary geneticists in order to understand the connection between genotype and phenotype and elucidate the mechanisms of evolutionary change. Despite considerable effort to address this question, there remain relatively few systems in which the genes shaping adaptations have been identified. 2. Here, we review the experimental tools that have been applied to document the molecular basis underlying evolution in several natural systems, in order to highlight their benefits, limitations and suitability. In most cases, a combination of DNA, RNA and functional methodologies with field experiments will be needed to uncover the genes and mechanisms shaping adaptation in nature.

Genetic evidence for ecological divergence in kokanee salmon

The evolution of locally adapted phenotypes among populations that experience divergent selective pressures is a central mechanism for generating and maintaining biodiversity. Recently, the advent of high-throughput DNA sequencing technology has provided tools for investigating the genetic basis of this process in natural populations of nonmodel organisms. Kokanee, the freshwater form of sockeye salmon (Oncorhynchus nerka), occurs as two reproductive ecotypes, which differ in spawning habitat (tributaries vs. shorelines); however, outside of the spawning season the two ecotypes co-occur in many lakes and lack diagnostic morphological characteristics. We used restriction site-associated DNA (RAD) sequencing to identify 6145 SNPs and genotype kokanee from multiple spawning sites in Okanagan Lake (British Columbia, Canada). Outlier tests revealed 18 loci putatively under divergent selection between ecotypes, all of which exhibited temporally stable allele frequencies within ecotypes. Six outliers were annotated to sequences in the NCBI database, two of which matched genes associated with early development. There was no evidence for neutral genetic differentiation; however, outlier loci demonstrated significant structure with respect to ecotype and had high assignment accuracy in mixed composition simulations. The absence of neutral structure combined with a small number of highly divergent outlier loci is consistent with theoretical predictions for the early stages of ecological divergence. These outlier loci were then applied to a realistic fisheries scenario in which additional RAD sequencing was used to genotype kokanee collected by trawl in Okanagan Lake, providing preliminary evidence that this approach may be an effective tool for conservation and management.

Do island plant populations really have lower genetic variation than mainland populations? Effects of selection and distribution range on genetic diversity estimates

Ecological and evolutionary studies largely assume that island populations display low levels of neutral genetic variation. However, this notion has only been formally tested in a few cases involving plant taxa, and the confounding effect of selection on genetic diversity (GD) estimates based on putatively neutral markers has typically been overlooked. Here, we generated nuclear microsatellite and plastid DNA sequence data in Periploca laevigata, a plant taxon with an island-mainland distribution area, to (i) investigate whether selection affects GD estimates of populations across contrasting habitats; and (ii) test the long-standing idea that island populations have lower GD than their mainland counterparts. Plastid data showed that colonization of the Canary Islands promoted strong lineage divergence within P. laevigata, which was accompanied by selective sweeps at several nuclear microsatellite loci. Inclusion of loci affected by strong divergent selection produced a significant downward bias in the GD estimates of the mainland lineage, but such underestimates were substantial (>14%) only when more than one loci under selection were included in the computations. When loci affected by selection were removed, we did not find evidence that insular Periploca populations have less GD than their mainland counterparts. The analysis of data obtained from a comprehensive literature survey reinforced this result, as overall comparisons of GD estimates between island and mainland populations were not significant across plant taxa (N = 66), with the only exception of island endemics with narrow distributions. This study suggests that identification and removal of markers potentially affected by selection should be routinely implemented in estimates of GD, particularly if different lineages are compared. Furthermore, it provides compelling evidence that the expectation of low GD cannot be generalized to island plant populations.

Selection on outlier loci and their association with adaptive phenotypes in Littorina saxatilis contact zones

A fundamental issue in speciation research is to evaluate phenotypic variation and the genomics driving the evolution of reproductive isolation between sister taxa. Above all, hybrid zones are excellent study systems for researchers to examine the association of genetic differentiation, phenotypic variation and the strength of selection. We investigated two contact zones in the marine gastropod Littorina saxatilis and utilized landmark-based geometric morphometric analysis together with amplified fragment length polymorphism (AFLP) markers to assess phenotypic and genomic divergence between ecotypes under divergent selection. From genetic markers, we calculated the cline width, linkage disequilibrium and the average effective selection on a locus. Additionally, we conducted an association analysis linking the outlier loci and phenotypic variation between ecotypes and show that a proportion of outlier loci are associated with key adaptive phenotypic traits.

Do the same genes underlie parallel phenotypic divergence in different Littorina saxatilis populations?

Parallel patterns of adaptive divergence and speciation are cited as powerful evidence for the role of selection driving these processes. However, it is often not clear whether parallel phenotypic divergence is underlain by parallel genetic changes. Here, we asked about the genetic basis of parallel divergence in the marine snail Littorina saxatilis, which has repeatedly evolved coexisting ecotypes adapted to either crab predation or wave action. We sequenced the transcriptome of snails of both ecotypes from three distant geographical locations (Spain, Sweden and United Kingdom) and mapped the reads to the L. saxatilis reference genome. We identified genomic regions potentially under divergent selection between ecotypes within each country, using an outlier approach based on F(ST) values calculated per locus. In line with previous studies indicating that gene reuse is generally common, we expected to find extensive sharing of outlier loci due to recent shared ancestry and gene flow between at least two of the locations in our study system. Contrary to our expectations, we found that most outliers were country specific, suggesting that much of the genetic basis of divergence is not shared among locations. However, we did find that more outliers were shared than expected by chance and that differentiation of shared outliers is often generated by the same SNPs. We discuss two mechanisms potentially explaining the limited amount of sharing we observed. First, a polygenic basis of divergent traits might allow for multiple distinct molecular mechanisms generating the same phenotypic patterns. Second, additional, location-specific axes of selection that we did not focus on in this study may produce distinct patterns of genetic divergence within each site.

Crater Lake Apoyo revisited--population genetics of an emerging species flock

The polytypic Nicaraguan Midas cichlids (Amphilophus cf. citrinellus) have been established as a model system for studying the mechanisms of speciation and patterns of diversification in allopatry and sympatry. The species assemblage in Crater Lake Apoyo has been accepted as a textbook example for sympatric speciation. Here, we present a first comprehensive data set of population genetic (mtDNA & AFLPs) proxies of species level differentiation for a representative set of individuals of all six endemic Amphilophus species occurring in Crater Lake Apoyo. AFLP genetic differentiation was partitioned into a neutral and non-neutral component based on outlier-loci detection approaches, and patterns of species divergence were explored with Bayesian clustering methods. Substantial levels of admixture between species were detected, indicating different levels of reproductive isolation between the six species. Analysis of neutral genetic variation revealed several A. zaliosus as being introgressed by an unknown contributor, hereby rendering the sympatrically evolving L. Apoyo flock polyphyletic. This is contrasted by the mtDNA analysis delivering a clear monophyly signal with Crater Lake Apoyo private haplotypes characterising all six described species, but also demonstrating different demographic histories as inferred from pairwise mismatch distributions.

Genomic distribution of AFLP markers relative to gene locations for different eukaryotic species

BACKGROUND

Amplified fragment length polymorphism (AFLP) markers are frequently used for a wide range of studies, such as genome-wide mapping, population genetic diversity estimation, hybridization and introgression studies, phylogenetic analyses, and detection of signatures of selection. An important issue to be addressed for some of these fields is the distribution of the markers across the genome, particularly in relation to gene sequences.

RESULTS

Using in-silico restriction fragment analysis of the genomes of nine eukaryotic species we characterise the distribution of AFLP fragments across the genome and, particularly, in relation to gene locations. First, we identify the physical position of markers across the chromosomes of all species. An observed accumulation of fragments around (peri) centromeric regions in some species is produced by repeated sequences, and this accumulation disappears when AFLP bands rather than fragments are considered. Second, we calculate the percentage of AFLP markers positioned within gene sequences. For the typical EcoRI/MseI enzyme pair, this ranges between 28 and 87% and is usually larger than that expected by chance because of the higher GC content of gene sequences relative to intergenic ones. In agreement with this, the use of enzyme pairs with GC-rich restriction sites substantially increases the above percentages. For example, using the enzyme system SacI/HpaII, 86% of AFLP markers are located within gene sequences in A. thaliana, and 100% of markers in Plasmodium falciparun. We further find that for a typical trait controlled by 50 genes of average size, if 1000 AFLPs are used in a study, the number of those within 1 kb distance from any of the genes would be only about 1-2, and only about 50% of the genes would have markers within that distance.

CONCLUSIONS

The high coverage of AFLP markers across the genomes and the high proportion of markers within or close to gene sequences make them suitable for genome scans and detecting large islands of differentiation in the genome. However, for specific traits, the percentage of AFLP markers close to genes can be rather small. Therefore, genome scans directed towards the search of markers closely linked to selected loci can be a difficult task in many instances.

Genetic structure and demographic history of Colletotrichum gloeosporioides sensu lato and C. truncatum isolates from Trinidad and Mexico

BACKGROUND

C. gloeosporioides sensu lato is one of the most economically important post-harvest diseases affecting papaya production worldwide. There is currently no information concerning the genetic structure or demographic history of this pathogen in any of the affected countries. Knowledge of molecular demographic parameters for different populations will improve our understanding of the biogeographic history as well as the evolutionary and adaptive potential of these pathogens. In this study, sequence data for ACT, GPDH, β-TUB and ITS gene regions were analyzed for C. gloeosporioides sensu lato and C. truncatum isolates infecting papaya in Trinidad and Mexico in order to determine the genetic structure and demographic history of these populations.

RESULTS

The data indicated that Mexico is the ancestral C. gloeosporioides sensu lato population with asymmetrical migration to Trinidad. Mexico also had the larger effective population size but, both Mexico and Trinidad populations exhibited population expansion. Mexico also had greater nucleotide diversity and high levels of diversity for each gene. There was significant sub-division of the Trinidad and Mexico populations and low levels of genetic divergence among populations for three of the four gene regions; β-TUB was shown to be under positive selection. There were also dissimilar haplotype characteristics for both populations. Mutation may play a role in shaping the population structure of C. gloeosporioides sensu lato isolates from Trinidad and from Mexico, especially with respect to the ACT and GPDH gene regions. There was no evidence of gene flow between the C. truncatum populations and it is possible that the Mexico and Trinidad populations emerged independently of each other.

CONCLUSIONS

The study revealed relevant information based on the genetic structure as well as the demographic history of two fungal pathogens infecting papaya, C. gloeosporioides sensu lato and C. truncatum, in Trinidad and Mexico. Understanding the genetic structure of pathogen populations will assist in determining the evolutionary potential of the pathogen and in identifying which evolutionary forces may have the greatest impact on durability of resistance. Intervention strategies that target these evolutionary forces would prove to be the most practical.

The population structure and recent colonization history of Oregon threespine stickleback determined using restriction-site associated DNA-sequencing

Understanding how genetic variation is partitioned across genomes within and among populations is a fundamental problem in ecological and evolutionary genetics. To address this problem, we studied the threespine stickleback fish, which has repeatedly undergone parallel phenotypic and genetic differentiation when oceanic fish have invaded freshwater habitats. While significant evolutionary genetic research has been performed using stickleback from geographic regions that have been deglaciated in the last 20 000 years, less research has focused on freshwater populations that predate the last glacial maximum. We performed restriction-site associated DNA-sequencing (RAD-seq) based population genomic analyses on stickleback from across Oregon, which was not glaciated during the last maximum. We sampled stickleback from coastal, Willamette Basin and central Oregon sites, analysed their genetic diversity using RAD-seq, performed structure analyses, reconstructed their phylogeographic history and tested the hypothesis of recent stickleback introduction into central Oregon, where incidence of this species was only recently documented. Our results showed a clear phylogeographic break between coastal and inland populations, with oceanic populations exhibiting the lowest levels of divergence from one another. Willamette Basin and central Oregon populations formed a clade of closely related populations, a finding consistent with a recent introduction of stickleback into central Oregon. Finally, genome-wide analysis of genetic diversity (π) and correlations of alleles within individuals in subpopulations (FIS) supported a role for introgressive hybridization in coastal populations and a recent expansion in central Oregon. Our results exhibit the power of next-generation sequencing genomic approaches such as RAD-seq to identify both historical population structure and recent colonization history.

Genomic signatures of relaxed disruptive selection associated with speciation reversal in whitefish

Background

Speciation reversal: the erosion of species differentiation via an increase in introgressive hybridization due to the weakening of previously divergent selection regimes, is thought to be an important, yet poorly understood, driver of biodiversity loss. Our study system, the Alpine whitefish (Coregonus spp.) species complex is a classic example of a recent postglacial adaptive radiation: forming an array of endemic lake flocks, with the independent origination of similar ecotypes among flocks. However, many of the lakes of the Alpine radiation have been seriously impacted by anthropogenic nutrient enrichment, resulting in a collapse in neutral genetic and phenotypic differentiation within the most polluted lakes. Here we investigate the effects of eutrophication on the selective forces that have shaped this radiation, using population genomics. We studied eight sympatric species assemblages belonging to five independent parallel adaptive radiations, and one species pair in secondary contact. We used AFLP markers, and applied F_ST outlier (BayeScan, Dfdist) and logistic regression analyses (MatSAM), to identify candidate regions for disruptive selection in the genome and their associations with adaptive traits within each lake flock. The number of outlier and adaptive trait associated loci identified per lake were then regressed against two variables (historical phosphorus concentration and contemporary oxygen concentration) representing the strength of eutrophication.

Results

Whilst we identify disruptive selection candidate regions in all lake flocks, we find similar trends, across analysis methods, towards fewer disruptive selection candidate regions and fewer adaptive trait/candidate loci associations in the more polluted lakes.

Conclusions

Weakened disruptive selection and a concomitant breakdown in reproductive isolating mechanisms in more polluted lakes has lead to increased gene flow between coexisting Alpine whitefish species. We hypothesize that the resulting higher rates of interspecific recombination reduce either the number or extent of genomic islands of divergence surrounding loci evolving under disruptive natural selection. This produces the negative trend seen in the number of selection candidate loci recovered during genome scans of whitefish species flocks, with increasing levels of anthropogenic eutrophication: as the likelihood decreases that AFLP restriction sites will fall within regions of heightened genomic divergence and therefore be classified as F_ST outlier loci. This study explores for the first time the potential effects of human-mediated relaxation of disruptive selection on heterogeneous genomic divergence between coexisting species.

Crossing the species barrier: genomic hotspots of introgression between two highly divergent Ciona intestinalis species

Inferring a realistic demographic model from genetic data is an important challenge to gain insights into the historical events during the speciation process and to detect molecular signatures of selection along genomes. Recent advances in divergence population genetics have reported that speciation in face of gene flow occurred more frequently than theoretically expected, but the approaches used did not account for genome-wide heterogeneity (GWH) in introgression rates. Here, we investigate the impact of GWH on the inference of divergence with gene flow between two cryptic species of the marine model Ciona intestinalis by analyzing polymorphism and divergence patterns in 852 protein-coding sequence loci. These morphologically similar entities are highly diverged molecular-wise, but evidence of hybridization has been reported in both laboratory and field studies. We compare various speciation models and test for GWH under the approximate Bayesian computation framework. Our results demonstrate the presence of significant extents of gene flow resulting from a recent secondary contact after >3 My of divergence in isolation. The inferred rates of introgression are relatively low, highly variable across loci and mostly unidirectional, which is consistent with the idea that numerous genetic incompatibilities have accumulated over time throughout the genomes of these highly diverged species. A genomic map of the level of gene flow identified two hotspots of introgression, that is, large genome regions of unidirectional introgression. This study clarifies the history and degree of isolation of two cryptic and partially sympatric model species and provides a methodological framework to investigate GWH at various stages of speciation process.

Microbial variome database: point mutations, adaptive or not, in bacterial core genomes

Analysis of genetic differences (gene presence/absence and nucleotide polymorphisms) among strains of a bacterial species is crucial to understanding molecular mechanisms of bacterial pathogenesis and selecting targets for novel antibacterial therapeutics. However, lack of genome-wide association studies on large and epidemiologically well-defined strain collections from the same species makes it difficult to identify the genes under positive selection and define adaptive polymorphisms in those genes. To address this need and to overcome existing limitations, we propose to create a "microbial variome"--a species-specific resource database of genomic variations based on molecular evolutionary analysis. Here, we present prototype variome databases of Escherichia coli and Salmonella enterica subspecies enterica (http://depts.washington.edu/sokurel/variome, last accessed March 26, 2013). The prototypes currently include the point mutations data of core protein-coding genes from completely sequenced genomes of 22 E. coli and 17 S. enterica strains. These publicly available databases allow for single- and multiple-field sorting, filtering, and searching of the gene variability data and the potential adaptive significance. Such resource databases would immensely help experimental research, clinical diagnostics, epidemiology, and environmental control of human pathogens.

Introgression of wing pattern alleles and speciation via homoploid hybridization in Heliconius butterflies: a review of evidence from the genome

The diverse Müllerian mimetic wing patterns of neotropical Heliconius (Nymphalidae) have been proposed to be not only aposematic signals to potential predators, but also intra- and interspecific recognition signals that allow the butterflies to maintain their specific identities, and which perhaps drive the process of speciation, as well. Adaptive features under differential selection that also serve as cues for assortative mating have been referred to as 'magic traits', which can drive ecological speciation. Such traits are expected to exhibit allelic differentiation between closely related species with ongoing gene flow, whereas unlinked neutral traits are expected to be homogenized to a greater degree by introgression. However, recent evidence suggests that interspecific hybridization among Heliconius butterflies may have resulted in adaptive introgression of these very same traits across species boundaries, and in the evolution of new species by homoploid hybrid speciation. The theory and data supporting various aspects of the apparent paradox of 'magic trait' introgression are reviewed, with emphasis on population genomic comparisons of Heliconius melpomene and its close relatives.

Outlier loci highlight the direction of introgression in oaks

Loci considered to be under selection are generally avoided in attempts to infer past demographic processes as they do not fit neutral model assumptions. However, opportunities to better reconstruct some aspects of past demography might thus be missed. Here we examined genetic differentiation between two sympatric European oak species with contrasting ecological dynamics (Quercus robur and Quercus petraea) with both outlier (i.e. loci possibly affected by divergent selection between species or by hitchhiking effects with genomic regions under selection) and nonoutlier loci. We sampled 855 individuals in six mixed forests in France and genotyped them with a set of 262 SNPs enriched with markers showing high interspecific differentiation, resulting in accurate species delimitation. We identified between 13 and 74 interspecific outlier loci, depending on the coalescent simulation models and parameters used. Greater genetic diversity was predicted in Q. petraea (a late-successional species) than in Q. robur (an early successional species) as introgression should theoretically occur predominantly from the resident species to the invading species. Remarkably, this prediction was verified with outlier loci but not with nonoutlier loci. We suggest that the lower effective interspecific gene flow at loci showing high interspecific divergence has better preserved the signal of past asymmetric introgression towards Q. petraea caused by the species' contrasting dynamics. Using markers under selection to reconstruct past demographic processes could therefore have broader potential than generally recognized.

Reproductive isolation and local adaptation quantified for a chromosome inversion in a malaria mosquito

Chromosome inversions have long been thought to be involved in speciation and local adaptation. We have little quantitative information, however, about the effects that inversion polymorphisms have on reproductive isolation and viability. Here we provide the first estimates from any organism for the total amount of reproductive isolation associated with an inversion segregating in natural populations. We sampled chromosomes from 751 mosquitoes of the malaria vector Anopheles funestus along a 1421 km transect in Cameroon that traverses savannah, highland, and rainforest ecological zones. We then developed a series of population genetic models that account for selection, migration, and assortative mating, and fit the models to the data using likelihood. Results from the best-fit models suggest there is strong local adaptation, with relative viabilities of homozygotes ranging from 25% to 130% compared to heterozygotes. Viabilities vary qualitatively between regions: the inversion is underdominant in the savannah, whereas in the highlands it is overdominant. The inversion is also implicated in strong assortative mating. In the savannah, the two homozygote forms show 92% reproductive isolation, suggesting that this one inversion can generate most of the genetic barriers needed for speciation.

Recombinant hybrids retain heterozygosity at many loci: new insights into the genomics of reproductive isolation in Populus

The maintenance of species barriers in the face of gene flow is often thought to result from strong selection against intermediate genotypes, thereby preserving genetic differentiation. Most speciation genomic studies thus aim to identify exceptionally divergent loci between populations, but divergence will be affected by many processes other than reproductive isolation (RI) and speciation. Through genomic studies of recombinant hybrids sampled in the wild, genetic variation associated with RI can be observed in situ, because selection against incompatible genotypes will leave detectable patterns of variation in the hybrid genomes. To better understand the mechanisms directly involved in RI, we investigated three natural 'replicate' hybrid zones between two divergent Populus species via locus-specific patterns of ancestry across recombinant hybrid genomes. As expected, genomic patterns in hybrids and their parental species were consistent with the presence of underdominant selection at several genomic regions. Surprisingly, many loci displayed greatly increased between-species heterozygosity in recombinant hybrids despite striking genetic differentiation between the parental genomes, the opposite of what would be expected with selection against intermediate genotypes. Only a limited, reproducible set of genotypic combinations was present in hybrid genomes across localities. In the absence of clearly delimited 'hybrid habitats', our results suggest that complex epistatic interactions within genomes play an important role in advanced stages of RI between these ecologically divergent forest trees. This calls for more genomic studies that test for unusual patterns of genomic ancestry in hybridizing species.

Do highly divergent loci reside in genomic regions affecting reproductive isolation? A test using next-generation sequence data in Timema stick insects

Background

Genetic divergence during speciation with gene flow is heterogeneous across the genome, with some regions exhibiting stronger differentiation than others. Exceptionally differentiated regions are often assumed to experience reduced introgression, i.e., reduced flow of alleles from one population into another because such regions are affected by divergent selection or cause reproductive isolation. In contrast, the remainder of the genome can be homogenized by high introgression. Although many studies have documented variation across the genome in genetic differentiation, there are few tests of this hypothesis that explicitly quantify introgression. Here, we provide such a test using 38,304 SNPs in populations of Timema cristinae stick insects. We quantify whether loci that are highly divergent between geographically separated (‘allopatric’) populations exhibit unusual patterns of introgression in admixed populations. To the extent this is true, highly divergent loci between allopatric populations contribute to reproductive isolation in admixed populations.

Results

As predicted, we find a substantial association between locus-specific divergence between allopatric populations and locus-specific introgression in admixed populations. However, many loci depart from this relationship, sometimes strongly so. We also report evidence for selection against foreign alleles due to local adaptation.

Conclusions

Loci that are strongly differentiated between allopatric populations sometimes contribute to reproductive isolation in admixed populations. However, geographic variation in selection and local adaptation, in aspects of genetic architecture (such as organization of genes, recombination rate variation, number and effect size of variants contributing to adaptation, etc.), and in stochastic evolutionary processes such as drift can cause strong differentiation of loci that do not always contribute to reproductive isolation. The results have implications for the theory of ‘genomic islands of speciation’.

Challenges and pitfalls in the characterization of anonymous outlier AFLP markers in non-model species: lessons from an ocellated lizard genome scan

In the last few years, dozens of studies have documented the detection of loci influenced by selection from genome scans in a wide range of non-model species. Many of those studies used amplified fragment length polymorphism (AFLP) markers, which became popular for being easily applicable to any organism. However, because they are anonymous markers, AFLPs impose many challenges for their isolation and identification. Most recent AFLP genome scans used capillary electrophoresis (CE), which adds even more obstacles to the isolation of bands with a specific size for sequencing. These caveats might explain the extremely low number of studies that moved from the detection of outlier AFLP markers to their actual isolation and characterization. We document our efforts to characterize a set of outlier AFLP markers from a previous genome scan with CE in ocellated lizards (Lacerta lepida). Seven outliers were successfully isolated, cloned and sequenced. Their sequences are noncoding and show internal indels or polymorphic repetitive elements (microsatellites). Three outliers were converted into codominant markers by using specific internal primers to sequence and screen population variability from undigested DNA. Amplification in closely related lizard species was also achieved, revealing remarkable interspecific conservation in outlier loci sequences. We stress the importance of following up AFLP genome scans to validate selection signatures of outlier loci, but also report the main challenges and pitfalls that may be faced during the process.

AFLP genome scan to detect genetic structure and candidate loci under selection for local adaptation of the invasive weed Mikania micrantha

Why some species become successful invaders is an important issue in invasive biology. However, limited genomic resources make it very difficult for identifying candidate genes involved in invasiveness. Mikania micrantha H.B.K. (Asteraceae), one of the world's most invasive weeds, has adapted rapidly in response to novel environments since its introduction to southern China. In its genome, we expect to find outlier loci under selection for local adaptation, critical to dissecting the molecular mechanisms of invasiveness. An explorative amplified fragment length polymorphism (AFLP) genome scan was used to detect candidate loci under selection in 28 M. micrantha populations across its entire introduced range in southern China. We also estimated population genetic parameters, bottleneck signatures, and linkage disequilibrium. In binary characters, such as presence or absence of AFLP bands, if all four character combinations are present, it is referred to as a character incompatibility. Since character incompatibility is deemed to be rare in populations with extensive asexual reproduction, a character incompatibility analysis was also performed in order to infer the predominant mating system in the introduced M. micrantha populations. Out of 483 AFLP loci examined using stringent significance criteria, 14 highly credible outlier loci were identified by Dfdist and Bayescan. Moreover, remarkable genetic variation, multiple introductions, substantial bottlenecks and character compatibility were found to occur in M. micrantha. Thus local adaptation at the genome level indeed exists in M. micrantha, and may represent a major evolutionary mechanism of successful invasion. Interactions between genetic diversity, multiple introductions, and reproductive modes contribute to increase the capacity of adaptive evolution.

Tissue sampling methods and standards for vertebrate genomics

The recent rise in speed and efficiency of new sequencing technologies have facilitated high-throughput sequencing, assembly and analyses of genomes, advancing ongoing efforts to analyze genetic sequences across major vertebrate groups. Standardized procedures in acquiring high quality DNA and RNA and establishing cell lines from target species will facilitate these initiatives. We provide a legal and methodological guide according to four standards of acquiring and storing tissue for the Genome 10K Project and similar initiatives as follows: four-star (banked tissue/cell cultures, RNA from multiple types of tissue for transcriptomes, and sufficient flash-frozen tissue for 1 mg of DNA, all from a single individual); three-star (RNA as above and frozen tissue for 1 mg of DNA); two-star (frozen tissue for at least 700 μg of DNA); and one-star (ethanol-preserved tissue for 700 μg of DNA or less of mixed quality). At a minimum, all tissues collected for the Genome 10K and other genomic projects should consider each species' natural history and follow institutional and legal requirements. Associated documentation should detail as much information as possible about provenance to ensure representative sampling and subsequent sequencing. Hopefully, the procedures outlined here will not only encourage success in the Genome 10K Project but also inspire the adaptation of standards by other genomic projects, including those involving other biota.

Incompatibility and competitive exclusion of genomic segments between sibling Drosophila species

The extent and nature of genetic incompatibilities between incipient races and sibling species is of fundamental importance to our view of speciation. However, with the exception of hybrid inviability and sterility factors, little is known about the extent of other, more subtle genetic incompatibilities between incipient species. Here we experimentally demonstrate the prevalence of such genetic incompatibilities between two young allopatric sibling species, Drosophila simulans and D. sechellia. Our experiments took advantage of 12 introgression lines that carried random introgressed D. sechellia segments in different parts of the D. simulans genome. First, we found that these introgression lines did not show any measurable sterility or inviability effects. To study if these sechellia introgressions in a simulans background contained other fitness consequences, we competed and genetically tracked the marked alleles within each introgression against the wild-type alleles for 20 generations. Strikingly, all marked D. sechellia introgression alleles rapidly decreased in frequency in only 6 to 7 generations. We then developed computer simulations to model our competition results. These simulations indicated that selection against D. sechellia introgression alleles was high (average s = 0.43) and that the marker alleles and the incompatible alleles did not separate in 78% of the introgressions. The latter result likely implies that most introgressions contain multiple genetic incompatibilities. Thus, this study reveals that, even at early stages of speciation, many parts of the genome diverge to a point where introducing foreign elements has detrimental fitness consequences, but which cannot be seen using standard sterility and inviability assays.

Determining the extent of genomic incompatibilities is a pivotal issue in understanding the process of speciation. A controversial topic that has recently sparked debate is whether there are few isolated genetic regions (so-called “genomic islands of speciation”) or extensive genetic regions (“genomic continents of speciation”) responsible for species divergence. To answer this question, most work has focused on species divergence at the DNA sequence level. Here, we present a new perspective by shifting the focus to the fitness and functional aspects of foreign genomic introgression. To illustrate our point, we performed an introgression experiment on two sibling species, D. sechellia and D. simulans. After introgressing random genomic segments of D. sechellia into D. simulans genetic background, a 20-generation competition experiment revealed that, even at the early stages of speciation, there are virtually always detrimental fitness consequences to introducing random foreign elements from one genome to another. This implies that incipient speciation may be characterized by widespread accumulation of genomic incompatibilities rather than a few isolated genes. This study shows that we should move beyond the sterility and inviability assays in order to understand the full extent of genetic incompatibilities during speciation.

Divergence hitchhiking and the spread of genomic isolation during ecological speciation-with-gene-flow

In allopatric populations, geographical separation simultaneously isolates the entire genome, allowing genetic divergence to accumulate virtually anywhere in the genome. In sympatric populations, however, the strong divergent selection required to overcome migration produces a genetic mosaic of divergent and non-divergent genomic regions. In some recent genome scans, each divergent genomic region has been interpreted as an independent incidence of migration/selection balance, such that the reduction of gene exchange is restricted to a few kilobases around each divergently selected gene. I propose an alternative mechanism, 'divergence hitchhiking' (DH), in which divergent selection can reduce gene exchange for several megabases around a gene under strong divergent selection. Not all genes/markers within a DH region are divergently selected, yet the entire region is protected to some degree from gene exchange, permitting genetic divergence from mechanisms other than divergent selection to accumulate secondarily. After contrasting DH and multilocus migration/selection balance (MM/SB), I outline a model in which genomic isolation at a given genomic location is jointly determined by DH and genome-wide effects of the progressive reduction in realized migration, then illustrate DH using data from several pairs of incipient species in the wild.