Mate Choice by Both Sexes Maintains Reproductive Isolation in a Species Flock of Pupfish (Cyprinodonspp) in the Bahamas

Seasonal stability and species specificity of environmentally acquired chemical mating signals in orchid bees

Traits that mediate reproductive isolation between species, such as those involved in mate choice and/or recognition, are predicted to experience stabilizing selection towards the species mean. Male orchid bees collect chemical compounds from many sources, such as plants and fungi, which they use as a perfume signal (pheromone) during courtship display, and are suggested to contribute to reproductive isolation between species. Environmentally acquired signals are more prone to variation as source availability can vary through space and time. If orchid bee perfumes are important for reproductive isolation between species, we expect them to exhibit stable species-specific differences in time and space. Here, we describe phenotypic patterns of inter- and intraspecific variation in the male perfumes of three sympatric species of Euglossa orchid bees across an entire year, investigating both their seasonality and species specificity. Our analysis revealed considerable within-species variation in perfumes. However, species specificity was maintained consistently throughout the year, supporting the idea that these perfumes could play an important role in reproductive isolation and are experiencing stabilizing selection towards a species mean. Our analysis also identified strong correlations in the abundance of some compounds, possibly due to shared collection sources between species. Our study suggests that orchid bee perfumes are robust in the face of environmental changes in resource availability and thus can maintain reproductive isolation between species.

Few Fixed Variants between Trophic Specialist Pupfish Species Reveal Candidate Cis-Regulatory Alleles Underlying Rapid Craniofacial Divergence

Investigating closely related species that rapidly evolved divergent feeding morphology is a powerful approach to identify genetic variation underlying variation in complex traits. This can also lead to the discovery of novel candidate genes influencing natural and clinical variation in human craniofacial phenotypes. We combined whole-genome resequencing of 258 individuals with 50 transcriptomes to identify candidate cis-acting genetic variation underlying rapidly evolving craniofacial phenotypes within an adaptive radiation of Cyprinodon pupfishes. This radiation consists of a dietary generalist species and two derived trophic niche specialists-a molluscivore and a scale-eating species. Despite extensive morphological divergence, these species only diverged 10 kya and produce fertile hybrids in the laboratory. Out of 9.3 million genome-wide SNPs and 80,012 structural variants, we found very few alleles fixed between species-only 157 SNPs and 87 deletions. Comparing gene expression across 38 purebred F1 offspring sampled at three early developmental stages, we identified 17 fixed variants within 10 kb of 12 genes that were highly differentially expressed between species. By measuring allele-specific expression in F1 hybrids from multiple crosses, we found that the majority of expression divergence between species was explained by trans-regulatory mechanisms. We also found strong evidence for two cis-regulatory alleles affecting expression divergence of two genes with putative effects on skeletal development (dync2li1 and pycr3). These results suggest that SNPs and structural variants contribute to the evolution of novel traits and highlight the utility of the San Salvador Island pupfish system as an evolutionary model for craniofacial development.

Asymmetric reinforcement in Lucania killifish: assessing reproductive isolation when both sexes choose

Reinforcement can occur when maladaptive hybridization in sympatry favors the evolution of conspecific preferences and target traits that promote behavioral isolation (BI). In many systems, enhanced BI is due to increased female preference for conspecifics. In others, BI is driven by male preference, and in other systems both sexes exert preferences. Some of these patterns can be attributed to classic sex-specific costs and benefits of preference. Alternatively, sex differences in conspecific preference can emerge due to asymmetric postzygotic isolation (e.g., hybrid offspring from female A × male B have lower fitness than hybrid offspring from female B × male A), which can lead to asymmetric BI (e.g., female A and male B are less likely to mate than female B and male A). Understanding reinforcement requires understanding how conspecific preferences evolve in sympatry. Yet, estimating conspecific preferences can be difficult when both sexes are choosy. In this study, we use Lucania killifish to test the hypothesis that patterns of reinforcement are driven by asymmetric postzygotic isolation between species. If true, we predicted that sympatric female Lucania goodei and sympatric male L. parva should have lower levels of BI compared with their sympatric counterparts, as they produce hybrid offspring with the highest fitness. To address the problem of measuring BI when both sexes are choosy, we inferred the contribution to BI of each partner using assays where one sex in the mating pair comes from an allopatric population with potentially low preference, whereas the other comes from a sympatric population with high preference. For one hybrid cross direction, we found that both female L. parva and male L. goodei have high contributions to BI in sympatry. In the other hybrid cross direction, we found that only female L. goodei contribute to BI. Sympatric male L. parva readily engaged in hybrid spawnings with allopatric L. goodei females. These results indicate that both asymmetric postzygotic isolation and the traditional sex-specific costs to preference likely affect the nature of selection on conspecific preferences and target traits.

