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.

Influence of host plant, geography and pheromone strain on genomic differentiation in sympatric populations of Ostrinia nubilalis.

Patterns of mating for the European corn borer (Ostrinia nubilalis) moth depend in part on variation in sex-pheromone blend. The ratio of (E)-11- and (Z)-11-tetradecenyl acetate (E11- and Z11-14:OAc) in the pheromone blend that females produce and males respond to differs between strains of O. nubilalis. Populations also vary in female oviposition preference for and larval performance on maize (C4) and nonmaize (C3) host plants. The relative contributions of sexual and ecological trait variation to the genetic structure of O. nubilalis remains unknown. Host-plant use (13 C/14 C ratios) and genetic differentiation were estimated among sympatric E and Z pheromone strain O. nubilalis males collected in sex-pheromone baited traps at 12 locations in Pennsylvania and New York between 2007 and 2010. Among genotypes at 65 single nucleotide polymorphism marker loci, variance at a position in the pheromone gland fatty acyl-reductase (pgfar) gene at the locus responsible for determining female pheromone ratio (Pher) explained 64% of the total genetic differentiation between males attracted to different pheromones (male response, Resp), providing evidence of sexual inter-selection at these unlinked loci. Principal coordinate, Bayesian clustering, and distance-based redundancy analysis (dbRDA) demonstrate that host plant history or geography does not significantly contribute to population variation or differentiation among males. In contrast, these analyses indicate that pheromone response and pgfar-defined strain contribute significantly to population genetic differentiation. This study suggests that behavioural divergence probably plays a larger role in driving genetic variation compared to host plant-defined ecological adaptation.

A combination of sexual and ecological divergence contributes to rearrangement spread during initial stages of speciation.

Chromosomal rearrangements between sympatric species often contain multiple loci contributing to assortative mating, local adaptation and hybrid sterility. When and how these associations arise during the process of speciation remains a subject of debate. Here, we address the relative roles of local adaptation and assortative mating on the dynamics of rearrangement evolution by studying how a rearrangement covaries with sexual and ecological trait divergence within a species. Previously, a chromosomal rearrangement that suppresses recombination on the Z (sex) chromosome was identified in European corn borer moths (Ostrinia nubilalis). We further characterize this recombination suppressor and explore its association with variation in sex pheromone communication and seasonal ecological adaptation in pairs of populations that are divergent in one or both of these characteristics. Direct estimates of recombination suppression in pedigree mapping families indicated that more than 39% of the Z chromosome (encompassing up to ~10 megabases and ~300 genes) resides within a nonrecombining unit, including pheromone olfactory receptor genes and a major quantitative trait locus that contributes to ecotype differences (Pdd). Combining direct and indirect estimates of recombination suppression, we found that the rearrangement was occasionally present between sexually isolated strains (E vs. Z) and between divergent ecotypes (univoltine vs. bivoltine). However, it was only consistently present when populations differed in both sexual and ecological traits. Our results suggest that independent of the forces that drove the initial establishment of the rearrangement, a combination of sexual and ecological divergence is required for rearrangement spread during speciation.

Behavioral Isolation due to Cascade Reinforcement in Lucania Killifish.

Reinforcement occurs when behavioral isolation is strengthened between species due to selection against hybridization in sympatry. Mate preferences and their target traits may change in sympatry as a consequence of reinforcement. This can potentially generate further behavioral isolation within species if sympatric populations evolve extreme preferences or traits that cause them to reject individuals from foreign populations as mates or be rejected as mates. This process is known as cascade reinforcement. We measured behavioral isolation between sympatric and allopatric populations of Lucania killifish to determine whether isolation evolves due to reinforcement between species and whether reinforcement affects preferences within species, consistent with the cascade reinforcement hypothesis. We measured mate preferences in both sexes between species (Lucania parva vs. Lucania goodei) and within species (among populations of L. parva). Between species, both male and female preferences for conspecifics were highest in sympatric populations. Within species, L. parva females from sympatric populations preferred their own native males over foreign males. Allopatric L. parva females and all L. parva males showed no preferences within species. Within species, behavioral isolation showed no association with ecological variables, such as salinity. Thus, reinforcement is a primary factor generating behavioral isolation in Lucania killifish, creating strong preferences in both sexes among species and leading to cascade reinforcement of female mate preference within species.

Predator experience overrides learned aversion to heterospecifics in stickleback species pairs.

