GENE GENEALOGIES REVEAL DIFFERENTIATION AT SEX PHEROMONE OLFACTORY RECEPTOR LOCI IN PHEROMONE STRAINS OF THE EUROPEAN CORN BORER,OSTRINIA NUBILALIS

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.

Pheromone Odorant Receptor Responses Reveal the Presence of a Cryptic, Redundant Sex Pheromone Component in the European Corn Borer, Ostrinia nubilalis

Two odorant receptors (ORs), OnubOR3 and OnubOR6, in the sex pheromone communication systems of E- and Z-strain European corn borers, Ostrinia nubilalis, were broadly receptive to analogs of their pheromone components. In addition to responding to their natural 14-carbon pheromone components, (Z)-11- and (E)-11-tetradecenyl acetates (Z11- and E11-14:OAc), these pheromone ORs responded to the longer-chain compounds, (Z)-11- and (E)-11-hexadecenyl acetate (Z11- and E11-16:OAc). Z11-16:OAc is a pheromone gland constituent of E-strain O. nubilalis females in Europe but has not previously been shown to have behavioral activity to males. Here, we demonstrate that Z11-16:OAc evokes high levels of upwind flight and source location in E-strain males when substituted for Z11-14:OAc (minor component) in the E-strain blend. Since Z11-16:OAc is found in the gland and has behavioral activity when Z11-14:OAc is missing, then it should be classified as a cryptic, redundant minor pheromone component in E-strain O. nubilalis. The opposite geometric isomer, E11-16:OAc, also functions in Z-strain O. nubilalis, substituting behaviorally for the E11-14:OAc minor component, but has not been found in Z-strain female glands. Single-sensillum recordings showed that sensory neurons of E- and Z-strain male antennae expressing OnubOR3 and OnubOR6 produced responses to these hexadecenyl acetates similar to those evoked by the natural (tetradecenyl acetate) pheromone components. We postulate that the wide responsiveness of these two ORs to the 16-carbon acetates could be a preadaptation for O. nubilalis to use these compounds as minor components in lieu of the respective 14-carbon acetates. Alternatively, the responsiveness of OnubOR3 to E11-16:OAc and OnubOR6 to Z11-16:OAc could represent a vestigial state of these receptors, with the 16-carbon acetates having previously acted as functional minor components in an ancestral blend.

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 (¹³ C/¹⁴ 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.

New WGS data and annotation of the heterosomal vs. autosomal localization of Ostrinia scapulalis (Lepidoptera, Crambidae) nuclear genomic scaffolds

Here, we introduce new whole-genome shotgun sequencing and annotation data describing the autosomal vs. Z-heterosomal localization of nuclear genomic scaffolds of the moth species Ostrinia scapulalis. Four WGS libraries (corresponding to 2 males and 2 females) were sequenced with an Illumina HiSeq2500 sequencing technology, and the so-called ‘AD-ratio’ method was applied to distinguish between autosomal and Z-heterosomal scaffolds based on sequencing depth comparisons between homogametic (male) and heterogametic (female) libraries. A total of 25,760 scaffolds (corresponding to 341.69 Mb) were labelled as autosomal and 1273 scaffolds (15.29 Mb) were labelled as Z-heterosomal, totaling about 357 Mb. Besides, 4874 scaffolds (29.07 Mb) remain ambiguous because of a lack of AD-ratio reproducibility between the two replicates. The annotation method was evaluated a posteriori, by comparing depth-based annotation with the exact localization of known genes. Raw genomic data have been deposited and made accessible via the EMBL ENA BioProject id PRJEB26557. Comprehensive annotation is made accessible via the LepidoDB database (http://bipaa.genouest.org/sp/ostrinia_scapulalis/download/genome/v1.2/).

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.

Genetic mapping of male pheromone response in the European corn borer identifies candidate genes regulating neurogenesis

The sexual pheromone communication system of moths is a model system for studies of the evolution of reproductive isolation. Females emit a blend of volatile components that males detect at a distance. Species differences in female pheromone composition and male response directly reinforce reproductive isolation in nature, because even slight variations in the species-specific pheromone blend are usually rejected by the male. The mechanisms by which a new pheromone signal-response system could evolve are enigmatic, because any deviation from the optimally attractive blend should be selected against. Here we investigate the genetic mechanisms enabling a switch in male response. We used a quantitative trait locus-mapping approach to identify the genetic basis of male response in the two pheromone races of the European corn borer, Ostrinia nubilalis Male response to a 99:1 vs. a 3:97 ratio of the E and Z isomers of the female pheromone is governed by a single, sex-linked locus. We found that the chromosomal region most tightly linked to this locus contains genes involved in neurogenesis but, in accordance with an earlier study, does not contain the odorant receptors expressed in the male antenna that detect the pheromone. This finding implies that differences in the development of neuronal pathways conveying information from the antenna, not differences in pheromone detection by the odorant receptors, are primarily responsible for the behavioral response differences among the males in this system. Comparison with other moth species reveals a previously unexplored mechanism by which male pheromone response can change in evolution.

The Genetic Basis of Pheromone Evolution in Moths

Moth sexual pheromones are widely studied as a fine-tuned system of intraspecific sexual communication that reinforces interspecific reproductive isolation. However, their evolution poses a dilemma: How can the female pheromone and male preference simultaneously change to create a new pattern of species-specific attraction? Solving this puzzle requires us to identify the genes underlying intraspecific variation in signals and responses and to understand the evolutionary mechanisms responsible for their interspecific divergence. Candidate gene approaches and functional analyses have yielded insights into large families of biosynthetic enzymes and pheromone receptors, although the factors controlling their expression remain largely unexplored. Intra- and interspecific crosses have provided tantalizing evidence of regulatory genes, although, to date, mapping resolution has been insufficient to identify them. Recent advances in high-throughput genome and transcriptome sequencing, together with established techniques, have great potential to help scientists identify the specific genetic changes underlying divergence and resolve the mystery of how moth sexual communication systems evolve.

Rapid evolution of chemosensory receptor genes in a pair of sibling species of orchid bees (Apidae: Euglossini)

Background

Insects rely more on chemical signals (semiochemicals) than on any other sensory modality to find, identify, and choose mates. In most insects, pheromone production is typically regulated through biosynthetic pathways, whereas pheromone sensory detection is controlled by the olfactory system. Orchid bees are exceptional in that their semiochemicals are not produced metabolically, but instead male bees collect odoriferous compounds (perfumes) from the environment and store them in specialized hind-leg pockets to subsequently expose during courtship display. Thus, the olfactory sensory system of orchid bees simultaneously controls male perfume traits (sender components) and female preferences (receiver components). This functional linkage increases the opportunities for parallel evolution of male traits and female preferences, particularly in response to genetic changes of chemosensory detection (e.g. Odorant Receptor genes). To identify whether shifts in pheromone composition among related lineages of orchid bees are associated with divergence in chemosensory genes of the olfactory periphery, we searched for patterns of divergent selection across the antennal transcriptomes of two recently diverged sibling species Euglossa dilemma and E. viridissima.

Results

We identified 3185 orthologous genes including 94 chemosensory loci from five different gene families (Odorant Receptors, Ionotropic Receptors, Gustatory Receptors, Odorant Binding Proteins, and Chemosensory Proteins). Our results revealed that orthologs with signatures of divergent selection between E. dilemma and E. viridissima were significantly enriched for chemosensory genes. Notably, elevated signals of divergent selection were almost exclusively observed among chemosensory receptors (i.e. Odorant Receptors).

Conclusions

Our results suggest that rapid changes in the chemosensory gene family occurred among closely related species of orchid bees. These findings are consistent with the hypothesis that strong divergent selection acting on chemosensory receptor genes plays an important role in the evolution and diversification of insect pheromone systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0451-9) contains supplementary material, which is available to authorized users.

Identification and Differential Expression of a Candidate Sex Pheromone Receptor in Natural Populations of Spodoptera litura

Olfaction is primarily mediated by highly specific olfactory receptors (ORs), a subfamily of which are the pheromone receptors that play a key role in sexual communication and can contribute to reproductive isolation. Here we cloned and identified an olfactory receptor, SlituOR3 (Genbank NO. JN835270), from Spodoptera litura, to be the candidate pheromone receptor. It exhibited male-biased expression in the antennae, where they were localized at the base of sensilla trichoidea. Conserved orthologues of these receptors were found amongst known pheromone receptors within the Lepidoptera, and SlituOR3 were placed amongst a clade of candidate pheromone receptors in a phylogeny tree of insect ORs. SlituOR3 is required for the EAG responses to both Z9E11-14:OAc and Z9E12-14:OAc SlituOR3 showed differential expression in S. litura populations attracted to traps baited with a series of sex pheromone blends composed of different ratios of (9Z,11E)-tetradecadienyl acetate (Z9E11-14:OAc) and (9Z,12E)-tetradecadienyl acetate (Z9E12-14:OAc). The changes in the expression level of SlitOR3 and antennal responses after SlitOR3 silencing suggested that SlitOR3 is required for the sex pheromone signaling. We infer that variation in transcription levels of olfactory receptors may modulate sex pheromone perception in male moths and could affect both of pest control and monitoring efficiency by pheromone application after long time mass trapping with one particular ratio of blend in the field.

A recombination suppressor contributes to ecological speciation in OSTRINIA moths

Despite unparalleled access to species' genomes in our post-genomic age, we often lack adequate biological explanations for a major hallmark of the speciation process-genetic divergence. In the presence of gene flow, chromosomal rearrangements such as inversions are thought to promote divergence and facilitate speciation by suppressing recombination. Using a combination of genetic crosses, phenotyping of a trait underlying ecological isolation, and population genetic analysis of wild populations, we set out to determine whether evidence supports a role for recombination suppressors during speciation between the Z and E strains of European corn borer moth (Ostrinia nubilalis). Our results are consistent with the presence of an inversion that has contributed to accumulation of ecologically adaptive alleles and genetic differentiation across roughly 20% of the Ostrinia sex chromosome (~4 Mb). Patterns in Ostrinia suggest that chromosomal divergence may involve two separate phases-one driving its transient origin through local adaptation and one determining its stable persistence through differential introgression. As the evolutionary rate of rearrangements in lepidopteran genomes appears to be one of the fastest among eukaryotes, structural mutations may have had a disproportionate role during adaptive divergence and speciation in Ostrinia and in other moths and butterflies.

Genetic mapping of two components of reproductive isolation between two sibling species of moths, Ostrinia nubilalis and O. scapulalis

We report the quantitative trait loci (QTL) mapping of reproductive isolation traits between Ostrinia nubilalis (the European corn borer) and its sibling species O. scapulalis (the Adzuki bean borer), focusing on two traits: mating isolation (mi) and pheromone production (Pher). Four genetic maps were generated from two backcross families, with two maps (one chromosomal map and one linkage map) per backcross. We located 165-323 AFLP markers on these four maps, resulting in the identification of 27-31 linkage groups, depending on the map considered. No-choice mating experiments with the offspring of each backcross led to the detection of at least two QTLs for mi in different linkage groups. QTLs underlying Pher were located in a third linkage group. The Z heterochromosome was identified by a specific marker (Tpi) and did not carry any of these QTLs. Finally, we considered the global divergence between the two sibling species, distortions of segregation throughout the genome, and the location and effect of mi and Pher QTLs in light of the known candidate genes for reproductive isolation within the genus Ostrinia and, more broadly, in phytophagous insects.

The language of speciation

The literature on speciation has expanded dramatically in recent years, catalyzed by the emergence of new conceptual frameworks, new theoretical approaches, and new methods for characterizing pattern and inferring process. As a consequence, the language used to describe the speciation process has become more complex. Increasing complexity may be an accurate reflection of current thinking with respect to how phenotypic differences limit gene flow, how selection results in the evolution of reproductive isolation, and genetic changes that contribute to speciation. However, increased language complexity has come at a cost; old definitions have been reconfigured and new terms have been introduced. In some instances, the introduction of new terminology has failed to recognize historical usage, leading to unnecessary ambiguity and redundancy. Although the writings of Mayr and Dobzhansky remain a reference point in the language of speciation, the last decades of the 20th century saw substantial changes in our thinking about the speciation process. During that period, the language of speciation remained relatively stable. In contrast, the first decade of the 21st century has witnessed a remarkable expansion of the language of speciation. Here, the origin and evolution of ideas about speciation are viewed through the lens of changing language use.

Large-scale candidate gene scan reveals the role of chemoreceptor genes in host plant specialization and speciation in the pea aphid

Understanding the drivers of speciation is critical to interpreting patterns of biodiversity. The identification of the genetic changes underlying adaptation and reproductive isolation is necessary to link barriers to gene flow to the causal origins of divergence. Here, we present a novel approach to the genetics of speciation, which should complement the commonly used approaches of quantitative trait locus mapping and genome-wide scans for selection. We present a large-scale candidate gene approach by means of sequence capture, applied to identifying the genetic changes underlying reproductive isolation in the pea aphid, a model system for the study of ecological speciation. Targeted resequencing enabled us to scale up the candidate gene approach, specifically testing for the role of chemosensory gene families in host plant specialization. Screening for the signature of divergence under selection at 172 candidate and noncandidate loci, we revealed a handful of loci that show high levels of differentiation among host races, which almost all correspond to odorant and gustatory receptor genes. This study offers the first indication that some chemoreceptor genes, often tightly linked together in the genome, could play a key role in local adaptation and reproductive isolation in the pea aphid and potentially other phytophagous insects. Our approach opens a new route toward the functional genomics of ecological speciation.

Sex-linked transcription factor involved in a shift of sex-pheromone preference in the silkmoth Bombyx mori

In the sex-pheromone communication systems of moths, odorant receptor (Or) specificity as well as higher olfactory information processing in males should be finely tuned to the pheromone of conspecific females. Accordingly, male sex-pheromone preference should have diversified along with the diversification of female sex pheromones; however, the genetic mechanisms that facilitated the diversification of male preference are not well understood. Here, we explored the mechanisms involved in a drastic shift in sex-pheromone preference in the silkmoth Bombyx mori using spli mutants in which the genomic structure of the gene Bmacj6, which encodes a class IV POU domain transcription factor, is disrupted or its expression is repressed. B. mori females secrete an ∼11:1 mixture of bombykol and bombykal. Bombykol alone elicits full male courtship behavior, whereas bombykal alone shows no apparent activity. In the spli mutants, the behavioral responsiveness of males to bombykol was markedly reduced, whereas bombykal alone evoked full courtship behavior. The reduced response of spli males to bombykol was explained by the paucity of bombykol receptors on the male antennae. It was also found that, in the spli males, neurons projecting into the toroid, a compartment in the brain where bombykol receptor neurons normally project, responded strongly to bombykal. The present study highlights a POU domain transcription factor, Bmacj6, which may have caused a shift of sex-pheromone preference in B. mori through Or gene choice and/or axon targeting.