Ecological divergence in sympatry causes gene misexpression in hybrids

Ecological speciation occurs when reproductive isolation evolves as a byproduct of adaptive divergence between populations. Selection favouring gene regulatory divergence between species could result in transgressive levels of gene expression in F1 hybrids that may lower hybrid fitness. We combined 58 resequenced genomes with 124 transcriptomes to identify patterns of hybrid gene misexpression that may be driven by adaptive regulatory divergence within a young radiation of Cyprinodon pupfishes, which consists of a dietary generalist and two trophic specialists-a molluscivore and a scale-eater. We found more differential gene expression between closely related sympatric specialists than between allopatric generalist populations separated by 1,000 km. Intriguingly, 9.6% of genes that were differentially expressed between sympatric species were also misexpressed in F1 hybrids. A subset of these genes were in highly differentiated genomic regions and enriched for functions important for trophic specialization, including head, muscle and brain development. These regions also included genes that showed evidence of hard selective sweeps and were significantly associated with oral jaw length-the most rapidly diversifying skeletal trait in this radiation. Our results indicate that divergent ecological selection in sympatry can contribute to hybrid gene misexpression which may act as a reproductive barrier between nascent species.

Coevolution of male and female mate choice can destabilize reproductive isolation

Sexual interactions play an important role in the evolution of reproductive isolation, with important consequences for speciation. Theoretical studies have focused on the evolution of mate preferences in each sex separately. However, mounting empirical evidence suggests that premating isolation often involves mutual mate choice. Here, using a population genetic model, we investigate how female and male mate choice coevolve under a phenotype matching rule and how this affects reproductive isolation. We show that the evolution of female preferences increases the mating success of males with reciprocal preferences, favouring mutual mate choice. However, the evolution of male preferences weakens indirect selection on female preferences and, with weak genetic drift, the coevolution of female and male mate choice leads to periodic episodes of random mating with increased hybridization (deterministic 'preference cycling' triggered by stochasticity). Thus, counterintuitively, the process of establishing premating isolation proves rather fragile if both male and female mate choice contribute to assortative mating.

John ME. How to Investigate the Origins of Novelty: Insights Gained from Genetic, Behavioral, and Fitness Perspectives

Biologists are drawn to the most extraordinary adaptations in the natural world, often referred to as evolutionary novelties, yet rarely do we understand the microevolutionary context underlying the origins of novel traits, behaviors, or ecological niches. Here we discuss insights gained into the origins of novelty from a research program spanning biological levels of organization from genotype to fitness in Caribbean pupfishes. We focus on a case study of the origins of novel trophic specialists on San Salvador Island, Bahamas and place this radiation in the context of other rapid radiations. We highlight questions that can be addressed about the origins of novelty at different biological levels, such as measuring the isolation of novel phenotypes on the fitness landscape, locating the spatial and temporal origins of adaptive variation contributing to novelty, detecting dysfunctional gene regulation due to adaptive divergence, and connecting behaviors with novel traits. Evolutionary novelties are rare, almost by definition, and we conclude that integrative case studies can provide insights into this rarity relative to the dynamics of adaptation to more common ecological niches and repeated parallel speciation, such as the relative isolation of novel phenotypes on fitness landscapes and the transient availability of ecological, genetic, and behavioral opportunities.

Hybrid gene misregulation in multiple developing tissues within a recent adaptive radiation of Cyprinodon pupfishes

Genetic incompatibilities constitute the final stages of reproductive isolation and speciation, but little is known about incompatibilities that occur within recent adaptive radiations among closely related diverging populations. Crossing divergent species to form hybrids can break up coadapted variation, resulting in genetic incompatibilities within developmental networks shaping divergent adaptive traits. We crossed two closely related sympatric Cyprinodon pupfish species–a dietary generalist and a specialized molluscivore–and measured expression levels in their F1 hybrids to identify regulatory variation underlying the novel craniofacial morphology found in this recent microendemic adaptive radiation. We extracted mRNA from eight day old whole-larvae tissue and from craniofacial tissues dissected from 17–20 day old larvae to compare gene expression between a total of seven F1 hybrids and 24 individuals from parental species populations. We found 3.9% of genes differentially expressed between generalists and molluscivores in whole-larvae tissues and 0.6% of genes differentially expressed in craniofacial tissue. We found that 2.1% of genes were misregulated in whole-larvae hybrids whereas 19.1% of genes were misregulated in hybrid craniofacial tissues, after correcting for sequencing biases. We also measured allele specific expression across 15,429 heterozygous sites to identify putative compensatory regulatory mechanisms underlying differential expression between generalists and molluscivores. Together, our results highlight the importance of considering misregulation as an early indicator of genetic incompatibilities in the context of rapidly diverging adaptive radiations and suggests that compensatory regulatory divergence drives hybrid gene misregulation in developing tissues that give rise to novel craniofacial traits.

