Stress Responses of Small Heat Shock Protein Genes in Lepidoptera Point to Limited Conservation of Function across Phylogeny

Bo Zhang; Jincheng Zheng; Yu Peng; Xiaoxia Liu; Ary A. Hoffmann; Chun-Sen Ma

doi:10.1371/journal.pone.0132700

Abstract

The small heat shock protein (sHsp) family is thought to play an important role in protein refolding and signal transduction, and thereby protect organisms from stress. However little is known about sHsp function and conservation across phylogenies. In the current study, we provide a comprehensive assessment of small Hsp genes and their stress responses in the oriental fruit moth (OFM), Grapholita molesta. Fourteen small heat shock proteins of OFM clustered with related Hsps in other Lepidoptera despite a high level of variability among them, and in contrast to the highly conserved Hsp11.1. The only known lepidopteran sHsp ortholog (Hsp21.3) was consistently unaffected under thermal stress in Lepidoptera where it has been characterized. However the phylogenetic position of the sHsps within the Lepidoptera was not associated with conservation of induction patterns under thermal extremes or diapause. These findings suggest that the sHsps have evolved rapidly to develop new functions within the Lepidoptera.

Citation: Zhang B, Zheng J, Peng Y, Liu X, Hoffmann AA, Ma C-S (2015) Stress Responses of Small Heat Shock Protein Genes in Lepidoptera Point to Limited Conservation of Function across Phylogeny. PLoS ONE 10(7): e0132700. https://doi.org/10.1371/journal.pone.0132700

Editor: Eric A. Shelden, Washington State University, UNITED STATES

Received: December 8, 2014; Accepted: June 17, 2015; Published: July 21, 2015

Copyright: © 2015 Zhang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: This research was supported by the National Special Fund for Agro-scientific Research in the Public Interest (No. 201103024) and Beijing Natural Science Foundation (5144032). CSM and BZ are the authors who received the fundings respectively. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Heat shock protein (HSP) families have been widely associated with physiological processes and stress resistance. The small heat shock proteins (sHSPs) (12–43 kDa) represent a particularly heterogeneous [1] family of HSPs, in contrast with other HSPs such as HSP70 and HSP90 that are highly conserved [2,3]. Many sHsp genes are upregulated under thermal stress and may protect insects against stressful conditions; for example Hsp26 is thought to protect Aedes aegypti larvae and pupae from heat stress [4], while Hsp20 may underlie a tradeoff between thermal protection at the cost of reproductive impairment in Liriomyza huidobrensis [5]. The sHsp genes may also play a role in cold tolerance, as suggested by the upregulation of Hsp20 expression in the recovery of Drosophila melanogaster from cold injury [6]. Hsp23 is strongly upregulated during pupal diapause in the flesh fly, Sarcophaga crassipalpis [7,8], while Hsp20.8 expression has been linked to diapause in Sesamia nonagrioides [9]. The sHsps therefore represent an important inducible chaperone family involved in stress responses.

The conservation of sHsps has rarely been considered in insects, in contrast to Hsp70 which shows a high degree of evolutionary conservation both in its sequence and in its induction pattern [10,11]. The functional significance of the sHps in many insect orders including the Lepidoptera remains unclear, despite an increasing amount of sequence information becoming available. Lepidopteran sHsps are thought to have evolved independently from sHsps in other orders, except that one ortholog, Hsp21.4-like protein, has been discovered in insect species from multiple orders [12–14]. The Hsp10 family also appears to be highly conserved across insects, and also in vertebrates where it acts as a co-chaperone with Hsp60 [15]. Some studies have considered sHsp expression patterns in Lepidoptera in an attempt to understand their potential function. Changes in expression patterns suggest that the sHsp genes in the armyworm, Spodoptera litura, may have roles in cold responses as well as moth metamorphosis [14,16]. In the codling moth, Cydia pomonella, three sHps were upregulated under heat stress and were persistently expressed during development [17], suggesting a role in both processes.

The Oriental fruit moth (OFM), Grapholita molesta, one of the world’s most invasive orchard pest insect of stone and pome fruits, causes severe economic loss in global fruit industry [18–20]. It is thought to have originated in China and spread to East Asia, Europe, Africa, Australia, and America during last century [20]. OFM has gained recognition as a global pest of peaches and other fruit crops since its introduction into America [21]. Its generation time depends on temperature and latitude, and varies considerably across regions; for example one or two generations occur in northeastern Asia, compared to six to seven annually in Georgia, USA [22,23]. To overwinter in cold environments, OFM diapause at a late larvae stage, which is induced by low temperature and a short photoperiod [18]. The developmental threshold of OFM ranges from 4°C to 11°C at the lower end, and from 30°C to 35°C at the upper end, suggesting a wide thermal performance breadth which has likely assisted their worldwide expansion [19].

Here we describe the sHsps genes in OFM and their expression patterns across developmental stages, across tissues and in response to heat stress, cold stress and diapause. We use these data to address a number of questions. (1) How do the sHsp and Hsp10 genes of OFM compare to those of other Lepidoptera? Do they fall into specific phylogenetic groupings? (2) Are there unique or common patterns of expression changes in the sHsps when they are under the different stresses or at different temporal/spatial stages? (3) Has the putative functional role of particular sHsps as defined by their patterns of induction been conserved across different lepidopteran species, such that related sHsps group together?

Materials and Methods

Insect preparation

We collected OFM larvae from an experimental station in Hebei province (N39°47′, E118°41′), where the pests could be sampled without permission. The field studies did not involve endangered or protected species. OFM were reared on fresh apples in the laboratory over two years at 24°C and 60% relative humidity under a photoperiod of 15: 9 (L: D). Six developmental stages of OFM, 3-day-old eggs, 1^st instar larvae, 3^rd instar larvae, 5^th instar larvae, pupae, and adults (equal numbers of males and females), were collected and frozen in liquid nitrogen before being stored at -80°C for further RNA extraction. The 5th instar larvae were dissected on ice, and each tissue including head, thorax, abdomen, fat body, and midgut were kept in RNAlater (Qiagen, Hilden, Germany) at -80°C.

To characterize expression patterns of Hsps in adults under thermal extremes, high and low temperatures were initially defined based on pilot experiments; these were the most severe conditions to which recently-eclosed adults could be exposed in the absence of mortality. Moths that had eclosed within 24 hours were therefore kept at -5°C or 40°C for 1 and 2 hours, and then left to recover at 24°C for 2 hours before being frozen in liquid nitrogen. Controls were treated identically but were not stressed. The thermal treatment is based on Drosophila Hsp studies [6,24] and a pilot experiment, in which different recovery and stress periods were considered. Those sHsps that were responsive under 2 h were also highly responsive under shorter or longer exposure times, while non-responsive sHsps failed to be expressed under all conditions tested. For the diapause treatment, diapause was induced by rearing moths at same temperature of 24°C under a photoperiod of 12 hours [25]. Each treatment consisted of five biological replicates, with 14–20 individuals (7–10 per sex) per replicate.

RNA extraction and cDNA synthesis

Total RNA was isolated with an RNeasy Mini kit (Qiagen, Valencia, CA, USA) and digested with DNase I (Invitrogen) following supplier’s instructions. The concentration of RNA was measured using a NanoDrop spectrophotometer (ND-1000, Thermo Fisher Scientific, Wilmington, DE). 2 μg of total RNA was used as a template in first strand cDNA synthesis using M-MLV reverse transcriptase following the manufacturer’s instruction with modified oligo (dT) (Promega, Madison, WI, USA). The final 25 μl volume cDNA templates were used for Real-time quantitative PCR of gene expression. The single strand cDNA from the whole body of OFM fifth instar larvae serving as the RACE template was synthesized by 1 μg of total RNA following instructions in the SMART RACE cDNA Amplification Kit (Clontech, Mountain View, CA, USA).

