The Evolution and Expression of the Moth Visual Opsin Family

Abstract

Because visual genes likely evolved in response to their ambient photic environment, the dichotomy between closely related nocturnal moths and diurnal butterflies forms an ideal basis for investigating their evolution. To investigate whether the visual genes of moths are associated with nocturnal dim-light environments or not, we cloned long-wavelength (R), blue (B) and ultraviolet (UV) opsin genes from 12 species of wild-captured moths and examined their evolutionary functions. Strong purifying selection appeared to constrain the functions of the genes. Dark-treatment altered the levels of mRNA expression in Helicoverpa armigera such that R and UV opsins were up-regulated after dark-treatment, the latter faster than the former. In contrast, B opsins were not significantly up-regulated. Diel changes of opsin mRNA levels in both wild-captured and lab-reared individuals showed no significant fluctuation within the same group. However, the former group had significantly elevated levels of expression compared with the latter. Consequently, environmental conditions appeared to affect the patterns of expression. These findings and the proportional expression of opsins suggested that moths potentially possessed color vision and the visual system played a more important role in the ecology of moths than previously appreciated. This aspect did not differ much from that of diurnal butterflies.

About the Authors:

Pengjun Xu

Contributed equally to this work with: Pengjun Xu, Bin Lu, Haijun Xiao

Affiliation: State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing, P.R. China

Bin Lu

Contributed equally to this work with: Pengjun Xu, Bin Lu, Haijun Xiao

Affiliation: Department of Herpetology, Chengdu Institute of Biology, Chinese Academy of Sciences, Chengdu, Sichuan, China

Haijun Xiao

Contributed equally to this work with: Pengjun Xu, Bin Lu, Haijun Xiao

Affiliation: State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing, P.R. China

Xiaowei Fu

Affiliation: State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing, P.R. China

Robert W. Murphy

Affiliation: State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, P.R. China

Kongming Wu

* E-mail: [email protected]

Affiliation: State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing, P.R. China

Introduction

The survival of species often involves adaptation to their environment. Vision is one of the most important senses of organisms and selection to the ambient light-environment may drive its evolution. The adaptive evolution between the visual systems of organisms and their ambient light-environment forms a model system for evolutionary research. For example, the opsins and eye morphology of fish are consistent with a possible adaptation to their ambient light environment [1]–[6].

Visual perception and the associated structural components are highly developed in insects. Lepidopterans, both butterflies and moths, depend on visual perception for foraging [7], [8], host plant identification [9], mate-choice [10] and long-distance migration [11]. Butterflies have evolved apposition eyes to better adapt their diurnal lifestyle, moreover, the opsin genes also adaptively evolved in color vision [12]–[15]. Unlike diurnal butterflies, most moths are primarily nocturnal. As with butterflies, they also require vision for feeding and long-distant migration [7], [16]–[18]. The superposition eyes of moths enhance light-capture and this may reflect an adaptation to dim-light vision [12]. Notwithstanding, little is known about how nocturnal lifestyle may have driven the evolution of opsin genes in moths.

Opsins are G-protein-coupled receptors characterized by seven transmembrane domain structures. A lysine residue in the seventh helix of opsins binds a light-sensitive vitamin A-derived chromophore and they determine the spectral sensitivity of the photopigments together [19]. Physiological and phylogenetic analyses indicate that early insects possessed trichromatic vision based on spectral peaks of pigments and subfamilies of opsins: long-wavelength (R) sensitive (>500 nm) proteins, blue (B) sensitive (400–500 nm) proteins and ultraviolet (UV) sensitive (325–400 nm) proteins [20], [21]. Although some insects retain trichromatic vision, many species appear to have adapted their range of spectral sensitivity to their ambient light environment. Gene duplications, losses and sequence variation of opsins occur, for example, in butterflies, mosquitoes and species of Drosophila [14], [15], [21]–[27]. Some members of Coleoptera lost B opsin [28], [29].

The dichotomy between closely related nocturnal moths and diurnal butterflies forms an ideal basis for investigating the evolution of visual genes in response to their ambient photic environment. We investigate the evolution of visual genes associated with nocturnal, dim-light environments by cloning 35 full-length cDNA opsin genes from 12 species of wild-captured moths. We analyze these data in concert with opsins from butterflies to investigate the evolutionary history of opsin and identify selection pressure genes in moths. We also seek to compare the patterns of diel-expression between wild-captured and lab-reared individuals of Helicoverpa armigera.

Materials and Methods

1. Ethics Statement

With permission, we captured the insects in Changdao Experiment Station (Shandong) of Institute of Plant Protection, Chinese Academy of Agricultural Sciences. The wild-captured moths used in this study were serious pests in China, therefore, no permits were required for the described insect collection and experimentation.

2. Insects

