Transcriptome analysis reveals key genes in response to high-temperature stress in Rhododendron molle

Abstract

Rhododendron molle is a deciduous rhododendron, a high-altitude plant prized for its medicinal and ornamental properties. A major challenge when introducing this plant to lower altitudes is understanding its response to high-temperature stress. Using transcriptome analysis, this study examined leaves under varying temperatures, identifying 344,593 transcripts, 124,901 Unigenes, and 12,089 differentially expressed genes (DEGs) at 36 °C high-temperature stress (ST). At 42 °C high-temperature stress (SY), 12,032 DEGs were found, indicating a significant impact of temperature on gene expression. A Gene Ontology analysis (GO) revealed that these DEGs are mostly involved in stress response, catalytic activity, binding, transporter activity, and immune processes. A Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis highlighted enrichment in pathways such as plant-pathogen interaction and spliceosomes, suggesting their key roles in the temperature stress response. Key genes such as Brassinosteroid-Insensitive 1-Associated Receptor Kinase (BAK1) and Pathogenesis-Related Gene Transcription Activation Factor (PTI6) were upregulated under ST stress, whereas heat shock proteins (HSP83A) and HSP90-1 were downregulated under SY. These findings offer insights into the molecular response of Rhododendron molle to high temperature, aiding further research in this area and potentially improving the plant’s cultivation and application.

Full Article

Transcriptome Analysis Reveals Key Genes in Response to High-Temperature Stress in Rhododendron molle

Linshi Wu,^a,b,* Yan Liu,^a,b,* Xinyun Liu,^a,b Qiaoyun Li,^a,b Xinyu Yi,^a,b Chan Chen,^a,b Ling Wang,^a,b and Juyang Liao ^a,b

Rhododendron molle is a deciduous rhododendron, a high-altitude plant prized for its medicinal and ornamental properties. A major challenge when introducing this plant to lower altitudes is understanding its response to high-temperature stress. Using transcriptome analysis, this study examined leaves under varying temperatures, identifying 344,593 transcripts, 124,901 Unigenes, and 12,089 differentially expressed genes (DEGs) at 36 °C high-temperature stress (ST). At 42 °C high-temperature stress (SY), 12,032 DEGs were found, indicating a significant impact of temperature on gene expression. A Gene Ontology analysis (GO) revealed that these DEGs are mostly involved in stress response, catalytic activity, binding, transporter activity, and immune processes. A Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis highlighted enrichment in pathways such as plant-pathogen interaction and spliceosomes, suggesting their key roles in the temperature stress response. Key genes such as Brassinosteroid-Insensitive 1-Associated Receptor Kinase (BAK1) and Pathogenesis-Related Gene Transcription Activation Factor (PTI6) were upregulated under ST stress, whereas heat shock proteins (HSP83A) and HSP90-1 were downregulated under SY. These findings offer insights into the molecular response of Rhododendron molle to high temperature, aiding further research in this area and potentially improving the plant’s cultivation and application.

DOI: 10.15376/biores.19.4.8238-8256

Keywords: Rhododendron molle; High temperature stress; Transcriptome analysis; Differentially expressed genes (DEGs)

Contact information: a: Hunan Botanical Garden, Box 4100004, Changsha P.R. China; b: Rhododendron Engineering Technology Research Center, Box 4100004, Changsha P.R. China; *These authors contributed equally to this work; Corresponding authors: Juyang Liao liaojuyang@163.com (Hunan Botanical Garden, Changsha, PR China)

INTRODUCTION

Rhododendron molle, distinctively known as the “yellow azalea” or “noisy sheep flower,” holds the singular distinction of being the only species native to China within the subgenus Pentanthera of the Rhododendron family (Ericaceae). Celebrated for its exceptional ornamental appeal, it plays a crucial role in enriching the color diversity of azaleas (Ureshino et al. 1998). This plant predominantly flourishes at altitudes of around one thousand meters, with its main habitat spread across central and southern China and a smaller presence in the southwestern region. Notably, Rhododendron molle is currently recognized as an endangered species (Wang et al. 2010), underscoring the need for conservation efforts to protect this unique and valuable botanical treasure.

Rhododendron molle, a plant with significant medicinal qualities, is enriched with compounds such as lignans in its flowers and flavonoids, rhododendron toxins, and methyl-lichenic acid in its leaves. Lignans are an important kind of natural phytoestrogen with various biological activities, including antioxidant, antiviral, and anti-tumor activities. The lignans in rhododendrons have little effect on animals, but the gray anotoxins have a strong insecticidal effect and have been widely used as a pesticide. When applied via contact or fumigation, they poison the stomach, interfere with feeding, and inhibit the growth and development of many kinds of pests. Rhododendron molle’s other constituents endow it with diverse therapeutic benefits. Traditionally, it has been used for its effectiveness in dispelling wind, alleviating pain, and reducing dampness, and it also serves as an anesthetic. It is particularly effective in treating conditions such as fractures, rheumatism, and persistent skin eczema (Liu et al. 2007; Cheng et al. 2011; Zhi et al. 2013; Li et al. 2018). Intriguingly, certain toxins in Rhododendron molle, especially those derived from the “noisy sheep flower,” have shown remarkable agricultural value. They demonstrate potent stomach-poisoning, contact-killing, and fumigation effects against various pests (Zhu et al. 2023). Notably, research has highlighted that these toxins are effective in deterring feeding in the larvae of pests like Leptinotarsa decemlineata and Spodoptera frugiperda, underscoring the plant’s potential in natural pest control (Shin and Yu 1993). Aside from its medicinal properties, Rhododendron molle also stands as a pivotal resource in azalea breeding. Azalea breeders have long recognized its excellent compatibility with other azalea species, making it an ideal candidate for hybridization (Geng et al. 2019). Many of its traits, typically recessive in the first generation, render it a superior parent for breeding programs. Despite its critical role in azalea cultivation, the challenges presented by Rhododendron molle such as its low seed germination rate and the difficult rooting and low survival rates of cuttings, have hindered its use in urban landscaping (Meijón et al. 2009). Consequently, intensifying research to improve Rhododendron molle cultivation techniques and promote their broader implementation is crucial for diversifying urban tree species and holds significant implications for agricultural production and medical research (Zhou and Zhu 2020).

