Genetic Coupling of Female Mate Choice with Polygenic Ecological Divergence Facilitates Stickleback Speciation

RIDGE, a tool tailored to detect gene flow barriers across species pairs

Characterizing the processes underlying reproductive isolation between diverging lineages is central to understanding speciation. Here, we present RIDGE-Reproductive Isolation Detection using Genomic polymorphisms-a tool tailored for quantifying gene flow barrier proportion and identifying the relevant genomic regions. RIDGE relies on an Approximate Bayesian Computation with a model-averaging approach to accommodate diverse scenarios of lineage divergence. It captures heterogeneity in effective migration rate along the genome while accounting for variation in linked selection and recombination. The barrier detection test relies on numerous summary statistics to compute a Bayes factor, offering a robust statistical framework that facilitates cross-species comparisons. Simulations revealed RIDGE's efficiency in capturing signals of ongoing migration. Model averaging proved particularly valuable in scenarios of high model uncertainty where no migration or migration homogeneity can be wrongly assumed, typically for recent divergence times <0.1 2Ne generations. Applying RIDGE to four published crow data sets, we first validated our tool by identifying a well-known large genomic region associated with mate choice patterns. Second, while we identified a significant overlap of outlier loci using RIDGE and traditional genomic scans, our results suggest that a substantial portion of previously identified outliers are likely false positives. Outlier detection relies on allele differentiation, relative measures of divergence and the count of shared polymorphisms and fixed differences. Our analyses also highlight the value of incorporating multiple summary statistics including our newly developed outlier ones that can be useful in challenging detection conditions.

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.

Does paternal care influence mate preference? Male and female mating behavior in Threespine Stickleback ecotypes that differ markedly in parental care

Reproductive isolation can occur due to divergence in sexual selection for particular traits. For example, differences in mate preference associated with body size can play an important role in divergence between groups. The importance of mate preference for population divergence may be influenced by other aspects of a mating system, such as the requirement for parental care. In Nova Scotia, Canada, two ecotypes of marine Threespine Stickleback occur sympatrically: a “common” ecotype wherein males provide parental care, and a “white” ecotype that does not exhibit paternal care. The goal of our study was to examine differences in male mate preference between white and common stickleback males to test the prediction that males who invest more in parental care may be more selective about their mates. Because of the link between size and fecundity in this species, we predict that males that invest in parental care should prefer large females, while males that do not provide care will not exhibit preference for larger female size. We found that common male stickleback preferred larger‐bodied females of both ecotypes, while white males showed a preference for larger‐bodied common females. Secondarily, we assessed whether females differed in their willingness to mate with males of different sizes and ecotypes. Common female stickleback had a higher response rate toward smaller white males, which may be associated with their relatively high courtship rates. Counter to previous studies on these ecotypes that suggest that mating is completely assortative, interecotype matings occurred in half of the observed spawning events. This observation, coupled with the results that males may prefer females based mainly on size and females respond to males who court more rigorously regardless of their ecotype, may lend insight into recent genetic evidence for hybridization in the wild.

Life on a beach leads to phenotypic divergence despite gene flow for an island lizard

Limited spatial separation within small islands suggests that observed population divergence may occur due to habitat differences without interruption to gene flow but strong evidence of this is scarce. The wall lizard Teira dugesii lives in starkly contrasting shingle beach and inland habitats on the island of Madeira. We used a matched pairs sampling design to examine morphological and genomic divergence between four beach and adjacent (<1 km) inland areas. Beach populations are significantly darker than corresponding inland populations. Geometric morphometric analyses reveal divergence in head morphology: beach lizards have generally wider snouts. Genotyping-by-sequencing allows the rejection of the hypothesis that beach populations form a distinct lineage. Bayesian analyses provide strong support for models that incorporate gene flow, relative to those that do not, replicated at all pairs of matched sites. Madeiran lizards show morphological divergence between habitats in the face of gene flow, revealing how divergence may originate within small islands.

Evolution of physical linkage between loci controlling ecological traits and mating preferences

Coupling of multiple barriers to gene-flow, such as divergent local adaptation and reproductive isolation, facilitates speciation. However, alleles at loci that contribute to barrier effects can be dissociated by recombination. Models of linkage between diverging alleles often consider elements that reduce recombination, such as chromosomal inversions and alleles that modify recombination rate between existing loci. In contrast, here, we consider the evolution of linkage due to the close proximity of loci on the same chromosome. Examples of such physical linkage exist in several species, but in other cases, strong associations are maintained without physical linkage. We use an individual-based model to study the conditions under which the physical linkage between loci controlling ecological traits and mating preferences might be expected to evolve. We modelled a single locus controlling an ecological trait that acts also as a mating cue. Mating preferences are controlled by multiple loci, formed by mutations that are randomly placed in the "genome", within varying distances from the ecological trait locus, allowing us to examine which genomic architectures spread across the population. Our model reveals that stronger physical linkage is favoured when mating preferences and selection are weaker. Under such conditions mating among divergent phenotypes is more frequent, and matching ecological trait and mating preference alleles are more likely to become dissociated by recombination, favouring the evolution of genetic linkage. While most theoretical studies on clustering of divergent loci focus on how physical linkage influences speciation, we show how physical linkage itself can arise, establishing conditions that can favour speciation.

