Next Article in Journal
Investigation of the Molecular Evolution of Treg Suppression Mechanisms Indicates a Convergent Origin
Next Article in Special Issue
Disarm The Bacteria: What Temperate Phages Can Do
Previous Article in Journal
“Losing the Brakes”—Suppressed Inhibitors Triggering Uncontrolled Wnt/ß-Catenin Signaling May Provide a Potential Therapeutic Target in Elderly Acute Myeloid Leukemia
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
CIMB
Volume 45
Issue 1
10.3390/cimb45010041
cimb-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 Analysis of the Heat-Shock Protein Gene in L. edodes and Expression Pattern Analysis under Heat Shock

by
Xu Zhao
1,2,3,
Kaiyong Yin
1,
Rencai Feng
1,2,
Renyun Miao
1,2,
Junbin Lin
1,2,
Luping Cao
4,
Yanqing Ni
5,
Wensheng Li
5 and
Qin Zhang
1,2,*
1
Institute of Urban Agriculture, Chinese Academy of Agricultural Sciences, Chengdu 610299, China
2
Chengdu National Agricultural Science and Technology Center, Chengdu 610299, China
3
Facility Agriculture and Equipment Research Institute, Gansu Academy of Agri-engineering Technology, Wuwei 733006, China
4
College of Life Science and Technology, Gansu Agricultural University, Lanzhou 730070, China
5
College of Food and Biological Engineering, Chengdu University, Chengdu 610106, China
*
Author to whom correspondence should be addressed.
Curr. Issues Mol. Biol. 2023, 45(1), 614-627; https://doi.org/10.3390/cimb45010041
Submission received: 30 November 2022 / Revised: 3 January 2023 / Accepted: 5 January 2023 / Published: 9 January 2023
(This article belongs to the Special Issue Microbial Engineering: Gene Expression Regulation and Its Application)
Browse Figures
Review Reports Versions Notes

Abstract

:
Lentinula edodes (L. edodes), one of the most popular edible mushrooms in China, is adversely affected by high temperature. Heat shock proteins (HSPs) play a crucial role in regulating the defense responses against the abiotic stresses in L. edodes. Some HSPs in L. edodes have been described previously, but a genome-wide analysis of these proteins is still lacking. Here, the HSP genes across the entire genome of the L. edodes mushroom were identified. The 34 LeHSP genes were subsequently classified into six subfamilies according to their molecular weights and the phylogenetic analysis. Sequence analysis showed that LeHSP proteins from the same subfamily have conserved domains and one to five similar motifs. Except for Chr 5 and 9, 34 LeHSPs genes were distributed on the other eight chromosomes. Three pairs of paralogs were identified because of sequence alignment and were confirmed as arising from segmental duplication. In LeHSPs’ promoters, different numbers of heat shock elements (HSEs) were predicted. The expression profiles of LeHSPs in 18N44 and 18 suggested that the thermo-tolerance of strain 18N44 might be related to high levels of LeHSPs transcript in response to heat stress. The quantitative real-time PCR (qRT-PCR) analysis of the 16 LeHSP genes in strains Le015 and Le027 verified their stress-inducible expression patterns under heat stress. Therefore, these comprehensive findings provide useful in-depth information on the evolution and function of LeHSPs and lay a theoretical foundation in breeding thermotolerant L. edodes varieties.

1. Introduction

