# Transcriptomic response of female adult moths to host and non-host plants in two closely related species

Marion Orsucci, Philippe Audiot, Sabine Nidelet, Franck Dorkeld, Alexandra Pommier, Marin Vabre, Dany Severac, Marine Rohmer, Bernhard Gschloessl, Rejane Streiff

► To cite this version:

Marion Orsucci, Philippe Audiot, Sabine Nidelet, Franck Dorkeld, Alexandra Pommier, et al. Transcriptomic response of female adult moths to host and non-host plants in two closely related species. BMC Evolutionary Biology, 2018, 18, 10.1186/s12862-018-1257-3. hal-01918051

HAL Id: hal-01918051

https://hal.science/hal-01918051

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

Open Access

<figure><img src="image_2.png" /></figure>

CrossMark

# Transcriptomic response of female adult moths to host and non-host plants in two closely related species

M. Orsucci¹,²,⁴* , P. Audiot¹, S. Nidelet¹, F. Dorkeld¹, A. Pommier¹, M. Vabre³, D. Severac⁵, M. Rohmer⁵, B. Gschloessl¹ and R. Streiff¹,²*

**Abstract**

**Background:** Divergent selection has been shown to promote speciation in many taxa and especially in phytophagous insects. In the *Ostrinia* species complex, the European corn borer (ECB) and adzuki bean borer (ABB) are two sibling species specialized to different host plants. The first is a well-known maize pest, whereas the second is a polyphagous species associated with various dicotyledons. Their specialization to host plants is driven by morphological, behavioral and physiological adaptations. In particular, previous studies have shown that ECB and ABB display marked behavior with regard to plant choice during oviposition, involving specific preference and avoidance mechanisms. In this study, our goal was to identify the mechanisms underlying this host-plant specialization in adult females through an analysis of their gene expression. We assembled and annotated a de novo reference transcriptome and measured differences in gene expression between ECB and ABB females, and between environments. We related differentially expressed genes to host preference behavior, and highlighted the functional categories involved. We also conducted a specific analysis of chemosensory genes, which are considered to be good candidates for host recognition before oviposition.

**Results:** We recorded more differentially expressed genes in ECB than in ABB samples, and noticed that the majority of genes potentially involved in the host preference were different between the two species. At the functional level, the response to plant environment in adult females involved many processes, including the metabolism of carbohydrates, lipids, proteins, and amino acids; detoxification mechanisms and immunity; and the chemosensory repertoire (as expected). Until now, most of the olfactory receptors described in *Ostrinia* spp. had been tested for their putative role in pheromone recognition by males. Here we observed that one specific olfactory receptor was clearly associated with ECB's discrimination between maize and mugwort conditions, highlighting a potential new candidate involved in plant odor discrimination in adult females.

**Conclusions:** Our results are a first step toward the identification of candidate genes and functions involved in chemosensory processes, carbohydrate metabolism, and virus and retrovirus dynamics. These candidates provide new avenues for research into understanding the role of divergent selection between different environments in species diversification.

**Keywords:** Specialization, Oviposition, Insect-plant interaction, Adaptation, Transcriptomics

<footnote>* Correspondence: marion.orsucci@gmail.com; rejane.streiff@inra.fr
¹CBGP, INRA, CIRAD, IRD, Montpellier SupAgro, Univ Montpellier, Montpellier, France
Full list of author information is available at the end of the article</footnote>
## Background

Speciation processes, by which populations diverge and become reproductively isolated, may be driven by divergent selection between different environments that affect various life history traits directly or indirectly related to fitness [12, 20, 22, 57, 72, 73]. In phytophagous insects, host races are relevant examples of population divergence driven or reinforced by host-plant specialization which may, in some cases, be a key component of reproductive isolation [14, 17, 21, 51]. Host plants directly affect the fitness of herbivores [2] due to their role in nutrition [14, 75], protection [26, 36], and reproduction [2, 83]. Indeed, plant metabolites are known to affect the rate of partner encounter and mating success [15, 62], the level of egg laying [18, 41], the choice of oviposition site of gravid females [67], as well as the development of eggs and larvae [56]. In addition to intrinsic host plant qualities, higher trophic level interactions such as plant-associated microorganisms and parasitoids may also impact realized fecundity by impairing larval performance [3]. The herbivorous host range results from these complex pluritrophic interactions, which occur throughout the entire life cycle of the insects [56].

In holometabolous species such as Lepidoptera, oviposition choice is crucial because, before the imago stage, individuals often have limited mobility and thus depend on a judicious choice of plant by the adult female [67]. The preference-performance hypothesis (or “mother-knows-best” hypothesis) suggests that the host preference of adult females should maximize the fitness of their offspring by ovipositing on the most suitable hosts [30]. In their literature survey, Gripenberg et al. [25] suggested that empirical data did support this hypothesis for a majority of the study cases. Molecular and functional studies highlighted the importance of chemosensory mechanisms for the recognition by ovipositing females of olfactory cues emitted by the plant during egg laying [10, 46]. Interestingly, Trona et al. [80] stressed that, in nature, sex pheromones and plant odors are perceived and coded as a whole, and suggested a dual role of plant signals in habitat selection and premating sexual communication.

European corn borer (ECB, *Ostrinia nubilalis* Hübner) and adzuki bean borer (ABB, *Ostrinia scapulalis* Walker), as defined by Frolov et al. [19], are two sibling species of phytophagous moths, co-occurring in a large part of the European continent. ABB occurs in the Palearctic region from Japan to Western Europe. ECB, which is thought to originate from Western Europe, occurs in Northern Africa, Western Europe, and Western Asia and was introduced in the Nearctic region (Northern America) at the end of the nineteenth century. To date, there is no simple phenotypic or molecular characteristics to distinguish these two species. A recent revision of the genus [19], which compiled dozens of years of genetic and phenotypic

data (morphology, pheromones, calling time and molecular markers), has concluded that the best criterion for discriminating between ECB and ABB regarding their reproductive isolation was the type of host plant. Indeed, these two species are specialized to different host ranges: ECB is mainly observed on an agricultural monoculture don, the maize *Zea mays* L., and ABB on several dicotyledons like mugwort (*Artemisia vulgaris* L.), hop (*Humulus lupulus* L.) and hemp (*Cannabis sativa* L.; [19]). In ECB and ABB, the preference-performance hypothesis has typically been observed through various sampling and independent studies [8, 49, 59], but both species exhibit contrasting behaviors and strategies. In a previous experiment, ECB and ABB adult moths were exposed to three plant environments during their mating and egg laying phases: a pure maize environment, a pure mugwort environment and a mixed environment with maize and mugwort plants [59]. The authors observed a strong preference for maize in ECB gravid females when they had the choice between maize and mugwort, coupled with a significant avoidance of mugwort when this plant was the only possible choice (see Additional file 1: Box S1 summary for details). In contrast, ABB showed a moderate preference for mugwort over maize during oviposition and no avoidance behavior. In addition to its oviposition behavior, higher survival and egg-laying rates were observed in ECB when maize was in its environment, with a preference for maize plants also as a resting site. In ABB, these latter traits were not affected by the plant environment (see Additional file 1: Box S1–2). These results suggested attractive and/or repellent molecules impacting oviposition choice [42] as well as synergetic effects on mating and egg laying success [34]. The olfactory responses of ECB females to plant odors have been previously characterized when they were presented with: i) injured and uninjured maize plants [74], ii) blends of maize volatiles in specific ratios [55, 78], and iii) maize vs. mugwort or hop odors [42]. However, the molecular mechanisms (e.g. genes) behind plant recognition and choice are still unknown.

We present here a comprehensive gene expression analysis of ECB and ABB adult females by means of high throughput RNA sequencing of females placed in various plant environments (sampled from the experiment in [59], see Additional file 1: Box S1–1 for details). We assembled and annotated a de novo reference transcriptome and measured differences in gene expression, first between ECB and ABB, then within each moth species between females confronted to contrasting plant environments: pure maize, pure mugwort, and a mixture of maize and mugwort. Furthermore, we linked differentially expressed (DE) genes to the host preference in ECB and ABB and investigated the functional categories involved. Alongside this exhaustive de novo transcriptomic survey, and because of the
importance of chemosensory genes, which are considered as good candidates for host recognition before oviposition, we conducted a specific analysis of a chemosensory gene subset previously characterized in ECB, ABB and other related *Ostrinia* species. As for the de novo transcripts, we measured gene expression differences in these chemosensory genes between ECB and ABB, and within each moth species, between contrasting plant environments. The discovery of candidate molecular mechanisms, combined with our knowledge of specialization patterns, offers an opportunity to better understand the role of the host plants in the divergence between ABB and ECB.

