Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.

Rapid sympatric ecological differentiation of crater lake cichlid fishes within historic times

Background: After a volcano erupts, a lake may form in the cooled crater and become an isolated aquatic ecosystem. This makes fishes in crater lakes informative for understanding sympatric evolution and ecological diversification in barren environments. From a geological and limnological perspective, such research offers insight about the process of crater lake ecosystem establishment and speciation. In the present study we use genetic and coalescence approaches to infer the colonization history of Midas cichlid fishes (Amphilophus cf. citrinellus) that inhabit a very young crater lake in Nicaragua-the ca. 1800 year-old Lake Apoyeque. This lake holds two sympatric, endemic morphs of Midas cichlid: one with large, hypertrophied lips (~20% of the total population) and another with thin lips. Here we test the associated ecological, morphological and genetic diversification of these two morphs and their potential to represent incipient speciation.
Results: Gene coalescence analyses [11 microsatellite loci and mitochondrial DNA (mtDNA) sequences] suggest that crater lake Apoyeque was colonized in a single event from the large neighbouring great lake Managua only about 100 years ago. This founding in historic times is also reflected in the extremely low nuclear and mitochondrial genetic diversity in Apoyeque. We found that sympatric adult thin- and thick-lipped fishes occupy distinct ecological trophic niches. Diet, body shape, head width, pharyngeal jaw size and shape and stable isotope values all differ significantly between the two lip-morphs. The eco-morphological features pharyngeal jaw shape, body shape, stomach contents and stable isotopes (d15N) all show a bimodal distribution of traits, which is compatible with the expectations of an initial stage of ecological speciation under disruptive selection. Genetic differentiation between the thin- and thick-lipped population is weak at mtDNA sequence (FST = 0.018) and absent at nuclear microsatellite loci (FST < 0.001).
Conclusions: This study provides empirical evidence of eco-morphological differentiation occurring very quickly after the colonization of a new and vacant habitat. Exceptionally low levels of neutral genetic diversity and inference from coalescence indicates that the Midas cichlid population in Apoyeque is much younger (ca. 100 years or generations old) than the crater itself (ca. 1 800 years old). This suggests either that the crater remained empty for many hundreds of years after its formation or that remnant volcanic activity prevented the establishment of a stable fish population during the early life of the crater lake. Based on our findings of eco-morphological variation in the Apoyeque Midas cichlids, and known patterns of adaptation in Midas cichlids in general, we suggest that this population may be in a very early stage of speciation (incipient species), promoted by disruptive selection and ecological diversification.

Convergent evolutionary processes driven by foraging opportunity in two sympatric morph pairs of Arctic charr with contrasting post-glacial origins

The expression of two or more discrete phenotypes amongst individuals within a species (morphs) provides multiple modes upon which selection can act semi-independently, and thus may be an important stage in speciation. In the present study, we compared two sympatric morph systems aiming to address hypotheses related to their evolutionary origin. Arctic charr in sympatry in Loch Tay, Scotland, exhibit one of two discrete, alternative body size phenotypes at maturity (large or small body size). Arctic charr in Loch Awe segregate into two temporally segregated spawning groups (breeding in either spring or autumn). Mitochondrial DNA restriction fragment length polymorphism analysis showed that the morph pairs in both lakes comprise separate gene pools, although segregation of the Loch Awe morphs is more subtle than that of Loch Tay. We conclude that the Loch Awe morphs diverged in situ (within the lake), whereas Loch Tay morphs most likely arose through multiple invasions by different ancestral groups that segregated before post-glacial invasion (i.e. in allopatry). Both morph pairs showed clear trophic segregation between planktonic and benthic resources (measured by stable isotope analysis) but this was significantly less distinct in Loch Tay than in Loch Awe. By contrast, both inter-morph morphological and life-history differences were more subtle in Loch Awe than in Loch Tay. The strong ecological but relatively weak morphological and life-history divergence of the in situ derived morphs compared to morphs with allopatric origins indicates a strong link between early ecological and subsequent genetic divergence of sympatric origin emerging species pairs. The emergence of parallel specialisms despite distinct genetic origins of these morph pairs suggests that the effect of available foraging opportunities may be at least as important as genetic origin in structuring sympatric divergence in post-glacial fishes with high levels of phenotypic plasticity. (c) 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, , .

Genome-wide evidence for speciation with gene flow in Heliconius butterflies

Most speciation events probably occur gradually, without complete and immediate reproductive isolation, but the full extent of gene flow between diverging species has rarely been characterized on a genome-wide scale. Documenting the extent and timing of admixture between diverging species can clarify the role of geographic isolation in speciation. Here we use new methodology to quantify admixture at different stages of divergence in Heliconius butterflies, based on whole-genome sequences of 31 individuals. Comparisons between sympatric and allopatric populations of H. melpomene, H. cydno, and H. timareta revealed a genome-wide trend of increased shared variation in sympatry, indicative of pervasive interspecific gene flow. Up to 40% of 100-kb genomic windows clustered by geography rather than by species, demonstrating that a very substantial fraction of the genome has been shared between sympatric species. Analyses of genetic variation shared over different time intervals suggested that admixture between these species has continued since early in speciation. Alleles shared between species during recent time intervals displayed higher levels of linkage disequilibrium than those shared over longer time intervals, suggesting that this admixture took place at multiple points during divergence and is probably ongoing. The signal of admixture was significantly reduced around loci controlling divergent wing patterns, as well as throughout the Z chromosome, consistent with strong selection for Müllerian mimicry and with known Z-linked hybrid incompatibility. Overall these results show that species divergence can occur in the face of persistent and genome-wide admixture over long periods of time.

