Skip to main content

HSP superfamily of genes in the malaria vector Anopheles sinensis: diversity, phylogenetics and association with pyrethroid resistance

Abstract

Background

Heat shock proteins (HSPs) are molecular chaperones that are involved in many normal cellular processes and various kinds of environmental stress. There is still no report regarding the diversity and phylogenetics research of HSP superfamily of genes at whole genome level in insects, and the HSP gene association with pyrethroid resistance is also not well known. The present study investigated the diversity, classification, scaffold location, characteristics, and phylogenetics of the superfamily of genes in Anopheles sinensis genome, and the HSP genes associated with pyrethroid resistance.

Methods

The present study identified the HSP genes in the An. sinensis genome, analysed their characteristics, and deduced phylogenetic relationships of all HSPs in An. sinensis, Anopheles gambiae, Culex quinquefasciatus and Aedes aegypti by bioinformatic methods. Importantly, the present study screened the HSPs associated with pyrethroid resistance using three field pyrethroid-resistant populations with RNA-seq and RT-qPCR, and looked over the HSP gene expression pattern for the first time in An. sinensis on the time-scale post insecticide treatment with RT-qPCR.

Results

There are 72 HSP genes in An. sinensis genome, and they are classified into five families and 11 subfamilies based on their molecular weight, homology and phylogenetics. Both RNA-seq and qPCR analysis revealed that the expression of AsHSP90AB, AsHSP70-2 and AsHSP21.7 are significantly upregulated in at least one field pyrethroid-resistant population. Eleven genes are significantly upregulated in different period after pyrethroid exposure. The HSP90, sHSP and HSP70 families are proposed to be involved in pyrethroid stress response based in expression analyses of three field pyrethroid-resistant populations, and expression pattern on the time scale post insecticide treatment. The AsHSP90AB gene is proposed to be the essential HSP gene for pyrethroid stress response in An. sinensis.

Conclusions

This study provides the information frame for HSP superfamily of genes, and lays an important basis for the better understanding and further research of HSP function in insect adaptability to diverse environments.

Background

The term “heat shock protein” (HSP) represents a superfamily of genes, and their proteins as “molecular chaperons” increase in synthesis in response of unfavorable environmental condition in both prokaryotic and eukaryotic cells [1]. Most of the HSPs have conserved sequences through bacteria to human, but their expression and function subject to environmental stress do vary from organism to organism [2]. The HSPs are traditionally classified into eight families based on their molecular weight (MW) from 10 to 110 kDa and homology: HSP110 (HSPH), HSP90 (HSPC), HSP70 (HSPA), HSP40 (DNAJ), small HSP (sHSP, HSPB), and the chaperonin families HSPD (HSP60), HSPE (HSP10) and CCT (TRiC) [3]. There have been a number of functional studies on HSPs in bacteria, algae, plant, amphibians, birds, and mammalian, especially in the model organisms of Anopheles gambiae, Drosophila melanogaster, Arabidopsis thaliana, Saccharomyces cerevisiae, Caenorhabditis elegans, Danio rerio and Mus musculus [4,5,6,7]. The HSP genes of insects encode molecular chaperones that help repair stress injuries via transportation and degradation of aggregated proteins, and they are susceptible to environment stresses, such as heat or cold stress, and insecticidal infection in insects [8, 9]. The HSP expression patterns have been identified in Drosophila species and some other insects [10,11,12,13]; however, the information regarding HSP whole-genome diversity and association with insecticide resistance is still quite limited for insects as well as mosquitoes.

Recent research showed that HSP90 was not only an important gene in Apolygus lucorum adults in response to extremely high temperature, but was also associated with the resistance or tolerance to cyhalothrin, imidacloprid, chlorpyrifos, and emamectin benzoate and cyhalothrin [14]. Two HSP70s were upregulated in a chlorpyrifos-resistant population of Plutella xylostella, whereas six sHSPs were downregulated [15]. The sHSPs expression in responses of indoxacarb and cantharidin varied, and the exposure to beta-cypermethrin and chlorfenapyr resulted in an increase of 13 sHSP transcripts and a reduction of 12 sHSP transcripts in Plutella xylostella, respectively, which indicates that different sHSPs might play distinct roles in the development and regulation of physiological activities [16]. One HSP70 was highly expressed in the captafol-exposed larvae of Drosophila melanogaster [17]. In the imidacloprid and deltamethrin treatment groups, only one HSP90 was upregulated more than two-fold in response to deltamethrin treatment in Sogatella furcifera [18]. Two HSPs were up-regulated in DDT-resistant field isolates in An. gambiae, supporting their link with insecticide resistance and/or stress response [19]. HSPs are induced in response to stress induced by environmental factors and may be involved in adverse reactions, including insecticide resistance [20].

Mosquitoes transmit disease resulting in around 700 million people of mosquito-borne illness each year and over one million deaths [21]. Mosquito control relies primarily on the use of insecticides through insecticide-impregnated bed nets and indoor residual spray. In the past decades, pyrethroids have become the preferred insecticide because of their low toxicity to humans, high efficacy against mosquito vectors and short residual action. However, there has been increasing number of mosquitoes to have developed resistance to pyrethroids [22]. Anopheles sinensis is a major malaria vector in China and other Southeast Asian countries [23]. Because of extensive and continued application, pyrethroid resistance in An. sinensis is now widespread in many regions of China [24, 25]. However, the molecular mechanisms of pyrethroid resistance in An. sinensis are not yet clearly understood and pose a challenge for the control of malaria in China. Recently, An. sinensis genome was deeply sequenced and assembled, and are comprehensively conducting the research on the molecular mechanism of pyrethroid resistance using the species as model. Some families of genes have been identified at whole-genome level, and their association with pyrethroid resistance have been investigated or summarized, such as cytochromes P450s [26] and carboxylesterases [27]. However, there is still no report for the genome-wide identification of the HSP superfamily of genes in insects as well as in An. sinensis. The association of HSPs with pyrethroid resistance has been little known in insects as well as in An. sinensis.

The present study identified and classified the HSPs superfamily of genes in An. sinensis genome, and analysed their phylogenetics and basic characteristics. More importantly, this study screened the HSP genes associated with pyrethroid resistance at whole-genome level through transcription comparison between resistant- and susceptible populations/strain from three geographical regions, quantitative PCR confirmation of candidate HSP genes, and expression changes of HSP genes in response to pyrethroid exposure. This study provides a comprehensive information frame for the HSP superfamily of genes in the An. sinensis and other insects, and lays a foundation for the functional study of HSPs in relation of pyrethroid resistance.

Methods

Whole-genome identification of Anopheles sinensis HSP genes