Genome-wide species delimitation analyses of a silverside fish species complex in central Mexico indicate taxonomic over-splitting

BACKGROUND

Delimiting species across a speciation continuum is a complex task, as the process of species origin is not generally instantaneous. The use of genome-wide data provides unprecedented resolution to address convoluted species delimitation cases, often unraveling cryptic diversity. However, because genome-wide approaches based on the multispecies coalescent model are known to confound population structure with species boundaries, often resulting in taxonomic over-splitting, it has become increasingly evident that species delimitation research must consider multiple lines of evidence. In this study, we used phylogenomic, population genomic, and coalescent-based species delimitation approaches, and examined those in light of morphological and ecological information, to investigate species numbers and boundaries comprising the Chirostoma "humboltianum group" (family Atherinidae). The humboltianum group is a taxonomically controversial species complex where previous morphological and mitochondrial studies produced conflicting species delimitation outcomes. We generated ddRADseq data for 77 individuals representing the nine nominal species in the group, spanning their distribution range in the central Mexican plateau.

RESULTS

Our results conflict with the morphospecies and ecological delimitation hypotheses, identifying four independently evolving lineages organized in three geographically cohesive clades: (i) chapalae and sphyraena groups in Lake Chapala, (ii) estor group in Lakes Pátzcuaro and Zirahuén, and (iii) humboltianum sensu stricto group in Lake Zacapu and Lerma river system.

CONCLUSIONS

Overall, our study provides an atypical example where genome-wide analyses delineate fewer species than previously recognized on the basis of morphology. It also highlights the influence of the geological history of the Chapala-Lerma hydrological system in driving allopatric speciation in the humboltianum group.

Three problems in the genetics of speciation by selection

Speciation is the process by which barriers to gene flow evolve between populations. Although we now know that speciation is largely driven by natural selection, knowledge of the agents of selection and the genetic and genomic mechanisms that facilitate divergence is required for a satisfactory theory of speciation. In this essay, we highlight three advances/problems in our understanding of speciation that have arisen from studies of the genes and genomic regions that underlie the evolution of reproductive isolation. First, we describe how the identification of “speciation” genes makes it possible to identify the agents of selection causing the evolution of reproductive isolation, while also noting that the link between the genetics of phenotypic divergence and intrinsic postzygotic reproductive barriers remains tenuous. Second, we discuss the important role of recombination suppressors in facilitating speciation with gene flow, but point out that the means and timing by which reproductive barriers become associated with recombination cold spots remains uncertain. Third, we establish the importance of ancient genetic variation in speciation, although we argue that the focus of speciation studies on evolutionarily young groups may bias conclusions in favor of ancient variation relative to new mutations.

What is adaptive radiation? Many manifestations of the phenomenon in an iconic lineage of Hawaiian spiders

Adaptive radiation provides the ideal context for identifying and testing the processes that drive evolutionary diversification. However, different adaptive radiations show a variety of different patterns, making it difficult to come up with universal rules that characterize all such systems. Diversification may occur via several mechanisms including non-adaptive divergence, adaptation to novel environments, or character displacement driven by competition. Here, we characterize the ways these different drivers contribute to present-day diversity patterns, using the exemplary adaptive radiation of Hawaiian long-jawed orbweaver (Tetragnatha) spiders. We present the most taxonomically comprehensive phylogenetic hypothesis to date for this group, using 10 molecular markers and representatives from every known species across the archipelago. Among the lineages that make up this remarkable radiation, we find evidence for multiple diversification modalities. Several clades appear to have diversified in allopatry under a narrow range of ecological conditions, highlighting the role of niche conservatism in speciation. Others have shifted into new environments and evolved traits that appear to be adaptive in those environments. Still others show evidence for character displacement by close relatives, often resulting in convergent evolution of stereotyped ecomorphs. All of the above mechanisms seem to have played a role in giving rise to the exceptional diversity of morphological, ecological and behavioral traits represented among the many species of Hawaiian Tetragnatha. Taking all these processes into account, and testing how they operate in different systems, may allow us to identify universal principles underlying adaptive radiation.

Non-parallel morphological divergence following colonization of a new host plant

Adaptation to new ecological niches is known to spur population diversification and may lead to speciation if gene flow is ceased. While adaptation to the same ecological niche is expected to be parallel, it is more difficult to predict whether selection against maladaptive hybridization in secondary sympatry results in parallel divergence also in traits that are not directly related to the ecological niches. Such parallelisms in response to selection for reproductive isolation can be identified through estimating parallelism in reproductive character displacement across different zones of secondary contact. Here, we use a host shift in the phytophagous peacock fly Tephritis conura, with both host races represented in two geographically separate areas East and West of the Baltic Sea to investigate convergence in morphological adaptations. We asked (i) if there are consistent morphological adaptations to a host plant shift and (ii) if the response to secondary sympatry with the alternate host race is parallel across contact zones. We found surprisingly low and variable, albeit significant, divergence between host races. Only one trait, the length of the female ovipositor, which serves an important function in the interaction with the hosts, was consistently different between host races. Instead, co-existence with the other host race significantly affected the degree of morphological divergence, but the divergence was largely driven by different traits in different contact zones. Thus, local stochastic fixation or reinforcement could generate trait divergence, and additional evidence is needed to conclude whether divergence is locally adaptive.