Opsin genes were cloned from 12 species of moths, including Agrotis segetum, Agrotis ypsilon, Argyrogramma agnate, Chilo supressalis, H. armigera, Loxostege sticticalis, Macroglossum stellatarum, Mamestra brassicae, Mythimna separata, Plutella xylostella, Spodoptera exigua, Spodoptera litura. Specimens species were captured and stored in liquid N₂ in July 2011 using a vertical pointing trap set up on Beihuang Island, Shandong province. This island, located in the Bohai Gulf (38° 23.200′ N, 120° 54,500′ E), had an area of about 2.5 km² [16]–[18].

We chose H. armigera to investigate the expression patterns of opsins. Diel changes of mRNA levels under natural light in July 2011 were investigated in adult moths captured by light trap every 3 h beginning at 6∶00 (Beijing time). Moths were directly dropped into liquid N₂. A lab-reared population from Langfang (Hebei province, 2005) was maintained at 25°C with a 14∶10, light:dark photoperiod. Different day-instar stages of adults, reared with a diet of 10% sugar and 2% vitamin complex, were collected daily at 9∶00. Diel changes of opsins mRNA level in the lab-reared population were determined from three day old adults (adults from the 3rd day after eclosion). Moths reared in a cage located at a window to sense the gradual change of light through dusk and dawn were collected every 3 h beginning at 6∶00 (Beijing time). The effect of dark-treatment on the expression pattern of opsins was investigated using moths three day old adults maintained in complete darkness and the control individuals reared at the same place while no dark-treatment.

3. Cloning and Sequencing of Opsins

Total RNA was isolated from individual adult moths using TRIzol reagent (Invitrogen, Carlsbad, CA, USA). Single stranded cDNA was synthesized using oligo(dT) and M-MLV Reverse Transcriptase (Promega, Madison, WL, USA). Genomic DNA was extracted from each individual using Easy Pure Genomic DNA Extraction Kit (TransGen, Beijing, China). The qualities of DNA, RNA and cDNA templates were ascertained by using conserved primers listed in Table S1 with the following PCR program: 4 min at 94°C; 30 s at 94°C, 30 s at 55°C, and 30 s at 72°C, 40 cycles; 10 min at 72°C.

Degenerate primers designed according to the known genes of other insects were used to amplify the 5′ends of opsins (Table S1), using the following PCR program: 4 min at 94°C; 30 s at 94°C, 30 s at 55°C, and 60 s at 72°C, 40 cycles; 10 min at 72°C. To obtain complete sequences of opsin, 5′ and 3′ RACE (Rapid Amplification of cDNA Ends) were performed using a SMART RACE cDNA Amplification Kit (Clontech, Palo Alto, CA, USA) with gene-specific primers according to the manufacture’s instructions (Table S1). To ensure the 5′/3′fragments were from the same gene, specific primers containing the full ORFs (open reading frames) were designed according to the 5′/3′ sequences of untranslated regions. In turn, these primers were used to amplify the entire ORF sequences with pfu enzyme (Table S1). The PCR program was as follows: 4 min at 94°C; 30 s at 94°C, 30 s at 55°C, and 150 s at 72°C, 40 cycles; 10 min at 72°C. All PCR products were purified, cloned into the pEASY-T Cloning Vector (TransGen, Beijing, China) and then used for the transformation of Escherichia coli DH5α. The clones were then sequenced.

4. Sequence Alignment, Phylogenetic and Evolutionary Analyses

All nucleotide sequences, including our sequences and others obtained from GenBank (Table S2), were aligned using Clustal W implemented in BioEdit [30] and then manually adjusted. The codon-matched sequences were phylogenetically analyzed using maximum likelihood (ML) and Bayesian inference (BI) with the same outgroup for each gene. ML analyses were performed using RAxML 7.3.2 [31] under the GTR+G substitution model. For estimating nodal support, nonparametric bootstrap proportions [32] with 1000 pseudoreplicates were used. BI analyses were performed using MrBayes 3.1.2 [33], [34]. For each gene, we partitioned the dataset according to codon position in each gene. The multiple partitions controlled for heterogeneity across dataset, such as variation in substitution rates. The best-fitting substitution-model for each partition (Table 1) was selected using JMODELTEST 0.1.1 [35]. The Bayesian information criterion (BIC) was used to select a model because of its high accuracy and precision [36]. For each partition-strategy, Metropolis-coupled Monte Carlo Markov chains were run for 10 million generations with two parallel searches using one cold and three heated chains, each starting with a random tree. Trees were sampled every 1000 generations using split frequencies <0.01 to indicate convergence. TRACER 1.5 [37] was used to determine when the log likelihood (lnL) of sampled trees reached a stationary. In all BI analyses, apparent stationarity was reached within 1 million generations; we conservatively discarded the first 5 million generations from each run as “burn-in” and used the sampled trees from the remaining 5 million generations (5001 trees) to calculate the frequency of nodal resolution in a 50% majority-rule consensus tree, termed Bayesian posterior probabilities (BPPs). Three replicate analyses were conducted to assess whether individual runs failed to converge upon the optimal posterior distribution and if likelihood values, branch lengths, tree topology, and posterior probabilities differed between runs. ML bootstrap proportions (MLBs) ≥70% and BPPs ≥0.95 were considered to indicate strong support for individual nodes [33], [38], [39]. The approach potentially assessed discordance among trees based on the implemented method of inference and provided a preferred topology for selection estimates.

