Original Research ARTICLE
Comparing the Expression of Olfaction-Related Genes in Gypsy Moth (Lymantria dispar) Adult Females and Larvae from One Flightless and Two Flight-Capable Populations
- 1Department of Environmental Systems Science, ETH Zürich, Zürich, Switzerland
- 2Department of Evolutionary Neuroethology, Max Planck Institute for Chemical Ecology, Jena, Germany
- 3Institute of Fundamental Sciences, College of Sciences, Massey University, Palmerston North, New Zealand
In insects, flight and sophisticated olfactory systems go hand in hand and are essential to survival and evolutionary success. Females of many Lepidopteran species have secondarily lost their flight ability, which may lead to changes in the olfactory capabilities of both larval and adult stages. The gypsy moth, Lymantria dispar, an important forest pest worldwide, is currently undergoing a diversification process with three recognized subspecies: the Asian gypsy moth (AGM), Lymantria dispar asiatica; the Japanese gypsy moth (JGM), Lymantria dispar japonica; and the European gypsy moth (EGM), Lymantria dispar dispar. Females of EGM populations from North America have lost their flight capacity whereas the JGM and AGM females are flight capable, making this an ideal system to investigate the relationship between flight and olfaction. We used next-generation sequencing to obtain female antennal and larval head capsule transcriptomes in order to (i) investigate the differences in expression of olfaction-related genes among populations; (ii) identify the most similar protein sequences reported for other organisms through a BLAST search, and (iii) establish the phylogenetic relationships of these sequences with respect to other insect species. Using this approach, we identified 115 putative chemosensory genes belonging to five families of olfaction-related genes. A principal component analysis (PCA) revealed that the gene-expression patterns of female antennal transcriptomes from different subspecies were more similar to one another than to the larval head capsules of their respective subspecies supporting strong chemosensory differences between the two developmental stages. An analysis of the shared and exclusively expressed genes for three populations shows no evidence that loss of flight affects the number or type of genes being expressed. These results indicate either (a) that loss of flight does not impact the olfactory gene repertoire or (b) that the secondary loss of flight in American EGM populations may be too recent to have caused major changes in the genes being expressed. However, we found higher expression values for most olfaction-related genes in EGM females, suggesting that differences in transcription rates could be an adaptation of flightless females to their chemical environment. Differences in olfactory genes and their expression in the larvae appear to be unrelated to the flight ability of adult females and are likely adaptations to different ecological pressures.
Flight is a leading factor contributing to the evolutionary success of insect species, enabling them to locate food and shelter, avoid predation and competition, and search for optimal oviposition sites for their offspring (Barbosa et al., 1989; Sattler, 1991; Hunter, 1995). Since host–plant location and oviposition in herbivorous insects are largely mediated by chemical cues (Bruce et al., 2005; Bruce and Pickett, 2011; Mescher and De Moraes, 2015), one would expect the evolution of flight to be accompanied by the development of sophisticated olfactory systems. New evidence even suggests that the odorant receptor family (OR), central to the olfactory systems of highly derived insects, emerged around the same time as flight (Missbach et al., 2014; Ioannidis et al., 2017). Furthermore, manipulation of OR-based odor detection in Drosophila also indicates that ORs play an important role in flight orientation (Getahun et al., 2016).
The females of many Lepidopteran species have secondarily lost their ability to fly, shifting the responsibility of host selection partly or entirely to the larvae (Barbosa et al., 1989; Sattler, 1991; Hunter, 1995). In this context, it is interesting to investigate whether the loss of flight has an impact on the olfaction of adults and larvae. The gypsy moth Lymantria dispar is one of the most important forest pest species worldwide, currently undergoing a diversification process involving the loss of flight by females of some populations (Schweitzer, 2004; Pogue and Schaeffer, 2007). These features make L. dispar an ideal model to explore changes in expression of olfaction-related genes that are associated with flight ability.
The first chemosensory proteins (CSPs) from adult L. dispar were identified as early as 1989 (Vogt et al., 1989, 1991). Thereafter, Plettner and coworkers have made great contributions to our understanding of olfaction in this species, in particular concerning the structure and function of its pheromone binding proteins (Kowcun et al., 2001; Honson et al., 2003; Honson and Plettner, 2006; Plettner and Gries, 2010; Gong and Plettner, 2011; Yu and Plettner, 2013). Recently, the L. dispar olfactory co-receptor (ORCO), a crucial component of olfactory receptor complexes, has been identified (Vosshall and Hansson, 2011; Lin et al., 2015). However, knowledge about olfaction-related proteins and the genes encoding them remains fragmentary for this species.
The gypsy moth is a highly polyphagous herbivore, capable of causing severe and widespread outbreaks in temperate Holarctic regions. At present, there are three recognized subspecies: the Asian Gypsy moth (AGM) Lymantria dispar asiatica, the Japanese Gypsy moth (JGM) Lymantria dispar japonica, and the European Gypsy moth (EGM) Lymantria dispar dispar (which encompasses both European and North American Gypsy moth populations). European Gypsy moth females from North American populations are flightless, possibly due to a founder effect associated with their introduction from Europe in the mid nineteenth century. In contrast, the Asian and Japanese females can fly and disperse over extended distances (Barlow, 2004; NBII, 2011; APHIS, 2013).
The loss of flight in the EGM females restricts their ability to make host-plant choices, transferring the responsibility to the larvae, which disperse either passively through ballooning in the early instars or actively by crawling in the late instars (Capinera and Barbosa, 1976; Lance and Barbosa, 1981, 1982). The extent to which flight capable females are involved in host-plant choices is not yet fully understood, but evidence suggests that both AGM and JGM actively disperse and display oviposition preferences under field conditions (Baranchikov, 1989; Sasaki et al., 2016).
Several efforts have been made to better understand the taxonomic and biogeographic distribution of female flight ability, as well as its heritability and phenotypic plasticity (Keena et al., 2001, 2007, 2008, 2010). However, no studies have yet documented variation in the odor perception systems of L. dispar subspecies, despite the likelihood that such differences may accompany the loss of female flight. Therefore, the aims of this study were to (i) Investigate the differences in expression of olfaction-related genes among populations, (ii) identify the most similar protein sequences reported for other organisms through a BLAST search, and (iii) establish the phylogenetic relationships of these sequences with respect to other model insect species, most of which have fully sequenced genomes.
To fulfill these aims we focused on five groups of chemosensory gene families: odorant receptors (ORs), odorant binding proteins (OBPs), CSPs, gustatory receptors (GRs), and ionotropic receptors (IRs). ORs are expressed in the cell membranes of olfactory sensory neurons (OSNs) and are responsible for the detection of odor molecules (Sanchez-Gracia et al., 2009). In general, OSNs will express either ORs or IRs, with the latter mostly tuned to compounds of lower molecular weight (Hallem et al., 2004, 2006; Benton et al., 2009; Silbering et al., 2011). All analyzed Lepidoptera species possess more OR than IR types (Croset et al., 2010; Koenig et al., 2015; van Schooten et al., 2016), and these play a role in the detection of plant volatiles as well as pheromones (Nakagawa et al., 2005; Grosse-Wilde et al., 2006, 2007; Tanaka et al., 2009). In insects, OSNs associated with basiconic or trichoid sensilla express one OR gene, along with the co-receptor ORCO, which is highly conserved and broadly expressed (Krieger et al., 2003; Touhara and Vosshall, 2009). Insect ORs are seven-transmembrane domain receptors with inverted membrane topology and are not phylogenetically related to vertebrate ORs (Benton et al., 2006).
OBPs contribute to the sensitivity of the olfactory system by binding, solubilizing and transporting odorants through the sensillar lymph (Leal, 2013). CSPs are likely to perform similar roles in chemical communication of insects as OBPs, but unlike these are also expressed in non-chemosensory tissues, and for this reason have been hypothesized to serve additional, as yet undiscovered, functions (Pelosi et al., 2005). Recent evidence suggests that OBPs are an adaptation to the detection of hydrophobic volatiles that became available as olfactory cues in the course of insect terrestrialization (Missbach et al., 2015); however, results in Drosophila suggest a different function for some OBPs (Larter et al., 2016). Structurally, insect OBPs and CSPs generally contain α-helical domains, but folded in two different patterns (Sandler et al., 2000; Lartigue et al., 2002; Tegoni et al., 2004).
GRs are typically expressed in gustatory receptor neurons (GRNs) within the taste sensillae in the mouthparts and are known to detect sugars, bitter compounds and non-volatile pheromones (Montell, 2013). However, some GR genes are also expressed in the antennae, suggesting that some members of this gene family may have an olfactory function (Hallem et al., 2006). This is further supported by the discovery of two GRs in Drosophila that act in the detection of CO2 (Yao and Carlson, 2010). GR proteins are highly divergent in sequence, sharing as little as 8% amino acid identity across insect species, and it has been hypothesized that the GR gene family is an ancient chemoreceptor family from which insect OR genes subsequently evolved (Robertson et al., 2003; Hallem et al., 2006; Benton, 2015).
IRs are also involved in chemoreception and comprise a large and highly diverse gene family closely related to ionotropic glutamate receptors (iGluR), typically present in the OSNs associated with the coeloconic sensillae in the antennae (Rytz et al., 2013). Recent reports suggest there are multiple variant IRs with different ligand-binding domains that lack the characteristic glutamate-interacting residues (Benton et al., 2009). Unlike ORs, which are exclusively found in pterygote insects, IRs are present in all protostome species studied so far and may have evolved as long as 550–850 million years ago (Croset et al., 2010; Missbach et al., 2014). Similar iGluR-like genes are also present in plants, animals and prokaryotes, indicating that this is an important and ancient group of chemoreceptors (Benton et al., 2009; Rytz et al., 2013).
Materials and Methods
Insects were provided as egg masses by Hannah Nadel, Supervisory Entomologist of the Animal and Plant Health Inspection Service (APHIS) of the United States Department of Agriculture (USDA). All egg masses came from laboratory cultures that had been maintained using carefully designed mating protocols to avoid the deleterious effects of inbreeding depression, details on the rearing system utilized for these colonies can be found in (Bell et al., 1981). Upon hatching larvae were fed ad libitum on artificial wheat germ diet prepared according to manufacturer's instructions (MP Biomedicals LLC, Illkirch, France) and food was replaced twice per week. Caterpillars, pupae, and adult moths were maintained in a climate chamber at 20°C, 60% relative humidity and 16/8 h photoperiod.
The European gypsy moth (EGM) culture (Lymantria dispar dispar) originated from flightless L. dispar populations collected in New Jersey (US). The Japanese gypsy moth (JGM) culture (Lymantria dispar japonica), originated from flight-capable populations coming from the Northern Iwate district and Takizawa, Morika, Nishine (Japan). The AGM culture (Lymantria dispar asiatica) originated from flight-capable populations coming from the Primorskiy Krai ports (Vostochnyy, Slavyanka, Vladivostok, Nadhodka) in Russia.
RNA extraction was performed following the same procedure as in Koenig et al. (2015), with minor changes as outlined below. Antennae of 50 adult female moths (1–2 days old) from each population were excised from the base of the antennal sclerite. Head capsules from 50 fifth instar larvae from each population were cut at the division point with the prothorax. Tissues were transferred to an Eppendorf tube, cooled with liquid nitrogen and stored at −86°C until extraction. For extraction, tissues were transferred into RL buffer (innuPREP RNA Mini Kit, Analytik Jena, Jena, Germany) and homogenized using a TissueLyser (Qiagen, Hilden, Germany). The resultant homogenate was used with the innuPREP RNA Mini Kit (Analytik Jena, Jena, Germany) following the manufacturers protocol.
Sequencing, Assembly, and Annotation
Total RNA was sent to the Max Planck Genome Centre Cologne (Germany) for construction of TruSeq libraries and subsequent sequencing on an Illumina HiSeq3000. Read data was trimmed and cleaned by the Genome Centre using standard protocols. The resulting Illumina reads were assembled with CLC Genomics Workbench 8 (CLCbio), using the de novo algorithm and default parameters. Annotation was performed using Blast2GO 3 (Conesa et al., 2005; Götz et al., 2008). Additionally, assembled transcripts belonging to target chemosensory families (OR, OBP, IR, GR, and CSP) were identified by comparison against custom, manually curated databases created using the available literature on other Lepidopteran species (Wanner and Robertson, 2008; Grosse-Wilde et al., 2011; Heliconius-Genome-Consortium, 2012; Briscoe et al., 2013; Koenig et al., 2015).
Each of the predicted protein sequences was further compared to available sequences using the blastp algorithm and the nr database (NCBI)1 to identify the most similar sequence, the organism expressing it and its putative function. We only report sequences yielding significant (E < 0.05) similarity values.
Alignments and Phylogenetic Trees
Protein sequences conceptually translated from the assembled transcripts were aligned with homologs from Bombyx mori, Danaus plexippus, Heliconius melpomene, and Manduca sexta (Wanner and Robertson, 2008; Grosse-Wilde et al., 2011; Heliconius-Genome-Consortium, 2012; Briscoe et al., 2013; Koenig et al., 2015). In the case of GRs and ORs, sequences from the waterflea Daphnia pulex were also included as an outgroup (Peñalva-Arana et al., 2009). For the CSPs and OBPs we included sequences from the Jumping Bristletail Lepismachilis y-signata and the Firebrat Thermobia domestica (Missbach et al., 2015). In the case of IRs (and Glu-Rs) sequences from Drosophila melanogaster have been included (Rytz et al., 2013).
For this purpose, we used MAFFT version 7 (Katoh et al., 2002; Katoh and Standley, 2013) with the “-auto” option. Phylogenetic trees were derived using the program FastTree-2, which uses the maximum likelihood method with a Shimodaira-Hasegawa test to estimate branch support values (Price et al., 2010). Figures were prepared for publication using the FigTree software 1.4.1 (Rambaut, 2007, 2012).
Some transcripts, corresponding to pheromone binding proteins (PBPs), general odorant binding proteins (GOBPs), glutamate receptors (Glu-Rs and Nmdars = N-methyl-D-aspartate receptors) and the ORCO were tentatively labeled following the naming code of closely related sequences.
Quantification of Gene Expression
For the quantification of gene expression levels in the respective tissues/subspecies, the annotated assemblies were used as a template, mapping the raw reads and performing RPKM analysis in CLC Genomics Workbench 8 using default settings. PCA plots were based on normalized count data that was transformed using the regularized log function implemented by the R package DESeq2 (doi: 10.1186/s13059-014-0550-8).
Gene Identification and Expression Patterns for the Three L. dispar Populations
We used next generation sequencing to obtain transcriptome assemblies of adult female antennae and larval head capsules from EGM, AGM, and JGM populations of L. dispar. The assemblies contained 28,004, 33,208, and 30,820 unique transcripts for EGM, AGM, and JGM populations, respectively. Blastx of the assembled transcripts to the NCBI refseq protein database revealed that that 46.3% (EGM), 52.6% (AGM), and 49.3% (JGM) had high homology (E < 1e-5) to previously characterized proteins at NCBI. To ascertain the transcript coverage of each assembly, we used Blastx to find the proportion of B. mori proteins that aligned in a high scoring alignment. We chose B. mori because it has one of the best characterized genomes of the Lepidoptera. This analysis showed that an average L. dispar transcript encodes just over half the expected protein sequence based on the best blastx hit to B. mori, possibly due to a high proportion of partial sequences (Supplementary Figure 1).
From the assembled transcripts we were able to identify 115 putative chemosensory transcripts belonging to the five families, 22 CSPs, 32 OBPs (including 2 GOBPs, and four pheromone binding proteins), 11 GRs, 33 ORs, and 16 IRs (Table 1). In addition we report 6 glutamate receptors (which are not chemosensory receptors) (Table 1). Our results show that 42 olfaction-related genes are found in at least one population in both female antennae and larval head capsules, 52 are exclusive to the female antennae, and 20 to the larval head capsules. A large contribution to the transcripts that are exclusive to the female antennae comes from the ORs (Table 1). Figure 1 depicts the differences and commonalities in gene expression (presence/absence) among the three populations for each chemosensory gene family.
Table 1. Normalized expression values (RPKM) for different families of chemosensory genes in the female antennae and larval head capsules of three populations of the gypsy moth Lymantria dispar.
Figure 1. Shared and exclusively expressed genes for three populations of the Gypsy moth for different classes of olfaction-related gene families in both female antennae and larval head capsules.
A Principal Component Analysis (PCA) comparing the gene expression patterns (Normalized gene expression values—RPKM) for the antennal and head capsule transcriptomes revealed that female antennal transcriptomes were clustered, being more similar to one another than to the larval transcriptomes of the same population. In contrast, larval transcriptomes were not clustered, but separated along the second component axis (Figure 2).
Figure 2. Principal component analysis (PCA), comparing the normalized gene expression patterns for the female antennae (ant) and the larval head capsule (head) transcriptomes of three subspecies of the Gypsy moth (L. dispar). The flight capable AGM (Asian gypsy moth) and JGM (Japanese gypsy moth), and the flightless EGM (European gypsy moth).
Best Match with Other Protein Sequences
After performing Blast searches with the individual protein sequences, we found that most putative L. dispar CSPs have a high sequence homology with those already published for a number of Lepidopteran species, the majority of which are Noctuid moths belonging to the genera Helicoverpa or Spodoptera (Table 2).
Table 2. List of L. dispar transcripts putatively involved in chemoreception, and characterization for the best hit after comparison with available protein sequences using the BlastP algorithm.
Phylogenetic Positioning of Putative Protein Sequences
We constructed phylogenetic trees from alignments of the L. dispar CSPs with other published sequences from model insect species (B. mori, H. melpomene, M. sexta, and D. plexippus, D. pulex, L. y-signata, T. domestica, and D. melanogaster). Fasta sequences used to construct the trees can be found in Supplementary Files 1–4.
The phylogenetic trees showed that CSPs aligned well within the published sequences, but in a few cases formed clusters containing only L. dispar sequences (e.g., CSPs 3, 16, 13, 21; CSPs 2, 8, 14, and 10; CSPs 5, 4, 1) (Figure 3). For OBPs most sequences were closely related to those reported for the model species, except OBPs 10, 3, 1, 2, 6, 8, 13, and 14 forming a branch unique to L. dispar and a few others forming single nodes (e.g., OBP11) (Figure 4).
Figure 3. Maximum likelihood dendrogram based on protein sequences of candidate chemosensory proteins (CSPs). Included are the putative sequences for Lymantria dispar (Ldis) plus those available for model Lepidoptera species Bombyx mori (Bmor), Danaus plexippus (Dple), Heliconius melpomene (Hmel), and Manduca sexta (Msex), we also included sequences from the Jumping Bristletail Lepismachilis y-signata (Lsig) and the Firebrat Thermobia domestica (Tdom).
Figure 4. Maximum likelihood dendrogram based on protein sequences of candidate odorant binding proteins (OBPs). Included are the putative sequences for Lymantria dispar (Ldis) plus those available for model Lepidoptera species Bombyx mori (Bmor), Danaus plexippus (Dple), Heliconius melpomene (Hmel), and Manduca sexta (Msex), we also included sequences from the Jumping Bristletail Lepismachilis y-signata (Lsig) and the Firebrat Thermobia domestica (Tdom).
In the case of the GRs and ORs, sequences are remarkably well nested within those of model species. Of particular interest is GR2 making a single node, and branch containing GR1 and ORs 2 and 4 unique to L. dispar (Figure 5). For most IRs, we found that the candidate gene sequences were partially aligned with those of the model species, with a few sequences (e.g., IR2) forming single nodes. Most Glu-Rs formed a branch unique to L. dispar (Figure 6).
Figure 5. Maximum likelihood dendrogram based on protein sequences of candidate gustatory receptors (GRs) and odorant receptors (ORs). Included are the putative sequences for Lymantria dispar (Ldis) plus those available for model Lepidoptera species Bombyx mori (Bmor), Danaus plexippus (Dple), Heliconius melpomene (Hmel), and Manduca sexta (Msex), we also included sequences from the waterflea Daphnia pulex (Dpu) as an outgroup.
Figure 6. Maximum likelihood dendrogram based on protein sequences of candidate ionotropic receptors (IRs). Included are the putative sequences for Lymantria dispar (Ldis) plus those available for model Lepidoptera species Bombyx mori (Bmor), Danaus plexippus (Dple), Heliconius melpomene (Hmel) and Manduca sexta (Msex), we also included sequences from the fruit fly Drosophila melanogaster (Dmel).
In recent years, considerable progress has been made in our understanding of insect olfaction. Antennal transcriptomes are available for insect species belonging to several orders, including Diptera, Coleoptera, and Lepidoptera (Grosse-Wilde et al., 2011; Andersson et al., 2013; Rinker et al., 2013; Leitch et al., 2015; Zhang et al., 2015). Within the Lepidoptera, the transcriptomes of model species such as B. mori, D. plexippus, H. melpomene, H. virescens, and M. sexta have been thoroughly investigated (Krieger et al., 2003, 2004; Nakagawa et al., 2005; Wanner et al., 2007; Wanner and Robertson, 2008; Tanaka et al., 2009; Briscoe et al., 2013; Koenig et al., 2015; van Schooten et al., 2016).
This knowledge is rapidly expanding to other economically important species like Helicoverpa armigera, Cydia pomonella, and Spodoptera littoralis, where it could greatly aid in improving already existing and developing new semiochemical-based management strategies (Bengtsson et al., 2012; Jacquin-Joly et al., 2012; Liu et al., 2012). This report represents an expansive characterization of the chemosensory transcripts and their encoded proteins of L. dispar, and increases the number of available olfactory-related sequences for Lepidopteran species of agricultural relevance. Our results may also help unveil how the expression of chemosensory genes changes throughout insect development and as a result of speciation processes, together with similar reports for other species (Poivet et al., 2013; Zhang et al., 2014, 2015; Walker et al., 2016).
We identified a total of 115 putative olfactory transcripts for L. dispar. This number is similar to the one reported on a previous study comparing S. littoralis adult antennae and larval head capsules (127) (Poivet et al., 2013), and another study investigating the adult antennal transcriptome of H. armigera (131), and H. assulta (129) (the latter did not include GRs and we excluded sensory neuron membrane proteins from the total count) (Zhang et al., 2014). The conserved number of olfaction-related genes suggests a core group of genes control olfaction in moth species belonging to the superfamily Noctuoidea (Kristensen et al., 2007; Zahiri et al., 2011). Given the similar number of genes across these species, we could speculate that olfactory differences emerge as a product of functional diversification, while the genes themselves are products of duplication. However, at this stage we can't rule out specific expansions of certain gene clusters balanced by contraction in others, and the evaluation of this possibility must await a more detailed understanding of olfactory differences in the Noctuoidea.
A study investigating expression patterns between adults and larvae of S. littoralis found that adults and larvae express similar numbers of OBPs and CSPs, while the caterpillar OR and IR repertoires were much smaller than the adult ones, and some GRs were found to be adult-specific (Poivet et al., 2013). We also encountered a similar number of CSPs being expressed in both stages and reduced IR, OR, and GR repertoires in the larval stages. However, in contrast to the previous study, we found that larvae had a higher CSP repertoire than adults including eight larval-specific genes. This pattern could reflect species-specific adaptations since in S. littoralis host-plant selection is mainly accomplished by adult females, who make suitable choices for the larvae as eclosion occurs rapidly after egg laying (Anderson and Alborn, 1999; Proffit et al., 2015). In contrast, L. dispar eggs of all populations undergo an overwintering process accompanied by changes in the distribution and quality of the resources from oviposition until larval hatching (Barbosa et al., 1989; Sattler, 1991; Hunter, 1995). Therefore, larval stages need to make host-choices to a greater or lesser extent, which may explain the observed differences in the number of CSP genes being expressed in the larval stages.
The strong reduction in the ORs in larvae vs. adults seems to be commonplace in insects and has been reported for a number of species, including D. melanogaster, Aedes aegypti, M. sexta, and B. mori (Hallem et al., 2004; Kreher et al., 2005; Bohbot et al., 2007; Tanaka et al., 2009; Koenig et al., 2015). In L. dispar, this reduction is quite dramatic, with only six ORs being expressed in the larvae vs. 35 in the female antennae. Results from the PCA analysis indicate that gene expression patterns of female antennal transcriptomes from different subspecies are more similar to one another than to the larval head capsules of their respective subspecies, further supporting strong differences in chemosensory perception between adult and larval stages.
After exploring the amount of shared and exclusively expressed genes for three populations (Figure 1), we observed that AGM and EGM populations share a high number of commonly expressed genes, whereas the JGM population appears to be more divergent, having a high number (14) of uniquely expressed genes. These results suggest that the observed differences are unrelated to flight capacity, indicating either that (a) loss of flight does not impact the olfactory gene repertoire or (b) the secondary loss of flight in the American EGM populations may be too recent to have caused major changes in the genes being expressed.
Interestingly, females from the flightless EGM population display higher gene expression values (RPKM) when compared with JGM and AGM females for most olfaction-related genes except CSPs (Table 1). This could indicate that changes in transcription rates could play an important role in the adaptation of flightless females to their chemical environment. The high variability in olfactory genes and their expression in the larvae suggest that these patterns are unrelated to loss of flight, and we speculate that they are rather adaptations to different ecological pressures.
A detailed comparison of the protein sequences with those reported for other Lepidopteran species through Blast searches and phylogenetic trees supports the common ancestry and high degree of conservation for most olfaction-related gene families within the Lepidoptera, and reveals a high sequence similarity between L. dispar and other members of the Noctuidae clade, particularly for ORs and GRs. A recent study investigating the evolution of these chemoreceptors in the Lepidoptera suggests that the common ancestor of this clade harbored only few OR and GR genes, and that while the number of genes increased greatly during the evolution of the clade, it remained relatively low in comparison to other insect groups. This high degree of conservation possibly occurred because olfaction-related gene expression in the Lepidoptera is under strict regulatory control, limiting the establishment of newly emerged genes (Engsontia et al., 2014).
Although most of our candidate sequences had close alignments with those reported for other model species, a few cases remain where L. dispar sequences were observed to form clusters or single nodes. Further studies are required to confirm the identity of these sequences and establish whether lineage-specific gene expansion occurs in the Lymantriinae clade (including closely related species such as the douglas-fir tussock moth Orgyia pseudotsugata and the nun moth Lymantria monacha) or the superfamily Noctuoidea (including more distantly related species such as Spodoptera spp. and Helicoverpa spp.).
This work represents the most complete description of chemosensory genes and proteins for L. dispar to date. Our results reveal differential gene expression between adult and larval stages characterized by fewer IR, OR, and GR genes being expressed in the larvae, but more CSP genes in comparison to the adults. Comparisons of protein sequences with those from other Lepidopteran species and organisms from different taxa support the common ancestry and high degree of conservation for most olfaction-related gene families. The gene expression patterns in female antennae are more similar to one another than they are to their respective larval stages, whereas larval gene expression patterns are highly divergent across populations. After exploring the number of unique and commonly expressed genes, AGM and EGM populations were found to share a high number of commonly expressed genes, whereas the JGM population appeared to be more divergent. These results indicate that either (a) loss of flight does not impact the olfactory gene repertoire or (b) the secondary loss of flight in American EGM populations may be too recent to cause major changes in the genes being expressed. Nevertheless, higher expression values for GRs, IRs, OBPs, and ORs in EGM females suggest that differences in transcription rates could be an adaptation of flightless females to their chemical environment. Differences in the larval olfactory-related gene expression, on the other hand, are likely responses to unique ecological pressures rather than to female flight ability. Further studies are required to understand the deeper evolutionary and ecological significance of these findings.
AC, EG, CD, MM, and BH designed research; AC and EG collected data; AC, EG, and DW analyzed data; all authors contributed to the writing process. All authors read and approved the submitted version.
This project was possible thanks to the funding provided to AC by the International Postdoctoral Fellowship Program in Plant Science “Plant Fellows,” A Marie Curie initiative managed by the Zurich-Basel Plant Science Center, and funds from ETH Zurich and the Max Planck Society.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We would like to express our gratitude to Hannah Nadel, Supervisory Entomologist of the Animal and Plant Health Inspection Service (APHIS) of the United States Department of Agriculture (USDA) for providing the insects used in this study. We would also like to acknowledge the reviewers for their constructive criticism and valuable suggestions.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fevo.2017.00115/full#supplementary-material
Supplementary Figure 1. Percent of Bombyx mori protein recovery in a blast search for three subspecies of the Gypsy moth (L. dispar). The flight capable AGM (Asian gypsy moth) and JGM (Japanese gypsy moth), and the flightless EGM (European gypsy moth).
APHIS (2013). Gypsy Moth. Available online at: http://www.aphis.usda.gov/plant_health/plant_pest_info/gypsy_moth/ (Accessed 02 Feb 2017).
Anderson, P., and Alborn, H. (1999). Effects on oviposition behaviour and larval development of Spodoptera littoralis by herbivore-induced changes in cotton plants. Entomol. Exp. Appl. 92, 45–51. doi: 10.1046/j.1570-7458.1999.00523.x
Andersson, M. N., Grosse-Wilde, E., Keeling, C. I., Bengtsonn, J. M., Yuen, M. M. S., Li, M., et al. (2013). Antennal transcriptome analysis of the chemosensory gene families in the tree killing bark beetles, Ips typographus and Dendroctonus ponderosae (Coleoptera: Curculionidae: Scolytinae). BMC Genomics 14:198. doi: 10.1186/1471-2164-14-198
Baranchikov, Y. N. (1989). “Ecological basis of the evolution of host relationships in Eurasian gypsy moth populations,” in Proceedings of the Lymantriidae: A Comparison of Features of New and Old World Tussock moths, eds W. E. Wallner, K. A. McManus (Broomall, PA: U.S. Department of Agriculture), 319–338.
Barlow, V. (2004). Gypsy Moth, Lymantria dispar. Avalable online at: http://northernwoodlands.org/articles/article/gypsy_moth_lymantria_dispar (Accessed 02 Feb 2017).
Bell, R. A., Owens, C. D., Shapiro, M., and Tardiff, J. R. (1981). “Mass rearing and virus production,” in The Gypsy Moth: Research Toward Integrated Pest Management, eds C. C. Doane and M. L. McManus (Washington, DC: USDA Forest Service; Technical Bulletin), 599–655.
Bengtsson, J. M., Trona, F., Montagné, N., Anfora, G., Ignell, R., Witzgall, P., et al. (2012). Putative chemosensory receptors of the codling moth, Cydia pomonella, identified by antennal transcriptome analysis. PLoS ONE 7:e31620. doi: 10.1371/journal.pone.0031620
Benton, R., Vannice, K. S., Gomez-Diaz, C., and Vosshall, L. B. (2009). Variant ionotropic glutamate receptors as chemosensory receptors in Drosophila. Cell 136, 149–162. doi: 10.1016/j.cell.2008.12.001
Benton, R., Sachse, S., Michnick, S. W., and Vosshall, L. B. (2006). Atypical membrane topology and heteromeric function of Drosophila odorant receptors in vivo. PLoS Biol. 4:e20. doi: 10.1371/journal.pbio.0040020
Bohbot, J., Pitts, R., Kwon, H. W., Rützler, M., Robertson, H., and Zwiebel, L. (2007). Molecular characterization of the Aedes aegypti odorant receptor gene family. Insect Mol. Biol. 16, 525–537. doi: 10.1111/j.1365-2583.2007.00748.x
Briscoe, A. D., Macias-Mu-oz, A., Kozak, K. M., Walters, J. R., Yuan, F., Jamie, G. A., et al. (2013). Female behaviour drives expression and evolution of gustatory receptors in butterflies. PLoS Genet. 9:e1003620. doi: 10.1371/journal.pgen.1003620
Conesa, A., Götz, S., García-Gómez, J. M., Terol, J., Talón, M., Robles, M., et al. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610
Croset, V., Rytz, R., Cummins, S. F., Budd, A., Brawand, D., Kaessmann, H., et al. (2010). Ancient protostome origin of chemosensory ionotropic glutamate receptors and the evolution of insect taste and olfaction. PLoS Genet. 6:e1001064. doi: 10.1371/journal.pgen.1001064
Engsontia, P., Sangket, U., Chotigeat, W., and Satasook, C. (2014). Molecular evolution of the odorant and gustatory receptor genes in lepidopteran insects: implications for their adaptation and speciation. J. Mol. Evol. 79, 21–39. doi: 10.1007/s00239-014-9633-0
Getahun, M. N., Thoma, M., Lavista-Llanos, S., Keesey, I., Fandino, R. A., Knaden, M., et al. (2016). Intracellular regulation of the insect chemoreceptor complex impacts odour localization in flying insects. J. Exp. Biol. 219, 3428–3438. doi: 10.1242/jeb.143396
Gong, Y. M., and Plettner, E. (2011). Effects of aromatic compounds on antennal responses and on the pheromone-binding proteins of the gypsy moth (Lymantria dispar). Chem. Senses 36, 291–300. doi: 10.1093/chemse/bjq130
Götz, S., García-Gómez, J. M., Terol, J., Williams, T. D., Nagaraj, S. H., Nueda, M. J., et al. (2008). High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 36, 3420–3435. doi: 10.1093/nar/gkn176
Grosse-Wilde, E., Gohl, T., Bouché, E., Breer, H., and Krieger, J. (2007). Candidate pheromone receptors provide the basis for the response of distinct antennal neurons to pheromonal compounds. Eur. J. Neurosci. 25, 2364–2373. doi: 10.1111/j.1460-9568.2007.05512.x
Grosse-Wilde, E., Kuebler, L. S., Bucks, S., Vogel, H., Wicher, D., Hansson, B. S., et al. (2011). Antennal transcriptome of manduca sexta. Proc. Natl. Acad. Sci. U.S.A. 108, 7449–7454. doi: 10.1073/pnas.1017963108
Grosse-Wilde, E., Svatoš, A., and Krieger, J. (2006). A pheromone-binding protein mediates the bombykol-induced activation of a pheromone receptor in vitro. Chem. Senses 31, 547–555. doi: 10.1093/chemse/bjj059
Honson, N., Johnson, M. A., Oliver, J. E., Prestwich, G. D., and Plettner, E. (2003). Structure–activity studies with pheromone-binding proteins of the gypsy moth, Lymantria dispar. Chem. Senses 28, 479–489. doi: 10.1093/chemse/28.6.479
Honson, N. S., and Plettner, E. (2006). Disulfide connectivity and reduction in pheromone-binding proteins of the gypsy moth, Lymantria dispar. Naturwissenschaften 93, 267–277. doi: 10.1007/s00114-006-0096-z
Ioannidis, P., Simao, F. A., Waterhouse, R. M., Manni, M., Seppey, M., Robertson, H. M., et al. (2017). Genomic features of the damselfly Calopteryx splendens representing a sister clade to most insect orders. Genome Biol. Evol. 9, 415–430. doi: 10.1093/gbe/evx006
Jacquin-Joly, E., Legeai, F., Montagné, N., Monsempes, C., François, M. C., Poulain, J., et al. (2012). Candidate chemosensory genes in female antennae of the noctuid moth Spodoptera littoralis. Int. J. Biol. Sci. 8, 1036–1050. doi: 10.7150/ijbs.4469
Katoh, K., Misawa, K., Kuma, K., and Miyata, T. (2002). MAFFT: a novel method for rapid multiple sequence alignment based on fast fourier transform. Nucleic Acids Res. 30, 3059–3066. doi: 10.1093/nar/gkf436
Keena, M. A., Côté, M. J., Grinberg, P. S., and Wallner, W. E. (2008). World distribution of female flight and genetic variation in Lymantria dispar (Lepidoptera: Lymantriidae). Environ. Entomol. 37, 636–649. doi: 10.1603/0046-225X(2008)37[636:WDOFFA]2.0.CO;2
Keena, M. A., Côté, M. J., Grinberg, P. S., and Wallner, W. E. (2010). “Predicting the female flight capability of Gypsy moth by using DNA markers,” in Proceedings 21st USDA Research Forum on Invasive Species, eds K. A. McManus and K. W. Gottschalk (Newtown Square, PA: U.S. Department of Agriculture, Forest Service, Northern Research Station), 35–37.
Keena, M. A., Wallner, W. E., Grinberg, P. S., and Cardé, R. T. (2001). Female flight propensity and capability in Lymantria dispar (Lepidoptera: Lymantriidae) from Russia, North America, and their reciprocal F1 hybrids. Environ. Entomol. 30, 380–387. doi: 10.1603/0046-225X-30.2.380
Koenig, C., Hirsh, A., Bucks, S., Klinner, C., Vogel, H., Shukla, A., et al. (2015). A reference gene set for chemosensory receptor genes of Manduca sexta. Insect Biochem. Mol. Biol. 66, 51–63. doi: 10.1016/j.ibmb.2015.09.007
Kowcun, A., Honson, N., and Plettner, E. (2001). Olfaction in the gypsy moth, Lymantria dispar- effect of Ph, ionic strength, and reductants on pheromone transport by pheromone-binding proteins. J. Biol. Chem. 276, 44770–44776. doi: 10.1074/jbc.M104688200
Krieger, J., Grosse-Wilde, E., Gohl, T., Dewer, Y. M., Raming, K., and Breer, H. (2004). Genes encoding candidate pheromone receptors in a moth (Heliothis virescens). Proc. Natl. Acad. Sci. U.S.A. 101, 11845–11850. doi: 10.1073/pnas.0403052101
Krieger, J., Klink, O., Mohl, C., Raming, K., and Breer, H. (2003). A candidate olfactory receptor subtype highly conserved across different insect orders. J. Comp. Physiol. A 189, 519–526. doi: 10.1007/s00359-003-0427-x
Lartigue, A., Campanacci, V., Roussel, A., Larsson, A. M., Jones, T. A., Tegoni, M., et al. (2002). X-ray structure and ligand binding study of a moth chemosensory protein. J. Biol. Chem. 277, 32094–32098. doi: 10.1074/jbc.M204371200
Leitch, O., Papanicolaou, A., Lennard, C., Kirkbride, K. P., and Anderson, A. (2015). Chemosensory genes identified in the antennal transcriptome of the blowfly Calliphora stygia. BMC Genomics 16:255. doi: 10.1186/s12864-015-1466-8
Lin, W., Yu, Y., Zhou, P., Zhang, J., Dou, L., Hao, Q., et al. (2015). Identification and knockdown of the olfactory receptor (OrCo) in gypsy moth, Lymantria dispar. Int. J. Biol. Sci. 11, 772–780. doi: 10.7150/ijbs.11898
Liu, Y., Gu, S., Zhang, Y., Guo, Y., and Wang, G. (2012). Candidate olfaction genes identified within the Helicoverpa armigera antennal transcriptome. PLoS ONE 7:e48260. doi: 10.1371/journal.pone.0048260
Missbach, C., Vogel, H., Hansson, B. S., and Groβe-Wilde, E. (2015). Identification of odorant binding proteins and chemosensory proteins in antennal transcriptomes of the jumping bristletail lepismachilis y-signata and the firebrat thermobia domestica: evidence for an independent OBP–OR Origin. Chem. Senses 40, 615–626. doi: 10.1093/chemse/bjv.050
Nakagawa, T., Sakurai, T., Nishioka, T., and Touhara, K. (2005). Insect sex-pheromone signals mediated by specific combinations of olfactory receptors. Science 307, 1638–1642. doi: 10.1126/science.1106267
NBII (2011). Global Invasive Species Database: Lymantria dispar. Availabe online at: http://www.issg.org/database/species/ecology.asp?si=96& (Accessed 02 Feb 2017).
Plettner, E., and Gries, R. (2010). Agonists and antagonists of antennal responses of gypsy moth (Lymantria dispar) to the pheromone (+)-disparlure and other odorants. J. Agric. Food Chem. 58, 3708–3719. doi: 10.1021/jf904139e
Pogue, M. G., and Schaeffer, P. W. (2007). A Review of Selected Species of Lymantria (Hubner, 1819). Washington, DC: US Department of Agriculture (USDA). Availabel online at: https://www.fs.fed.us/foresthealth/technology/pdfs/Pogue_Lymantria_web.pdf (Accessed 02 Feb 2017).
Poivet, E., Gallot, A., Montagné, N., Glaser, N., Legeai, F., Jacquin-Joly, E., et al. (2013). A comparison of the olfactory gene repertoires of adults and larvae in the noctuid moth Spodoptera littoralis. PLoS ONE 8:e60263. doi: 10.1371/journal.pone.0060263
Proffit, M., Khallaf, M. A., Carrasco, D., Larsson, M. C., and Anderson, P. (2015). ‘Do you remember the first time?’ Host plant preference in a moth is modulated by experiences during larval feeding and adult mating. Ecol. Lett. 18, 365–374. doi: 10.1111/ele.12419
Rambaut, A. (2007). FigTree, A Graphical Viewer of Phylogenetic Trees. Available online at: http://tree.bio.ed.ac.uk/software/figtree/ (Accessed 02 Feb 2017).
Rambaut, A. (2012). FigTree version 1.4. Available online at: http://tree.bio.ed.ac.uk/software/figtree/ (Accessed 02 Feb 2017).
Rinker, D. C., Zhou, X., Pitts, R. J., Rokas, A., and Zwiebel, L. J. (2013). Antennal transcriptome profiles of anopheline mosquitoes reveal human host olfactory specialization in Anopheles gambiae. BMC Genomics 14:749. doi: 10.1186/1471-2164-14-749
Robertson, H. M., Warr, C. G., and Carlson, J. R. (2003). Molecular evolution of the insect chemoreceptor gene superfamily in Drosophila melanogaster. Proc. Natl. Acad. Sci. U.S.A. 100(Suppl. 2), 14537–14542. doi: 10.1073/pnas.2335847100
Rytz, R., Croset, V., and Benton, R. (2013). Ionotropic receptors (IRs): chemosensory ionotropic glutamate receptors in Drosophila and beyond. Insect Biochem. Mol. Biol. 43, 888–897. doi: 10.1016/j.ibmb.2013.02.007
Sandler, B. H., Nikonova, L., Leal, W. S., and Clardy, J. (2000). Sexual attraction in the silkworm moth: structure of the pheromone-binding-protein–bombykol complex. Chem. Biol. 7, 143–151. doi: 10.1016/S1074-5521(00)00078-8
Sasaki, T., Jikumaru, S., Azuma, W., Kuroda, K., and Ishii, H. (2016). Oviposition site selection by Japanese gypsy moth (Lymatria dispar japonica) in a warm-temperate secondary forest in western Japan. Forest Sci. Technol. 12, 130–136. doi: 10.1080/21580103.2015.1132782
Silbering, A. F., Rytz, R., Grosjean, Y., Abuin, L., Ramdya, P., Jefferis, G. S., et al. (2011). Complementary function and integrated wiring of the evolutionarily distinct Drosophila olfactory subsystems. J. Neurosci. 31, 13357–13375. doi: 10.1523/JNEUROSCI.2360-11.2011
Tanaka, K., Uda, Y., Ono, Y., Nakagawa, T., Suwa, M., Suwa, M., et al. (2009). Highly selective tuning of a silkworm olfactory receptor to a key mulberry leaf volatile. Curr. Biol. 19, 881–890. doi: 10.1016/j.cub.2009.04.035
van Schooten, B., Jiggins, C. D., Briscoe, A. D., and Papa, R. (2016). Genome-wide analysis of ionotropic receptors provides insight into their evolution in Heliconius butterflies. BMC Genomics 17:254. doi: 10.1186/s12864-016-2572-y
Vogt, R. G., Prestwich, G. D., and Lerner, M. R. (1991). Odorant-binding-protein subfamilies associate with distinct classes of olfactory receptor neurons in insects. J. Neurobiol. 22, 74–84. doi: 10.1002/neu.480220108
Walker, W. B. III., Gonzalez, F., Garczynski, S. F., and Witzgall, P. (2016). The chemosensory receptors of codling moth Cydia pomonella–expression in larvae and adults. Sci. Rep. 6:23518. doi: 10.1038/srep23518
Wanner, K., Anderson, A., Trowell, S., Theilmann, D., Robertson, H., Robertson, H. M., et al. (2007). Female-biased expression of odourant receptor genes in the adult antennae of the silkworm, Bombyx mori. Insect Mol. Biol. 16, 107–119. doi: 10.1111/j.1365-2583.2007.00708.x
Wanner, K., and Robertson, H. (2008). The gustatory receptor family in the silkworm moth Bombyx mori is characterized by a large expansion of a single lineage of putative bitter receptors. Insect Mol. Biol. 17, 621–629. doi: 10.1111/j.1365-2583.2008.00836.x
Yu, Y., and Plettner, E. (2013). Enantiomer and conformer recognition of (+) and (-)-disparlure and their analogs by the pheromone binding proteins of the gypsy moth, Lymantria dispar. Bioorg. Med. Chem. 21, 1811–1822. doi: 10.1016/j.bmc.2013.01.043
Zahiri, R., Kitching, I. J., Lafontaine, J. D., Mutanen, M., Kaila, L., Holloway, J. D., et al. (2011). A new molecular phylogeny offers hope for a stable family level classification of the Noctuoidea (Lepidoptera). Zool. Scr. 40, 158–173. doi: 10.1111/j.1463-6409.2010.00459.x
Zhang, J., Wang, B., Dong, S., Cao, D., Dong, J., Walker, W. B., et al. (2015). Antennal transcriptome analysis and comparison of chemosensory gene families in two closely related noctuidae moths, Helicoverpa armigera and H. assulta. PLoS ONE 10:e0117054. doi: 10.1371/journal.pone.0117054
Zhang, S., Zhang, Z., Wang, H., and Kong, X. (2014). Antennal transcriptome analysis and comparison of olfactory genes in two sympatric defoliators, Dendrolimus houi and Dendrolimus kikuchii (Lepidoptera: Lasiocampidae). Insect Biochem. Mol. Biol. 52, 69–81. doi: 10.1016/j.ibmb.2014.06.006
Keywords: Lymantria dispar, transcriptome, odorant receptor, ionotropic receptor, gustatory receptor, odorant binding protein, chemosensory protein
Citation: Clavijo McCormick A, Grosse-Wilde E, Wheeler D, Mescher MC, Hansson BS and De Moraes CM (2017) Comparing the Expression of Olfaction-Related Genes in Gypsy Moth (Lymantria dispar) Adult Females and Larvae from One Flightless and Two Flight-Capable Populations. Front. Ecol. Evol. 5:115. doi: 10.3389/fevo.2017.00115
Received: 09 June 2017; Accepted: 12 September 2017;
Published: 29 September 2017.
Edited by:Stefano Colazza, University of Palermo, Italy
Reviewed by:Julien Pelletier, Keele University, United Kingdom
Nicolas Montagné, Université Pierre et Marie Curie, France
Copyright © 2017 Clavijo McCormick, Grosse-Wilde, Wheeler, Mescher, Hansson and De Moraes. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Consuelo M. De Moraes, email@example.com
†Present Address: Andrea Clavijo McCormick, Massey University, Institute of Agriculture and Environment, Palmerston North, New Zealand