## Results

### ECB- and ABB-ref assemblies and annotation

We sequenced RNA from adult ECB and ABB moths sampled from a previous experiment [59], where phenotypic and behavioral traits in moths faced with different plant environments were measured. Three plant environments were tested: i) “maize” conditions where only maize plants were available, they being the usual host for ECB; ii) “mugwort” conditions where only mugwort plants were available, they being the usual host for ABB; and iii) “choice” conditions where both plants were available (the entire experimental framework is given in Additional file 1: Box S1–1.) Twenty-four RNA samples corresponding to two repetitions, three experimental conditions, two tissues (head-thorax: HT, and abdomen: ABDO) and two moth species (ABB and ECB) were sequenced with HiSeq technology, and raw reads were assembled to build de novo transcripts. Details on raw reads and contaminants per library are given in Additional file 2: Table S1, and the main features of the de novo assemblies (HiSeq only or HiSeq-454 combined) are detailed in Table 1. We observed that combining published 454 transcripts and de novo HiSeq transcripts doubled the length of transcripts (Table 1) compared to de novo HiSeq assembly alone. Taking the longest sequence as the only representative for each unigene resulted in 9415 ECB genes and 9992 ABB

genes (Table 1). These sets of genes are hereafter referred as ECB-ref and ABB-ref respectively. Homology searches in genomic and protein databases yielded significant hits for 61% of the ECB-ref genes, of which 80% had an associated GO (Gene Ontology) term. More than 97% of the annotated genes were homologous to Lepidoptera genes, followed by homologs to other insect orders (Additional file 2: Table S2). To a much lesser extent, hits with genes of plants, viruses and bacteria were observed. Among these microorganisms, viruses were most prevalent. Details about their detection and annotations are given in Additional file 2: Table S5-A.

### DE genes in ECB- and ABB-ref: Global pattern

DE gene detection was done in both ECB-ref and ABB-ref by mapping sequencing reads for each experimental condition on de novo transcripts. Qualitative trends were observed via multivariate analyses, and quantitative statistics (e.g. log₂ fold change, Log2FC hereafter) and models were designed with the “EdgeR” function [70]. Generalized linear models (GLMs) were in particular designed to test for the effects of the moth species (ABB vs. ECB) and experimental conditions (maize, mugwort and choice) on gene expression.

Quantitative and qualitative trends were highly correlated between ECB-ref and ABB-ref analyses (data not shown). Thus for simplicity we will present the results for ECB-ref only. Rather than constructing a mixed reference by pooling ECB and ABB reads, we chose to conduct two analyses to avoid chimeric transcripts between species. On a global analytic scale, we did not observe noticeable differences in the number or nature (function and GO term) of the DE genes. Nonetheless, we specifically searched for ABB-ref genes strongly mapped by ABB reads but not by ECB reads. These genes, which probably did not appear in the ECB-ref, may account for some ABB peculiarities like less stringent host choice. Genes with such patterns were very rare and, to date, they have not hit with a known protein or function (Additional file 2: Table S3).

<table><caption>**Table 1** Sequencing data and principal features of ECB and ABB assemblies</caption><thead><tr><th></th><th></th><th></th><th>ECB</th><th>ABB</th></tr></thead><tbody><tr><td rowspan="2">Number of raw reads</td><td rowspan="2">HiSeq</td><td>after sequencing</td><td>632,618,343</td><td>655,684,618</td></tr><tr><td>after filtering</td><td>466,145,250</td><td>498,030,304</td></tr><tr><td>Number of unigenes</td><td>HiSeq</td><td></td><td>41,155</td><td>42,913</td></tr><tr><td rowspan="3">Number of transcripts</td><td>HiSeq</td><td></td><td>41,212</td><td>42,969</td></tr><tr><td>HiSeq-454</td><td>all transcripts</td><td>11,456</td><td>12,024</td></tr><tr><td>HiSeq-454</td><td>longest transcript per unigene</td><td>9415</td><td>9992</td></tr><tr><td rowspan="2">Transcript N50 length (bp)</td><td>HiSeq</td><td></td><td>530.41</td><td>544.16</td></tr><tr><td>HiSeq-454</td><td></td><td>1221</td><td>1260</td></tr><tr><td>Assembled transcriptome (Mbp)</td><td>HiSeq-454</td><td></td><td>13.7</td><td>15.2</td></tr></tbody></table>
The overall percentage of mapped reads on the
ECB-ref transcripts varied from 64 to 71% in the 12 HT
RNA libraries, and from 58 to 65% in the 12 ABDO li-
braries (Additional file 2: Table S4). The proportion of
multiple mapped reads was less than 1% in HT (0.51%)
and ABDO (0.39%) libraries. After filtering genes with
coverage of 1X or less, we identified a total of 9326 and
9186 genes expressed across all experimental conditions
in the HT and ABDO samples, respectively. Most of
them (around 98% in HT and 95% in ABDO) were
shared across all experimental conditions (Additional file
3: Figure S1). We detected 24 and 17 genes found to be
expressed only in ECB HT and ABDO samples, respect-
ively, and 19 and 118 genes found to be expressed only
in ABB HT and ABDO samples, respectively. Fisher's
exact test did not reveal any significant enrichment in
GO categories for these lists of genes specific to the ECB
and ABB libraries.

DE genes and source of variation: Species and
experimental effects

Apart from the tissues (Fig. 1a), the main source of differ-
ential expression was due to contrasts between ECB and
ABB samples (Fig. 1b). In HT samples, we detected 337 and
268 transcripts significantly over- and under-expressed in
ECB vs. ABB. In ABDO samples, 828 and 132 transcripts
were significantly over- and under-expressed in ECB vs.
ABB (Table 2, model 2).

Annotated DE genes in HT tissues between ECB and ABB (Additional file 2: Table S6) were mostly enriched in GO terms (Additional file 2: Table S7) relating to: i) DNA/RNA metabolic processes (e.g. reverse transcriptases, transposases, pol proteins and transferases) and ii) organic compound metabolic processes (e.g. fructofuranosidases and UDP-glucuronosyltransferases). In most DE genes (|Log2FC| > 2), we also observed several cytochrome P450s involved in detoxification processes. One had previously been characterized in *O. furnacalis* (gb|ACF17813.2, Additional file 2: Table S6) to the DIMBOA toxic compound of maize, and another was an immune-induced protein previously characterized in *Ostrinia* species (gb|AGV28583.1, Additional file 2: Table S6). In ABDO tissues, we observed enrichment in GO terms relating to DNA metabolic processes and molecular binding and transport (Additional file 2: Table S7). In most DE genes (|Log2FC| > 5, Additional file 2: Table S6), we identified apolipoporphins potentially involved in immune response [87] as well as the DIMBOA-induced cytochrome P450 previously described in *O. furnacalis* (gb|ACF17813.2, |log2FC| > 5, Additional file 2: Table S6; [88]).

The experimental condition had the least influence on
the variation in gene expression (Table 2, model 2, Fig. 1c).
Between 122 and 187 DE genes were observed through
comparisons between the experimental conditions (choice

vs. mugwort, choice vs. maize, and mugwort vs. maize) in
HT samples, and even fewer were observed in ABDO sam-
ples (between 9 and 35). These contrasting numbers of DE
genes reflected a higher between-repetition dispersion and
between-experimental condition overlap in ECB ABDO
tissues compared to ECB HT tissues (Fig. 2a-b), as well as
in ABB samples vs. ECB samples (Fig. 2c-d vs. Fig. 2a-b).
What is more, DE gene lists between the experimental
conditions (mugwort, maize and choice) mostly differed
from each other (Fig. 3).

Annotated DE genes between experimental conditions (Additional file 2: Table S8) were enriched in a wide variety of GO categories for biological processes (Additional file 2: Table S7) involving feeding and mating behavior, sensing, molecule transport and binding, organic compounds and DNA/RNA metabolism. In most DE genes (|Log2FC| > 5, Additional file 2: Table S8) of the HT and ABDO samples, we identified various enzymes potentially involved in digestion, detoxification and immunity (e.g. glucose dehydrogenase, valine-rich midgut proteins [VMP], heat shock proteins, and cytochrome P450) as well as one protein potentially involved in sensing (transient receptor potential protein).