[Figure omitted. See PDF.]

Table 1. Model selection for each codon position of genes in MrBayes analysis.

https://doi.org/10.1371/journal.pone.0078140.t001

5. Branch and Branch-Sites Test of Selection

Maximum Likelihood [40] was employed to test for differences in selection pressure, using the CODEML program of PAML 4.5 [41]. Branch models and branch-site models were employed to detect positive selection acting on the particular lineages. These tree-based tests of selection required the unrooted input-trees. RAxML and Mrbayes trees were converted to unrooted trees by introducing opsin genes of Drosophila_melanogaster. To test the hypothesis, for branch models, five hypotheses were evaluated: 1) one d_n/d_s ratio for all branches; 2) d_n/d_s ratio = 1 for all branches; 3) moth lineage and butterfly lineage have different d_n/d_s ratio (ω₂ and ω₁); 4) neutral evolution for moth (ω₂ = 1); 5) the free-ratio model with free d_n/d_s ratio for each branch. For branch-site modes, the moth lineage was defined as the foreground branch, and remaining lineages were defined as background branches, which were specified in the tree file by using branch labels. Likelihood ratio test (LRT) was used to investigate if the alternative model, indicating positive selection, was superior to the null model.

6. Real-time PCR Analysis

Total RNA was isolated and cDNA was synthesized from whole individual adults. Using β-actin as the reference gene, real-time PCR (qPCR) was carried out with the TaqMan method (FAM) in 20 µl reaction agent comprised of 1 µl of template cDNA, 2*Premix Ex Taq™ (Takara, Japan), 0.2 µM of each primer and 0.4 µM probe (Table S1) on a 7500 Fast Real-time PCR System (Applied Biosystems). Thermal cycling conditions were: 95°C, 2 min and 45 cycles of 95°C for 15 s, 60°C for 34 s. After PCR, a melting curve analysis was performed to demonstrate the specificity of the PCR product, as displayed by a single peak. The cDNA sample of each group was replicated three times. Fold differences of the genes were calculated according to the 2^−ΔΔCT method [42]. Four repeats were used for each data-point.

For comparing proportional changes of mRNA levels of opsins, the absolute quantification methodology using standard curve was performed [43]. Fragments containing the primers and probes of qPCR of opsin genes were amplified with our de novo primers (Table S1) and cloned into the pEASY-T Cloning Vector (TransGen, Beijing, China). These plasmids were used for quantification assay.

Statistical analyses were conducted using STATA v.9.0. Student’s t-test or ANOVA with Bonferroni multiple comparisons were used to determine the level of significance in the relative levels of mRNA expression.

Results

1. Moth Opsins

Highly conserved regions of aligned lepidopteran opsin sequences occurred in the 5′ untranslated region (UTR) of R, B and UV opsins, respectively (figure S1). After locating the ORFs of R, B and UV opsins, pairs of primers were designed based on conserved regions and there were used to amplify the 5′ end fragments of moth opsin genes containing the start codon. Taken together with the 3′ RACE, we successfully obtained complete full-length cDNA of 35 opsins from 12 species of moths. A Blast search at NCBI (http://www.ncbi.nih.gov/BLAST/) revealed that these genes were highly similar to opsins of other insects, especially Lepidoptera. No pseudogenes were found. Alignment of deduced amino acid sequences showed that these genes exhibited typical opsin characteristics, such as seven transmembrane domains, two sites with disulphide cysteine, and a lysine residue that binds the retinal chromophore through a retinylidine Schiff-base linkage in the seventh transmembrane domain (figure S2).

2. Phylogenetic Inference

ML and BI analysis generated quite different topologies depending on the gene and method of tree construction. One possible reason was the absence of sufficient information to resolve the phylogenetic relationships among moths using each gene alone. After concatenating all opsin genes, the topologies from ML and BI still differed. Consequently, both the ML (figure 1) and BI (figure S3) trees for each opsin genes were used for evaluating selection only; this study did not seek to resolve the phylogeny among these moths.

[Figure omitted. See PDF.]

Figure 1. Phylogeny of opsin genes in Lepidoptera based on maximum likelihood.

Values on the nodes are the nonparametric bootstrap proportions (MLBPs). Moths denoted in blue and butterflies in red. Branch-specific ω values are shown on nodes of the common ancestors. a) R opsin; b) B opsin; c) UV opsin. Figures at right show the moths.

https://doi.org/10.1371/journal.pone.0078140.g001

3. Natural Selection Analysis