The HSP genes were identified from the genome and two sets of transcriptome sequences of An. sinensis. The former was achieved by Institute of Entomology and Molecular Biology, Chongqing Normal University, China (publication in preparation), and the later were downloaded from the GenBank [24, 28]. Three approaches were used to the identification as previously described for the gene searching at whole-genome level in other species. First, this study searched the An. sinensis amino acid (a.a.) database and nucleotide database using BLAST (Basic Local Alignment Search Tool) [29] with E-value cut off of 1 × 10−5 using HSP sequences of Homo sapiens and An. gambiae as query. Second, An. sinensis HSP sequences retrieved by five hidden Markov model (HMM) for the HSP superfamily (HSP90, PF00183; HSP70, PF00012; HSP60, PF00118; HSP40, PF01556; HSP20, PF00011) to search against An. sinensis a.a. database using Pfam (version 27.0) (http://pfam.sanger.ac.uk/) [30]. Third, this study combined the sequences retrieved from two approaches above as query to search against the An. sinensis genome assembly with an E-value cut-off of 1 × 10−5. Then putative HSP sequences were predicted by Fgenesh+ (http://www.softberry.com/) to determine whether they are full-length or partial sequences of genes. Full-length sequences of the partial sequences of genes were determined by extending the flanking regions (1 kb or longer), and manually corrected by comparison with other homologous HSP genes. Finally, the HSP genes and their sequences were determined in combining the results from these steps above. The scaffold position of each HSP gene was mapped to the An. sinensis genome using BLAST searching, and the scaffold distribution of HSPs was drawn manually. The scaffolds containing only one gene were not displayed, and the gene clusters were defined by the presence of four or more genes in a region of less than 200 kb [31].

Phylogenetics analysis and characterization of An. sinensis HSPs

The deduced HSP a.a. sequences of An. sinensis, and the HSP a.a. sequences of An. gambiae, Aedes aegypti, and Culex quinquefasciatus downloaded from NCBI (http://www.ncbi.nlm.nih.gov/database) were used to refer the phylogenetic relationships of HSP genes. The alignment of full amino acid sequence was conducted using ClustalW2.0 program with the default Needle algorithm [32] and corrected manually using GeneDoc [33]. Phylogenetic trees were constructed using the neighbor-joining (NJ) algorithm with the Poisson model, p-distance and pairwise deletion of gaps using MEGA5.0 [34]. The maximum likelihood (ML) method was also used to infer the phylogenetic relationships using MEGA5.0 with the best-fit model (JTT + G+I + F) predicted by the Modeltest 3.7 [34]. Tree topology was assessed by bootstrap analysis with 1000 resampling replicates. The deduced a.a. sequences were further confirmed by the analysis of its characteristics. The theoretical molecular weights and isoelectric points of the predicted a.a. were calculated using the ExpaSy server (http://www.expasy.ch/tools/pi tool.html). The signal peptide sequences were predicted using SignalP 4.1 server (http://www.cbs.dtu.dk/services/SignalP/), and the subcellular localization using the Wolfpsort (http://www.genscript.com/psort/wolf_psort.html). The conserved domains were examined using the program Pfam (version 27.0) (http://pfam.sanger.ac.uk/) [30], and the gene structures were determined from the Gene Structure Display Server (http://gsds.cbi.pku.edu.cn/).

Screening of HSP genes associated with pyrethroid resistance using RNA-seq

Four populations/strain of An. sinensis were used to screen HSP genes associated with pyrethroid resistance using RNA-seq. They are three field pyrethroid-resistant populations from Anhui Province (AH-FR), Chongqing Municipality (CQ-FR), and Yunnan Province (YN-FR), and one laboratory pyrethroid-susceptible strain originally collected from Wuxi, Jiangsu Province (WX-LS). The mosquito larvae and pupae from fields were collected from rice fields, and were reared to adults in an insectary under 27 ± 1 °C and 70 ± 10% RH. The female adults 3 day post emergence, which were preliminarily identified to An. sinensis using morphological characters, were tested for pyrethroid susceptibility using the standard WHO tube bioassay using test papers with a diagnostic concentration of deltamethrin (0.05%) (CAS: 52918-63-5, Sigma-Aldrich, St. Louis, MO, USA) [35]. The individuals that were still alive post 24 h of recovery after 1 h 0.05% deltamethrin treatment were considered to be deltamethrin resistant. All individuals in the An. sinensis susceptible laboratory strain died in the first 1 h 0.05% deltamethrin treatment. Prior to RNA extraction, genomic DNA was extracted from two or three legs of each individual mosquito for molecular identification using the ribosomal DNA internal transcribed spacer 2 (rDNA-ITS2)-based method to confirm the morphological identification [36].

At least 100 female adults 3 day post emergence from each of these three pyrethroid-resistant populations and the laboratory strain were preserved in RNAlater (Qiagen Shanghai, China) for RNA extraction and RNA-seq library construction performed as previously reported [28], one pool of 15 individuals was used for the RNA-seq, respectively. After RNA quality was assessed using an Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA), the cDNA for each library was synthesized and amplified using the Mint-2 cDNA synthesis kit (Evrogen, Moscow, Russia). Illumina HiSeq™ 2000 was used for cDNA library sequencing at Beijing Genomics Institute (BGI), following to manufacturer’s instructions. Reads obtained from Illumina sequencing were mapped to the An. sinensis genome using TopHat [37], and expression were determined in term of fragment per kb per million reads (FPKM) using Cufflinks [37, 38]. Differential accumulation of transcripts between deltamethrin-resistant and -susceptible mosquitoes was assessed by the Cuffdiff program within Cufflinks. To minimize the impact of sequencing length and nucleotide composition, FPKM for each gene of each sample were calculated to determine the expression quantity [38]. The gene comparison pair (field-resistant/library-susceptible sample) with the FPKM fold change ≥ 2 [39, 40] and the p value ≤ 0.05 [− Log10(p-value) ≥ 1.301] were considered to significantly up-regulated [41]. In contrast, the pair with FPKM fold change ≤ 2 and the p-value ≤ 0.05 ≥ 1.301] were treated as significantly down-regulated [41].

qPCR verification of potential HSP genes associated with pyrethroid resistance

The genes significantly up-regulated at least one resistant population in the RNA-seq analysis and those genes earlier reported associated with insecticide/deltamethrin resistance were chosen to carry out reverse-transcription quantitative-PCR (RT-qPCR) verification. These three pyrethroid resistant populations and the laboratory susceptible strain were investigated for the expression confirmation of selected genes. The RT-qPCRs were conducted with three biological replicates (three mosquitoes per sample) and three technique replicates. Total RNA was isolated from dehydrated/rehydrated samples using Trizol Reagent (Invitrogen, Shanghai, China) following the supplier’s instructions. RNA concentration was measured using a Nanodrop-1000 spectrophotometer (Thermo scientific). For each sample, complementary DNAs (cDNAs) were synthesized from 1.0 μg RNA, treated with PrimScript ™ RT Reagent Kit with gDNA Eraser (TaKaRa, Dalian, China) and stored at − 20 °C. Gene-specific RT-qPCR primers were designed using Primer Premier 5.0 against the An. sinensis genome sequence (Table 1). Real-time reactions were conducted on a thermal cycler (CFX, Bio-Rad, USA) in a 20 μL reaction system containing 10.0 μL of 2 × qPCR mix (Bio-Rad, USA), 0.8 μL each of gene-specific primers (10 mM each) and 1 μL the cDNAs templates, and 7.4 μL of double distilled water. Thermal cycling conditions were 94 °C for 3 min; 40 cycles of 95 °C for 5 s, 60 °C for 15 s, and 72 °C for 15 s, and followed by melting curve analysis (60 to 95 °C). The relative normalized expression was calculated using Bio-Rad CFX96 software. Samples were normalized using the ribosomal protein S7 (RPS7) and the ribosomal protein L49 (RPL49) Ct values. Basal expression levels were represented as folds over the expression levels of RPS7 and RPL49. Expression folds were calculated with the 2−∆∆Ct method [42] between treatment and control samples for each biological replicate. All data were presented as mean ± SD (standard deviation) of three biological and three technical replicates. Significant differences between treatment groups and the control group were analysed using Student’s t test with p ≤ 0.05 pairs. One-way ANOVA was used for multiple comparisons.