DE genes and behavioral traits

We then focused on DE genes potentially related to specific
traits and behaviors observed at the phenotypic level
(Additional file 1: Box S1–2), using a specific contrast ana-
lysis in the GLM procedure of the EdgeR function (Table 3).
We categorized DE genes identified by model 3 (Table 3)
for maize preference in ECB (ECB-pref), mugwort avoid-
ance in ECB (ECB-avoid) and mugwort preference in ABB
(ABB-pref). Remarkably, many more DE genes were ob-
served in ECB (260 and 786 for ECB-pref, and 132 and 149
for ECB-avoid; Tables 3 and 4, and more details in
Additional file 2: Table S9) than in ABB (56 and 4 for
ABB-pref, Tables 3 and 4, and more details in Additional
file 2: Table S9). Moreover, the preference and avoidance
gene lists mostly differed from each other (Fig. 4 and
Additional file 2: Table S10).

We observed enrichment in various GO categories for the DE ECB-pref, ABB-pref and ECB-avoid genes involving the metabolism of carbohydrates, lipids, proteins, and amino acids, the detoxification mechanisms, immunity, and the chemosensory repertoire (Table 4). These functions were also observed in the gene lists specific to ECB-pref, ABB-pref, and ECB-avoid contrasts (Additional file 2: Table S10).

Lastly, we considered the DE gene lists in which the ex-
pression also differed between ECB and ABB. A reverse tran-
scriptase, various enzymes involved in digestion and
detoxification (e.g. mucin, Dzip1, Beta-fructofuranosidase,
phospholipase, and serine protease), an immune protein
(lysozyme precursor), a heat shock protein, and an egg yolk
protein (Table 4) were included. All these DE genes were in-
volved in ECB-pref and ECB-avoid, making them the best
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 1** PCA individual factor map of mean raw counts of 9413 genes in 24 RNA libraries. **a** plot on 1–2 axes, ellipses group the libraries’ mean raw counts per tissue (ABDO vs. HT); **b** plot on 2–3 axes, ellipses group the libraries’ mean raw counts per moth species (ECB vs. ABB); **c** plot on 4–5 axes, ellipses group the libraries’ mean raw counts per experimental condition (choice vs. maize vs. mugwort). Species are represented in each plot: ABB (circle) and ECB (diamond). Different colors were used to distinguish between the samples according to the experimental conditions: choice (blue), maize (red) and mugwort (yellow). The percentage of explained variation is indicated on each dimension axis</figcaption></figure>

candidates for the divergence driven by specialization to
maize.

DE transcripts in candidate sensory genes

We compiled a set of genes known to be involved in sen-
sory perception which had been previously described in
ECB, ABB and related species of the *Ostrinia* genus
(Additional file 2: Table S11). This gene set, hereafter called
“Senso-ref,” was used as a gene reference independent of
ABB-ref and ECB-ref in order to focus on sensory genes.
Some 109 chemosensory genes of the 139 Senso-ref
candidates were expressed in at least one library (Additional
file 2: Table S12). As for the de novo transcriptome, most
of the variation in expression for these genes was due to
tissue type (ABDO vs. HT) with 65 DE genes, followed by
moth species (ECB vs. ABB) with 36 DE genes. Overall, 28
genes were differentially expressed depending on the ex-
perimental conditions, of which a majority (24) differed be-
tween pure mugwort and pure maize conditions. When we
focused on DE candidate genes potentially involved in traits
and behaviors (ABB-pref, ECB-pref and ECB-avoid) ob-
served at the phenotypic level (Additional file 1: Box S1–2,
Additional file 2: Table S12) using specific contrast analysis
in the GLM procedure, we highlighted eight genes involved
in maize preference in ECB samples: four odorant receptors
(ORs), one aldehyde oxidase, one carboxylesterase, and one
ionotropic receptor (IR). Only two genes (ORs) were found
in mugwort avoidance in ECB samples, and three genes
(ORs) in mugwort preference in ABB samples. Notably,
genes differentially expressed in ECB-pref differed from
those detected in ABB-pref and ECB-avoid (Additional
file 2: Table S9). Finally, of the 28 genes differentially
expressed depending on the plant environment, 11 (5 ORs,
5 IRs and 1 carboxylesterase) were also differentially
expressed between ABB and ECB samples.

Discussion

Using an RNA-Seq approach without a priori and a specific
analysis targeting chemosensory genes, we searched for
genes expressed in ECB and ABB adult females exposed to
various environments with or without their usual host
plant. We identified genes and functions that differentiated
ECB from ABB whatever, or depending on, the
<table><caption>**Table 2** GLM model 1</caption><thead><tr><th rowspan="2">Models</th><th rowspan="2"># libraries</th><th rowspan="2">Comparison between</th><th rowspan="2">Contrast used<br>(1; -1)</th><th colspan="2"># DE transcripts</th></tr><tr><th>HT tissue</th><th>ABDO tissue</th></tr></thead><tbody><tr><td rowspan="12">Model 1: ti = modalities<sup>a</sup></td><td rowspan="12">12</td><td rowspan="3">moth species</td><td rowspan="3">ECB vs. ABB</td><td>total = 605</td><td>total = 960</td></tr><tr><td>under-expressed = 268</td><td>under-expressed = 132</td></tr><tr><td>over-expressed = 337</td><td>over-expressed = 828</td></tr><tr><td rowspan="3">experimental condition</td><td rowspan="3">Choice vs. Maize</td><td>total = 122</td><td>total = 9</td></tr><tr><td>under-expressed = 66</td><td>under-expressed = 4</td></tr><tr><td>over-expressed = 56</td><td>over-expressed = 5</td></tr><tr><td rowspan="3">Choice vs. Mugwort</td><td>total = 187</td><td>total = 25</td></tr><tr><td>under-expressed = 118</td><td>under-expressed = 1</td></tr><tr><td>over-expressed = 69</td><td>over-expressed = 24</td></tr><tr><td rowspan="3">Maize vs. Mugwort</td><td>total = 129</td><td>total = 35</td></tr><tr><td>under-expressed = 86</td><td>under-expressed = 1</td></tr><tr><td>over-expressed = 43</td><td>over-expressed = 34</td></tr></tbody></table>

<figure><img src="image_6.png" /></figure>

<footnote><sup>a</sup>"modalities" is a fusion of two variables: moth species with experimental condition</footnote>
<figure><img src="image_7.png" /><figcaption>**Fig. 3** Venn diagram of DE genes between experimental conditions (mugwort vs. choice, maize vs. choice, and mugwort vs. maize) in HT (a) and ABDO (b) tissues</figcaption></figure>

environment. We also classified these DE genes according to their expression pattern to relate them to preference and avoidance behaviors in ECB and ABB gravid females. In this way we identified candidate genes for species divergence by host specialization.

### Response to plant environment in ECB and ABB adult females: Global pattern

On average across all transcripts, we observed first of all that more genes were up- or down-regulated in HT tissues than in ABDO tissues. This observation corroborated the expectation that sensory and neuronal signals, which are mainly active in HT organs, would be major factors in oviposition [66, 80]. Second, we recorded more DE genes in ECB samples than in ABB samples. We may assume that this result is a signature of the strength and variability of host specialization in ECB (combining host preference and alternative host avoidance), while ABB appeared less specialized to its own host in terms of physiological and behavioral traits ([59], Additional file 1: Box S1). Lastly, we found that the

genes involved in the preference for their hosts differed between ECB and ABB, as well as those involved in mugwort rejection in ECB, suggesting that the ECB shift onto maize was characterized by genetic novelties.

Response to plant environment in ECB and ABB adult females involved many processes, including metabolism of carbohydrates, lipids, proteins, and amino acids, detoxification mechanisms, immunity and chemosensory repertoire. These functions are key processes of host specialization in phytophagous insects [77], but were until now mainly illustrated at the larval stage [16, 29, 32, 43, 63, 76]. In the present study, we provide novel results for the adult stage, revealing some other candidates for host specialization like for example the Takeout-like proteins, previously reported as involved in circadian rhythms and feeding behavior in *Drosophila* [71].

### Chemosensory genes in host choice: New candidate genes and divergence between ECB and ABB

Targeted analyses of chemosensory genes supported their implication in host discrimination in ECB, ABB, or