Heterosis counteracts hybrid breakdown to forestall speciation by parallel natural selection

In contrast to ecological speciation, where reproductive isolation evolves as a consequence of divergent natural selection, speciation by parallel natural selection has been less thoroughly studied. To test whether parallel evolution drives speciation, we leveraged the repeated evolution of benthic and limnetic ecotypes of threespine stickleback fish and estimated fitness for pure crosses and within-ecotype hybrids in semi-natural ponds and in laboratory aquaria. In ponds, we detected hybrid breakdown in both ecotypes but this was counterbalanced by heterosis and the strength of post-zygotic isolation was nil. In aquaria, we detected heterosis in limnetic crosses and breakdown in benthic crosses, which is suggestive of process- and ecotype-specific environment-dependence. In ponds, heterosis and breakdown were three times greater in limnetic crosses than in benthic crosses, contrasting the prediction that the fitness consequences of hybridization should be greater in crosses among more derived ecotypes. Consistent with a primary role for stochastic processes, patterns differed among crosses between populations from different lakes. Yet, the observation of qualitatively similar patterns of heterosis and hybrid breakdown for both ecotypes when averaging the lake pairs indicates that the outcome of hybridization is repeatable in a general sense.

Imbalanced segregation of recombinant haplotypes in hybrid populations reveals inter- and intrachromosomal Dobzhansky-Muller incompatibilities

Dobzhansky-Muller incompatibilities (DMIs) are a major component of reproductive isolation between species. DMIs imply negative epistasis and are exposed when two diverged populations hybridize. Mapping the locations of DMIs has largely relied on classical genetic mapping. Approaches to date are hampered by low power and the challenge of identifying DMI loci on the same chromosome, because strong initial linkage of parental haplotypes weakens statistical tests. Here, we propose new statistics to infer negative epistasis from haplotype frequencies in hybrid populations. When two divergent populations hybridize, the variance in heterozygosity at two loci decreases faster with time at DMI loci than at random pairs of loci. When two populations hybridize at near-even admixture proportions, the deviation of the observed variance from its expectation becomes negative for the DMI pair. This negative deviation enables us to detect intermediate to strong negative epistasis both within and between chromosomes. In practice, the detection window in hybrid populations depends on the demographic scenario, the recombination rate, and the strength of epistasis. When the initial proportion of the two parental populations is uneven, only strong DMIs can be detected with our method unless migration prevents parental haplotypes from being lost. We use the new statistics to infer candidate DMIs from three hybrid populations of swordtail fish. We identify numerous new DMI candidates, some of which are inferred to interact with several loci within and between chromosomes. Moreover, we discuss our results in the context of an expected enrichment in intrachromosomal over interchromosomal DMIs.

Analysis of ancestry heterozygosity suggests that hybrid incompatibilities in threespine stickleback are environment dependent

Hybrid incompatibilities occur when interactions between opposite ancestry alleles at different loci reduce the fitness of hybrids. Most work on incompatibilities has focused on those that are "intrinsic," meaning they affect viability and sterility in the laboratory. Theory predicts that ecological selection can also underlie hybrid incompatibilities, but tests of this hypothesis using sequence data are scarce. In this article, we compiled genetic data for F2 hybrid crosses between divergent populations of threespine stickleback fish (Gasterosteus aculeatus L.) that were born and raised in either the field (seminatural experimental ponds) or the laboratory (aquaria). Because selection against incompatibilities results in elevated ancestry heterozygosity, we tested the prediction that ancestry heterozygosity will be higher in pond-raised fish compared to those raised in aquaria. We found that ancestry heterozygosity was elevated by approximately 3% in crosses raised in ponds compared to those raised in aquaria. Additional analyses support a phenotypic basis for incompatibility and suggest that environment-specific single-locus heterozygote advantage is not the cause of selection on ancestry heterozygosity. Our study provides evidence that, in stickleback, a coarse-albeit indirect-signal of environment-dependent hybrid incompatibility is reliably detectable and suggests that extrinsic incompatibilities can evolve before intrinsic incompatibilities.

Flexibility of nuptial colouration in a unique ecotype of threespine stickleback (Gasterosteus aculeatus)

Nuptial colouration in animals may serve as a signal of competitor and (or) mate quality during breeding. In many temperate fishes, nuptial colouration develops during discrete breeding seasons and is a target of sexual selection. We examine nuptial colouration and behaviour of a unique ecotype of threespine stickleback (Gasterosteus aculeatus Linnaeus, 1758), wherein males turn from dull brown-grey to pearlescent white during the breeding season. The main goal of this work was to determine the relative role of white colouration in intersexual competition and mate choice. In a combination of field and laboratory work, we found that males are brightest white when engaging in courtship activities in the presence of a female; this indicates that white colouration may be primarily related to enhancing signalling during mate attraction. White colouration intensity increased as the breeding season progressed and may be related to an influx of conspecifics. Colour change from cryptic grey to bright white occurred rapidly (<90 s) and may be deployed to enhance behavioural signals. We conclude that bright white colouration in the white ecotype is a potential signal of mate quality and may have evolved from a previously existing capacity for colour plasticity in common threespine stickleback.

