Next Article in Journal
Plastic Fly: What Drosophila melanogaster Can Tell Us about the Biological Effects and the Carcinogenic Potential of Nanopolystyrene
Previous Article in Journal
Effect of Lanthanide Ions and Triazole Ligands on the Molecular Properties, Spectroscopy and Pharmacological Activity
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 25
Issue 14
10.3390/ijms25147963
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Identification and Interaction Analysis of Turbot Heat Shock Protein 40 and 70 Families Suggest the Mechanism of Chaperone Proteins Involved in Immune Response after Bacterial Infection

by
Yuanwei Geng
1,†,
Yuxuan Gai
1,2,3,†,
Yanping Zhang
4,
Shengwei Zhao
1,
Anlan Jiang
1,
Xueqing Li
1,
Kaiqing Deng
1,
Fuxuan Zhang
1,
Lingling Tan
1,* and
Lin Song
1,2,3,*
1
School of Life Science, Qingdao Agricultural University, Qingdao 266109, China
2
Key Laboratory of Applied Mycology, Qingdao Agricultural University, Qingdao 266109, China
3
Qingdao International Center on Microbes Utilizing Biogas, Qingdao Agricultural University, Qingdao 266109, China
4
College of Entrepreneurship and Innovation, Qingdao Agricultural University, Qingdao 266109, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2024, 25(14), 7963; https://doi.org/10.3390/ijms25147963
Submission received: 3 July 2024 / Revised: 17 July 2024 / Accepted: 19 July 2024 / Published: 21 July 2024
(This article belongs to the Section Biochemistry)
Browse Figures
Versions Notes

Abstract

:
Hsp40–Hsp70 typically function in concert as molecular chaperones, and their roles in post-infection immune responses are increasingly recognized. However, in the economically important fish species Scophthalmus maximus (turbot), there is still a lack in the systematic identification, interaction models, and binding site analysis of these proteins. Herein, 62 Hsp40 genes and 16 Hsp70 genes were identified in the turbot at a genome-wide level and were unevenly distributed on 22 chromosomes through chromosomal distribution analysis. Phylogenetic and syntenic analysis provided strong evidence in supporting the orthologies and paralogies of these HSPs. Protein–protein interaction and expression analysis was conducted to predict the expression profile after challenging with Aeromonas salmonicida. dnajb1b and hspa1a were found to have a co-expression trend under infection stresses. Molecular docking was performed using Auto-Dock Tool and PyMOL for this pair of chaperone proteins. It was discovered that in addition to the interaction sites in the J domain, the carboxyl-terminal domain of Hsp40 also plays a crucial role in its interaction with Hsp70. This is important for the mechanistic understanding of the Hsp40–Hsp70 chaperone system, providing a theoretical basis for turbot disease resistance breeding, and effective value for the prevention of certain diseases in turbot.

1. Introduction

The turbot, Scophthalmus maximus, is a significant economic species in aquaculture [1,2]. Stressors during turbot cultivation can be categorized into abiotic factors, such as high temperature, salinity, and water quality, and biotic factors, including bacterial, viral, and parasitic diseases, with bacterial infections being particularly common. Recent years have seen an increase in outbreaks of infections caused by Aeromonas salmonicida, emerging as a major disease threat [3]. Heat shock proteins (HSPs) have long been considered molecular chaperones that assist in the proper folding of proteins under high temperature stress [4,5,6]; however, increasing research has highlighted the pivotal role of HSPs in post-infection immunity [7,8,9,10,11]. Hsp40 and Hsp70, as a pair of molecular chaperones, typically interact structurally to function effectively [12]. Nevertheless, in the immune response to infection in turbot, the focus has predominantly been on the function of Hsp70, neglecting the interaction between this chaperone pair. This study will systematically identify the Hsp40 and Hsp70 families in turbot and analyze their interaction following A. salmonicida infection. This study aims to fill the gap in our understanding of the Hsp40 and Hsp70 interaction, offering a deeper comprehension of their mechanistic roles in fighting infections and managing stress in aquaculture environments.
Members of the Hsp40 family are also known as J proteins (DnaJ) as they all contain a highly conserved J domain [13]. Hsp40 can be classified into three groups: group A possesses the J domain, the Gly/Phe-rich region, and the cysteine repeats; group B possess the J domain and the Gly/Phe-rich region, but lack the cysteine repeats; and group C only possesses a helical J domain [5,14,15]. The J domain is a key site for the interaction of Hsp40 and Hsp70, which contains four α-helical structures. The second helix is enriched with lysines, a region rich in positively charged amino acid residues that are necessary for the ATPase activity of Hsp70. There is a tripeptide between the second and third helix consisting of histidine, proline, and aspartic acid (HPD) that constitutes its core structure [16,17]. It was recently found that E. coli DnaJ (Hsp40) cannot activate ATP hydrolysis of DnaK (Hsp70) by the J domain alone and that the ATPase stimulation requires a G/F-rich region [18]. Faust et al. revealed the NMR structure of the JD-GF domain of Hsp40 and found that upon binding to Hsp70, the G/F-rich region could form a stabilized α-helical structure [19]. Hsp70 are highly conserved in amino acid sequences and structures, and the main structure can be divided into four regions from the N-terminal to the C-terminal: a 40 kDa N-terminal ATPase domain (NBD), an 18 kDa substrate binding domain (SBD), a 15 kDa conserved peptide-binding functional domain, and a 10 kDa variable region at the C-terminal end, including an EEVD motif [20,21].
Hsp40 typically binds to heat shock protein 70 to facilitate ATP hydrolysis in order to perform its molecular chaperone function and participates in the next cycle upon completion. In the Hsp40–Hsp70 chaperone cycle, Hsp40 carries the substrate and stimulates the ATP hydrolysis to ADP on Hsp70, causing a conformational change in Hsp70 that allows it to bind tightly with the substrate. Subsequently, the nucleotide exchange factor (NEF) in turn stimulates the release of ADP from Hsp70 and facilitates the binding of a new round of ATP to the ATP-binding site on Hsp70, triggering the release of the substrate [12]. In addition, these processes are dependent on the ATPase activity of Hsp70, and Hsp40 is in turn thought to regulate how tightly Hsp70 binds to Hsp40-carrying substrates by modulating Hsp70 enzyme activity [22,23,24,25,26,27]. However, the mechanism by which Hsp40 regulates the ATPase activity of Hsp70 remains unclear. Faust et al. reported that the J domain as well as the GF domain of Hsp40 have important roles in the binding process with Hsp70 and that the binding of Hsp40 to Hsp70 requires the support of other components in addition to the roles of the J domain and the GF domain [19,28]. Therefore, we speculate that the binding site of Hsp40 to Hsp70 is not only present on the J domain. The specific binding sites of Hsp40 and Hsp70 as well as the sequence characteristics need to be further explored.
Since turbot is extremely demanding in terms of temperature and susceptible to a wide range of environmental stresses, the important functions of the HSPs’ family under different stresses make the regulation of their genes more pronounced in turbot, which can be investigated as a molecular biomarker of stress in fish. In addition, given the important functions of Hsp70 genes in response to environmental stresses and the abundant genome sequence resources, more and more Hsp70 gene families have been identified at the genome-wide level in a wide range of scleractinian fishes, and in recent years, the potential roles of these genes in responding to abiotic and biotic stresses have also been described [29,30,31,32]. Moreover, the high-quality genome [33] assemblies and multiple published RNA-seq dataset [34,35] resources about turbot in recent years have made it possible to systematically conduct identification, phylogenetic analyses, and functional research of Hsp40 and Hsp70 genes in turbot, which facilitates subsequent research.
In this study, genome-wide identification of the Hsp40 and Hsp70 families of heat shock proteins was performed, and the protein interaction between one pair of chaperones Dnajb1b and Hspa1a, which is involved in the immune response after bacterial infection, was analyzed.

