Widespread yet heterogeneous genomic divergence

Genetics and the Evolution of Prezygotic Isolation

The significance of prezygotic isolation for speciation has been recognized at least since the Modern Synthesis. However, fundamental questions remain. For example, how are genetic associations between traits that contribute to prezygotic isolation maintained? What is the source of genetic variation underlying the evolution of these traits? And how do prezygotic barriers affect patterns of gene flow? We address these questions by reviewing genetic features shared across plants and animals that influence prezygotic isolation. Emerging technologies increasingly enable the identification and functional characterization of the genes involved, allowing us to test established theoretical expectations. Embedding these genes in their developmental context will allow further predictions about what constrains the evolution of prezygotic isolation. Ongoing improvements in statistical and computational tools will reveal how pre- and postzygotic isolation may differ in how they influence gene flow across the genome. Finally, we highlight opportunities for progress by combining theory with appropriate data.

Compensatory Base Changes Reveal Sexual Incompatibility among Members of the Anopheles subpictus Sensu Lato (Diptera: Culicidae) Species Complex in Sri Lanka

The mosquito Anopheles (Cellia) subpictus sensu lato (s.l.) is a major secondary vector of malaria in Sri Lanka. The sibling species composition in this species complex in Sri Lanka remains debatable. Compensatory base changes (CBCs) in the secondary structures of internal transcribed spacer 2 (ITS2) are reliable sources to predict sexual incompatibility among closely related species. The objective of the present study was to investigate the An. subpictus s.l. populations in Sri Lanka using the CBC analysis. Mosquito DNA was amplified and sequenced for the ITS2 region. The sequences were annotated using ITS2 Database. ITS2 secondary structures were constructed and analyzed for CBCs using various bioinformatics tools. The ITS2 regions consisted of two different lengths, 575 bp and 480 bp. The two CBCs and three hemi CBCs identified in the present study suggest that there may be at least two sexually incompatible sibling species. In conclusion, it is likely that there may be only two reproductively isolated sibling species in the An. subpictus species complex in Sri Lanka. However, due to high divergence of ITS2 in these species, it is reasonable to assume that they may be undergoing a speciation event to separate as a distinct species.

Structure-behaviour correlations between two genetically closely related snail species

Species, through their structure and composition, have evolved to respond to environmental constraints. Predator-prey interactions are among environmental pressures that can lead to speciation, but it remains unclear how this pressure can be related to the material structure and performance. Recently, two land snails, Karaftohelix editha and Karaftohelix gainesi, were found to exhibit divergent phenotypes and responses to predation despite sharing the same habitat and most of their genome. Indeed, under attack from a beetle, K. editha snails retract into their shell whereas K. gainesi snails swing their shell. In this paper, we looked at the microstructure, composition, morphology and mechanics of the shells of those two species and discuss potential relationships between material structure and the snail defence behaviour. The results of this study provide additional arguments for the role of predator-prey interactions on speciation, as well as an unusual approach for the design of biomimetic structures adapted to a particular function.

Environmental drivers and genomic architecture of trait differentiation in fire-adapted Banksia attenuata ecotypes

Trait divergence between populations is considered an adaptive response to different environments, but to what extent this response is accompanied by genetic differentiation is less clear since it may be phenotypic plasticity. In this study, we analyzed phenotypic variation between two Banksia attenuata growth forms, lignotuberous (shrub) and epicormic resprouting (tree), in fire-prone environments to identify the environmental factors that have driven this phenotypic divergence. We linked genotype with phenotype and traced candidate genes using differential gene expression analysis. Fire intervals determined the phenotypic divergence between growth forms in B. attenuata. A genome-wide association study identified 69 single nucleotide polymorphisms, putatively associated with growth form, whereas no growth form- or phenotype-specific genotypes were identified. Genomic differentiation between the two growth forms was low (F_st  = 0.024). Differential gene expression analysis identified 37 genes/transcripts that were differentially expressed in the two growth forms. A small heat-shock protein gene, associated with lignotuber presence, was differentially expressed in the two forms. We conclude that different fire regimes induce phenotypic polymorphism in B. attenuata, whereas phenotypic trait divergence involves the differential expression of a small fraction of genes that interact strongly with the disturbance regime. Thus, phenotypic plasticity among resprouters is the general strategy for surviving varying fire regimes.

