The molecular basis of phenotypic evolution: beyond the usual suspects.

It has been well documented that mutations in coding DNA or cis-regulatory elements underlie natural phenotypic variation in many organisms. However, the development of sophisticated functional tools in recent years in a wide range of traditionally non-model systems have revealed many 'unusual suspects' in the molecular bases of phenotypic evolution, including upstream open reading frames (uORFs), cryptic splice sites, and small RNAs. Furthermore, large-scale genome sequencing, especially long-read sequencing, has identified a cornucopia of structural variation underlying phenotypic divergence and elucidated the composition of supergenes that control complex multi-trait polymorphisms. In this review article we highlight recent studies that demonstrate this great diversity of molecular mechanisms producing adaptive genetic variation and the panoply of evolutionary paths leading to the 'grandeur of life'.

Heterosis in crop improvement.

Heterosis, also known as hybrid vigor, is the phenomenon wherein a progeny exhibits superior traits relative to one or both parents. In terms of crop breeding, this usually refers to the yield advantage of F1 hybrids over both inbred parents. The development of high-yielding hybrid cultivars across a wider range of crops is key to meeting future food demands. However, conventional hybrid breeding strategies are proving to be exceptionally challenging to apply commercially in many self-pollinating crops, particularly wheat and barley. Currently in these crops, the relative performance advantage of hybrids over inbred line cultivars does not outweigh the cost of hybrid seed production. Here, we review the genetic basis of heterosis, discuss the challenges in hybrid breeding, and propose a strategy to recruit multiple heterosis-associated genes to develop lines with improved agronomic characteristics. This strategy leverages modern genetic engineering tools to synthesize supergenes by fusing multiple heterotic alleles across multiple heterosis-associated loci. We outline a plan to assess the feasibility of this approach to improve line performance using barley (Hordeum vulgare) as the model self-pollinating crop species, and a few heterosis-associated genes. The proposed method can be applied to all crops for which heterotic gene combinations can be identified.

Radiation and hybridization underpin the spread of the fire ant social supergene.

Supergenes are clusters of tightly linked genes that jointly produce complex phenotypes. Although widespread in nature, how such genomic elements are formed and how they spread are in most cases unclear. In the fire ant Solenopsis invicta and closely related species, a "social supergene controls whether a colony maintains one or multiple queens. Here, we show that the three inversions constituting the Social b (Sb) supergene emerged sequentially during the separation of the ancestral lineages of S. invicta and Solenopsis richteri. The two first inversions arose in the ancestral population of both species, while the third one arose in the S. richteri lineage. Once completely assembled in the S. richteri lineage, the supergene first introgressed into S. invicta, and from there into the other species of the socially polymorphic group of South American fire ant species. Surprisingly, the introgression of this large and important genomic element occurred despite recent hybridization being uncommon between several of the species. These results highlight how supergenes can readily move across species boundaries, possibly because of fitness benefits they provide and/or expression of selfish properties favoring their transmission.

Evolutionary genomics of socially polymorphic populations of Pogonomyrmex californicus.

BACKGROUND: Social insects vary considerably in their social organization both between and within species. In the California harvester ant, Pogonomyrmex californicus (Buckley 1867), colonies are commonly founded and headed by a single queen (haplometrosis, primary monogyny). However, in some populations in California (USA), unrelated queens cooperate not only during founding (pleometrosis) but also throughout the life of the colony (primary polygyny). The genetic architecture and evolutionary dynamics of this complex social niche polymorphism (haplometrosis vs pleometrosis) have remained unknown. RESULTS: We provide a first analysis of its genomic basis and evolutionary history using population genomics comparing individuals from a haplometrotic population to those from a pleometrotic population. We discovered a recently evolved (< 200 k years), 8-Mb non-recombining region segregating with the observed social niche polymorphism. This region shares several characteristics with supergenes underlying social polymorphisms in other socially polymorphic ant species. However, we also find remarkable differences from previously described social supergenes. Particularly, four additional genomic regions not in linkage with the supergene show signatures of a selective sweep in the pleometrotic population. Within these regions, we find for example genes crucial for epigenetic regulation via histone modification (chameau) and DNA methylation (Dnmt1). CONCLUSIONS: Altogether, our results suggest that social morph in this species is a polygenic trait involving a potential young supergene. Further studies targeting haplo- and pleometrotic individuals from a single population are however required to conclusively resolve whether these genetic differences underlie the alternative social phenotypes or have emerged through genetic drift.