Table 1 Primers used for qPCR verification for the potential genes associated with pyrethroid resistance, detected by RNA-seq

Expression analysis of HSP90 and sHSP genes under pyrethroid stress

All four HSP90 family of genes and all ten sHSP family of genes identified in the present study were assessed for the expression response on time-scale subject to pyrethroid treatment, because these two families of genes had more reports to be associated with insecticide resistance. The fourth-instar stage larvae of laboratory pyrethroid-susceptible strain were employed for the assessment. The 50 percent lethal (LC50) concentration with deltamethrin was first determined as the pyrethroid treatment concentration of the larval samples. Five pyrethroid concentrations (0.00059, 0.00118, 0.00256, 0.00512, 0.01000 mg/L) were applied for the determination under the standard insectary condition indicated above. For each concentration, 30 larvae without feeding were placed in 200 ml of treated or untreated water in trays with four replicates, larval mortality was recorded 24 h later, and the LC50 was calculated based on the mortality using WHO standard method [43]. Thirty-fourth larvae were subsequently treated with the LC50 of deltamethrin or untreated water as control for 24 h, and then moved to rearing-conditioned water with four replicates. The samples were collected at 1 h, 3 h, 6 h, 12 h, 24 h, 36 h and 48 h post moving to water, and were subject to expression analysis using qPCR technique with three biological and three technical replications as stated above.

Results

HSP superfamily of genes and their classification in Anopheles sinensis genome

A total of 72 HSP genes are identified from An. sinensis genome and transcriptome database (Additional file 1: Table S1). Among these genes, 71 genes have complete open reading frame (ORF) and are supported by transcripts, and one gene (AsDNAJB9) also has complete ORF sequence but no transcript has been found. The deduced a.a. sequence of AsDNAJB9 shares 46% sequence identity to other reported insect HSP40 s and has the conserved domain for HSP40 recognition; therefore, it is also considered as a true HSP40 gene. These 72 genes were classified into five families (HSP90, HSP70, chaperonins, HSP40 and sHSP) and 11 subfamilies based on the molecular weight and homology analysis. The nomenclature of HSP90 and HSP70 family also refers the nomenclature for An. gambiae HSP genes retrieved from VectorBase, and the nomenclature of chaperonins and HSP40 follows the nomenclature of the HUGO Gene Nomenclature Committee. There are 69, 69 and 88 HSP genes in An. gambiae, Cx. quinquefasciatus and Aedes aegypti genomes in the present study, respectively, and they are also classified into five families and 11 subfamilies (Table 2). There are 40, 36, 25 and 35 genes in HSP40 family, 11, 11, 14 and 13 genes in Chaperonins family, 10, 10, 14 and 24 genes in sHSP family, 7, 8, 11 and 8 genes in HSP70 family, and 4, 4, 5 and 8 genes in HSP90 family, separately in An. sinensis, An. gambiae, Culex quinquefasciatus and Aedes aegypti. It appears that the HSP40 family is the largest family in term of gene number, and the HSP90 family is the smallest family for all of these four species investigated.

Table 2 Family, subfamily and number of the HSP genes in An. sinensis, An. gambiae, Cx. quinquefasciatus and Ae. aegypti genomes

Gene structure and genomic location of Anopheles sinensis HSP genes

The 72 HSP genes have a total of 220 exons with each gene having 1–15 exons, and their exon sizes range from 3 bp in chaperonins family to 2593 bp in HSP70 family (Additional file 1: Table S1). In these genes, 17 genes (23.6% of total) have only one exon (i.e. no intron), 36 genes (50.0%) each have 2–3 exons (i.e. 1–2 introns), 15 genes (20.8%) each have 4–6 exons, and four genes have especially bigger number of exons (7 for HSPH1 and DNAJC23, nine for DNAJC20, and 15 for DNAJC13) (Fig. 1a, Additional file 2: Fig S1). Only HSP40 family has genes with exon number larger than five (except for HSPH1 in HSP70 family), which might be due to large diversity in the family (40 genes, 55.6% of total). The 17 genes with only one exon exist in the HSP90 (1 gene), HSP70 (1), chaperonins (2), HSP40 (9) and sHSP (4) family, respectively. There are a total 148 introns in the 72 HSP genes, with sizes ranging from 19 to 4790 bp and 37.4% intron having a length of 60–69 bp (Fig. 1b).

Fig. 1
figure 1

Number/frequency distribution of exon number (a) and intron length (b) of the 72 HSP genes in Anopheles sinensis

The 72 HSP genes are located on 33 different scaffolds, in which 24 genes are on scaffold6, scaffold15, scaffold25, scaffold63 and scaffold116 (Fig. 2), and 68 genes scattered on 28 other scaffolds with intergenic distance larger than 200 kb. Nine genes in sHSP family are tandemly located in one cluster on scaffold25 with intergenic distance less than 200 Kb, which suggests that these nine genes might be derived from a single gene through a series of gene duplication events. Similarly, AsHSP90AA and AsHSP90AB in HSP 90 family on scaffold15, AsDNAJC5 and AsDNAJB9 in HSP40 family on scaffold63, and AsDNAJB1 and AsDNAJB2 in HSP40 family on scaffold116 each cluster together with intergenic distance less than 200 Kb, which suggests that each of these three pairs might experience one gene duplication event. The results reveal that the gene expansion in An. sinensis mainly happened in the sHSP family for the function that needs to be further investigated. Most of sHSP genes are also located on one chromosome (chromosome 2) of An. gambiae (6 of 10 genes) [47]. The genome of An. sinensis is comparative with An. gambiae in gene location and possible gene expansion in the sHSP family. The gene expansions in other families need to be further investigated with more insect species.

Fig. 2
figure 2

The location of Anopheles sinensis HSP genes on Scaffolds. The left side on scaffold (S) denotes physical locations, and right side is for gene names. One individual gene cluster is shown in Scaffolds25

Phylogenetic relationships of Anopheles sinensis HSP genes

The phylogenetic analysis with NJ method generally classifies the 72 HSP genes in An. sinensis into five families as expected: HSP90, HSP70, Chaperonins, HSP40 and sHSP (Fig. 3). The HSP90, HSP70 and sHSP families of clades are each supported by a 100% bootstrap values, which suggests that these three families each be a monophyly; however, the Chaperonins and HSP40 families of branches are only supported by 33% and 20% bootstrap values, respectively, which suggests that these two families be not monophyly. The phylogenetic relationships well fit the modern classification at family level [3]. This study further investigated and discussed the phylogenetic relationships of each family of genes with inclusion of the HSP genes in An. gambiae, Culex quinquefasciatus and Aedes aegypti using amino acid sequences with ML method.

Fig. 3
figure 3

Phylogeny of the HSP genes in Anopheles sinensis. The phylogenetic relationships were inferred based on the amino acids with neighbor-joining method using MEGA5.0. The bootstrap values calculated from 1000 replicates were marked on each corresponding node. HSP90, HSP70, Chaperonins, HSP40 (DNAJ) and sHSP represent five different families in the HSP superfamily