Predation risk can alter female mating decisions because the costs of mate searching and selecting attractive mates increase when predators are present. In response to predators, females have been found to plastically adjust mate preference within species, but little is known about how predators alter sexual isolation and hybridization among species. We tested the effects of predator exposure on sexual isolation between benthic and limnetic threespine sticklebacks (Gasterosteus spp.). Female discrimination against heterospecific mates was measured before and after females experienced a simulated attack by a trout predator or a control exposure to a harmless object. In the absence of predators, females showed increased aversion to heterospecifics over time. We found that predator exposure made females less discriminating and precluded this learned aversion to heterospecifics. Benthic and limnetic males differ in coloration, and predator exposure also affected sexual isolation by weakening female preferences for colourful males. Predator effects on sexual selection were also tested but predators had few effects on female choosiness among conspecific mates. Our results suggest that predation risk may disrupt the cognitive processes associated with mate choice and lead to fluctuations in the strength of sexual isolation between species.

Temperature impacts how sugar resources alter reproductive investment in the European corn borer moth.

Investment of resources in reproduction can be based on individual state, environmental conditions, and perceived mate quality. Changing climates impact many aspects of the environment by increasing temperature, decreasing precipitation, and altering resource availability. Access to high-quality resources is known to improve survival under elevated temperatures, but its effects on reproduction in warming environments are largely unexplored. Here, we investigate the effects of elevated temperature and sugar resources on reproductive output within and between E- and Z-pheromone strains of the European corn borer moth [Ostrinia nubilalis (Hübner) (Lepidoptera: Crambidae)]. Corn borers prefer mates from their own strain, with reproductive output being highest for within-strain pairs. In this experiment, mating pairs were provided with a 20% sugar solution while exposed to either ambient (23 °C) or elevated (28 °C) temperatures. We measured reproductive investment as the total number of egg clusters laid 3 days after pairing. We found that at ambient temperature, sugar supplementation resulted in high investment across all pairs, including with usually unpreferred mates. However, when sugar was provided at elevated temperature, more egg clusters were laid in pairs with preferred (within-strain) mates as compared to less preferred (between-strain) mates. These results differ from temperature effects in the absence of sugar and suggest that the effects of sugar on reproductive investment in less preferred mates depend on temperature. Changes in investment may be due to differences in the allocation of extra resources to thermoregulation at elevated temperatures. Our results suggest the possibility of interactive effects of temperature and resources on sexual selection.

Plastic responses to parents and predators lead to divergent shoaling behaviour in sticklebacks.

Population divergence in antipredator defence and behaviour occurs rapidly and repeatedly. Genetic differences, phenotypic plasticity or parental effects may all contribute to divergence, but the relative importance of each of these mechanisms remains unknown. We exposed juveniles to parents and predators to measure how induced changes contribute to shoaling behaviour differences between two threespine stickleback species (benthics and limnetics: Gasterosteus spp). We found that limnetics increased shoaling in response to predator attacks, whereas benthics did not alter their behaviour. Care by limnetic fathers led to increased shoaling in both limnetic and benthic offspring. Shoaling helps limnetics avoid trout and avian predation; our results suggest that this adaptive behaviour is the result of a combination of paternal effects, predator-induced plasticity and genetic differences between species. These results suggest that plasticity substantially contributes to the rapid divergence in shoaling behaviour across the post-Pleistocene radiation of sticklebacks.

Genomic architecture of supergenes: connecting form and function.

Supergenes are tightly linked sets of loci that are inherited together and control complex phenotypes. While classical supergenes-governing traits such as wing patterns in Heliconius butterflies or heterostyly in Primula-have been studied since the Modern Synthesis, we still understand very little about how they evolve and persist in nature. The genetic architecture of supergenes is a critical factor affecting their evolutionary fate, as it can change key parameters such as recombination rate and effective population size, potentially redirecting molecular evolution of the supergene in addition to the surrounding genomic region. To understand supergene evolution, we must link genomic architecture with evolutionary patterns and processes. This is now becoming possible with recent advances in sequencing technology and powerful forward computer simulations. The present theme issue brings together theoretical and empirical papers, as well as opinion and synthesis papers, which showcase the architectural diversity of supergenes and connect this to critical processes in supergene evolution, such as polymorphism maintenance and mutation accumulation. Here, we summarize those insights to highlight new ideas and methods that illuminate the path forward for the study of supergenes in nature. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Consequences of coupled barriers to gene flow for the build-up of genomic differentiation.