Light environment influences mating behaviours during the early stages of divergence in tropical butterflies

Speciation is facilitated when traits under divergent selection also act as mating cues. Fluctuations in sensory conditions can alter signal perception independently of adaptation to the broader sensory environment, but how this fine-scale variation may constrain or promote behavioural isolation has received little attention. The warning patterns of Heliconius butterflies are under selection for aposematism and act as mating cues. Using computer vision, we extracted behavioural data from 1481 h of video footage, for 387 individuals. We show that the putative hybrid species H. heurippa and its close relative H. timareta linaresi differ in their response to divergent warning patterns, but that these differences are strengthened with increased local illuminance. Trials with live individuals reveal low-level assortative mating that is sufficiently explained by differences in visual attraction. Finally, results from hybrid butterflies are consistent with linkage between a major warning pattern gene and the corresponding behaviour, though the differences in behaviour we observe are unlikely to cause rapid reproductive isolation as predicted under a model of hybrid trait speciation. Overall, our results reveal that the contribution of ecological mating cues to reproductive isolation may depend on the immediate sensory conditions during which they are displayed to conspecifics.

Extensive Linkage and Genetic Coupling of Song and Preference Loci Underlying Rapid Speciation in Laupala Crickets

In nature, closely related species commonly display divergent mating behaviors, suggesting a central role for such traits in the origin of species. Elucidating the genetic basis of divergence in these traits is necessary to understand the evolutionary process leading to reproductive barriers and speciation. The rapidly speciating Hawaiian crickets of the genus Laupala provides an ideal system for dissecting the genetic basis of mating behavior divergence. In Laupala, closely related species differ markedly in male song pulse rate and female preference for pulse rate. These behaviors play an important role in determining mating patterns. Previous studies identified a genetic architecture consisting of numerous small to moderate effect loci causing interspecific differences in pulse rate and preference, including colocalizing pulse rate and preference QTL on linkage group one (LG1). To further interrogate these QTL, we conduct a fine mapping study using high-density SNP linkage maps. With improved statistical power and map resolution, we provide robust evidence for genetic coupling between song and preference, along with two additional pulse rate QTL on LG1, revealing a more resolved picture of the genetic architecture underlying mating behavior divergence. Our sequence-based genetic map, along with dramatically narrowed QTL confidence intervals, allowed us to annotate genes within the QTL regions and identify several exciting candidate genes underlying variation in pulse rate and preference divergence. Such knowledge suggests potential molecular mechanisms underlying the evolution of behavioral barriers.

Courtship behavior, nesting microhabitat, and assortative mating in sympatric stickleback species pairs

The maintenance of reproductive isolation in the face of gene flow is a particularly contentious topic, but differences in reproductive behavior may provide the key to explaining this phenomenon. However, we do not yet fully understand how behavior contributes to maintaining species boundaries. How important are behavioral differences during reproduction? To what extent does assortative mating maintain reproductive isolation in recently diverged populations and how important are "magic traits"? Assortative mating can arise as a by-product of accumulated differences between divergent populations as well as an adaptive response to contact between those populations, but this is often overlooked. Here we address these questions using recently described species pairs of three-spined stickleback (Gasterosteus aculeatus), from two separate locations and a phenotypically intermediate allopatric population on the island of North Uist, Scottish Western Isles. We identified stark differences in the preferred nesting substrate and courtship behavior of species pair males. We showed that all males selectively court females of their own ecotype and all females prefer males of the same ecotype, regardless of whether they are from species pairs or allopatric populations. We also showed that mate choice does not appear to be driven by body size differences (a potential "magic trait"). By explicitly comparing the strength of these mating preferences between species pairs and single-ecotype locations, we were able to show that present levels of assortative mating due to direct mate choice are likely a by-product of other adaptations between ecotypes, and not subject to obvious selection in species pairs. Our results suggest that ecological divergence in mating characteristics, particularly nesting microhabitat may be more important than direct mate choice in maintaining reproductive isolation in stickleback species pairs.

Clustering of loci controlling species differences in male chemical bouquets of sympatric Heliconius butterflies