HSP90 family (Additional file 3: Fig S2A)

Twenty-one HSP90 genes from An. sinensis (4 genes), An. gambiae (4), Culex quinquefasciatus (5) and Aedes aegypti (8) are grouped into three subfamilies, HSP90A (13 genes), HSP90B (4) and TRAP (4). All these three subfamilies are supported by 100% bootstrap value, which suggests that these all be monophyly.

HSP70 family (Additional file 3: Fig S2B)

The HSP70 family of genes are currently classified into two subfamilies, HSP70 and HSP110 [3], and these genes function in cytoplasm, mitochondrion and endoplasmic reticulum (ER) [50]. The present research based 34 HSP70 genes suggests that the HSP70 subfamily be a monophyly with 100% bootstrap value, whereas the HSP110 subfamily be a paraphyly. The HSP70 subfamily of genes are further divided into three groups based on their functional sites, cytoplasm, mitochondrion and ER group, which are all supported by 100% bootstrap value. Similarly, the HSP110 subfamily of genes are further divided into cytoplasm and ER group, which are both supported by 100% bootstrap value.

Chaperonins family (Additional file 3: Fig S2C)

The chaperonins family of genes are nowadays classified into three subfamilies, CCT, HSPD and HSPE [3]. This research based on 49 genes suggests that the CCT, HSPD and HSPE subfamilies be monophyly, supported by 94%, 100% and 100% bootstrap values, respectively. The CCT subfamily of genes is further divided into eight groups (A to H), which are all supported at least 96% bootstrap value.

HSP40 (DNAJ) family (Additional file 3: Fig S2D)

The HSP40 family of genes is used to be classified into three subfamilies, DNAJA, DNAJB and DNAJC [3]. The domain structure of all HSP40 genes identified in An. sinensis (Additional file 4: Table S2). The present study based on 136 genes suggest that the DNAJA, DNAJB and DNAJC be monophyly supported by 83%, 99% and 100% bootstrap values, respectively.

sHSP gene family (Additional file 3: Fig S2E)

A total of 58 sHSP genes are classified into eight clusters in the phylogenetic analysis, four orthologous clusters and four species-specific clusters. There are 1–2 genes in each of these four orthologous clusters for the four mosquito species investigated. For the four species-specific clusters, there are six, four, ten and ten genes in An. sinensis, An. gambiae, Culex quinquefasciatus and Aedes aegypti, respectively.

HSP genes associated with pyrethroid resistance

Transcript accumulation levels in a RNA-seq experiment are expected to reflect gene transcription. Comparing each of these three field pyrethroid-resistant populations against the laboratory susceptible strain (WX-LS), four HSP genes (AsHSP90AB, AsHSP70-2, AsHSP21.7 and AsDNAJB4) are significantly overexpressed in at least one population with p-values ≤ 0.05 in T-test (Fig. 4). The AsHSP90AB is significantly upregulated in all three populations with the folds of 1.58 (in AH-FR), 3.63 (CQ-FR) and 3.74 (YN-FR), the AsHSP70-2 significantly upregulated in both AH-FR (2.56 folds) and YN-FR (3.56), and the AsHSP21.7 and AsDNAJB4 alone in YN-FR with the folds of 1.77 and 1.55, respectively. All of the four HSP genes are selected to perform RT-qPCR for validation using the samples consistent with those for RNA-seq analysis.

Fig. 4
figure 4

Expression profile of HSP genes in three pyrethroid-resistant populations compared to pyrethroid-susceptible strain (WX-LS) Anopheles sinensis, detected by RNA-seq. The three resistant populations are Anhui (AH-FR, marked in blue), Chongqing (CQ-FR, red) and Yunnan (YN-FR, green), respectively. Genes with FPKM value ≥ 2 and p-value < 0.05 are considered to be significantly upregulated, whereas those with FPKM value < 0.5 and p-value < 0.05 are regarded as significantly downregulated. Vertical dotted lines illustrate the 1 (FPKM = 2) and − 1 (FPKM = 0.5) of Log2 (fold changes of FPKM value) values on the X-axis, and horizontal dotted lines denote the p-value = 0.05 of − Log10 (p-value) on the Y-axis. AsHSP90AB is significantly upregulated in all three resistant populations, AsHSP70-2 is significantly upregulated in both AH-FR and CQ-FR and AsHSP21.7 and AsDNAJB4 are significantly upregulated only in YN-FR

The RT-qPCR verification results are largely consistent with those from RNA-seq, and some differences might stem from the different batch of samples, target gene splicing and/or method difference (Fig. 5). Most importantly, the AsHSP90AB in HSP90 family is also significantly up-regulated in all three resistant populations in the RT-qPCR analysis. The AsHSP70-2 is significantly up-regulated also in AH-FR but not in YN-FR in the RT-qPCR verification. The AsHSP21.7 in sHSP family is significantly upregulated in CQ-FR but not in YN-FR in the RT-qPCR verification. The AsDNAJB4 in HSP40 family has no significant expression difference for any of the three resistant populations in the RT-qPCR verification in comparison of the susceptible strain.

Fig. 5
figure 5

qPCR verification of four genes significantly differentially expressed in RNA-seq analysis. The relative expression levels of these genes are shown as the mean ± SD of three biological and three technical replicates in qPCR analysis. The RPS7 and RPL49 genes were used as internal reference for expression normalization. The population/strain pairs marked as different letters are significantly different in expression (p-value ≤ 0.05), and those marked as same letters are not (p-value ≥ 0.05), determined by one-way ANOVA analysis

Expression response of HSP90 and sHSP genes under pyrethroid stresses

The fourth-instar stage larvae of laboratory pyrethroid-susceptible strain were treated with 0.00160 mg/L of deltamethrin, based on the determination of LC50 concentration. The obvious expression pattern on the time scale post pyrethroid treatment was detected for the four HSP90 genes and ten sHSP genes. The AsHSP90AB, AsHSP25.6, AsHSP19 and AsHSP21 are significantly up-regulated by 1.33 to 9.47 folds through 1 h to 48 h pyrethroid treatment (Fig. 6); the AsHSP21.7 is significantly upregulated by 1.44 to 2.47 folds through 1 h to 24 h; the AsHSP13.7 is significantly upregulated by 1.99 to 3.38 folds through 6 h to 48 h; the AsHSP90AA, AsTRAP, AsHSP23.5 and AsHSP20.5 are significantly upregulated by 1.67 to 6.90 folds through 12 h to 48 h; and the AsHSP24.5 is significantly upregulated by 1.48 to 2.99 folds through 24 h to 48 h. However, the AsHSP24.7 is significantly downregulated through 3 h to 24 h, and the AsHSP90B1 and AsHSP17.8 show no significant expression variations through 1 h to 48 h. This is the first time to look over the HSP gene expression pattern on the time scale post-insecticide treatment, and the study suggests that both HSP90 and sHSP families genes be overtranscribed across the detoxification process after pyrethroid exposure with specific expression pattern. There is a need to further investigate the involvement of other HSP families of genes.

Fig. 6
figure 6