Rhododendron molle grows under shrubs or mixed forests on mountain ridges in mountainous grasslands or hilly areas at an altitude of 1000 meters, with the most suitable growth temperature being 15 to 20 ℃; these conditions adversely affect the introduction and garden cultivation of Rhododendron molle (Arisumi et al. 1989). Additionally, global climate change has led to more frequent occurrences of hot weather and water shortages across the nation, impacting the normal growth of plants. This situation of high temperatures is expected to worsen in the future (Schlenker and Maximilian 2018; Yuan et al. 2019). Therefore, it is crucial to understand and explore the molecular regulatory mechanisms that enable Rhododendron molle to withstand high-temperature stress, which is particularly important for breeding purposes.

RNA-Seq can be used to comprehensively obtain the sequence and expression information of the transcriptome of a particular organism under different environmental conditions and analyze the differential expression of genes at the RNA level to identify key DEGs. This, in turn, allows an in-depth study of the regulatory mechanisms and functions of differential genes in different environments (Yan et al. 2004). For example, 18 genes encoding heat stress proteins (HSPs) were identified in Paeonia suffruticosa by means of transcriptome analysis, and all of them showed significant up-regulation under high-temperature treatment. A transcriptome analysis of Narcissus tazetta var. chinensis in response to high-temperature stress identified 96 known miRNAs, 25 novel miRNAs and 52,286 predicted siRNAs that may play important roles in response to heat stress in daffodils (Yan et al. 2004). In this study, we analyzed the transcripts of Rhododendron molle under high-temperature stress using RNA-seq technology to obtain the differentially expressed genes. Our results can provide a reference for the later molecular breeding of Rhododendron molle variants that are robust against high temperatures.

EXPERIMENTAL

Study Sites and Plant Materials

The experimental materials were kept in Hunan Provincial Botanical Garden (113°01′30″E, 28°06′40″N) Two-year-old cuttings of Rhododendron truncatulum were transplanted into pots (h≥20cm), with one cutting planted in each pot. A total of 2 high-temperature treatments (36 ℃ and 42 ℃ referred to as ST and SY, respectively) were set up, and a 25 ℃ treatment was used as the control (CK). The experimental design included lighting at 2000 lx(14 h day/10 h night) in an artificial climate chamber. Three biological replicates were planted in the control and each treatment. Under continuous high-temperature stress in the SY treatment, Rhododendron molle displayed symptoms of yellowing and blackening of its leaves after 10 days, with signs of death observed after 14 days. Therefore, a high-temperature stress duration of 15 days was designed (Jiang et al. 2023). On the 15th day, healthy leaves without pests or diseases were collected from the target plants, rapidly frozen in liquid nitrogen, and stored at -80 °C. These samples were then sent to Nanjing Jisi Huiyuan Biotechnology Co., Ltd. for further analysis.

RNA Extraction and Sequencing Library Construction

The mRNA was extracted from the biological samples using the Poly (A) structure, and the RNA quality was determined using a Nanodrop ultra-micro UV-Vis spectrophotometer and Agilent 2100 Bioanalyzer. Eukaryotic mRNA was enriched using magnetic beads with Oligo (dT), and the first cDNA strand was synthesized by adding Fragmentation Buffer to randomly interrupt the mRNA, followed by the synthesis of the double-stranded cDNA strand, purification of the cDNA, and end-repair, etc. Finally, the cDNA library was enriched via PCR.

Assembly and Functional Annotation of the Transcriptome