The free-ratio model provided a better fit to these three opsin genes. Similar conditions occurred among the three opsin genes. For branch model, significant differences between moths and butterflies in d_n/d_s were not detected (Table 2). The branch-specific ω values of moths were 0.0311, 0.0302 and 0.0326 for R, B and UV opsins, respectively. The corresponding ω values of butterflies were higher than those of moths: 0.0797, 0.0522 and 0.0422, respectively (figure 1). Therefore, the hypothesis of strong purifying selection for moths was not rejected. For the branch-site model, the null model also could not be rejected. Similar results were generated for each gene regardless of employment of a ML or BI tree. All analyses indicated that the three types of opsin genes experienced strong purifying selection. The different topologies did not affect the selection-analyses. Therefore, we only showed the selection-analyses results using the ML tree.

[Figure omitted. See PDF.]

Table 2. Selective patterns for opsin genes.

https://doi.org/10.1371/journal.pone.0078140.t002

4. The Expression Level of Opsins mRNA

Lab-reared H. armigera was used to investigate variation among different instar and dark-treated changes in opsin-expression. The results of qPCR indicated that the expression of opsins was dynamic and differed among day-instar-stage of adults but there were no significant difference among different days except for UV in female between day 1 and day 5 (Bonferroni multiple comparisons P = 0.033) (figure 2a). Dark-treatment changed the expression of opsins. The highest expression of UV opsin occurred at 15 min of darkness and then expression was down-regulated to a level similar to untreated individuals. R opsin was significantly up-regulated at 30 min then down-regulated at 1 h and 2 h followed significantly by the highest expression at 4 h. B opsin mRNA was significantly down-regulated after 1 h of darkness, and then up-regulated at 4 h albeit not significantly so (figure 2b).

[Figure omitted. See PDF.]

Figure 2. Relative expression level of R, B and UV opsins at different day instar stage of adults (a) and time of dark-treatment (b).

Mean ± SE. n = 4 for each point. No significant difference occurs among different day instar stage of adults except for UV opsin in females, which is shown using different letters. The “*” and “**” denote statistical significance of the expression levels at P<0.05 and P<0.01, respectively.

https://doi.org/10.1371/journal.pone.0078140.g002

The qPCR analysis did not reveal significant fluctuations of the three opsins within the same group (for wild-captured individuals, B: F = 1.90, d.f. = 7, 24, P = 0.1142; R: F = 1.50, d.f. = 7, 24, P = 0.2143; UV: F = 1.00, d.f. = 7, 21, P = 0.4596; for lab-reared individuals, B: F = 1.78, d.f. = 7, 24, P = 0.1377; R: F = 0.37, d.f. = 7, 24, P = 0.9082; UV: F = 2.32, d.f. = 7, 24, P = 0.0598). The wild-captured moths had significantly higher measures of relative opsin-expression than lab-reared moths at almost all time-points in H. armigera (figure 3, students t test P<0.05 and P<0.01, respectively). The standard curves for each opsin gene were constructed (figure S4) and these served as proportional measures of opsin-expression. The analysis indicated that the expression of UV opsins was proportionally higher in wild-captured than in lab-reared H. armigera (figure 4, figure S5). The proportional change of B opsin was stable in both wild-captured and lab-reared individuals (figure 4).

[Figure omitted. See PDF.]

Figure 3. Relative expression level of R, B and UV opsins between indoor and outdoor Helicoverpa armigera.

Mean ± SE. “Indoor” stands for lab-reared individuals, n = 4 for each point. “Outdoor” stands for the wild-captured individuals, n = 4 for each point except for UV at 3∶00 (n = 3) and 21∶00 (n = 2). The “*” and “**” denote statistical significance of the expression levels at P<0.05 and P<0.01, respectively.

https://doi.org/10.1371/journal.pone.0078140.g003

[Figure omitted. See PDF.]

Figure 4. Proportional (expression relative to total opsin pool) expression of R, B and UV opsins in H. armigera including indoor and outdoor individuals.

Mean ± SE.

https://doi.org/10.1371/journal.pone.0078140.g004

Discussion

1. The Evolution of Moth Opsin

Selection pressure promotes the evolution of sensory systems by adapting organisms to local ecological conditions [44]. The cloning of opsins from several species of moths facilitates the identification of potential evolutionary adaptations of moths to their ambient light environments. The amino acid sequences of the moth pigments show a high percentage of identity with other arthropod opsins and possess a retinal binding site at site 296 (K296), a feature typical of all opsins. In vertebrate visual opsins, E113 acts as a counterion to stabilise the protonated Schiff base, however, in many vertebrate non-visual pigments and invertebrate opsins, the counterion has shifted to site 181 (E181) with a subsequent E113Y substitution [19]. In moths, Y113 is present in blue- and red-sensitive pigments, F113 in UV-sensitive opsins, and E181 in all three opsin classes. Strong signals of purifying selection exist in the three opsin genes. This suggests that these genes have conserved functions in moths. The rod-only retina of deep-sea fishes is sensitive to wavelengths of approximate 480 nm, which matches the ambient dim-light [45]–[48]. Surely, species living in the dim-light environment still need to capture available light. Therefore, opsin genes of species living in dim-light environments should retain functions.