Expression of 14 HSP genes subject to pyrethroid exposure, detected by qPCR. The 14 genes include four from HSP90 family and ten from sHSP family. The fourth-instar larvae at 1 h, 3 h, 6 h, 12 h, 24 h, 36 h and 48 h post pyrethroid treatment were investigated using qPCR analysis with three biological and three technical replications. The expression levels (mean ± SD) of each gene, represented with bars, are normalized in reference of RPL49 and RPS7. The samples marked with “*” were significantly differently expressed (p-value ≤ 0.05) compared with the samples without pyrethroid treatment at the corresponding time, and the samples without “*” were not significantly different in expression determined by one-way ANOVA analysis (p ≤ 0.05)

Discussion

In history, the HSP superfamily of genes across all organisms was classified into several families based on their molecular weights, HSP105/110, HSP90, HSP70, HSP60, HSP40, small HSP (sHSP) and HSP10 [44, 45]. This is the first report for whole-genome identification of HSP superfamily of genes in the four species as well in insects. However, the HSP110 family does not exist in insect [45], and the HSP10 (HSPE) and HSP60 (HSPD) families have been degraded as subfamilies in insects [46]. The HSP gene number (72) of An. sinensis is comparative with that of An. gambiae (69) and Culex quinquefasciatus (69), but much less than that of Aedes aegypti (88). HSPs have been shown to be markedly associated with the resistance to heavy metals, pesticides and oxidative stress in insects, and the difference of gene number might result from the adaptation to different environment [40, 46]. The much larger number in Aedes aegypti appears due to the expansion of sHSP family (24 genes in Aedes aegypti, and 10–14 in other three species) and HSP90 family (8 genes in Aedes aegypti, and 4–5 in other three species).

The established phylogenetic relationship of HSP90 family genes are consistent with earlier studies, in which the HSP90A is a sister with the HSP90B, and the TRAP originated earlier than both HSP90A and HSP90B [48, 49, 69]. Both the HSP90A and HSP90B exist in all eukaryotic kingdoms, and the TRAP also in bacteria [49]. In this study, the four mosquito species investigated are lack of HTPG and HSP90C subfamilies of genes, like in Drosophila melanogaster [49]. In the HSP70 family, the C-terminal motif contains with diverse subcellular localizations signatures [50]. In HSP110 subfamily have similar domains as canonical HSP70 subfamily but have long insertions and C-terminal extensions [50]. The ATPase domain and the C-terminal helical lid of HSP110 subfamily were thought to mediate the interaction with HSP70 subfamily [50]. In the chaperonins family, the CCT subfamily of genes in cytoplasm is a multi-subunit protein complex which functions in cytoskeletal protein folding in all eukaryotes [51], and the HSPD and HSPE subfamilies in the mitochondria are orthologs of the Escherichia coli GroEL (HSP60) and GroES (HSP10), respectively [3]. It has been reported that HSPE (HSP10) proteins serve as the co-factor of HSPD (HSP60) to assist in the folding of newly synthesized proteins imported into mitochondria [3]. In the HSP40 family, classified into three subfamilies, DNAJA, DNAJB and DNAJC are in accordance with those of previous studies, Bombyx mori [52]. For sHSP family, four orthologous clusters and four species-specific clusters was found. This clustering pattern is similar to those in Bombyx mori and An. gambiae [47].

Previous studies showed that HSP superfamily members were differently expressed in diverse insecticides of insects. This may be correlated with the facts that HSP genes were significant in response to insecticides resistance of insects. This study checked the expression profiles of HSP genes in three field pyrethroid-resistant populations against the laboratory susceptible strain of An. sinensis, and the expression patterns of HSP genes verified by qPCR. Similarly, the gene has been reported to be significantly overexpressed in chlorpyrifos-resistant resistance strains in Plutella xylostella (HSP90) [53], and in response of DDE induces in Ruditapes decussatus (HSP90) [54] and abamectin treatment in Tetranychus cinnabarinus (HSP90) [65]. In addition, the increased expression of HSP90 has been associated with pesticide exposure in Apis mellifera [55]. The transcriptional expression profile results also indicated that the expression of AlHSP90 in female adults treated with chlorpyrifos and emamectin benzoate and in male adults treated with cyhalothrin were higher than that with other treatments [14]. Earlier studies show that one HSP70 gene was significantly upregulated in the DDT-resistant strains of Aedes aegypti [56], two HSP70s in a chlorpyrifos-resistant population of Plutella xylostella [15], one HSP70 in organophosphorus insecticide resistant Chironomus yoshimatsu [20], and two HSP70s in chlorpyrifos-resistant Laodelphax striatella [57]. One HSP70 gene has been reported to involve in cellular damage in reproductive tissues induced by cypermethrin insecticide in Drosophila melanogaster [58], and Colorado potato beetles produce more HSP70 in response of imidacloprid [59]. The AsHSP70-2 is significantly up-regulated in AH-FR, and a number of HSP70 family of genes were significantly upregulated in populations/strains in a number of species. These findings suggest that the HSP70 families of genes might be also involved in stress response process. These results suggest that HSP70 expression be a sensitive indicator of exposure to certain insecticides and in conjunction with other biomarkers, and it may be useful for assessing exposure to environmental stressors ecosystems [20]. However, HSP70 genes involved might be different along species, geographical populations and insecticides. It is significantly upregulated in response of induce of cypermethrin insecticide in Caenorhabditis elegans (HSP16) [60]. The expression levels of sHSP19.7 and sHSP20.7 in cultured cells of Mamestra brassicae were significantly up-regulated in response to high concentrations of chlorfenapyr [61]. However, six sHSPs were downregulated in a chlorpyrifos-resistant population of Plutella xylostella [15], and one sHSP was downregulated in response of imidacloprid treatment of Sogatella furcifera [18]. It appears that some sHSP genes are responsible for the defense against insecticide stress, but the gene response vary upon insecticide and insect species. The AccDnaJB12 was upregulated from 1 to 1.5 h in response of lambda-cyhalothrin and paraquat treatment in Apis cerana cerana [62]. The OcHsp40 mRNA levels had no significant difference observed with Cd concentrations of Oxya chinensis [63].

In arthropods, HSP90 proteins have been shown to be involved in tolerance and resistance to pesticides [64, 65]. HSP90 family of genes has been known to play a role in protein folding and posttranslational regulation, in particular in steroid hormone targeting and cell death and apoptosis regulating [66]. This study reveals that the AsHSP90AB is significantly upregulated in the all three pyrethroid-resistant populations investigated and through 1 h to 48 h post pyrethroid exposure, and this result suggests that the AsHSP90AB be the essential HSP gene for pyrethroid stress response. The AsHSP90AA and AsTRAP in HSP90 family might also be involved in pyrethroid stress response. However, three HSP90 genes and two sHSP genes down-regulated during permethrin exposure on Anopheles stephensi third instar larvae [70]. The AsHSP21.7 is significantly upregulated in CQ-FR, and a number of sHSP family of genes are significantly over-expressed over the detoxification process post pyrethroid exposure. Twelve of 14 sHSPs genes were significantly up-regulated in the fourth instar larvae of Plutella xylostella after beta-cypermethrin exposure [15]. These findings suggest that the sHSPs families of genes might be also associated with insecticide stress process. HSP90 and sHSP families genes expression patterns were different in different insecticides and different concentration in different insects.