High-throughput sequencing (e.g., Illumina HiSeq 2500) was used for experimental processing. The Unigene sequence were compared with the NR (Marchler et al. 2015), Swiss-Prot (Gasteiger et al. 2001), GO (Gene Ontology Consortium, 2004), COG (Tatusov et al. 2003), KOG (Jensen et al. 2007) and KEGG (Kanehisa et al. 2012) databases using BLAST v2.3.0+ (http://blast.ncbi.nlm.nih.gov/Blast.cgi) software, and the Unigene amino acid sequence were predicted and compared with the Pfam database (Mistry et al. 2021) using HMMER v3.4 (Finn et al. 2011) software to obtain the Unigene annotation information.

Functional Analysis of the DEGs

The sequenced reads were compared with the Unigene library using Bowtie2 v2.4.4 (Langmead and Salzberg 2012), and the expression level was estimated based on the comparison results in combination with RSEM (Li and Dewey 2011). The expressing abundance of the corresponding Unigene is expressed using FPKM values.

RESULTS AND DISCUSSION

Sequencing Data Assembly Results

Each sequencing sample was repeated once, and after sequencing quality control, the accuracy of each sample base was greater than 96% (corresponding to Q20 ), and the GC content was basically below 50%. Thus the sequencing data could be used for assembly (Moreton et al. 2016). The data output and assembly statistics for each sample are shown in Table 1. For the de novo assembly of the transcriptome, Trinity assembly software was adopted. The assembly resulted in a total of 124,901 Unigenes, with a total length of 530,057,171 nt. and an average length of 1538.21 nt, N50. All the Unigenes were arranged in order from longest to shortest, and their lengths were added. When a fragment’s length is 50% of the total number of fragments, the length and number of the fragment, i.e., the length and number of N50, indicate the quality of the assembly: the lower these values are for the Unigenes, the better the quality.

Table 1. Sequencing Sata Quality

Transcriptome Unigene Cluster Results

Since there were no reference genomic data for Rhododendron molle, Clean reads were used in Trinity software for hybrid splicing. A total of 344,593 transcripts and 124,901 Unigenes were obtained. The average length of the transcripts was 1538.21, and the N50 values of the transcripts and Unigenes were 2466 and 1725, respectively (Table 2), indicating high assembly integrity. This data met the requirements for the next step of gene function annotation and classification.

Table 2. Data Assembly Statistics

The assembled transcriptome sequences were compared with COG, GO, KEGG, KOG, Swiss-Prot, Pfam, NR, and other databases to obtain the annotation information of each database. The numbers of Unigenes annotated in the above databases were 10,973, 32,259, 51,672, 29,049, 23,701, 35,661, 62,917, respectively. The gene function annotation statistics are shown in Table 3.

Table 3. Annotation Statistics of Single Gene Clusters

Analysis of Differentially Expressed Genes Under Different Heat Stresses

Based on the RPKM values, the transcript gene expression differences under the different treatments were assessed, with screening criteria of FDR≤0.05 and FC≥2. In Fig. 1, the horizontal coordinate indicates the logarithmic value of the difference multiplier between the two treatment groups, and the vertical coordinate indicates the negative log10 value of the FDR of the difference between the two subgroups. In these volcano plots of differentially expressed genes, red (up-regulated expression) and green (down-regulated expression) dots indicate differential gene expression, while gray dots indicate no difference. By comparing the differential expression of genes between high-temperature stress and the control, the volcano plot analysis (Fig. 1) showed a total of 12,089 differentially expressed genes, of which 6004 were up-regulated and 6,085 were down-regulated, in the control compared with the 36 °C high-temperature stress. Additionally there were 12,032 differentially expressed genes in the control compared to the 42 °C high-temperature stress, of which 5,431 were up-regulated and 6,601 were down-regulated. There were 4594 common differentially expressed genes between the two groups of differential genes (Fig. 2), of which 1224 were co-up-regulated and 3007 were co-down-regulated, in addition to 363 common reverse differentially expressed genes between the two. Some of the common differentially expressed genes responded to both 36 °C and 42 °C high-temperature stress. The number of down-regulated genes in Rhododendron molle was higher than the number of up-regulated genes with increasing temperature under the high-temperature stress treatments, indicating that high-temperature stress significantly suppressed some of the physiological activities of Rhododendron molle.

Fig. 1. Volcano plot of the number of differentially expressed genes. Each point in the differential expression volcano plot represents a gene, and the horizontal coordinate indicates the logarithmic value of the expression ploidy of a gene in the two samples; the vertical coordinate indicates the negative logarithmic value of FDR. The larger the absolute value of the horizontal coordinate, the greater the difference in expression ploidy between the two samples; the larger the value of the vertical coordinate, the more significant the differential expression and the more reliable the differentially expressed genes obtained from the screening. The green dots represent down-regulated differentially expressed genes, the red dots represent up-regulated differentially expressed genes, and the gray dots represent non-differentially expressed genes.

Fig. 2. Venn diagram of the number of differentially expressed genes

GO Classification and Enrichment Analysis of Significantly Different Genes

The obtained target genes were subjected to GO enrichment analysis. The biological functions of the differentially expressed genes were mainly analyzed in terms of their molecular functions, cellular components, and biological processes. The biological processes enriched on metabolic processes, cellular processes, stimulus responses, and biological regulation; the cellular component part enriched on cells, cellular components, organelles, organelle components, membranes, membrane components, and protein complexes; and the molecular function part enriched on binders and catalytic activity. The GO enrichment results for Rhododendron molle under different high-temperature stresses are displayed in Fig. 3. The significantly enriched GO term categories did not change, but the numbers of enriched genes all increased, indicating that high-temperature stress led to increased plant cell destruction.

Fig. 3. Differential gene GO functions. Cellular component 1: Extracellular region; 2: Cell; 3: Nucleoid; 4: Membrane; 5: Cell junction; 6: Membrane-enclosed lumen; 7: Protein-containing complex; 8: Organelle; 9: Other organism; 10: Other organism part; 11: Extracellular region part; 12: Organelle part; 13: Membrane part; 14: Cell part; 15: Supramolecular complex . molecular function 1: Catalytic activity; 2: Structural molecule activity; 3: Transporter activity; 4: Binding; 5: Antioxidants activity; 6: Nutrient reservoir activity; 7: molecular transducer activity; 8: molecular function activity; 9: molecular carrier activity; 10: Transcription regulator activity. Biological process: 1: Reproduction; 2: Cell killing; 3: Immune system process; 4: Metabolic process; 5: Cell proliferation; 6: Cellular process; 7: Carbon utilization; 8: Nitrogen utilization; 9: Reproductive process; 10: Biological adhesion; 11: Signaling; 12: Multicellular organismal process; 13: Developmental process；14: Growth；15: Locomotion；16: Pigmentation；17: Rhythmic process; 18: Response to stimulus; 19: Localization; 20: multi-organism process; 21: Biological regulation; 22: Cellular component organization; 23: Detoxification.

KEGG Enrichment Analysis of Significantly Different Genes

The KEGG database was used to analyze the enrichment pathways of the differentially expressed genes detected in Rhododendron molle under different high-temperature stress treatments. As a result 4594 differentially expressed genes (DEGs) were found to be enriched in 125 pathways. The KEGG metabolic pathways for the control in comparison to 36 °C high-temperature stress mainly involved sulfur metabolism, astragalus, diarylheptane and gingerol biosynthesis, starch and sucrose metabolism, plant-pathogen interactions, phytohormone signaling, photosynthesis tentacle protein, phenylpropanoid biosynthesis, and other pathways. Among the top 20 metabolic pathways with the most reliable enrichment significance, as shown in Table 4, the plant-pathogen interaction pathway was annotated to a maximum of 88 genes; followed by spliceosomes (annotated to 60 differentially expressed genes), plant hormone signaling (annotated to 59 differentially expressed genes), ribosomes (annotated to 58 differentially expressed genes), and other pathways. The KEGG metabolic pathways for the control in comparison to 42 °C high-temperature stress mainly involved in plant-pathogen interactions, ribosomes, spliceosomes, endoplasmic reticulum protein processing, and other pathways. Among the top 20 metabolic pathways with the most reliable enrichment significance, the plant-pathogen interaction pathway was annotated to the most genes (100), followed by the ribosomal (annotated to 83 differentially expressed genes), ribosomal (annotated to 73 differentially expressed genes), and endoplasmic reticulum protein processing (annotated to 73 differentially expressed genes) pathways.

Table 4. KEGG Pathway Enrichment Information for the Top 10 Differentially Expressed Genes

Notes: ST: 36 °C high-temperature stress; SY: 42 °C high-temperature stress

Fig. 4. KEGG functional enrichment analysis of differentially genes

Functional Identification of Differentially Expressed Genes

From the enrichment results for differentially expressed genes in the KEGG signaling pathways, the related genes that were up-regulated or down-regulated in all metabolic pathways of lambda R. molle under high-temperature stress were identified. Two significantly up-regulated genes were screened from among the most enriched genes in the plant-pathogen interactions pathway under the ST high-temperature treatment(36 °C): oleuropein receptor protein-associated receptor kinase (BAK1) and disease process-associated gene transcriptional activator (PTI6).Four genes were significantly down-regulated, including mitogen-activated protein kinase kinase 1a (MKK1a), mitogen-activated protein kinase kinase kinase 5 (MKK5), and retroviral-associated Pol polyprotein (RE1) (Tables 5 through 7). Screening under the spliceosome pathway showed significant upregulation of the sulfate glycoprotein (RH14) and RNA unwinding helicase (SMD1A) and significant downregulation of a small nuclear ribonucleoprotein (STA1). Under the phytohormone signal transduction pathway, indole-3-acetic acid amide synthase (GH3.10), two-component response regulator (RR22), and DELLA protein (GAI1) were the three up-regulated genes and gene erythromycin receptor (GID1B) and oleuropein receptor protein-associated receptor kinase (BAK1) were the two down-regulated genes. Oleuropein receptor protein-related receptor kinase (BAK1), WRKY transcription factor (WRKY22), and transcriptional activator of disease process-related genes (PTI6) were significantly up-regulated in the phytopathogen interactions pathway under the SY high-temperature treatment (42 °C), while heat shock protein HSP83A, heat kinin (HSP90-1), and mitogen-activated protein kinase kinase 5 (MKK5) were significantly downregulated. The ribosomal proteins RPS3C and RPL24 showed up-regulation while RPS20B RPL10A and showed down-regulation under the ribosomal pathway. Protein homologs (MAGO1) showed up-regulation and serine-arginine-rich splicing factors (SC35, SR34, SR34B, SCL33) all showed down-regulation under the spliceosomal pathway (Tables 8-10). Genes in the metabolic pathway exhibited different levels of coordination involving multiple pathways such as photosynthesis, respiration, redox, phytohormone signaling, and gluconeogenesis, which work together to combat high-temperature stress.

Table 5. Information on Relevant DEGs Under the Interactions Pathway in Plants

Table 6. Information on Relevant DEGs Under the Spliceosome Pathway

Table 7. Information on Relevant DEGs Under the Phytohormone Signaling Pathway

Table 8. Information on Relevant DEGs Under the Plant-Pathogens Interactions Pathway

Table 9. Information on Relevant DEGs Under the Ribosomal Pathway

Table 10. Information on Relevant DEGs Under the Ribosomal Pathway

Discussion

High-temperature stress alters plant morphology, physiological and biochemical processes, and gene expression, impeding normal growth and development (Wang et al. 2023). This not only alters physiological metabolic processes such as photosynthesis, osmoregulation, antioxidant metabolism, phytohormones, and sugar metabolism but also affects a variety of signaling factors, transcription factors, and functional genes (Hasanuzzaman et al. 2013). The response of plants to high-temperature stress is a complex biological process (Kotak et al. 2007). In a previous study, the heat tolerance of five species of rhododendron was studied, and the catalase (CAT) activities, ascorbate peroxidase (APX), proline (Pro) content, hydrogen peroxide (H₂O₂), and malondialdehyde (MDA) in the leaves of R. molle increased with increasing stress temperature. The superoxide dismutase (SOD) activity increased slightly under mild stress at 30 °C but decreased under severe stress at 38 °C (Zhang et al. 2011). In a study on the effect of ERM fungi on the heat resistance of Rhododendron molle, in the control group without any fungi added, the electrical conductivity and MDA of R. molle increased with the temperature, while proline (Pro), POD, and SOD decreased with increasing stress (Luo 2016).

Rhododendrons grow best in cool environments and do not tolerate high temperatures. Thus, the high temperatures of China’s summer climate make it difficult for transplanted Rhododendron molle to survive. Under the same watering conditions, rhododendrons with high heat resistance were found to have a higher water utilization rate, CO₂ utilization rate, and light energy utilization rate; stronger photosynthetic capacity and antioxidant enzyme activity; and higher Pro and chlorophyll contents than heat-resistant rhododendrons in the same period (Li 2019). In our next steps, we will continue to study the differences in the expression levels of relevant genes in rhododendrons under water stress conditions, in combination with the gene changes that occur under high-temperature stress, to comprehensively analyze the key genes regulating stress resistance in rhododendrons in hot summers, so as to provide references for later rhododendron research.

In this study, high-temperature stress affected the physiological metabolism and growth of Rhododendron molle, and the transcriptome sequences of Rhododendron molle leaves under high-temperature stress and a control treatment were compared to reveal the response mechanism of Rhododendron molle to high-temperature stress at the molecular level. A study carried out in Hunan Provincial Botanical Garden showed that with an increase in the temperature and duration of high-temperature stress, the physiological and biochemical responses of Rhododendron molle leaves were enhanced, including a gradual decrease in SOD, MDA, and POD activities (Ying et al. 2018). To elucidate the biological mechanisms that occur in Rhododendron in response to high-temperature stress, the present study was conducted to analyze the transcriptome of Rhododendron molle leaves treated under high-temperature stress by means of Illumina high-throughput sequencing technology. The total number of genes in the sample was 124901, with a total length of 530057171 nt and an average length of 1538.21 nt. A total of 344593 transcripts and 124901 Unigenes were obtained after mixed splicing. A GO functional annotation analysis showed that the DEGs under high-temperature stress were mainly concentrated in metabolic processes, cellular processes, cells, cellular components, binders, and catalytic activity, and the most enriched DEGs were found in metabolic processes, This indicates that a large number of transcription factors were activated in Rhododendron molle to regulated its response to high-temperature stress by inducing the expression of high-temperature stress genes; in particular, cell destruction was severely damaged (Yu et al. 2011). An analysis of KEGG metabolic pathways revealed that the most enriched metabolic pathways for the DEGs were plant-pathogen interactions, spliceosomes, ribosomes, phytohormone signaling, and endoplasmic reticulum protein processing. This indicates that the metabolic pathway plays an important role in the response of Rhododendron molle to high-temperature stress. The enrichment pathways screened for transcription factors such as WRKY were previously found to be involved in the synthesis of plant terpenoids to withstand adversity (Singh et al. 2017; Alfieri et al. 2018; Zhan et al. 2022).

Brassinolide (BR), also known as brassinosteroid, is a naturally occurring compound with biological activity, characterized by a steroidal structure. It represents the sixth type of plant hormone, following gibberellins, auxins, cytokinins, abscisic acid, and ethylene (Sasse 1997). The discovery of BR was one of the most significant findings in the field of plant hormones, following that of gibberellins. This hormone is involved in regulating various processes in plants, including seed germination, nutritional growth, reproductive growth, and responses to external environmental stressors (Zheng et al. 2018; Rozhon et al. 2019). Cao and Zhao (2008) found that BR can enhance the activity of superoxide dismutase (SOD) and peroxidase (POD), increase the expression of their isoenzymes, eliminate free radicals generated within cells under high-temperature conditions, stabilize membrane structures, reduce electrolyte leakage, and increase the protein content in plants. This enables normal plant growth in high-temperature stress conditions. A study by Jin et al. (2015) confirmed that BR can alleviate the damage caused by high-temperature stress in banyan tree seedlings. Furthermore, a study conducted by Yang et al. (2019) suggested that brassinosteroid (BR) treatment can alleviate the inhibitory effects of high-temperature stress on the growth of Toona sinensis seedlings. In our research, we found that brassinosteroid (BR) receptor protein-associated receptor kinase (BAK1) was significantly upregulated during the heat-stress treatment (ST). This suggests that brassinosteroids play a crucial role in the response of Rhododendron molle to high-temperature stress.

A particular group of plant-coded proteins induced by different stress stimuli, named “pathogenesis-related proteins” (PRs), is assigned an important role in plant defense against pathogenic constraints and in general adaptation to stressful environments, Upon investigation, it was discovered that these genes play a pivotal role in responding to both biotic and abiotic stressors, underscoring their critical importance in the survival of plants (Edreva 2005). In previous studies, there was no mention of PRs being involved in plant responses to heat stress. In our research, we discovered that Pathogenesis-Related gene Transcriptional Activator (PTI6) was significantly upregulated during heat treatment (ST), indicating that PRs may play a crucial role in the response of sheep sorrel to high-temperature stress. This finding provides important theoretical support for the involvement of PRs in the response to heat stress. Therefore, PTI6 emerged as a novel and advantageous candidate gene in our study, laying a solid foundation for future in-depth investigations into Rhododendron molle‘s adaptation to high-temperature stress.

Key heat-stimulated transcription factors (Hsf) play an important role in the response to high temperature (Pirkkala et al. 2001; Mittal et al. 2009; Liu et al. 2015). Hsf genes, as heat shock transcription factors, often bind to heat-excited elements in the heat shock protein (HSP) promoter and play important roles in the regulation of plant development and various stress responses (Li et al. 2023). In our study, we observed a significant downregulation of heat shock proteins (HSP83A and HSP90-1) under the SY treatment. This contradicts the typical plant response to combat high-temperature stress, suggesting that this downregulation may be one of the reasons for temperature sensitivity in sheep sorrel.

CONCLUSIONS

1. A transcriptome analysis was conducted on the leaves of Rhododendron molle. Leaves from plants kept at 25 °C were used as the control group and were compared with leaves from plants kept at 36 °C or 42 °C.
2. According to the GO enrichment, KEGG enrichment, and expression trend analyses of the DEGs, four candidate genes expressed in response to HTS were identified. BAK1, PTI6, HSP83A, and HSP90-1. BAK1 and PTI6 exhibited significant upregulation under the ST treatment, while HSP83A and HSP90-1 showed significant downregulation under the SY treatment.
3. The results of this study will be helpful for further studies of the molecular mechanisms of Rhododendron molle in response to HTS. The results provide a theoretical basis for exploring these molecular mechanisms and selecting high-quality heat-tolerant varieties.

ACKNOWLEDGMENTS

The authors are grateful for the support of Hunan Province Forestry Science and Technology Innovation Program (XLKY202212, XLKY202313, XLKJ202203, XLKT202303).

REFERENCES CITED

Alfieri, M., Vaccaro, M. C., Cappetta, E., Ambrosone, A., De Tommasi, N., and Leone, A. (2018). “Coactivation of MEP-biosynthetic genes and accumulation of abietane diterpenes in Salvia sclarea by heterologous expression of WRKY and MYC2 transcription factors,” Scientific Reports 8(1), 11009. DOI: 10.1038/s41598-018-29389-4

Arisumi, K., Matsuo, E., and Sakata, Y. (1989). “Breeding for the heat-resistant rhododendrons: IV. The reliability of early selection of seedlings on to the heat resistance,” Memoirs of the Faculty of Agriculture Kagoshima University 25, 75-82.

Cao, Y. Y., and Zhao, H. (2008). “Protective roles of Brassinolide on rice seedlings under high temperature stress,” Chinese Journal of Rice Science 21(5), 525-529.

Cheng, X. A., Xie, J. J., Hu, M. Y., Zhang, Y. B., and Huang, J. F. (2011). “Induction of intracellular Ca²⁺ and pH changes in Sf9 insect cells by rhodojaponin-III, a natural botanic insecticide isolated from Rhododendron molle,” Molecules 16(4), 3179-3196. DOI: 10.3390/molecules16043179

Edreva, A. (2005). “Pathogenesis-related proteins: Research progress in the last 15 years,” Gen Appl Plant Physiol, 31(12), 105-24.

Finn, R. D., Clements, J., and Eddy, S. R. (2011). “HMMER web server: Interactive sequence similarity searching,” Nucleic Acids Research 39 (suppl_2), W29-W37. DOI: 10.1093/nar/gkr367

Gasteiger, E., Jung, E., and Bairoch, A. (2001). “SWISS-PROT: connecting biomolecular knowledge via a protein database,” Current issues in Molecular Biology 3(3), 47-55.

Gene Ontology Consortium. (2004). “The Gene Ontology (GO) database and informatics resource,” Nucleic Acids Research 32, 258-261. DOI: 10.1093/nar/gkh036

Geng, X. M., Yang, Q. Y., Yue, Y., and Ozaki, Y. (2019). “Effects of high temperature on photosynthesis, membrane lipid peroxidation and osmotic adjustment in four Rhododendron species,” Journal of the Faculty of Agriculture, Kyushu University 64(1), 33-38. DOI: 10.5109/2231631

Hasanuzzaman, M., Nahar, K., Alam, M. M., Roychowdhury, R., and Fujita, M. (2013). “Physiological, biochemical, and molecular mechanisms of heat stress tolerance in plants,” International Journal of Molecular Sciences 14(5), 9643-9684. DOI: 10.3390/ijms14059643

Jensen, L. J., Julien, P., Kuhn, M., von Mering, C., Muller, J., Doerks, T., and Bork, P. (2007). “eggNOG: Automated construction and annotation of orthologous groups of genes,” Nucleic Acids Research 36(suppl_1), D250-D254. DOI: 10.1093/nar/gkm796

Jiang, J. C., Liao, J. Y., Cao, S. J., Xu, K., and Wu, L. S. (2023). “Effects of high temperature stress on the growth and physiology of three Rhododendron species,” Northern Horticulture 9(18), 54-62. DOI: 10.11937/bfyy.2022.5177

Jin, S. H., Li, X. Q., Wang, G. G., and Zhu, X. T. (2015). “Brassinosteroids alleviate high-temperature injury in Ficus concinna seedlings via maintaining higher antioxidant defence and glyoxalase systems,” Aob Plants 7, 16-21. DOI: 10.1093/aobpla/plv009

Kanehisa, M., Goto, S., Sato, Y., Furumichi, M., and Tanabe, M. (2012). “KEGG for integration and interpretation of large-scale molecular data sets,” Nucleic Acids Research 40(D1), D109-D114. DOI: 10.1093/nar/gkr988

Kotak, S., Larkindale, J., Lee, U., von Koskull-Döring, P., Vierling, E., and Scharf, K. D. (2007). “Complexity of the heat stress response in plants,” Current Opinion in Plant Biology 10(3), 310-316. DOI: 10.1016/j.pbi.2007.04.011

Langmead, B., and Salzberg, S. L. (2012). “Fast gapped-read alignment with Bowtie 2,” Nature Methods 9(4), 357-359. DOI: 10.1038/nmeth.1923

Li, B., and Dewey, C. N. (2011). “RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome,” BMC Bioinformatics 12, 1-16. DOI: 10.1186/1471-2105-12-323

Li, H. (2019) Study on Evaluation of Heat Resistance Resources and Physiological Mechanism of Rhododendron, Master’s thesis, Yangzhou University. DOI:10.27441/d.cnki.gyzdu. 2019.000835

Li, S., Chen, H., Yu, H., Li, Y., and Wang, L. (2023). “Responses and adaptations of fruit trees to high temperatures,” Fruit Research 3(1), 15-19. DOI: 10.48130/FruRes-2023-0023

Li, Y., Zhu, Y. X., Zhang, Z. X., Liu, Y. L., Liu, Y. B., Qu, J., and Yu, S. S. (2018). “Diterpenoids from the fruits of Rhododendron molle, potent analgesics for acute pain,” Tetrahedron 74(7), 693-699. DOI: 10.1016/j.tet.2017.12.017

Liu, Z. L., Goh, S. H., and Ho, S. H. (2007). “Screening of Chinese medicinal herbs for bioactivity against Sitophilus zeamais Motschulsky and Tribolium castaneum (Herbst),” Journal of Stored Products Research 43(3), 290-296. DOI: 10.1016/j.jspr.2006.06.010

Liu, Z., Xin, M., Qin, J., Peng, H., Ni, Z., Yao, Y., and Sun, Q. (2015). “Temporal transcriptome profiling reveals expression partitioning of homeologous genes contributing to heat and drought acclimation in wheat (Triticum aestivum L.),” BMC Plant Biology 15, 1-20. DOI: 10.1186/s12870-015-0511-8

Lou, Q. (2016). Study on the Effect of ERM Fungi on Introduction and Heat-Resistance of Rhododendron molle, Master’s thesis, South China Agricultural University.

Marchler-Bauer, A., Derbyshire, M. K., Gonzales, N. R., Lu, S., Chitsaz, F., Geer, L. Y., and Bryant, S. H. (2015). “CDD: NCBI’s conserved domain database,” Nucleic Acids Research 43(D1), D222-D226. DOI: 10.1093/nar/gku1221

Meijón, M., Rodríguez, R., Canal, M. J., and Feito, I. (2009). “Improvement of compactness and floral quality in azalea by means of application of plant growth regulators,” Scientia Horticulturae 119(2), 169-176. DOI: 10.1016/j.scienta.2008.07.023

Mistry, J., Chuguransky, S., Williams, L., Qureshi, M., Salazar, G. A., Sonnhammer, E. L., and Bateman, A. (2021). “Pfam: The protein families database in 2021,” Nucleic Acids Research 49(D1), D412-D419. DOI: 10.1093/nar/gkaa913

Mittal, D., Chakrabarti, S., Sarkar, A., Singh, A., and Grover, A. (2009). “Heat shock factor gene family in rice: genomic organization and transcript expression profiling in response to high temperature, low temperature and oxidative stresses,” Plant Physiology and Biochemistry 47(9), 785-795. DOI: 10.1016/j.plaphy.2009.05.003

Moreton, J., Izquierdo, A., and Emes, R. D. (2016). “Assembly, assessment, and availability of de novo generated eukaryotic transcriptomes,” Frontiers in Genetics 6, 169488. DOI: 10.3389/fgene.2015.00361

Pirkkala, L., Nykänen, P., and Sistonen, L. E. A. (2001). “Roles of the heat shock transcription factors in regulation of the heat shock response and beyond,” The FASEB Journal 15(7), 1118-1131. DOI: 10.1096/fj00-0294rev

Rozhon, W., Akter, S., Fernandez, A., and Poppenberger, B. (2019). “Inhibitors of brassinosteroid biosynthesis and signal transduction,” Molecules 24(23), 4372. DOI: 10.3390/molecules24234372

Sasse, J. M. (1997). “Recent progress in brassinosteroid research,” Physiologia Plantarum 100(3), 696-701. DOI: 10.1111/j.1399-3054.1997.tb03076.x

Schlenker, W., and Auffhammer, M. (2018). “The cost of a warming climate,” Nature 26: 498-499. DOI: 10.1038/d41586-018-05198-7

Shin, F, C., and Yu‐Tong, Q. (1993). “Experiments on the application of botanical insecticides for the control of diamondback moth in South China,” Journal of Applied Entomology 116(1‐5), 479-486. DOI: 10.1111/j.1439-0418.1993.tb01224.x

Singh, A. K., Kumar, S. R., Dwivedi, V., Rai, A., Pal, S., Shasany, A. K., and Nagegowda, D. A. (2017). “A WRKY transcription factor from Withania somnifera regulates triterpenoid with anolide accumulation and biotic stress tolerance through modulation of phytosterol and defense pathways,” New Phytologist 215(3), 1115-1131. DOI: 10.1111/nph.14663

Tatusov, R. L., Fedorova, N. D., Jackson, J. D., Jacobs, A. R., Kiryutin, B., Koonin, E. V., and Natale, D. A. (2003). “The COG database: An updated version includes eukaryotes,” BMC Bioinformatics 4, 1-14. DOI: 10.1186/1471-2105-4-41

Ureshino, K., Miyajima, I., and Akabane, M. (1998). “Effectiveness of three-way crossing for the breeding of yellow-flowered evergreen azalea,” Euphytica 104, 113-118. DOI: 10.1023/A:1018661004467

Wang, S., Lin, S., Zhu, C., Yang, Y., Li, S., Zhang, J., and Shi, J. (2010). “Highly acylated diterpenoids with a new 3, 4-secograyanane skeleton from the flower buds of Rhododendron molle,” Organic Letters 12(7), 1560-1563. DOI: 10.1021/ol1002797

Wang, X., Altaf, M. A., Hao, Y., Wang, Z., and Zhu, G. (2023). “Effect of heat stress on root architecture, photosynthesis, and antioxidant profile of water spinach (Ipomoea aquatica Forsk) seedlings,” Horticulturae 9(8), 923. DOI: 10.3390/horticulturae9080923

Yan, D., Saito, K., Ohmi, Y., Fujie, N., and Ohtsuka, K. (2004). “Paeoniflorin, a novel heat shock protein–inducing compound,” Cell Stress and Chaperones 9(4), 378. DOI: 10.1379/CSC-51R.1

Yang, C. C., Li, Z., Li, J. A., Liu, S. T., Zhang, F. H., and Li, Y. L. (2019). “Effects of spraying Brassinolide on growth and photosynthesis of seedlings of Toona sinensis in different temperatures,” Journal of Central South University of Forestry and Technology 39(1), 25-31.

Ying, Z., Yu, W., Hu, X., Shi, Y., Liu, Y., Zhong, Y., and Yu, X. (2018). “Physiological and transcriptomic analysis revealed the involvement of crucial factors in heat stress response of Rhododendron hainanense,” Gene 660, 109-119. DOI: 10.1016/j.gene.2018.03.082

Yu, Q., Jin, X., Liu, C., and Wen, Y. (2022). “An integrated analysis of transcriptome and miRNA sequencing provides insights into the dynamic regulations during flower morphogenesis in petunia,” Horticulturae 8(4), 284. DOI: 10.3390/horticulturae8040284

Yuan, C., Zhu, G., Yang, S., Xu, G., Li, Y., Gong, H., and Wu, C. (2019). “Soil warming increases soil temperature sensitivity in subtropical forests of SW China,” PeerJ 7, e7721. DOI: 10.7717/peerj.7721

Zhan, X., Chen, Z., Chen, R., and Shen, C. (2022). “Environmental and genetic factors involved in plant protection-associated secondary metabolite biosynthesis pathways,” Frontiers in Plant Science 13, 877304. DOI: 10.3389/fpls.2022.877304

Zhang, L. H., Sun , B. T., Zhou, G .,Wang, S. S., Li, X. H., and Dan, W. (2011). “Physiological changes and heat tolerance comparison of five Rhododendron species under high temperature stress,” Guihaia 31(5), 651-658. DOl:10.3969/i.issn.1000-3142.2011.05.016

Zheng, L., Ma, J., Zhang, L., Gao, C., Zhang, D., Zhao, C., and Han, M. (2018). “Revealing critical mechanisms of BR-mediated apple nursery tree growth using iTRAQ-based proteomic analysis,” Journal of Proteomics 173, 139-154. DOI: 10.1016/j.jprot.2017.12.007

Zhi, X., Xiao, L., Liang, S. et al. (2013). “Chemical constituents of Rhododendron molle,” Chemistry of Natural Compounds 49, 454-456. DOI: 10.1007/s10600-013-0637-6

Zhou, G. L., and Zhu, P. (2020). “De novo transcriptome sequencing of Rhododendron molle and identification of genes involved in the biosynthesis of secondary metabolites,” BMC Plant Biology 20, 1-19. DOI: 10.1186/s12870-020-02586-y

Zhu, Y. L., Deng, L., Tang, Y., Fan, X. Z., Han, Y., Pan, M., and Liao, H. B. (2023). “New polychlorinated bibenzyls from Rhododendron minutiflorum,” Natural Products and Bioprospecting 13(1), 2-8. DOI: 10.1007/s13659-022-00364-x

Article submitted: March 25, 2024; Peer review completed: June 29, 2024; Revisions accepted: August 23, 2024; Published: September 12, 2024.

DOI: 10.15376/biores.19.4.8238-8256