Theory predicts that when different barriers to gene flow become coincident, their joint effects enhance reproductive isolation and genomic divergence beyond their individual effects, but empirical tests of this "coupling" hypothesis are rare. Here, we analyze patterns of gene exchange among populations of European corn borer moths that vary in the number of acting barriers, allowing for comparisons of genomic variation when barrier traits or loci are in coincident or independent states. We find that divergence is mainly restricted to barrier loci when populations differ by a single barrier, whereas the coincidence of temporal and behavioral barriers is associated with divergence of two chromosomes harboring barrier loci. Furthermore, differentiation at temporal barrier loci increases in the presence of behavioral divergence and differentiation at behavioral barrier loci increases in the presence of temporal divergence. Our results demonstrate how the joint action of coincident barrier effects leads to levels of genomic differentiation that far exceed those of single barriers acting alone, consistent with theory arguing that coupling allows indirect selection to combine with direct selection and thereby lead to a stronger overall barrier to gene flow. Thus, the state of barriers-independent or coupled-strongly influences the accumulation of genomic differentiation.

Sexual imprinting on ecologically divergent traits leads to sexual isolation in sticklebacks.

During sexual imprinting, offspring learn parental phenotypes and then select mates who are similar to their parents. Imprinting has been thought to contribute to the process of speciation in only a few rare cases; this is despite imprinting's potential to generate assortative mating and solve the problem of recombination in ecological speciation. If offspring imprint on parental traits under divergent selection, these traits will then be involved in both adaptation and mate preference. Such 'magic traits' easily generate sexual isolation and facilitate speciation. In this study, we show that imprinting occurs in two ecologically divergent stickleback species (benthics and limnetics: Gasterosteus spp.). Cross-fostered females preferred mates of their foster father's species. Furthermore, imprinting is essential for sexual isolation between species; isolation was reduced when females were raised without fathers. Daughters imprinted on father odour and colour during a critical period early in development. These traits have diverged between the species owing to differences in ecology. Therefore, we provide the first evidence that imprinting links ecological adaptation to sexual isolation between species. Our results suggest that imprinting may facilitate the evolution of sexual isolation during ecological speciation, may be especially important in cases of rapid diversification, and thus play an integral role in the generation of biodiversity.

Elevated temperature increases reproductive investment in less preferred mates in the invasive European corn borer moth.

Rapidly changing environments may weaken sexual selection and lead to indiscriminate mating by interfering with the reception of mating signals or by increasing the costs associated with mate choice. If temperature alters sexual selection, it may impact population response and adaptation to climate change. Here, we examine how differences in temperature of the mating environment influence reproductive investment in the European corn borer moth (Ostrinia nubilalis). Mate preference in this species is known to be related to pheromone usage, with assortative mating occurring between genetically distinct E and Z strains that differ in the composition of female and male pheromones. We compared egg production within and between corn borer lines derived from four different populations that vary in pheromone composition and other traits. Pairs of adults were placed in a mating environment that matched the pupal environment (ambient temperature) or at elevated temperature (5°C above the pupal environment). At ambient temperature, we found that within-line pairs produced eggs sooner and produced more egg clusters than between-line pairs. However, at elevated temperature, between-line pairs produced the same number of egg clusters as within-line pairs. These results suggest that elevated temperature increased investment in matings with typically less preferred, between-line mates. This increased investment could result in changes in gene flow among corn borer populations in warming environments.

bric à brac controls sex pheromone choice by male European corn borer moths.

The sex pheromone system of ~160,000 moth species acts as a powerful form of assortative mating whereby females attract conspecific males with a species-specific blend of volatile compounds. Understanding how female pheromone production and male preference coevolve to produce this diversity requires knowledge of the genes underlying change in both traits. In the European corn borer moth, pheromone blend variation is controlled by two alleles of an autosomal fatty-acyl reductase gene expressed in the female pheromone gland (pgFAR). Here we show that asymmetric male preference is controlled by cis-acting variation in a sex-linked transcription factor expressed in the developing male antenna, bric à brac (bab). A genome-wide association study of preference using pheromone-trapped males implicates variation in the 293 kb bab intron 1, rather than the coding sequence. Linkage disequilibrium between bab intron 1 and pgFAR further validates bab as the preference locus, and demonstrates that the two genes interact to contribute to assortative mating. Thus, lack of physical linkage is not a constraint for coevolutionary divergence of female pheromone production and male behavioral response genes, in contrast to what is often predicted by evolutionary theory.

