The coincidence of ecological opportunity with hybridization explains rapid adaptive radiation in Lake Mweru cichlid fishes

The process of adaptive radiation was classically hypothesized to require isolation of a lineage from its source (no gene flow) and from related species (no competition). Alternatively, hybridization between species may generate genetic variation that facilitates adaptive radiation. Here we study haplochromine cichlid assemblages in two African Great Lakes to test these hypotheses. Greater biotic isolation (fewer lineages) predicts fewer constraints by competition and hence more ecological opportunity in Lake Bangweulu, whereas opportunity for hybridization predicts increased genetic potential in Lake Mweru. In Lake Bangweulu, we find no evidence for hybridization but also no adaptive radiation. We show that the Bangweulu lineages also colonized Lake Mweru, where they hybridized with Congolese lineages and then underwent multiple adaptive radiations that are strikingly complementary in ecology and morphology. Our data suggest that the presence of several related lineages does not necessarily prevent adaptive radiation, although it constrains the trajectories of morphological diversification. It might instead facilitate adaptive radiation when hybridization generates genetic variation, without which radiation may start much later, progress more slowly or never occur.

Genomic architecture of parallel ecological divergence: Beyond a single environmental contrast

The study of parallel ecological divergence provides important clues to the operation of natural selection. Parallel divergence often occurs in heterogeneous environments with different kinds of environmental gradients in different locations, but the genomic basis underlying this process is unknown. We investigated the genomics of rapid parallel adaptation in the marine snail Littorina saxatilis in response to two independent environmental axes (crab-predation versus wave-action and low-shore versus high-shore). Using pooled whole-genome resequencing, we show that sharing of genomic regions of high differentiation between environments is generally low but increases at smaller spatial scales. We identify different shared genomic regions of divergence for each environmental axis and show that most of these regions overlap with candidate chromosomal inversions. Several inversion regions are divergent and polymorphic across many localities. We argue that chromosomal inversions could store shared variation that fuels rapid parallel adaptation to heterogeneous environments, possibly as balanced polymorphism shared by adaptive gene flow.

Admixture between old lineages facilitated contemporary ecological speciation in Lake Constance stickleback

Ecological speciation can sometimes rapidly generate reproductively isolated populations coexisting in sympatry, but the origin of genetic variation permitting this is rarely known. We previously explored the genomics of very recent ecological speciation into lake and stream ecotypes in stickleback from Lake Constance. Here, we reconstruct the origin of alleles underlying ecological speciation by combining demographic modelling on genome-wide single nucleotide polymorphisms, phenotypic data and mitochondrial sequence data in the wider European biogeographical context. We find that parallel differentiation between lake and stream ecotypes across replicate lake-stream ecotones resulted from recent secondary contact and admixture between old East and West European lineages. Unexpectedly, West European alleles that introgressed across the hybrid zone at the western end of the lake, were recruited to genomic islands of differentiation between ecotypes at the eastern end of the lake. Our results highlight an overlooked outcome of secondary contact: ecological speciation facilitated by admixture variation.

Ecological speciation can proceed rapidly, but the origin of genetic variation facilitating it has remained elusive. Here, the authors show that secondary contact and introgression between deeply diverged lineages of stickleback fish facilitated rapid ecological speciation into lake and stream ecotypes in Lake Constance.

Geographic variation in opsin expression does not align with opsin genotype in Lake Victoria cichlid populations

Sensory adaptation to the local environment can contribute to speciation. Aquatic environments are well suited for studying this process: The natural attenuation of light through water results in heterogeneous light environments, to which vision-dependent species must adapt for communication and survival. Here, we study visual adaptation in sympatric Pundamilia cichlids from southeastern Lake Victoria. Species with blue or red male nuptial coloration co-occur at many rocky islands but tend to be depth-differentiated, entailing different visual habitats, more strongly at some islands than others. Divergent visual adaptation to these environments has been implicated as a major factor in the divergence of P. pundamilia and P. nyererei, as they show consistent differentiation in the long-wavelength-sensitive visual pigment gene sequence (LWS opsin). In addition to sequence variation, variation in the opsin gene expression levels may contribute to visual adaptation. We characterized opsin gene expression and LWS genotype across Pundamilia populations inhabiting turbid and clear waters, to examine how different mechanisms of visual tuning contribute to visual adaptation. As predicted, the short-wavelength-sensitive opsin (SWS2b) was expressed exclusively in a population from clear water. Contrary to prediction however, expression levels of the other opsins were species- and island-dependent and did not align with species differences in LWS genotype. Specifically, in two locations with turbid water, the shallow-water dwelling blue species expressed more LWS and less RH2A than the deeper-dwelling red species, while the opposite pattern occurred in the two locations with clear water. Visual modeling suggests that the observed distribution of opsin expression profiles and LWS genotypes does not maximize visual performance, implying the involvement of additional visual tuning mechanisms and/or incomplete adaptation.