The degree to which loci promoting reproductive isolation cluster in the genome-that is, the genetic architecture of reproductive isolation-can influence the tempo and mode of speciation. Tight linkage between these loci can facilitate speciation in the face of gene flow. Pheromones play a role in reproductive isolation in many Lepidoptera species, and the role of endogenously produced compounds as secondary metabolites decreases the likelihood of pleiotropy associated with many barrier loci. Heliconius butterflies use male sex pheromones to both court females (aphrodisiac wing pheromones) and ward off male courtship (male-transferred antiaphrodisiac genital pheromones), and it is likely that these compounds play a role in reproductive isolation between Heliconius species. Using a set of backcross hybrids between H. melpomene and H. cydno, we investigated the genetic architecture of putative male pheromone compound production. We found a set of 40 significant quantitative trait loci (QTL) representing 33 potential pheromone compounds. QTL clustered significantly on two chromosomes, chromosome 8 for genital compounds and chromosome 20 for wing compounds, and chromosome 20 was enriched for potential pheromone biosynthesis genes. There was minimal overlap between pheromone QTL and known QTL for mate choice and color pattern. Nonetheless, we did detect linkage between a QTL for wing androconial area and optix, a color pattern locus known to play a role in reproductive isolation in these species. This tight clustering of putative pheromone loci might contribute to coincident reproductive isolating barriers, facilitating speciation despite ongoing gene flow.

The effects of the sex chromosomes on the inheritance of species-specific traits of the copulatory organ shape in Drosophila virilis and Drosophila lummei

The shape of the male genitalia in many taxa is the most rapidly evolving morphological structure, often driving reproductive isolation, and is therefore widely used in systematics as a key character to distinguish between sibling species. However, only a few studies have used the genital arch of the male copulatory organ as a model to study the genetic basis of species-specific differences in the Drosophila copulatory system. Moreover, almost nothing is known about the effects of the sex chromosomes on the shape of the male mating organ. In our study, we used a set of crosses between D. virilis and D. lummei and applied the methods of quantitative genetics to assess the variability of the shape of the male copulatory organ and the effects of the sex chromosomes and autosomes on its variance. Our results showed that the male genital shape depends on the species composition of the sex chromosomes and autosomes. Epistatic interactions of the sex chromosomes with autosomes and the species origin of the Y-chromosome in a male in interspecific crosses also influenced the expression of species-specific traits in the shape of the male copulatory system. Overall, the effects of sex chromosomes were comparable to the effects of autosomes despite the great differences in gene numbers between them. It may be reasonably considered that sexual selection for specific genes associated with the shape of the male mating organ prevents the demasculinization of the X chromosome.

Tadpoles of hybridising fire-bellied toads (B. bombina and B. variegata) differ in their susceptibility to predation

The role of adaptive divergence in the formation of new species has been the subject of much recent debate. The most direct evidence comes from traits that can be shown to have diverged under natural selection and that now contribute to reproductive isolation. Here, we investigate differential adaptation of two fire-bellied toads (Anura, Bombinatoridae) to two types of aquatic habitat. Bombina bombina and B. variegata are two anciently diverged taxa that now reproduce in predator-rich ponds and ephemeral aquatic sites, respectively. Nevertheless, they hybridise extensively wherever their distribution ranges adjoin. We show in laboratory experiments that, as expected, B. variegata tadpoles are at relatively greater risk of predation from dragonfly larvae, even when they display a predator-induced phenotype. These tadpoles spent relatively more time swimming and so prompted more attacks from the visually hunting predators. We argue in the discussion that genomic regions linked to high activity in B. variegata should be barred from introgression into the B. bombina gene pool and thus contribute to gene flow barriers that keep the two taxa from merging into one.

Putative chromosomal rearrangements are associated primarily with ecotype divergence rather than geographic separation in an intertidal, poorly dispersing snail

Littorina saxatilis is becoming a model system for understanding the genomic basis of ecological speciation. The parallel formation of crab-adapted ecotypes that exhibit partial reproductive isolation from wave-adapted ecotypes has enabled genomic investigation of conspicuous shell traits. Recent genomic studies suggest that chromosomal rearrangements may enable ecotype divergence by reducing gene flow. However, the genomic architecture of traits that are divergent between ecotypes remains poorly understood. Here, we use 11,504 single nucleotide polymorphism (SNP) markers called using the recently released L. saxatilis genome to genotype 462 crab ecotype, wave ecotype and phenotypically intermediate Spanish L. saxatilis individuals with scored phenotypes. We used redundancy analysis to study the genetic architecture of loci associated with shell shape, shape corrected for size, shell size and shell ornamentation, and to compare levels of co-association among different traits. We discovered 341 SNPs associated with shell traits. Loci associated with trait divergence between ecotypes were often located inside putative chromosomal rearrangements recently characterized in Swedish L. saxatilis. In contrast, we found that shell shape corrected for size varied primarily by geographic site rather than by ecotype and showed little association with these putative rearrangements. We conclude that genomic regions of elevated divergence inside putative rearrangements were associated with divergence of L. saxatilis ecotypes along steep environmental axes-consistent with models of adaptation with gene flow-but were not associated with divergence among the three geographical sites. Our findings support predictions from models indicating the importance of genomic regions of reduced recombination allowing co-association of loci during ecological speciation with ongoing gene flow.

Genomic landscape of reproductive isolation in Lucania killifish: The role of sex loci and salinity