<table><caption>**Table 3** GLM model 2</caption><thead><tr><th rowspan="2">Targeted trait</th><th rowspan="2"># libraries</th><th rowspan="2">Models</th><th rowspan="2">Contrast used<br>(0.5; 0.5; -1)</th><th colspan="2"># DE transcripts</th></tr><tr><th>HT tissue</th><th>ABDO tissue</th></tr></thead><tbody><tr><td rowspan="3">Maize preference in ECB</td><td rowspan="3">6</td><td rowspan="3">Model 3:<br>ti = experimental condition</td><td>Choice and Maize vs. Mugwort</td><td>total = 260</td><td>total = 786</td></tr><tr><td>under-expressed = 148</td><td>under-expressed = 608</td></tr><tr><td>over-expressed = 118</td><td>over-expressed = 178</td></tr><tr><td rowspan="3">Mugwort avoidance in ECB</td><td rowspan="3"></td><td rowspan="3">Choice and Mugwort vs. Maize</td><td>Choice and Maize vs. Mugwort</td><td>total = 132</td><td>total = 149</td></tr><tr><td>under-expressed = 84</td><td>under-expressed = 128</td></tr><tr><td>over-expressed = 48</td><td>over-expressed = 21</td></tr><tr><td rowspan="3">Mugwort preference in ABB</td><td rowspan="3">6</td><td rowspan="3">Choice and Mugwort vs. Maize</td><td>Choice and Maize vs. Mugwort</td><td>total = 56</td><td>total = 4</td></tr><tr><td>under-expressed = 15</td><td>under-expressed = 3</td></tr><tr><td>over-expressed = 41</td><td>over-expressed = 1</td></tr></tbody></table>
<table><caption>**Table 4** DE genes involved in ECB-pref, ECB-avoid and ABB-pref</caption><thead><tr><th rowspan="2">Transcript ID</th><th colspan="2">ECB</th><th colspan="2">ECB avoid</th><th colspan="2">ABB</th><th colspan="2">Species</th><th colspan="2">nr BLAST hit</th><th rowspan="2">annotation</th></tr><tr><th>pref</th><th>ABDO</th><th>pref</th><th>HT</th><th>pref</th><th>ABDO</th><th>pref</th><th>HT</th><th>pref</th><th>ABDO</th></tr></thead><tbody><tr><td>isotig08066</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>x</td><td></td><td>mucin-17-like [Amyelois transitella]</td></tr><tr><td>isotig04815</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>juvenile hormone esterase-like [Amyelois transitella]</td></tr><tr><td>isotig06592</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>glucose dehydrogenase [acceptor], partial [Zootermopsis nevadensis]</td></tr><tr><td>isotig07523</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>protein takeout-like [Bombyx mori]</td></tr><tr><td>isotig01371</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>x</td><td></td><td></td><td>protein split ends-like [Amyelois transitella]</td></tr><tr><td>isotig08378</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>coagulation factor-like protein 3 [Hypantria cunea]</td></tr><tr><td>isotig09172</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VMP25 [Manduca sexta]</td></tr><tr><td>isotig04519</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>x</td><td></td><td></td><td>lysozyme precursor [Ostrinia furnacalis]</td></tr><tr><td>isotig09335</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>x</td><td></td><td>serine protease K12H4.7 [Amyelois transitella]</td></tr><tr><td>isotig03611</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>VMP93 [Manduca sexta]</td></tr><tr><td>isotig04272</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>lipase 3-like [Amyelois transitella]</td></tr><tr><td>isotig07502</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>follicular epithelium yolk protein subunit [Plodia interpunctella]</td></tr><tr><td>isotig09064</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>zinc finger protein Dzip1 [Papilio polytes]</td></tr><tr><td>isotig05796</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>x</td><td></td><td></td><td>zinc finger protein Dzip1-like [Plutella xylostella]</td></tr><tr><td>isotig08253</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>x</td><td></td><td></td><td>follicle-specific yolk polypeptide-4 [Galleria mellonella]</td></tr><tr><td>isotig04271</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>x</td><td></td><td></td><td>heat shock protein 70-1 [Melitaea cinxia]</td></tr><tr><td>isotig08527</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>x</td><td></td><td></td><td>Jerky protein [Operophtera brumata]</td></tr><tr><td>isotig01095</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>homeobox protein SMOX-1-like [Papilio polytes]</td></tr><tr><td>isotig08752</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>transient receptor potential protein-like [Amyelois transitella]</td></tr><tr><td>isotig03905</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>peroxidase-like [Papilio polytes]</td></tr><tr><td>isotig10705</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>transient receptor potential protein-like [Amyelois transitella]</td></tr><tr><td>isotig08302</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>transient receptor potential protein-like [Amyelois transitella]</td></tr><tr><td>isotig10768</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>x</td><td></td><td></td><td>RNA-directed DNA polymerase [Operophtera brumata]</td></tr><tr><td>isotig01012</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>x</td><td></td><td></td><td>Beta-fructofuranosidase 2 [Operophtera brumata]</td></tr><tr><td>isotig06078</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>plexin domain-containing protein 2-like [Plutella xylostella]</td></tr><tr><td>isotig09499</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>x</td><td></td><td></td><td>zinc finger CCHC domain-containing protein 24-like [Bombyx mori]</td></tr><tr><td>isotig06672</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>x</td><td></td><td>calcium-independent phospholipase A2-gamma-like [Papilio xuthus]</td></tr><tr><td>isotig09126</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>elongation of very long chain fatty acids protein [Papilio xuthus]</td></tr><tr><td>isotig09158</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>RNA-directed DNA polymerase from transposon BS [Papilio machaon]</td></tr></tbody></table>

For clarity, transcripts with no hits or with anonymous annotation (e.g. 'uncharacterized protein') were removed. In red: up-regulated; in blue: down-regulated. Crosses are for genes differentially regulated between species. A more detailed table is given in the Additional file 2: (Table S9) Genes with a |Log2FC| > 2 are presented with their best "nr" hit annotation

both. Moreover, five ORs, five IRs and one carboxyles-
terase appeared as suitable candidates for their putative
role in species divergence by specialization. Indeed, their
expression responded to plant environments and, at the
same time, differed between moth species. Among these
genes, one OR (LC002733.1, Additional file 2: Table S12)

was particularly relevant, because it clearly discriminated
between maize and mugwort conditions in ECB, either for
maize preference or mugwort avoidance. In contrast to
many other ORs involved in pheromone recognition by
males, Wanner et al. [82] and Yang et al. [85] observed
that this OR was expressed higher in ECB females. We

a HT tissue

<figure><img src="image_8.png" /></figure>

b ABDO tissue

<figure><img src="image_9.png" /><figcaption>**Fig. 4** Venn diagram of DE genes involved in ECB-pref, ECB-avoid and ABB-pref in HT (a) and ABDO (b) tissues. The ECB samples are represented by warm colors and the ABB samples by cool colors</figcaption></figure>
suggest that the DE analysis may lead to the discovery of new potential candidates involved in plant odor discrimination in females during oviposition. Lastly, among the genes differentially expressed between moth species and potentially involved in ECB-pref, we identified one OR (JN169130, [40]) previously reported as more expressed in males than in females in ECB and the Asian corn borer (ACB), another sibling species of the *Ostrinia* species complex. This receptor may thus be an example of pleiotropy or synergetic activity for pheromone recognition in males and plant recognition in females. An essential requirement for speciation via specialization is the existence of a link between reproductive isolation and specialization due to divergent environment selection [51]. In ECB and ABB, two pheromone lineages coexist and are partially reproductively isolated due to the discrimination by males of the pheromone blends emitted by females [23, 37]. A gene implicated in both pheromone recognition by males and host plant discrimination for oviposition by females would be a direct link between mating isolation and host specificity [81].

## Major role of metabolism of organic compounds and divergence between ECB and ABB

Besides the recognition of plant chemical cues, the analysis based on whole de novo assembly highlighted a wide variety of other functional categories involved in either moth species divergence, host plant choice, or both.

The organic metabolic processes represented by key genes previously identified and involved in the differentiation between ECB and ABB [60] as well as in various phytophagous lepidopterans [77] highlighted the importance of the digestion and detoxification of plant metabolites, in particular at the larval stage. At the adult stage, chemosensory processes have been emphasized, because ovipositing females use plant secondary metabolites as cues to locate and select suitable plants [79]. Our study suggests that the digestion and detoxification of non-volatile compounds are also important factors in host choice. This result may be unexpected because ECB is often considered as a “capital breeder,” i.e. as an insect that feeds only during the larval stage, as is often the case for lepidopterans [7]. On the contrary, our results and previous work [13, 39] show that nutrition (and in particular carbohydrates) may play a key role in oviposition behavior, egg weight, fecundity and female longevity.