Don't throw out the sympatric speciation with the crater lake water: fine-scale investigation of introgression provides equivocal support for causal role of secondary gene flow in one of the clearest examples of sympatric speciation

Genomic data has revealed complex histories of colonization and repeated gene flow previously unrecognized in some of the most celebrated examples of sympatric speciation and radiation. However, much of the evidence for secondary gene flow into these radiations comes from summary statistics calculated from sparse genomic sampling without knowledge of which specific genomic regions introgressed. This tells us little about how gene flow potentially influenced sympatric diversification. Here, we investigated whole genomes of Barombi Mbo crater lake cichlids for fine-scale patterns of introgression with neighboring riverine cichlid populations. We found evidence of secondary gene flow into the radiation scattered across <0.24% of the genome; however, from our analyses, it is not clear if the functional diversity in these regions contributed to the ecological, sexual, and morphological diversity found in the lake. Unlike similar studies, we found no obvious candidate genes for adaptive introgression and we cannot rule out that secondary gene flow was predominantly neutral with respect to the diversification process. We also found evidence for differential assortment of ancestral polymorphisms found in riverine populations between sympatric sister species, suggesting the presence of an ancestral hybrid swarm. Although the history of gene flow and colonization is more complicated than previously assumed, the lack of compelling evidence for secondary gene flow's role in species diversification suggests that we should not yet rule out one of the most celebrated examples of sympatric speciation in nature without a more thorough investigation of the timing and functional role of each introgressed region.

Chromosome arm-specific patterns of polymorphism associated with chromosomal inversions in the major African malaria vector, Anopheles funestus

Chromosomal inversions facilitate local adaptation of beneficial mutations and modulate genetic polymorphism, but the extent of their effects within the genome is still insufficiently understood. The genome of Anopheles funestus, a malaria mosquito endemic to sub-Saharan Africa, contains an impressive number of paracentric polymorphic inversions, which are unevenly distributed among chromosomes and provide an excellent framework for investigating the genomic impacts of chromosomal rearrangements. Here, we present results of a fine-scale analysis of genetic variation within the genome of two weakly differentiated populations of Anopheles funestus inhabiting contrasting moisture conditions in Cameroon. Using population genomic analyses, we found that genetic divergence between the two populations is centred on regions of the genome corresponding to three inversions, which are characterized by high values of FST , absolute sequence divergence and fixed differences. Importantly, in contrast to the 2L chromosome arm, which is collinear, nucleotide diversity is significantly reduced along the entire length of three autosome arms bearing multiple overlapping chromosomal rearrangements. These findings support the idea that interactions between reduced recombination and natural selection within inversions contribute to sculpt nucleotide polymorphism across chromosomes in An. funestus.

Genomic insights into adaptive divergence and speciation among malaria vectors of the Anopheles nili group

Ongoing speciation in the most important African malaria vectors gives rise to cryptic populations, which differ remarkably in their behavior, ecology, and capacity to vector malaria parasites. Understanding the population structure and the drivers of genetic differentiation among mosquitoes is crucial for effective disease control because heterogeneity within vector species contributes to variability in malaria cases and allow fractions of populations to escape control efforts. To examine population structure and the potential impacts of recent large-scale control interventions, we have investigated the genomic patterns of differentiation in mosquitoes belonging to the Anopheles nili group-a large taxonomic group that diverged ~3 Myr ago. Using 4,343 single nucleotide polymorphisms (SNPs), we detected strong population structure characterized by high-FST values between multiple divergent populations adapted to different habitats within the Central African rainforest. Delineating the cryptic species within the Anopheles nili group is challenging due to incongruence between morphology, ribosomal DNA, and SNP markers consistent with incomplete lineage sorting and/or interspecific gene flow. A very high proportion of loci are fixed (FST = 1) within the genome of putative species, which suggests that ecological and/or reproductive barriers are maintained by strong selection on a substantial number of genes.

A history of hybrids? Genomic patterns of introgression in the True Geese

Background