HSPs are highly demanded by the insects to combat environmental stresses, and are associated with excess expressions of apoptotic genes under insecticide stress, which results in higher apoptosis [53]. In the presence of abiotic and biotic stressors, HSPs upregulated are thought to participate in stress tolerance and promote cell survival mainly through refolding proteins and preventing their denaturation [67, 68]. The difference in expression pattern of these HSPs may be due to a compensation effect among the HSP genes [53]. The insect HSPs response to insecticide stress has received increasing attention [49, 53, 62]. In addition, inducible HSPs as playing a vital role in the insecticide resistance phenotype [9]. Previous studies showed that blood feeding increases HSP expression in mosquito, that may suggest a role for HSPs in blood meal-induced reduction of insecticide toxicity [4]. As molecular chaperones, HSPs may be induced by insecticides, and they also contribute to insecticide resistance [20].

Conclusions

This is the first study to look over the whole HSP superfamily of genes in insects. The study provides useful insights into the diversity, classification, scaffold location, characteristics, and phylogenetics of HSP superfamily of genes in An. sinensis genome, and the HSP genes associated with pyrethroid resistance. There are a total of 72 HSP genes in An. sinensis, and they are classified into five families (HSP90, HSP70, Chaperonins, HSP40 and sHSP) and 11 subfamilies based on the molecular weight, homology and phylogenetics analyses. The HSP90, sHSPs and HSP70 families of genes are proposed to be involved in pyrethroid stress response based on expression analyses of three field pyrethroid-resistant populations, and expression pattern on the time scale post insecticide treatment. The AsHSP90AB is proposed to be the essential HSP gene for pyrethroid stress response in An. sinensis. This study provides the information frame for HSP superfamily of genes, and lays an important basis for the better understanding and further research of HSP function in adaptability of insects to diverse environments.

Abbreviations

HSPs:

heat shock proteins

a.a.:

amino acid

ML:

maximum likelihood method

NJ:

neighbour joining

FPKM:

fragment per kb per million reads

HSPH:

heat shock protein 110

HSPC:

heat shock protein 90

HSPA:

heat shock protein 70

CCT:

chaperonin-containing tailless complex polypeptide 1

HSPD:

heat shock protein 60

HSPE:

heat shock protein 10

HMM:

Hidden Markov model

WHO:

World Health Organization