Most species of moths are either nocturnal or crepuscular. Twilight is blue-shifted relative to daylight and may also contain a relatively strong red light [21]. Interestingly, our results show that B and R opsins of moths, which respond to blue and red light, are under strong purifying selection. In moths, these two genes seem to have functional constraints. Unlike moths, the B and R opsins of many diurnal butterflies are highly diverse. This likely reflects adaptations to variable light-niches [14], [22], [49]. However, the loss of B in the darkling beetle (Tribolium castaneum) may be an evolutionary sensory modality tradeoff, most likely driven by the dramatic increase of chemoreceptor rather than a nocturnal lifestyle [28].

UV opsin is the short-wavelength opsin that functions in color vision and strong UV light exposure, which might be harmful to nocturnal insects [50], [51]. Generally, UV opsin might not be so important to nocturnal moths because UV light is weakest at night. However, our results indicate that this gene continues to experience functional constraints. Two non-exclusive observations may explain the function of UV opsin in moths. First, moths have limited diurnal activities though they are far more active at night. In this scenario, UV opsin functions to protect them from UV damage. Second, in addition to offering protection from UV light, UV opsin may still play an important role in color vision, and perhaps other functions. In this case, we do not expect a functional shift in moths that differs from that of diurnal butterflies.

2. The Expression Pattern of Moth Opsin

Adaptation of the visual system usually involves two mechanisms: spectral tuning by amino acid sequence variation and expression-regulation [52]–[55]. The levels of opsin-expression could be a function of visual pigment sensitivity. In this case, regulation is a response to conditions such as light, age and heredity [2], [5], [56]–[59]. Our results show that the expression of opsins in H. armigera also changes after dark-treatment yet uniquely. UV and R opsins are immediately up-regulated, although the latter occurs more slowly than the former, In contrast, dark-treatment does not significantly up-regulate B opsin. Taken with the evidence that strong purifying selection acts on all three kinds of opsin genes, we infer that the opsin genes not only work well for moths, they also facilitate nocturnal activities. Electroretinogram (ERG) analyses of H. armigera suggest that photosensitivity gets stronger with dark-treatment ≤2 h, especially in the UV region [60]. Variation occurs in the significant up-regulation of R and UV opsins but not at each point of detection, yet in other cases significant up-regulation of B opsin is not detected within the first 2 h of dark-treatment. During this time, the screening pigments rapidly moved making the compound eye transform from being light in appearance to becoming black (data not shown). This observation suggests that the superposing optics and movement of screening pigments also play important roles in the nocturnal lifestyle of moths. Some nocturnal butterflies also evolved superposing optics [61] and this observation supports our conjecture.

The expression-patterns of opsin genes correspond significantly with environmental light [2], [4], [5], [59]. Phototaxis has been well-established in moths [62], [63]. However, no studies report diel changes of mRNA levels. Our comparison of wild-captured and lab-reared (>5 yr) H. armigera fails to detect significantly different levels of RNA between day and night within the same group, yet opsins are down-regulated significantly in lab-reared H. armigera. This observation most likely results from the weaker light in the laboratory as opposed to the wild. Consequently, the expression-patterns of moth opsins could vary in different environments. The relative measures of opsin-expression to the housekeeping gene may indicate differences in diel opsin-regulation both within and among species. However, proportional measures of opsin-expression might serve better for making inferences about color vision than relative measure [5]. Diurnal insects have elaborated optic systems that use ambient light for color-vision and navigation [11], [21]. Vision is involved in inter-individual communication. In contrast, nocturnal insects usually communicate using pheromones or songs [64]. Nocturnal insects also have the ability to use light for color vision and flight control [65]–[67]. Our proportional measures identify changes in expression in time in different test populations, suggesting that H. armigera has potentially nocturnal color vision. In particular, the proportion of UV opsin is lower in lab-reared than wild-captured moths, mostly due to small differences in exposure to UV light. This result indicates that UV opsin has an important function in wild moths.

Conclusion

The three kinds of moth opsins are under strong purifying selection. The expression of opsin genes is plastic and varies with the environmental light. Relative and proportional measures of opsin-expression vary with time of day and environment, suggesting that the vision system plays a more important role in the sensory ecology of moths than previously appreciated.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Alignments of 5′ untranslated region (UTR) of B (a), R (b) and UV (c) opsin genes from insects of Lepidoptera. “.” indicates sites identical with the first row sequence. “−” signifies a gap. Highly conserved regions shown with shading.

https://doi.org/10.1371/journal.pone.0078140.s001

(TIF)

Figure S2.