2. Results

2.1. Identification of Hsp40 and Hsp70 Genes in the Turbot Genome

A total of 62 Hsp40 genes and 16 Hsp70 genes were found in the turbot genome based on the HMM search and domain analysis. As shown in Figure 1, the Hsp40 protein sequences vary from 118 AA to 4580 AA in length, and nearly 27.4% of the members were between 302 AA and 404 AA in length (Figure 1A). The protein molecular weights range from 12.78 kDa to 519.1 kDa and are mainly concentrated between 40.73 kD to 52.04 kDa (Figure 1B). The lengths of the Hsp70 protein sequence range from 439 AA to 986 AA, and the protein molecular weights extend from 47.64 kDa to 110.91 kDa, with a major concentration between 70.23 kDa and 79.47 kDa.
The distribution of conserved domains on the Hsp40 and Hsp70 proteins is illustrated in Figure 1C,D. The members of Hsp40 were classified into three subfamilies based on the Conserved Domains Database (CDD): DNAJA, DNAJB, and DNAJC. The DNAJA subfamily has J domain (JD), zinc finger domain, and carboxyl-terminal domain (CTD). The DNAJB subfamily has J domains, while not all members of DNAJB possess CTD. The DNAJC subfamily has a more diverse range of conserved domains; however, they all have J domains in common. In the Hsp70 protein family, several structurally similar and functionally related nucleotide binding domains (NBDs) are observed in each protein (Figure 1C,D). As shown in Figure 1C,D, the number of exons in Hsp40 and Hsp70 genes varied greatly from 1 to 57 and from 2 to 25. However, paralogous genes derived from the same group were shown to possess similar gene structures. For instance, both Dnajb5 and Dnajb5l contained three introns and four exons; both hspa1a and hspa1b in group HSPA1 contained only one intron and two exons, which was a composition obviously different from the members in the other groups (Table S1).

2.2. Chromosome Distribution Analysis

The 62 and 16 members of the Hsp40 and Hsp70 gene families are distributed on 22 chromosomes of the turbot reference genome (Figure 2). The results showed an average of three Hsp40 genes per chromosome. Among these chromosomes, the highest number of Hsp40 genes (five) was found on chromosomes Chr3 and Chr14, whereas the lowest number was found on chromosome Chr4, with only one Hsp40 gene. Hsp70 gene distribution also showed clustering on some chromosomes, such as six Hsp70 genes: hspa8a.1, hspa8a.2, hspa4a, hspa12b, hspa14, and hyou1, located on Chr2. In addition, two Hsp70 genes were observed on Chr9 and Chr11, respectively. The remaining six Hsp70 genes were distributed on different chromosomes, such as Chr1, Chr3, Chr10, Chr15, Chr17, and Chr20.

2.3. Motif Analysis of the Identified Hsp40 and Hsp70 Proteins

To identify the conserved motifs within the turbot Hsp40 and Hsp70 families, MEME analysis was conducted. Five conserved motifs (named Motif1 to Motif5) were identified in each of the A, B, and C subfamily proteins of the Hsp40 family, and ten motifs (named Motif1 to Motif10) were identified in the Hsp70 family.
The detailed information of Hsp40 and Hsp70 family motifs were listed. The results showed that the sequence lengths of the five motifs in subfamily A ranged from a minimum of 12 amino acids (Motif3) to a maximum of 50 amino acids (Motif4) (Figure 3A). The sequence lengths of the five motifs in subfamily B ranged from a minimum of 23 amino acids (Motif3) to a maximum of 50 amino acids (Motif4, 5) (Figure 3B). The sequence lengths of the five motifs in subfamily C ranged from a minimum of 15 amino acids (Motif4) to a maximum of 40 amino acids (Motif5) (Figure 3C). By comparing the three Hsp40 subfamilies, we found similar Motif1s containing the conserved tripeptide HPD in all three subfamilies, which is coincident with the conserved polypeptide sequence between the second and third alpha helices of the J domain. In addition, the motif corresponding to the GF-rich region was found in both subfamilies A and B (Motif3). No such motif was found in subfamily C. In Hsp70, only Motif4 was found in all the members, suggesting that Motif4 is a key sequence in the Hsp70 family (Figure 3D). This variable motif distribution may provide potent support for the phylogenetic classification and contribute to the functional diversity of the Hsp40 and Hsp70 protein families in turbot.

2.4. Phylogenetic Analysis

To understand the structural classification of the Hsp40 and Hsp70 gene families and clarify the phylogenetic relationships of Hsp40 and Hsp70 gene families among different vertebrate species, neighbor-joining (NJ) trees were constructed. Accession ID of the protein used in the tree is listed in Table S2. There are three subfamilies in the Hsp40 family, subfamily A, subfamily B, and subfamily C. Due to the large number of members in subfamily C, subfamily C was further divided into three groups (dnajc1–10, dnajc11–20, dnajc21–30) according to the order of sequence number for the phylogenetic tree construction. Hsp70 proteins of the selected vertebrates were constructed to build the phylogenetic tree. The results showed that subfamily A clustered into four clusters; subfamily B clustered into nine clusters; and subfamily C clustered into 31 clusters in the Hsp40 family (Figure 4A–E). The Hsp70 family clustered into nine clusters (Figure 4F). The results showed that turbot was more closely related to half-smooth tongue sole and rainbow trout.

2.5. Synteny and Gene Duplication Analysis

We counted the number of Hsp40 and Hsp70 genes in zebrafish, large yellow croaker, damselfish, rainbow trout, half-smooth tongue sole, and turbot and found that damselfish had the highest number of Hsp70 genes, which was likely due to the fact that they live in the tropics. The number of Hsp70 genes in turbot is similar to that of rainbow trout and large yellow croaker, but it differs in the number of half-smooth tongue sole, which is also a Pleuronectiformes. From the evolutionary perspective, it could be due to the fact that turbot have acquired more anti-stress genes. As shown in Figure 5B, 13 paralogues of Hsp40 genes were found in the turbot genome, indicating duplication events might be involved. Moreover, four duplication events involving the Hsp70 gene was found. Among these duplication events, two pairs of duplicated genes in Hsp40 occurred on adjacent chromosomes, for example, dnajc7l on chromosome 16 and dnajc7 on chromosome 17.
To further understand the possible evolutionary events involving the Hsp40 and Hsp70 gene families in different teleost species, three comparative syntenic maps of turbot with half-smooth tongue sole and rainbow trout were constructed. The results showed that 59 Hsp40 orthologues and 13 Hsp70 orthologues were found between turbot and half-smooth tongue sole as well as 127 Hsp40 orthologues and 27 Hsp70 orthologues between turbot and rainbow trout. The number of homologous gene pairs between turbot and rainbow trout was significantly more than half-smooth tongue sole, probably because both turbot and rainbow trout are more favored to survive in colder temperatures.

2.6. Protein-Protein Interaction Analysis

The 62 Hsp40 and 16 Hsp70 proteins were imported into the STRING database to obtain protein interaction network and screen for proteins with co-expression relationships. The analysis results were analyzed by Cytoscape 3.7.2, as shown in Figure 6 and Table S3 (Supplemental Information); nodes represent proteins, and edges represent relationships between proteins. The colors of nodes from yellow to red and the colors of edges from light blue to dark blue represent small to large degree values, respectively. The degree values of Hsp70 members were higher than those of Hsp40 by showing more edges and darker node colors in Hsp70 members.

2.7. Expression Analysis