Interestingly, the activation of genes belonging to specific metabolic pathways differed between ECB and ABB, with for example mucin, serine protease and phospholipase genes involved in ECB-pref and ABB-pref. Hence these genes require further investigation to test whether they reflect differences in the metabolic strategies employed by ECB and ABB to survive, develop, and reproduce on their respective hosts.

## Unexpected involvement: Virus and retrovirus-like transposons

Besides the metabolism of organic compounds, in most DE genes we also observed functions relating to detoxification (cytochrome P450) and an immune-induced protein previously characterized in *Ostrinia* species (gb|AGV28583.1, Additional file 2: Table S6). These immune and detoxification genes may be mobilized to detoxify plant metabolites as well as being involved in the defense against bacteria and viruses. Indeed, we observed the genes of plants, viruses and bacteria, which were probably due to remnants of the food bolus or the presence of environmental, commensal, symbiotic, or pathogenic microorganisms. The RNA sequencing was not designed to target the microorganism community, so a fine and quantitative analysis would be unwarranted. Yet the presence of some of this alien fauna impacted the gene expression of ECB and ABB females, as illustrated by the activity of DNA-directed RNA polymerases (Additional file 2: Table S6) implicated in RNA virus replication. In particular for viruses, we observed contrasting patterns between ECB and ABB and between experimental conditions, with a higher prevalence of cypovirus-like genes in ABB HT samples in maize conditions than in ECB or other conditions (Additional file 2: Table S5-A). Thus we suggest that immune and defense genes differentially expressed between ECB and ABB, or between maize and mugwort conditions, may have been activated to contend pathogens present in the experimental environment. Interestingly, the prevalence of these microorganisms as well as the defense response of female moths differed between ECB and ABB samples, and between environments. These results, and previous similar observations at the larval stage [60], warrant further studies in other moth populations and plant species, to i) evaluate the variability and structuration of the microorganism communities in situ, ii) measure the immune response and its divergence between ECB and ABB, and iii) characterize the interactions between host plants, microorganisms and species divergence. Such multitrophic interactions have been demonstrated in different studies [64]. They merit further investigation to gauge their importance in the context of ECB-ABB divergence, since we cannot exclude opportunistic infections in experimental populations through chance or stress in a single treatment, rather than a generalizable biological reality.

Additionally, we observed the expression of retroviral-like transposons (Additional file 2: Table S5-B) and various genes associated with their transposition and spread, in particular in HT tissues of ECB. A putative role for transposable elements in ECB and ABB divergence has not yet been reported, but these elements are known to interfere with the adaptive response to various environmental stresses in
various other taxa [11]. Bursts of transposable elements may thus in some circumstances be a driver of speciation events [5]. In ECB, Lillehoj [45] proposed a scenario in which aflatoxin (a toxin produced by the microbial inhabitants in the insect's digestive tract) induced the activation of mobile genetic elements such as transposons, plasmids and chromosomal mutation, in response to environmental changes. According to this author, aflatoxin-producing fungi may be symbiotic in certain equilibrium conditions and pathogenic in conditions of nutritional imbalance, notably when elevated levels of digestible carbon compounds are present. The ecological disequilibrium caused by crop monoculture is supposed to have favored ECB populations with two or more life cycles per crop rotation over univoltine populations prone to a diapause period. Moreover, aflatoxin properties in host DNA and transposon activation may have facilitated this novelty in ECB. No similar mechanism has been described in ABB, and further studies are required to directly test the role of transposon activity in ECB-ABB divergence and maize colonization by ECB.

## Conclusions

Our results highlighted some potential candidate genes and functions involved in specialization and in species divergence between ECB and ABB. Through ad-hoc categorization, we discriminated between genes in ECB and ABB that responded to the environment with similar or divergent patterns. Independent approaches such as functional and population genomics will be required to confirm their role in the history of ECB host shift onto maize. Besides, we confirmed the role of sensing processes and identified new candidates for the chemosensory repertoire. More precise studies must now be conducted to characterize the role of these candidates in specific traits of the female adult stage such as oviposition, mating, or metabolism, since gene expression is potentially highly sensitive to many environmental factors. These candidates should also be tested in different populations and contexts, since they could help with the design of effective control solutions for crop pest management [65]. Finally, we highlighted other avenues for research on carbohydrate metabolism and virus and retrovirus dynamics, which hold the promise of enhancing our understanding of divergence between cryptic sibling species.

## Methods

### Moth sampling and rearing

In this study, we sequenced RNA from adult moths sampled from a previous experiment [59], where phenotypic and behavioral traits in moths faced with different plant environments were measured. The moths in the experiment originated from a single generation of laboratory-reared larvae initially sampled in maize and mugwort stands near Versailles (northern France). In this region, several genetic

studies had verified that the host plant (maize vs. mugwort) was the best criterion for discriminating between ECB and ABB samples [47, 48, 52]. Hence we then applied the no-menclature proposed by Frolov et al. [19] for specimens collected on maize and mugwort respectively.

The experimental framework is summarized in Additional file 1: Box S1 under "Experimental design". Briefly, around 20 males and females were released into three types of conditions: (i) "choice" cages, with maize and mugwort plants, (ii) "maize" cages (the usual host for ECB) and (iii) "mugwort" cages (the usual host for ABB). After 72 h, several life history traits, and particularly egg laying, were quantified. After the measurements, all adults were recaptured and the females were used for the acquisition of expression data by RNA sequencing (see figure in Additional file 1: Box S1).

At the time of sampling, we expected most females to be searching for favorable oviposition sites, because in the female adult lifetime (ca. 10–20 days in the laboratory) most females begin to mate and lay eggs 2 days after emergence ([24, 58] and pers. obs.). They continue to lay eggs for several days (with or without additional mating) with egg laying activity peaking at about 10 days. Sampling after 72 h was thus the best compromise between sufficient egg laying and low mortality. This sampling method aimed at capturing the gene activity underlying oviposition behavior and physiology (host searching, egg maturation and laying) by maximizing the number of females ready to oviposit, although some females in the sample were probably not in that target state.

### Dissection, RNA extraction and sequencing

As we primarily targeted genes involved in host preference and host avoidance during oviposition, we focused on females only. Moreover, we assumed that these traits certainly involved chemosensory genes. Chemical reception is mediated by specialized sensory neurons located in the antennae, mouthparts or legs, but chemosensory genes are also involved in non-sensory functions and expressed in other tissues [84]. Thus, we chose to split by tissue samples into head-thorax (HT) and abdomen (ABDO) tissues. HT tissues contained the legs, antennae, and mouthparts while ABDO tissues contained the digestive and genital organs. In addition, the use of separate pools increased the sequencing coverage per pool of tissue. In total, we dissected 220 females on ice to prevent thawing and RNA degradation. Total RNA was extracted from each tissue and individual with the AllPrep DNA/RNA 96 Kit (Qiagen), according to the manufacturer's protocol. RNA quality and quantity were assessed using a NanoDrop (ThermoFisher), and equimolar quantities of individual extracts were pooled to obtain 24 samples corresponding to 2 repetitions × 3 experimental conditions × 2 tissues × 2 moth species. Each pool
contained 7 to 11 females depending on the RNA samples available (Additional file 1: Table S4).

Twenty-four cDNA libraries were constructed with the TruSeq Stranded mRNA Sample Preparation Kit (Illumina) according to the manufacturer's protocol and sent to the MGX platform (Montpellier, France) for single-end 1x50bp sequencing on six lanes of a HiSeq 2000 (Illumina).

## De novo transcriptome assembly and annotation

After sequencing, all reads were subjected to quality control and trimming using Trimmomatic (version 0.25) to remove Illumina sequencing adapters and low quality reads [6]. Ribosomal and bacterial contaminants were identified using sortmeRNA ([33] version 2.0, default parameters), and Bowtie2 ([35]; version 2.2.4) mapping against a bacteria subset dataset extracted from the NCBI database (October 2014).

First, we applied default parameters and the normalization procedure of the Trinity package (version r20140413p1 [28]) to assemble high quality reads of the 12 ECB libraries and of the 12 ABB libraries separately into two comprehensive transcriptomes. Then, for each of the resulting ECB and ABB assemblies, we included a 454 de novo transcript set previously obtained by Gschloessl et al. [27] from four ECB and four ABB adults, respectively. Such meta-assemblies that combine different technologies of sequencing reads have been proved to increase contig length and coverage in other studies [50, 61]. To achieve this goal, we used the EMBOSS [68] splitter program (version 6.6.0) to create short overlapping subsequences (length 300 bp, overlap 200 bp) for each of the 454 transcripts and HiSeq transcript sets, respectively. All 454 and HiSeq subsequences were then assembled with Newbler (454 Life Sciences Corporation, version 2.9) set at 40 bp for minimum overlap length, 97% for minimum overlap identity, and 'Het' (for heterozygote). Finally, we selected the longest transcript as the representative sequence for each gene to circumvent any putative alternative transcripts that might bias the expression analyses.

For functional and GO-term annotations, a custom script was developed to compare ECB-ref and ABB-ref genes against NCBI nr (version May 2016) and UniProtKb ([4], version June 2016) using *blastx* BLAST (blast+, version 2.2.28, [9]) and InterProScan [31].

## Candidate sensory genes

In addition to the reconstruction of the de novo transcripts, we compiled a set of 141 reference genes involved in sensory perception previously described in ECB, ABB and related species of the *Ostrinia* genus [40, 53, 54, 82, 85, 86]. They included odorant receptors (ORs), gustatory receptors (GRs), ionotropic receptors (IRs), sensory neuron membrane proteins, and odorant degrading enzymes. As some

of them may be orthologous between species, or similar between studies, we removed the redundant sequences as identified by reciprocal BLAST analyses. After this filtering step, the sensory dataset was reduced to 139 unique sequences (Additional file 2: Table S11). This gene set, "Senso-ref", was used as an independent gene reference from ABB-ref and ECB-ref to focus on sensory genes.

## Raw read mapping and normalization

Raw reads were mapped onto ABB-ref, ECB-ref and Senso-ref using Bowtie2 ([35], version 2.2.4), set to the default parameters and "sensitive." After mapping, we produced mapped read counts for each gene using the SAMtools program [44]. Transcript-specific variations (the dispersion and logarithmic fold-changes of read counts across the experimental conditions) were estimated with the Bioconductor EdgeR (Empirical Analysis of Digital Gene Expression Data in R) package ([69], version 3.8.6). A negative binomial distribution was chosen as the model for count variations, and transcripts with less than one count per million for at least half of the libraries were removed. Read counts were normalized across the conditions (i.e. across the sequenced libraries) using the *calcNormFactors* EdgeR function and TMM option [70]. Log₂ fold changes (Log2FC) were calculated for each gene between conditions.

In addition, a principal component analysis (PCA) was conducted on the normalized read counts, as implemented in the FactoMineR library of R [38]. The following categorical variables were inputted as supplementary variables: the moth species (ABB vs. ECB), the experimental condition (choice, mugwort, maize), and the tissue (ABDO vs. HT).

## Identification of differentially expressed genes: response to environment and divergence between ECB and ABB

The main source of variation in expression was due, as expected, to the HT vs. ABDO tissues rather than to the species or experimental conditions (Fig. 1). Thus, we analyzed the ABDO and HT libraries separately. Variations in read counts by transcript ($t_i$) were analyzed using the EdgeR linear model to estimate the statistical evidence for the experimental effects. We combined the two explanatory variables *moth species* and *experimental condition* into a so-called "modalities" factor, and applied the following model:

Model 1: $ti$ = modalities.

We then used a contrast analysis to identify *i*) transcripts whose expression varied significantly between ABB and ECB whatever the experimental condition, and *ii*) transcripts whose expression varied significantly between the experimental conditions (maize and mugwort, maize and choice, mugwort and choice) whatever the moth species. Transcripts were considered differentially expressed (DE) across the entire experiment if they had a false discovery rate (FDR) adjusted *p*-value < 0.01.
Gene ontology (GO) enrichment was calculated in R using topGO package (version 2.26.0) from the Bioconductor project [1]. We used the GO annotation described previously as a custom input. We removed GO terms with fewer than five annotated genes (using the nodeSize parameter set to 5). We used the weight algorithm and Fisher's exact test. We reported enriched GO categories with *p-values* < 0.01 for DE genes between ECB and ABB species, and for DE genes between the experimental environments (maize and mugwort, maize and choice, mugwort and choice).