The alignments of amino acid sequences of B (a), UV (b) and R (c) opsins amplified in this study. “.” indicates sites identical to the first row sequence. Transmembrane domains are shaded.  = the site of the chromophore Schiff-base linkage (K). • = the counterion (Glu) (E). ▾ = two disulphide Cys (C) residues. ▪ =  Position 90 in bovine rodopsin, the residues for blue and UV. ○ =  Position 113 in bovine rodopsin. AA =  Argyrogramma agnata, AS =  Agrotis segetum, AY =  Agrotis ypsilon, CS =  Chilo supressalis, HA =  Helicoverpa armigera, LS =  Loxostege sticticalis, MB =  Mamestra brassicae, MS =  Mythimna separata, MSt =  Macroglossum stellatarum, PX =  Plutella xylostella, SL =  Spodoptera litura, SE =  Spodoptera exigua.

https://doi.org/10.1371/journal.pone.0078140.s002

(TIF)

Figure S3.

Phylogenetic reconstruction of opsin genes in Lepidopteran based on Bayesian inference. Values on the nodes are the Bayesian posterior probabilities (BPPs). Moths denoted in blue and butterflies in red. a) R opsin; b) B opsin; c) UV opsin.

https://doi.org/10.1371/journal.pone.0078140.s003

(TIF)

Figure S4.

The standard curves for opsins of H. armigera determined by triplicate sampling. The primer efficiency is 103.7%, 98.7% and 100.7% for R, B and UV opsins, respectively.

https://doi.org/10.1371/journal.pone.0078140.s004

(TIF)

Figure S5.

Proportional (expression relative to total opsin pool) expression of UV opsins in H. armigera between indoor and outdoor individuals. Mean ± SE. The “*” and “**” denote statistical significance of the expression levels at P<0.05 and P<0.01, respectively.

https://doi.org/10.1371/journal.pone.0078140.s005

(TIF)

Table S1.

Primers used in this study. The abbreviation of species name was shown in bracket.

https://doi.org/10.1371/journal.pone.0078140.s006

(PDF)

Table S2.

The GenBank accession number of genes used in this study.

https://doi.org/10.1371/journal.pone.0078140.s007

(PDF)

Author Contributions

Conceived and designed the experiments: PJX HJX KMW. Performed the experiments: PJX HJX XWF. Analyzed the data: PJX BL HJX KMW. Contributed reagents/materials/analysis tools: PJX HJX KMW. Wrote the paper: PJX BL HJX RWM KMW.

Citation: Xu P, Lu B, Xiao H, Fu X, Murphy RW, Wu K (2013) The Evolution and Expression of the Moth Visual Opsin Family. PLoS ONE 8(10): e78140. https://doi.org/10.1371/journal.pone.0078140

References

1. Yokoyama S, Zhang H, Radlwimmer FB, Blow NS (1999) Adaptive evolution of color vision of the Comoran coelacanth (Latimeria chalumnae). Proc Natl Acad Sci U S A 96: 6279–6284.

2. Hofmann CM, O’Quin KE, Smith AR, Carleton KL (2010) Plasticity of opsin gene expression in cichlids from Lake Malawi. Mol Ecol 19: 2064–2074.

3. O’Quin KE, Hofmann CM, Hofmann HA, Carleton KL (2010) Parallel evolution of opsin gene expression in African cichlid fishes. Mol Biol Evol 27: 2839–2854.

4. Tobler M, Coleman SW, Perkins BD, Rosenthal GG (2010) Reduced opsin gene expression in a cave-dwelling fish. Biol Lett 6: 98–101.

5. Fuller RC, Claricoates KM (2011) Rapid light-induced shifts in opsin expression: finding new opsins, discerning mechanisms of change, and implications for visual sensitivity. Mol Ecol 20: 3321–3335.

6. Schmitz L, Wainwright PC (2011) Nocturnality constrains morphological and functional diversity in the eyes of reef fishes. BMC Evol Biol 11: 338.

7. Cutler DE, Bennett RR, Stevenson RD, White RH (1995) Feeding behavior in the nocturnal moth Manduca sexta is mediated mainly by violet receptors, but where are they located in the retina? J Exp Biol 198: 1909–1917.

8. Weiss MR (2001) Vision and learning in some neglected pollinators: beetles, flies, moths and butterflies. Cambridge, Cambridge University Press.

9. Kelber A (1999) Ovipositing butterflies use a red receptor to see green. J Exp Biol 202: 2619–2630.

10. Jiggins CD, Naisbit RE, Coe RL, Mallet J (2001) Reproductive isolation caused by colour pattern mimicry. Nature 411: 302–305.

11. Sauman I, Briscoe AD, Zhu H, Shi D, Froy O, et al. (2005) Connecting the navigational clock to sun compass input in monarch butterfly brain. Neuron 46: 457–467.

12. Stavenga DG, Arikawa K (2006) Evolution of color and vision of butterflies. Arthropod Struct Dev 35: 307–318.

13. Frentiu FD, Bernard GD, Cuevas CI, Sison-Mangus MP, Prudic KL, et al. (2007a) Adaptive evolution of color vision as seen through the eyes of butterflies. Proc Natl Acad Sci U S A 104: 8634–8640.