OPEN RESEARCH BADGE

This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://hdl.handle.net/10411/I1IUUQ.

Admixture between ancient lineages, selection, and the formation of sympatric stickleback species-pairs

Ecological speciation has become a popular model for the development and maintenance of reproductive isolation in closely related sympatric pairs of species or ecotypes. An implicit assumption has been that such pairs originate (possibly with gene flow) from a recent, genetically homogeneous ancestor. However, recent genomic data has revealed that currently sympatric taxa are often a result of secondary contact between ancestrally allopatric lineages. This has sparked an interest in the importance of initial hybridization upon secondary contact, with genomic re-analysis of classic examples of ecological speciation often implicating admixture in speciation. We describe a novel occurrence of unusually well-developed reproductive isolation in a model system for ecological speciation: the three-spined stickleback (Gasterosteus aculeatus), breeding sympatrically in multiple lagoons on the Scottish island of North Uist. Using morphological data, targeted genotyping and genome-wide single nucleotide polymorphism (SNP) data we show that lagoon resident and anadromous ecotypes are strongly reproductively isolated with an estimated hybridization rate of only ∼1%. We use palaeoecological and genetic data to test three hypotheses to explain the existence of these species-pairs. Our results suggest that recent, purely ecological speciation from a genetically homogeneous ancestor is probably not solely responsible for the evolution of species-pairs. Instead we reveal a complex colonisation history with multiple ancestral lineages contributing to the genetic composition of species-pairs, alongside strong disruptive selection. Our results imply a role for admixture upon secondary contact and are consistent with the recent suggestion that the genomic underpinning of ecological speciation often has an older, allopatric origin.

What can be learned by scanning the genome for molecular convergence in wild populations?

Convergent evolution, where independent lineages evolve similar phenotypes in response to similar challenges, can provide valuable insight into how selection operates and the limitations it encounters. However, it has only recently become possible to explore how convergent evolution is reflected at the genomic level. The overlapping outlier approach (OOA), where genome scans of multiple independent lineages are used to find outliers that overlap and therefore identify convergently evolving loci, is becoming popular. Here, we present a quantitative analysis of 34 studies that used this approach across many sampling designs, taxa, and sampling intensities. We found that OOA studies with increased biological sampling power within replicates have increased likelihood of finding overlapping, "convergent" signals of adaptation between them. When identifying convergent loci as overlapping outliers, it is tempting to assume that any false-positive outliers derived from individual scans will fail to overlap across replicates, but this cannot be guaranteed. We highlight how population demographics and genomic context can contribute toward both true convergence and false positives in OOA studies. We finish with an exploration of emerging methods that couple genome scans with phenotype and environmental measures, leveraging added information from genome data to more directly test hypotheses of the likelihood of convergent evolution.

Evidence that DNA repair genes, a family of tumor suppressor genes, are associated with evolution rate and size of genomes

Adaptive radiation and evolutionary stasis are characterized by very different evolution rates. The main aim of this study was to investigate if any genes have a special role to a high or low evolution rate. The availability of animal genomes permitted comparison of gene content of genomes of 24 vertebrate species that evolved through adaptive radiation (representing high evolutionary rate) and of 20 vertebrate species that are considered as living fossils (representing a slow evolutionary rate or evolutionary stasis). Mammals, birds, reptiles, and bony fishes were included in the analysis. Pathway analysis was performed for genes found to be specific in adaptive radiation or evolutionary stasis respectively. Pathway analysis revealed that DNA repair and cellular response to DNA damage are important (false discovery rate = 8.35 × 10⁻⁵; 7.15 × 10⁻⁶, respectively) for species evolved through adaptive radiation. This was confirmed by further genetic in silico analysis (p = 5.30 × 10⁻³). Nucleotide excision repair and base excision repair were the most significant pathways. Additionally, the number of DNA repair genes was found to be linearly related to the genome size and the protein number (proteome) of the 44 animals analyzed (p < 1.00 × 10⁻⁴), this being compatible with Drake’s rule. This is the first study where radiated and living fossil species have been genetically compared. Evidence has been found that cancer-related genes have a special role in radiated species. Linear association of the number of DNA repair genes with the species genome size has also been revealed. These comparative genetics results can support the idea of punctuated equilibrium evolution.

Disentangling structural genomic and behavioural barriers in a sea of connectivity

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

Inter-chromosomal coupling between vision and pigmentation genes during genomic divergence