Adaptive divergence and speciation among sexual and pseudoviviparous populations of Festuca

Pseudovivipary is an environmentally induced flowering abnormality in which vegetative shoots replace seminiferous (sexual) inflorescences. Pseudovivipary is usually retained in transplantation experiments, indicating that the trait is not solely induced by the growing environment. Pseudovivipary is the defining characteristic of Festuca vivipara, and arguably the only feature separating this species from its closest seminiferous relative, Festuca ovina. We performed phylogenetic and population genetic analysis on sympatric F. ovina and F. vivipara samples to establish whether pseudovivipary is an adaptive trait that accurately defines the separation of genetically distinct Festuca species. Chloroplast and nuclear marker-based analyses revealed that variation at a geographical level can exceed that between F. vivipara and F. ovina. We deduced that F. vivipara is a recent species that frequently arises independently within F. ovina populations and has not accumulated significant genetic differentiation from its progenitor. We inferred local gene flow between the species. We identified one amplified fragment length polymorphism marker that may be linked to a pseudovivipary-related region of the genome, and several other markers provide evidence of regional local adaptation in Festuca populations. We conclude that F. vivipara can only be appropriately recognized as a morphologically and ecologically distinct species; it lacks genetic differentiation from its relatives. This is the first report of a ‘failure in normal flowering development’ that repeatedly appears to be adaptive, such that the trait responsible for species recognition constantly reappears on a local basis.

Cryptic speciation and recombination in the fungus Paracoccidioides brasiliensis as revealed by gene genealogies

Paracoccidioides brasiliensis is the etiologic agent of paracoccidioidomycosis, a disease confined to Latin America and of marked importance in the endemic areas due to its frequency and severity. This species is considered to be clonal according to mycological criteria and has been shown to vary in virulence. To characterize natural genetic variation and reproductive mode in this fungus, we analyzed P. brasiliensis phylogenetically in search of cryptic species and possible recombination using concordance and nondiscordance of gene genealogies with respect to phylogenies of eight regions in five nuclear loci. Our data indicate that this fungus consists of at least three distinct, previously unrecognized species: S1 (species 1 with 38 isolates), PS2 (phylogenetic species 2 with six isolates), and PS3 (phylogenetic species 3 with 21 isolates). Genealogies of four of the regions studied strongly supported the PS2 clade, composed of five Brazilian and one Venezuelan isolate. The second clade, PS3, composed solely of 21 Colombian isolates, was strongly supported by the alpha-tubulin genealogy. The remaining 38 individuals formed S1. Two of the three lineages of P. brasiliensis, S1 and PS2, are sympatric across their range, suggesting barriers to gene flow other than geographic isolation. Our study provides the first evidence for possible sexual reproduction in P. brasiliensis S1, but does not rule it out in the other two species.

Interploidy hybridization in sympatric zones: the formation of Epidendrum fulgens x E. puniceoluteum hybrids (Epidendroideae, Orchidaceae)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Mating system variation and assortative mating of sympatric bromeliads (Pitcairnia spp.) endemic to neotropical inselbergs

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Community Genetics Reveal Elevated Levels of Sympatric Gene Flow among Morphologically Similar but Not among Morphologically Dissimilar Species of Lake Victoria Cichlid Fish

We examined genetic structure among five species of Lake Victoria haplochromine cichlids in four island communities, using a full factorial sampling design that compared genetic differentiation between pairs of species and populations of varying morphological similarity and geographical proximity. We found that allopatric conspecific populations were on average significantly more strongly differentiated than sympatric heterospecific populations of morphologically similar species. Allopatric heterospecific populations of morphologically dissimilar species were most differentiated. Our work demonstrates that phenotypic divergence can be maintained and perhaps even evolve in sympatry despite considerable gene flow between species. Conversely, phenotypic resemblance among conspecific populations can be maintained despite geographical isolation. Additionally we show that anthropogenically increased hybridization does not affect all sympatric species evenly but predominantly affects morphologically similar and closely related species. This has important implications for the evolution of reproductive isolation between species These findings are also consistent with the hypothesis of speciation reversal due to weakening of divergent selection and reproductive isolation as a consequence of habitat homogenization and offers an evolutionary mechanistic explanation for the observation that species poor assemblages in turbid areas of the lake are characterized by just one or two species in each of a few morphologically distinct genera.

Speciation reversal and biodiversity dynamics with hybridization in changing environments