14. Frentiu FD, Bernard GD, Sison-Mangus MP, Brower AV, Briscoe AD (2007b) Gene duplication is an evolutionary mechanism for expanding spectral diversity in the long-wavelength photopigments of butterflies. Mol Biol Evol 24: 2016–2028.

15. Briscoe AD, Bybee SM, Bernard GD, Yuan F, Sison-Mangus MP, et al. (2010) Positive selection of a duplicated UV-sensitive visual pigment coincides with wing pigment evolution in Heliconius butterflies. Proc. Natl. Acad. Sci. U S A 107: 3628–3633.

16. Feng HQ, Fu XW, Wu B, Wu KM (2009) Seasonal migration of Helicoverpa armigera (Lepidoptera: Noctuidae) over the Bohai Sea. J Econ Entomol 102: 95–104.

17. Feng HQ, Wu KM, Cheng DF, Guo YY (2004) Northward migration of Helicoverpa armigera (Lepidoptera: Noctuidae) and other moths in early summer observed with radar in northern China. J Econ Entomol 97: 1874–1883.

18. Feng HQ, Wu KM, Ni YX, Cheng DF, Guo YY (2005) Return migration of Helicoverpa armigera (Lepidoptera: Noctuidae) during autumn in northern China. Bull Entomol Res 95: 361–370.

19. Terakita A (2005) The opsins. Genome Biol 6: 213.

20. Chittka L (1996) Dose bee color vision predate the evolution of flower color? Naturwossenschaften 83: 136–138.

21. Briscoe AD, Chittka L (2001) The evolution of color vision in insects. Annu Rev Entomol 46: 471–510.

22. Briscoe AD (2000) Six opsins from the butterfly Papilio glaucus: molecular phylogenetic evidence for paralogous origins of red-sensitive visual pigments in insects. J Mol Evol 51: 110–121.

23. Spaethe J, Briscoe AD (2004) Early duplication and functional diversification of the opsin gene family in insects. Mol Biol Evol 21: 1583–1594.

24. Sison-Mangus MP, Bernard GD, Lampel J, Briscoe AD (2006) Beauty in the eye of the beholder: the two blue opsins of lycaenid butterflies and the opsin gene-driven evolution of sexually dimorphic eyes. J Exp Biol 209: 3079–3090.

25. Sison-Mangus MP, Briscoe AD, Zaccardi G, Knuttel H, Kelber A (2008) The lycaenid butterfly Polyommatus icarus uses a duplicated blue opsin to see green. J Exp Biol 211: 361–369.

26. Bernard GD, Stavenga DG (1979) Spectral sensitivities of retinular cells measured in intact living flies by an optical method. J Comp Physiol A 134: 95–107.

27. Salcedo E, Huber A, Henrich S, Chadwell LV, Chou WH, et al. (1999) Blueand green-absorbing visual pigments of Drosophila: ectopic expression and physiological characterization of the R8 photoreceptor cell-specific Rh5 and Rh6 rhodopsins. J Neurosci 19: 10716–10726.

28. Jackowska M, Bao R, Liu Z, McDonald EC, Cook TA, et al. (2007) Genomic and gene regulatory signatures of cryptozoic adaptation: Loss of blue sensitive photoreceptors through expansion of long wavelength-opsin expression in the red flour beetle Tribolium castaneum. Front Zool 4 24.

29. Oba Y, Kainuma T (2009) Diel changes in the expression of long wavelength-sensitive and ultraviolet-sensitive opsin genes in the Japanese firefly, Luciola cruciata. Gene 436: 66–70.

30. Thompson JD, Higgins DG (1994) Gibson (1994) CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucl Acids Res 22: 4673–4680.

31. Stamatakis A (2006) RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22: 2688–2690.

32. Felsenstein J (1985) Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39: 783–791.

33. Huelsenbeck JP, Rannala B (2004) Frequentist properties of Bayesian posterior probabilities of phylogenetic trees under simple and complex substitution models. Syst Biol 53: 904–913.

34. Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574.

35. Posada D (2008) jModelTest: phylogenetic model averaging. Mol Biol Evol 25: 1253–1256.

36. Luo A, Qiao H, Zhang Y, Shi W, Ho SY, et al. (2010) Performance of criteria for selecting evolutionary models in phylogenetics: a comprehensive study based on simulated datasets. BMC Evol Biol 10: 242.

37. Drummond AJ, Rambaut A (2007) Tracer v1.4, available from http://beast.bio.ed.ac.uk/Tracer.

38. Hillis DM, Bull JJ (1993) An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analysis. Syst Biol 42: 182–192.

39. Erixon P, Svennblad B, Britton T, Oxelman B (2003) Reliability of Bayesian posterior probabilities and bootstrap frequencies in phylogenetics. Syst Biol 52: 665–673.

40. Nielsen R, Yang Z (1998) Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics 148: 929–936.

41. Yang Z (2007) PAML 4: a program package for phylogenetic analysis by maximum likelihood. Mol Biol Evol 24: 1586–1591.

42. Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 25: 402–408.