Recombination between loci underlying mate choice and ecological traits is a major evolutionary force acting against speciation with gene flow. The evolution of linkage disequilibrium between such loci is therefore a fundamental step in the origin of species. Here, we show that this process can take place in the absence of physical linkage in hamlets-a group of closely related reef fishes from the wider Caribbean that differ essentially in colour pattern and are reproductively isolated through strong visually-based assortative mating. Using full-genome analysis, we identify four narrow genomic intervals that are consistently differentiated among sympatric species in a backdrop of extremely low genomic divergence. These four intervals include genes involved in pigmentation (sox10), axial patterning (hoxc13a), photoreceptor development (casz1) and visual sensitivity (SWS and LWS opsins) that develop islands of long-distance and inter-chromosomal linkage disequilibrium as species diverge. The relatively simple genomic architecture of species differences facilitates the evolution of linkage disequilibrium in the presence of gene flow.

Patterns of genomic differentiation between two Lake Victoria cichlid species, Haplochromis pyrrhocephalus and H. sp. 'macula'

BACKGROUND

The molecular basis of the incipient stage of speciation is still poorly understood. Cichlid fish species in Lake Victoria are a prime example of recent speciation events and a suitable system to study the adaptation and reproductive isolation of species.

RESULTS

Here, we report the pattern of genomic differentiation between two Lake Victoria cichlid species collected in sympatry, Haplochromis pyrrhocephalus and H. sp. 'macula,' based on the pooled genome sequences of 20 individuals of each species. Despite their ecological differences, population genomics analyses demonstrate that the two species are very close to a single panmictic population due to extensive gene flow. However, we identified 21 highly differentiated short genomic regions with fixed nucleotide differences. At least 15 of these regions contained genes with predicted roles in adaptation and reproductive isolation, such as visual adaptation, circadian clock, developmental processes, adaptation to hypoxia, and sexual selection. The nonsynonymous fixed differences in one of these genes, LWS, were reported as substitutions causing shift in absorption spectra of LWS pigments. Fixed differences were found in the promoter regions of four other differentially expressed genes, indicating that these substitutions may alter gene expression levels.

CONCLUSIONS

These diverged short genomic regions may have contributed to the differentiation of two ecologically different species. Moreover, the origins of adaptive variants within the differentiated regions predate the geological formation of Lake Victoria; thus Lake Victoria cichlid species diversified via selection on standing genetic variation.

Visual adaptation and microhabitat choice in Lake Victoria cichlid fish

When different genotypes choose different habitats to better match their phenotypes, genetic differentiation within a population may be promoted. Mating within those habitats may subsequently contribute to reproductive isolation. In cichlid fish, visual adaptation to alternative visual environments is hypothesized to contribute to speciation. Here, we investigated whether variation in visual sensitivity causes different visual habitat preferences, using two closely related cichlid species that occur at different but overlapping water depths in Lake Victoria and that differ in visual perception (Pundamilia spp.). In addition to species differences, we explored potential effects of visual plasticity, by rearing fish in two different light conditions: broad-spectrum (mimicking shallow water) and red-shifted (mimicking deeper waters). Contrary to expectations, fish did not prefer the light environment that mimicked their typical natural habitat. Instead, we found an overall preference for the broad-spectrum environment. We also found a transient influence of the rearing condition, indicating that the assessment of microhabitat preference requires repeated testing to control for familiarity effects. Together, our results show that cichlid fish exert visual habitat preference but do not support straightforward visual habitat matching.

Genomic insights into the vulnerability of sympatric whitefish species flocks

The erosion of habitat heterogeneity can reduce species diversity directly but can also lead to the loss of distinctiveness of sympatric species through speciation reversal. We know little about changes in genomic differentiation during the early stages of these processes, which can be mediated by anthropogenic perturbation. Here, we analyse three sympatric whitefish species (Coregonus spp) sampled across two neighbouring and connected Swiss pre-alpine lakes, which have been differentially affected by anthropogenic eutrophication. Our data set comprises 16,173 loci genotyped across 138 whitefish using restriction-site associated DNA sequencing (RADseq). Our analysis suggests that in each of the two lakes, the population of a different, but ecologically similar, whitefish species declined following a recent period of eutrophication. Genomic signatures consistent with hybridization are more pronounced in the more severely impacted lake. Comparisons between sympatric pairs of whitefish species with contrasting ecology, where one is shallow benthic and the other one more profundal pelagic, reveal genomic differentiation that is largely correlated along the genome, while differentiation is uncorrelated between pairs of allopatric provenance with similar ecology. We identify four genomic loci that provide evidence of parallel divergent adaptation between the shallow benthic species and the two different more profundal species. Functional annotations available for two of those loci are consistent with divergent ecological adaptation. Our genomic analysis indicates the action of divergent natural selection between sympatric whitefish species in pre-alpine lakes and reveals the vulnerability of these species to anthropogenic alterations of the environment and associated adaptive landscape.