Genomic Basis of Circannual Rhythm in the European Corn Borer Moth.

Synchronizing the annual timing of physiological, morphological, and behavioral transitions with seasons enables survival in temperate environments [1]. The capacity to adjust life history timing and track local seasonal cycles can facilitate geographic expansion [2], adaptation [3], and tolerance [4-6] during rapid environmental change. Understanding the proximate causes of variation in seasonal timing improves prediction of future response and persistence [7, 8]. However, relatively little is known about the molecular basis generating this diversity [9], particularly in Lepidoptera, a group with many species in decline [10, 11]. In insects, the stress-tolerant physiological state of diapause enables coping with seasonal challenges [1, 12-15]. Seasonal changes in photoperiod and temperature are used to synchronize diapause with winter, and timing of diapause transitions varies widely within and among species [9, 16]. Changes in spring diapause termination in the European corn borer moth (Ostrinia nubilalis) have allowed populations to respond to shorter winters and emerge ∼3 weeks earlier in the year [17]. Multiple whole-genome approaches suggest two circadian clock genes, period (per) and pigment-dispersing factor receptor (Pdfr), underlie this polymorphism. Per and Pdfr are within interacting quantitative trait loci (QTL) and differ in allele frequency among individuals that end diapause early or late, with alleles maintained in high linkage disequilibrium. Our results provide testable hypotheses about the physiological role of circadian clock genes in the circannual timer. We predict these gene candidates will be targets of selection for future adaptation under continued global climate change [18].

Non-Pleiotropic Coupling of Daily and Seasonal Temporal Isolation in the European Corn Borer.

Speciation often involves the coupling of multiple isolating barriers to produce reproductive isolation, but how coupling is generated among different premating barriers is unknown. We measure the degree of coupling between the daily mating time and seasonal mating time between strains of European corn borer (Ostrinia nubilalis) and evaluate the hypothesis that the coupling of different forms of allochrony is due to a shared genetic architecture, involving genes with pleiotropic effects on both timing phenotypes. We measure differences in gene expression at peak mating times and compare these genes to previously identified candidates that are associated with changes in seasonal mating time between the corn borer strains. We find that the E strain, which mates earlier in the season, also mates 2.7 h earlier in the night than the Z strain. Earlier daily mating is correlated with the differences in expression of the circadian clock genes cycle, slimb, and vrille. However, different circadian clock genes associate with daily and seasonal timing, suggesting that the coupling of timing traits is maintained by natural selection rather than pleiotropy. Juvenile hormone gene expression was associated with both types of timing, suggesting that circadian genes activate common downstream modules that may impose constraint on future evolution of these traits.

Learning to speciate: The biased learning of mate preferences promotes adaptive radiation.

Bursts of rapid repeated speciation called adaptive radiations have generated much of Earth's biodiversity and fascinated biologists since Darwin, but we still do not know why some lineages radiate and others do not. Understanding what causes assortative mating to evolve rapidly and repeatedly in the same lineage is key to understanding adaptive radiation. Many species that have undergone adaptive radiations exhibit mate preference learning, where individuals acquire mate preferences by observing the phenotypes of other members of their populations. Mate preference learning can be biased if individuals also learn phenotypes to avoid in mates, and shift their preferences away from these avoided phenotypes. We used individual-based computational simulations to study whether biased and unbiased mate preference learning promotes ecological speciation and adaptive radiation. We found that ecological speciation can be rapid and repeated when mate preferences are biased, but is inhibited when mate preferences are learned without bias. Our results suggest that biased mate preference learning may play an important role in generating animal biodiversity through adaptive radiation.

Insight into genomic changes accompanying divergence: genetic linkage maps and synteny of Lucania goodei and L. parva reveal a Robertsonian fusion.

Linkage maps are important tools in evolutionary genetics and in studies of speciation. We performed a karyotyping study and constructed high-density linkage maps for two closely related killifish species, Lucania parva and L. goodei, that differ in salinity tolerance and still hybridize in their contact zone in Florida. Using SNPs from orthologous EST contigs, we compared synteny between the two species to determine how genomic architecture has shifted with divergence. Karyotyping revealed that L. goodei possesses 24 acrocentric chromosomes (1N) whereas L. parva possesses 23 chromosomes (1N), one of which is a large metacentric chromosome. Likewise, high-density single-nucleotide polymorphism-based linkage maps indicated 24 linkage groups for L. goodei and 23 linkage groups for L. parva. Synteny mapping revealed two linkage groups in L. goodei that were highly syntenic with the largest linkage group in L. parva. Together, this evidence points to the largest linkage group in L. parva being the result of a chromosomal fusion. We further compared synteny between Lucania with the genome of a more distant teleost relative medaka (Oryzias latipes) and found good conservation of synteny at the chromosomal level. Each Lucania LG had a single best match with each medaka chromosome. These results provide the groundwork for future studies on the genetic architecture of reproductive isolation and salinity tolerance in Lucania and other Fundulidae.