43. Wong ML, Medrano JF (2005) Real-time PCR for mRNA quantitation. Biotechniques 39: 75–85.

44. Endler JA, Basolo AL (1998) Sensory ecology, receiver biases and sexual selection. Trends Ecol Evol 13: 415–420.

45. Hope AJ, Partridge JC, Dulai KS, Hunt DM (1997) Mechanisms of wavelength tuning in the rod opsins of deep-sea fishes. Proc R Soc Lond B 264: 155–163.

46. Douglas RH, Partridge JC, Marshall NJ (1998) The eyes of deep-sea fish. I: Lens pigmentation, tapeta and visual pigments. Prog Retin Eye Res 17: 597–636.

47. Wagner HJ, Frohlich E, Negishi K, Collin SP (1998) The eyes of deep-sea fish. II. Functional morphology of the retina. Prog Retin Eye Res 17: 637–685.

48. Hunt DM, Dulai KS, Partridge JC, Cottrill P, Bowmaker JK (2001) The molecular basis for spectral tuning of rod visual pigments in deep-sea fish. J Exp Biol 204: 3333–3344.

49. Awata H, Wakakuwa M, Arikawa K (2009) Evolutiaon of color vision in pierid butterflies: blue opsin duplication, ommatidial heterogeneity and eye regionalization in Colias erate. J Comp Physiol A Neuroethol Sens Neural Behav Physiol 195: 401–408.

50. Koyanagi M, Nagata T, Katoh K, Yamashita S, Tokunaga F (2008) Molecular evolution of arthropod color vision deduced from multiple opsin genes of jumping spiders. J Mol Evol 66: 130–137.

51. Kashiyama K, Seki T, Numata H, Goto SG (2009) Molecular characterization of visual pigments in Branchiopoda and the evolution of opsins in Arthropoda. Mol Biol Evol 26: 299–311.

52. Yokoyama S, Shi Y (2000) Genetics and evolution of ultraviolet vision in vertebrates. FEBS Lett 486: 167–172.

53. Fuller RC, Carleton KL, Fadool JM, Spady TC, Travis J (2004) Population variation in opsin expression in the bluefin killifish, Lucania goodei: a real-time PCR study. J Comp Physiol A Neuroethol Sens Neural Behav Physiol 190: 147–154.

54. Parry JW, Carleton KL, Spady T, Carboo A, Hunt DM, et al. (2005) Mix and match color vision: tuning spectral sensitivity by differential opsin gene expression in Lake Malawi cichlids. Curr Biol 15: 1734–1739.

55. Carleton K (2009) Cichlid fish visual systems: mechanisms of spectral tuning. Integr Zool 4: 75–86.

56. Sasagawa H, Narita R, Kitagawa Y, Kadowaki T (2003) The expression of genes encoding visual components is regulated by a circadian clock, light environment and age in the honeybee (Apis mellifera). Eur J Neurosci 17: 963–970.

57. Hofmann CM, Carleton KL (2009) Gene duplication and differential gene expression play an important role in the diversification of visual pigments in fish. Integr Comp Biol 49: 630–643.

58. Larmuseau MH, Vancampenhout K, Raeymaekers JA, Van Houdt JK, Volckaert FA (2010) Differential modes of selection on the rhodopsin gene in coastal Baltic and North Sea populations of the sand goby, Pomatoschistus minutus.. Mol Ecol 19: 2256–2268.

59. Smith AR, D’Annunzio L, Smith AE, Sharma A, Hofmann CM, et al. (2010) Intraspecific cone opsin expression variation in the cichlids of Lake Malawi. Mol Ecol 20: 299–310.

60. Wei GS, Zhang QW, Zhou M, Wu WG (2002) Characteristic response of the compound eyes of Helicoverpa armigera to light. Acta Entomologica Sinica 45: 323–328 (In chinese)..

61. Yack JE, Johnson SE, Brown SG, Warrant EJ (2007) The eyes of Macrosoma sp. (Lepidoptera: Hedyloidea): a nocturnal butterfly with superposition optics. Arthropod Struct Dev 36: 11–22.

62. Danthanarayana W (1976) Diel and lunar flight periodicities in the light brown apple moth, Epiphyas postvittana (Walker) (Tortricidae) and their possible adaptive significance. Aust J Zool 24: 65–73.

63. Frank DK (1988) Impact of outdoor lighting on moths: an assessment. J Lepidopt Soc 42: 63–93.

64. Lewis T (1984) The elements and frontiers of insect communication. London: Academic Press.

65. Kelber A, Balkenius A, Warrant EJ (2003) Colour vision in diurnal and nocturnal hawkmoths. Integr Comp Biol 43: 571–579.

66. Baird E, Kreiss E, Wcislo W, Warrant E, Dacke M (2011) Nocturnal insects use optic flow for flight control. Biol Lett 7: 499–501.

67. Warrant E, Dacke M (2010) Vision and visual navigation in nocturnal insects. Annu Rev Entomol 56: 239–254.