## Categorization of DE genes according to behavioral traits

We further refined the categorization of DE genes to relate them to the preference and avoidance traits and behaviors previously observed in Orsucci et al. [59] (Additional file 1: Box S1). To do this, we split the data set and analyzed the ECB and ABB data separately using a simplified model:

Model 2: ti = experimental condition

We compared *i*) gene expression between mugwort samples and the mean expression of the two other experimental conditions in the ECB sub-data to detect genes potentially involved in maize preference during oviposition in ECB (hereafter referred to as ECB-pref); *ii*) gene expression between maize samples and the mean expression of the two other experimental conditions in the ABB sub-data to detect genes involved in mugwort preference during oviposition in ABB (ABB-pref); and *iii*) gene expression between maize-only and the mean expression of the two other experimental conditions in the ECB sub-data to detect genes involved in mugwort avoidance during oviposition in the mugwort environment in ECB (ECB-avoid).

For all models, the *p-values* were adjusted for multiple testing using the Benjamini-Hochberg procedure (FDR < 0.01). GO enrichment was calculated as previously for the ECB-pref, ABB-pref and ECB-avoid gene lists.

## Additional files

**Additional file 1: Box S1.** Experimental framework of RNAseq sequencing: material, methods and main results on the behavior and life history traits of ECB and ABB during the oviposition phase (derived from [59]). (PDF 574 kb)

**Additional file 2: Table S1.** Details on sequencing reads and contaminants per library. **Table S2.** Top hit species after BLAST analysis, for the ECB-ref transcripts (sheet 1) and ABB-ref transcripts (sheet 2). **Table S3.** Read mapping against ABB-ref transcripts for ABB and ECB samples. **Table S4.** Sequencing characteristics of the RNA pools. **Table S5.** Detailed information for DE genes homologous with virus (A) and retrovirus-like transposon sequences (B). **Table S6.** Detailed information for genes DE between ECB and ABB species. **Table S7.** Fisher's test of enrichment in GO term categories between the DE transcripts (test set) and the whole transcriptome (reference set). **Table S8.** Detailed information for genes DE between experimental environments (maize, choice, mugwort). **Table S9.** Detailed information for genes involved in ECB-pref, ECB-avoid and ABB-pref.

**Table S10.** Detailed information for genes specific to ABB-pref, ECB-pref and ECB-avoid contrasts in HT and ABDO tissues. **Table S11.** Published chemosensory *Ostrinia* spp. genes. **Table S12.** Detailed information for DE chemosensory genes.

**Additional file 3: Figure S1.** Venn diagram for ECB-ref (A and B) and ABB-ref (C and D) transcripts shared between the different experimental conditions for the HT samples (A and C) and ABDO samples (B and D). (PDF 183 kb)

## Abbreviations

ABB: Adzuki bean borer; ABDO: Abdomen; avoid: Avoidance; bp: Base pair; DE: Differentially expressed; ECB: European corn borer; FC: Fold change; FDR: False discovery rate; GLM: Generalized linear model; GO: Gene ontology; GR: Gustatory receptor; HT: Head-thorax; IR: Ionotropic receptor; L4: Fourth larval instar; LRT: Likelihood ratio test; Mbp: Mega base pair; OR: Odorant receptor; PCA: Principal component analysis; pref: Preference; TMM: Trimmed mean M values

## Acknowledgments

We wish to thank Béatrice Ramora, Anne Zanetto and Christophe Reynaud for their technical assistance in plant cultivation at Diascope (INRA experimental research station, Mauguio, France). Our thanks also to Denis Bourguet, Pascal Milesi and Carole Smadja for their critical and helpful comments on the manuscript, and to Elodie de Roux-Labbé for proofreading the manuscript.

## Funding

For this work, RS received support from ANR Adapt-Ome (ANR-13-BSV7-0012-01) for the sampling, study design, sequencing and data analysis of RNA-Seq data, and manuscript writing. In addition, University of Montpellier has funded the PhD grant of MO.