An embryonic staging series up to hatching for Cyprinodon variegatus: An emerging fish model for developmental, evolutionary, and ecological research

Using multiple taxa to research development is necessary for making general conclusions about developmental patterns and mechanisms. We present a staging series for Cyprinodon variegatus as a basis for further study of the developmental biology of fishes in the genus Cyprinodon and for comparative work on teleost fishes beyond the standard models. Cyprinodon are small, euryhaline fishes, widely distributed in fresh, brackish, and hypersaline waters of southern and eastern North America. Cyprinodontids are closely related to fundulids, providing a comparative reference point to the embryological model, Fundulus heteroclitus. Ecologists and evolutionary biologists commonly study Cyprinodon, and we have been using Cyprinodon to study skull variation and its genetic basis among closely related species. We divided embryonic development of C. variegatus into 34 morphologically identifiable stages. We reference our staging series to that already defined for a related model species, Oryzias latipes (medaka) that is studied by a large community of researchers. We provide a description of the early chondrogenesis and ossification of skull and caudal fin bones during the latter stages of embryonic development. We show that Cyprinodon are tractable for studying development. Eggs can be obtained easily from breeding pairs and our study provides a staging system to facilitate future developmental studies.

Sex roles and sexual selection: lessons from a dynamic model system

Our understanding of sexual selection has greatly improved during the last decades. The focus is no longer solely on males, but also on how female competition and male mate choice shape ornamentation and other sexually selected traits in females. At the same time, the focus has shifted from documenting sexual selection to exploring variation and spatiotemporal dynamics of sexual selection, and their evolutionary consequences. Here, I review insights from a model system with exceptionally dynamic sexual selection, the two-spotted goby fish Gobiusculus flavescens. The species displays a complete reversal of sex roles over a 3-month breeding season. The reversal is driven by a dramatic change in the operational sex ratio, which is heavily male-biased at the start of the season and heavily female-biased late in the season. Early in the season, breeding-ready males outnumber mature females, causing males to be highly competitive, and leading to sexual selection on males. Late in the season, mating-ready females are in excess, engage more in courtship and aggression than males, and rarely reject mating opportunities. With typically many females simultaneously courting available males late in the season, males become selective and prefer more colorful females. This variable sexual selection regime likely explains why both male and female G. flavescens have ornamental colors. The G. flavescens model system reveals that sexual behavior and sexual selection can be astonishingly dynamic in response to short-term fluctuations in mating competition. Future work should explore whether sexual selection is equally dynamic on a spatial scale, and related spatiotemporal dynamics.

The roles of inter- and intra-sexual selection in behavioral isolation between native and invasive pupfishes

Male-male competition and female mate choice may both play important roles in driving and maintaining reproductive isolation between species. When previously allopatric species come into secondary contact with each other due to introductions, they provide an opportunity to evaluate the identity and strength of reproductive isolating mechanisms. If reproductive isolation is not maintained, hybridization may occur. We examined how reproductive isolating mechanisms mediate hybridization between endemic populations of the Red River pupfish Cyprinodon rubrofluviatilis and the recently introduced sheepshead minnow C. variegatus. In lab-based dominance trials, males of both species won the same number of competitions. However, male C. rubrofluviatilis that won competitions were more aggressive than C. variegatus winners, and more aggression was needed to win against competitor C. variagatus than allopatric C. rubrofluviatilis. Duration of fights also differed based on the relatedness of the competitor. In dichotomous mate choice trials, there were no conspecific or heterospecific preferences expressed by females of either species. Our findings that male-male aggression differs between closely and distantly related groups, but female choice does not suggest that male-male competition may be the more likely mechanism to impede gene flow in this system.

Adaptive introgression from distant Caribbean islands contributed to the diversification of a microendemic adaptive radiation of trophic specialist pupfishes

Rapid diversification often involves complex histories of gene flow that leave variable and conflicting signatures of evolutionary relatedness across the genome. Identifying the extent and source of variation in these evolutionary relationships can provide insight into the evolutionary mechanisms involved in rapid radiations. Here we compare the discordant evolutionary relationships associated with species phenotypes across 42 whole genomes from a sympatric adaptive radiation of Cyprinodon pupfishes endemic to San Salvador Island, Bahamas and several outgroup pupfish species in order to understand the rarity of these trophic specialists within the larger radiation of Cyprinodon. 82% of the genome depicts close evolutionary relationships among the San Salvador Island species reflecting their geographic proximity, but the vast majority of variants fixed between specialist species lie in regions with discordant topologies. Top candidate adaptive introgression regions include signatures of selective sweeps and adaptive introgression of genetic variation from a single population in the northwestern Bahamas into each of the specialist species. Hard selective sweeps of genetic variation on San Salvador Island contributed 5 times more to speciation of trophic specialists than adaptive introgression of Caribbean genetic variation; however, four of the 11 introgressed regions came from a single distant island and were associated with the primary axis of oral jaw divergence within the radiation. For example, standing variation in a proto-oncogene (ski) known to have effects on jaw size introgressed into one San Salvador Island specialist from an island 300 km away approximately 10 kya. The complex emerging picture of the origins of adaptive radiation on San Salvador Island indicates that multiple sources of genetic variation contributed to the adaptive phenotypes of novel trophic specialists on the island. Our findings suggest that a suite of factors, including rare adaptive introgression, may be necessary for adaptive radiation in addition to ecological opportunity.

The Cyprinodon variegatus genome reveals gene expression changes underlying differences in skull morphology among closely related species

BACKGROUND

Understanding the genetic and developmental origins of phenotypic novelty is central to the study of biological diversity. In this study we identify modifications to the expression of genes at four developmental stages that may underlie jaw morphological differences among three closely related species of pupfish (genus Cyprinodon) from San Salvador Island, Bahamas. Pupfishes on San Salvador Island are trophically differentiated and include two endemic species that have evolved jaw morphologies unlike that of any other species in the genus Cyprinodon.

RESULTS

We find that gene expression differs significantly across recently diverged species of pupfish. Genes such as Bmp4 and calmodulin, previously implicated in jaw diversification in African cichlid fishes and Galapagos finches, were not found to be differentially expressed among species of pupfish. Instead we find multiple growth factors and cytokine/chemokine genes to be differentially expressed among these pupfish taxa. These include both genes and pathways known to affect craniofacial development, such as Wnt signaling, as well as novel genes and pathways not previously implicated in craniofacial development. These data highlight both shared and potentially unique sources of jaw diversity in pupfish and those identified in other evolutionary model systems such as Galapagos finches and African cichlids.

CONCLUSIONS

We identify modifications to the expression of genes involved in Wnt signaling, Igf signaling, and the inflammation response as promising avenues for future research. Our project provides insight into the magnitude of gene expression changes contributing to the evolution of morphological novelties, such as jaw structure, in recently diverged pupfish species.

Context dependence in complex adaptive landscapes: frequency and trait-dependent selection surfaces within an adaptive radiation of Caribbean pupfishes

The adaptive landscape provides the foundational bridge between micro- and macroevolution. One well-known caveat to this perspective is that fitness surfaces depend on ecological context, including competitor frequency, traits measured, and resource abundance. However, this view is based largely on intraspecific studies. It is still unknown how context-dependence affects the larger features of peaks and valleys on the landscape which ultimately drive speciation and adaptive radiation. Here, I explore this question using one of the most complex fitness landscapes measured in the wild in a sympatric pupfish radiation endemic to San Salvador Island, Bahamas by tracking survival and growth of laboratory-reared F2 hybrids. I present new analyses of the effects of competitor frequency, dietary isotopes, and trait subsets on this fitness landscape. Contrary to expectations, decreasing competitor frequency increased survival only among very common phenotypes, whereas less common phenotypes rarely survived despite few competitors, suggesting that performance, not competitor frequency, shapes large-scale features of the fitness landscape. Dietary isotopes were weakly correlated with phenotype and growth, but did not explain additional survival variation. Nonlinear fitness surfaces varied substantially among trait subsets, revealing one-, two-, and three-peak landscapes, demonstrating the complexity of selection in the wild, even among similar functional traits.