Adaptation to different environments can directly and indirectly generate reproductive isolation between species. Bluefin killifish (Lucania goodei) and rainwater killifish (L. parva) are sister species that have diverged across a salinity gradient and are reproductively isolated by habitat, behavioural, extrinsic and intrinsic post-zygotic isolation. We asked if salinity adaptation contributes indirectly to other forms of reproductive isolation via linked selection and hypothesized that low recombination regions, such as sex chromosomes or chromosomal rearrangements, might facilitate this process. We conducted QTL mapping in backcrosses between L. parva and L. goodei to explore the genetic architecture of salinity tolerance, behavioural isolation and intrinsic isolation. We mapped traits relative to a chromosome that has undergone a centric fusion in L. parva (relative to L. goodei). We found that the sex locus appears to be male determining (XX-XY), was located on the fused chromosome and was implicated in intrinsic isolation. QTL associated with salinity tolerance were spread across the genome and did not overly co-localize with regions associated with behavioural or intrinsic isolation. This preliminary analysis of the genetic architecture of reproductive isolation between Lucania species does not support the hypothesis that divergent natural selection for salinity tolerance led to behavioural and intrinsic isolation as a by-product. Combined with previous studies in this system, our work suggests that adaptation as a function of salinity contributes to habitat isolation and that reinforcement may have contributed to the evolution of behavioural isolation instead, possibly facilitated by linkage between behavioural isolation and intrinsic isolation loci on the fused chromosome.

Visual mate preference evolution during butterfly speciation is linked to neural processing genes

Many animal species remain separate not because their individuals fail to produce viable hybrids but because they "choose" not to mate. However, we still know very little of the genetic mechanisms underlying changes in these mate preference behaviours. Heliconius butterflies display bright warning patterns, which they also use to recognize conspecifics. Here, we couple QTL for divergence in visual preference behaviours with population genomic and gene expression analyses of neural tissue (central brain, optic lobes and ommatidia) across development in two sympatric Heliconius species. Within a region containing 200 genes, we identify five genes that are strongly associated with divergent visual preferences. Three of these have previously been implicated in key components of neural signalling (specifically an ionotropic glutamate receptor and two regucalcins), and overall our candidates suggest shifts in behaviour involve changes in visual integration or processing. This would allow preference evolution without altering perception of the wider environment.

Genes and speciation: is it time to abandon the biological species concept?

The biological species concept (BSC) is the cornerstone of neo-Darwinian thinking. In BSC, species do not exchange genes either during or after speciation. However, as gene flow during speciation is increasingly being reported in a substantial literature, it seems time to reassess the revered, but often doubted, BSC. Contrary to the common perception, BSC should expect substantial gene flow at the onset of speciation, not least because geographical isolation develops gradually. Although BSC does not stipulate how speciation begins, it does require a sustained period of isolation for speciation to complete its course. Evidence against BSC must demonstrate that the observed gene flow does not merely occur at the onset of speciation but continues until its completion. Importantly, recent genomic analyses cannot reject this more realistic version of BSC, although future analyses may still prove it wrong. The ultimate acceptance or rejection of BSC is not merely about a historical debate; rather, it is about the fundamental nature of species - are species (and, hence, divergent adaptations) driven by a relatively small number of genes, or by thousands of them? Many levels of biology, ranging from taxonomy to biodiversity, depend on this resolution.

On the logic of Fisherian sexual selection

In Fisher's model of sexual selection, a female preference for a male trait spreads together with the trait because their genetic bases become correlated. This can be interpreted as a "greenbeard" system: a preference gene, by inducing a female to mate with a trait-bearing male, favors itself because the male is disproportionately likely also to carry the preference gene. Here, we use this logic to argue that Fisherian sexual selection in diploids proceeds via two channels: (i) trait-bearing males are disproportionately the product of matings between preference-bearing mothers and trait-bearing fathers, and thus trait and preference genes are correlated "in trans"; (ii) trait and preference genes come into gametic phase disequilibrium, and thus are correlated "in cis." Gametic phase disequilibrium is generated by three distinct mechanisms that we identify. The trans channel does not operate when sexual selection is restricted to the haploid phase, and therefore represents a fundamental difference between haploid and diploid models of sexual selection. We show that the cis and trans channels contribute equally to the spread of the preference when recombination between the preference and trait loci is free, but that the trans channel is substantially more important when linkage is tight.

Chemical signals act as the main reproductive barrier between sister and mimetic Heliconius butterflies

Colour pattern is the main trait that drives mate recognition between Heliconius species that are phylogenetically close. However, when this cue is compromised such as in cases of mimetic, sympatric and closely related species, alternative mating signals must evolve to ensure reproductive isolation and species integrity. The closely related species Heliconius melpomene malleti and H. timareta florencia occur in the same geographical region, and despite being co-mimics, they display strong reproductive isolation. In order to test which cues differ between species, and potentially contribute to reproductive isolation, we quantified differences in the wing phenotype and the male chemical profile. As expected, the wing colour pattern was indistinguishable between the two species, while the chemical profile of the androconial and genital males' extracts showed marked differences. We then conducted behavioural experiments to study the importance of these signals in mate recognition by females. In agreement with our previous results, we found that chemical blends and not wing colour pattern drive the preference of females for conspecific males. Also, experiments with hybrid males and females suggested an important genetic component for both chemical production and preference. Altogether, these results suggest that chemicals are the major reproductive barrier opposing gene flow between these two sister and co-mimic species.

