Yao Chen1
Yuanhang Xiang1
Zhengrong Hu2Yang Gao1Youxin Zhang1
Minghui Chen1A. B. M. Khaldun3
Xuebing Yan1*
Jibiao Fan1*- 1College of Animal Science and Technology, Yangzhou University, Yangzhou, China
- 2Hunan Tobacco Research Institute, Changsha, China
- 3Bangladesh Agricultural Research Institute, Dhaka, Bangladesh
Soil salinization is a major problem all over the world. The accumulation of salt in soil reduces the root water uptake and directly affects plant growth and metabolic activities. Brassinosteroid is a plant hormone that plays an important role in regulation of plant growth and physiological process, including promotion of cell expansion and elongation, signal transduction and stress response. Exogenous 24-epibrassinolide (EBL) has been proved to alleviate various environmental stress in plants. However, the role that EBL plays in salt stress response is still unknown in tall fescue (Festuca arundinacea). In this study, the physiology and molecular mechanisms regulated by exogenous EBL of salt stress response in tall fescue was investigated. Tall fescue plants were divided into four groups, including control (CK), NaCl solution (SALT), 24-epibrassinolide (EBL), NaCl solution + 24-epibrassinolide (SE). During the growth period of tall fescue, we found that electrolyte leakage (EL) and malondialdehyde (MDA) were decreased, chlorophyll (Chl) content and antioxidant enzyme activity were increased in leaves of tall fescue in SE group compared with SALT group, indicating that EBL improved the salt tolerance in grasses. Transcriptomic profiling analysis showed that after 12 h of treatments, 10,265, 13,830 and 10,537 differential genes were expressed in EBL, SALT, and SE groups compared with control, respectively. These differentially expressed genes (DEGs) mainly focused on binding, catalytic activity, cellular process, metabolic process, cellular anatomical entity. Moreover, most of the differential genes were expressed in the plant hormone signal transduction pathway. These results helped us to better understand the mechanism of exogenous 24-epibrassinolide to improve the salt tolerance of tall fescue.
Introduction
Soil salinization is a major problem worldwide. Around 1 billion hectares of land throughout the world are affected by salinity, and these soils are distributed in more than 100 countries (Food and Agriculture Organization of the United Nations, 2015). Soil salinity is induced by natural causes and human behaviors such as improper irrigation (Munns and Tester, 2008). And the accumulation of salt in the soil can seriously affect growth and yield of the plant (Zelm et al., 2020). Soil salinization makes it difficult for roots to absorb water and further disrupt the water potential and ion balance within plant cells. The accumulation of salt in the aerial parts will reduce the photosynthetic capacity of plants and destroy the antioxidant enzyme system, which resulting in a large increase in oxygen free radicals and continuous reduction in growth. And with the extension of time, the damage caused by salt stress is more and more severe, which eventually leads to the stagnation of plant growth or even death (Zhu, 2001; Munns and Tester, 2008; Zelm et al., 2020). The extent and severity of salt-affected lands will worsen due to issues such as global warming (Munns and Gilliham, 2015).
To defense salt stress, plants have developed various strategies to optimize growth conditions under salinity. In addition to regulating plant growth and development, phytohormones also participate in various environmental stresses to improve stress resistance (Yu et al., 2020). Brassinosteroids (BRs) are a new class of plant hormones that are essential for plant growth and development (Kim et al., 2009). Furthermore, it exhibits excellent growth-promoting activity at very low concentrations (Avalbaev et al., 2021). 24-Epibrassinolide (EBL), an active by-product of brassinolide biosynthesis, can stimulate plant metabolic processes such as protein and nucleic acid biosynthesis (Bajguz, 2000). It can significantly improve the ability of plants to resist abiotic stresses, including salt, heavy metal, drought, and low-temperature (Arora et al., 2012; Wu et al., 2017; Anwar et al., 2018; Tanveer et al., 2019).
Tall fescue (Festuca arundinacea Schreb) is a broad-leaved outbreeding allohexaploid grass (Hand et al., 2010). As one of the most important cool-season grass species in the world, tall fescue is widely used as turfgrass and forage (Wang et al., 2004). On account of its strong adaptability, tall fescue can grow in temperate and subtropical regions (Cao et al., 2009). Nonetheless, regional salinization limits the growth of tall fescue, which in turn affects turf persistence and forage yield (Wang et al., 2017).
Salinity can cause the increase of reactive oxygen species (ROS) and the content of malondialdehyde (MDA) (Miller et al., 2010; Hu et al., 2012). Plants have a complex antioxidant defense system to alleviate this oxidative damage, including catalase (CAT), peroxidase (POD), superoxide dismutase (SOD) (Lou et al., 2017). In addition, transcriptome analysis is also an available method to study salt stress responses and has been applied to many plants such as Arabidopsis, rice, alfalfa, bermudagrass (Shen et al., 2014; Shi et al., 2015; Wang et al., 2016; Bhattarai et al., 2021). However, there are few studies on the transcriptomic data on the effects of exogenous 24-epibrassinolide in tall fescue under salt stress. Therefore, this study elucidates the mechanism of action of 24-epibrassinolide in alleviating salt stress from both the physiological and transcriptomic levels.
Materials and methods
Plant material and growth conditions
Tall fescue (F. arundinacea) cultivar Houndog 6 was used for this study and quartz sand was used as growth material. All materials were grown in a climatic chamber (ROX-1000CH, Changzhou, China) at 25°C temperature, 14/10 h (light/dark), 60% relative humidity.
Treatments
All plants were divided into four groups: control (CK), 1 μmol/L 24-epibrassinolide (EBL), 300 mmol/L NaCl solution (SALT), 300 mmol/L NaCl solution + 1 μmol/L 24-epibrassinolide (SE). Three repetitions were set up for each group, and 100 seeds were evenly sown in each pot. The plants were watered daily until all seeds germinated, then irrigated with 50 ml 1/2 Hoagland nutrient solution per pot. After the plants had grown for 20 days, NaCl and 24-epibrassinolide treatments were conducted. The salt treatment was increased by daily increments of 100 mmol/L NaCl until a final salinity level (300 mmol/L) was reached. Plants were watered with NaCl solution until the liquid flowed from the drain holes at the bottom of the pot. The leaves of plants in the EBL and SE treatments were sprayed with 24-epibrassinolide evenly every day.
Measurement of chlorophyll content
0.1 g of fresh leaves were immersed in a 15 ml centrifuge tube containing 10 ml of dimethylsulfoxide and were placed in the dark for 48 h. Then, the absorbance of the extraction was measured at 645 and 663 nm with an ultraviolet spectrophotometer (Hiscox and Israelstam, 1979).
Chlorophyll content was calculated with the following formula:
OD645 and OD663 indicate the absorbance of the extraction at 645 and 663 nm, respectively.
Measurement of electrolyte leakage
The determination of electrolyte leakage (EL) referred to the method of Feng et al. (2021) with slight modifications. In detail, 0.1 g of fresh leaves were chopped and placed in a 50 ml centrifuge tube containing 25 ml of deionized water. The tubes were shaken at room temperature for 24 h at 200 rpm. The initial conductivity (Ci) was measured with a conductivity meter (DDS-11A conductivity Meter, Hangzhou, China). Then, the leaf tissue in the tube was autoclaved at 121°C for 20 min to release the electrolytes completely. The conductivity (Cmax) was measured again after the solution had been cooled to room temperature. The EL was calculated by the formula:
Crude enzyme extraction
0.2 g of fresh leaves were quickly ground into powder with liquid nitrogen. 4 mL of pre-chilled sodium phosphate buffer (pH 7.8) was added to the powder and the homogenate was transferred to a 15 mL centrifuge tube. After centrifugation at 4°C at 12,000 rpm for 15 min, the supernatant was extracted as crude enzyme solution.
Measurement of malonaldehyde content
Malonaldehyde (MDA) content was determined by thiobarbituric acid (TBA) method (Fan et al., 2015). Add 1 ml of crude enzyme solution to 2 ml of MDA reaction buffer containing 0.5% (v/v) TBA and 20% (v/v) trichloroacetic acid. The mixed solution was water-bathed at 95°C for 30 min and then rapidly cooled to room temperature in ice. Then the tubes were centrifuged for 10 min at 25°C with 12,000 rpm. The solution was measured at 532 and 600 nm with an ultraviolet spectrophotometer. MDA content was calculated with following formula:
Where L is the volume of the extraction, l is the thickness of the cuvette, ε is the molar absorption coefficient of 155 mM–1 cm–1, and FW is the fresh weight of the leaves.
Quantification of antioxidant enzymes
To determine CAT activity, 0.1 ml of crude enzyme solution was added to a mixed solution containing 2 ml phosphate buffer (100 mM, pH 7.0), 0.5 ml 100 mM hydrogen peroxide (H2O2). Absorbance at 240 nm was measured three times with an ultraviolet spectrophotometer and the values were recorded per minute. A reduction of one unit of the absorbance per minute was defined as one unit CAT activity.
To determine POD activity, 20 μL of crude enzyme solution was added to a mixed solution containing 2.5 mL phosphate buffer (100 mM, pH 6.0), 1.4 μL guaiacol, 0.95 μL 30% hydrogen peroxide (30% H2O2). Absorbance at 470 nm was measured three times with an ultraviolet spectrophotometer and the values were recorded per minute. Increment of one unit of the absorbance per minute was defined as one unit POD activity.
To determine SOD activity, 0.1 mL of crude enzyme solution was added to a mixed solution containing 1.5 ml phosphate buffer (50 mM, pH 7.8), 0.3 ml 130 mM methionine (Met), 0.3 ml 750 μM nitroblue tetrazolium (NBT), 0.3 ml 100 μM disodium ethylenediamine tetraacetic acid (EDTA-Na2), and 0.3 mL 20 μM riboflavin. Take the mixed solution without crude enzyme solution as control. Another new tube was regarded as blank, in which phosphate buffer (50 mM, pH 7.8) was used instead of the crude enzyme solution. The control and the experimental groups were placed under the condition of 4,000 Lux for 30 min, and the blank was placed in the dark for 30 min. Absorbance at 560 nm was measured with an ultraviolet spectrophotometer and the blank was used as zero adjustment. Taking NBT photoreduction by 50% as one unit of SOD activity.
RNA extraction and library construction
Total RNA was extracted from plant tissues using the ethanol precipitation protocol and CTAB-PBIOZOL reagent according to the manufacturer’s instructions. The concentration and purity of RNA were tested by the Agilent 2100 Bioanalyzer (Thermo Fisher Scientific, MA, USA) and NanoDrop. Use DNase I to digest the DNA fragments presenting in total RNA samples. The magnetic beads were purified to recover the reaction products. RNase H (Illumina, CA, USA) was used to remove the rRNA. Purified mRNA from previous steps was fragmented into small pieces with fragment buffer at appropriate temperature. Then, First-strand cDNA was generated in First Strand Reaction System by PCR, and the second-strand cDNA was generated as well. The reaction product was purified by magnetic beads. A-Tailing Mix and RNA Index Adapters were added by incubating to carry out end repair. The cDNA fragments with adapters were amplified by PCR. The PCR products were heat denatured into single-stranded DNA, and a single-stranded circular DNA library was obtained by circularizing the single-stranded DNA with a bridge primer.
Sequencing data filtering and de novo assembly
Raw reads were filtered via SOAPnuke (v1.4.0) (Chen et al., 2018) by removing reads containing adaptors, poly-N or low quality, and then clean reads were obtained. Trinity (v2.0.6) (Grabherr et al., 2011) was used to assemble the clean reads and Tgicl (v2.0.6) (Pertea et al., 2003) to perform clustering and eliminate redundant data in the assembled transcripts to obtain unique genes. Unigenes were divided into two parts. One part were clusters, which were the results of further de-redundancy. And there were several unigenes with high similarity (greater than 70%) in the same cluster, starting with CL, followed by the number of the gene family. The rest were singletons, starting with Unigene, referring to a single unigene without clustering.
Coding DNA sequence forecast and unigene annotation
The candidate coding regions in Unigene were identified using TransDecoder (v3.0.1) (Kim et al., 2015), the longest open reading frame was extracted. Based on sequence similarity, Diamond Blastp was used to align to SwissProt database, and then Hmmscan (v3.0) was used to screen the Blast results in the Pfam database. Trans Decoder. Predict was used to predict coding regions finally. The assembled unigenes were annotated with seven functional databases (KEGG, GO, NR, NT, SwissProt, Pfam, and KOG).
Gene expression analysis
Gene expression levels in each sample were calculated by RSEM (v1.2.8) (Li and Dewey, 2011) and expression analysis of differential genes was performed by the DESeq2 with Qvalue (adjusted p-value) ≤ 0.05 (Anders and Huber, 2010; Love et al., 2014).
Quantitative real-time PCR
cDNA was synthesized using the HiScript® III RT SuperMix for qPCR (+ gDNA wiper) (Vazyme, Nanjing, China). Quantitative real-time PCR (qRT-pCR) reactions were performed on QuantStudio™ 3 Real-Time PCR System (Thermo Fisher Scientific, MA, USA) using the AceQ® Universal SYBR qPCR Master Mix (Vazyme, Nanjing, China). The thermal cycle program was as following: 95°C for 5 min, 40 cycles of 95°C for 10 s and 60°C for 30 s; 95°C for 15 s, 60°C for 1 min, 95°C for 15 s. The information of primers for the genes for qRT-PCR are listed in Supplementary Table 1. The expression data were calculated by the 2–ΔΔCt method (Pfaffl, 2001).
Results
Effects of exogenous 24-epibrassinolide on growth and physiological responses under salt stress
In order to investigate whether exogenous 24-epibrassinolide could improve the growth of tall fescue under salt stress, the plant height was measured. The results showed that the plant height of the four groups increased at 10 days after different treatments. Among them, the two groups (SALT and SE) treated with NaCl solution were significantly lower than the CK and the EBL. After application of 24-epibrassinolide, the plant heights of EBL and SE groups were increased by 4.5% and 5.2% compared with the CK and SALT, respectively (Figure 1A and Supplementary Figure 1). These results indicated that 24-epibrassinolide can enhance the growth of tall fescue under salt stress remarkably.
Figure 1. Effects of exogenous 24-epibrassinolide on growth and physiological characteristics of tall fescue under salt stress. (A) Plant height of tall fescue; (B) Chlorophyll content of tall fescue; (C) MDA content; (D) Electrolyte Leakage; (E) SOD activity; (F) POD activity; (G) CAT activity. Bars with different letters indicate significant difference at P ≤ 0.05 (Duncan test).
To investigate the effect of 24-epibrassinolide on the chlorophyll of the leaves under salt stress, the chlorophyll content was determined. After 10 days of treatment, the chlorophyll content in the leaves after SE treatment was the highest, followed by SALT and EBL, and the content of CK was the lowest. With the application of 24-epibrassinolide, the chlorophyll content in the leaves after EBL and SE treatments increased by 9.9% and 6.1%, respectively compared with CK and SALT (Figure 1B). These results implied that exogenous 24-epibrassinolide could maintain the chlorophyll stability of tall fescue under salt stress.
To investigate whether the exogenous 24-epibrassinolide played a positive role in maintaining cell membrane stability of tall fescue under salt stress, MDA content and EL alterations were determined. The contents of MDA and EL in tall fescue after 10 days of salt treatment were much higher than those of the CK. However, there were no significant differences observed in MDA and EL contents between CK and EBL groups. After exogenous application of 24-epibrassinolide, MDA and EL contents of SE treatment decreased by 1.09% and 27.8% compared with the SALT, respectively (Figures 1C,D). These results suggested that exogenous 24-epibrassinolide could maintain cell membrane stability.
To investigate the effects of 24-epibrassinolide on antioxidant enzymes, activities of SOD, POD and CAT were determined (Figures 1E–G). After 10 days of salt treatment, the activities of SOD and POD were enhanced, but CAT had no significant difference compared with control. After the application of 24-epibrassinolide, activities of SOD, POD and CAT increased by 10.7, 5.5, and 3.4% in the leaves of EBL treatment compared with the CK. The SOD, POD, and CAT activities in plants of SE treatment were 34.6, 2.5, and 7.5% higher than those of SALT, respectively. These results demonstrated that 24-epibrassinolide plays an important role in enhancing antioxidative enzyme activities of tall fescue in response to salt stress.
Transcriptome sequencing analysis and unigene annotation
The transcriptome sequencing obtained 129,866 unigenes after de novo assembly and de-redundancy. The total length, average length, N50 and GC content were 135,692,978 bp, 1,044 bp, 1,466 bp, and 48.89%, respectively (Supplementary Table 2). Among them, 19,111 unigenes with a length of 200–300 bp accounted for the largest number, and the number of unigenes with a length of 2,900-3,000 bp accounted for the least, with a total of 541. It is worth noting that there were 3,293 genes with length greater than or equal to 3,000 bp (Figure 2A). It can be seen from the expression stacking chart that the Fragments Per Kilobase of exon model per Million mapped fragments (FPKM) values between 1 and 10 in each group accounted for the most, and the number greater than 10 was the least (Figure 2B).
Figure 2. Overview of the tall fescue transcriptome dataset. (A) Unigene length distribution; (B) Number of genes in different FPKM intervals for each sample; (C) Venn diagram of KEGG, NT, NR, KOG, and Swissprot annotation results of the bermudagrass reference transcriptome; (D) Distribution and quantity statistics of KEGG functions.
Unigenes were compared to seven functional databases for annotation, there were 85,000 in NR (65.45%), 76,023 in NT (58.54%), 58,513 in SwissProt (45.06%), 58,312 in KOG (44.90%), 61,897 in KEGG (47.66%), 61,289 in GO (47.19%) and 58,434 in Pfam (45.00%) (Figure 2C and Supplementary Table 3). We used the KEGG database to analyze pathway enrichment for all unigenes. As a result, 61,897 unigenes were annotated into five categories: cellular processes, environmental information processing, genetic information processing, metabolism, and organismal system. In this classification, genes were further divided into 19 functional pathways. Metabolic categories had the most types of enriched pathways, among which, carbohydrate metabolism accounts for the largest proportion except the global and overview maps. Genetic information processing was the second functional category with four subcategories. Although cellular processes and environmental information processing both have only one subclass, 2,479 and 3,288 single genes were enriched respectively (Figure 2D). In addition, the transcription factor families in which the unigenes belong were classified and counted. The results showed that the most genes belonged to the MYB family, followed by the bHLH family (Supplementary Table 4).
Identification and gene ontology analysis of differentially expressed genes
We performed volcano plots for the differential gene expression of the four comparison groups. Compared with CK, EBL group had 6,618 genes up-regulated and 3,647 genes down-regulated, respectively. In the CK vs. SALT comparison, there were 13,830 different genes, including 7,968 up-regulated genes and 5,862 down-regulated genes. There were 10,537 differential genes in the groups of CK and SE, including 6,409 up-regulated genes and 4,128 down-regulated genes. In the SALT vs. SE comparison, 117 DEGs were obtained, containing 58 upregulated genes and 59 downregulated genes (Figure 3).
Figure 3. Differential genes between groups in volcano maps. (A) The distribution of differential genes between CK and EBL; (B) The distribution of differential genes between CK and SALT; (C) The distribution of differential genes between CK and SE; (D) The distribution of differential genes between SALT and SE. Up-regulated genes were shown in red and down-regulated genes were shown in blue.
Based on the results, we selected the top 20 up-regulated and 20 down-regulated genes with the smallest Q-value in each comparison group for Gene Ontology (GO) analysis. According to GO annotation, differentially expressed genes (DEGs) are divided into three functional categories: molecular function (MF), cellular component (CC), and biological process (BP). And each functional category contains multiple levels of subcategories. In the CK vs. EBL comparison and the CK vs. SALT comparison, the most enriched GO term among the up-regulated genes were DNA-binding transcription factor activity in MF and protein-chromophore linkage in BP, respectively. The term with the highest number of down-regulated genes for both comparison groups was chloroplast in CC (Figures 4A,B). In the CK vs. SE comparison, chloroplast thylakoid membrane in CC was the most in up-regulated genes. Among down-regulated genes, carbohydrate metabolic process in BP and chloroplast in CC occupied the most by the same count (Figure 4C). In the SALT vs. SE comparison, the most abundant term in up-regulated genes was quinone binding in MF. The terms in down-regulated genes were actin filament depolymerization, sterol biosynthetic process in BP, actin cytoskeleton in CC, actin binding and oxidoreductase activity, acting on the CH-OH group of donors, NAD or NADP as acceptor in MF, respectively (Figure 4D).
Figure 4. Gene ontology (GO) classification of genes with significant differences between groups. (A) GO classification of differential genes between CK and EBL; (B) GO classification of differential genes between CK and SALT; (C) GO classification of differential genes between CK and SE; (D) GO classification of differential genes between SALT and SE. The top 20 up-regulated genes between two groups were shown in orange and the top 20 down-regulated genes were shown in blue.
Trend analysis of differentially expressed genes among different comparisons
A bidirectional clustering heatmap was created based on the 117 DEGs of SALT vs. SE comparison. DEGs with similar expression patterns were divided into multiple clusters using the k-means clustering method. Each cluster contained a heatmap and a trend graph. There were CL1355.Contig9, Unigene25288, Unigene5787, CL4441.Contig1, CL873.Contig36 in cluster 1 and the expression trends in the four groups were increased, and then decreased and finally increased again. There were 15 genes in cluster 2. Its trend was the same as cluster 1, but the fluctuation was smoother. Seven genes were in cluster 3 and the expression level of SALT was the highest among the four groups. There were 90 genes in cluster 4. It had the lowest gene expression among the four clusters, and the fluctuation of its trend was the smallest (Figure 5 and Supplementary Figure 2).
Figure 5. Trend analysis of DEG expression in four groups. (A) Heatmap of gene expression in cluster 1; (B) Trend graph of gene expression in cluster 1; (C) Heatmap of gene expression in cluster 2; (D) Trend graph of gene expression in cluster 2; (E) Heatmap of gene expression in cluster 3; (F) Trend graph of gene expression in cluster 3; (G) Heatmap of gene expression in cluster 4; (H) Trend graph of gene expression in cluster 4.
Plant hormone signal transduction
Plant hormones play important roles in growth and development of plants. Here, genes involved in plant hormone signal transduction pathways were differentially expressed in EBL vs. SE comparison (Figure 6). Auxin is an important regulatory hormone in plants. Several components of auxin signal transduction pathway were down-regulated after SE treatment compared to EBL, including one auxin-responsive protein IAA (Aux/IAA), one SAUR family protein (SAUR). There were three DELLA proteins (DELLA) down-regulated in the gibberellin pathway. Similarly, the salicylic acid pathway involved in disease resistance also contains only down-regulated genes. There was an up-regulated Arabidopsis histidine-containing phosphotransfer protein (AHP), two down-regulated genes of two-component response regulator ARR-B family (B-ARR) and two-component response regulator ARR-A family (A-ARR) in the cytokinin transduction pathway. In the ethylene transduction pathway related to fruit ripening and senescence, genes were all up-regulated. Most genes were also up-regulated in the abscisic acid signaling pathway. For jasmonic acid, three jasmonate ZIM domain-containing protein (JAZ) were down-regulated and three transcription factor MYC2 (MYC2) were up-regulated after SE treatment. In the brassinosteroid pathway, both brassinosteroid insensitive 1-associated receptor kinase 1 (BAK1) and protein brassinosteroid insensitive 1 (BRI1) contain one up-regulated and one down-regulated gene, respectively. Compared with SALT, SE had up-regulated genes in auxin and cytokinin signal transduction pathways, which indicated that the application of 24-epibrassinolide could improve plant growth (Supplementary Figure 3).
Figure 6. Plant hormone signal transduction in EBL vs. SE comparison. Up-regulated DEGs were indicated in red, down-regulated DEGs were indicated in green.
Verification of differentially expressed genes
We selected eight DEGs for transcriptome sequencing validation (Figure 7 and Supplementary Figure 4). Two genes were up-regulated in the three treated groups compared to the control group, and two genes were down-regulated in the three treated groups compared with the control group. Among the remaining four genes, two were expressed in the EBL and the SE, and the other two were expressed in the SALT or the SE. The expression levels of Unigene30391 and CL2028.Contig1 in the three treatment groups were higher than those of CK, and the expression levels of these two genes were the highest in EBL and SALT groups, respectively. CL1847.Contig5 and CL1052.Contig9 were both down-regulated in the three treatment groups, and the expression levels of these two genes were the lowest in EBL and SALT groups, respectively. Unigene4922 and Unigene7275 were only expressed in EBL and SE groups, and were significantly higher than CK and SALT. The expression of CL1002.Contig34 was the largest in SALT group, SE was the same as CK, and EBL was slightly larger than CK. CL7400.Contig4 was expressed in both SALT and SE, and their expression levels were two times and one time higher than that of CK, respectively.
Figure 7. Quantitative real-time PCR (qRT-pCR) analysis of several differentially expressed genes in tall fescue. (A) Relative expression of Unigene30391; (B) Relative expression of CL2028.Contig1; (C) Relative expression of CL1847.Contig5; (D) Relative expression of CL1052.Contig9; (E) Relative expression of Unigene4922; (F) Relative expression of Unigene7275; (G) Relative expression of CL1002.Contig34; (H) Relative expression of CL7400.Contig4. Bars with different letters indicate significant difference at P ≤ 0.05 (Duncan test).
Discussion
In saline environments, plants develop multiple strategies at the physiological and molecular levels to mitigate the damage caused by salt stress. The results of this study showed that NaCl treatment inhibited the growth and physiology of tall fescue, while the addition of EBL could confer increased tolerance to salt stress. Under the influence of sodium chloride, the plant height and chlorophyll content of tall fescue decreased, while the MDA content and electrolyte leakage increased. These growth and physiological changes were improved when 24-epibrassinolide was used. We found that SOD and CAT activities of tall fescue increased under salt stress, which was consistent with the results of Xu et al. (2013). After adding EBL, the activity of the antioxidant enzyme system was enhanced to improve the resistance of plants to salt stress. Wu et al. (2017) also found that exogenous 24-epibrassinolide increased the chlorophyll content of perennial ryegrass, enhanced the activities of SOD and CAT, and decreased the EL and MDA content. Similarly, Alam et al. (2019) found that EBL reduced the electrolyte leakage and hydrogen peroxide content, as well as the accumulation of Na+ of soybean. Moreover, Kolomeichuk et al. (2020) demonstrated that 24-epibrassinolide could increase peroxidase activity and improve the K+/Na+ ratio in potato leaves under NaCl treatment. These results suggest that 24-epibrassinolide may be an effective molecule for enhancing plant salt tolerance.
In order to understand the response mechanism of 24-epibrassinolide in alleviating salt stress, DEGs were subjected to GO analysis. In the CK vs. SALT comparison (Supplementary Table 5), most DEGs were annotated to cellular anatomical entity (6,045) in cellular component, catalytic activity (4,429) and binding (4,423) in molecular function, cellular process (3,990) and metabolic process (3,613) in biological process. The results were accorded with the GO annotation results of differential genes in tall fescue after 24 h of salt treatment by Amombo et al. (2018). Interestingly, the top five annotation categories in the CK vs. EBL comparison (Supplementary Table 5) were the same as CK vs. SALT. We randomly selected several genes for validation, and these genes were up-regulated or down-regulated under salt or 24-epibrassinolide treatment compared to the control. The GO annotations of these genes in the molecular function were calcium ion binding, chlorophyll binding, DNA binding and ATP binding, respectively. In addition to salt stress, most of the molecular functional subcategories of GO analysis are binding under drought, heat and lead stress (Li et al., 2017; Shu et al., 2020; Liu et al., 2021; Ouertani et al., 2021). This suggests that binding-related genes play a significant role in plant resistance to various environmental stresses.
There are various signals and pathways in plants to enhance salt tolerance (Yang and Guo, 2018; Zelm et al., 2020). Salt stress induces ionic and osmotic stress, leading to an increase in cytoplasmic Ca2+ content, and then the SOS pathway is triggered to export excess Na+ to reduce salt toxicity (Zhao et al., 2021). In addition, the high-affinity potassium transporter (HKT1), the mitogen-activated protein kinases (MAPKs) and various plant hormones are also important anti-stress signals (Pan et al., 2011; Fujita et al., 2013; Zhang et al., 2013; Khan et al., 2019; Nolan et al., 2020; Zhou et al., 2020). In this study, the plant hormone signal transduction, MAPK signaling pathway-plant and plant-pathogen interaction were the most enriched signal transduction pathways in the signal transduction of DEGs. Salt stress is regulated by the interaction of multiple plant hormones, and auxin is a key mediator in this process (Ribba et al., 2020). Auxin is perceived by Transport Inhibitor Response 1 (TIR1), which ubiquitinates Aux/IAA. Then the inhibition of AUXIN RESPONSE FACTORS (ARFs) is released, and auxin-induced genes expression are activated (Leyser, 2018). The results showed that genes related with auxin were down-regulated after SE treatment compared to EBL. In addition, among the other three hormones related to plant growth, gibberellin genes were all down-regulated, most genes of cytokinins and half of brassinosteroids were also down-regulated. However, there were up-regulated genes in the ethylene and jasmonic acid pathways associated with senescence. Growth-related genes were reduced, while defense-related genes were increased. This may be a protective mechanism adopted by plants to resist the damage caused by external stress (He et al., 2022).
Most transcription factors such as MYB, WRKY, bHLH, and NAC play important roles in plant adaptation to salt stress (Baillo et al., 2019). Many MYB proteins have been implicated in salt tolerance. MYB74 and MYB20 can be induced and the interaction between MPK4 and MYB42 can be enhanced by NaCl treatment (Wang et al., 2021). Among the transcription factor families related to plant stress resistance, the bHLH family is the second largest family after the MYB family. bHLH is an indispensable factor in the salt tolerance of various plants (Sun et al., 2018). The top six transcription factors that we annotated in the TF database with the highest number were MYB, bHLH, AP2-EREBP, WRKY, FAR1, NAC. Moreover, the genes annotated to these six families were both up-regulated and down-regulated after NaCl and 24-epibrassinolide treatment, indicating that these transcription factors were actively involved in the salt response of tall fescue.
Conclusion
To better understand the effect of 24-epibrassinolide on the salt tolerance of tall fescue, we performed studies at the physiological and transcriptomic levels. Through RNA sequencing, a total of 129,866 unigenes were acquired in four groups of CK, EBL, SALT, and SE. These unigenes were annotated in seven functional databases. A large number of DEGs were obtained in different comparison groups. There were 10,265 and 13,830 DEGs in EBL and NaCl treatments compared with control, respectively. Most of them were annotated in metabolic process, binding, catalytic activity in GO database. Most of the DEGs were involved in plant signal transduction pathways, indicating that hormones actively mediated the salt response mechanism of tall fescue. At the physiological level, the plant height of tall fescue was significantly decreased and the level of membrane lipid peroxidation was increased under NaCl treatment. After the application of 24-epibrassinolide, the contents of MDA and EL in leaves were remarkably reduced, and the content of chlorophyll and the activities of antioxidant enzymes were enhanced. These results suggest that 24-epibrassinolide has non-negligible effects on the enhancement of plant salt tolerance.
Data availability statement
The data presented in the study are deposited in the China National Center repository, accession number PRJCA010198 (https://ngdc.cncb.ac.cn/bioproject/browse/PRJCA010198).
Author contributions
YC, JF, XY, and ZH conceived and designed the research. YC, YX, YG, YZ, and MC performed the experiments. YC analyzed the data and drafted the manuscript. YC, JF, and AK revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the Project of Forestry Science and Technology Innovation and Promotion of Jiangsu (Grant No. LYKJ [2021]09).
Acknowledgments
We thank the CAS Key Laboratory of Plant Germplasm Enhancement and Specialty Agriculture (Wuhan, China) for the seeds of Houndog 6 of this study.
Conflict of interest
ZH was employed by Hunan Tobacco Research Institute.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.976341/full#supplementary-material
Supplementary Figure 1 | Growth of tall fescue after 10 days of 24-epibrassinolide and salt stress treatments.
Supplementary Figure 2 | A bidirectional clustering heatmap of DEGs in the SALT vs. SE comparison.
Supplementary Figure 3 | Plant hormone signal transduction in SALT vs. SE comparison.
Supplementary Figure 4 | FPKM values of genes used for qRT-PCR analysis. (A) FPKM values of Unigene30391; (B) FPKM values of CL2028.Contig1; (C) FPKM values of CL1847.Contig5; (D) FPKM values of CL1052.Contig9; (E) FPKM values of Unigene4922; (F) FPKM values of Unigene7275; (G) FPKM values of CL1002.Contig34; (H) FPKM values of CL7400.Contig4.
References
Alam, P., Albalawi, T. H., Altalayan, F. H., Bakht, M. A., Ahanger, M. A., Raja, V., et al. (2019). 24-Epibrassinolide (EBR) confers tolerance against NaCl stress in soybean plants by up-regulating antioxidant system, ascorbate-glutathione cycle, and glyoxalase system. Biomolecules 9:640. doi: 10.3390/biom9110640
Amombo, E., Li, X. N., Wang, G. Y., An, S., Wang, W., and Fu, J. M. (2018). Comprehensive transcriptome profiling and identification of potential genes responsible for salt tolerance in tall fescue leaves under salinity stress. Genes 9:466. doi: 10.3390/genes9100466
Anders, S., and Huber, W. (2010). Differential expression analysis for sequence count data. Genome Biol. 11:R106. doi: 10.1186/gb-2010-11-10-r106
Anwar, A., Bai, L. Q., Miao, L., Liu, Y. M., Li, S. Z., Yu, X. C., et al. (2018). 24-Epibrassinolide ameliorates endogenous hormone levels to enhance low-temperature stress tolerance in cucumber seedlings. Int. J. Mol. Sci. 19:2497. doi: 10.3390/ijms19092497
Arora, P., Bhardwaj, R., and Kanwar, M. K. (2012). Effect of 24-epibrassinolide on growth, protein content and antioxidative defense system of Brassica juncea L. subjected to cobalt ion toxicity. Acta Physiol. Plant. 34, 2007–2017. doi: 10.1007/s11738-012-1002-2
Avalbaev, A., Yuldashev, R., Fedorova, K., Petrova, N., Fedina, E., Gilmanova, R., et al. (2021). 24-epibrassinolide-induced growth promotion of wheat seedlings is associated with changes in the proteome and tyrosine phosphoproteome. Plant Biol. 23, 456–463. doi: 10.1111/plb.13233
Baillo, E. H., Kimotho, R. N., Zhang, Z. B., and Xu, P. (2019). Transcription factors associated with abiotic and biotic stress tolerance and their potential for crops improvement. Genes 10:771.
Bajguz, A. (2000). Effect of brassinosteroids on nucleic acids and protein content in cultured cells of Chlorella vulgaris. Plant Physiol. Biochem. 38, 209–215. doi: 10.1016/S0981-9428(00)00733-6
Bhattarai, S., Fu, Y. B., Coulman, B., Tanino, K., Karunakaran, C., and Biligetu, B. (2021). Transcriptomic analysis of differentially expressed genes in leaves and roots of two alfalfa (Medicago sativa L.) cultivars with different salt tolerance. BMC Plant Biol. 21:446. doi: 10.1186/s12870-021-03201-4
Cao, Y. J., Wei, Q., Liao, Y., Song, H. L., Li, X., Xiang, C. B., et al. (2009). Ectopic overexpression of AtHDG11 in tall fescue resulted in enhanced tolerance to drought and salt stress. Plant Cell Rep. 28, 579–588. doi: 10.1007/s00299-008-0659-x
Chen, Y. X., Chen, Y. S., Shi, C. M., Huang, Z. B., Zhang, Y., Li, S. K., et al. (2018). SOAPnuke: A MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data. Gigascience 7, 1–6. doi: 10.1093/gigascience/gix120
Fan, J. B., Hu, Z. R., Xie, Y., Chan, Z. L., Chen, K., Amombo, E., et al. (2015). Alleviation of cold damage to photosystem II and metabolisms by melatonin in bermudagrass. Front. Plant Sci. 6:925. doi: 10.3389/fpls.2015.00925
Feng, Q. J., Song, S. R., Yang, Y., Amee, M., Chen, L., and Xie, Y. (2021). Comparative physiological and metabolic analyzes of two Italian ryegrass (Lolium multiflorum) cultivars with contrasting salinity tolerance. Physiol. Plant 172, 1688–1699. doi: 10.1111/ppl.13374
Food and Agriculture Organization of the United Nations (2015). Status of the world’s soil resources (SWSR) – main report. Available online at: http://www.fao.org/3/i5199e/i5199e.pdf (Accessed June 1, 2020)
Fujita, Y., Yoshida, T., and Yamaguchi-Shinozaki, K. (2013). Pivotal role of the AREB/ABF-SnRK2 pathway in ABRE-mediated transcription in response to osmotic stress in plants. Physiol. Plant. 147, 15–27. doi: 10.1111/j.1399-3054.2012.01635.x
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Trinity: Reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883
Hand, M. L., Cogan, N. O., Stewart, A. V., and Forster, J. W. (2010). Evolutionary history of tall fescue morphotypes inferred from molecular phylogenetics of the Lolium-Festuca species complex. BMC Evol. Biol. 10:303. doi: 10.1186/1471-2148-10-303
He, Z. H., Webster, S., and He, S. Y. (2022). Growth–defense trade-offs in plants. Curr. Biol. 32, R634–R639. doi: 10.1016/j.cub.2022.04.070
Hiscox, J., and Israelstam, G. (1979). A method for the extraction of chlorophyll from leaf tissue without maceration. Can. J. Bot. 57, 1332–1334. doi: 10.1139/b79-163
Hu, L. X., Li, H. Y., Pang, H. C., and Fu, J. M. (2012). Responses of antioxidant gene, protein and enzymes to salinity stress in two genotypes of perennial ryegrass (Lolium perenne) differing in salt tolerance. J. Plant Physiol. 169, 146–156. doi: 10.1016/j.jplph.2011.08.020
Khan, A., Kamran, M., Imran, M., Al-Harrasi, A., Al-Rawahi, A., Al-Amri, I., et al. (2019). Silicon and salicylic acid confer high-pH stress tolerance in tomato seedlings. Sci. Rep. 9:19788. doi: 10.1038/s41598-019-55651-4
Kim, H. S., Lee, B. Y., Won, E. J., Han, J., Hwang, D., Park, H. G., et al. (2015). Identification of xenobiotic biodegradation and metabolism-related genes in the copepod Tigriopus japonicus whole transcriptome analysis. Mar. Genomics 3, 207–208. doi: 10.1016/j.margen.2015.05.011
Kim, T. W., Guan, S. H., Sun, Y., Deng, Z. P., Tang, W. Q., Shang, J. X., et al. (2009). Brassinosteroid signal transduction from cell surface receptor kinases to nuclear transcription factors. Nat. Cell Biol. 11, 1254–1260. doi: 10.1038/ncb1970
Kolomeichuk, L. V., Efimova, M. V., Zlobin, I. E., Kreslavski, V. D., Murgan, O. K., Kovtun, I. S., et al. (2020). 24-Epibrassinolide alleviates the toxic effects of NaCl on photosynthetic processes in potato plants. Photosynth Res. 146, 151–163. doi: 10.1007/s11120-020-00708-z
Li, B., and Dewey, C. N. (2011). RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics 12:323. doi: 10.1186/1471-2105-12-323
Li, H. Y., Hu, T., Amombo, E., and Fu, J. M. (2017). Transcriptome profilings of two tall fescue (Festuca arundinacea) cultivars in response to lead (Pb) stress. BMC Genomics 18:145. doi: 10.1186/s12864-016-3479-3
Liu, F. W., Ali, T., and Liu, Z. (2021). Comparative transcriptomic analysis reveals the effects of drought on the biosynthesis of methyleugenol in Asarum sieboldii Miq. Biomolecules 11:1233. doi: 10.3390/biom11081233
Lou, Y. H., Zhao, P., Wang, D. L., Amombo, E., Sun, X., Wang, H., et al. (2017).)Germination, physiological responses and gene expression of tall fescue (Festuca arundinacea Schreb.) growing under Pb and Cd. PLoS One 12:e0169495. doi: 10.1371/journal.pone.0169495
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8
Miller, G., Suzuki, N., Ciftci-Yilmaz, S., and Mittler, R. (2010). Reactive oxygen species homeostasis and signalling during drought and salinity stresses. Plant Cell Environ. 33, 453–467. 02041.x doi: 10.1111/j.1365-3040.2009.02041.x
Munns, R., and Gilliham, M. (2015). Salinity tolerance of crops – what is the cost? New Phytol. 208, 668–673. doi: 10.1111/nph.13519
Munns, R., and Tester, M. (2008). Mechanisms of salinity tolerance. Annu. Rev. Plant Biol. 59, 651–681. doi: 10.1146/annurev.arplant.59.032607.092911
Nolan, T. M., Vukašinoviíc, N., Liu, D., Russinova, E., and Yin, Y. (2020). Brassinosteroids: Multidimensional regulators of plant growth, development, and stress responses. Plant Cell 32, 295–318. doi: 10.1105/tpc.19.00335
Ouertani, R. H., Arasappan, D., Abid, G., Chikha, M. B., Jardak, R., Mahmoudi, M., et al. (2021). Transcriptomic analysis of salt-stress-responsive genes in barley roots and leaves. Int. J. Mol. Sci. 22:8155. doi: 10.3390/ijms22158155
Pan, J., Zhang, M., Kong, X., Xing, X., Liu, Y., Zhou, Y., et al. (2011). ZmMPK17, a novel maize group D MAP kinase gene, is involved in multiple stress responses. Planta 235, 661–676. doi: 10.1007/s00425-011-1510-0
Pertea, G., Huang, X., Liang, F., Antonescu, V., Sultana, R., Karamycheva, S., et al. (2003). TIGR gene indices clustering tools (TGICL): A software system for fast clustering of large EST datasets. Bioinformatics 19, 651–652. doi: 10.1093/bioinformatics/btg034
Pfaffl, M. W. (2001). A new mathematical model for relative quantification in real-time RT–PCR. Nucleic Acids Res. 29:e45. doi: 10.1093/nar/29.9.e45
Ribba, T., Garrido-Vargas, F., and O’Brien, J. A. (2020). Auxin-mediated responses under salt stress: From developmental regulation to biotechnological applications. J. Exp. Bot. 71, 3843–3853. doi: 10.1093/jxb/eraa241
Shen, X. Y., Wang, Z. L., Song, X. F., Xu, J. J., Jiang, C. Y., Zhao, Y. X., et al. (2014). Transcriptomic profiling revealed an important role of cell wall remodeling and ethylene signaling pathway during salt acclimation in Arabidopsis. Plant Mol Biol. 86, 303–317. doi: 10.1007/s11103-014-0230-9
Shi, H. T., Jiang, C., Ye, T. T., Tan, D. X., Reiter, R. J., Zhang, H., et al. (2015). Comparative physiological, metabolomic, and transcriptomic analyses reveal mechanisms of improved abiotic stress resistance in bermudagrass [Cynodon dactylon (L). Pers.] by exogenous melatonin. J. Exp. Bot. 66, 681–694. doi: 10.1093/jxb/eru373
Shu, B. S., Wu, Y. X., Qu, M. Q., Pu, X. H., Wu, Z. Z., and Lin, J. T. (2020). Comparative transcriptomic analyses revealed genes and pathways responsive to heat stress in Diaphorina citri. Gene 727:144246. doi: 10.1016/j.gene.2019.144246
Sun, X., Wang, Y., and Sui, N. (2018). Transcriptional regulation of bHLH during plant response to stress. Biochem. Biophys. Res. Commun. 503, 397–401. doi: 10.1016/j.bbrc.2018.07.123
Tanveer, M., Shahzad, B., Sharma, A., and Khan, E. A. (2019). 24-Epibrassinolide application in plants: An implication for improving drought stress tolerance in plants. Plant Physiol. Biochem. 135, 295–303. doi: 10.1016/j.plaphy.2018.12.013
Wang, G. Y., Bi, A. Y., Amombo, E., Li, H. Y., Zhang, L., Cheng, C., et al. (2017). Exogenous calcium enhances the photosystem II photochemistry response in salt stressed tall fescue. Front. Plant. Sci. 8:2032. doi: 10.3389/fpls.2017.02032
Wang, W. S., Zhao, X. Q., Li, M., Huang, L. Y., Xu, J. L., Zhang, F., et al. (2016). Complex molecular mechanisms underlying seedling salt tolerance in rice revealed by comparative transcriptome and metabolomic profiling. J. Exp. Bot. 67, 405–419. doi: 10.1093/jxb/erv476
Wang, X. P., Niu, Y. L., and Zheng, Y. (2021). Multiple functions of MYB transcription factors in abiotic stress responses. Int. J. Mol. Sci. 22:6125. doi: 10.3390/ijms22116125
Wang, Z. Y., Hopkins, A., Lawrence, R., Bell, J., and Scott, M. (2004). Field evaluation and risk assessment of transgenic tall fescue (Festuca arundinacea) plants. Mol. Breed. Forage Turf 11, 367–379. doi: 10.1007/1-4020-2591-2_32
Wu, W. L., Zhang, Q., Ervin, E. H., Yang, Z. P., and Zhang, X. Z. (2017). Physiological mechanism of enhancing salt stress tolerance of perennial ryegrass by 24-epibrassinolide. Front. Plant Sci. 8:1017. doi: 10.3389/fpls.2017.01017
Xu, R., Yamada, M., and Fujiyama, H. (2013). Lipid Peroxidation and Antioxidative Enzymes of Two Turfgrass Species Under Salinity Stress. Pedosphere 23, 213–222. doi: 10.1016/S1002-0160(13)60009-0
Yang, Y. Q., and Guo, Y. (2018). Unraveling salt stress signaling in plants. J. Integr. Plant Biol. 60, 796–804. doi: 10.1111/jipb.12689
Yu, Z. P., Duan, X. B., Luo, L., Dai, S. J., Ding, Z. J., and Xia, G. M. (2020). How plant hormones mediate salt stress responses. Trends Plant Sci. 25, 1117–1130. doi: 10.1016/j.tplants.2020.06.008
Zelm, E. V., Zhang, Y. X., and Testerink, C. (2020). Salt tolerance mechanisms of plants. Annu. Rev. Plant Biol. 71, 403–433. doi: 10.1146/annurev-arplant-050718-100005
Zhang, D., Jiang, S., Pan, J., Kong, X., Zhou, Y., Liu, Y., et al. (2013). The overexpression of a maize mitogen-activated protein kinase gene (ZmMPK5) confers salt stress tolerance and induces defence responses in tobacco. Plant Biol. 16, 558–570. doi: 10.1111/plb.12084
Zhao, S. S., Zhang, Q. K., Liu, M. Y., Zhou, H. P., Ma, C. L., and Wang, P. P. (2021). Regulation of plant responses to salt stress. Int J Mol Sci. 22:4609. doi: 10.3390/ijms22094609
Zhou, J. H., Li, Z. Y., Xiao, G. Q., Zhai, M. J., Pan, X. W., Huang, R. F., et al. (2020). CYP71D8L is a key regulator involved in growth and stress responses by mediating gibberellin homeostasis in rice. J. Exp. Bot. 71, 1160–1170. doi: 10.1093/jxb/erz491
Keywords: 24-epibrassinolide, salt stress, tall fescue, physiological responses, transcriptomic analysis, gene expression
Citation: Chen Y, Xiang Y, Hu Z, Gao Y, Zhang Y, Chen M, Khaldun ABM, Yan X and Fan J (2022) Transcriptomic profiling revealed the role of 24-epibrassinolide in alleviating salt stress damage in tall fescue (Festuca arundinacea). Front. Plant Sci. 13:976341. doi: 10.3389/fpls.2022.976341
Received: 23 June 2022; Accepted: 29 August 2022;
Published: 22 September 2022.
Edited by:
Maofeng Chai, Qingdao Agricultural University, ChinaCopyright © 2022 Chen, Xiang, Hu, Gao, Zhang, Chen, Khaldun, Yan and Fan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xuebing Yan, yxbbjzz@163.com; Jibiao Fan, fanjibiao11@126.com