To study Hsp40 and Hsp70 gene expression in turbot after infection with A. salmonicida, we analyzed their transcript abundance changes in transcriptome and RT–PCR data (Figure 7). After infection with A. salmonicida, five differentially expressed Hsp40 and Hsp70 genes were found in turbot spleen tissue, and fourteen differentially expressed Hsp40 and Hsp70 genes were discovered in turbot head kidney tissues. In the head kidney, only the expression of hsc70 was down-regulated, while the expression of the other 13 genes was up-regulated. In the spleen, dnajb1b, hspa1a, and hspa1b showed an overall trend of up-regulation of expression, and the overall expression of dnajc22 and hspa12a was down-regulated. In addition, dnajb1b and hspa1a show a more consistent expression trend compared to hspa1b, as their expression levels both increase with infection time. Therefore, we consider them co-expressed genes and hypothesize that they have interacting sites in the context of stress response between Hsp40 and Hsp70.

2.8. Molecular Docking Analysis

To investigate the protein interaction sites of the co-expressed genes dnajb1b and hspa1a, we conducted molecular docking analysis. Hydrogen bonds in the secondary structure of protein are important for maintaining structural stability. The shorter the bond distance, the greater the bond energy. We found the amino acid residue sites and hydrogen bond distances where Hsp40 and Hsp70 interact with each other by PDBePISA, and screened the sites that can form relatively stable interactions according to the hydrogen bond distances. As shown in Figure 8 and Table S4 (Supplemental Information), molecular chaperones Dnajb1b and Hspa1a formed a few groups of hydrogen bonds through amino acid residue sites, such as LYS 21-THR 297 (distance = 2.25 Å) and GLU 105-ALA 282 (distance = 2.51 Å), revealing that proteins Dnajb1b and Hspa1a formed a stable protein docking model.
When combining domain and motif analysis (Figure 9), certain amino acids in the JD and GF-rich regions of Dnajb1b can bind to specific amino acid sites within the same motif on the NBD of Hspa1a. D18, K21, S56, and D57 are binding sites in Hsp40 for Hsp70NBD. Among them, D18 and K21 are within a motif containing the conserved tripeptide HPD, and these two amino acids both bind to T297 of the NBD of Hspa1a. E105 and F107 are amino acids in the GF-rich region of Dnajb1b that bind to the NBD of Hspa1a. E105 can bind to three amino acids in the NBD, namely S280, A282, and S283, while F107 can only bind to S283. The CTD domain of Dnajb1b contains amino acids that bind to all three domains of Hspa1a. K191, L199, N200, R207, K211, Q216, and E232 are binding sites in the CTD of Dnajb1b for the NBD of Hspa1a. K197 is a binding site in the CTD of Dnajb1b for the SBD of Hspa1a, and S195 is a binding site in the CTD of Hsp40 for the CTD of Hspa1a.

3. Discussion

Hsp40 and Hsp70, as a pair of molecular chaperones, play a crucial role in regulating cellular protein homeostasis and participating in the immune activities of organisms [36,37]. However, in turbot, systematic analysis of Hsp40 and Hsp70 is limited, and their interaction sites remain unclear. In our study, we identified 62 Hsp40 genes and 16 Hsp70 genes. Based on the conserved domain database (CDD), we classified Hsp40 into three subfamilies: A, B, and C, and Hsp70 into nine groups. Additionally, we constructed six phylogenetic trees and identified five conserved motifs in the three subfamilies of Hsp40 and ten conserved motifs in Hsp70. In the syntenic analysis, we studied the synteny of genes within turbot and the synteny between turbot and other species, such as half-smooth tongue sole and rainbow trout. In the expression analysis, we constructed a protein–protein interaction (PPI) network for Hsp40 and Hsp70, identifying differentially expressed Hsp40 and Hsp70 genes in the spleen and head kidney of turbot following infection by Aeromonas salmonicida. We selected a pair of Hsp40 and Hsp70 proteins that exhibited interaction and co-expression for molecular docking analysis to further explore their role in the immune response. This study focuses on addressing the gap in the prediction of interaction sites between Hsp40 and Hsp70 in turbot, potentially aiding in the development of new immunoregulatory strategies or treatment. These findings provide essential foundational information for research on the structure and function of heat shock proteins and may guide related genetic engineering and drug development projects.
Based on the constructed Hsp40 and Hsp70 phylogenetic trees, it is evident that fish and terrestrial vertebrates are evolutionarily related to each other. However, due to their distinct living environments, the affinities between fish and terrestrial vertebrates differ significantly, with fish species being more closely related to each other. The phylogenetic tree reveals a close relationship between turbot and zebrafish, rainbow trout, and half-smooth tongue sole. This, combined with the synteny analysis, further confirms the close affinity between turbot and these species. Notably, the number of orthologous gene pairs between turbot and rainbow trout was significantly higher than that of half-smooth tongue sole and zebrafish for the Hsp40 and Hsp70 families. Therefore, when researching relevant vaccine adjuvants for turbot, we might consider referring to the research findings related to half-smooth tongue sole to provide new insights for the development of vac-cines against bacterial diseases in turbot.
In addition to homologues, the number of genes varies from species to species. Differences in the number of genes may be the result of gene duplication or gene loss during evolution. Gene duplication is the primary mechanism for gene family expansion and provides the raw material for the evolution of new gene functions. For example, in the constructed developmental-evolutionary tree, it can be found that dnajb1 has only one in humans, mice, chickens, three-toed turtles, and tropical clawed frogs but contains two genes (dnajb1a, dnajb1b) in zebrafish, rainbow trout, and half-smooth tongue sole, which are fishes. It was further confirmed that dnajb1a and dnajb1b are a set of duplicate genes in the synteny analysis. In addition, a significant up-regulation of dnajb1b expression was found in turbot after infection with A. salmonicida, and it was hypothesized that the duplication of dnajb1b was for better adaptation to environmental stresses in the aquatic environment.
In expression analysis, it can be found that dnajb1b, hspa1a, and dnaja1b show very similar changes in expression over time after infection with A. salmonicida and they all have a tendency toward up-regulated expression. In addition, upon combining the protein interaction analysis graphs (Figure 6) of Hsp40 and Hsp70, it can be found there is a strong interaction between dnajb1b and hspa1a. Therefore, based on the co-expression trend between them, dnajb1b and hspa1a are presumed to be one pair of chaperone molecules.
Based on previous studies on the interaction between Hsp40 and Hsp70, researchers usually focus only on the interaction between the JD domain of Hsp40 and Hsp70 [38,39,40,41,42,43,44,45,46,47,48]. However, through molecular docking analysis, we have also discovered that the CTD of Dnajb1b can bind to Hspa1a and even interact with three domains of Hspa1a. Surprisingly, most of the binding sites between these domains are located within the conserved motifs we have identified. These conserved motifs may provide potential information for identifying molecular targets of interaction. Additionally, we identified a highly conserved peptide motif in the CTD domain of Dnajb1b, where certain amino acids can bind to amino acids in all three domains of Hspa1a (Figure 9). A study has shown that the subfamily B of Hsp40 has a self-inhibition mechanism different from subfamily A [14]. Under normal conditions, the JD of Hsp40 is in a self-inhibited state and does not bind to Hsp70. When the CTD of Hsp40 binds to the CTD of Hsp70, the JD is released from self-inhibition and binds to Hsp70. In our analysis, there are 22 pairs of interacting amino acids between Dnajb1b and Hspa1a (Table S2), with half of them being pairs of amino acids in the CTD of Dnajb1b interacting with amino acids in Hspa1a. This highlights the significant role of the CTD of Hsp40. Moreover, according to Table S2, the hydrogen bond lengths between certain amino acids in the CTD of Dnajb1b and amino acids in the NBD of Hspa1a (e.g., LYS 211-ALA 56, distance = 2.38 Å; LYS 211-GLU 50, distance = 2.49 Å) and those between amino acids in the JD of Dnajb1b and amino acids in the NBD of Hspa1a (e.g., LYS 21-THR 297, Distance = 2.25 Å) are relatively short, forming stable binding models. This further demonstrates that the CTD of Hsp40 plays a crucial role in its interaction with Hsp70, and certain conserved amino acid motifs in the CTD may function similarly to the HPD conserved tripeptide in the JD.
Hsp40 and Hsp70 are widely present and highly conserved protein families in evolution from bacteria to mammals. Recently, Hsp40 and Hsp70 have been increasingly used in the treatment and prevention of fish diseases [49,50,51,52,53,54,55] and even in the human diseases, such as cancer [56,57,58,59], neurodegeneration [60,61,62], and autoimmune diseases [63,64,65,66,67]. Most of these diseases are caused by the inefficiency of HSPs to correct the misfolded proteins. The conserved motifs on Hsp40 interacting with Hsp70 identified in this paper may provide some molecular targets for improving the efficiency of Hsp40–Hsp70 interaction and the theory basis for disease treatment.