Within-generation and transgenerational plasticity of mate choice in oceanic stickleback under climate change

Plasticity, both within and across generations, can shape sexual traits involved in mate choice and reproductive success, and thus direct measures of fitness. Especially, transgenerational plasticity (TGP), where parental environment influences offspring plasticity in future environments, could compensate for otherwise negative effects of environmental change on offspring sexual traits. We conducted a mate choice experiment using stickleback ( Gasterosteus aculeatus) with different thermal histories (ambient 17°C or elevated 21°C) within and across generations under simulated ocean warming using outdoor mesocosms. Parentage analysis of egg clutches revealed that maternal developmental temperature and reproductive (mesocosm) environment affected egg size, with females that developed at 17°C laying smaller eggs in 21°C mesocosms, likely owing to metabolic costs at elevated temperature. Paternal developmental temperature interacted with the reproductive environment to influence mating success, particularly under simulated ocean warming, with males that developed at 21°C showing lower overall mating success compared with 17°C males, but higher mating success in 21°C mesocosms. Furthermore, mating success of males was influenced by the interaction between F1 developmental temperature and F0 parent acclimation temperature, demonstrating the potential role of both TGP and within-generation plasticity in shaping traits involved in sexual selection and mate choice, potentially facilitating rapid responses to environmental change. This article is part of the theme issue 'The role of plasticity in phenotypic adaptation to rapid environmental change'.

Genetic coupling of signal and preference facilitates sexual isolation during rapid speciation

The divergence of sexual signals is ultimately a coevolutionary process: while signals and preferences diverge between lineages, they must remain coordinated within lineages for matings to occur. Divergence in sexual signals makes a major contribution to evolving species barriers. Therefore, the genetic architecture underlying signal-preference coevolution is essential to understanding speciation but remains largely unknown. In Laupala crickets where male song pulse rate and female pulse rate preferences have coevolved repeatedly and rapidly, we tested two contrasting hypotheses for the genetic architecture underlying signal-preference coevolution: linkage disequilibrium between unlinked loci and genetic coupling (linkage disequilibrium resulting from pleiotropy of a shared locus or tight physical linkage). Through selective introgression and quantitative trait locus (QTL) fine mapping, we estimated the location of QTL underlying interspecific variation in both female preference and male pulse rate from the same mapping populations. Remarkably, map estimates of the pulse rate and preference loci are as close as 0.06 cM apart, the strongest evidence to date for genetic coupling between signal and preference loci. As the second pair of colocalizing signal and preference loci in the Laupala genome, our finding supports an intriguing pattern, pointing to a major role for genetic coupling in the quantitative evolution of a reproductive barrier and rapid speciation in Laupala. Owing to its effect on suppressing recombination, a coupled, quantitative genetic architecture offers a powerful and parsimonious genetic mechanism for signal-preference coevolution and the establishment of positive genetic covariance on which the Fisherian runaway process of sexual selection relies.

Major Histocompatibility Complex (MHC) Genes and Disease Resistance in Fish

Fascinating about classical major histocompatibility complex (MHC) molecules is their polymorphism. The present study is a review and discussion of the fish MHC situation. The basic pattern of MHC variation in fish is similar to mammals, with MHC class I versus class II, and polymorphic classical versus nonpolymorphic nonclassical. However, in many or all teleost fishes, important differences with mammalian or human MHC were observed: (1) The allelic/haplotype diversification levels of classical MHC class I tend to be much higher than in mammals and involve structural positions within but also outside the peptide binding groove; (2) Teleost fish classical MHC class I and class II loci are not linked. The present article summarizes previous studies that performed quantitative trait loci (QTL) analysis for mapping differences in teleost fish disease resistance, and discusses them from MHC point of view. Overall, those QTL studies suggest the possible importance of genomic regions including classical MHC class II and nonclassical MHC class I genes, whereas similar observations were not made for the genomic regions with the highly diversified classical MHC class I alleles. It must be concluded that despite decades of knowing MHC polymorphism in jawed vertebrate species including fish, firm conclusions (as opposed to appealing hypotheses) on the reasons for MHC polymorphism cannot be made, and that the types of polymorphism observed in fish may not be explained by disease-resistance models alone.

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.

Genetic dissection of assortative mating behavior

The evolution of new species is made easier when traits under divergent ecological selection are also mating cues. Such ecological mating cues are now considered more common than previously thought, but we still know little about the genetic changes underlying their evolution or more generally about the genetic basis for assortative mating behaviors. Both tight physical linkage and the existence of large-effect preference loci will strengthen genetic associations between behavioral and ecological barriers, promoting the evolution of assortative mating. The warning patterns of Heliconius melpomene and H. cydno are under disruptive selection due to increased predation of nonmimetic hybrids and are used during mate recognition. We carried out a genome-wide quantitative trait locus (QTL) analysis of preference behaviors between these species and showed that divergent male preference has a simple genetic basis. We identify three QTLs that together explain a large proportion (approximately 60%) of the difference in preference behavior observed between the parental species. One of these QTLs is just 1.2 (0-4.8) centiMorgans (cM) from the major color pattern gene optix, and, individually, all three have a large effect on the preference phenotype. Genomic divergence between H. cydno and H. melpomene is high but broadly heterogenous, and admixture is reduced at the preference-optix color pattern locus but not the other preference QTLs. The simple genetic architecture we reveal will facilitate the evolution and maintenance of new species despite ongoing gene flow by coupling behavioral and ecological aspects of reproductive isolation.