References

  1. Tissieres A, Mitchell HK, Tracy UM. Protein synthesis in salivary glands of Drosophila melanogaster: relation to chromosome puffs. J Mol Biol. 1974;84:389–98.

    Article  CAS  Google Scholar 

  2. Pratt WB, Toft DO. Regulation of signaling protein function and trafficking by the hsp90/hsp70-based chaperone machinery. Exp Biol Med. 2003;228:111–33.

    Article  CAS  Google Scholar 

  3. Kampinga HH, Hageman J, Vos MJ, Kubota H, Tanguay RM, Bruford EA, et al. Guidelines for the nomenclature of the human heat shock proteins. Cell Stress Chaperones. 2009;14:105–11.

    Article  CAS  Google Scholar 

  4. Benoit JB, Giancarlo LM, Patrick KR, Phillips ZP, Krause TB, Denlinger DL. Drinking a hot blood meal elicits a protective heat shock response in mosquitoes. Proc Natl Acad Sci USA. 2011;108:8026–9.

    Article  CAS  Google Scholar 

  5. Beckmann RP, Mizzen LE, Welch WJ. Interaction of Hsp 70 with newly synthesized proteins: implications for protein folding and assembly. Science. 1990;248:850–4.

    Article  CAS  Google Scholar 

  6. Kelley WL. The J-domain family and the recruitment of chaperone power. Trends Biochem Sci. 1998;23:222–7.

    Article  CAS  Google Scholar 

  7. Walsh P, Bursac D, Law YC, Cyr D, Lithgow T. The J-protein family: modulating protein assembly, disassembly and translocation. EMBO Rep. 2004;5:567–71.

    Article  CAS  Google Scholar 

  8. Rinehart JP, Li A, Yocum GD, Robich RM, Hayward SA, Denlinger DL. Up-regulation of heat shock proteins is essential for cold survival during insect diapause. Proc Natl Acad Sci USA. 2007;104:11130–7.

    Article  CAS  Google Scholar 

  9. Oliver SV, Brooke BD. The effect of elevated temperatures on the life history and insecticide resistance phenotype of the major malaria vector Anopheles arabiensis (Diptera: Culicidae). Malar J. 2017;16:73.

    Article  Google Scholar 

  10. Feder ME, Hofmann GE. Heat-shock proteins, molecular chaperones, and the stress response: evolutionary and ecological physiology. Annu Rev Physiol. 1999;61:243–82.

    Article  CAS  Google Scholar 

  11. Kregel KC. Heat shock proteins: modifying factors in physiological stress responses and acquired thermotolerance. J Appl Physiol. 2002;92:2177–86.

    Article  CAS  Google Scholar 

  12. Zheng H, Nagaraja GM, Kaur P, Asea EE, Asea A. Chaperokine function of recombinant Hsp72 produced in insect cells using a baculovirus expression system is retained. J Biol Chem. 2010;285:349–56.

    Article  CAS  Google Scholar 

  13. Tower J. Heat shock proteins and Drosophila aging. Exp Gerontol. 2011;46:355–62.

    Article  CAS  Google Scholar 

  14. Sun Y, Sheng Y, Bai L, Zhang Y, Xiao Y, Xiao L, et al. Characterizing heat shock protein 90 gene of Apolygus lucorum (Meyer-Dür) and its expression in response to different temperature and pesticide stresses. Cell Stress Chaperone. 2014;19:725–39.

    Article  CAS  Google Scholar 

  15. Xia XF, Lin HL, Zheng DD, Yang G, You MS. Identification and expression patterns of heat shock protein genes in the diamondback moth, Plutella xylostella (Lepidoptera: Yponomeutidae). Acta Entomol Sin. 2013;56:457–64.

    CAS  Google Scholar 

  16. Chen XE, Zhang YL. Identification of multiple small heat-shock protein genes in Plutella xylostella (L.) and their expression profiles in response to abiotic stresses. Cell Stress Chaperone. 2015;20:23–35.

    Article  CAS  Google Scholar 

  17. Nazir A, Saxena DK, Chowdhuri DK. Induction of hsp70 in transgenic Drosophila: biomarker of exposure against phthalimide group of chemicals. Biochim Biophys Acta. 2003;1621:218–25.

    Article  CAS  Google Scholar 

  18. Zhou C, Yang H, Wang Z, Long G, Jin D. Comparative transcriptome analysis of Sogatella furcifera (Horváth) exposed to different insecticides. Sci Rep. 2018;8:8773.

    Article  Google Scholar 

  19. Tene BF, Poupardin R, Costantini C, Awono-Ambene P, Wondji CS, Ranson H, et al. Resistance to DDT in an urban setting: common mechanisms implicated in both M and S forms of Anopheles gambiae in the city of Yaoundé, Cameroon. PLoS ONE. 2013;8:e61408.

    Article  Google Scholar 

  20. Yoshimi T, Minowa K, Karouna-Renier NK, Watanabe C, Sugaya Y, Miura T. Activation of a stress-induced gene by insecticides in the midge, Chironomus yoshimatsui. J Biochem Mol Toxic. 2010;16:10–7.

    Article  Google Scholar 

  21. Caraballo H, King K. Emergency department management of mosquito-borne illness: malaria, dengue, and West Nile virus. Emerg Med Pract. 2014;16:1–23.

    PubMed  Google Scholar 

  22. Coleman M, Hemingway J, Gleave KA, Wiebe A, Gething PW, Moyes CL. Developing global maps of insecticide resistance risk to improve vector control. Malar J. 2017;16:86.

    Article  Google Scholar 

  23. Sleigh AC, Liu XL, Jackson S, Li P, Shang LY. Resurgence of vivax malaria in Henan Province, China. Bull World Health Organ. 1998;76:265–70.

    CAS  PubMed  PubMed Central  Google Scholar 

  24. Zhu G, Zhong D, Cao J, Zhou H, Li J, Liu Y, et al. Transcriptome profiling of pyrethroid resistant and susceptible mosquitoes in the malaria vector, Anopheles sinensis. BMC Genomics. 2014;15:448.

    Article  Google Scholar 

  25. Chang X, Zhong D, Fang Q, Hartsel J, Zhou G, Shi L, et al. Multiple resistances and complex mechanisms of Anopheles sinensis mosquito: a major obstacle to mosquito-borne diseases control and elimination in China. PLoS Neglect Trop Dis. 2014;8:e2889.

    Article  Google Scholar 

  26. Yan ZW, He ZB, Yan ZT, Si FL, Zhou Y, Chen B. Genome-wide and expression-profiling analyses suggest the main cytochrome P450 genes related to pyrethroid resistance in the malaria vector, Anopheles sinensis (Diptera Culicidae). Pest Manag Sci. 2018;74:1810–20.

    Article  CAS  Google Scholar 

  27. Wu XM, Xu BY, Si FL, Li J, Yan ZT, Yan ZW, He X, Chen B. Identification of carboxylesterase genes associated with pyrethroid resistance in the malaria vector Anopheles sinensis (Diptera: Culicidae). Pest Manag Sci. 2018;74:159–69.

    Article  CAS  Google Scholar 

  28. Chen B, Zhang YJ, He Z, Li W, Si F, Tang Y, et al. De novo transcriptome sequencing and sequence analysis of the malaria vector Anopheles sinensis (Diptera: Culicidae). Parasit Vectors. 2014;7:314.

    Article  Google Scholar 

  29. Lobo I. Basic Local Alignment Search Tool (BLAST). J Mol Bio. 2012;215:403–10.

    Google Scholar 

  30. Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, et al. Pfam: the protein families database. Nucleic Acids Res. 2014;42:222–30.

    Article  Google Scholar 

  31. Holub EB. The arms race is ancient history in Arabidopsis, the wildflower. Nat Rev Genet. 2001;2:516–27.

    Article  CAS  Google Scholar 

  32. Thompson JD, Gibson TJ, Higgins DG. Multiple sequence alignment using ClustalW and ClustalX. Curr Protoc Bioinformatics. 2002, Chapter 2:Unit 2.3.

  33. Nicholas KB, Nicholas HBJ, Deerfield DWI. GeneDoc: analysis and visualization of genetic variation. EMBnet News. 1997;14:3.

    Google Scholar 

  34. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.

    Article  CAS  Google Scholar 

  35. WHO. Test procedures for insecticide resistance monitoring in malaria vector mosquitoes. Geneva: World Health Organization; 2013.

    Google Scholar 

  36. Gao Q, Beebe NW, Cooper RD. Molecular identification of the malaria vectors Anopheles anthropophagus and Anopheles sinensis (Diptera: Culicidae) in central China using polymerase chain reaction and appraisal of their position within the Hyrcanus group. J Med Entomol. 2004;41:5–11.

    Article  CAS  Google Scholar 

  37. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–11.

    Article  CAS  Google Scholar 

  38. Cole T, Williams BA, Geo P, Ali M, Gordon K, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511–5.

    Article  Google Scholar 

  39. Emelyanov VV. Phylogenetic relationships of organellar Hsp90 homologs reveal fundamental differences to organellar Hsp70 and Hsp60 evolution. Gene. 2002;299:125–33.

    Article  CAS  Google Scholar 

  40. Porcelli D, Butlin RK, Gaston KJ, Joly D, Snook RR. The environmental genomics of metazoan thermal adaptation. Heredity. 2015;114:502–14.

    Article  CAS  Google Scholar 

  41. Mortazavi A, Williams BA, Mccue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.

    Article  CAS  Google Scholar 

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

    Article  Google Scholar 

  43. World Health Organization. Instruction for determining the susceptibility or resistance of mosquito larvae to insecticides. Geneva: World Health Organization; 1981.

    Google Scholar 

  44. Chen JF, Gao T, Wan SQ, Zhang YH, Yang JK, Yu YB, Wang WD. Genome-wide identification, classification and expression analysis of the HSP gene superfamily in tea plant (Camellia sinensis). Insect J Mol Sci. 2018;19:2633–52.

    Article  Google Scholar 

  45. Wang XR, Wang C, Ban FX, Zhu DT, Liu SS, Wang XW. Genome-wide identification and characterization of HSP gene superfamily in whitefly (Bemisia tabaci) and expression profiling analysis under temperature stress. Insect Sci. 2017;17:1–14.

    Article  CAS  Google Scholar 

  46. Wu S, Huang Z, Rebeca CL, Zhu X, Guo Y, Lin Q, et al. De novo characterization of the pine aphid Cinara pinitabulaeformis Zhang et Zhang transcriptome and analysis of genes relevant to pesticides. PLoS ONE. 2017;12:e0178496.

    Article  Google Scholar 

  47. Li ZW, Li X, Yu QY, Xiang ZH, Kishino H, Zhang Z. The small heat shock protein (sHSP) genes in the silkworm, Bombyx mori, and comparative analysis with other insect sHSP genes. BMC Evol Biol. 2009;9:215.

    Article  Google Scholar 

  48. Chen B, Zhong D, Monteiro A. Comparative genomics and evolution of the HSP90 family of genes across all kingdoms of organisms. BMC Genomics. 2006;7:156.

    Article  Google Scholar 

  49. Gupta RS. Phylogenetic analysis of the 90 kD heat shock family of protein sequences and an examination of the relationship among animals, plants, and fungi species. Mol Biol Evol. 1995;12:1063–73.

    CAS  PubMed  Google Scholar 

  50. Guy CL, Li QB. The organization and evolution of the spinach stress 70 molecular chaperone gene family. Plant Cell. 1998;10:539–56.

    Article  CAS  Google Scholar 

  51. Valpuesta JM, Martı́N-Benito J, Gómez-Puertas P, Carrascosa JL, Willison KR. Structure and function of a protein folding machine: the eukaryotic cytosolic chaperonin CCT. FEBS Lett. 2002;529:11–6.

    Article  CAS  Google Scholar 

  52. Li Y, Bu C, Li T, Wang S, Jiang F, Yi Y, et al. Cloning and analysis of DnaJ family members in the silkworm, Bombyx mori. Gene. 2016;576:88–98.

    Article  CAS  Google Scholar 

  53. Zhang LJ, Wu ZL, Wang KF, Liu Q, Zhuang HM, Wu G. Trade-off between thermal tolerance and insecticide resistance in Plutella xylostella. Ecol Evol. 2015;5:515–30.

    Article  Google Scholar 

  54. Dowling V, Hoarau PC, Romeo M, O’Halloran J, Pelt FV, O’Brien N, et al. Protein carbonylation and heat shock response in Ruditapes decussatus following p, p′-dichlorodiphenyldichloroethylene (DDE) exposure: a proteomic approach reveals that DDE causes oxidative stress. Aquat Toxicol. 2006;77:11–8.

    Article  CAS  Google Scholar 

  55. Škerl MIS, Gregorc A. Heat shock proteins and cell death in situ localisation in hypopharyngeal glands of honeybee (Apis mellifera carnica) workers after imidacloprid or coumaphos treatment. Apidologie. 2010;41:73–86.

    Article  Google Scholar 

  56. Yadav P, Barde PV, Gokhale MD, Vipat V, Mishra AC, Pal JK, et al. Effect of temperature and insecticide stresses on Aedes aegypti larvae and their influence on the susceptibility of mosquitoes to dengue-2 virus. Southeast Asian J Trop Med Public Health. 2005;36:1139–44.

    CAS  PubMed  Google Scholar 

  57. Wang L, Shan D, Zhang Y, Liu X, Sun Y, Zhang Z, et al. Effects of high temperature on life history traits and heat shock protein expression in chlorpyrifos-resistant Laodelphax striatella. Pestic Biochem Phys. 2017;136:64–9.

    Article  CAS  Google Scholar 

  58. Mukhopadhyay I, Siddique HR, Bajpai VK, Saxena DK, Chowdhuri DK. Synthetic pyrethroid cypermethrin induced cellular damage in reproductive tissues of Drosophila melanogaster: Hsp70 as a marker of cellular damage. Arch Environ Cont Tox. 2006;51:673–80.

    Article  CAS  Google Scholar 

  59. Chen J, Kitazumi A, Alpuerto J, Alyokhin A, Reyes BL. Heat-induced mortality and expression of heat shock proteins in Colorado potato beetles treated with imidacloprid. Acta Entomol Sin. 2016;23:548–54.

    CAS  Google Scholar 

  60. Shashikumar S, Rajini PS. Cypermethrin-induced alterations in vital physiological parameters and oxidative balance in Caenorhabditis elegans. Pestic Biochem Phys. 2010;97:235–42.

    Article  CAS  Google Scholar 

  61. Sonoda S, Tsumuki H. Induction of heat shock protein genes by chlorfenapyr in cultured cells of the cabbage armyworm, Mamestra brassicae. Arch Insect Biochem Phys. 2010;65:210–22.

    Article  Google Scholar 

  62. Li G, Zhao H, Zhang X, Zhang Y, Zhao H, Yang X, et al. Environmental stress responses of DnaJA1, DnaJB12 and DnaJC8 in Apis cerana cerana. Front Genet. 2018;9:445–58.

    Article  Google Scholar 

  63. Zhang Y, Liu Y, Zhang J, Guo YP, Ma EB. Molecular cloning and mRNA expression of heat shock protein genes and their response to cadmium stress in the grasshopper Oxya chinensis. PLoS ONE. 2015;10:e0131244.

    Article  Google Scholar 

  64. Wei DD, Chen EH, Ding TB, Chen SC, Dou W, Wang JJ. De novo assembly, gene annotation, and marker discovery in stored-product pest Liposcelis entomophila (Enderlein) using transcriptome sequences. PLoS ONE. 2013;8:e80046.

    Article  Google Scholar 

  65. Feng H, Wang L, Liu Y, He L, Li M, Lu W, et al. Molecular characterization and expression of a heat shock protein gene (HSP90) from the carmine spider mite, Tetranychus cinnabarinus (Boisduval). J Insect Sci. 2010;10:112–26.

    Article  Google Scholar 

  66. Zou J, Guo Y, Guettouche T, Smith DF, Voellmy R. Repression of heat shock transcription factor HSF1 activation by HSP90 (HSP90 complex) that forms a stress-sensitive complex with HSF1. Cell. 1998;94:471–80.

    Article  CAS  Google Scholar 

  67. Chen XE, Zhang Y. Identification of multiple small heat-shock protein genes in Plutella xylostella (L.) and their expression profiles in response to abiotic stresses. Cell Stress Chaperones. 2015;20:23–35.

    Article  CAS  Google Scholar 

  68. Rand EED, Smit S, Beukes M, Apostolides Z, Pirk CWW, Nicolson SW. Detoxification mechanisms of honey bees (Apis mellifera) resulting in tolerance of dietary nicotine. Sci Rep. 2015;5:11779.

    Article  Google Scholar 

  69. Chen B, Piel WH, Gui L, Bruford E, Monteiro A. The HSP90 family of genes in the human genome: insights into their divergence and evolution. Genomics. 2005;86:627–37.

    Article  CAS  Google Scholar 

  70. Marco LD, Sassera D, Epis S, Mastrantonio V, Ferrari M, Ricci I, et al. The choreography of the chemical defensome response to insecticide stress: insights into the Anopheles stephensi transcriptome using RNA-Seq. Sci Rep. 2017;7:41312.

    Article  Google Scholar 