4. Materials and Methods

4.1. Identification and Characterization of Hsp40 and Hsp70 Genes in Turbot

The turbot reference genome sequence, annotation files, and protein sequences were collected from NCBI (https://www.ncbi.nlm.nih.gov, accessed on 1 March 2023). The turbot Hsp40 and Hsp70 protein sequences were used as a query database, and the turbot amino acid sequences were searched using the BLASTp program with default parameters. Moreover, the candidate protein sequences were submitted to the Conserved Domains Database (CDD) to further confirm the existence of conserved domains with an e-value of 10−4. The exon–intron structure of turbot Hsp40 and Hsp70 genes was visualized using the Gene Structure Display Server (GSDS, version 2.0, http://gsds.cbi.pku.edu.cn, accessed on 6 August 2022) with the genome GTF file.

4.2. Chromosomal Distribution Analysis of Hsp40 and Hsp70 Genes

The turbot genome file and genome Generic Feature Format (GFF) file were collected from NCBI. Information on chromosome length and gene position was obtained from the genome Generic Feature Format (GFF) file, and then the genes were visualized using the gene location visualize function in TBtools (v1.098761) [68].

4.3. Motif Analysis of Hsp40 and Hsp70

Motif analysis of Hsp40 member classifications was performed using the MEME program (http://meme-suite.org/tools/meme, accessed on 30 March 2023). Conserved motifs of Hsp40 and Hsp70, with the maximum number of motifs, was set to five and ten, respectively.

4.4. Phylogenetic Analysis of Hsp40 and Hsp70 Proteins

Hsp40 and Hsp70 proteins from the turbot and other selected vertebrate species were aligned using ClustalW [69] and MUSCLE [70]. Then, we constructed the neighbor-joining (NJ) phylogenetic trees with the pairwise deletion option based on the Jones–Taylor–Thornton (JTT) amino acid substitution model, using MEGA11.0 software [71] with bootstrap test of 1000 replicates. EvolView (v2.0) (https://www.evolgenius.info/evolview, accessed on 28 December 2023) was used to visualize the phylogenetic tree.

4.5. Synteny Analysis of Hsp40 and Hsp70 Genes

Gene duplication analysis of Hsp40 and Hsp70 genes was performed using TBtools software to identify duplication events among these Hsp40 and Hsp70 genes using the TBtools with the Graphics–Comparative Genomics–One Step MCScanX function [68].
Using the MCScanX program of Tbtools, we identified regions associated with Hsp40 and Hsp70 members with synteny among zebrafish, rainbow trout, and half-smooth tongue sole, and the related figures were generated by using TBtools software (Version 1.098761) [68,72]

4.6. Interaction Analysis of Hsp40 and Hsp70

Interaction analysis of Hsp40 and Hsp70 proteins was performed by String (http://cn.string-db.org, accessed on 22 April 2023). The protein interaction relationships were retrieved. The combined score from the export file was imported into Cytoscape to build a PPI network diagram [73].

4.7. Expression Analysis of Hsp40 and Hsp70

Two sets of transcriptomic data for Hsp40 and Hsp70 genes after infection with Aeromonas salmonicida were collected from the researches of Ting Xue et al. [34] and A. Romero et al. [37]. The 0 h (control) samples and the experiment groups of 12 h, 24 h, and 96 h samples35ere collected from the researches of Ting Xue et al. [34]. The 0 h (control) samples and the experiment group of 12 h samples were collected from the researches of A. Romero et al. [35].
The FPKM value was used to estimate the gene expression level, and the differentially expressed genes were obtained by DESeq with p-value < 0.05 and |log2FC| ≥ 1.

4.8. Molecular Docking of Dnajb1b-Hspa1a

Molecular Docking Analysis: Rigid protein–protein docking was performed between Dnajb1b and Hspa1a to investigate the relationships by using GRAMM-X (http://gramm.compbio.ku.edu/, accessed on 20 October 2023). The protein domains of Dnajb1b and Hspa1a were obtained from the Protein Data Bank PDB database (http://www.rcsb.org/, accessed on 20 October 2023), which was sourced from the Alpha fold [74]. PyMOL (Version 2.4), and PDBePISA (https://www.ebi.ac.uk/pdbe/pisa/, accessed on 20 October 2023) was used to investigate protein–protein interactions and further visual analysis [75]. Hydrogen bonds between residues with a bond distance of less than 5 Å are selected as stable interaction pairs.

5. Conclusions

We have identified 62 Hsp40 genes and 16 Hsp70 genes in the turbot genome. Their regulation pattern was related to tissue and time. Most of the Hsp40 and Hsp70 genes were up-regulated, and only a few of them were down-regulated when infected with A. salmonicida. According to the expression pattern, we built a relatively stable molecular model between Hsp40 and Hsp70, which revealed their interaction sites among each domain. The carboxyl-terminal domain of Hsp40 plays a crucial role in its interaction with Hsp70, and some conserved amino acid motifs in the carboxyl-terminal domain may be functionally similar to the conserved tripeptide of HPD in the J domain. This is important for the mechanistic understanding of the Hsp40–Hsp70 chaperone system and the development of small molecule inhibitors targeting Hsp70 or Hsp40, providing a theoretical basis for turbot disease resistance breeding and effective value for the prevention of certain diseases in turbot.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms25147963/s1.

Author Contributions

Conceptualization, L.S.; methodology, Y.G. (Yuanwei Geng) and Y.G. (Yuxuan Gai); software, Y.G. (Yuanwei Geng) and Y.G. (Yuxuan Gai); validation, A.J., F.Z. and K.D.; formal analysis, Y.G. (Yuanwei Geng) and S.Z.; data curation, Y.G. (Yuxuan Gai) and X.L.; writing—original draft preparation, Y.G. (Yuanwei Geng); writing—review and editing, L.S. and Y.G. (Yuxuan Gai); visualization, Y.Z. and S.Z.; supervision, L.S. and L.T.; funding acquisition, L.S. and L.T. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by Qingdao Agricultural University Doctoral Start-Up Fund, grant number: 663-1117004; Horizontal project: Development of Antarctic krill feed for turbot and research on nutrition and disease resistance, project number: H 20220193.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Pereiro, P.; Figueras, A.; Novoa, B. Turbot (Scophthalmus maximus) vs. VHSV (Viral Hemorrhagic Septicemia Virus): A Review. Front. Physiol. 2016, 7, 192. [Google Scholar] [CrossRef] [PubMed]
  2. Bjørndal, T.; Øiestad, V. The Development of a New Farmed Species: Production Technology and Markets for Turbot; Working Paper; Institute for Research in Economics and Business Administration: Bergen, Norway, 2010. [Google Scholar]
  3. Dallaire-Dufresne, S.; Tanaka, K.H.; Trudel, M.V.; Lafaille, A.; Charette, S.J. Virulence, genomic features, and plasticity of Aeromonas salmonicida subsp. salmonicida, the causative agent of fish furunculosis. Vet. Microbiol. 2014, 169, 1–7. [Google Scholar] [CrossRef] [PubMed]
  4. Hartl, F.U.; Bracher, A.; Hayer-Hartl, M. Molecular chaperones in protein folding and proteostasis. Nature 2011, 475, 324–332. [Google Scholar] [CrossRef]
  5. Li, J.; Qian, X.; Sha, B. Heat shock protein 40: Structural studies and their functional implications. Protein Pept. Lett. 2009, 16, 606–612. [Google Scholar] [CrossRef]
  6. Xu, H. Cochaperones enable Hsp70 to use ATP energy to stabilize native proteins out of the folding equilibrium. Sci. Rep. 2018, 8, 13213. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  7. Kumada, K.; Fuse, N.; Tamura, T.; Okamori, C.; Kurata, S. HSP70/DNAJA3 chaperone/cochaperone regulates NF-κB activity in immune responses. Biochem. Biophys. Res. Commun. 2019, 513, 947–951. [Google Scholar] [CrossRef] [PubMed]
  8. Ghazaei, C. Role and mechanism of the Hsp70 molecular chaperone machines in bacterial pathogens. J. Med. Microbiol. 2017, 66, 259–265. [Google Scholar] [CrossRef] [PubMed]
  9. Iyer, K.; Chand, K.; Mitra, A.; Trivedi, J.; Mitra, D. Diversity in heat shock protein families: Functional implications in virus infection with a comprehensive insight of their role in the HIV-1 life cycle. Cell Stress Chaperones 2021, 26, 743–768. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  10. Liu, J.Z.; Whitham, S.A. Overexpression of a soybean nuclear localized type-III DnaJ domain-containing HSP40 reveals its roles in cell death and disease resistance. Plant J. 2013, 74, 110–121. [Google Scholar] [CrossRef] [PubMed]
  11. Li, J.; Zhang, Y.; Liu, Y.; Zhang, Y.; Xiao, S.; Yu, Z. Co-expression of heat shock protein (HSP) 40 and HSP70 in Pinctada martensii response to thermal, low salinity and bacterial challenges. Fish Shellfish Immunol. 2016, 48, 239–243. [Google Scholar] [CrossRef] [PubMed]
  12. Liu, Q.; Liang, C.; Zhou, L. Structural and functional analysis of the Hsp70/Hsp40 chaperone system. Protein Sci. 2020, 29, 378–390. [Google Scholar] [CrossRef] [PubMed]
  13. Kampinga, H.H.; Craig, E.A. The HSP70 chaperone machinery: J proteins as drivers of functional specificity. Nat. Rev. Mol. Cell Biol. 2010, 11, 579–592. [Google Scholar] [CrossRef] [PubMed]
  14. Fan, C.Y.; Lee, S.; Cyr, D.M. Mechanisms for regulation of Hsp70 function by Hsp40. Cell Stress Chaperones 2003, 8, 309–316. [Google Scholar] [CrossRef] [PubMed]
  15. Alderson, T.R.; Kim, J.H.; Markley, J.L. Dynamical structures of Hsp70 and Hsp70-Hsp40 complexes. Structure 2016, 24, 1014–1030. [Google Scholar] [CrossRef] [PubMed]
  16. Karzai, A.W.; McMacken, R. A bipartite signaling mechanism involved in DnaJ-mediated activation of the Escherichia coli DnaK protein. J. Biol. Chem. 1996, 271, 11236–11246. [Google Scholar] [CrossRef] [PubMed]
  17. Kelley, W.L. The J-domain family and the recruitment of chaperone power. Trends Biochem. Sci. 1998, 23, 222–227. [Google Scholar] [CrossRef] [PubMed]
  18. Wu, S.; Hong, L.; Wang, Y.; Yu, J.; Yang, J.; Yang, J.; Zhang, H.; Perrett, S. Kinetics of the conformational cycle of Hsp70 reveals the importance of the dynamic and heterogeneous nature of Hsp70 for its function. Proc. Natl. Acad. Sci. USA 2020, 117, 7814–7823. [Google Scholar] [CrossRef] [PubMed]
  19. Faust, O.; Abayev-Avraham, M.; Wentink, A.S.; Maurer, M.; Nillegoda, N.B.; London, N.; Bukau, B.; Rosenzweig, R. HSP40 proteins use class-specific regulation to drive HSP70 functional diversity. Nature 2020, 587, 489–494. [Google Scholar] [CrossRef]
  20. Freeman, B.C.; Myers, M.P.; Schumacher, R.; Morimoto, R.I. Identification of a regulatory motif in Hsp70 that affects ATPase activity, substrate binding and interaction with HDJ-1. EMBO J. 1995, 14, 2281–2292. [Google Scholar] [CrossRef]
  21. Allan, R.K.; Ratajczak, T. Versatile TPR domains accommodate different modes of target protein recognition and function. Cell Stress Chaperones 2011, 16, 353–367. [Google Scholar] [CrossRef]
  22. Kuppner, M.C.; Gastpar, R.; Gelwer, S.; Nössner, E.; Ochmann, O.; Scharner, A.; Issels, R.D. The role of heat shock protein (hsp70) in dendritic cell maturation: Hsp70 induces the maturation of immature dendritic cells but reduces DC differentiation from monocyte precursors. Eur. J. Immunol. 2001, 31, 1602–1609. [Google Scholar] [CrossRef] [PubMed]
  23. Chen, T.; Cao, X. Stress for maintaining memory: HSP70 as a mobile messenger for innate and adaptive immunity. Eur. J. Immunol. 2010, 40, 1541–1544. [Google Scholar] [CrossRef] [PubMed]
  24. Rosenzweig, R.; Nillegoda, N.B.; Mayer, M.P.; Bukau, B. The Hsp70 chaperone network. Nat. Rev. Mol. Cell Biol. 2019, 20, 665–680. [Google Scholar] [CrossRef] [PubMed]
  25. Borges, T.J.; Wieten, L.; Van Herwijnen, M.J.; Broere, F.; Van Der Zee, R.; Bonorino, C.; Van Eden, W. The anti-inflammatory mechanisms of Hsp70. Front. Immunol. 2012, 3, 95. [Google Scholar] [CrossRef] [PubMed]
  26. Gao, X.-C.; Zhou, C.-J.; Zhou, Z.-R.; Wu, M.; Cao, C.-Y.; Hu, H.-Y. The C-terminal Helices of Heat Shock Protein 70 Are Essential for J-domain Binding and ATPase Activation*. J. Biol. Chem. 2012, 287, 6044–6052. [Google Scholar] [CrossRef]
  27. Ko, S.-K.; Kim, J.; Na, D.C.; Park, S.; Park, S.-H.; Hyun, J.Y.; Baek, K.-H.; Kim, N.D.; Kim, N.-K.; Park, Y.N.; et al. A Small Molecule Inhibitor of ATPase Activity of HSP70 Induces Apoptosis and Has Antitumor Activities. Chem. Biol. 2015, 22, 391–403. [Google Scholar] [CrossRef] [PubMed]
  28. Faust, O.; Rosenzweig, R. Structural and Biochemical Properties of Hsp40/Hsp70 Chaperone System. Adv. Exp. Med. Biol. 2020, 1243, 3–20. [Google Scholar] [CrossRef] [PubMed]
  29. Zhao, H.; Wang, Q.; Zhao, H.; Chen, C. Transcriptome profiles revealed high- and low-salinity water altered gill homeostasis in half-smooth tongue sole (Cynoglossus semilaevis). Comp. Biochem. Physiol. Part D Genom. Proteom. 2022, 42, 100989. [Google Scholar] [CrossRef]
  30. Huang, J.; Li, Y.; Liu, Z.; Kang, Y.; Wang, J. Transcriptomic responses to heat stress in rainbow trout Oncorhynchus mykiss head kidney. Fish Shellfish Immunol. 2018, 82, 32–40. [Google Scholar] [CrossRef] [PubMed]
  31. Song, L.; Li, C.; Xie, Y.; Liu, S.; Zhang, J.; Yao, J.; Jiang, C.; Li, Y.; Liu, Z. Genome-wide identification of Hsp70 genes in channel catfish and their regulated expression after bacterial infection. Fish Shellfish Immunol. 2016, 49, 154–162. [Google Scholar] [CrossRef]
  32. Du, Y.; Tang, X.; Sheng, X.; Xing, J.; Zhan, W. The influence of concentration of inactivated Edwardsiella tarda bacterin and immersion time on antigen uptake and expression of immune-related genes in Japanese flounder (Paralichthys olivaceus). Microb. Pathog. 2017, 103, 19–28. [Google Scholar] [CrossRef] [PubMed]
  33. Xu, X.-W.; Zheng, W.; Meng, Z.; Xu, W.; Liu, Y.; Chen, S. Identification of stress-related genes by co-expression network analysis based on the improved turbot genome. Sci. Data 2022, 9, 374. [Google Scholar] [CrossRef] [PubMed]
  34. Xue, T.; Liu, Y.; Cao, M.; Li, J.; Tian, M.; Zhang, L.; Wang, B.; Liu, X.; Li, C. Transcriptome analysis reveals deep insights into the early immune response of turbot (Scophthalmus maximus) induced by inactivated Aeromonas salmonicida vaccine. Fish Shellfish Immunol. 2021, 119, 163–172. [Google Scholar] [CrossRef] [PubMed]
  35. Romero, A.; Rey-Campos, M.; Pereiro, P.; Librán-Pérez, M.; Figueras, A.; Novoa, B. Transcriptomic analysis of turbot (Scophthalmus maximus) treated with zymosan a reveals that lncRNAs and inflammation-related genes mediate the protection conferred against Aeromonas salmonicida. Fish Shellfish Immunol. 2024, 147, 109456. [Google Scholar] [CrossRef] [PubMed]
  36. Roufayel, R.; Kadry, S. Molecular Chaperone HSP70 and Key Regulators of Apoptosis—A Review. Curr. Mol. Med. 2019, 19, 315–325. [Google Scholar] [CrossRef] [PubMed]
  37. Nunes, K.P.; de Oliveira, A.A. HSP70: From Signaling Mechanisms to Therapeutics. Biomolecules 2023, 13, 1141. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  38. Cyr, D.M.; Ramos, C.H. Specification of Hsp70 Function by Hsp40 Co-chaperones. Subcell. Biochem. 2023, 101, 127–139. [Google Scholar] [CrossRef] [PubMed]
  39. He, Y.; Wang, Z. The roles of HSP40/DNAJ protein family in neurodegenerative diseases. Zhejiang Da Xue Xue Bao Yi Xue Ban 2022, 51, 640–646. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  40. Chang, Y.L.; Yang, C.C.; Huang, Y.Y.; Chen, Y.A.; Yang, C.W.; Liao, C.Y.; Li, H.; Wu, C.S.; Lin, C.H.; Teng, S.C. The HSP40 family chaperone isoform DNAJB6b prevents neuronal cells from tau aggregation. BMC Biol. 2023, 21, 293. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  41. Cyr, D.M.; Ramos, C.H. Specification of Hsp70 function by Type I and Type II Hsp40. Subcell. Biochem. 2015, 78, 91–102. [Google Scholar] [CrossRef] [PubMed]
  42. Malinverni, D.; Jost Lopez, A.; De Los Rios, P.; Hummer, G.; Barducci, A. Modeling Hsp70/Hsp40 interaction by multi-scale molecular simulations and coevolutionary sequence analysis. eLife 2017, 6, e23471. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  43. Velasco-Carneros, L.; Cuéllar, J.; Dublang, L.; Santiago, C.; Maréchal, J.D.; Martín-Benito, J.; Maestro, M.; Fernández-Higuero, J.Á.; Orozco, N.; Moro, F.; et al. The self-association equilibrium of DNAJA2 regulates its interaction with unfolded substrate proteins and with Hsc70. Nat. Commun. 2023, 14, 5436, Erratum in Nat. Commun. 2024, 15, 99. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  44. Bascos, N.A.D.; Mayer, M.P.; Bukau, B.; Landry, S.J. The Hsp40 J-domain modulates Hsp70 conformation and ATPase activity with a semi-elliptical spring. Protein Sci. 2017, 26, 1838–1851. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  45. Ciesielski, S.J.; Young, C.; Ciesielska, E.J.; Ciesielski, G.L. The Hsp70 and JDP proteins: Structure-function perspective on molecular chaperone activity. Enzymes 2023, 54, 221–245. [Google Scholar] [CrossRef] [PubMed]
  46. Wickramaratne, A.C.; Liao, J.Y.; Doyle, S.M.; Hoskins, J.R.; Puller, G.; Scott, M.L.; Alao, J.P.; Obaseki, I.; Dinan, J.C.; Maity, T.K.; et al. J-domain Proteins form Binary Complexes with Hsp90 and Ternary Complexes with Hsp90 and Hsp70. J. Mol. Biol. 2023, 435, 168184. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  47. Malinverni, D.; Zamuner, S.; Rebeaud, M.E.; Barducci, A.; Nillegoda, N.B.; De Los Rios, P. Data-driven large-scale genomic analysis reveals an intricate phylogenetic and functional landscape in J-domain proteins. Proc. Natl. Acad. Sci. USA 2023, 8, e2218217120. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  48. Horne, B.E.; Li, T.; Genevaux, P.; Georgopoulos, C.; Landry, S.J. The Hsp40 J-domain stimulates Hsp70 when tethered by the client to the ATPase domain. J. Biol. Chem. 2010, 285, 21679–21688. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  49. Hou, G.; Zhang, Q.; Li, C.; Ding, G.; Hu, L.; Chen, X.; Lv, Z.; Fan, Y.; Zou, J.; Xiao, T.; et al. An Aquareovirus Exploits Membrane-Anchored HSP70 To Promote Viral Entry. Microbiol. Spectr. 2023, 11, e0405522. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  50. Pham, T.H.; Cheng, T.C.; Wang, P.C.; Chen, S.C. Protective efficacy of four heat-shock proteins as recombinant vaccines against photobacteriosis in Asian seabass (Lates calcarifer). Fish Shellfish Immunol. 2021, 111, 179–188. [Google Scholar] [CrossRef] [PubMed]
  51. Li, X.; Fan, K.; Liu, Y.; Liu, Y.; Liu, P.F. Administration of a recombinant ALDH7A1 (rA7) indicates potential regulation of the metabolite and immunology pathways in Atlantic salmon infected with Aeromonas salmonicida. J. Fish Dis. 2021, 44, 961–977. [Google Scholar] [CrossRef] [PubMed]
  52. Pham, P.H.; Sokeechand, B.S.H.; Hamilton, M.E.; Misk, E.; Jones, G.; Lee, L.E.J.; Lumsden, J.S.; Bols, N.C. VER-155008 induced Hsp70 proteins expression in fish cell cultures while impeding replication of two RNA viruses. Antivir. Res. 2019, 162, 151–162. [Google Scholar] [CrossRef] [PubMed]
  53. Rao, Y.; Wan, Q.; Su, H.; Xiao, X.; Liao, Z.; Ji, J.; Yang, C.; Lin, L.; Su, J. ROS-induced HSP70 promotes cytoplasmic translocation of high-mobility group box 1b and stimulates antiviral autophagy in grass carp kidney cells. J. Biol. Chem. 2018, 293, 17387–17401. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  54. Shan, L.P.; Chen, X.H.; Ling, F.; Zhu, B.; Wang, G.X. Targeting Heat Shock Protein 70 as an antiviral strategy against grass carp reovirus infection. Virus Res. 2018, 247, 1–9. [Google Scholar] [CrossRef] [PubMed]
  55. Vahdatiraad, L.; Heidari, B.; Zarei, S.; Sohrabi, T.; Ghafouri, H. Protective effects of HSP Inducer on diazinon-exposed stellate sturgeon (Acipenser stellatus) fry: Insights on HSP70 gene expression, immune response, and enzyme indices. PLoS ONE 2023, 18, e0294188. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  56. Nitzsche, B.; Höpfner, M.; Biersack, B. Synthetic Small Molecule Modulators of Hsp70 and Hsp40 Chaperones as Promising Anticancer Agents. Int. J. Mol. Sci. 2023, 24, 4083. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  57. Theivanthiran, B.; Yarla, N.; Haykal, T.; Nguyen, Y.V.; Cao, L.; Ferreira, M.; Holtzhausen, A.; Al-Rohil, R.; Salama, A.K.S.; Beasley, G.M.; et al. Tumor-intrinsic NLRP3-HSP70-TLR4 axis drives premetastatic niche development and hyperprogression during anti-PD-1 immunotherapy. Sci. Transl. Med. 2022, 14, eabq7019. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  58. Zhao, K.; Zhou, G.; Liu, Y.; Zhang, J.; Chen, Y.; Liu, L.; Zhang, G. HSP70 Family in Cancer: Signaling Mechanisms and Therapeutic Advances. Biomolecules 2023, 13, 601. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  59. Elmallah, M.I.Y.; Cordonnier, M.; Vautrot, V.; Chanteloup, G.; Garrido, C.; Gobbo, J. Membrane-anchored heat-shock protein 70 (Hsp70) in cancer. Cancer Lett. 2020, 469, 134–141. [Google Scholar] [CrossRef] [PubMed]
  60. Perez, V.A.; Sanders, D.W.; Mendoza-Oliva, A.; Stopschinski, B.E.; Mullapudi, V.; White, C.L.; Joachimiak, L.A.; Diamond, M.I. DnaJC7 specifically regulates tau seeding. eLife 2023, 12, e86936. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  61. Dong, Y.; Li, T.; Ma, Z.; Zhou, C.; Wang, X.; Li, J. HSPA1A, HSPA2, and HSPA8 Are Potential Molecular Biomarkers for Prognosis among HSP70 Family in Alzheimer’s Disease. Dis. Markers 2022, 2022, 9480398. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  62. Chanteloup, G.; Cordonnier, M.; Moreno-Ramos, T.; Pytel, V.; Matías-Guiu, J.; Gobbo, J.; Cabrera-Martín, M.N.; Gómez-Pinedo, U.; Garrido, C.; Matías-Guiu, J.A. Exosomal HSP70 for Monitoring of Frontotemporal Dementia and Alzheimer’s Disease: Clinical and FDG-PET Correlation. J. Alzheimer’s Dis. 2019, 71, 1263–1269. [Google Scholar] [CrossRef] [PubMed]
  63. Niu, L.; Lou, F.; Sun, Y.; Sun, L.; Cai, X.; Liu, Z.; Zhou, H.; Wang, H.; Wang, Z.; Bai, J.; et al. A micropeptide encoded by lncRNA MIR155HG suppresses autoimmune inflammation via modulating antigen presentation. Sci. Adv. 2020, 6, eaaz2059. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  64. van Eden, W. Immune tolerance therapies for autoimmune diseases based on heat shock protein T-cell epitopes. Philos. Trans. R. Soc. B Biol. Sci. 2018, 373, 20160531. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  65. Rodriguez-Iturbe, B.; Lanaspa, M.A.; Johnson, R.J. The role of autoimmune reactivity induced by heat shock protein 70 in the pathogenesis of essential hypertension. Br. J. Pharmacol. 2019, 176, 1829–1838. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  66. Tukaj, S.; Mantej, J.; Sitko, K.; Zillikens, D.; Ludwig, R.J.; Bieber, K.; Kasperkiewicz, M. Pathological Relevance of Anti-Hsp70 IgG Autoantibodies in Epidermolysis Bullosa Acquisita. Front. Immunol. 2022, 13, 877958. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  67. Wang, L.; Fu, Y.; Yu, B.; Jiang, X.; Liu, H.; Liu, J.; Zha, B.; Chu, Y. HSP70, a Novel Regulatory Molecule in B Cell-Mediated Suppression of Autoimmune Diseases. J. Mol. Biol. 2021, 433, 166634. [Google Scholar] [CrossRef] [PubMed]
  68. Chen, C.J.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.H.; Xia, R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef] [PubMed]
  69. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; McGettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23, 2947–2948. [Google Scholar] [CrossRef] [PubMed]
  70. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  71. Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular evolutionary genetics analysis version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef]
  72. Wang, Y.; Tang, H.; DeBarry, J.D.; Tan, X.; Li, J.; Wang, X.; Lee, T.H.; Jin, H.; Marler, B.; Guo, H.; et al. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012, 40, e49. [Google Scholar] [CrossRef] [PubMed]
  73. Doncheva, N.T.; Morris, J.H.; Gorodkin, J.; Jensen, L.J. Cytoscape StringApp: Network Analysis and Visualization of Proteomics Data. J. Proteome Res. 2019, 18, 623–632. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  74. Abramson, J.; Adler, J.; Dunger, J.; Evans, R.; Green, T.; Pritzel, A.; Ronneberger, O.; Willmore, L.; Ballard, A.J.; Bambrick, J.; et al. Accurate structure prediction of biomolecular interactions with AlphaFold 3. Nature 2024, 630, 493–500. [Google Scholar] [CrossRef] [PubMed]
  75. DeLano, W.L. The PyMOL Molecular Graphics System. 2002. Available online: http://www.pymol.org (accessed on 20 October 2023).
Ijms 25 07963 g001aIjms 25 07963 g001b
Figure 1. The characteristics and domains of Hsp40 and Hsp70 members in S. maximus. (A) Protein sequence length distribution; AA, amino acid; (B) the protein molecular weight distribution; (C) distribution of conserved domains and gene structure on the Hsp40; (D) distribution of conserved domains and gene structure on the Hsp70.
Figure 1. The characteristics and domains of Hsp40 and Hsp70 members in S. maximus. (A) Protein sequence length distribution; AA, amino acid; (B) the protein molecular weight distribution; (C) distribution of conserved domains and gene structure on the Hsp40; (D) distribution of conserved domains and gene structure on the Hsp70.
Ijms 25 07963 g001aIjms 25 07963 g001b
Ijms 25 07963 g002
Figure 2. The genome positions of Hsp40 and Hsp70 members in S. maximus. The Hsp40 genes and the Hsp70 genes are labeled in red and blue respectively.
Figure 2. The genome positions of Hsp40 and Hsp70 members in S. maximus. The Hsp40 genes and the Hsp70 genes are labeled in red and blue respectively.
Ijms 25 07963 g002
Ijms 25 07963 g003
Figure 3. Conserved motif information for Hsp40 and Hsp70 proteins in S. maximus, and each colored rectangular box represents a motif. (A) Motif information in turbot Hsp40 family subfamily A; (B) motif information in turbot Hsp40 family subfamily B; (C) motif information in turbot Hsp40 family subfamily C; (D) motif information in turbot Hsp70 family.
Figure 3. Conserved motif information for Hsp40 and Hsp70 proteins in S. maximus, and each colored rectangular box represents a motif. (A) Motif information in turbot Hsp40 family subfamily A; (B) motif information in turbot Hsp40 family subfamily B; (C) motif information in turbot Hsp40 family subfamily C; (D) motif information in turbot Hsp70 family.
Ijms 25 07963 g003
Ijms 25 07963 g004
Figure 4. The phylogenetic tree constructed from Hsp40 and Hsp70 proteins in S. maximus, C. semilaevis, O. mykiss, D. rerio, X. tropicalis, T. c. triunguis, G. gallus, M. musculus, and H. sapiens. Different colors indicate different groups. The genes for turbot, half-smooth tongue sole, and rainbow trout are labeled with red dots and black and gray triangles, respectively. (A) Phylogenetic tree of subfamily A of the Hsp40 family; (B) phylogenetic tree of subfamily B of the Hsp40 family; (C) phylogenetic tree of dnajc1~dnajc10 in subfamily C of the Hsp40 family; (D) phylogenetic tree of dnajc11~dnajc20 in subfamily C of the Hsp40 family; (E) phylogenetic tree of dnajc21~dnajc30 in subfamily C of the Hsp40 family; (F) phylogenetic tree of Hsp70 family.
Figure 4. The phylogenetic tree constructed from Hsp40 and Hsp70 proteins in S. maximus, C. semilaevis, O. mykiss, D. rerio, X. tropicalis, T. c. triunguis, G. gallus, M. musculus, and H. sapiens. Different colors indicate different groups. The genes for turbot, half-smooth tongue sole, and rainbow trout are labeled with red dots and black and gray triangles, respectively. (A) Phylogenetic tree of subfamily A of the Hsp40 family; (B) phylogenetic tree of subfamily B of the Hsp40 family; (C) phylogenetic tree of dnajc1~dnajc10 in subfamily C of the Hsp40 family; (D) phylogenetic tree of dnajc11~dnajc20 in subfamily C of the Hsp40 family; (E) phylogenetic tree of dnajc21~dnajc30 in subfamily C of the Hsp40 family; (F) phylogenetic tree of Hsp70 family.
Ijms 25 07963 g004
Ijms 25 07963 g005
Figure 5. Duplication event analysis for the Hsp40 and Hsp70 gene families in the turbot genome and the synteny analysis between turbot and other fishes. (A) Number of Hsp40 and Hsp70 genes in six fish species; (B) the paralogues of Hsp40 and Hsp70 in the turbot genome—the blue and red lines indicate the Hsp40 and Hsp70 families members, respectively; (C) sytenic analysis of Hsp40 genes between S. maximus, C. semilaevis and O. mykiss—the blue lines indicate the Hsp40 family members in different species; (D) sytenic analysis of Hsp70 genes between S. maximus, C. semilaevis, and O. mykiss—the violet lines indicate the Hsp70 family members in different species.
Figure 5. Duplication event analysis for the Hsp40 and Hsp70 gene families in the turbot genome and the synteny analysis between turbot and other fishes. (A) Number of Hsp40 and Hsp70 genes in six fish species; (B) the paralogues of Hsp40 and Hsp70 in the turbot genome—the blue and red lines indicate the Hsp40 and Hsp70 families members, respectively; (C) sytenic analysis of Hsp40 genes between S. maximus, C. semilaevis and O. mykiss—the blue lines indicate the Hsp40 family members in different species; (D) sytenic analysis of Hsp70 genes between S. maximus, C. semilaevis, and O. mykiss—the violet lines indicate the Hsp70 family members in different species.
Ijms 25 07963 g005
Ijms 25 07963 g006
Figure 6. PPI network of the Hsp40 and Hsp70 proteins.
Figure 6. PPI network of the Hsp40 and Hsp70 proteins.
Ijms 25 07963 g006
Ijms 25 07963 g007
Figure 7. Expression analysis of differentially expressed genes of Hsp40 and Hsp70 families after infection with A. salmonicida in turbot.
Figure 7. Expression analysis of differentially expressed genes of Hsp40 and Hsp70 families after infection with A. salmonicida in turbot.
Ijms 25 07963 g007
Ijms 25 07963 g008
Figure 8. Dnajb1b and Hspa1a protein molecular docking model. Surface diagram of the docking model and their interfacing residues between Dnajb1b and Hspa1a protein (Dnajb1b, yellow; Hspa1a, blue; hydrogen bond interaction, dotted line).
Figure 8. Dnajb1b and Hspa1a protein molecular docking model. Surface diagram of the docking model and their interfacing residues between Dnajb1b and Hspa1a protein (Dnajb1b, yellow; Hspa1a, blue; hydrogen bond interaction, dotted line).
Ijms 25 07963 g008
Ijms 25 07963 g009
Figure 9. Dnajb1b and Hspa1a protein molecular docking model combined with motif analysis (Dnajb1b, yellow; Hspa1a, blue; hydrogen bonding interactions, dashed lines). (A) Molecular structure of the JD of Dnajb1b interacting with the NBD of Hspa1a; (B) molecular structure of the GF-rich regions of Dnajb1b interacting with the NBD of Hspa1a; (C) molecular structure of the CTD of Dnajb1b interacting with the NBD of Hspa1a; (D) molecular structure of the CTD of Dnajb1b interacting with both the SBD and the CTD of Hspa1a; (E) domain interactions between Dnajb1b and Hspa1a and the motif where the interactions amino acid sites are located.
Figure 9. Dnajb1b and Hspa1a protein molecular docking model combined with motif analysis (Dnajb1b, yellow; Hspa1a, blue; hydrogen bonding interactions, dashed lines). (A) Molecular structure of the JD of Dnajb1b interacting with the NBD of Hspa1a; (B) molecular structure of the GF-rich regions of Dnajb1b interacting with the NBD of Hspa1a; (C) molecular structure of the CTD of Dnajb1b interacting with the NBD of Hspa1a; (D) molecular structure of the CTD of Dnajb1b interacting with both the SBD and the CTD of Hspa1a; (E) domain interactions between Dnajb1b and Hspa1a and the motif where the interactions amino acid sites are located.
Ijms 25 07963 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Geng, Y.; Gai, Y.; Zhang, Y.; Zhao, S.; Jiang, A.; Li, X.; Deng, K.; Zhang, F.; Tan, L.; Song, L. Genome-Wide Identification and Interaction Analysis of Turbot Heat Shock Protein 40 and 70 Families Suggest the Mechanism of Chaperone Proteins Involved in Immune Response after Bacterial Infection. Int. J. Mol. Sci. 2024, 25, 7963. https://doi.org/10.3390/ijms25147963

AMA Style

Geng Y, Gai Y, Zhang Y, Zhao S, Jiang A, Li X, Deng K, Zhang F, Tan L, Song L. Genome-Wide Identification and Interaction Analysis of Turbot Heat Shock Protein 40 and 70 Families Suggest the Mechanism of Chaperone Proteins Involved in Immune Response after Bacterial Infection. International Journal of Molecular Sciences. 2024; 25(14):7963. https://doi.org/10.3390/ijms25147963

Chicago/Turabian Style

Geng, Yuanwei, Yuxuan Gai, Yanping Zhang, Shengwei Zhao, Anlan Jiang, Xueqing Li, Kaiqing Deng, Fuxuan Zhang, Lingling Tan, and Lin Song. 2024. "Genome-Wide Identification and Interaction Analysis of Turbot Heat Shock Protein 40 and 70 Families Suggest the Mechanism of Chaperone Proteins Involved in Immune Response after Bacterial Infection" International Journal of Molecular Sciences 25, no. 14: 7963. https://doi.org/10.3390/ijms25147963

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop