The small heat shock protein (sHSP) genes in the silkworm, Bombyx mori, and comparative analysis with other insect sHSP genes
© Li et al; licensee BioMed Central Ltd. 2009
Received: 20 May 2009
Accepted: 28 August 2009
Published: 28 August 2009
Small heat shock proteins (sHSPs) are products of heat shock response and of other stress responses, and ubiquitous in all three domains of life, archaea, bacteria, and eukarya. They mainly function as molecular chaperones to protect proteins from being denatured in extreme conditions. Study on insect sHSPs could provide some insights into evolution of insects that have adapted to diverse niches in the world.
Taking advantage of the newly assembled genome sequence, we performed a genome-wide analysis of the candidate sHSP genes in the silkworm, Bombyx mori. Based on known silkworm sHSP sequences, we identified 16 silkworm sHSP genes. Most of them are distributed on two silkworm chromosomes 5 and 27, respectively. 15 of 16 silkworm sHSPs have expression evidence. The comparative analysis of insect sHSPs from B. mori, Drosophila melanogaster, Apis mellifera, Tribolium castaneum, and Anopheles gambiae revealed that there is only one orthologous cluster whereas remaining clusters are species-specific on the phylogenetic tree. This suggested that most of sHSPs might have diverged in function across insects investigated. In addition, the data presented in this study also revealed that sHSPs in the insect orthologous cluster are highly conserved in both sequence and expression pattern. In sum, insect sHSPs show a completely different evolutionary pattern from that found in vertebrate sHSPs.
B. mori has the largest number of insect sHSP genes characterized to date, including 16 genes. The inference that most species-specific sHSPs might have diverged in function across insects investigated will help us understand the adaptability of these insects to diverse environments.
It was in Drosophila that researchers observed the phenomenon of heat shock response (HSR) at the first time . Then these proteins highly expressed in thermal stress were defined as heat shock proteins (HSPs). Previous studies suggest that the small heat shock proteins (sHSPs) are abundant and ubiquitous in almost all organisms, from bacteria to algae with single cell to the higher organisms including human [2–5]. HSPs can be divided into five families, including HSP100, HSP90, HSP70, HSP60, and sHSP. This classification is mainly based on the molecular weight as well as the homologous relationship of HSPs. All of these families are conserved, except for sHSP family that is more diverse than other four. However, there are some common characteristics in the function and structure of sHSP. The molecular weights of sHSPs range from 12 kDa to 42 kDa, usually bellow 30 kDa [2, 4]. The sHSPs have an α-crystalling domain comprising about 100 amino acid residues, which is the conserved structure of all sHSP sequences [6–8]. There is a conserved β-sheet sandwich in sHSP secondary structure, and it is these β-sheet structures that help several subunits of sHSP to form a large oligomer [2, 9, 10]. This stable multimeric structure formed by sHSPs has the function of molecular chaperone, which binds to the proteins and prevents them from thermal denaturation [10–12]. The low molecular weight, the conserved secondary structure, and protein domain are principal features used to identify a protein belonging to sHSP family.
In addition to functioning as molecular chaperones to protect proteins from being denatured in high temperature stress [10, 13], sHSPs can also develop the protection function in the conditions of other stresses, such as cold, drought, oxidation, hypertonic stress, UV, and heavy metals [4, 14], even high population density of organisms . The HSPs including sHSPs are playing the part of chaperone function not only in stress conditions, but also in normal development . Although the sHSP family has relatively conserved function of molecular chaperone and the C-terminal of these proteins harbors the conserved α-crystalling domain, the N-terminal sequences of these proteins are variable. This indicates that the conserved C-terminal sequence is a significant part for sustaining the chaperone and other functions of sHSP in cells or organisms whereas the diverse N-terminal sequences may be associated with the diverse expressions, functions, and evolutionary patterns among sHSPs.
There have been numerous studies on HSPs in bacteria, algae, plant, amphibians, birds, and mammalian, especially in the model organisms of Arabidopsis thaliana, Saccharomyces cerevisiae, Caenorhabditis elegans, Drosophila melanogaster, zebra fish, Mus musculus[3–5]. Although the heat shock response was fist discovered in insect , the study of HSP in insects, especially for sHSP, is not as extensive and penetrating as in other organisms. sHSPs primarily have chaperone activity and reflect the response machine of organisms to some extreme stresses existing in environment. Insects are one of the most successful organisms in the world and have a strong ability to adapt to various habitats. So the study of insect sHSP is necessary.
In this study, we first identified the sHSP genes of Bombyx mori based on the new assembly of the silkworm genome sequence http://silkworm.swu.edu.cn/silkdb/. Then we performed comparative analyses with the sHSPs from other four insects whose complete genome sequences are available. It was found that there is an orthologous cluster including five sHSP genes that come from each of five insects. This group of genes is so conserved in sequence among insects, which is in sharp contrast to the other sHSP genes that are species-specific and show the evolutionary pattern of lineage-specific expansion. In addition, the expression patterns of the silkworm sHSPs were investigated by available microarray data. Our data provide some new insights into evolution and functions of the insect sHSPs.
sHSP genes in silkworm and other insects
The sHSP genes in B. mori
Scaffold and Interval
Predicted molecular weight (Da)
Sakano et al.2006
Sakano et al.2006
Sakano et al.2006
Li et al.2005
Sakano et al.2006
Liu et al.2008
Sakano et al.2006
B. mori has 28 chromosome pairs. It can be known from Table 1 that the genes encoding sHSPs were distributed on at least four chromosomes 27, 5, 19, and 13. Because the scaffolds nscaf2803 and scaffold2784 could not be located on the chromosomes in current version of the silkworm genome sequence, the locations of Bm21.6 (in nscaf2803) and Bm20.1 (in scaffold2784) were not determined. It is clear that the distribution of the silkworm sHSP genes in the genome is not homogeneous. There are six sHSP genes located on the silkworm chromosome 27, five genes on chromosome 5, two genes on chromosome 19, and only one gene on chromosome 13. In addition, we found that genes distributed on the chromosomes 27 and 5 have no intron. One gene, Bm22.6 distributed on chromosome 19, contains one intron (the intron length is 2982 bp), and another one, Bm21.4 distributed on chromosome 13, contains two introns (the lengths are 5263 bp and 620 bp, respectively). Compared with the full-length (400~750 bp) of most of the silkworm sHSP genes, these two genes containing introns are relatively too long, especially for Bm21.4. In the previous report , the complete gene sequence of Bm21.4 could not be amplified through PCR method using silkworm genomic DNA as template. Probably, this is due to existence of a relatively big intron. Bm22.6, Bm21.4, and Bm42.3 are different from the other silkworm sHSP genes in chromosome location and the existence of intron. It can be seen from Table 1 that Bm20.8, Bm20.4, Bm19.1, and Bm19.9 are tandem arranged on the middle region of chromosome 27 and there is no any other gene interrupting this genomic region. Bm 23.8 is located on the near end of chromosome 27. However, only two of five sHSP genes on chromosome 5 are tightly tandem arranged (Table 1 and Additional file 1).
The sHSP genes in other four insects
Predicted molecular weight (Da)
Silkworm has the largest number of sHSP genes (16 genes) whereas An. gambiae has the least number of sHSP genes (7 genes) among the insects investigated. With the information of chromosome location and intron number, we can roughly divide the sHSP genes into two types for each of the species. One is that the genes clustered on the same chromosome (two chromosomes for silkworm) tend to have no intron and the other is that the genes (1 to 3 genes) dispersed on chromosomes tend to have intron.
We also performed homology search in the genome of S. cerevisia and found two genes that have high identity with silkworm sHSP genes. One is Scer23.9, located on the chromosome 2. This gene has been reported previously and its NCBI accession number is NP_009628.1. We used it as the outgroup in the phylogenetic analyses.
Phylogenetic analysis of insect sHSP genes
Figure 2 shows that there are two types of clusters on the phylogenetic tree. One type is an orthologous cluster with a very high supporting value measured by the posterior probability, which includes five members that come from each of five insects, respectively. These genes all have at least two introns and dispersed in genomes without any linkage relationships to other sHSP genes. The sHSP genes of D. melanogaster (Dmel20.8) and An. gambiae (Agam20.9) are located on the chromosome X, which is the sex chromosome. The branching pattern of the orthologous cluster shows that Dmel20.8 and Agam20.9 are closer in five orthologous genes and Amel25.6 and Tcas21.8b are closer. Surprisingly, Bm21.4 did not cluster closely to Diptera but was placed on the basal of the cluster, which is not consistent with the result of a recent study on insect molecular systematics . This may imply that Bm21.4 might have experienced a specific evolutionary process.
The other type on the phylogenetic tree is species-specific clusters as a result of lineage-specific expansion. Silkworm-specific cluster includes 12 sHSP genes but has a relative low support value. Most of these genes are located on two chromosomes 5 and 27, however, Bm27.4 and Bm26.6 located on chromosome 5 are not included in this cluster. Bm42.3 is also contained in this cluster.
The posterior probabilities of the species-specific clusters in other four insects are higher (Figure 2). In Ap. mellifera, the species-specific cluster contains 7 sHSP genes and all these genes are distributed on chromosome 2. Additionally, 6 of these genes have no intron except Amel21.3 contains one intron. It is interesting that Amel22.0 located on chromosome 2 and with three introns is not the member of this cluster. In D. melanogaster, sHSP genes distributed on the left arm of chromosome 3 are all included in the species-specific cluster. The gene of Dmel24.5 located on chromosome 2 is also in this cluster. All these genes have no intron. Most of members of the species-specific cluster of T. castaneum sHSP are located on chromosome 8 and only Tcas22.2 is on chromosome 3. They are all single exon genes. Three sHSP genes on the same chromosome of An. gambiae formed a cluster and they all have no intron. In addition, there are other sHSP specific genes in An. gambiae, including Agam21.7, Agam21.6, and Agam23.5b. These three genes have two or three introns and are also located on chromosome 2. All tandem arranged genes of each species except An. gambiae are included in the species-specific cluster (Table 1, Table 2 and Figure 2).
Most strikingly, most members in each of species-specific clusters are usually distributed on the same chromosome in a tandem manner. This gene expansion may be the result of gene duplication events created by unequal crossing over.
The orthologous sHSP genes in insects
We performed blast searching Genbank with these orthologous sHSP sequences to look for orthologous sHSPs in other insects or distant species. We did find that many other insects have this orthologous gene, such as Acyrthosiphon pisum (Accession no. XP_001949402), Aedes aegypti (XP_001657982), Culex quinquefasciatus (XM_001870731), Heliconius erato (ABS57447), Locusta migratoria (ABC84493), and Nasonia vitripennis (XP_001607669). Furthermore, an arthropods of Ixodes scapularis (EEC06453) belonging to Arachnoidea and four nematoda organisms of Caenorhabditis briggsae (XP_001676833), C. elegans (NP_001024374), Brugia malayi (U48407) and Dirofilaria immitis (AAB08736) also contain the orthologous sHSP genes which are slightly divergent from the insect orthologs. Compared with the five orthologous sHSP genes in Figure 2, the homologous genes mentioned above have not only high identity but also a conserved N-terminal, which is the principal structure characteristic of the genes in the orthologous cluster (Additional file 3). So, this kind of orthologous genes exists not only in insects but also in other invertebrates.
To confirm whether the orthologous group genes exist in vertebrate, we reconstructed a larger phylogenetic tree including sHSP genes from above five insects, nematode, amphioxus, mouse and human (Additional file 4). The sequences of sHSP genes in nematode, mouse and human came from previous studies [5, 21, 22]. 14 amphioxus sHSP genes were obtained by searching the genome of Branchiostoma floridae  using blast method. The phylogenetic analysis results indicated that, as expectation, a nematode sHSP gene Ce25 was included in the highly supported orthologous cluster but the human and mouse sHSP genes were not in it. However, the HsapHspB11 from Human and MmusHSPB11 from mouse as well as two sHSP genes from amphioxus formed a cluster although the support value is low. Furthermore, this cluster seems to be close to the orthologous group (Additional file 4). Therefore, it appears that a kind of highly conserved sHSP genes may exist even in vertebrates. In addition, we observed that there are many orthologous sHSP clusters between human and mouse, which is consistent with previous studies [21, 22].
The species specific sHSP genes in insects, especially in the silkworm
Our study revealed that most insect sHSP genes are species-specific (Figure 2). The silkworm contains the most number of the species-specific sHSPs. They all have a conserved α-crystalling domain and a diversified N-terminal region. Whether there are divergent expression patterns in these paralogous genes is worthy to be detected for this domesticated insect, B. mori.
We searched for EST evidences of sHSP genes through SilkDB . 14 silkworm sHSP genes have complete EST sequences except Bm21.6 and Bm42.3. But Bm21.6 had been cloned in previous studies, as well as the genes of Bm19.9 [17, 18], Bm23.8, Bm20.8, Bm20.4, Bm21.4, and Bm20.1 , Bm24.2 , Bm19.1, Bm15.7, Bm 22.6, Bm27.4, Bm20.2, and Bm19.5 (not published). So, 15 silkworm sHSP genes are certainly transcribed. In fact, there is no EST sequence associated with Bm42.3. Several possible reasons may explain the absence of EST evidence for Bm42.3. One is that this candidate silkworm sHSP gene is not expressed in usual conditions or its expression is too weak to be detected. Another is that it might have lost its regulatory elements. All of these should be confirmed in future.
In this study, we first identified 16 silkworm sHSP genes based on the new assembly of the B. mori genome sequence , which is the largest number of the sHSPs of five insects in different orders that whole genome sequences are available. Previous study based on six silkworm sHSP genes suggested that the sHSP genes in the silkworm may be classified into two groups . Similarly, the comparative analysis of sHSPs in five insects revealed that there are two types of clusters of insect sHSPs: one is orthologous and the other is species-specific. Interestingly, there is one orthologous cluster that includes one sHSP from each insect whereas most insect sHSPs are species-specific (Figure 2). These results suggested that most insect sHSPs might have experienced lineage-specific expansion and diverged in function at least amongst insect orders. Nevertheless, a few insect sHSPs such as those in orthologous cluster are highly conserved in whole sequence as well as in gene structure, suggesting that they may also have highly conserved function. The phylogenetic relation of the orthologous cluster may be interpreted as the speciation history.
Most members of species-specific cluster are tandem arranged genes in each insect and there is no orthologous relationship between tandem genes of different insects. This case may occur if a primitive tandem set with several genes existed in the ancient of the five insects, subsequently, frequent and general gene conversion events homogenize the tandem sHSP genes in each insect after the speciation, finally the originally orthologous sHSP genes became species-specific. So we estimated the gene conversion events of sHSP genes in each insect with GeneConv program http://www.math.wustl.edu/~sawyer/geneconv/. Surprisingly, although the P values are below 0.01, the predicted gene conversion nucleotide residues are no longer than 40 bp in our tests (most are below 25 bp) except a gene pair in the silkworm of Bm20.8 and Bm 20.4 (P value, 0.001; length, 128 bp). These results imply that gene conversion is not the major mechanism explaining the topology in Figure 2. Therefore, our observations do not support the assumption above. Members of species-specific cluster or tandem arranged genes might be the result of independent gene expansion occurred during evolution of each insect.
The gene expansion process may associate with gene duplication and loss events, which are known as a birth-and-death model. This model was successfully used to interpret some evolutionary process of most gene families including HSP70 and cytochrome P450 [27, 28]. If this is the case, Figure 2 suggests that some gene duplication events of sHSPs occurred before the radiation of insects, and that one or two copies experienced lineage-specific expansions through successive gene duplications after the radiation of insects.
The tandem arranged pattern of sHSP genes is not exclusive in insects. For example, sHSP16 genes in C. elegans and C. briggsae also show similar pattern . The tandem repeats of sHSP genes may facilitate organisms to rapidly respond to changing environmental conditions because of regulation reason. All the five insects have independently developed this strategy for sHSP genes, which is similar to the convergent duplication of the proneural genes in insects . However, the 10 sHSP genes in human are separately distributed on eight chromosomes . The different arranged forms of sHSP genes may reflect different evolutionary processes that organisms experienced. Tandem sHSP genes may be a better way for insects to regulate gene expressions in diverse environments.
It is known that the sHSPs are members of HSP family with function of molecular chaperone. They are ubiquitous in all organisms and extremely conserved in the sequence structure of α-crystalling domain resided in C-terminal. The conserved C-terminal region might be the result of strongly functional constraint of sHSPs as a molecular chaperone. sHSPs bind to other cellular proteins in thermal and other extremely devastating stresses, and protect them from denaturation. They also have functions in protein folding and transportation, the embryo development and immunization mechanism of organisms. The conservative degrees of sHSP genes are heterogeneous in different protein regions. In general, the N-terminal coding sequences are more variable than C-terminal with the α-crystalling domain. However, the sHSPs in the orthologous cluster showed high conservation along whole sequence. The orthologous sHSP gene (Bm21.4) in silkworm was cloned and found to be not induced by thermic stress (45°C), and it was expressed constitutively in normal conditions (room temperature) in fat body and other tissues of silkworm . The microarray data also revealed the similar expression patterns (the probe number is sw01265). So Bm21.4 is special in sequence and expression compared with other silkworm sHSP genes . These observations suggest that insect orthologous sHSPs may not be associated with response of insects to environmental stresses. In contrast, they may be involved in basically metabolic processes in insects.
One previous study based on the limited number of sHSP sequences suggested that the phylogenetic relationship of sHSP genes in insects is quite different from that of other species or taxonomical groups . Mammals have at least ten sHSP subfamilies and they also have orthologous genes in zebrafish . The identified sHSP genes in amphibians are also orthologous to those of mammals . In general, there are multiple orthologous sHSP clusters in vertebrate and angiosperm sHSPs also have similar pattern [4, 32]. However, insect sHSP genes principally show one orthologous cluster that is composed of each sHSP from each species investigated, and most of the genes are species-specific (Figure 2). Thus, insect sHSPs display completely different evolutionary pattern from that of vertebrates or plants.
To determine whether the evolutionary pattern observed in five insects is universal for other arthropods and nematode, we searched the GenBank for the orthologous sHSPs in corresponding species. Indeed, we found orthologous sHSP genes in many other arthropods and nematode (Additional file 3). Most of these species are insects, such as L. migratoria, He. erato , Ac. pisum, and Ae. aegypti. These insects cover the orders Lepidoptera, Hymenoptera, Diptera, Hemiptera, Coleoptera, and Orthoptera. This ubiquitous existence in insects suggests that all of the species of insecta may possess one copy of this conserved gene in their genomes. There is a species of Arachnida, I. scapularis, that contains this orthologous gene. This gene is also detected within each of nine nematode species . A species of chordate, Branchiostoma lanceolatum, also contains this orthologous sHSP gene named B1Hsp24.1 (GenBank accession number: AJ616737) with relatively lower identity . Insecta and Arachnida belong to the phylum Arthropoda. Furthermore, Arthopods, Branchiostoma, and phylum Nematoda are all invertebrates. So, these orthologous sHSP genes do exist in Arthropods and Nematode and probably in all invertebrates.
That insects harbor orthologous sHSP genes is interesting. To know how these gens are related to the vertebrate sHSPs, we performed the phylogenetic analysis on the sHSPs from five insects, nematode, amphioxus, mouse and human. Although no highly supported orthologous cluster was found on the phylogenetic tree, it seems that vertebrates also have the highly conserved sHSPs that are close to insect orthologs. Therefore, it is important to understand functions of these conserved sHSPs of different organisms in future.
The orthologous sHSP genes have not only the conserved sequences but also the conserved gene structure across insect orders probably and the sequences of α-crystalling domains in these genes are more conserved than that in species-specific genes, all of which is the result of strongly purifying selection during evolution. P27, another orthologous gene of Bm21.4 in Di. immitis, is also constitutively expressed in male and female adult worms and developmental phase such as L3–L4 molt and is located in hypodermal tissues. Additionally, it is not up-regulated by heat shock at 43°C . In C. elegans, Ce25 is also constitutively expressed in normal conditions and this protein is principally located in body wall muscle specially bonding to vinculin and α-actinin . However, the orthologous gene in migratory locust, hsp20.6, can regulate its expression associated with the population density, for example, high density corresponding to an increased expression . We infer that these orthologous sHSP genes may be constitutively expressed in organisms and may not be involved in the heat shock response. It is likely that they are house-keeping genes and important in the development, adaptability and metabolism of organisms.
Another important feature of insect sHSPs is that most of them are species-specific in five insects investigated. This suggested that functions of most sHSPs might have diverged among the insects, reflecting that the sHSPs might have contributed to adaptability of insects to different niches. Silkworm has the largest cluster of species-specific sHSPs which includes 12 genes although the cluster has a relatively low supporting value measured by posterior probability. 12 silkworm species-specific sHSPgenes diverged in sequences especially in N-terminal region (Figure 1). There are several previous studies focused on these proteins [17–19, 36–38]. Silkworm sHSP genes Bm20.4, Bm20.8, Bm19.9, Bm20.1, Bm24.2, Bm23.8, and Bm19.5 have increased expression by heat shock stress [17–19, 38]. It is the sHSP genes belonging to species-specific cluster that principally participate in heat shock response of silkworm and other insects. However, these genes are also expressed in usual conditions without extreme stress and have distinct expression levels in different tissues of silkworm (Figure 4). The expression of Bm19.9 is significantly high in testis and ovary , and Bm20.1 is highly expressed in testis of the fifth instar . There are other sHSP genes having high expression intensities in gonad through the microarray data (Figure 4). Silkworm sHSP genes highly expressed in gonad may play an important role in the development of germocyte. Bm19.1, Bm15.7, and Bm23.8 are increasedly expressed in midgut. Furthermore, Bm22.6, which does not belong to the species specific cluster, is also highly expressed in midgut (Figure 4). Additionally, Bm22.6 is one of the 11 genes significantly up-regulated in midgut following the infection of BmDNV-Z, a densonucleosis virus . Bm22.6 and Bm19.1 are also remarkably expressed in integument (Figure 4). Integument and midgut are tissues with immune function. Previous studies suggest that HSPs may participate in the innate immune responses of silkworm . So the silkworm sHSP genes highly expressed in midgut or integument may have functions in immune defense mechanism.
There are some changes in cellular physiology before and after the heat stress. Bm20.1 principally located in nucleus and Bm20.8 in cytoplasm, are all transferred to the membrane after heat shock, and then recovered to the previous subcellular localization in normal conditions. This implies that sHSP have functions in the membrane structure and metabolism in extremely conditions [41, 42]. Silkworm sHSP genes also have been studied in embryonic development of diapause and non-diapause eggs , silkworm parthenogenesis , and different silkworm breeds of bivoltine and multivoltine . The studies in heat stress induced pattern, expression diversity in tissues, and other functions without chaperone of silkworm sHSP genes suggest that there are divergences in the biological functions among these genes. These species-specific sHSP genes may resist harmful stimulus for organisms and participate in the reproduction, development, and the normal metabolism activities. It is believed that the diversity and divergence of sHSP genes might have facilitated the insects adapting to various environments.
The data presented in this study revealed that the majority of insect sHSP genes in each species are located on one chromosome except for that most silkworm sHSP genes are located on two chromosomes 5 and 27. It is possible that the distribution of most silkworm sHSP genes on chromosomes 5 and 27 might be the result of split of an ancient chromosome. However, there is opposed evidence against this scenario, which is that silkworm chromosomes 5 and 28 correspond to one chromosome in He. erato and silkworm chromosomes 23 and 27 correspond to another one . So, the expanding process of sHSP genes in silkworm is complicated. The silkworm has the most number of sHSP genes, and An. gambiae has the least. Whether there is relationship between the number of the sHSP genes and the resistance to thermal or other stresses is unknown. The multivoltine silkworms are more resistant than bivoltine in thermal stress, and the expressions of some sHSP genes in multivoltine breeds are lower compared with the bivoltine breeds after heat shock at 45°C . The gregarious locust expresses more sHSP genes than solitary locust due to different population density, and the former have a smaller body size and a weaker fecundity . Therefore, the sensitivity to heat and other stresses may influence the development and reproduction of various voltinism breeds in silkworm.
Based on known silkworm sHSP sequences and the newly assembled genome sequence, we identified 16 silkworm sHSP genes. 15 of 16 silkworm sHSPs have expression evidence. The comparative analysis with other insect sHSPs revealed that there is only one orthologous cluster whereas remaining clusters are species-specific. Our data demonstrated that the evolutionary pattern of insect sHSPs which is different from that found in vertebrates and will provide some new insight into evolution of insect sHSPs. It is believed that studies on the evolution, function and mechanism of stress response genes in insects will help us understand the adaptability of insects to diverse environments.
Identification of insect sHSP genes
We first retrieved 6 known sHSP sequences of B. mori from Genbank  and then used these sequences as queries to search the new assembly of the silkworm genome sequence in SilkDB http://silkworm.swu.edu.cn/silkdb for a complete list of the silkworm sHSPs. The E value for evaluating all sequences in homology search is 10-4. Finally, we obtained 16 sequences as candidate genes for the silkworm sHSPs. The other four insect annotated gene and protein sequences were downloaded from NCBI website http://www.ncbi.nlm.nih.gov/. The four insects are D. melanogaster, Ap. mellifera, T. castaneum, and An. gambiae. We used above 16 sHSP sequences of B. mori as queries to perform the same homology search for other insect sHSPs as above.
One sHSP gene sequences of S. cerevisiae was used as the outgroup. This sequence was retrieved from Genbank (Accession no. NP_009628).
The secondary structure and structural domains and chromosome locations of the insect sHSPs
The information of protein secondary structure and domain has advantage for helping researchers to understand the function of these proteins. So the structure information of B. mori candidate sHSP genes was predicted by PredictProtein http://www.predictprotein.org/, one online service tool. The prediction of protein domain was performed on InterProScan online service http://www.ebi.ac.uk/Tools/InterProScan/.
The location information of B. mori sHSP genes on the chromosomes was obtained by SilkMap online tool http://silkworm.swu.edu.cn/silksoft/silkmap.html. The location information of the other four insect sHSP genes was determined by Map View tool http://www.ncbi.nlm.nih.gov/mapview/ at NCBI website.
The amino acid sequences of sHSP genes of five insects were aligned by three methods, clustal X v.1.8, clustal W v.1.83 , and muscle v.3.6 . The resulting alignments were almost identical and one of them was selected for further analyses. The phylogenetic tree for the insect sHSPs was first reconstructed by Neighbor-Joining (NJ) method in which distance was estimated by JTT amino acid matrix implemented in the program MEGA 4 . The pairwise deletion option was used in the NJ tree reconstruction and the accuracy of the tree topology was assessed by bootstrap analysis with 1,000 resampling replicates.
Using the amino acid sequences of the conserved α-crystalling domain with 100 residues, we also reconstructed the phylogenetic relationship of the insect sHSP genes by the program MrBayes v.3.1.2 . MrBayes is a program of reconstructing phylogenetic relationships based on Bayesian and maximum likelihood methods. This program simulates the posterior distribution of the phylogenetic trees through the method of Metropolis-coupled Markov Chain Monte Carlo . WAG model was demonstrated to be the best model used in the MrBayes program to reconstruct the phylogenetic tree for the amino acid sequences of insect sHSPs. After MrBayes program with WAG model run 6500000 generations, the value of average standard deviation of split frequencies became stable below 0.01 and finally reached to 0.007103 at the end of the running. So, the early 3500000 generations data were seen as the burnin, which were deleted in the analysis of the results. All the trees are rooted by the sHSP gene sequence of S. cerevisiae, scer23.9.
Expression information based on microarray database
B. mori microarray database http://silkworm.swu.edu.cn/microarray/ is composed of expression data for 22987 probes of 9 silkworm tissues on day 3 of the fifth instar in silkworm . We used blastn search to obtain expression information of each silkworm sHSP gene in several tissues at normal or some stress conditions. We also searched the EST data for the silkworm sHSP genes in SilkDB to confirm whether these genes identified from the genome search approach are expressed in the organism or not http://silkworm.swu.edu.cn/silkdb.
We thank one anonymous reviewer for constructive comments that improved the manuscript. This study was supported by the Hi-Tech Research and Development (863) Program of China (2006AA10A117), a grant from National Science Foundation of China (No. 30671587), and the Programme of Introducing Talents of Discipline to Universities (B07045).
- Ritossa F: A new puffing pattern induced by a temperature shock and DNP in Drosophila. Experientien. 1962, 18: 571-573.View ArticleGoogle Scholar
- Kim KK, Kim R, Kim SH: Crystal structure of a small heat-shock protein. Nature. 1998, 394 (6693): 595-599.View ArticlePubMedGoogle Scholar
- Waters ER, Rioflorido I: Evolutionary analysis of the small heat shock proteins in five complete algal genomes. J Mol Evol. 2007, 65 (2): 162-174.View ArticlePubMedGoogle Scholar
- Waters ER, Aevermann BD, Sanders-Reed Z: Comparative analysis of the small heat shock proteins in three angiosperm genomes. Cell Stress Chaperones. 2008, 13 (2): 127-142.PubMed CentralView ArticlePubMedGoogle Scholar
- Aevermann BD, Waters ER: A comparative genomic analysis of the small heat shock proteins in Caenorhabditis elegans and briggsae. Genetica. 2008, 133 (3): 307-319.View ArticlePubMedGoogle Scholar
- Caspers GJ, Leunissen JAM, de Jong WW: The expanding small heat-shock protein family, and structure predictions of the conserved "alpha-crystallin domain". J Mol Evol. 1995, 40 (3): 238-248.View ArticlePubMedGoogle Scholar
- de Jong WW, Caspers GJ, Leunissen JA: Genealogy of the alpha-crystallin-small heat-shock protein superfamily. Int J Biol Macromol. 1998, 22 (3–4): 151-162.View ArticlePubMedGoogle Scholar
- Fu X, Jiao W, Chang Z: Phylogenetic and biochemical studies reveal a potential evolutionary origin of small heat shock proteins of animals from bacterial class A. J Mol Evol. 2006, 62 (3): 257-266.View ArticlePubMedGoogle Scholar
- van Montfort RLM, Basha E, Friedrich KL, Slingsby C, Vierling E: Crystal structure and assembly of a eukaryotic small heat shock protein. Nat Struct Biol. 2001, 8 (21): 1025-1030.View ArticlePubMedGoogle Scholar
- van Montfort RLM, Slingsby C, Vierling E: Structure and function of the small heat shock protein/alpha-crystallin family of molecular chaperones. Adv Protein Chem. 2001, 59: 105-156.View ArticlePubMedGoogle Scholar
- Lee GJ, Pokala N, Vierling E: Structure and in vitro molecular chaperone activity of cytosolic small heat shock proteins from pea. J Biol Chem. 1995, 270 (18): 10432-10438.View ArticlePubMedGoogle Scholar
- Nakamoto H, Vigh L: The small heat shock proteins and their clients. Cell Mol Life Sci. 2007, 64 (3): 294-306.View ArticlePubMedGoogle Scholar
- Sun Y, MacRae TH: Small heat shock proteins: molecular structure and chaperone function. Cell Mol Life Sci. 2005, 62 (21): 2460-2476.View ArticlePubMedGoogle Scholar
- Dasgupta S, Hohman TC, Carper D: Hypertonic stress induces αB-crystallin expression. Exp Eye Res. 1992, 54 (3): 461-470.View ArticlePubMedGoogle Scholar
- Wang HS, Wang XH, Zhou CS, Huang LH, Zhang SF, Guo W, Kang L: cDNA cloning of heat shock proteins and their expression in the two phases of the migratory locust. Insect Mol Biol. 2007, 16 (2): 207-219.View ArticlePubMedGoogle Scholar
- The International Silkworm Genome Consortium: The genome of a lepidopteran model insect, the silkworm Bombyx mori. Insect Biochem Mol Biol. 2008, 38 (12): 1036-1045.View ArticleGoogle Scholar
- Li B, Xia QY, Fujii H, Banno Y, Lu C: Expression of the small heat-shock protein Bmhspl9.9 gene in silkworm (Bombyx mori). J Agr Biotechnol. 2005, 13 (2): 195-201. (in Chinese with English abstract)Google Scholar
- Sakano D, Li B, Xia QY, Yamamoto K, Fujii H, Aso Y: Genes encoding small heat shock proteins of the silkworm, Bombyx mori. Bios Biotech Bioch. 2006, 70 (10): 2443-2450.View ArticleGoogle Scholar
- Liu BB, Chai CL, Wang XY, Huang FF, Lu C: Cloning and functional study of Bombyx mori hsp24 3 gene. Scientia Agricultura Sinica. 2008, 41 (12): 4201-4208. (in Chinese with English abstract)Google Scholar
- Savard J, Tautz D, Richards S, Weinstock GM, Gibbs RA, Werren JH, Tettelin H, Lercher MJ: Phylogenomic analysis reveals bees and wasps (Hymenoptera) at the base of the radiation of Holometabolous insects. Genome Res. 2006, 16 (11): 1334-1338.PubMed CentralView ArticlePubMedGoogle Scholar
- Kappé G, Franck E, Verschuure P, Wilbert C, Boelens WC, Leunissen JAM, de Jong WW: The human genome encodes 10 α-crystallin-related small heat shock proteins:HspB1–10. Cell Stress Chaperones. 2003, 8 (1): 53-61.PubMed CentralView ArticlePubMedGoogle Scholar
- Kampinga HH, Hageman J, Vos MJ, Kubota H, Tanguay RM, Bruford EA, Cheetham ME, Chen B, Hightower LE: Guidelines for the nomenclature of the human heat shock proteins. Cell Stress Chaperones. 2009, 14 (1): 105-111.PubMed CentralView ArticlePubMedGoogle Scholar
- Putnam NH, Butts T, Ferrier DE, Furlong RF, Hellsten U, Kawashima T, Robinson-Rechavi M, Shoguchi E, Terry A, Yu JK, Benito-Gutierrez EL, Dubchak I, Garcia-Fernandez J, Gibson-Brown JJ, Grigoriev IV, Horton AC, de Jong PJ, Jurka J, Kapitonov VV, Kohara Y, Kuroki Y, Lindquist E, Lucas S, Osoegawa K, Pennacchio LA, Salamov AA, Satou Y, Sauka-Spengler T, Schmutz J, Shin-I T, Toyoda A, Bronner-Fraser M, Fujiyama A, Holland LZ, Holland PW, Satoh N, Rokhsar DS: The amphioxus genome and the evolution of the chordate karyotype. Nature. 2008, 453 (7198): 1064-1071.View ArticlePubMedGoogle Scholar
- Wang J, Xia QY, He XM, Dai MT, Ruan J, Chen J, Yu G, Yuan HF, Hu YF, Li RQ, Feng T, Ye C, Lu C, Wang J, Li SG, Wong GKS, Yang HM, Wang J, Xiang ZH, Zhou ZY, Yu J: SilkDB: a knowledgebase for silkworm biology and genomics. Nucl Acids Res. 2005, 33: D339-D402.Google Scholar
- Xia QY, Cheng DJ, Duan J, Wang GH, Cheng TC, Zha XF, Liu C, Zhao P, Dai FY, Zhang Z, He NJ, Zhang L, Xiang ZH: Microarray-based gene expression profiles in multiple tissues of the domesticated silkworm, Bombyx mori. Genome Biol. 2007, 8 (8): R162-PubMed CentralView ArticlePubMedGoogle Scholar
- Sawyer SA: Statistical tests for detecting gene conversion. Mol Biol Evol. 1989, 6 (5): 526-538.PubMedGoogle Scholar
- Nikolaidis N, Nei M: Concerted and nonconcerted evolution of the Hsp70 gene superfamily in two sibling species of nematodes. Mol Biol Evol. 2004, 21 (3): 498-505.View ArticlePubMedGoogle Scholar
- Thomas JH: Rapid birth-death evolution specific to xenobiotic cytochrome P450 genes in vertebrates. PLoS Genet. 2007, 3 (5): 720-728.View ArticleGoogle Scholar
- Negre B, Simpson P: Evolution of the achaete-scute complex in insects: convergent duplication of proneural genes. Trends Genet. 2009, 25 (4): 147-152.View ArticlePubMedGoogle Scholar
- Elicker KS, Hutson LD: Genome-wide analysis and expression profiling of the small heat shock proteins in zebrafish. Gene. 2007, 403 (1–2): 60-69.PubMed CentralView ArticlePubMedGoogle Scholar
- Tuttle AM, Gauley J, Chan N, Heikkila JJ: Analysis of the expression and function of the small heat shock protein gene, hsp27, in Xenopus laevis embryos. Comp Biochem Physiol A Mol Integr Physiol. 2007, 147 (1): 112-21.View ArticlePubMedGoogle Scholar
- Franck E, Madsen O, van Rheede T, Ricard G, Huynen MA, de JongW: Evolutionary diversity of vertebrate small heat shock proteins. J Mol Evol. 2004, 59 (6): 792-805.View ArticlePubMedGoogle Scholar
- Pringle EG, Baxter SW, Webster CL, Papanicolaou A, Lee SF, Jiggins CD: Synteny and chromosome evolution in the Lepidoptera: Evidence from mapping in Heliconius melpomene. Genetics. 2007, 177 (1): 417-426.PubMed CentralView ArticlePubMedGoogle Scholar
- Lillibridge CD, Rudin W, Philipp MT: Dirofilaria immitis: Ultrastructural localization, molecular characterization, and analysis of the expression of p27, a small heat shock protein homolog of Nematodes. Exp Parasitol. 1996, 83 (1): 30-45.View ArticlePubMedGoogle Scholar
- Ding L, Candido EPM: HSP25, a small heat shock protein associated with dense bodies and M-lines of body wall muscle in Caenorhabditis elegans. J Biol Chem. 2000, 275 (13): 9510-9517.View ArticlePubMedGoogle Scholar
- Sun WZ, Li B, Wang YW, Chai CL, Liu BB, Lu C: Expression in vitro and fluorescence in situ hybridization of the HSP20. 8 Gene of Bombyx mori. Science of Sericulture. 2007, 33 (1): 14-19.Google Scholar
- Zhao EW, Sheng Q, Lv ZB, Chen J, Nie ZM, Wang D, Liu LL, Shen HD, Shu JH, Chen JQ, Wu XF, Zhang RZ: Expression and functional research of BmHSP20.8 from silkworm, Bombyx mori. Sciencepaper Online. 2009, (in Chinese with English abstract), [http://www.paper.edu.cn/paper.php?serial_number=200902-314]Google Scholar
- Mao J, Lv ZB, Chen J, Zhang YZ: Analysis of the expression and function of the small heat shock protein HSP20.1 in Bombyx mori. Sciencepaper Online. 2009, (in Chinese with English abstract), [http://www.paper.edu.cn/paper.php?serial_number=200902-546]Google Scholar
- Bao YY, Li MW, Zhao YP, Ge JQ, Wang CS, Huang YP, Zhang CX: Differentially expressed genes in resistant and susceptible Bombyx mori strains infected with a densonucleosis virus. Insect Biochem Mol Biol. 2008, 38 (9): 853-861.View ArticlePubMedGoogle Scholar
- Song KH, Jung SJ, Seo YR, Kang SW, Han SS: Identification of up-regulated proteins in the hemolymph of immunized Bombyx mori larvae. Comp Biochem Phys Part D. 2006, 1 (2): 260-266.Google Scholar
- Tsvetkova NM, Horváth I, Török Z, Wolkers WF, Balogi Z, Shigapova N, Crowe LM, Tablin F, Vierling E, Crowe JH, Vigh L: Small heat-shock proteins regulate membrane lipid polymorphism. Proc Natl Acad Sci USA. 2002, 99 (21): 13504-13509.PubMed CentralView ArticlePubMedGoogle Scholar
- Horváth I, Multhoff G, Sonnleitner A, Vá L: Membrane-associated stress proteins: More than simply chaperones. BBA-Biomembranes. 2008, 1778 (7–8): 1653-1664.View ArticlePubMedGoogle Scholar
- Saravanakumar R, Ponnuvel KM, Qadri SMH: Expression of metabolic enzyme genes and heat-shock protein genes during embryonic development in diapause and non-diapause egg of multivoltine silkworm Bombyx mori. Biol. 2008, 63 (5): 737-744.Google Scholar
- Xia JY, Long XH, Chen J, Nie ZM, Wang D, Lu ZB, Xu J, He PA, Zhang YZ: Bioinformatics analysis of heat shock related proteins from differential protein spots associated with parthenogenesis in Bombyx mori. Science of Sericulture. 2007, 33 (4): 568-573. (in Chinese with English abstract)Google Scholar
- Hosseini Moghaddam SH, Du X, Li J, Cao JR, Zhong BX, Chen YY: Proteome analysis on differentially expressed proteins of the fat body of two silkworm breeds, Bombyx mori, exposed to heat shock exposure. Biotechnol Bioproc Eng. 2008, 13: 624-631.View ArticleGoogle Scholar
- Higgins DG, Thompson JD, Gibson TJ: Using CLUSTAL for multiple sequence alignments. Methods Enzymol. 1996, 226: 383-402.View ArticleGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucl Acids Res. 2004, 32 (5): 1792-1797.PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.