Clone of full length OFM sHsps

The conserved α-crystalline domain of lepidopteran sHsps was used to design the degenerated primers for full length sHsps in OFM. 5’- and 3’ RACE was done separately by universal primers and gene specific primers with the following cycles: five cycles consisting of 94°C for 5 s, 72°C for 3 min, then five cycles consisting of 94°C for 5 s, 70°C for 10s and 72°C for 3 min, and then 25 cycles consisting of 94°C for 5 s, 68°C for 10 s and 72°C for 180 s. A final extension was performed under 72°C for 10 min. The PCR products were visualized by electrophoresis in a 1.5% agarose gel and purified using an ABgene Ultra PCR Clean-Up Kit (Thermo Scientific) before they were cloned into a pGEM-Teasy vector (Promega, Madison, WI, USA). The products were sequenced on a 3730xl DNA Analyser (Applied Biosystems).

Real-time quantitative PCR (RT-qPCR)

The sHSP RT qPCR primers were designed based on the full length cDNA transcripts (S1 Table). All cDNA templates were reversely transcribed from 2 μg total RNA in the final 25 μl volume, and 1 μl was used in the RT-qPCR reaction. β-actin and GAPDH were both chosen as the reference genes for normalizing the mRNA expression level, however data were also normalized without using reference genes through the Norma-gene normalization algorithm [26] which involves normalization based on the experimental data rather than reference genes to avoid the issue of error in reference gene expression, and the assumption that the reference genes are not influenced by the different conditions tested. The sHsp transcript standards were produced by cloning the PCR products into plasmids. Seven 10-fold serial dilutions of stock plasmid served as quantitative standards to estimate the relative expression of each sHsp gene transcript. RT-qPCR was conducted in a 20 μl volume comprising of 10 μl 2×SYBR Green PCR Master Mix (Takara, Dalian, China), 1μl cDNA template/plasmid, 1μl RT-qPCR primers (10mM working solution), and 8 μl dd H₂O. A Stratagene Mx3000P thermal cycler was used. The PCR cycling parameters were as follows: initial denaturation of 10 s at 95°C, and 40 cycles of 5 s at 95°C, 20 s at 55°C, and 20 s at 72°C. Each experimental group included five biological replicates, each of which contained three technical repeats.

Differences between treatments for normalized expression were compared either by t-test (for comparison of two means), or by one-way analysis of variance (ANOVA) followed by a Tukey B test for posthoc comparisons among means, run with SPSS 18.0 software. Treatment differences were considered significant at p < 0.01. Values are denoted as means ± SE (standard error). Unsupervised hierarchical clustering was performed with Cluster v3.0 software using uncentered Pearson correlations and complete linkage, and run with Java TreeView software [27,28].

Phylogenetic analysis