The impacts of hybridization on the process of speciation are manifold, leading to distinct patterns across the genome. Genetic differentiation accumulates in certain genomic regions, while divergence is hampered in other regions by homogenizing gene flow, resulting in a heterogeneous genomic landscape. A consequence of this heterogeneity is that genomes are mosaics of different gene histories that can be compared to unravel complex speciation and hybridization events. However, incomplete lineage sorting (often the outcome of rapid speciation) can result in similar patterns. New statistical techniques, such as the D-statistic and hybridization networks, can be applied to disentangle the contributions of hybridization and incomplete lineage sorting. We unravel patterns of hybridization and incomplete lineage sorting during and after the diversification of the True Geese (family Anatidae, tribe Anserini, genera Anser and Branta) using an exon-based hybridization network approach and taking advantage of discordant gene tree histories by re-sequencing all taxa of this clade. In addition, we determine the timing of introgression and reconstruct historical effective population sizes for all goose species to infer which demographic or biogeographic factors might explain the observed patterns of introgression.

Results

We find indications for ancient interspecific gene flow during the diversification of the True Geese and were able to pinpoint several putative hybridization events. Specifically, in the genus Branta, both the ancestor of the White-cheeked Geese (Hawaiian Goose, Canada Goose, Cackling Goose and Barnacle Goose) and the ancestor of the Brent Goose hybridized with Red-breasted Goose. One hybridization network suggests a hybrid origin for the Red-breasted Goose, but this scenario seems unlikely and it not supported by the D-statistic analysis. The complex, highly reticulated evolutionary history of the genus Anser hampered the estimation of ancient hybridization events by means of hybridization networks. The reconstruction of historical effective population sizes shows that most species showed a steady increase during the Pliocene and Pleistocene. These large effective population sizes might have facilitated contact between diverging goose species, resulting in the establishment of hybrid zones and consequent gene flow.

Conclusions

Our analyses suggest that the evolutionary history of the True Geese is influenced by introgressive hybridization. The approach that we have used, based on genome-wide phylogenetic incongruence and network analyses, will be a useful procedure to reconstruct the complex evolutionary histories of many naturally hybridizing species groups.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1048-2) contains supplementary material, which is available to authorized users.

Pollutants and Insecticides Drive Local Adaptation in African Malaria Mosquitoes

The Anopheles gambiae complex contains a number of highly anthropophilic mosquito species that have acquired exceptional ability to thrive in complex human habitats. Thus, examining the evolutionary history of this Afrotropical mosquito may yield vital information on the selective processes that occurred during the adaptation to human-dominated environments. We performed reduced representation sequencing on 941 mosquitoes of the Anopheles gambiae complex collected across four ecogeographic zones in Cameroon. We find evidence for genetic and geographic subdivision within An. coluzzii and An. gambiae sensu stricto-the two most significant malaria vectors in the region. Importantly, in both species, rural and urban populations are genetically differentiated. Genome scans reveal pervasive signatures of selection centered on genes involved in xenobiotic resistance. Notably, a selective sweep containing detoxification enzymes is prominent in urban mosquitoes that exploit polluted breeding sites. Overall, our study suggests that recent anthropogenic environmental modifications and widespread use of insecticides are driving population differentiation and local adaptation in vectors with potentially significant consequences for malaria epidemiology.

Evolution of stickleback in 50 years on earthquake-uplifted islands

How rapidly can animal populations in the wild evolve when faced with sudden environmental shifts? Uplift during the 1964 Great Alaska Earthquake abruptly created freshwater ponds on multiple islands in Prince William Sound and the Gulf of Alaska. In the short time since the earthquake, the phenotypes of resident freshwater threespine stickleback fish on at least three of these islands have changed dramatically from their oceanic ancestors. To test the hypothesis that these freshwater populations were derived from oceanic ancestors only 50 y ago, we generated over 130,000 single-nucleotide polymorphism genotypes from more than 1,000 individuals using restriction site-associated DNA sequencing (RAD-seq). Population genomic analyses of these data support the hypothesis of recent and repeated, independent colonization of freshwater habitats by oceanic ancestors. We find evidence of recurrent gene flow between oceanic and freshwater ecotypes where they co-occur. Our data implicate natural selection in phenotypic diversification and support the hypothesis that the metapopulation organization of this species helps maintain a large pool of genetic variation that can be redeployed rapidly when oceanic stickleback colonize freshwater environments. We find that the freshwater populations, despite population genetic analyses clearly supporting their young age, have diverged phenotypically from oceanic ancestors to nearly the same extent as populations that were likely founded thousands of years ago. Our results support the intriguing hypothesis that most stickleback evolution in fresh water occurs within the first few decades after invasion of a novel environment.