Low Mutation Load in a Supergene Underpinning Alternative Male Mating Strategies in Ruff (Calidris pugnax).

A paradox in evolutionary biology is how supergenes can maintain high fitness despite reduced effective population size, the suppression of recombination, and the expected accumulation of mutational load. The ruff supergene involves 2 rare inversion haplotypes (satellite and faeder). These are recessive lethals but with dominant effects on male mating strategies, plumage, and body size. Sequence divergence to the wild-type (independent) haplotype indicates that the inversion could be as old as 4 million years. Here, we have constructed a highly contiguous genome assembly of the inversion region for both the independent and satellite haplotypes. Based on the new data, we estimate that the recombination event(s) creating the satellite haplotype occurred only about 70,000 yr ago. Contrary to expectations for supergenes, we find no substantial expansion of repeats and only a modest mutation load on the satellite and faeder haplotypes despite high sequence divergence to the non-inverted haplotype (1.46%). The essential centromere protein N (CENPN) gene is disrupted by the inversion and is as well conserved on the inversion haplotypes as on the noninversion haplotype. These results suggest that the inversion may be much younger than previously thought. The low mutation load, despite recessive lethality, may be explained by the introgression of the inversion from a now extinct lineage.

Recent Evolution of a Maternally Acting Sex-Determining Supergene in a Fly with Single-Sex Broods.