The NCBI BLAST tool was used to detect sHsps and ensure the fragments sequenced involved sHsp conserved domains. The open reading frames (ORFs), translations, and predicted secondary structure of sHsps were obtained online (http://www.bioinformatics.org/sms2). The molecular weights of sHsps were also estimated by an online tool (http://web.expasy.org/protparam). For comparisons of sHsps across diverse groups, the sHsp gene sequences derived from other insect species were downloaded from GenBank, (S2 Table) and aligned in MEGA6 [29]. All phylogenetic analyses were run using Bayesian Inference (BI), Maximum Likelihood (ML) and Neighbor Joining (NJ) methods. BI analysis was performed using MrBayes, v 3.1.2 [30]. Two sets of four chains were allowed to run simultaneously for 1,000,000 generations. Each set was sampled every 100 generations with a burn-in of 25%. Stationarity was considered to have been reached when the average standard deviation of split frequencies was less than 0.01. Bayesian posterior probabilities (BPP) were estimated based on a 50% majority rule consensus tree of the remaining trees. ML analysis was conducted using an online program (http://www.atgc-montpellier.fr/phyml/) with a GTR substitution model. Bootstrap values were obtained for the nodes based on 1000 replicates. NJ trees were built in MEGA5 based on p-distance and bootstrap values for nodes were based on 1000 replicates.

Phylogenetic reconstruction and signal detection

To examine the possible evolutionary conservation of sHsp functions in response to temperature and diapause, the expression traits were analyzed to trace phylogenetic signal based on a tree generated through BI. We collected data on temperature responses of lepidopteran species from published references, and characterized sHsp expression into three categories (upregulated, no change, or unknown) under the three conditions considered (heat, cold and diapause). Downregulated patterns were not considered because these were absent from our study and the published literature. We then evaluated the number of evolutionary steps between genes with similar functions (same pattern of expression under a condition) based on parsimony reconstruction, and compared these with the number of steps expected by chance under a null model. Expression characteristics were reshuffled 1000 times across the tips of the phylogeny to obtain confidence intervals and test for significance under p = 0.01. These comparisons were run in Mesquite 1.12 [31].

Results

Cloning and characterization of GmHsps

The full-length Hsps were obtained by using 5’ and 3’ rapid amplification of cDNA ends from the cDNA library of the whole body of OFM. Fourteen GmHsp genes were sequenced and deposited in GenBank (KP843895- KP843908). Based on their predicted molecular weight, they were named as GmHsp11.1, 18.9, 19.6, 19.8a, 19.8b, 19.9, 20.4, 21.3, 21.4, 21.7, 22.1, 22.5, 24.8, and 31.8. The ORFs spanned between 312 and 849 bp, coding from 103 to 282 amino acid proteins.

GmHsps were classified into two families, 13 of them being sHsps and one being Hsp10 (also named as chaperonin 10, Cpn10) (Figs 1 and 2). A conserved α-crystalline domain comprising about 100 amino acids was detected in all OFM sHsps, containing 9 β-sheet sandwich structures numbered β2 to β10 (Fig 1). The typical 12 putative dimer interfaces were found to have highly conserved sites of the metazoan α-crystalline domain (ACD), while GmHsp22.1 only presented eight putative dimer interfaces. While GmHsp20.4 was classified as a sHsp, no dimer interface was found, as its conserved domain was ACD sHsp p23-like. Due to the high divergence of the sHsp sequences, genetic distance comparisons were based on the 105 bp conserved ACD characterized by NCBI. GmHsp20.4 showed higher average genetic distances based on the number of amino acid substitutions per site between sequences, ranging from 1.881 substitutions with GmHsp24.8 to 2.708 with GmHsp19.8a. While C and N terminals showed significant differentiation between the full-length gene of GmHsp21.4 and GmHsp21.7, their distance between conserved ACD domains was the lowest (0.00). The only GmHsp11.1 in OFM showed typical Cpn10 superfamily characteristics (NCBI–Conserved domains accession: cl09113), with two roof hairpins, 13 putative oligomerisation interfaces, and one mobile loop (Fig 2).

Download:

Fig 1. 13 sHsp gene alignment of G. molesta.

The predicted secondary structure showing the α-crystallin domain (ACD) is underlined by a blue solid line. Yellow arrows indicate β - strands. Red boxes indicate putative dimer interfaces. The green highlighted fragment is the ACD p23-like superfamily.

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

Download:

Fig 2. Alignment of Hsp11.1 of G. molesta with thirteen other insect species.

The predicted secondary structure of the Cpn10 superfamily is underlined by a blue solid line. Red arrows show two roof hairpins. Yellow boxes indicate a mobile loop, and green dots indicate 13 putative oligomerisation interfaces.

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

Phylogenetic analysis and alignment of GmHsps

To analyze the relationships of OFM sHsps to those of other insects, 71 lepidopteran sHsps including the 13 GmsHsps were collected, and about 100 bp of ACD were truncated to construct Bayesian phylogenetic trees (Fig 3A), which were similar with other trees obtained by different approaches (S1 Fig). The sHsps from lepidopteran species presented two clusters, one with a metazoan ACD and the other one with an ACD p23-like domain, though only with 50% and 73% supporting values respectively. The high level of variation among metazoan ACDs was also evident among lepidopteran sHsps. Three GmHsps showed a close relationship with those from the related species Cydia pomonella, from which only three sHsps have been identified so far [17]. The phylogenetic tree also showed that not all the sHsps from the same species clustered in the same branch, as in the case of Bombyx mori and Danaus plexippus. The orthologous clusters were prevalent in lepidopteran sHsps, including GmHsp19.6-like, GmHsp21.3-like, and GmHsp31.8-like clusters. An unrooted neighbor-joining tree of 13 sHsps was constructed to show the relationship within the sHsp family in this species (Fig 3B). The conserved domain (ACD) of sHsps was clustered with high supporting values, while GmHsp20.4 containing the p23-like domain was isolated. This analysis also indicated that the sHsps with a similar molecular weight did not follow phylogeny, suggesting that the molecular weights of sHsps depended on the variable N- and C-terminal extensions rather than overall similarity.

Download:

Fig 3. Phylogenetic trees of OFM sHsps.

A. Phylogenetic tree of the lepidopteran sHsps. A total of 71 lepidopteran sHsps and two outgroup sequences were used in this Bayesian phylogenetic tree. Only the conserved a-crystalline domains were aligned and used for tree construction. Percentage bootstrap values above 50% were indicated on each cluster. B. An unrooted neighbor-joining tree was constructed with 13 sHsp amino acid sequences from G. molesta. Bootstrap supports of > 40% are shown on the nodes of the tree. The G. molesta sHsps are labeled in bold. Thin arrows indicated the close relationship between codling moth and OFM.

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

Compared with the complex relationship and high level of divergence of OFM sHsps, GmHsp11.1 of the Cpn10 family exhibited a high level of conservation among insects from different orders (Fig 4). Three lepidopteran Cpn10 genes were clustered within a branch with 100% support. Within the Diptera, the suborders Nematocera (mosquitos) and Brachycera (fruit flies) were separated.

Download:

Fig 4. Phylogeny of Cpn10 amino acid sequences of G. molesta and 14 other insect species.

The tree was constructed using the Bayesian method. Bootstrap support above 50% is indicated.

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

Expression in response to thermal treatments and diapause

Expression of most GmHsp genes was significantly upregulated under thermal stress, with average fold changes relative to controls ranging from 5-fold to 270-fold. In response to heat and cold stress, eight genes (GmHsp19.6, 19.8b, 19.9, 21.4, 21.7, 24.8, 31.8, 11.1) were consistently upregulated, particularly GmHsp19.6, whose expression significantly increased with treatment duration (Tables 1 and 2). Ten GmHsp genes were upregulated under heat shock compared with eight genes under the cold treatment. Three genes (GmHsp19.8a, GmHsp20.4 and GmHsp21.3) did not change in expression during the thermal treatments. The hierarchical cluster analysis for all differentially expressed genes revealed changing patterns of expression for most genes in response to heat and cold stress (S2 Fig). Expression changes of 14 GmHsps in diapausing and nondiapausing OFM are shown in Fig 5. GmHsp 21.3 exhibited the highest fold change in diapause individuals which was highly significant (t = -35.74, df = 8, p < 0.001), and together with five other genes showed upregulated expression during diapause.

Download:

Fig 5. Mean normalized expression (with standard errors) of fourteen GmHsps in diapausing G. molesta.

The values were based on log2 transformation of RT-qPCR of the assayed Hsps. Asterisks indicate mean values that are significantly (P < 0.01) different in 5^th instar larvae between diapausing and control groups. A value of 0 indicates no difference in expression, whereas positive and negative values indicate upregulation and downregulation respectively.

https://doi.org/10.1371/journal.pone.0132700.g005

Download:

Table 1. Comparison of the overall expression of GmHsp genes under control conditions, after exposure to heat for one hour (H1h), or two hours (H2h) followed by a two-hour recovery period in each case.

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

Download:

Table 2. Comparison of the overall expression of GmHsp genes under control conditions, after exposure to cold for one hour (C1h) or two hours (C2h) followed by two-hour recovery period in each case.

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

All GmHsps were detected in each dissected tissue and developmental stage but expressed at different levels (S3 and S4 Figs). GmHsp19.9 and GmHsp21.3 exhibited higher expression levels in the head, with a 4.26 and 4.27 fold change in expression respectively. GmHsp31.8 was the only gene with relatively higher (2.71-fold) expression in the midgut. GmHsp11.1 exhibited high expression during development except at the final larval instar stage and in pupae. GmHsp21.3 showed a 24-fold increase at the 5^th instar stage as development proceeded. In contrast, GmHsp20.4 and GmHsp22.1 showed relatively low expression (1-fold and 2-fold on average respectively) across all developmental stages.

Low phylogenetic signal in lepidopteran sHsp expression patterns

The lepidopteran sHsp genes were tested for phylogenetic signal of expression patterns in Mesquite (S3 Table). Expression responses to the three stresses of heat, cold, and diapause showed no evidence of being conserved within phylogenetic lineages (S5 Fig). This was based on the number of steps expected for each character when compared to the number of steps expected from randomly generated trees. A character was considered to possess phylogenetic signal if there were fewer steps than in 99% of the randomly generated trees. There were four evolutionary steps between two traits in the heat response character, falling out of the range of the null model (mean step = 8.4, range 6–9). For expression changes under cold and diapause conditions there were 15 and 14 evolutionary steps respectively, which did not differ significantly from the 18.1 (range 15–21) and 13.1 (range 11–14) steps predicted under a null model. These findings suggest that there was no detectable phylogenetic signature in functional patterns of the sHsps based on the criteria given.

Discussion

In the present study, fourteen new members from the sHsp and Hsp10 families have been identified from OFM. The expression levels of nearly half of these sHsp genes changed significantly in response to thermal stresses. These sHsp are likely to be important Hsp members involved in development and stress tolerance in OFM. However, no phylogenetic signal was found in sHsp responses to stress or diapause, suggesting an inconsistency between similarities in sequence, structure and function of the lepidopteran sHsps. This may reflect rapid evolution of sHsp genes towards new functions under various environments.

High variation of sHsp in sequence feature

All thirteen sHsps in OFM can be separated into two different branches, ACD and p23-like sHsps (Figs 1 and 3). Among them, 12 sHsps were found to contain characteristic ACD, which shows the same number of paralogues of typical ACD in the D. melanogaster genome [32]. One orthologous group was reported in Lepidoptera [13], and GmHsp21.3 is a member of this common ortholog in OFM. It is very different from sHsps of humans in which all 10 orthologous sHsps can be identified in other mammals [33]. Indeed, most of the sHSPs of Lepidoptera show extensive sequence variation because of the N-terminal arm of divergent sequence and variable length and a C-terminal extension [34,35]. This high variation of sHsp sequences contrasts markedly to conserved patterns for the other Hsp families, including Hsp70, Hsp90, and Hsp10 [36]. GmHsp11.1 was identified as a member of Hsp10 which is highly conserved among insects. The conservation of Hsp10 is as high as for Hsp90 (41–76% vs. 25–88%) [37,38].

The potential functions of sHsp in temporal and spatial expression

Gene functions are often deduced from known gene families with similar conserved sequence structure [39,40]. Four out of 12 sHsp members of Drosophila have been well-characterized for their expression patterns and functions [41], but this information cannot necessarily be used to deduce function for sHsp genes in other insects. The function of Hsp21.3 also has not yet been established even though it is commonly found in Lepidoptera. We thus inferred the potential functions of sHsps in OFM from temporal and spatial patterns of expression in our study and from previous research suggesting that various sHsp may be involved in insect development as well as stress responses. A developmental function for Hsps has been suggested by tissue-specific and tightly regulated patterns of expression during development in Drosophila [42–44]. The expression profiles of GmHsp11.1 may point to a putative function in reproduction, embryo development and maternal effects, because of its high expression at the adult and egg stage, in contrast to a steep decrease at the late larval stage. Hps10 is known as a chaperone in mitochondria where it is involved with cell proliferation and differentiation [45]. The increased sHsps (GmHsp18.9, GmHsp19.8a, and GmHsp31.8) from the muscles/cuticle of the abdomen might link to metamorphosis.

GmHsp21.3 may also play a role in development given its increased expression from the neo-larval to pupal stages. The 5^th instar stage is key for OFM entering pupal metamorphosis or diapause. High expression of three GmHsp genes (GmHsp19.9, 21.3, and 21.7) in the head of 5^th instars may point to a possible function in signal reception and transmission. Many genes are potentially involved with diapause [46,47], and in our study GmHsp21.3 is the most markedly upregulated together with five other GmHsps during early diapause. Hsp23 has previously been shown upregulated in diapausing flesh fly [7,8], while ArHsp21 was shown to be upregulated in diapausing embryos of Artemia franciscana [48]. The upregulation of sHsp genes during diapause of OFM may point to a protection function in cold environments, contributing to cold-hardiness of overwintering insects in combination with other Hsps such as Hsp70 [49]. However, there was no common feature of sHsps that were upregulated in diapausing insects.

Our expression results suggest that sHsps may be important in sub-lethal temperature tolerance in OFM, acting as molecular chaperones under extreme conditions. In OFM, sub-lethal temperatures frequently arise in different seasons. For example, low temperatures at night often occur in early spring or late autumn when OFM adults might regulate expression of genes such as GmHsp19.6 and GmHsp21.7. Hot days in summer that are increasing in frequency under climate warming might induce expression of genes like GmHsp19.9 and GmHsp24.8, particularly in adults that cannot easily find a refuge unlike larvae inside plant tissue.

Lack of phylogenetic signal of sHsps expression for function divergence

We analyzed sequence and expression data available for sHsps to test for potential functional conservation (heat, cold and diapause) among lepidopteran insects. We found that similarities in sequence were inconsistent with functional conservation in sHsp genes, as assessed by the association between stress responses and the expression of sHsps across a lepidopteran phylogeny covering five families. Our findings do not fit in well with the notion that conserved domains can be used to extrapolate gene function from model species to target species under investigation. This may reflect the fact that sHsps with the unified structural ACD do not display conserved functional binding sites affecting expression patterns. Those sHsps with relatively constant expression patterns may have “housekeeping” roles in maintaining molecular structures at certain stages [50] rather than acting as conventional stress chaperones.

The inconsistent sHsp response patterns suggest that sHsp responses are evolutionarily labile. While sHsps are induced at high temperatures in many organisms [51], not all described sHsps show increased expression when organisms are placed under temperature stress [12,14,52]. Given the large number of sHsps in organisms, it is possible that only a subset is required for chaperone functions to maintain cellular viability. This may leave sHsps free to evolve additional functions apart from needing to act as stress chaperons and protective agents [42]. The absence of strong phylogenetic signal points to functional constraints being restricted to closely related groups rather than higher taxonomic levels, even though sHsps might not necessarily always diverge after speciation, as in the case of the ortholog hsp21.3.

Genome sequence data is continuing to expand knowledge about the structure of sHSPs [34,35], but the most sHsp families in insects remain poorly characterized. Thus, more information is needed to further assess the presence/absence of phylogenetic signal and links to functional roles. This involves testing more species with divergent traits (particularly for cold and diapause responses) and also characterizing expression patterns in more detail. Signal might be present but not evident from whole organism comparisons; for instance lepidopteran species might exhibit temporal and/or spatial variation in expression. In addition, the classification of sHsp by molecular weight and conserved domains should be further characterized, and functional tests need to also focus on manipulations of gene expression such as through RNAi [53].

Supporting Information

S1 Fig. Phylogenetic trees of OFM sHsps based on the conserved a-crystalline domains (about 105 amino acids) for the tree construction.

A. The NJ tree with bootstrap support > 40% on the nodes. B. The ML tree with bootstrap support > 40% on the nodes. The G. molesta sHsps are labeled in bold. Thin arrows indicated the close relationship between codling moth and OFM.

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

(DOCX)

S2 Fig. Hierarchical clustering of differentially expressed genes in response to heat and cold treatments (1h and 2h, respectively).

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

(DOCX)

S3 Fig. Normalized mRNA expression levels of the fourteen GmHsps in different tissues.

Error bars are standard errors.

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

(DOCX)

S4 Fig. Normalized mRNA expression of the 14 GmHsp genes in OFM at different stages.

Error bars are standard errors.

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

(DOCX)

S5 Fig. Phylogenetic tree of the lepidopteran species considered, including the state of change in the branches and tips in three characters involving heat (A), cold (B), and diapause (C) responses.

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

(DOCX)

S1 Table. The GmsHsp real-time PCR primers, annealing temperatures, and expected fragment lengths.

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

(DOCX)

S2 Table. List of sHsps of Lepidopteran species used in phylogenetic analysis.

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

(DOCX)

S3 Table. List of sHsp expression information used in lepidopteran phylogenetic signal detection.

https://doi.org/10.1371/journal.pone.0132700.s008

(DOCX)

Acknowledgments

We thank Dr. Gang Ma, Dr. Fei Zhao and Dr. Wei Zhang for discussions. We also appreciate valuable comments for manuscript improvement from the editor and two anonymous reviewers.

Author Contributions

Conceived and designed the experiments: BZ CSM. Performed the experiments: BZ JZ YP. Analyzed the data: BZ JZ YP AAH CSM. Contributed reagents/materials/analysis tools: XL AAH CSM. Wrote the paper: BZ JZ YP XL AAH CSM.

References

1. Taylor RP, Benjamin IJ (2005) Small heat shock proteins: a new classification scheme in mammals. J Mol Cell Cardiol 38: 433–444. pmid:15733903
- View Article
- PubMed/NCBI
- Google Scholar
2. Brocchieri L, Conway de Macario E, Macario AJ (2008) Hsp70 genes in the human genome: Conservation and differentiation patterns predict a wide array of overlapping and specialized functions. BMC Evol Biol 8: 19. pmid:18215318
- View Article
- PubMed/NCBI
- Google Scholar
3. Johnson JL (2012) Evolution and function of diverse Hsp90 homologs and cochaperone proteins. Biochim Biophys Acta 1823: 607–613. pmid:22008467
- View Article
- PubMed/NCBI
- Google Scholar
4. Zhao LM, Becnel JJ, Clark GG, Linthicum KJ (2010) Expression of AeaHsp26 and AeaHsp83 in Aedes aegypti (Diptera: Culicidae) larvae and pupae in response to heat shock stress. J Medi Entomol 47: 367–375.
- View Article
- Google Scholar
5. Huang LH, Kang L (2007) Cloning and interspecific altered expression of heat shock protein genes in two leafminer species in response to thermal stress. Insect Mol Biol 16: 491–500. pmid:17651238
- View Article
- PubMed/NCBI
- Google Scholar
6. Colinet H, Lee SF, Hoffmann A (2010) Temporal expression of heat shock genes during cold stress and recovery from chill coma in adult Drosophila melanogaster. FEBS J 277: 174–185. pmid:19968716
- View Article
- PubMed/NCBI
- Google Scholar
7. Rinehart JP, Denlinger DL (2000) Heat-shock protein 90 is down-regulated during pupal diapause in the flesh fly, Sarcophaga crassipalpis, but remains responsive to thermal stress. Insect Mol Biol 9: 641–645. pmid:11122473
- View Article
- PubMed/NCBI
- Google Scholar
8. Hayward SA, Pavlides SC, Tammariello SP, Rinehart JP, Denlinger DL (2005) Temporal expression patterns of diapause-associated genes in flesh fly pupae from the onset of diapause through post-diapause quiescence. J Insect Physiol 51: 631–640. pmid:15993127
- View Article
- PubMed/NCBI
- Google Scholar
9. Gkouvitsas T, Kontogiannatos D, Kourti A (2008) Differential expression of two small Hsps during diapause in the corn stalk borer Sesamia nonagrioides (Lef.). J Insect Physiol 54: 1503–1510. pmid:18804477
- View Article
- PubMed/NCBI
- Google Scholar
10. Mukhopadhyay I, Saxena DK, Chowdhuri DK (2003) Hazardous effects of effluent from the chrome plating industry: 70 kDa heat shock protein expression as a marker of cellular damage in transgenic Drosophila melanogaster (hsp70-lacZ). Environ Health Perspect 111: 1926–1932. pmid:14644668
- View Article
- PubMed/NCBI
- Google Scholar
11. Mayer MP, Bukau B (2005) Hsp70 chaperones: Cellular functions and molecular mechanism. Cell Mol Life Sci 62: 670–684. pmid:15770419
- View Article
- PubMed/NCBI
- Google Scholar
12. Sakano D, Li B, Xia Q, Yamamoto K, Banno Y, Fujii H, et al. (2006) Genes encoding small heat shock proteins of the silkworm, Bombyx mori. Biosci Biotechnol Biochem 70: 2443–2450. pmid:17031055
- View Article
- PubMed/NCBI
- Google Scholar
13. Li ZW, Li X, Yu QY, Xiang ZH, Kishino H, Zhang Z (2009) The small heat shock protein (sHSP) genes in the silkworm, Bombyx mori, and comparative analysis with other insect sHSP genes. BMC Evol Biol 9: 215. pmid:19715580
- View Article
- PubMed/NCBI
- Google Scholar
14. Shen Y, Gu J, Huang LH, Zheng SC, Liu L, Xu WH, et al. (2011) Cloning and expression analysis of six small heat shock protein genes in the common cutworm, Spodoptera litura. J Insect Physiol 57: 908–914. pmid:21510953
- View Article
- PubMed/NCBI
- Google Scholar
15. Tazir Y, Steisslinger V, Soblik H, Younis AE, Beckmann S, Grevelding CG, et al. (2009) Molecular and functional characterisation of the heat shock protein 10 of Strongyloides ratti. Mol Biochem Parasit 168: 149–157.
- View Article
- Google Scholar
16. Gu J, Huang LX, Shen Y, Huang LH, Feng QL (2012) Hsp70 and small Hsps are the major heat shock protein members involved in midgut metamorphosis in the common cutworm, Spodoptera litura. Insect Mol Biol 21: 535–543. pmid:22957810
- View Article
- PubMed/NCBI
- Google Scholar
17. Garczynski SF, Unruh TR, Guedot C, Neven LG (2011) Characterization of three transcripts encoding small heat shock proteins expressed in the codling moth, Cydia pomonella (Lepidoptera: Tortricidae). Insect Sci 18: 473–483.
- View Article
- Google Scholar
18. Russell DA (1986) Ecology of the oriental fruit moth, Grapholita molesta (Busck) in New Zealand: University of Auckland. PhD thesis.
19. Jones MM (2010) Susceptibility of Oriental fruit moth, (Grapholita molesta (Busck)) to selected insecticides and mixtures: University of Illinois at Urbana-Champaign. PhD thesis.
20. Kirk H, Dorn S, Mazzi D (2013) Worldwide population genetic structure of the oriental fruit moth (Grapholita molesta), a globally invasive pest. BMC Ecol 13: 12. pmid:23531126
- View Article
- PubMed/NCBI
- Google Scholar
21. Jones MM, Robertson JL, Weinzierl RA (2010) Susceptibility of oriental fruit moth (Lepidoptera: Tortricidae) larvae to selected reduced-risk insecticides. J Econ Entomol 103: 1815–1820. pmid:21061985
- View Article
- PubMed/NCBI
- Google Scholar
22. Rothschild G, Vickers R (1991) Biology, ecology and control of the oriental fruit moth; In Tortricid pests: their biology, natural enemies and control. Edited by van der Geest LPS, Evenhuis HH. Amsterdam, The Netherlands: Elsevier.
23. Tao W, Sheng C, Chen R, Zhang Y, Zhang S, Chang T, et al. (2011) Investigation on generations of Grapholita molesta in Changchun region by sex pheromen traps and studies on related techniques. J Jilin Agri Sci 36: 43–45.
- View Article
- Google Scholar
24. Colinet H, Hoffmann AA (2012) Comparing phenotypic effects and molecular correlates of developmental, gradual and rapid cold acclimation responses in Drosophila melanogaster. Funct Ecol 26: 84–93.
- View Article
- Google Scholar
25. Dickson R (1949) Factors govering the induction of diapause in the oriental fruit moth. Ann Entomol Soc Am 42: 511–537.
- View Article
- Google Scholar
26. Heckmann LH, Sorensen PB, Krogh PH, Sorensen JG (2011) NORMA-Gene: A simple and robust method for qPCR normalization based on target gene data. Bmc Bioinformatics 12: 250. pmid:21693017
- View Article
- PubMed/NCBI
- Google Scholar
27. Eisen MB, Spellman PT, Brown PO, Botstein D (1998) Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci Unit States Am 95: 14863–14868.
- View Article
- Google Scholar
28. Saldanha AJ (2004) Java Treeview—extensible visualization of microarray data. Bioinformatics 20: 3246–3248. pmid:15180930
- View Article
- PubMed/NCBI
- Google Scholar
29. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S, et al. (2011) MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 28: 2731–2739. pmid:21546353
- View Article
- PubMed/NCBI
- Google Scholar
30. Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574. pmid:12912839
- View Article
- PubMed/NCBI
- Google Scholar
31. Maddison WP, Maddison DR (2006) Mesquite: A modular system for evolutionary analysis. Ver. 1.1. Available via http://www.mesquiteproject.org.
32. Michaud S, Morrow G, Marchand J, Tanguay RM (2002) Drosophila small heat shock proteins: cell and organelle-specific chaperones? Prog Mol Subcell Biol 28: 79–101. pmid:11908067
- View Article
- PubMed/NCBI
- Google Scholar
33. Kampinga HH, Hageman J, Vos MJ, Kubota H, Tanguay RM, Bruford EA, et al. (2009) Guidelines for the nomenclature of the human heat shock proteins. Cell Stress Chaperon 14: 105–111.
- View Article
- Google Scholar
34. Poulain P, Gelly JC, Flatters D (2010) Detection and architecture of small heat shock protein monomers. PLoS One 5: e9990. pmid:20383329
- View Article
- PubMed/NCBI
- Google Scholar
35. Kriehuber T, Rattei T, Weinmaier T, Bepperling A, Haslbeck M, Buchner J (2010) Independent evolution of the core domain and its flanking sequences in small heat shock proteins. FASEB J 24: 3633–3642. pmid:20501794
- View Article
- PubMed/NCBI
- Google Scholar
36. Huang LH, Wang HS, Kang L (2008) Different evolutionary lineages of large and small heat shock proteins in eukaryotes. Cell Res 18: 1074–1076. pmid:18766171
- View Article
- PubMed/NCBI
- Google Scholar
37. Zhu JY, Zhao N, Yang B (2012) Global transcriptome profiling of the pine shoot beetle, Tomicus yunnanensis (Coleoptera: Scolytinae). PLoS One 7: e32291. pmid:22384206
- View Article
- PubMed/NCBI
- Google Scholar
38. Hull JJ, Geib SM, Fabrick JA, Brent CS (2013) Sequencing and de novo assembly of the western tarnished plant bug (Lygus hesperus) transcriptome. PLoS One 8: e55105. pmid:23357950
- View Article
- PubMed/NCBI
- Google Scholar
39. Levine MT, McCoy C, Vermaak D, Lee YC, Hiatt MA, Matsen FA, et al. (2012) Phylogenomic analysis reveals dynamic evolutionary history of the Drosophila heterochromatin protein 1 (HP1) gene family. PLoS Genet 8: e1002729. pmid:22737079
- View Article
- PubMed/NCBI
- Google Scholar
40. Mi H, Muruganujan A, Thomas PD (2013) PANTHER in 2013: modeling the evolution of gene function, and other gene attributes, in the context of phylogenetic trees. Nucleic Acids Res 41: D377–386. pmid:23193289
- View Article
- PubMed/NCBI
- Google Scholar
41. Michaud S, Tanguay RM (2003) Expression of the Hsp23 chaperone during Drosophila embryogenesis: association to distinct neural and glial lineages. BMC Dev Biol 3: 9. pmid:14617383
- View Article
- PubMed/NCBI
- Google Scholar
42. Basha E, O'Neill H, Vierling E (2012) Small heat shock proteins and alpha-crystallins: dynamic proteins with flexible functions. Trends Biochem Sci 37: 106–117. pmid:22177323
- View Article
- PubMed/NCBI
- Google Scholar
43. Morrow G, Tanguay RM (2012) Small heat shock protein expression and functions during development. Int J Biochem Cell Biol 44: 1613–1621. pmid:22502646
- View Article
- PubMed/NCBI
- Google Scholar
44. Wheeler JC, Bieschke ET, Tower J (1995) Muscle-specific expression of Drosophila hsp70 in response to aging and oxidative stress. Proc Natl Acad Sci Unit States Am 92: 10408–10412.
- View Article
- Google Scholar
45. Jia H, Halilou AI, Hu L, Cai W, Liu J, Huang B (2011) Heat shock protein 10 (Hsp10) in immune-related diseases: one coin, two sides. Int J Biochem Mol Biol 2: 47–57. pmid:21969171
- View Article
- PubMed/NCBI
- Google Scholar
46. Denlinger DL (2002) Regulation of diapause. Annu Rev Entomol 47: 93–122. pmid:11729070
- View Article
- PubMed/NCBI
- Google Scholar
47. Zhang Q, Denlinger DL (2010) Molecular characterization of heat shock protein 90, 70 and 70 cognate cDNAs and their expression patterns during thermal stress and pupal diapause in the corn earworm. J Insect Physiol 56: 138–150. pmid:19782689
- View Article
- PubMed/NCBI
- Google Scholar
48. Qiu Z, Macrae TH (2008) ArHsp21, a developmentally regulated small heat-shock protein synthesized in diapausing embryos of Artemia franciscana. Biochem J 411: 605–611. pmid:18095938
- View Article
- PubMed/NCBI
- Google Scholar
49. Rinehart JP, Li A, Yocum GD, Robich RM, Hayward SA, Denlinger DL (2007) Up-regulation of heat shock proteins is essential for cold survival during insect diapause. Proc Natl Acad Sci Unit States Am 104: 11130–11137.
- View Article
- Google Scholar
50. Mymrikov EV, Seit-Nebi AS, Gusev NB (2011) Large potentials of small heat shock proteins. Physiol Rev 91: 1123–1159. pmid:22013208
- View Article
- PubMed/NCBI
- Google Scholar
51. Jakob U, Gaestel M, Engel K, Buchner J (1993) Small heat shock proteins are molecular chaperones. J Biol Chem 268: 1517–1520. pmid:8093612
- View Article
- PubMed/NCBI
- Google Scholar
52. Lu MX, Hua J, Cui YD, Du YZ (2014) Five small heat shock protein genes from Chilo suppressalis: characteristics of gene, genomic organization, structural analysis, and transcription profiles. Cell Stress Chaperon 19: 91–104.
- View Article
- Google Scholar
53. Bousquet F, Nojima T, Houot B, Chauvel I, Chaudy S, Dupas S, et al (2012) Expression of a desaturase gene, desat1, in neural and nonneural tissues separately affects perception and emission of sex pheromones in Drosophila. Proc Natl Acad Sci Unit States Am 109: 249–254.
- View Article
- Google Scholar

[ref1] 1. Taylor RP, Benjamin IJ (2005) Small heat shock proteins: a new classification scheme in mammals. J Mol Cell Cardiol 38: 433–444. pmid:15733903
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Brocchieri L, Conway de Macario E, Macario AJ (2008) Hsp70 genes in the human genome: Conservation and differentiation patterns predict a wide array of overlapping and specialized functions. BMC Evol Biol 8: 19. pmid:18215318
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Johnson JL (2012) Evolution and function of diverse Hsp90 homologs and cochaperone proteins. Biochim Biophys Acta 1823: 607–613. pmid:22008467
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. Zhao LM, Becnel JJ, Clark GG, Linthicum KJ (2010) Expression of AeaHsp26 and AeaHsp83 in Aedes aegypti (Diptera: Culicidae) larvae and pupae in response to heat shock stress. J Medi Entomol 47: 367–375.
View Article
Google Scholar

[14] View Article

[15] Google Scholar

[ref5] 5. Huang LH, Kang L (2007) Cloning and interspecific altered expression of heat shock protein genes in two leafminer species in response to thermal stress. Insect Mol Biol 16: 491–500. pmid:17651238
View Article
PubMed/NCBI
Google Scholar

[17] View Article

[18] PubMed/NCBI

[19] Google Scholar

[ref6] 6. Colinet H, Lee SF, Hoffmann A (2010) Temporal expression of heat shock genes during cold stress and recovery from chill coma in adult Drosophila melanogaster. FEBS J 277: 174–185. pmid:19968716
View Article
PubMed/NCBI
Google Scholar

[21] View Article

[22] PubMed/NCBI

[23] Google Scholar

[ref7] 7. Rinehart JP, Denlinger DL (2000) Heat-shock protein 90 is down-regulated during pupal diapause in the flesh fly, Sarcophaga crassipalpis, but remains responsive to thermal stress. Insect Mol Biol 9: 641–645. pmid:11122473
View Article
PubMed/NCBI
Google Scholar

[25] View Article

[26] PubMed/NCBI

[27] Google Scholar

[ref8] 8. Hayward SA, Pavlides SC, Tammariello SP, Rinehart JP, Denlinger DL (2005) Temporal expression patterns of diapause-associated genes in flesh fly pupae from the onset of diapause through post-diapause quiescence. J Insect Physiol 51: 631–640. pmid:15993127
View Article
PubMed/NCBI
Google Scholar

[29] View Article

[30] PubMed/NCBI

[31] Google Scholar

[ref9] 9. Gkouvitsas T, Kontogiannatos D, Kourti A (2008) Differential expression of two small Hsps during diapause in the corn stalk borer Sesamia nonagrioides (Lef.). J Insect Physiol 54: 1503–1510. pmid:18804477
View Article
PubMed/NCBI
Google Scholar

[33] View Article

[34] PubMed/NCBI

[35] Google Scholar

[ref10] 10. Mukhopadhyay I, Saxena DK, Chowdhuri DK (2003) Hazardous effects of effluent from the chrome plating industry: 70 kDa heat shock protein expression as a marker of cellular damage in transgenic Drosophila melanogaster (hsp70-lacZ). Environ Health Perspect 111: 1926–1932. pmid:14644668
View Article
PubMed/NCBI
Google Scholar

[37] View Article

[38] PubMed/NCBI

[39] Google Scholar

[ref11] 11. Mayer MP, Bukau B (2005) Hsp70 chaperones: Cellular functions and molecular mechanism. Cell Mol Life Sci 62: 670–684. pmid:15770419
View Article
PubMed/NCBI
Google Scholar

[41] View Article

[42] PubMed/NCBI

[43] Google Scholar

[ref12] 12. Sakano D, Li B, Xia Q, Yamamoto K, Banno Y, Fujii H, et al. (2006) Genes encoding small heat shock proteins of the silkworm, Bombyx mori. Biosci Biotechnol Biochem 70: 2443–2450. pmid:17031055
View Article
PubMed/NCBI
Google Scholar

[45] View Article

[46] PubMed/NCBI

[47] Google Scholar

[ref13] 13. Li ZW, Li X, Yu QY, Xiang ZH, Kishino H, Zhang Z (2009) The small heat shock protein (sHSP) genes in the silkworm, Bombyx mori, and comparative analysis with other insect sHSP genes. BMC Evol Biol 9: 215. pmid:19715580
View Article
PubMed/NCBI
Google Scholar

[49] View Article

[50] PubMed/NCBI

[51] Google Scholar

[ref14] 14. Shen Y, Gu J, Huang LH, Zheng SC, Liu L, Xu WH, et al. (2011) Cloning and expression analysis of six small heat shock protein genes in the common cutworm, Spodoptera litura. J Insect Physiol 57: 908–914. pmid:21510953
View Article
PubMed/NCBI
Google Scholar

[53] View Article

[54] PubMed/NCBI

[55] Google Scholar

[ref15] 15. Tazir Y, Steisslinger V, Soblik H, Younis AE, Beckmann S, Grevelding CG, et al. (2009) Molecular and functional characterisation of the heat shock protein 10 of Strongyloides ratti. Mol Biochem Parasit 168: 149–157.
View Article
Google Scholar

[57] View Article

[58] Google Scholar

[ref16] 16. Gu J, Huang LX, Shen Y, Huang LH, Feng QL (2012) Hsp70 and small Hsps are the major heat shock protein members involved in midgut metamorphosis in the common cutworm, Spodoptera litura. Insect Mol Biol 21: 535–543. pmid:22957810
View Article
PubMed/NCBI
Google Scholar

[60] View Article

[61] PubMed/NCBI

[62] Google Scholar

[ref17] 17. Garczynski SF, Unruh TR, Guedot C, Neven LG (2011) Characterization of three transcripts encoding small heat shock proteins expressed in the codling moth, Cydia pomonella (Lepidoptera: Tortricidae). Insect Sci 18: 473–483.
View Article
Google Scholar

[64] View Article

[65] Google Scholar

[ref18] 18. Russell DA (1986) Ecology of the oriental fruit moth, Grapholita molesta (Busck) in New Zealand: University of Auckland. PhD thesis.

[ref19] 19. Jones MM (2010) Susceptibility of Oriental fruit moth, (Grapholita molesta (Busck)) to selected insecticides and mixtures: University of Illinois at Urbana-Champaign. PhD thesis.

[ref20] 20. Kirk H, Dorn S, Mazzi D (2013) Worldwide population genetic structure of the oriental fruit moth (Grapholita molesta), a globally invasive pest. BMC Ecol 13: 12. pmid:23531126
View Article
PubMed/NCBI
Google Scholar

[69] View Article

[70] PubMed/NCBI

[71] Google Scholar

[ref21] 21. Jones MM, Robertson JL, Weinzierl RA (2010) Susceptibility of oriental fruit moth (Lepidoptera: Tortricidae) larvae to selected reduced-risk insecticides. J Econ Entomol 103: 1815–1820. pmid:21061985
View Article
PubMed/NCBI
Google Scholar

[73] View Article

[74] PubMed/NCBI

[75] Google Scholar

[ref22] 22. Rothschild G, Vickers R (1991) Biology, ecology and control of the oriental fruit moth; In Tortricid pests: their biology, natural enemies and control. Edited by van der Geest LPS, Evenhuis HH. Amsterdam, The Netherlands: Elsevier.

[ref23] 23. Tao W, Sheng C, Chen R, Zhang Y, Zhang S, Chang T, et al. (2011) Investigation on generations of Grapholita molesta in Changchun region by sex pheromen traps and studies on related techniques. J Jilin Agri Sci 36: 43–45.
View Article
Google Scholar

[78] View Article

[79] Google Scholar

[ref24] 24. Colinet H, Hoffmann AA (2012) Comparing phenotypic effects and molecular correlates of developmental, gradual and rapid cold acclimation responses in Drosophila melanogaster. Funct Ecol 26: 84–93.
View Article
Google Scholar

[81] View Article

[82] Google Scholar

[ref25] 25. Dickson R (1949) Factors govering the induction of diapause in the oriental fruit moth. Ann Entomol Soc Am 42: 511–537.
View Article
Google Scholar

[84] View Article

[85] Google Scholar

[ref26] 26. Heckmann LH, Sorensen PB, Krogh PH, Sorensen JG (2011) NORMA-Gene: A simple and robust method for qPCR normalization based on target gene data. Bmc Bioinformatics 12: 250. pmid:21693017
View Article
PubMed/NCBI
Google Scholar

[87] View Article

[88] PubMed/NCBI

[89] Google Scholar

[ref27] 27. Eisen MB, Spellman PT, Brown PO, Botstein D (1998) Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci Unit States Am 95: 14863–14868.
View Article
Google Scholar

[91] View Article

[92] Google Scholar

[ref28] 28. Saldanha AJ (2004) Java Treeview—extensible visualization of microarray data. Bioinformatics 20: 3246–3248. pmid:15180930
View Article
PubMed/NCBI
Google Scholar

[94] View Article

[95] PubMed/NCBI

[96] Google Scholar

[ref29] 29. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S, et al. (2011) MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 28: 2731–2739. pmid:21546353
View Article
PubMed/NCBI
Google Scholar

[98] View Article

[99] PubMed/NCBI

[100] Google Scholar

[ref30] 30. Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574. pmid:12912839
View Article
PubMed/NCBI
Google Scholar

[102] View Article

[103] PubMed/NCBI

[104] Google Scholar

[ref31] 31. Maddison WP, Maddison DR (2006) Mesquite: A modular system for evolutionary analysis. Ver. 1.1. Available via http://www.mesquiteproject.org.

[ref32] 32. Michaud S, Morrow G, Marchand J, Tanguay RM (2002) Drosophila small heat shock proteins: cell and organelle-specific chaperones? Prog Mol Subcell Biol 28: 79–101. pmid:11908067
View Article
PubMed/NCBI
Google Scholar

[107] View Article

[108] PubMed/NCBI

[109] Google Scholar

[ref33] 33. Kampinga HH, Hageman J, Vos MJ, Kubota H, Tanguay RM, Bruford EA, et al. (2009) Guidelines for the nomenclature of the human heat shock proteins. Cell Stress Chaperon 14: 105–111.
View Article
Google Scholar

[111] View Article

[112] Google Scholar

[ref34] 34. Poulain P, Gelly JC, Flatters D (2010) Detection and architecture of small heat shock protein monomers. PLoS One 5: e9990. pmid:20383329
View Article
PubMed/NCBI
Google Scholar

[114] View Article

[115] PubMed/NCBI

[116] Google Scholar

[ref35] 35. Kriehuber T, Rattei T, Weinmaier T, Bepperling A, Haslbeck M, Buchner J (2010) Independent evolution of the core domain and its flanking sequences in small heat shock proteins. FASEB J 24: 3633–3642. pmid:20501794
View Article
PubMed/NCBI
Google Scholar

[118] View Article

[119] PubMed/NCBI

[120] Google Scholar

[ref36] 36. Huang LH, Wang HS, Kang L (2008) Different evolutionary lineages of large and small heat shock proteins in eukaryotes. Cell Res 18: 1074–1076. pmid:18766171
View Article
PubMed/NCBI
Google Scholar

[122] View Article

[123] PubMed/NCBI

[124] Google Scholar

[ref37] 37. Zhu JY, Zhao N, Yang B (2012) Global transcriptome profiling of the pine shoot beetle, Tomicus yunnanensis (Coleoptera: Scolytinae). PLoS One 7: e32291. pmid:22384206
View Article
PubMed/NCBI
Google Scholar

[126] View Article

[127] PubMed/NCBI

[128] Google Scholar

[ref38] 38. Hull JJ, Geib SM, Fabrick JA, Brent CS (2013) Sequencing and de novo assembly of the western tarnished plant bug (Lygus hesperus) transcriptome. PLoS One 8: e55105. pmid:23357950
View Article
PubMed/NCBI
Google Scholar

[130] View Article

[131] PubMed/NCBI

[132] Google Scholar

[ref39] 39. Levine MT, McCoy C, Vermaak D, Lee YC, Hiatt MA, Matsen FA, et al. (2012) Phylogenomic analysis reveals dynamic evolutionary history of the Drosophila heterochromatin protein 1 (HP1) gene family. PLoS Genet 8: e1002729. pmid:22737079
View Article
PubMed/NCBI
Google Scholar

[134] View Article

[135] PubMed/NCBI

[136] Google Scholar

[ref40] 40. Mi H, Muruganujan A, Thomas PD (2013) PANTHER in 2013: modeling the evolution of gene function, and other gene attributes, in the context of phylogenetic trees. Nucleic Acids Res 41: D377–386. pmid:23193289
View Article
PubMed/NCBI
Google Scholar

[138] View Article

[139] PubMed/NCBI

[140] Google Scholar

[ref41] 41. Michaud S, Tanguay RM (2003) Expression of the Hsp23 chaperone during Drosophila embryogenesis: association to distinct neural and glial lineages. BMC Dev Biol 3: 9. pmid:14617383
View Article
PubMed/NCBI
Google Scholar

[142] View Article

[143] PubMed/NCBI

[144] Google Scholar

[ref42] 42. Basha E, O'Neill H, Vierling E (2012) Small heat shock proteins and alpha-crystallins: dynamic proteins with flexible functions. Trends Biochem Sci 37: 106–117. pmid:22177323
View Article
PubMed/NCBI
Google Scholar

[146] View Article

[147] PubMed/NCBI

[148] Google Scholar

[ref43] 43. Morrow G, Tanguay RM (2012) Small heat shock protein expression and functions during development. Int J Biochem Cell Biol 44: 1613–1621. pmid:22502646
View Article
PubMed/NCBI
Google Scholar

[150] View Article

[151] PubMed/NCBI

[152] Google Scholar

[ref44] 44. Wheeler JC, Bieschke ET, Tower J (1995) Muscle-specific expression of Drosophila hsp70 in response to aging and oxidative stress. Proc Natl Acad Sci Unit States Am 92: 10408–10412.
View Article
Google Scholar

[154] View Article

[155] Google Scholar

[ref45] 45. Jia H, Halilou AI, Hu L, Cai W, Liu J, Huang B (2011) Heat shock protein 10 (Hsp10) in immune-related diseases: one coin, two sides. Int J Biochem Mol Biol 2: 47–57. pmid:21969171
View Article
PubMed/NCBI
Google Scholar

[157] View Article

[158] PubMed/NCBI

[159] Google Scholar

[ref46] 46. Denlinger DL (2002) Regulation of diapause. Annu Rev Entomol 47: 93–122. pmid:11729070
View Article
PubMed/NCBI
Google Scholar

[161] View Article

[162] PubMed/NCBI

[163] Google Scholar

[ref47] 47. Zhang Q, Denlinger DL (2010) Molecular characterization of heat shock protein 90, 70 and 70 cognate cDNAs and their expression patterns during thermal stress and pupal diapause in the corn earworm. J Insect Physiol 56: 138–150. pmid:19782689
View Article
PubMed/NCBI
Google Scholar

[165] View Article

[166] PubMed/NCBI

[167] Google Scholar

[ref48] 48. Qiu Z, Macrae TH (2008) ArHsp21, a developmentally regulated small heat-shock protein synthesized in diapausing embryos of Artemia franciscana. Biochem J 411: 605–611. pmid:18095938
View Article
PubMed/NCBI
Google Scholar

[169] View Article

[170] PubMed/NCBI

[171] Google Scholar

[ref49] 49. Rinehart JP, Li A, Yocum GD, Robich RM, Hayward SA, Denlinger DL (2007) Up-regulation of heat shock proteins is essential for cold survival during insect diapause. Proc Natl Acad Sci Unit States Am 104: 11130–11137.
View Article
Google Scholar

[173] View Article

[174] Google Scholar

[ref50] 50. Mymrikov EV, Seit-Nebi AS, Gusev NB (2011) Large potentials of small heat shock proteins. Physiol Rev 91: 1123–1159. pmid:22013208
View Article
PubMed/NCBI
Google Scholar

[176] View Article

[177] PubMed/NCBI

[178] Google Scholar

[ref51] 51. Jakob U, Gaestel M, Engel K, Buchner J (1993) Small heat shock proteins are molecular chaperones. J Biol Chem 268: 1517–1520. pmid:8093612
View Article
PubMed/NCBI
Google Scholar

[180] View Article

[181] PubMed/NCBI

[182] Google Scholar

[ref52] 52. Lu MX, Hua J, Cui YD, Du YZ (2014) Five small heat shock protein genes from Chilo suppressalis: characteristics of gene, genomic organization, structural analysis, and transcription profiles. Cell Stress Chaperon 19: 91–104.
View Article
Google Scholar

[184] View Article

[185] Google Scholar

[ref53] 53. Bousquet F, Nojima T, Houot B, Chauvel I, Chaudy S, Dupas S, et al (2012) Expression of a desaturase gene, desat1, in neural and nonneural tissues separately affects perception and emission of sex pheromones in Drosophila. Proc Natl Acad Sci Unit States Am 109: 249–254.
View Article
Google Scholar

[187] View Article

[188] Google Scholar

Figures

Abstract

Introduction

Materials and Methods

Insect preparation

RNA extraction and cDNA synthesis

Clone of full length OFM sHsps

Real-time quantitative PCR (RT-qPCR)

Phylogenetic analysis

Phylogenetic reconstruction and signal detection

Results

Cloning and characterization of GmHsps

Phylogenetic analysis and alignment of GmHsps

Expression in response to thermal treatments and diapause

Low phylogenetic signal in lepidopteran sHsp expression patterns

Discussion

High variation of sHsp in sequence feature

The potential functions of sHsp in temporal and spatial expression

Lack of phylogenetic signal of sHsps expression for function divergence

Supporting Information

S1 Fig. Phylogenetic trees of OFM sHsps based on the conserved a-crystalline domains (about 105 amino acids) for the tree construction.

S2 Fig. Hierarchical clustering of differentially expressed genes in response to heat and cold treatments (1h and 2h, respectively).

S3 Fig. Normalized mRNA expression levels of the fourteen GmHsps in different tissues.

S4 Fig. Normalized mRNA expression of the 14 GmHsp genes in OFM at different stages.

S5 Fig. Phylogenetic tree of the lepidopteran species considered, including the state of change in the branches and tips in three characters involving heat (A), cold (B), and diapause (C) responses.

S1 Table. The GmsHsp real-time PCR primers, annealing temperatures, and expected fragment lengths.

S2 Table. List of sHsps of Lepidopteran species used in phylogenetic analysis.

S3 Table. List of sHsp expression information used in lepidopteran phylogenetic signal detection.

Acknowledgments

Author Contributions

References