## Availability of data and materials

Raw RNA-Seq reads and the assembled ABB and ECB transcriptomes have been deposited with links to the BioProject accession number PRJNA396660 in the NCBI BioProject database (https://www.ncbi.nlm.nih.gov/bioproject/). The assembled ABB and ECB transcript sets have also been made available on the LepidoDB web site. (http://bipaa.genouest.org/sp/ostrinia_scapulalis/download/transcriptome/ECB_ABB_femaleAdultTranscripts/).

## Authors' contributions

MO, PA and RS conceived the study. MO, PA and AP performed the population sampling, moth rearing, RT experiment, and data collection; SN performed the RNA extractions; and RS regularly assisted in these steps. MV was in charge of plant rearing and set up the experimental conditions at the Diascope site. DS and MR constructed the RNA libraries and performed the HiSeq sequencing. FD and BG generated de novo transcriptome assemblies and annotations. MO and RS performed DE gene identification and statistical analyses and wrote the manuscript, with input from the co-authors. All authors read and approved the final manuscript.

## Ethics approval and consent to participate

No ethics approval was required for the conduct of the experiments in this study.

## Consent for publication

Not applicable.

## Competing interests

The authors declare that they have no competing interests.

## Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Author details

¹CBGP, INRA, CIRAD, IRD, Montpellier SupAgro, Univ Montpellier, Montpellier, France. ²DGIMI, INRA, Univ Montpellier, Montpellier, France. ³DIASCOPE, INRA, Mauguio, France. ⁴Present address: Department of Ecology and Genetics, EBC, Uppsala University, Norbyvägen 18D, 75236 Uppsala, Sweden.

⁵MGX-Montpellier GenomiX, c/o Institut de Génomique Fonctionnelle, 34094 Montpellier Cedex 5, France.
Received: 2 March 2018 Accepted: 30 August 2018

Published online: 20 September 2018

## References

1. Alexa, A. & Rahnenfuhrer, J., 2010. topGO: topGO: Enrichment analysis for Gene Ontology. R package version 2.18.0. October.

2. Awmack CS, Leather SR. Host plant quality and fecundity in herbivorous insects. Annu Rev Entomol. 2006;51(1):309–30.

3. Babikova Z, et al. Arbuscular mycorrhizal fungi and aphids interact by changing host plant quality and volatile emission. Func. 2014;28:375–85.

4. Bateman A, et al. UniProt: a hub for protein information. Nucleic Acids Res. 2015;43(D1):D204–12.

5. Belyayev A. Bursts of transposable elements as an evolutionary driving force. J Evol Biol. 2014;27(12):2573–84.

6. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.

7. Bonoan RE, Al-Wathiqui N, Lewis S. Linking larval nutrition to adult reproductive traits in the European corn borer *Ostrinia nubilalis*. Physiol Entomol. 2015;40(4):309–16.

8. Calcagno V, Thomas Y, Bourguet D. Sympatric host races of the European corn borer: adaptation to host plants and hybrid performance. J Evol Biol. 2007;20(5):1720–9.

9. Camacho C, et al. BLAST+: architecture and applications. BMC Bioinform. 2009;10(1):421.

10. Carrasco D, Larsson MC, Anderson P. Insect host plant selection in complex environments. Curr Opin Insect Sci. 2015;8:1–7.

11. Casacuberta E, González J. The impact of transposable elements in environmental adaptation. Mol Ecol. 2013;22(6):1503–17.

12. Coyne JA, Orr HA. Speciation in Sinauer. Sunderland: Sinauer Associates; 2004.

13. Derridj S, et al. Role of nutrients found in the phylloplane, in the insect host-plant selection for oviposition. In: 8th International Symposium on Insect-plant Relationships. Dordrecht: Kluwer Academic Publishers; 1992. p. 139–40.

14. Drès M, Mallet J. Host races in plant-feeding insects and their importance in sympatric speciation. Philos Trans R Soc Lond Biol Sci. 2002;357(1420):471–92.

15. Dupuy F, et al. A background of a volatile plant compound alters neural and behavioral responses to the sex pheromone blend in a moth. Front Physiol. 2017;8:79.

16. Eyres I, et al. Differential gene expression according to race and host plant in the pea aphid. Mol Ecol. 2016;25(17):4197–215.

17. Forister ML, et al. Revisiting the evolution of ecological specialization, with emphasis on insect–plant interactions. Ecology. 2012;93(5):981–91.

18. Foster SP, Howard AJ. The effects of mating, age at mating, and plant stimuli, on the lifetime fecundity and fertility of the generalist herbivore *Epiphyas postvittana*. Entomol Exp Applicata. 1999;91(2):287–95.

19. Frolov AN, Bourguet D, Ponsard S. Reconsidering the taxonomy of several *Ostrinia* species in the light of reproductive isolation: a tale for Ernst Mayr. Biol J Linn Soc. 2007;91(1):49–72.

20. Funk DJ. Isolating a role for natural selection in speciation: host adaptation and sexual isolation in *Neochlamisus bebbianae* leaf beetles. Evolution. 1998; 52(6):1744–59.

21. Futuyma DJ, Moreno G. The evolution of ecological specialization. Annu Rev Ecol Evol Syst. 1988;19:207–33.

22. Gavrillets S. Fitness landscapes and the origin of species. Princeton and Oxford: P. U. Press; 2004.

23. Glover TJ, Tang XH, Roelofs WL. Sex pheromone blend discrimination by male moths from E and Z strains of European corn borer. J Chem Ecol. 1987; 13(1):143–51.

24. Gohari H, Hawlitzky N. Activité reproductrice de la pyrale du maïs, *Ostrinia nubilalis* Hbn. (Lep. *Pyralidae*) à basses températures constantes. Agronomie. 1986;6(10):911–7.

25. Grippenberg S, et al. A meta-analysis of preference-performance relationships in phytophagous insects. Ecol Lett. 2010;13(3):383–93.

26. Gross P. Insect behavioral and morphological defenses against parasitoids. Annu Rev Entomol. 1993;38:251–73.

27. Gschloessl B, et al. De novo transcriptomic resources for two sibling species of moths: *Ostrinia nubilalis* and *O. scapulalis*. BMC Res Notes. 2013;6:73.

28. Haas BJ, et al. Reference generation and analysis with trinity. Nat Protoc. 2014;8(8):1–43.

29. Hoang K, Matzkin LM, Bono JM. Transcriptional variation associated with cactus host plant adaptation in *Drosophila mettleri* populations. Mol Ecol. 2015;24(20):5186–99.

30. Jaenike J. Host specialization in phytophagous insects. Annu Rev Ecol Syst. 1990;21:243–73.

31. Jones P, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30(9):1236–40.

32. Koenig C, et al. The plastic response of *Manduca sexta* to host and non-host plants. Insect Biochem Mol Biol. 2015;63:72–85.

33. Kopylova E, Noé L, Touzet H. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics. 2012;28:3211–7.

34. Landolt PJ, Phillips TW. Host plant influences on sex pheromone behavior of phytophagous insects. Annu Rev Entomol. 1997;42:371–91.

35. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie2. Nat Methods. 2012;9(4):357–9.

36. Larsson S, Häggström H, Denno R. Preference for protected feeding sites by larvae of the willow-feeding leaf beetle *Galerucella lineola*. Ecol Entomol. 1997;22:445–52.

37. Lassance JM. Journey in the *Ostrinia* world: from pest to model in chemical ecology. J Chem Ecol. 2010;36(10):1155–69.

38. Lê S, Josse J, Husson F. FactoMineR: an R package for multivariate analysis. J Stat Softw. 2008;25(1):1–18.

39. Leahy TC, Andow DA. Egg weight, fecundity, and longevity are increased by adult feeding in *Ostrinia nubilalis* (Lepidoptera: Pyralidae). Ann Entomol Soc Am. 1994;87(3):342–9.

40. Leary GP, et al. Single mutation to a sex pheromone receptor provides adaptive specificity between closely related moth species. Proc Natl Acad Sci U S A. 2012;109(35):14081–6.

41. Leather SR, Watt AD, Barbour D. The effect of host-plant and delayed mating on the fecundity and lifespan of the pine beauty moth, *Panolus flammea* (Denis and Schiffermuller) (Lepidoptera: Noctuidae): their influence on population dynamics and relevance to pest management. Bull Entomol Res. 1985;75(4):641–51.

42. Leppik E, Frérot B. Volatile organic compounds and host-plant specialization in European corn borer E and Z pheromone races. Chemoecology. 2012; 22(2):119–29.

43. Li H, et al. Identification of differential expression genes associated with host selection and adaptation between two sibling insect species by transcriptional profile analysis. BMC Genomics. 2013;14(1):582.

44. Li H, et al. The sequence alignment / map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.

45. Lillehoj EB. Aflatoxin: an ecologically elicited genetic activation signal. In: Smith JE, Henderson RS, editors. Mycotoxins and animal foods: CRC Press, Inc.; 1991. p. 1–36.

46. Ma L, et al. Identification and comparative study of chemosensory genes related to host selection by legs transcriptome analysis in the tea geometrid *Ectropis obliqua*. PLoS One. 2016;11(3):1–25.

47. Malausa T, et al. Molecular differentiation at nuclear loci in French host races of the European corn borer (*Ostrinia nubilalis*). Genetics. 2007a;176(4): 2343–55.

48. Malausa T, et al. Genetic structure and gene flow in French populations of two *Ostrinia* taxa: host races or sibling species? Mol Ecol. 2007b;16(20):4210–22.

49. Malausa T, et al. Differences in oviposition behaviour of two sympatric sibling species of the genus *Ostrinia*. Bull Entomol Res. 2008;98(2):193–201.

50. Marchant A, et al. Comparing de novo and reference-based transcriptome assembly strategies by applying them to the blood-sucking bug Rhodnius prolixus. Insect Biochem Mol Biol. 2016;69:25–33.

51. Matsubayashi, K.W., Ohshima, I. & Nosil, P., 2010. Ecological speciation in phytophagous insects. Entomol Exp Applicata. Matsubayashi, K.W., Ohshima, I. & Nosil, P., 2010. Ecological speciation in phytophagous insects. Entomologia Experimentalis et Applicata. 134(1), pp.1–27.

52. Midamegbe A, et al. Scanning the European corn borer (*Ostrinia spp.*) genome for adaptive divergence between host-affiliated sibling species. Mol Ecol. 2011;20(7):1414–30.

53. Miura N, et al. A male-specific odorant receptor conserved through the evolution of sex pheromones in *Ostrinia* moth species. Int J Biol Sci. 2009; 5(4):319–30.

54. Miura N, et al. Broadly and narrowly tuned odorant receptors are involved in female sex pheromone reception in *Ostrinia* moths. Insect Biochem Mol Biol. 2010;40(1):64–73.

55. Molnár BP, et al. Electrophysiologically-active maize volatiles attract gravid female European corn borer, *Ostrinia nubilalis*. J Chem Ecol. 2015;41(11):997–1005.
56. Moreau J, Benrey B, Thiery D. Assessing larval food quality for phytophagous insects: are the facts as simple as they appear? Funct Ecol. 2006;20(4):592–600.

57. Nosil, P. et al. Natural selection and divergence in mate preference during speciation. *Genetica*. 2007.

58. Orsucci, M., 2015. Rôle de la spécialisation à la plante hôte et de l'isolement reproducteur dans la divergence de lépidoptères ravageurs de cultures. PhD Thesis, Université de Montpellier.

59. Orsucci M, et al. Host specialisation in two phytophagous moth species: a multiform pattern involving attraction avoidance and performance. J Evol Biol. 2016;29:114–25.

60. Orsucci M, et al. Larval gene expression and host plant specialization in two sibling moth species (Ostrinia spp.). BMC Genomics. 2018;19:265.

61. Perera OP, et al. Generation of a transcriptome in a model lepidopteran pest, *Heliothis virescens*, using multiple sequencing strategies for profiling midgut gene expression. PLoS One. 2015;10(6):e0128563.

62. Pregitzer P, et al. Plant odorants interfere with detection of sex pheromone signals by male *Heliothis virescens*. *Front Cell Neurosci*. 2012;6:42.

63. Ragland GJ, et al. Differences in performance and transcriptome-wide gene expression associated with *Rhagoletis* (Diptera: Tephritidae) larvae feeding in alternate host fruit environments. Mol Ecol. 2015;24(11): 2759–76.

64. Raymond B, et al. Host plant species can influence the fitness of herbivore pathogens: the winter moth and its nucleopolyhedrovirus. Oecologia. 2002; 131(4):533–41.

65. Reisenman CE, Lei H, Guerenstein PG. Neuroethology of olfactory-guided behavior and its potential application in the control of harmful insects. *Front Physiol. 2016;7*(June):1–21.

66. Reisenman CE, Riffell JA. The neural bases of host plant selection in a Neuroecology framework. *Front Physiol. 2015;6*(Aug):1–17.

67. Renwick J. Oviposition behavior in Lepidoptera. Annu Rev Entomol. 1994; 39(1):377–400.

68. Rice P, Longden I, Bleasby A. EMBOSS: The European Molecular Biology Open Software Suite. Trends Genet. Elsevier. 2000;16:276–7.

69. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. *Bioinformatics*. 2010;26(1):139–40.

70. Robinson MD, Oshlack A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010;11.

71. Sarov-Blat L, et al. The Drosophila takeout gene is a novel molecular link between circadian rhythms and feeding behavior. Cell. 2000;101(6): 647–56.

72. Schlüter D. Ecological character displacement in adaptive radiation. Am Nat. 2000;156(s4):S4–S16.

73. Schlüter D. Ecology and the origin of species. *Trends Ecol Evol*. 2001;16(7): 372–80.

74. Schurr K, Holdaway FG. Olfactory responses of female *Ostrinia nubilalis* (Lepidoptera: Pyraustinae). Entomol Exp Appl. 1970;13(6278):455–61.

75. Scriber JM, Slansky F. The nutritional ecology of immature insects. Annu Rev Entomol. 1981;26:183–211.

76. Silva-Brandão KL, et al. Transcript expression plasticity as a response to alternative larval host plants in the speciation process of corn and rice strains of *Spodoptera frugiperda*. BMC Genomics. 2017;18(1):792.

77. Simon, J.-C. et al. Genomics of adaptation to host-plants in herbivorous insects. Brief Funct Genom. 2015. 1–11.

78. Solé J, et al. Behavioural and electrophysiological responses of the European corn borer *Ostrinia nubilalis* to host-plant volatiles and related chemicals. *Physiol Entomol. 2010;35*(4):354–63.

79. Städler E. Plant chemical cues important for egg deposition by herbivorous insects. In: Chemoecology of Insect Eggs and Egg Deposition. Hilker M, Meiners T. eds. Berlin: Blackwell Publishing; 2002. p. 171-204.

80. Trona F, et al. Neural coding merges sex and habitat chemosensory signals in an insect herbivore. Proc. R Soc B Biol Sci. 2013;280:20130267.

81. Via S, Hawthorne DJ. The genetic architecture of ecological Specialization : correlated gene effects on host use and. Am Nat. 2002; 159(3):S76–88.

82. Wanner KW, et al. Sex pheromone receptor specificity in the European corn borer moth, *Ostrinia nubilalis*. PLoS One. 2010;5(1):e8685.

83. Wood TK, Keese MC. Host-plant-induced assortative mating in *Enchenopa* treehoppers. Evolution. 1990;44(443):619–28.

84. Wu Z, et al. Antennal and abdominal transcriptomes reveal chemosensory genes in the Asian Citrus psyllid, Diaphorina citri. PLoS One. 2016;11(7):1–23.

85. Yang B, et al. Identification of candidate odorant receptors in asian corn borer *Ostrinia furnacalis*. PLoS One. 2015;10(3):1–19.

86. Yasukochi Y, et al. Sex-linked pheromone receptor genes of the European corn borer, *Ostrinia nubilalis*, are in tandem arrays. PLoS One. 2011;6(4):1–10.

87. Zdybicka-Barabas A, Cytryńska M. Apolipophorins and insects immune response. Invertebr Surviv J. 2013;10:58–68.

88. Zhang T, He K, Wang Z. Transcriptome comparison analysis of *Ostrinia furnacalis* in four developmental stages. Sci Rep. 2016;6(1):35008.

Ready to submit your research? Choose BMC and benefit from:

• fast, convenient online submission

• thorough peer review by experienced researchers in your field

• rapid publication on acceptance

• support for research data, including large and complex data types

• gold Open Access which fosters wider collaboration and increased citations

• maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_10.png" /></figure>