Differential introgression of a female competitive trait in a hybrid zone between sex-role reversed species

Mating behavior between recently diverged species in secondary contact can impede or promote reproductive isolation. Traditionally, researchers focus on the importance of female mate choice and male-male competition in maintaining or eroding species barriers. Although female-female competition is widespread, little is known about its role in the speciation process. Here, we investigate a case of interspecific female competition and its influence on patterns of phenotypic and genetic introgression between species. We examine a hybrid zone between sex-role reversed, Neotropical shorebird species, the northern jacana (Jacana spinosa) and wattled jacana (J. jacana), in which female-female competition is a major determinant of reproductive success. Previous work found that females of the more aggressive and larger species, J. spinosa, disproportionately mother hybrid offspring, potentially by monopolizing breeding territories in sympatry with J. jacana. We find a cline shift of female body mass relative to the genetic center of the hybrid zone, consistent with asymmetric introgression of this competitive trait. We suggest that divergence in sexual characteristics between sex-role reversed females can influence patterns of gene flow upon secondary contact, similar to males in systems with more typical sex roles.

Prevalence of disruptive selection predicts extent of species differentiation in Lake Victoria cichlids

Theory suggests that speciation with gene flow is most likely when both sexual and ecological selection are divergent or disruptive. Divergent sexual and natural selection on the visual system have been demonstrated before in sympatric, morphologically similar sister species of Lake Victoria cichlids, but this does not explain the subtle morphological differences between them. To investigate the significance of natural selection on morphology during speciation, we here ask whether the prevalence of disruptive ecological selection differs between sympatric sister species that are at different stages of speciation. Some of our species pairs do (Pundamilia) and others do not (Neochromis) differ distinctively in sexually selected male nuptial coloration. We find that (i) evidence for disruptive selection, and for evolutionary response to it, is prevalent in traits that are differentiated between sister species; (ii) prevalence of both predicts the extent of genetic differentiation; and (iii) genetic differentiation is weaker in species pairs with conserved male nuptial coloration. Our results speak to the existence of two different mechanisms of speciation with gene flow: speciation mainly by sexual selection tightly followed by ecological character displacement in some cases and speciation mainly by divergent ecological selection in others.

The onset of ecological diversification 50 years after colonization of a crater lake by haplochromine cichlid fishes

Adaptive radiation research typically relies on the study of evolution in retrospective, leaving the predictive value of the concept hard to evaluate. Several radiations, including the cichlid fishes in the East African Great Lakes, have been studied extensively, yet no study has investigated the onset of the intraspecific processes of niche expansion and differentiation shortly after colonization of an adaptive zone by cichlids. Haplochromine cichlids of one of the two lineages that seeded the Lake Victoria radiation recently arrived in Lake Chala, a lake perfectly suited for within-lake cichlid speciation. Here, we infer the colonization and demographic history, quantify phenotypic, ecological and genomic diversity and diversification, and investigate the selection regime to ask if the population shows signs of diversification resembling the onset of adaptive radiation. We find that since their arrival in the lake, haplochromines have colonized a wide range of depth habitats associated with ecological and morphological expansion and the beginning of phenotypic differentiation and potentially nascent speciation, consistent with the very early onset of an adaptive radiation process. Moreover, we demonstrate evidence of rugged phenotypic fitness surfaces, indicating that current ecological selection may contribute to the phenotypic diversification.

A Dense Linkage Map of Lake Victoria Cichlids Improved the Pundamilia Genome Assembly and Revealed a Major QTL for Sex-Determination

Genetic linkage maps are essential for comparative genomics, high quality genome sequence assembly and fine scale quantitative trait locus (QTL) mapping. In the present study we identified and genotyped markers via restriction-site associated DNA (RAD) sequencing and constructed a genetic linkage map based on 1,597 SNP markers of an interspecific F2 cross of two closely related Lake Victoria cichlids (Pundamilia pundamilia and P sp. 'red head'). The SNP markers were distributed on 22 linkage groups and the total map size was 1,594 cM with an average marker distance of 1.01 cM. This high-resolution genetic linkage map was used to anchor the scaffolds of the Pundamilia genome and estimate recombination rates along the genome. Via QTL mapping we identified a major QTL for sex in a ∼1.9 Mb region on Pun-LG10, which is homologous to Oreochromis niloticus LG 23 (Ore-LG23) and includes a well-known vertebrate sex-determination gene (amh).