The basidiomycete Lentinula edodes (Berk.) Singer (Xianggu in Chinese or Shiitake in Japanese) is an edible mushroom that originates from China. It was first cultivated in China approximately 800 years ago and is currently the second most commonly produced edible mushroom in the world, with applications as decoctions, essences, and alternative medicines [1,2]. In 2020, the L. edodes production reached 11.64 million tons in China, accounting for 90% of the total global production (https://www.chyxx.com/industry/202201/992805.html (accessed on 10 Octorber 2022)).
During the cultivation of L. edodes, the environmental temperature is extremely rigorous and the optimum temperature for mycelia growth is 24~27 °C. Growth temperature higher than 32 °C will result in poor mycelia, and growth temperature higher than 35 °C can lead to severe retardation or even death. Hence, heat stress is considered one of the most important environmental stressors challenging L. edodes production.
To reduce the adverse effects under heat stress, fungus has evolved a series of self-defense mechanisms against heat shock to protect its body from damage [3,4]. Heat shock proteins (HSPs), also referred to as heat stress proteins, are a kind of highly conserved proteins across primitive prokaryotes (bacteria) to higher level eukaryotes [5]. The first HSP was discovered in Drosophila melanogaster under heat stress [6]. Later, many members of HSP families were found in almost all living organisms. HSPs are produced in response to thermal stress as well as various stressful conditions including drought, salinity, extreme temperatures, chemical toxicity, and oxidative stress [7,8,9], thereby making HSPs extremely important “moonlight” protein families [10]. The expression of HSP genes is regulated by heat shock factors (HSFs). These HSFs can recognize and bind to specific DNA sequences such as 5′-nGAAC-3′ [11].
Members of HSP families share a highly conserved domain named “heat-shock” and are classified into six major categories according to their molecular weights and sequence homologies: HSP100, HSP90, HSP70, HSP60, HSP40, and small HSP (SHSP) [12].The cellular mechanisms underlying the function of HSPs under abiotic stress are not fully understood, but many of them were found to act as chaperones or co-chaperones involved in protein synthesis and folding [7], in the degradation of abnormal protein [13], and in cellular localization [14] during stress conditions. HSPs in a few model plants have been well characterized. Arabidopsis 21, 27, 18, 7, and 8 genes have been identified for encoding, including the HSF, SHSP, HSP70, HSP90, and HSP100 families, respectively [15]. Genome-wide identifications of HSPs have been performed in wheat, soybean, and rice [12,16,17,18]. To date, limited information is available on genome-wide identification of the entire HSP families in L. edodes. The response of HSPs to heat stress is unclear and clarifying the mechanism of heat stress tolerance in L. edodes is vital.
HSP families play a crucial role in response to abiotic stress, but the relationship between heat stress resistance and HSP expression is ambiguous in L. edodes. Therefore, a total of 34 HSP genes are identified in the L. edodes genome with bioinformatics methods. Then, the physicochemical characteristics, chromosomal location, phylogenetic tree, gene structure, conserved motifs, homologous gene pairs, expression in response to heat stress, and heat shock elements (HSEs) in the promoter of LeHSP genes are analyzed. The LeHSP expression profiles under heat stress indicate that different combinations of LeHSPs are involved in the thermotolerance of strains 18N44 and Le027. This study systematically elucidates the evolution of HSPs in L. edodes and relevant role of LeHSP genes in response to heat stress. Thus, this work not only lays a foundation for further research on the biological function of LeHSP genes, but also provides a reference base for the molecular mechanism of heat resistance breeding in L. edodes.

2. Materials and methods

2.1. Identification of HSPs in L. edodes

The protein sequences of L. edodes were downloaded from the NCBI database (https://www.ncbi.nlm.nih.gov/genome/?term=lentinula+edodes (accessed on 1 July 2022)) and used for the construction of the local database [19]. To identify the HSP family, the hidden Markov models of the HSP20 (PF00011), HSP40 (PF00226 and PF01556), HSP60 (PF00118), HSP70 (PF00012), HSP90 (PF00183), and HSP100 domains (PF02861 and PF10431) were downloaded from the Pfam database (https://pfam.xfam.org/ (accessed on 4 July 2022)) [20] and used as queries (p < 0.0001) to search against the local database of L. edodes by using the Simple HMMER search plug-in implemented in the TBtools software (v1.098773) [21]. To ascertain all HSP family members, a blastp search was performed by employing all HSP amino acid sequences from the HMMER search as queries against the non-redundant protein sequences of L. edodes with an e-value of 1e−4. After the removal of the redundant HSP proteins, the putative HSP proteins were verified for the presence of HSP related domains using other databases including the NCBICDD (https://www.ncbi.nlm.nih.gov/cdd (accessed on 18 September 2022)), HAMMER scan (https://www.ebi.ac.uk/Tools/hmmer/search/hmmscan (accessed on 18 September 2022)), and SMART (http://smart.embl-heidelberg.de/ (accessed on 18 September 2022)). The identified genes were named with a prefix ‘Le’ for HSP40, HSP60, HSP70, HSP90, HSP100 and ‘LeS’ for small HSP, and numbered in increasing order with their position on the chromosome proceeding from the short to the long arm.

2.2. Chromosomal Localization and Gene Structure

The chromosomal localizations and intron–exon structures of HSP genes were retrieved from the annotation files of L. edodes and visualized by Map Genes on Genome and Gene Structure View (Advanced), respectively. Both plug-ins were implemented in the TBtools software. The motifs of HSP proteins were identified by the MEME program (Version 5.4.1, https://meme-suite.org/meme/tools/meme (accessed on 4 July 2022)) [22] with the parameter of a maximum of 5 motifs and other parameters as default. The number of amino acids, molecular weights, and theoretical isoelectric points (pI) were analyzed by the Geneious primer software (https://www.geneious.com/prime/ (accessed on 6 May 2022)) [23].

2.3. Duplication of LeHSPs in L. edodes

Genes were defined as paralogues if they meet the following criteria: (a) the length of the alignable sequence covers more than 75% of the longer gene; and (b) the aligned regions have more than 75% similarity. In addition, paralogous pairs were defined as tandemly duplicated genes if they are separated by five or fewer genes within a 100-kb region; if the paralogous pairs were mapped onto duplicated chromosomal fragment, then they were designated as segmental duplications [24].

2.4. Phylogenetic Reconstruction

The full amino acid sequences of the HSPs of members from Agaricusbisporus, Pleurotuspulmonarius, Pleurotuseryngii, and Volvariellavolvacea were downloaded from the NCBI.
Phylogenetic analysis was performed using the Geneious primer software (Version 2022.2). In brief, all LeHSP protein sequences were first aligned according to MUSCLE Alignment with default parameters in Geneious primer. A Geneious tree was subsequently constructed based on the aligned sequences using the UPGMA method with the Jukes–Cantor model. The final phylogenetic tree was embellished and visualized by iTOL (https://itol.embl.de/ (accessed on 1 July 2022)) [25] and Adobe Illustrator 2020.

2.5. Heat Shock Element Analysis in the Promoter Regions of HSP Genes

The promoter sequences (1 kb upstream of the initiation ATG start codon) of each HSP gene were extracted using the Gtf/Gff3 Sequences Extrator plug-in TBtools. The locations and numbers of the heat stock element (HSE, 5′nGAAC-3′) were predicted and visualized by TBtools.

2.6. Strain Materials and High-Temperature Treatments

L. edodes strains 18 and 18N44 were used for RNA-seq analysis. Both strains were obtained from the Institute of Edible Fungi, Shanghai Academy of Agricultural Sciences. Strain 18 is an elite L. edodes in agronomic traits but is sensitive to heat stress, and strain 18N44 was selected from the protoplasts of the original strain 18 which was exposed to ultra-violet rays. Strain 18N44 showed good performance in both agronomic traits and thermotolerance [26]. A similar method described by Zhao et al. [27] with some modification was used for high-temperature treatment. In brief, the mycelia of both strains 18 and 18N44 were incubated at 25 °C for 14 days on PDA medium and then cultured in the PDB medium. After incubating at 25 °C for 14 days with shaking at 150 r/min, the mycelia were treated with high-temperature stress at 37 °C for 0, 4, 12, and 24 h. The mycelia were collected and immediately frozen in liquid nitrogen and stored at −80 °C for RNA isolation and RNA-seq. Three biological replications were carried out for each sample. Two other L. edodes strains Le015 and Le027 were used for expression analysis. Strain Le015 is sensitive to heat-stress, and Le027 is a thermo-resistant strain. The same cultivation conditions and heat stress treatments were applied for strains Le015 and Le027. The mycelia were collected for RNA isolation, and the obtained RNAs were used for qPCR analysis.

2.7. Extraction and Purification of Total RNA from Mycelia

Total RNA was isolated using the Fungal Total RNA Isolation Kit (B518629, Sangon Biotech, China) and purified following the manufacturer’s protocol. The purity and concentration of RNA was measured by NanoDrop 2000 C (Thermo Fisher Scientific, Waltham, MA, USA).

2.8. RNA-Sequencing and Analysis of HSP Genes

The qualified RNA samples from strains 18 and 18N44 were sent to Oebiotech (Shanghai, China) for transcriptome sequencing. The FPKM [28] and read counts value of each unigene was calculated using bowtie2 [29] and eXpress [30]. The heatmap of the gene expression pattern was generated according to the FPKM values using the R (version 4.0.4) package heatmap (1.0.8).

2.9. Gene Expression by RT-qPCR

Quantitative reverse transcription-PCR was carried out using II Green One-Step qRT-PCR SuperMix (TransScrip, China) according to the manufacturer’s instructions. The reference genes of TUB and UBI [27,31,32] were used to normalize the differences between samples. The relative expressions of HSP genes were calculated using the EΔΔCt method [33] against the geometric mean of two internal reference genes. The primers were designed by Primer 5.0. The primers for RT-qPCR are listed in Table S1. The results were calculated by Excel and displayed as means ± standard deviation (SD).

3. Results

3.1. Identification and Analysis of HSP Genes in L. edodes

A total of 34 HSP genes were identified in the L. edodes genome, including eight small HSP (SHSP), six HSP40, nine HSP60, six HSP70, two HSP90, and three HSP100 genes (Table 1). These 34 HSP genes were further named based on their chromosomal locations. The protein sequence lengths of LeHSPs varied from 153 to 890 amino acids with molecular weights ranging from 17.41 kDa to 98.436 kDa. The theoretical isoelectric point (pI) values of the LeHSPs ranged from 4.4 (HSP70.6) to 10.3 (LeSHSP.1).

3.2. Gene Structure of HSPs

The intron–exon structures were determined by the alignment of the genomic DNA sequences with the coding sequences of HSPs (Figure 1). The number of exons of LeHSP genes varied from 2 to 19 exons (Figure 1, Table 1). LeSHSP genes had 2 to 4 exons. The number of exons in LeHSP40 genes ranged from 4 to 12, with 3 in LeHSP40.5 and 11 introns in LeHSP40.4. The number of exons of LeHSP60s were between 4 (LeHSP60.5) and 11 (LeHSP50.7). The exons number in LeHSP70s ranged from 4 to 11. Two LeHSP90 genes were identified from L. edodes. LeHSP90.1 had 8 exons, and LeHSP90.2 had 19. Three LeHSP100 genes were identified, within two of them (LeHSP100.2 and LeHSP100.3) had 11 exons, and the remainder (LeHSP100.1) had only 8 exons.
Domain analysis can help clarify the functions of the HSPs. The conserved domains of LeSHSP, LeSHP40, LeHSP60, LeHSP70, LeHSP90, and LeHSP100 were identified by Bath Web CD-Search tool (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi (accessed on 18 September 2022)) against the CDD, Pfam, and SMART databases. LeSHSP.1, LeSHSP.3, and LeSHSP.8 had a conserved ACD domain consisting of 80 to 100 amino acids (Table S2). The other LeSHSPs contained the conserved IbpA domain (cl00175). Two conserved domains, DnaJ_C (PF1556) and DnaJ (PF00226), existed in all LeHSP40 proteins. LeHSP60s contained a compound domain Cpn60_TCP1 (PF00118). Two LeHSP70s (LeHSP70.4 and LeHSP70.6) had a conserved domain HSP70 superfamily (cl37948). The other four LeHSP70s contained another HSP70 (PF00012) domain. LeHSP90.1 and LeHSP90.2 contained two conserved domains, HSP90 (PF000183) and ATPase_c. LeHSP100s contained four conserved domains including AAA_2 (PF077724), AAA_lid_9 (PF17871), AAA (PF00004), and ClpB_D2_small (PF10431). In addition to the four conserved domains mentioned above, LeHSP100.1 had another specific domain, Clp_N (PF02861).
The sequence motifs of each LeHSP member were individually predicted by the MEME tool (Table S3). The maximum number of motifs for each protein was set to five. The sequence motifs were conserved within but not between subfamilies, as evidenced by the absence of common motifs in all HSPs (Figure 1). The number of conserved motifs in each HSP subfamily also varied. For instance, of the five motifs predicted in members of LeSHSPs, only one motif was shared by all LeSHSP members. By contrast, all five predicted motifs in the LeHSP100 subfamily were shared by each member.
Except for LeHSP40.4 which possessed four conserved motifs, the other LeHSP40s contained five conserved motifs. LeHSP60s had five conserved motifs, except for LeHSP60.5 which lacked Motif 14. LeHSP70.4 only had two conserved motifs which were Motifs 28 and 29. Compared to LeHSP70.4, LeHSP70.6 had extra Motif 27. Five conserved motifs were found in LeHSP70.1, LeHSP70.2, LeHSP70.3, and LeHSP70.5. LeHSP90s and LHSP100s contained five conserved motifs. Detailed information on motifs is presented in Table S3.

3.3. Phylogenetic Analysis of HSP Proteins

The full-length amino acid sequences of HSP proteins from L. edodes, A. bisporus, P. pulmonarius, P. eryngii, and V. volvaceas were used for phylogenetic tree construction (Figure 2). The 186 HSP sequences were categorized into three distinct groups according to the sequence characteristics. Group I was subdivided into three subgroups representing members from SHSPs, HSP90s, and HSP100s. The SHSP subfamily was an extensive subgroup which included 41 protein sequences. Group II consisted of 43 HSP70 sequences. All members from the HSP40 and HSP60 subfamilies were clustered into Group III. Within Group III, all HSP60 sequences and one HSP40 sequence (Accession number: KAF8652705) from V. volvaceas were clustered into a subgroup (Subgroup I). The other 31 HSP40 sequences were clustered into a separate subgroup (Subgroup II).

3.4. Chromosomal Location and Duplication of LeHSP Genes

Chromosomal location analysis revealed that the 34 LeHSP genes were unevenly distributed on all chromosomes of L. edodes with the exemption of Chromosomes 5 and 9 (Figure 3A). Chromosome 7 contained the largest number (10) of LeHSPs followed by Chromosome 1 with nine genes. Only one gene was positioned on Chromosome 10 (LeHSP70.6). Chromosomes 6 and 8 contained four LeHSP genes. Chromosomes 2, 3, and 4 each possessed two LeHSP genes. No rule was discerned for the distribution of different types of LeHSP genes. Eight LeSHSP genes were, respectively, located on Chromosomes 2 (1 gene), 3 (1 gene), 4 (1 gene), 6 (2 genes), and 7 (3 genes). Six LeHSP40 genes were distributed on Chromosomes 1 (2 genes) and 7 (4 genes). Chromosomes 1, 4, 7, and 8 contained three, one, three, and two LeHSP60 genes, respectively. Six LeHSP70 genes were mapped on Chromosomes 1 (2 genes), 3 (1 gene), 8 (2 genes), and 10 (1 gene). Two LeHSP90 genes were located on Chromosomes 1 and 6. Three LeHSP100 genes were scattered on Chromosomes 1, 2, and 6.
One the basis of the evolutionary relationships among the coding sequences of the 34 LeHSP genes, one and two paralogous pairs were observed in LeHSP100 and LeSHSP genes, respectively (Figure 3B). All pairs were segmental duplications according to Rudd et al. [34].

3.5. The Binding Site of HSF in the Promoters of LeHSP Genes

The HSE 5′-nGAAC-3′ was specifically recognized by the HSF, and the affinity between HSF and HSE was affected by the number of HSEs displayed on the promoter. To further analyze the combining capacity between LeHSF and the promoter of LeHSP, the 1 kb upstream sequences from the initiation codon (ATG) of LeHSP genes were extracted and the number of HSEs was identified (Figure 4). All LeHSP genes’ promoters possessed at least one HSE, suggesting that HSEs might play pivotal roles in response to heat stress. However, the number of HSEs in individual genes highly varied. For instance, the promoters of LeHSP60.2 contained only one HSE, and the promoters of LeSHSP.2 had up to nine HSEs.

3.6. Expression Profile of LeHSP Genes in Response to Heat Stress

To understand the function of LeHSP genes in response to heat stress, the expression patterns of LeHSP genes in both Strains 18 and 18N44 (thermo-resistant) under heat-stress were investigated (Figure 5). However, the FPKM values of the other four LeHSP genes (LeSHSP.6, LeSHSP.8, LeHSP40.3, and LeHSP40.6) were zero or undetectable in the RNA-seq data, so this part of the expression data was discarded during the analysis. Thirty LeHSP genes had FPKM values and were up-expressed in response to heat stress. These LeHSP genes were divided into three groups according to the resulting expression pattern. In Group I, 12 LeHSP genes were highly expressed at 4 or 12 h under heat stress in Strain 18. In Group II, 11 LeHSP genes were clustered and showed a high level of expression at 24 h in Strain 18N44. In Group III, 6 of 7 LeHSP genes were highly expressed at 24 h in Strains 18 and 18N44. However, the expression pattern of LeHSP.4 was unusual as it was highly expressed at 24 h in Strain 18 and low expressed at 24 h in Strain 18N44.
To confirm the expressions of LeHSPs in response to heat stress, RT-qPCR was performed at the same time-points of heat stress in two other strains (Le015 and Le027) by selecting 16 LeHSP genes which had representative expression profiles (Figure 6). According to the RT-qPCR results, most LeHSPs were sensitive to heat stress and showed different expression patterns during the 24 h heat stress. The expressions of LeSHSP.1, LeSHSP.4, and LeSHSP.5 were up-regulated under heat stress and peaked at 4 h. The relative expressions of LeSHSP.4 and LeSHSP.5 increased dramatically by more than 400-fold under heat stress compared to the control (without heat stress treatment) in both Le015 and Le027. An increase in expression was observed for the LeHSP40 genes (LeHSP40.1, LeHSP40.2, and LeHSP40.4) at 4 h of heat stress, and LeHSP40.1 and LeHSP40.2 peaked after 4 h of heat stress. The expressions of LeHSP60s (LeHSP60.3, LeHSP60.4, and LeHSP60.6) were down-regulated in response to heat stress. The expression of LeHSP70.2 was unchanged under heat stress. LeHSP70.3 was up-regulated only after 4 h of heat stress, peaked at 4 h (Le027) or 12 h (Le015), and then dropped subsequently. The expression of LeHSP70.4 was increased after 12 h of heat stress and maintained a high level after that period of heat stress. LeHSP90s were extremely induced by heat stress in Strain Le027, but the expression of LeHSP90.2 in Strain Le015 was repressed under heat stress. The expression levels of LeHSP100s were up-regulated and peaked at 4 h of heat shock, then dropped subsequently. Compared to Le015, the expression of LeHSP100 was significantly higher in Le027.

4. Discussion

HSPs play a crucial role in prokaryotes and eukaryotes in response to various stress conditions, including heat stress [35]. To date, the identifications and functions of HSPs have been well investigated in animals and plants but largely lagged in edible fungi. In the present study, 34 HSP genes, with 8 LeSHSPs, 6 LeHSP40s, 9 LeHSP60s, 2 LeHSP90s, and 3 LeHSP100s, were identified from the entire genome of L. edodes (Le(Bin) 0899 ss11). Previous research [36] identified 11 LeSHSPs, 1 LeHSP60, 13 LeHSP70s, 2 LeHSP90s, 1 LeHSP98 and 1 LeHSP104 in Strains YS606 and YS3357 of L. edodes through RNA-seq data. We believed that two reasons may be responsible for the different number of LeHSPs in the two studies. Firstly, the variation in the number of HSPs found in L. edodes may be caused by the different genome sources used in these two studies. Secondly, some LeHSP genes expressed under the limit of detection using RNA sequencing and these genes were not analyzed in Wang’s study.
Compared to other species, L. edodes had fewer HSP genes. A total of 64, 79, and 753 HSP genes were identified in A. thaliana [16], S. italica [37], and T. aestivum [12], respectively. The low number of HSP genes in L. edodes may arise from its small genome. For instance, A. thaliana, whose genome size is 125 Mb [38], possesses few HSP genes, but T. aestivum with 6.3 Gb of genome size [12] has numerous HSP genes.
In this study, we identified 34 HSP genes in L. edodes, then analyzed the structure, chromosomal location, gene duplication, binding site of HSF, and expression pattern under heat stress. Variation occurred in terms of the number of introns and exons, but the motifs were conserved in different classifications of HSP proteins. The phylogenetic and motif analysis of HSPs in L. edodes confirm that members from the same sub-group share the same features of conserved domains, indicating they possess a similar biological function in response to different stresses. The different numbers of HSEs laid in the LeHSPs’ promoter revealed that all LeHSP genes have the potential to be recognized by HSFs [39]. Gene duplications play a key role in expansions of gene family members in many species, such as tomato [40], pepper [41], and cucumber [42]. In this study, a total of three gene duplications were found from the L.edodes genome and these duplications belonged to segmental duplication events, implying that the three duplication events contributed to the expansion of the LeHSP family.
According to the expression patterns, most LeHSPs were expressed in response to heat stress. In Strains 18 and 18N44, the transcript level of all LeHSP genes showed an increase after 4 h to 24 h of heat stress at 37 °C (Figure 5). In the thermo-tolerant Strain 18N44, the expression level of LeHSPs, including LeHSP40.2, LeHSP60s, LeHSP70.2 and LeHSP70.6, were significantly up-regulated under heat stress. In previous study, they confirmed that the relative expression of LeHSP100, LeHSP90, and LeHSP60 is far higher in 18N44 than that in strain 18 [43]. The RT-qPCR results reveal that except for LeHSP60 (LeHSP60.3, LeHSP60.4, and LeHSP60.6) whose expression showed a decrease in response to heat stress, the expression levels of the 13 other LeHSPs (LeSHSP, LeHSP40, LeHSP70, LeHSP90 and LeHSP100) were higher under heat stress than under the normal condition (25 °C) (Figure 6). Given the results, we deduced that the combination of different types and numbers of up-regulated HSP expression can contribute to thermotolerance. Previous studies have successfully examined the function of HSP40 (LeDnaJ07) in L. edodes and found that overexpression of LeDanJ07 could enhance mycelial resistance in response to heat stress in YS55 [2]. Identification and expression of DnaJ protein (HSP40) was well studied in L. edodes, and the expression of LeDnaJ genes was induced by cadmium, Trichoderma atroviride, and heat stress, thereby indicating that LeDnaJ proteins were involved in multiple stresses [44]. In addition, the proteome analysis exhibited that the contents of HSP40 (DnaJ), HSP70, and HSP98 were significantly up-regulated in the thermotolerant Strain S606 after heat stress [36]. Up-regulation of members of various HSP genes in response to heat stress are reported in other species [45,46]. For example, 2-fold to 20-fold increases in expressions were detected for 11 members of the A. thaliana HSP70 family under heat stress [16]. In finger millet, a HSP gene, HSP17.8, showed a 40-fold up-regulation under heat stress, was cloned, and proved to have major contributions to thermo-tolerance in the crop [47]. The results also confirmed the association of HSP genes (HSP70, HSP21, and SHSP21) with the high temperature tolerance of Agasicles hygrophila [48]. The expression of HSP genes were regulated by the HSFs. When cells were exposed to abiotic stress, HSFs promptly bind to the promoter of HSPs to regulate gene transcription. Then, the newly synthesized HSPs act as molecular chaperone to prevent the substrate protein from aggregating irreversibly [49,50]. Subsequently, the different signaling pathways related to enhancing tolerance to heat stress will be activated by HSPs. In this study, we confirmed that thermotolerance was associated with up-regulated HSP expression, but the mechanism by which HSPs enhance thermo-resistance in L. edodes requires further investigation.

5. Conclusions

In this study, 34 HSPs were identified and characterized after scanning for the L. edodes genome. LeHSPs were classified according to their molecular weight and verified by phylogenetic tree construction and conserved motifs analysis. The gene structure, distribution, motif analysis, duplication, and phylogenetic analysis revealed a clear evolutionary history for the HSP family in L. edodes. Three gene duplication events were found in the LeHSP gene family, and they belong to segmental duplications. The different numbers of HSEs in the LeHSPs’ promoters indicated that HSEs could be induced by heat stress.
The expression analysis indicated that the combination of different types and numbers of up-regulated LeHsp genes mainly contributed to thermo-tolerance. This study not only lays a scientific foundation for in-depth understanding of the regulation mechanism of HSP genes, but also provides a novel direction for breeding a mycelium of L. edodes that is tolerant to high temperature.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cimb45010041/s1, Table S1. Primer pairs used for gene expression. Table S2. Information of Domains in HSPs. Table S3. Amino sequence of motifs in HSPs.

Author Contributions

X.Z. and Q.Z. were involved in the design and coordination of the study, conceptualized the experiments, analyzed data, and produced the manuscript. K.Y., R.F. and R.M. participated in the design and coordination of the study, conceived of the experiments, analyzed data, and contributed to manuscript writing. J.L. conceived of experiments and contributed to manuscript writing. L.C., Y.N. and W.L. performed experiments and analyzed data related to the RT-qPCR. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Scientific and Technological Innovation Talents of Sichuan Province (Grant No. 2022JDRC0034), the Gansu Provincial Youth Science and Technology Fund Program (Grant No. 20JR5RA064), the Central Public-interest Scientific Institution Basal Research Fund (Grant No. S2022007), and the Local Financial Funds of the National Agricultural Science and Technology Center, Chengdu (Grant No. NASC2021KR06).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors acknowledge Xiaojiang Guo for his input in bioinformatics analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Finimundy, T.C.; Dillon, A.J.P.; Henriques, J.A.P.; Ely, M.R. A Review on General Nutritional Compounds and Pharmacological Properties of the Lentinula edodes Mushroom. Food Nutr. Sci. 2014, 05, 1095–1105. [Google Scholar] [CrossRef] [Green Version]
  2. Wang, G.; Luo, Y.; Wang, C.; Zhou, Y.; Mou, C.; Kang, H.; Xiao, Y.; Bian, Y.; Gong, Y.H. Hsp40 Protein LeDnaJ07 Enhances the Thermotolerance of Lentinula edodes and Regulates IAA Biosynthesis by Interacting LetrpE. Front. Microbiol. 2020, 11, 707. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Kotak, S.; Larkindale, J.; Lee, U.; von Koskull-Doring, P.; Vierling, E.; Scharf, K.D. Complexity of the heat stress response in plants. Curr. Opin. Plant Biol. 2007, 10, 310–316. [Google Scholar] [CrossRef] [PubMed]
  4. Luo, L.; Zhang, S.; Wu, J.; Sun, X.; Ma, A. Heat stress in macrofungi: Effects and response mechanisms. Appl. Microbiol. Biotechnol. 2021, 105, 7567–7576. [Google Scholar] [CrossRef] [PubMed]
  5. Roy, A.; Tamuli, R. Heat shock proteins and the calcineurin-crz1 signaling regulate stress responses in fungi. Arch. Microbiol. 2022, 204, 240. [Google Scholar] [CrossRef] [PubMed]
  6. Ritossa, F. A new puffing pattern induced by temperature shock and DNP in drosophila. Experientia 1962, 18, 571–573. [Google Scholar] [CrossRef]
  7. Wang, W.; Vinocur, B.; Shoseyov, O.; Altman, A. Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response. Trends Plant Sci. 2004, 9, 244–252. [Google Scholar] [CrossRef]
  8. Li, Z.; Long, R.; Zhang, T.; Yang, Q.; Kang, J. Research Progress on Plant Heat Shock Protein. Biotechnol. Bull. 2016, 32, 7–13. [Google Scholar] [CrossRef]
  9. Xu, J.; Xue, C.; Xue, D.; Zhao, J.; Gai, J.; Guo, N.; Xing, H. Overexpression of GmHsp90s, a heat shock protein 90 (Hsp90) gene family cloning from soybean, decrease damage of abiotic stresses in Arabidopsis thaliana. PLoS ONE 2013, 8, e69810. [Google Scholar] [CrossRef] [PubMed]
  10. Park, J.C.; Kim, D.H.; Lee, Y.; Lee, M.C.; Kim, T.K.; Yim, J.H.; Lee, J.S. Genome-wide identification and structural analysis of heat shock protein gene families in the marine rotifer Brachionus spp.: Potential application in molecular ecotoxicology. Comp. Biochem. Physiol. Part D Genom. Proteom. 2020, 36, 100749. [Google Scholar] [CrossRef]
  11. Jaeger, A.M.; Makley, L.N.; Gestwicki, J.E.; Thiele, D.J. Genomic heat shock element sequences drive cooperative human heat shock factor 1 DNA binding and selectivity. J. Biol. Chem. 2014, 289, 30459–30469. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Kumar, A.; Sharma, S.; Chunduri, V.; Kaur, A.; Kaur, S.; Malhotra, N.; Kumar, A.; Kapoor, P.; Kumari, A.; Kaur, J.; et al. Genome-wide Identification and Characterization of Heat Shock Protein Family Reveals Role in Development and Stress Conditions in Triticum aestivum L. Sci. Rep. 2020, 10, 7858. [Google Scholar] [CrossRef] [PubMed]
  13. Xu, S.; Sun, W.; Tian, J.; Wang, C. Plant Heat Shock Protein HSP100/ClpB and Its Applications in Improvement of Heat and Cold Resistances in Plants. Plant Physiol. Commun. 2008, 44, 804–810. [Google Scholar] [CrossRef] [Green Version]
  14. Lee, J.; An, Y.S.; Kim, M.R.; Kim, Y.A.; Lee, J.K.; Hwang, C.S.; Chung, E.; Park, I.C.; Yi, J.Y. Heat Shock Protein 90 Regulates Subcellular Localization of Smads in Mv1Lu Cells. J. Cell. Biochem. 2016, 117, 230–238. [Google Scholar] [CrossRef]
  15. Hu, W.; Hu, G.; Han, B. Genome-wide survey and expression profiling of heat shock proteins and heat shock factors revealed overlapped and stress specific response under abiotic stresses in rice. Plant Sci. 2009, 176, 583–590. [Google Scholar] [CrossRef]
  16. Sung, D.Y.; Vierling, E.; Guy, C.L. Comprehensive expression profile analysis of the Arabidopsis Hsp70 gene family. Plant Physiol. 2001, 126, 789–800. [Google Scholar] [CrossRef] [Green Version]
  17. Wang, Y.; Lin, S.; Song, Q.; Li, K.; Tao, H.; Huang, J.; Chen, X.; Que, S.; He, H. Genome-wide identification of heat shock proteins (Hsps) and Hsp interactors in rice: Hsp70s as a case study. BMC Genom. 2014, 15, 344. [Google Scholar] [CrossRef] [Green Version]
  18. Zhang, L.; Zhao, H.K.; Dong, Q.L.; Zhang, Y.Y.; Wang, Y.M.; Li, H.Y.; Xing, G.J.; Li, Q.Y.; Dong, Y.S. Genome-wide analysis and expression profiling under heat and drought treatments of HSP70 gene family in soybean (Glycine max L.). Front. Plant Sci. 2015, 6, 773. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, J.; Shen, N.; Li, C.; Xiang, X.; Liu, G.; Gui, Y.; Patev, S.; Hibbett, D.S.; Barry, K.; Andreopoulos, W.; et al. Population genomics provides insights into the genetic basis of adaptive evolution in the mushroom-forming fungus Lentinula edodes. J. Adv. Res. 2021, 38, 91–106. [Google Scholar] [CrossRef]
  20. El-Gebali, S.; Mistry, J.; Bateman, A.; Eddy, S.R.; Luciani, A.; Potter, S.C.; Qureshi, M.; Richardson, L.J.; Salazar, G.A.; Smart, A.; et al. The Pfam protein families database in 2019. Nucleic Acids Res. 2019, 47, 427–432. [Google Scholar] [CrossRef]
  21. Chen, C.; Chen, H.; Zhang, Y.; R. Thomas, H.; H. Frank, M.; He, Y.; Xia, R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  22. Bailey, T.L.; Williams, N.; Misleh, C.; Li, W.W. MEME: Discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 2006, 34, 369–373. [Google Scholar] [CrossRef] [PubMed]
  23. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C.; et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Panzade, K.P.; Kale, S.S.; Chavan, N.R.; Hatzade, B. Genome-wide analysis of Hsp70 and Hsp100 gene families in Ziziphus jujuba. Cell Stress Chaperones 2020, 26, 341–353. [Google Scholar] [CrossRef]
  25. Letunic, I.; Bork, P. Interactive Tree Of Life (iTOL) v5: An online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021, 49, W293–W296. [Google Scholar] [CrossRef]
  26. Wang, L.; Zhao, Y.; Zhang, B.; Chen, M. Breeding thermo-tolerant strains of Lentinula edodes by UV induced protoplast mutagenesis. Microbiol. China 2014, 41, 1350–1357. [Google Scholar] [CrossRef]
  27. Zhao, X.; Yang, H.; Chen, M.; Song, X.; Yu, C.; Zhao, Y.; Wu, Y. Reference Gene Selection for Quantitative Real-Time PCR of Mycelia from Lentinula edodes under High-Temperature Stress. Biomed. Res. Int. 2018, 2018, 1670328. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Trapnell, C.; Williams, B.A.; Pertea, G.; Mortazavi, A.; Kwan, G.; van Baren, M.J.; Salzberg, S.L.; Wold, B.J.; Pachter, L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 2010, 28, 511–515. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [Green Version]
  30. Roberts, A.; Pachter, L. Streaming fragment assignment for real-time analysis of sequencing experiments. Nat. Methods 2013, 10, 71–73. [Google Scholar] [CrossRef]
  31. Tao, Y.; van Peer, A.F.; Huang, Q.; Shao, Y.; Zhang, L.; Xie, B.; Jiang, Y.; Zhu, J.; Xie, B. Identification of novel and robust internal control genes from Volvariella volvacea that are suitable for RT-qPCR in filamentous fungi. Sci. Rep. 2016, 6, 29236. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Glare, E.M.; Divjak, M.; Bailey, M.J.; Walters, E.H. β-Actin and GAPDH housekeeping gene expression in asthmatic airways is variable and not suitable for normalising mRNA levels. Thorax 2002, 57, 765–770. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Pfaffl, M.W. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29, 2002–2007. [Google Scholar] [CrossRef] [PubMed]
  34. Rudd, M.K.; Keene, J.; Bunke, B.; Kaminsky, E.B.; Adam, M.P.; Mulle, J.G.; Ledbetter, D.H.; Martin, C.L. Segmental duplications mediate novel, clinically relevant chromosome rearrangements. Hum. Mol. Genet 2009, 18, 2957–2962. [Google Scholar] [CrossRef] [Green Version]
  35. Antomio, D.M. Heat shock proteins: Facts, thoughts, and dreams. Shock 1999, 11, 1–12. [Google Scholar] [CrossRef]
  36. Wang, G.Z.; Ma, C.J.; Luo, Y.; Zhou, S.S.; Zhou, Y.; Ma, X.L.; Cai, Y.L.; Yu, J.J.; Bian, Y.B.; Gong, Y.H. Proteome and Transcriptome Reveal Involvement of Heat Shock Proteins and Indoleacetic Acid Metabolism Process in Lentinula Edodes Thermotolerance. Cell Physiol. Biochem. 2018, 50, 1617–1637. [Google Scholar] [CrossRef]
  37. Singh, R.K.; Jaishankar, J.; Muthamilarasan, M.; Shweta, S.; Dangi, A.; Prasad, M. Genome-wide analysis of heat shock proteins in C4 model, foxtail millet identifies potential candidates for crop improvement under abiotic stress. Sci. Rep. 2016, 6, 32641. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Scharf, K.D.; Siddique, M.; Vierling, E. The expanding family of Arabidopsis thaliana small heat stress proteins and a new family of proteins containing alpha-crystallin domains (Acd proteins). Cell Stress Chaperones 2001, 6, 225–237. [Google Scholar] [CrossRef]
  39. Santoro, N.; Johansson, N.; Thiele, D.J. Heat Shock Element Architecture Is an Important Determinant in the Temperature and Transactivation Domain Requirements for Heat Shock Transcription Factor. Mol. Cell. Biol. 1998, 18, 6340–6352. [Google Scholar] [CrossRef] [Green Version]
  40. Yu, J.; Cheng, Y.; Feng, K.; Ruan, M.; Ye, Q.; Wang, R.; Li, Z.; Zhou, G.; Yao, Z.; Yang, Y.; et al. Genome-Wide Identification and Expression Profiling of Tomato Hsp20 Gene Family in Response to Biotic and Abiotic Stresses. Front. Plant Sci. 2016, 7, 1215. [Google Scholar] [CrossRef] [PubMed]
  41. Guo, M.; Liu, J.H.; Lu, J.P.; Zhai, Y.F.; Wang, H.; Gong, Z.H.; Wang, S.B.; Lu, M.H. Genome-wide analysis of the CaHsp20 gene family in pepper: Comprehensive sequence and expression profile analysis under heat stress. Front. Plant Sci. 2015, 6, 806. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Wang, X.; Zheng, Y.; Chen, B.; Zhi, C.; Qiao, L.; Liu, C.; Pan, Y.; Cheng, Z. Genome-wide identification of small heat shock protein (HSP20) homologs in three cucurbit species and the expression profiles of CsHSP20s under several abiotic stresses. Int J. Biol. Macromol. 2021, 190, 827–836. [Google Scholar] [CrossRef] [PubMed]
  43. Miaomiao, X. The Gene Expression and DFunction Study of Hydrophobin Heat Shock Proteins From Lentinula edodes; Nanjing Agricultural University: Nanjing, China, 2016. [Google Scholar]
  44. Wang, G.Z.; Luo, Y.; Chen, W.; Gong, Y.H.; Zhou, Y.; Bian, Y.B. The DnaJ Gene Family in Shiitake Culinary-Medicinal Mushroom, Lentinus edodes (Agaricomycetes): Comprehensive Identification, Characterization, and Expression Profiles under Different Conditions. Int. J. Med. Mushrooms 2019, 21, 909–919. [Google Scholar] [CrossRef]
  45. Alshameri, A.; Al-Qurainy, F.; Gaafar, A.R.; Khan, S.; Nadeem, M.; Alansi, S. Identification of Heat-Responsive Genes in Guar [Cyamopsis tetragonoloba (L.) Taub]. Int. J. Genomics 2020, 2020, 3126592. [Google Scholar] [CrossRef]
  46. Zhao, W.; Ma, Z.; Liu, S.; Yang, W.; Ma, J. Transcriptome Profiling Reveals Potential Genes and Pathways Supporting Ananas comosus L. Merr’s High Temperature Stress Tolerance. Trop. Plant Biol. 2021, 14, 132–142. [Google Scholar] [CrossRef]
  47. Chopperla, R.; Singh, S.; Tomar, R.S.; Mohanty, S.; Khan, S.; Reddy, N.; Padaria, J.C.; Solanke, A.U. Isolation and allelic characterization of finger millet (Eleusine coracana L.) small heat shock protein EcHSP17.8 for stress tolerance. Indian J. Genet. Plant Breed. 2017, 78, 95. [Google Scholar] [CrossRef]
  48. Jin, J.; Zhao, M.; Wang, Y.; Zhou, Z.; Wan, F.; Guo, J. Induced Thermotolerance and Expression of Three Key Hsp Genes (Hsp70, Hsp21, and sHsp21) and Their Roles in the High Temperature Tolerance of Agasicles hygrophila. Front. Physiol. 2019, 10, 1593. [Google Scholar] [CrossRef] [Green Version]
  49. Mani, N.; Ramakrishna, K.; Suguna, K. Characterization of rice small heat shock proteins targeted to different cellular organelles. Cell Stress Chaperones 2015, 20, 451–460. [Google Scholar] [CrossRef] [Green Version]
  50. Morimoto, R.I. Cells in stress: Transcriptional activation of heat shock genes. Science 1993, 259, 1409–1410. [Google Scholar] [CrossRef]
Cimb 45 00041 g001 550
Figure 1. Phylogenetic relationship, gene structure, and conserved motif analysis of LeHSP genes. (A) Phylogenetic tree of 34 LeHSP proteins. Azure, brown, red, purple, blue, and green colors represent the LeSHSP, LeHSP40, LeHSP60, LeHSP90, LeHSP100, and LeHSP70 genes, respectively. (B) Exon/intron organization of LeHSP genes. Green boxes represent exons and black lines with the same length represent introns. (C) Distributions of conserved motifs in LeHSPs. The conserved motifs of 34 proteins were identified and visualized by MEME.
Figure 1. Phylogenetic relationship, gene structure, and conserved motif analysis of LeHSP genes. (A) Phylogenetic tree of 34 LeHSP proteins. Azure, brown, red, purple, blue, and green colors represent the LeSHSP, LeHSP40, LeHSP60, LeHSP90, LeHSP100, and LeHSP70 genes, respectively. (B) Exon/intron organization of LeHSP genes. Green boxes represent exons and black lines with the same length represent introns. (C) Distributions of conserved motifs in LeHSPs. The conserved motifs of 34 proteins were identified and visualized by MEME.
Cimb 45 00041 g001
Cimb 45 00041 g002 550
Figure 2. Phylogenetic tree of HSPs from L. edodes, A. bisporus, P. pulmonarius, P. eryngii, and V. volvaceas. The phylogenetic tree was constructed using the UPGMA method and Jukes-Cantor model with the Geneious primer software. Different subfamilies are shaded with different colors.
Figure 2. Phylogenetic tree of HSPs from L. edodes, A. bisporus, P. pulmonarius, P. eryngii, and V. volvaceas. The phylogenetic tree was constructed using the UPGMA method and Jukes-Cantor model with the Geneious primer software. Different subfamilies are shaded with different colors.
Cimb 45 00041 g002
Cimb 45 00041 g003 550
Figure 3. Chromosomal location and gene duplication of LeHSP genes. (A) 34 LeHSP genes were mapped on 8 chromosomes. Chromosome numbers are indicated at the top of each bar. The gene names are labelled in different colors. The scale on the left is in megabases. (B) Relationship of LeHSPs with one another in different subfamilies. The white color represents shared zero percent identity, and deep blue indicates shared 100% identity.
Figure 3. Chromosomal location and gene duplication of LeHSP genes. (A) 34 LeHSP genes were mapped on 8 chromosomes. Chromosome numbers are indicated at the top of each bar. The gene names are labelled in different colors. The scale on the left is in megabases. (B) Relationship of LeHSPs with one another in different subfamilies. The white color represents shared zero percent identity, and deep blue indicates shared 100% identity.
Cimb 45 00041 g003
Cimb 45 00041 g004 550
Figure 4. Predicted binding site of HSF in LeHSP promoters. Promoter sequences (−1000 bp) of 34 LeHSP genes were used to predict the HSEs (5−nGAAC−3′) using TBtools. The location of HSEs can be inferred according to the scale at the bottom.
Figure 4. Predicted binding site of HSF in LeHSP promoters. Promoter sequences (−1000 bp) of 34 LeHSP genes were used to predict the HSEs (5−nGAAC−3′) using TBtools. The location of HSEs can be inferred according to the scale at the bottom.
Cimb 45 00041 g004
Cimb 45 00041 g005 550
Figure 5. Heat map of the expression profiles of LeHSP genes in response to heat stress. The FPKM values of LeHSP genes were Z-scores normalized and centralized by a unit variance scaling method.
Figure 5. Heat map of the expression profiles of LeHSP genes in response to heat stress. The FPKM values of LeHSP genes were Z-scores normalized and centralized by a unit variance scaling method.
Cimb 45 00041 g005
Cimb 45 00041 g006 550
Figure 6. Expression profiles of LeHSP genes under heat stress in the strains of Le015 and Le027. Transcript level of different LeHSP genes (relative to geometric mean of LeTUB and LeUBI) in response to heat stress in Le015 (red) and Le027 (blue). X-axes showed heat stress time. Y-axes indicated scales of relative expression level. One way ANOVA was used for statistical analysis. p-values are indicated as follows ** p < 0.01.
Figure 6. Expression profiles of LeHSP genes under heat stress in the strains of Le015 and Le027. Transcript level of different LeHSP genes (relative to geometric mean of LeTUB and LeUBI) in response to heat stress in Le015 (red) and Le027 (blue). X-axes showed heat stress time. Y-axes indicated scales of relative expression level. One way ANOVA was used for statistical analysis. p-values are indicated as follows ** p < 0.01.
Cimb 45 00041 g006
Table
Table 1. Molecular features of Hsp genes in L. edodes.
Table 1. Molecular features of Hsp genes in L. edodes.
NameProtein IDExonIntronLocationChrAAMW(kDa)PI
HSP100LeHSP100.1XP_046091847.1986279738–6282864189098.4365.48
LeHSP100.2XP_046082612.111102522518–2526440277184.9595.59
LeHSP100.3XP_046080071.111103743879–3748058674983.1085.54
HSP90LeHSP90.1XP_046088585.1872407716–2410082170680.1854.62
LeHSP90.2XP_046090891.119181815650–1819056581791.7124.29
HSP70LeHSP70.1XP_046088664.1542657999–2660171165170.9624.78
LeHSP70.2XP_046088714.1982800610–28029171613675.24
LeHSP70.3XP_046086623.111104172952–4175420362768.0465.46
LeHSP70.4XP_046083696.1981532467–1535557888697.5975.52
LeHSP70.5XP_046079422.1981824230–1826698867873.7424.74
LeHSP70.6XP_046090274.143288215–2903461063166.3244.41
HSP60LeHSP60.1XP_046082828.1981726567–1728602153957.887.33
LeHSP60.2XP_046088531.1872303108–2305148155459.7226.24
LeHSP60.3XP_046092210.1654944807–4946672156060.8815.98
LeHSP60.4XP_046084390.1981792382–1794487455059.8336.63
LeHSP60.5XP_046084477.1431298157–1300114759862.5925.26
LeHSP60.6XP_046084950.1871721053–1723408757562.3775.79
LeHSP60.7XP_046088311.111102236637–2238719754259.2435.29
LeHSP60.8XP_046083662.1871633717–1635690856160.1385.91
LeHSP60.9XP_046079390.11091750215–1752273852456.4625.16
HSP40LeHSP40.1XP_046088819.1653097898–3099296139843.5476.44
LeHSP40.2XP_046091483.1657116980–7118740149352.3949.53
LeHSP40.3XP_046084554.111101053764–1055625742246.4174.59
LeHSP40.4XP_046084878.112111936586–1938308736841.5016.18
LeHSP40.5XP_046079519.1432621252–2622535737239.8579.66
LeHSP40.6XP_046082726.1542888022–2889557742948.31710.02
SHSPLeSHSP.1XP_046080632.1433183573–3184343219021.26110.33
LeSHSP.2XP_046088441.1325037923–5037741315417.6716.54
LeSHSP.3XP_046084164.1211241615–1242240418820.374.83
LeSHSP.4XP_046088091.1322943783–2944370615617.6665.65
LeSHSP.5XP_046088092.1322944818–2945400615517.4095.45
LeSHSP.6XP_046082167.132382364–382938715317.4915.57
LeSHSP.7XP_046082168.143383366–383941716218.4926.45
LeSHSP.8XP_046080378.1322094217–2095081725429.0545.22
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

Zhao, X.; Yin, K.; Feng, R.; Miao, R.; Lin, J.; Cao, L.; Ni, Y.; Li, W.; Zhang, Q. Genome-Wide Identification and Analysis of the Heat-Shock Protein Gene in L. edodes and Expression Pattern Analysis under Heat Shock. Curr. Issues Mol. Biol. 2023, 45, 614-627. https://doi.org/10.3390/cimb45010041

AMA Style

Zhao X, Yin K, Feng R, Miao R, Lin J, Cao L, Ni Y, Li W, Zhang Q. Genome-Wide Identification and Analysis of the Heat-Shock Protein Gene in L. edodes and Expression Pattern Analysis under Heat Shock. Current Issues in Molecular Biology. 2023; 45(1):614-627. https://doi.org/10.3390/cimb45010041

Chicago/Turabian Style

Zhao, Xu, Kaiyong Yin, Rencai Feng, Renyun Miao, Junbin Lin, Luping Cao, Yanqing Ni, Wensheng Li, and Qin Zhang. 2023. "Genome-Wide Identification and Analysis of the Heat-Shock Protein Gene in L. edodes and Expression Pattern Analysis under Heat Shock" Current Issues in Molecular Biology 45, no. 1: 614-627. https://doi.org/10.3390/cimb45010041

APA Style

Zhao, X., Yin, K., Feng, R., Miao, R., Lin, J., Cao, L., Ni, Y., Li, W., & Zhang, Q. (2023). Genome-Wide Identification and Analysis of the Heat-Shock Protein Gene in L. edodes and Expression Pattern Analysis under Heat Shock. Current Issues in Molecular Biology, 45(1), 614-627. https://doi.org/10.3390/cimb45010041

Article Metrics

Back to TopTop