Sex determination is a key developmental process, yet it is remarkably variable across the tree of life. The dipteran family Sciaridae exhibits one of the most unusual sex determination systems in which mothers control offspring sex through selective elimination of paternal X chromosomes. Whereas in some members of the family females produce mixed-sex broods, others such as the dark-winged fungus gnat Bradysia coprophila are monogenic, with females producing single-sex broods. Female-producing females were previously found to be heterozygous for a large X-linked paracentric inversion (X'), which is maternally inherited and absent from male-producing females. Here, we assembled and characterized the X' sequence. As close sequence homology between the X and X' made identification of the inversion challenging, we developed a k-mer-based approach to bin genomic reads before assembly. We confirmed that the inversion spans most of the X' chromosome (∼55 Mb) and encodes ∼3,500 genes. Analysis of the divergence between the inversion and the homologous region of the X revealed that it originated very recently (<0.5 Ma). Surprisingly, we found that the X' is more complex than previously thought and is likely to have undergone multiple rearrangements that have produced regions of varying ages, resembling a supergene composed of evolutionary strata. We found functional degradation of ∼7.3% of genes within the region of recombination suppression, but no evidence of accumulation of repetitive elements. Our findings provide an indication that sex-linked inversions are driving turnover of the strange sex determination system in this family of flies.

The Snapdragon Genomes Reveal the Evolutionary Dynamics of the S-Locus Supergene.

The genus Antirrhinum has been used as a model to study self-incompatibility extensively. The multi-allelic S-locus, carrying a pistil S-RNase and dozens of S-locus F-box (SLF) genes, underlies the genetic control of self-incompatibility (SI) in Antirrhinum hispanicum. However, there have been limited studies on the genomic organization of the S-locus supergene due to a lack of high-quality genomic data. Here, we present the chromosome-level reference and haplotype-resolved genome assemblies of a self-incompatible A. hispanicum line, AhS7S8. For the first time, 2 complete A. hispanicum S-haplotypes spanning ∼1.2 Mb and containing a total of 32 SLFs were reconstructed, whereas most of the SLFs derived from retroelement-mediated proximal or tandem duplication ∼122 Mya. Back then, the S-RNase gene and incipient SLFs came into linkage to form the pro-type of type-1 S-locus in the common ancestor of eudicots. Furthermore, we detected a pleiotropic cis-transcription factor (TF) associated with regulating the expression of SLFs, and two miRNAs may control the expression of this TF. Interspecific S-locus and intraspecific S-haplotype comparisons revealed the dynamic nature and polymorphism of the S-locus supergene mediated by continuous gene duplication, segmental translocation or loss, and TE-mediated transposition events. Our data provide an excellent resource for future research on the evolutionary studies of the S-RNase-based self-incompatibility system.

Hybrid Mimulus flowers attract a new pollinator.

Hybridization is common in flowering plants and is believed to be an important force driving adaptation and speciation. The flowers of hybrids often exhibit new trait combinations, which, theoretically, could attract new species of pollinators. In this study, we found that the hybrids between a hummingbird-pollinated species Mimulus cardinalis and a self-pollinated species Mimulus parishii attract bumblebees (Bombus impatiens), a pollinator not attracted to either of the progenitor species. This novel attraction is explained by new combinations of floral traits in hybrids, including, most importantly, petal color, in addition to nectar concentration and corolla size. To understand how petal color variation is perceived by bumblebees, we performed reflectance spectroscopy and multispectral imaging to model the flower appearance in bee vision. This analysis showed that color variation would impact the ease of detection. We also found that YUP, the genetic locus responsible for a large portion of floral color variation and previously shown to be important in bee interactions with other Mimulus species, also played an important role in this novel attraction. These results together suggest that the attraction of new pollinators to hybrid plants could be an underexplored avenue for pollinator shift and speciation.

The interplay of local adaptation and gene flow may lead to the formation of supergenes.

Supergenes are genetic architectures resulting in the segregation of alternative combinations of alleles underlying complex phenotypes. The co-segregation of alleles at linked loci is often facilitated by polymorphic chromosomal rearrangements suppressing recombination locally. Supergenes are involved in many complex polymorphisms, including sexual, colour or behavioural polymorphisms in numerous plants, fungi, mammals, fish, and insects. Despite a long history of empirical and theoretical research, the formation of supergenes remains poorly understood. Here, using a two-island population genetic model, we explore how gene flow and the evolution of overdominant chromosomal inversions may jointly lead to the formation of supergenes. We show that the evolution of inversions in differentiated populations, both under disruptive selection, leads to an increase in frequency of poorly adapted, immigrant haplotypes. Indeed, rare allelic combinations, such as immigrant haplotypes, are more frequently reshuffled by recombination than common allelic combinations, and therefore benefit from the recombination suppression generated by inversions. When an inversion capturing a locally adapted haplotype spreads but is associated with a fitness cost hampering its fixation (e.g. a recessive mutation load), the maintenance of a non-inverted haplotype in the population is enhanced; under certain conditions, the immigrant haplotype persists alongside the inverted local haplotype, while the standard local haplotype disappears. This establishes a stable, local polymorphism with two non-recombining haplotypes encoding alternative adaptive strategies, that is, a supergene. These results bring new light to the importance of local adaptation, overdominance, and gene flow in the formation of supergenes and inversion polymorphisms in general.

Longer or shorter spines: Reciprocal trait evolution in stickleback via triallelic regulatory changes in Stanniocalcin2a.

Vertebrates have repeatedly modified skeletal structures to adapt to their environments. The threespine stickleback is an excellent system for studying skeletal modifications, as different wild populations have either increased or decreased the lengths of their prominent dorsal and pelvic spines in different freshwater environments. Here we identify a regulatory locus that has a major morphological effect on the length of stickleback dorsal and pelvic spines, which we term Maser (major spine enhancer). Maser maps in a closely linked supergene complex that controls multiple armor, feeding, and behavioral traits on chromosome IV. Natural alleles in Maser are differentiated between marine and freshwater sticklebacks; however, alleles found among freshwater populations are also differentiated, with distinct alleles found in short- and long-spined freshwater populations. The distinct freshwater alleles either increase or decrease expression of the bone growth inhibitor gene Stanniocalcin2a in developing spines, providing a simple genetic mechanism for either increasing or decreasing spine lengths in natural populations. Genomic surveys suggest many recurrently differentiated loci in sticklebacks are similarly specialized into three or more distinct alleles, providing multiple ancient standing variants in particular genes that may contribute to a range of phenotypes in different environments.

High-quality haploid genomes corroborate 29 chromosomes and highly conserved synteny of genes in Hyles hawkmoths (Lepidoptera: Sphingidae).

BACKGROUND: Morphological and traditional genetic studies of the young Pliocene genus Hyles have led to the understanding that despite its importance for taxonomy, phenotypic similarity of wing patterns does not correlate with phylogenetic relationship. To gain insights into various aspects of speciation in the Spurge Hawkmoth (Hyles euphorbiae), we assembled a chromosome-level genome and investigated some of its characteristics. RESULTS: The genome of a male H. euphorbiae was sequenced using PacBio and Hi-C data, yielding a 504 Mb assembly (scaffold N50 of 18.2 Mb) with 99.9% of data represented by the 29 largest scaffolds forming the haploid chromosome set. Consistent with this, FISH analysis of the karyotype revealed n = 29 chromosomes and a WZ/ZZ (female/male) sex chromosome system. Estimates of chromosome length based on the karyotype image provided an additional quality metric of assembled chromosome size. Rescaffolding the published male H. vespertilio genome resulted in a high-quality assembly (651 Mb, scaffold N50 of 22 Mb) with 98% of sequence data in the 29 chromosomes. The larger genome size of H. vespertilio (average 1C DNA value of 562 Mb) was accompanied by a proportional increase in repeats from 45% in H. euphorbiae (measured as 472 Mb) to almost 55% in H. vespertilio. Several wing pattern genes were found on the same chromosomes in the two species, with varying amounts and positions of repetitive elements and inversions possibly corrupting their function. CONCLUSIONS: Our two-fold comparative genomics approach revealed high gene synteny of the Hyles genomes to other Sphingidae and high correspondence to intact Merian elements, the ancestral linkage groups of Lepidoptera, with the exception of three simple fusion events. We propose a standardized approach for genome taxonomy using nucleotide homology via scaffold chaining as the primary tool combined with Oxford plots based on Merian elements to infer and visualize directionality of chromosomal rearrangements. The identification of wing pattern genes promises future understanding of the evolution of forewing patterns in the genus Hyles, although further sequencing data from more individuals are needed. The genomic data obtained provide additional reliable references for further comparative studies in hawkmoths (Sphingidae).

Comparative Genomics Elucidates the Origin of a Supergene Controlling Floral Heteromorphism.

Supergenes are nonrecombining genomic regions ensuring the coinheritance of multiple, coadapted genes. Despite the importance of supergenes in adaptation, little is known on how they originate. A classic example of supergene is the S locus controlling heterostyly, a floral heteromorphism occurring in 28 angiosperm families. In Primula, heterostyly is characterized by the cooccurrence of two complementary, self-incompatible floral morphs and is controlled by five genes clustered in the hemizygous, ca. 300-kb S locus. Here, we present the first chromosome-scale genome assembly of any heterostylous species, that of Primula veris (cowslip). By leveraging the high contiguity of the P. veris assembly and comparative genomic analyses, we demonstrated that the S-locus evolved via multiple, asynchronous gene duplications and independent gene translocations. Furthermore, we discovered a new whole-genome duplication in Ericales that is specific to the Primula lineage. We also propose a mechanism for the origin of S-locus hemizygosity via nonhomologous recombination involving the newly discovered two pairs of CFB genes flanking the S locus. Finally, we detected only weak signatures of degeneration in the S locus, as predicted for hemizygous supergenes. The present study provides a useful resource for future research addressing key questions on the evolution of supergenes in general and the S locus in particular: How do supergenes arise? What is the role of genome architecture in the evolution of complex adaptations? Is the molecular architecture of heterostyly supergenes across angiosperms similar to that of Primula?

A supergene in seaweed flies modulates male traits and female perception.

Supergenes, tightly linked sets of alleles, offer some of the most spectacular examples of polymorphism persisting under long-term balancing selection. However, we still do not understand their evolution and persistence, especially in the face of accumulation of deleterious elements. Here, we show that an overdominant supergene in seaweed flies, Coelopa frigida, modulates male traits, potentially facilitating disassortative mating and promoting intraspecific polymorphism. Across two continents, the Cf-Inv(1) supergene strongly affected the composition of male cuticular hydrocarbons (CHCs) but only weakly affected CHC composition in females. Using gas chromatography-electroantennographic detection, we show that females can sense male CHCs and that there may be differential perception between genotypes. Combining our phenotypic results with RNA-seq data, we show that candidate genes for CHC biosynthesis primarily show differential expression for Cf-Inv(1) in males but not females. Conversely, candidate genes for odorant detection were differentially expressed in both sexes but showed high levels of divergence between supergene haplotypes. We suggest that the reduced recombination between supergene haplotypes may have led to rapid divergence in mate preferences as well as increasing linkage between male traits, and overdominant loci. Together this probably helped to maintain the polymorphism despite deleterious effects in homozygotes.

Whole-genome analyses disentangle reticulate evolution of primroses in a biodiversity hotspot.

Biodiversity hotspots, such as the Caucasus mountains, provide unprecedented opportunities for understanding the evolutionary processes that shape species diversity and richness. Therefore, we investigated the evolution of Primula sect. Primula, a clade with a high degree of endemism in the Caucasus. We performed phylogenetic and network analyses of whole-genome resequencing data from the entire nuclear genome, the entire chloroplast genome, and the entire heterostyly supergene. The different characteristics of the genomic partitions and the resulting phylogenetic incongruences enabled us to disentangle evolutionary histories resulting from tokogenetic vs cladogenetic processes. We provide the first phylogeny inferred from the heterostyly supergene that includes all species of Primula sect. Primula. Our results identified recurrent admixture at deep nodes between lineages in the Caucasus as the cause of non-monophyly in Primula. Biogeographic analyses support the 'out-of-the-Caucasus' hypothesis, emphasizing the importance of this hotspot as a cradle for biodiversity. Our findings provide novel insights into causal processes of phylogenetic discordance, demonstrating that genome-wide analyses from partitions with contrasting genetic characteristics and broad geographic sampling are crucial for disentangling the diversification of species-rich clades in biodiversity hotspots.

Genomic evidence supports the genetic convergence of a supergene controlling the distylous floral syndrome.

Heterostyly, a plant sexual polymorphism controlled by the S-locus supergene, has evolved numerous times among angiosperm lineages and represents a classic example of convergent evolution in form and function. Determining whether underlying molecular convergence occurs could provide insights on constraints to floral evolution. Here, we investigated S-locus genes in distylous Gelsemium (Gelsemiaceae) to determine whether there is evidence of molecular convergence with unrelated distylous species. We used several approaches, including anatomical measurements of sex-organ development and transcriptome and whole-genome sequencing, to identify components of the S-locus supergene. We also performed evolutionary analysis with candidate S-locus genes and compared them with those reported in Primula and Turnera. The candidate S-locus supergene of Gelsemium contained four genes, of which three appear to have originated from gene duplication events within Gelsemiaceae. The style-length genes GeCYP in Gelsemium and CYP734A50 in Primula likely arose from duplication of the same gene, CYP734A1. Three out of four S-locus genes in Gelsemium elegans were hemizygous, as previously reported in Primula and Turnera. We provide genomic evidence on the genetic convergence of the supergene underlying distyly among distantly related angiosperm lineages and help to illuminate the genetic architecture involved in the evolution of heterostyly.

A terpene synthase supergene locus determines chemotype in Melaleuca alternifolia (tea tree).

Leaf oil terpenes vary categorically in many plant populations, leading to discrete phenotypes of adaptive and economic significance, but for most species, a genetic explanation for the concerted fluctuation in terpene chemistry remains unresolved. To uncover the genetic architecture underlying multi-component terpene chemotypes in Melaleuca alternifolia (tea tree), a genome-wide association study was undertaken for 148 individuals representing all six recognised chemotypes. A number of single nucleotide polymorphisms in a genomic region of c. 400 kb explained large proportions of the variation in key monoterpenes of tea tree oil. The region contained a cluster of 10 monoterpene synthase genes, including four genes predicted to encode synthases for 1,8-cineole, terpinolene, and the terpinen-4-ol precursor, sabinene hydrate. Chemotype-dependent null alleles at some sites suggested structural variants within this gene cluster, providing a possible basis for linkage disequilibrium in this region. Genotyping in a separate domesticated population revealed that all alleles surrounding this gene cluster were fixed after artificial selection for a single chemotype. These observations indicate that a supergene accounts for chemotypes in M. alternifolia. A genetic model with three haplotypes, encompassing the four characterised monoterpene synthase genes, explained the six terpene chemotypes, and was consistent with available biparental cross-segregation data.

Direct and indirect genetic effects of a social supergene.

Indirect genetic effects describe phenotypic variation that results from differences in the genotypic composition of social partners. Such effects represent heritable sources of environmental variation in eusocial organisms because individuals are typically reared by their siblings. In the fire ant Solenopsis invicta, a social supergene exhibits striking indirect genetic effects on worker regulation of colony queen number, such that the genotypic composition of workers at the supergene determines whether colonies contain a single or multiple queens. We assessed the direct and indirect genetic effects of this supergene on gene expression in brains and abdominal tissues from laboratory-reared workers and compared these with previously published data from field-collected prereproductive queens. We found that direct genetic effects caused larger gene expression changes and were more consistent across tissue types and castes than indirect genetic effects. Indirect genetic effects influenced the expression of many loci but were generally restricted to the abdominal tissues. Further, indirect genetic effects were only detected when the genotypic composition of social partners differed throughout the development and adult life of focal workers, and were often only significant with relatively lenient statistical cutoffs. Our study provides insight into direct and indirect genetic effects of a social supergene on gene regulatory dynamics across tissues and castes in a complex society.

Supergene evolution via stepwise duplications and neofunctionalization of a floral-organ identity gene.

Heterostyly represents a fascinating adaptation to promote outbreeding in plants that evolved multiple times independently. While l-morph individuals form flowers with long styles, short anthers, and small pollen grains, S-morph individuals have flowers with short styles, long anthers, and large pollen grains. The difference between the morphs is controlled by an S-locus "supergene" consisting of several distinct genes that determine different traits of the syndrome and are held together, because recombination between them is suppressed. In Primula, the S locus is a roughly 300-kb hemizygous region containing five predicted genes. However, with one exception, their roles remain unclear, as does the evolutionary buildup of the S locus. Here we demonstrate that the MADS-box GLOBOSA2 (GLO2) gene at the S locus determines anther position. In Primula forbesii S-morph plants, GLO2 promotes growth by cell expansion in the fused tube of petals and stamen filaments beneath the anther insertion point; by contrast, neither pollen size nor male incompatibility is affected by GLO2 activity. The paralogue GLO1, from which GLO2 arose by duplication, has maintained the ancestral B-class function in specifying petal and stamen identity, indicating that GLO2 underwent neofunctionalization, likely at the level of the encoded protein. Genetic mapping and phylogenetic analysis indicate that the duplications giving rise to the style-length-determining gene CYP734A50 and to GLO2 occurred sequentially, with the CYP734A50 duplication likely the first. Together these results provide the most detailed insight into the assembly of a plant supergene yet and have important implications for the evolution of heterostyly.

Maternal effect killing by a supergene controlling ant social organization.

Supergenes underlie striking polymorphisms in nature, yet the evolutionary mechanisms by which they arise and persist remain enigmatic. These clusters of linked loci can spread in populations because they captured coadapted alleles or by selfishly distorting the laws of Mendelian inheritance. Here, we show that the supergene haplotype associated with multiple-queen colonies in Alpine silver ants is a maternal effect killer. All eggs from heterozygous queens failed to hatch when they did not inherit this haplotype. Hence, the haplotype specific to multiple-queen colonies is a selfish genetic element that enhances its own transmission by causing developmental arrest of progeny that do not carry it. At the population level, such transmission ratio distortion favors the spread of multiple-queen colonies, to the detriment of the alternative haplotype associated with single-queen colonies. Hence, selfish gene drive by one haplotype will impact the evolutionary dynamics of alternative forms of colony social organization. This killer hidden in a social supergene shows that large nonrecombining genomic regions are prone to cause multifarious effects across levels of biological organization.

Characterization of Queen Supergene Pheromone in the Red Imported Fire Ant Using Worker Discrimination Assays.

Ants use chemical signals to communicate for various purposes related to colony function. Social organization in the red imported fire ant, Solenopsis invicta, is determined by the Sb supergene, with colonies of the monogyne (single-queen) form lacking the element and colonies of the polygyne (multiple-queen) form possessing it. Polygyne workers accept new reproductive queens in their nest, but only those carrying Sb; young winged queens lacking this genetic element are executed as they mature sexually in their natal nest or as they attempt to enter a foreign nest to initiate reproduction after mating and shedding their wings. It has been suggested that queen supergene genotype status is signaled to workers by unsaturated cuticular hydrocarbons, while queen reproductive status is signaled by piperidines (venom alkaloids). We used high-throughput behavioral assays to study worker acceptance of paper dummies dosed with fractions of extracts of polygyne queens, or blends of synthetic counterparts of queen cuticular compounds. We show that the queen supergene pheromone comprises a blend of monoene and diene unsaturated hydrocarbons. Our assays also reveal that unsaturated hydrocarbons elicit discrimination by polygyne workers only when associated with additional compounds that signal queen fertility. This synergistic effect was obtained with a polar fraction of queen extracts, but not by the piperidine alkaloids, suggesting that the chemical(s) indicating queen reproductive status are compounds more polar than cuticular hydrocarbons but are not the piperidine alkaloids. Our results advance understanding of the role of chemical signaling that is central to the regulation of social organization in an important invasive pest and model ant species.