Download references

Authors’ contributions

BC and FLS conceived and designed the study. FLS, BC, LQ, QYH, YZ and ZTY performed the data analysis or experiment, and FLS and BC drafted the manuscript. All authors read and approved the final manuscript.

Acknowledgements

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets generated and/or analysed during the current study are available from the corresponding author on reasonable request.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Funding

This research was supported by the following, the National Natural Science Foundation of China (31672363, 31872262), the National Key Program of Science and Technology Foundation Work of China (2015FY210300), and the Science and Technology Research Program of Chongqing Municipal Education Commission (KJQN201800527).

Publisher’s Note

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

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Bin Chen.

Additional files

Additional file 1: Table S1.

Detailed information and intron–exon organization of An. sinensis HSPs genes.

Additional file 2: Fig. S1.

Gene structure of predicted Anopheles sinensis HSP genes. Blue boxes represent exons and black lines represent introns. The numbers indicate the splicing phases of the HSP transporter genes: 0 refers to phase 0; 1 to phase 1; and 2 to phase 2.

Additional file 3: Fig. S2.

Phylogenetic relationships of HSP genes from four insect species. The phylogenetic tree was constructed using the maximum-likelihood method based on predicted amino acid sequences. Analysis was performed with the program package MEGA5.0. The number at the branch point of the node represents the bootstrap values calculated from 1000 replications, and the gaps were deleted with pairwise deletion method. As, An. sinensis; Ag, An. gambiae; Cq, Cx. quinquefasciatus and Ae, Ae. aegypti. (A) HSP90 family; (B) HSP70 family; (C) Chaperonins family; (D) HSP40 family; (E) sHSP family.

Additional file 4: Table S2.

Type and domain structure of 40 HSP40 (DNAJ) genes identified on An. sinensis genome. All these genes have complete open reading frame (ORF) sequences, with all supported by transcripts identified.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Si, FL., Qiao, L., He, QY. et al. HSP superfamily of genes in the malaria vector Anopheles sinensis: diversity, phylogenetics and association with pyrethroid resistance. Malar J 18, 132 (2019). https://doi.org/10.1186/s12936-019-2770-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-019-2770-6

Keywords