Functional and population genomic divergence within and between two species of killifish adapted to different osmotic niches.

Adaptation to salinity affects species distributions, promotes speciation, and guides many evolutionary patterns in fishes. To uncover the basis of a complex trait like osmoregulation, genome-level analyses are sensible. We combine population genomic scans with genome expression profiling to discover candidate genes and pathways associated with divergence between osmotic environments. We compared transcriptome sequence divergence between multiple freshwater and saltwater populations of the rainwater killifish, Lucania parva. We also compared sequence divergence between L. parva and its sister species, Lucania goodei, a freshwater specialist. We found highly differentiated single nucleotide polymorphisms (SNPs) between freshwater and saltwater L. parva populations in cell junction and ion transport genes, including V-type H(+) ATPase. Between species, we found divergence in reproduction and osmotic stress genes. Genes that were differentially expressed between species during osmotic acclimation included genes involved in ion transport and cell volume regulation. Gene sets that were divergent in coding sequence and divergent in expression did not overlap, although they did converge in function. Like many studies using genomic scans, our approach may miss some loci that contribute to adaptation but have complicated patterns of allelic variation. Our study suggests that gene expression and coding sequence may evolve independently as populations adapt to a complex physiological challenge.

Female mate preferences for male body size and shape promote sexual isolation in threespine sticklebacks.

Female mate preferences for ecologically relevant traits may enhance natural selection, leading to rapid divergence. They may also forge a link between mate choice within species and sexual isolation between species. Here, we examine female mate preference for two ecologically important traits: body size and body shape. We measured female preferences within and between species of benthic, limnetic, and anadromous threespine sticklebacks (Gasterosteus aculeatus species complex). We found that mate preferences differed between species and between contexts (i.e., within vs. between species). Within species, anadromous females preferred males that were deep bodied for their size, benthic females preferred larger males (as measured by centroid size), and limnetic females preferred males that were more limnetic shaped. In heterospecific mating trials between benthics and limnetics, limnetic females continued to prefer males that were more limnetic like in shape when presented with benthic males. Benthic females showed no preferences for size when presented with limnetic males. These results show that females use ecologically relevant traits to select mates in all three species and that female preference has diverged between species. These results suggest that sexual selection may act in concert with natural selection on stickleback size and shape. Further, our results suggest that female preferences may track adaptation to local environments and contribute to sexual isolation between benthic and limnetic sticklebacks.

Postzygotic Isolation Evolves before Prezygotic Isolation between Fresh and Saltwater Populations of the Rainwater Killifish, Lucania parva.

Divergent natural selection has the potential to drive the evolution of reproductive isolation. The euryhaline killifish Lucania parva has stable populations in both fresh water and salt water. Lucania parva and its sister species, the freshwater L. goodei, are isolated by both prezygotic and postzygotic barriers. To further test whether adaptation to salinity has led to the evolution of these isolating barriers, we tested for incipient reproductive isolation within L. parva by crossing freshwater and saltwater populations. We found no evidence for prezygotic isolation, but reduced hybrid survival indicated that postzygotic isolation existed between L. parva populations. Therefore, postzygotic isolation evolved before prezygotic isolation in these ecologically divergent populations. Previous work on these species raised eggs with methylene blue, which acts as a fungicide. We found this fungicide distorts the pattern of postzygotic isolation by increasing fresh water survival in L. parva, masking species/population differences, and underestimating hybrid inviability.

The impact of learning on sexual selection and speciation.

Learning is widespread in nature, occurring in most animal taxa and in several different ecological contexts and, thus, might play a key role in evolutionary processes. Here, we review the accumulating empirical evidence for the involvement of learning in mate choice and the consequences for sexual selection and reproductive isolation. We distinguish two broad categories: learned mate preferences and learned traits under mate selection (such as bird song). We point out that the context of learning, namely how and when learning takes place, often makes a crucial difference to the predicted evolutionary outcome. Factors causing biases in learning and when one should expect the evolution of learning itself are also explored.