The Genetics of Mating Song Evolution Underlying Rapid Speciation: Linking Quantitative Variation to Candidate Genes for Behavioral Isolation

Differences in mating behaviors evolve early during speciation, eventually contributing to reproductive barriers between species. Knowledge of the genetic and genomic basis of these behaviors is therefore integral to a causal understanding of speciation. Acoustic behaviors are often part of the mating ritual in animal species. The temporal rhythms of mating songs are notably species-specific in many vertebrates and arthropods and often underlie assortative mating. Despite discoveries of mutations that disrupt the temporal rhythm of these songs, we know surprisingly little about genes affecting naturally occurring variation in the temporal pattern of singing behavior. In the rapidly speciating Hawaiian cricket genus Laupala, the striking species variation in song rhythms constitutes a behavioral barrier to reproduction between species. Here, we mapped the largest-effect locus underlying interspecific variation in song rhythm between two Laupala species to a narrow genomic region, wherein we find no known candidate genes affecting song temporal rhythm in Drosophila Whole-genome sequencing, gene prediction, and functional annotation of this region reveal an exciting and promising candidate gene, the putative cyclic nucleotide-gated ion channel-like gene, for natural variation in mating behavior. Identification and molecular characterization of the candidate gene reveals a nonsynonymous mutation in a conserved binding domain, suggesting that ion channels are important targets of selection on rhythmic signaling during establishment of behavioral isolation and rapid speciation.

Role of sexual imprinting in assortative mating and premating isolation in Darwin's finches

Global biodiversity is being degraded at an unprecedented rate, so it is important to preserve the potential for future speciation. Providing for the future requires understanding speciation as a contemporary ecological process. Phylogenetically young adaptive radiations are a good choice for detailed study because diversification is ongoing. A key question is how incipient species become reproductively isolated from each other. Barriers to gene exchange have been investigated experimentally in the laboratory and in the field, but little information exists from the quantitative study of mating patterns in nature. Although the degree to which genetic variation underlying mate-preference learning is unknown, we provide evidence that two species of Darwin's finches imprint on morphological cues of their parents and mate assortatively. Statistical evidence of presumed imprinting is stronger for sons than for daughters and is stronger for imprinting on fathers than on mothers. In combination, morphology and species-specific song learned from the father constitute a barrier to interbreeding. The barrier becomes stronger the more the species diverge morphologically and ecologically. It occasionally breaks down, and the species hybridize. Hybridization is most likely to happen when species are similar to each other in adaptive morphological traits, e.g., body size and beak size and shape. Hybridization can lead to the formation of a new species reproductively isolated from the parental species as a result of sexual imprinting. Conservation of sufficiently diverse natural habitat is needed to sustain a large sample of extant biota and preserve the potential for future speciation.

Notes from a snail island: Littorinid evolution and adaptation

The most successful study systems are built on a foundation of decades of research on the basic biology, ecology and life history of the organisms in question. Combined with new technologies, this can provide a formidable means to address important issues in evolutionary biology and molecular ecology. Littorinid marine snails are a good example of this, with a rich literature on their taxonomy, speciation, thermal tolerance and behavioural adaptations. In August 2017, an international meeting on Littorinid evolution was held at the Tjärnö Marine Research Laboratory in Western Sweden. In this meeting review, I provide a summary of some of the exciting work on parallel evolution, sexual selection and adaptation to environmental stress presented there. I argue that newly available genomic resources present an opportunity for integrating the traditionally divergent fields of speciation and environmental adaptation in Littorinid research.

Female Choice Undermines the Emergence of Strong Sexual Isolation between Locally Adapted Populations of Atlantic Mollies (Poecilia mexicana)

Divergent selection between ecologically dissimilar habitats promotes local adaptation, which can lead to reproductive isolation (RI). Populations in the Poecilia mexicana species complex have independently adapted to toxic hydrogen sulfide and show varying degrees of RI. Here, we examined the variation in the mate choice component of prezygotic RI. Mate choice tests across drainages (with stimulus males from another drainage) suggest that specific features of the males coupled with a general female preference for yellow color patterns explain the observed variation. Analyses of male body coloration identified the intensity of yellow fin coloration as a strong candidate to explain this pattern, and common-garden rearing suggested heritable population differences. Male sexual ornamentation apparently evolved differently across sulfide-adapted populations, for example because of differences in natural counterselection via predation. The ubiquitous preference for yellow color ornaments in poeciliid females likely undermines the emergence of strong RI, as female discrimination in favor of own males becomes weaker when yellow fin coloration in the respective sulfide ecotype increases. Our study illustrates the complexity of the (partly non-parallel) pathways to divergence among replicated ecological gradients. We suggest that future work should identify the genomic loci involved in the pattern reported here, making use of the increasing genomic and transcriptomic datasets available for our study system.