Can we continue to neglect genomic variation in introgression rates when inferring the history of speciation? A case study in a Mytilus hybrid zone

The use of molecular data to reconstruct the history of divergence and gene flow between populations of closely related taxa represents a challenging problem. It has been proposed that the long-standing debate about the geography of speciation can be resolved by comparing the likelihoods of a model of isolation with migration and a model of secondary contact. However, data are commonly only fit to a model of isolation with migration and rarely tested against the secondary contact alternative. Furthermore, most demographic inference methods have neglected variation in introgression rates and assume that the gene flow parameter (Nm) is similar among loci. Here, we show that neglecting this source of variation can give misleading results. We analysed DNA sequences sampled from populations of the marine mussels, Mytilus edulis and M. galloprovincialis, across a well-studied mosaic hybrid zone in Europe and evaluated various scenarios of speciation, with or without variation in introgression rates, using an Approximate Bayesian Computation (ABC) approach. Models with heterogeneous gene flow across loci always outperformed models assuming equal migration rates irrespective of the history of gene flow being considered. By incorporating this heterogeneity, the best-supported scenario was a long period of allopatric isolation during the first three-quarters of the time since divergence followed by secondary contact and introgression during the last quarter. By contrast, constraining migration to be homogeneous failed to discriminate among any of the different models of gene flow tested. Our simulations thus provide statistical support for the secondary contact scenario in the European Mytilus hybrid zone that the standard coalescent approach failed to confirm. Our results demonstrate that genomic variation in introgression rates can have profound impacts on the biological conclusions drawn from inference methods and needs to be incorporated in future studies.

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.

Heterogeneous Genomic Differentiation in marine threespine sticklebacks: adaptation along an environmental gradient

Evolutionary divergence among populations occupying ecologically distinct environments can occur even in the face of on-going gene flow. However, the genetic underpinnings, as well as the scale and magnitude at which this differentiation occurs in marine habitats are not well understood. We investigated the patterns and degree of genomic heterogeneity in threespine sticklebacks (Gasterosteus aculeatus) by assessing genetic variability in 20 nongenic and 20 genic (associated with genes important for freshwater adaptation) microsatellite loci in samples collected from 38 locations spanning the entire Baltic Sea coast to the North Sea boundary. Population divergence (F(ST) ≈ 0.026) and structuring (five genetic clusters) was significantly more pronounced in the genic as compared to nongenic markers (F(ST) ≈ 0.008; no genetic clusters). Patterns of divergence in the genic markers--45% of which were identified as outliers--correlated with local differences in salinity. Yet, a strong positive correlation between divergence in genic and nongenic markers, and their association with environmental factors suggests that adaptive divergence is reducing gene flow across the genome. Apart from providing a clear demonstration of heterogeneous genomic patterns of differentiation in a marine species, the results are indicative of adaptive population structuring across the relatively young Baltic Sea in spite of ample opportunities for gene flow.

Genetic hitchhiking and the dynamic buildup of genomic divergence during speciation with gene flow

A major issue in evolutionary biology is explaining patterns of differentiation observed in population genomic data, as divergence can be due to both direct selection on a locus and genetic hitchhiking. "Divergence hitchhiking" (DH) theory postulates that divergent selection on a locus reduces gene flow at physically linked sites, facilitating the formation of localized clusters of tightly linked, diverged loci. "Genome hitchhiking" (GH) theory emphasizes genome-wide effects of divergent selection. Past theoretical investigations of DH and GH focused on static snapshots of divergence. Here, we used simulations assessing a variety of strengths of selection, migration rates, population sizes, and mutation rates to investigate the relative importance of direct selection, GH, and DH in facilitating the dynamic buildup of genomic divergence as speciation proceeds through time. When divergently selected mutations were limiting, GH promoted divergence, but DH had little measurable effect. When populations were small and divergently selected mutations were common, DH enhanced the accumulation of weakly selected mutations, but this contributed little to reproductive isolation. In general, GH promoted reproductive isolation by reducing effective migration rates below that due to direct selection alone, and was important for genome-wide "congealing" or "coupling" of differentiation (F(ST)) across loci as speciation progressed.