A considerable fraction of the world's biodiversity is of recent evolutionary origin and has evolved as a by-product of, and is maintained by, divergent adaptation in heterogeneous environments. Conservationists have paid attention to genetic homogenization caused by human-induced translocations (e.g. biological invasions and stocking), and to the importance of environmental heterogeneity for the ecological coexistence of species. However, far less attention has been paid to the consequences of loss of environmental heterogeneity to the genetic coexistence of sympatric species. Our review of empirical observations and our theoretical considerations on the causes and consequences of interspecific hybridization suggest that a loss of environmental heterogeneity causes a loss of biodiversity through increased genetic admixture, effectively reversing speciation. Loss of heterogeneity relaxes divergent selection and removes ecological barriers to gene flow between divergently adapted species, promoting interspecific introgressive hybridization. Since heterogeneity of natural environments is rapidly deteriorating in most biomes, the evolutionary ecology of speciation reversal ought to be fully integrated into conservation biology.

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 FST 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 FST 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.

Divergent selection during speciation of Lake Malawi cichlid fishes inferred from parallel radiations in nuptial coloration

Repeated evolution of the same phenotypic difference during independent episodes of speciation is strong evidence for selection during speciation. More than 1,000 species of cichlids, >10% of the world's freshwater fish species, have arisen within the past million years in Lakes Malawi and Victoria in eastern Africa. Many pairs of closely related sympatric species differ in their nuptial coloration in very similar ways. Nuptial coloration is important in their mate choice, and speciation by sexual selection on genetically or ecologically constrained variation in nuptial coloration had been proposed, which would repeatedly produce similar nuptial types in different populations, a prediction that was difficult to test in the absence of population-level phylogenies. We measured genetic similarity between individuals within and between populations, species, and lake regions by typing 59 individuals at >2,000 polymorphic genetic loci. From these data, we reconstructed, to our knowledge, the first larger species level phylogeny for the most diverse group of Lake Malawi cichlids. We used the genetic and phylogenetic data to test the divergent selection scenario against colonization, character displacement, and hybridization scenarios that could also explain diverse communities. Diversity has arisen by replicated radiations into the same color types, resulting in phenotypically very different, yet closely related, species within and phenotypically highly similar yet unrelated sets of species between regions, which is consistent with divergent selection during speciation and is inconsistent with colonization and character displacement models.

Ecological speciation and phenotypic plasticity affect ecosystems

Phenotypic differences among closely related populations and species can cause contrasting effects on ecosystems; however, it is unknown whether such effects result from genetic divergence, phenotypic plasticity, or both. To test this, we reared sympatric limnetic and benthic species of whitefish from a young adaptive radiation in a common garden, where the benthic species was raised on two distinct food types. We then used these fish in a mesocosm experiment to test for contrasting ecosystem effects of closely related species and of plastically induced differences within a species. We found that strong contrasting ecosystem effects resulted more frequently from genetic divergence, although they were not stronger overall than those resulting from phenotypic plasticity. Overall, our results provide evidence that genetically based differences among closely related species that evolved during a young adaptive radiation can affect ecosystems, and that phenotypic plasticity can modify the ecosystem effects of such species.

Direct mate choice maintains diversity among sympatric cichlids in Lake Victoria

Mate choice may play an important role in animal speciation. The haplochromine cichlids of Lake Victoria are suitable to test this hypothesis. Diversity in ecology, coloration and anatomy evolved in these fish faster than postzygotic barriers to gene flow, and little is known about how this diversity is maintained. It was tested whether recognizable forms are selection-maintained morphs or reproductively isolated species by investigating in the field reproductive timing, location of spawning sites, and mate choice behaviour. There was a large interspecific overlap in timing of breeding and location of spawning sites, which was largest in members of the same genus. Behavioural mate choice of such closely related taxa was highly assortative, such that it is likely that they are sexually isolated species and that direct mate choice is the major force that directs gene flow and maintains form diversity. The results differ from what is known about recent radiations of other lacustrine fish groups where speciation seems to be driven by diverging microhabitat preferences or diverging timing of reproduction, but are in agreement with predictions from models of speciation by diverging mate preferences.

Speciation of Antechinus stuartii and A. subtropicus (Marsupialia : Dasyuridae) in eastern Australia: molecular and morphological evidence

This paper evaluates the systematic status of the Antechinus populations of northern New South Wales and southern Queensland, using a combined morphological and molecular (allozymes and mitochondrial DNA) approach. Analysis of the d-loop section of the mitochondrial DNA control region revealed two highly supported clades within A. stuartii sensu lato that were sympatric in the Border Ranges of northern New South Wales. However, genetic distances between these clades were small ( approximately 3%), indicating that time of divergence was probably relatively recent. Allozyme electrophoresis also showed very small differences between clades/ species. Analyses of cranial and dental characters showed that the members of each of these clades differed morphologically and that the clades corresponded to A. stuartii and the recently described A. subtropicus. The combined results support the species status of A. stuartii and A. subtropicus, and suggest that speciation was likely a result of a recent vicariant event.