Physiological and transcriptome analyses for assessing the effects of exogenous uniconazole on drought tolerance in hemp (Cannabis sativa L.)
grid.412064.50000 0004 1808 3449College of Agronomy, Heilongjiang Bayi Agricultural University, Daqing, 163316 China
grid.494628.50000 0004 1760 1486Daqing Branch of Heilongjiang Academy of Sciences, Daqing, 163319 China
grid.411846.e0000 0001 0685 868XCollege of Coastal Agricultural Sciences, Guangdong Ocean University, Zhanjiang, 524088 China
grid.411846.e0000 0001 0685 868XShenzhen Research Institute of Guangdong Ocean University, Shenzhen, 518108 China
grid.495707.80000 0001 0627 4537Institute of Cereal Crops, Henan Academy of Agricultural Sciences, Zhengzhou, 450002 China
Abstract
Uniconazole (S-(+)-uniconazole), a plant growth retardant, exerts key roles in modulating growth and development and increasing abiotic stress tolerance in plants. However, the underlying mechanisms by which uniconazole regulates drought response remain largely unknown. Here, the effects of exogenous uniconazole on drought tolerance in hemp were studied via physiological and transcriptome analyses of the drought-sensitive industrial hemp cultivar Hanma No. 2 grown under drought stress. Exogenous uniconazole treatment increased hemp tolerance to drought-induced damage by enhancing chlorophyll content and photosynthesis capacity, regulating activities of enzymes involved in carbon and nitrogen metabolism, and altering endogenous hormone levels. Expression of genes associated with porphyrin and chlorophyll metabolism, photosynthesis-antenna proteins, photosynthesis, starch and sucrose metabolism, nitrogen metabolism, and plant hormone signal transduction were significantly regulated by uniconazole compared with that by control (distilled water) under drought stress. Numerous genes were differentially expressed to increase chlorophyll content, enhance photosynthesis, regulate carbon–nitrogen metabolism-related enzyme activities, and alter endogenous hormone levels. Thus, uniconazole regulated physiological and molecular characteristics of photosynthesis, carbon–nitrogen metabolism, and plant hormone signal transduction to enhance drought resistance in industrial hemp.
Introduction
Drought stress, an abiotic stress, has more important effects on crop yield and quality compared to other abiotic stresses, inducing a highly vulnerable state in plants because of high temperature and low water content1. The photosynthetic system is greatly damaged owing to drought-induced low leaf water content, membrane lipid peroxidation, damaged leaf microstructure, higher leaf cell membrane permeability, and decreased photosynthetic pigment content2. Photosynthesis, as an index of the photosynthetic system, regulates key processes in carbon metabolism. Drought stress reduces net photosynthetic rate and carbon fixation capacity, which affects carbon absorption and metabolism in plants3. Therefore, carbon metabolism must be tightly regulated to enhance drought resistance in plants by changing starch content, with extensive accumulation of soluble sugar and sucrose4, and inverting sucrose catalyzed by sucrose synthase (SS) and sucrose phosphate synthase (SPS)5. Nitrogen metabolism is closely related to photosynthesis and carbon metabolism and is significant for plants adapting to environmental changes. The activities of key enzymes involved in N metabolism may play a major role in plant photosynthetic adaptation under drought stress6. Tolerance to drought stress is also manifested as increased levels of drought-responsive proteins and the expression of genes involved in signal transduction7. Plants initiate a series of morphological, molecular, and physiological and biochemical changes such as hormone regulation and gene expression to tolerate drought stress8. It is thus vital to study the mechanisms underlying drought stress resistance to develop effective strategies for alleviating drought-induced damage and improving agricultural production.
Hemp (Cannabis sativa L.) has one of the longest histories among fibre and traditional economic crops. It provides major raw materials such as fibre, seeds, xyloid stems, and floral leaves for numerous traditional and innovative industries globally9–11. Industrial hemp has a delta-9-tetrahydrocannabinol content of ≤ 0.3%, according to international standards. Drought stress during the growth period of hemp can lead to an increased incidence of diseases and insect pests, slow growth, and delay fibre and seed maturation, which ultimately affects yield and quality11. Few studies have focused on the mechanisms underlying drought stress resistance in hemp12. Currently, physiological and transcriptome analyses are widely used in plant science to assess changes that occur under drought stress but have been rarely utilized to study such changes in hemp.
Plant growth regulators can be used to regulate plant growth and stress tolerance under conditions of abiotic stress. Uniconazole, as one of the most efficient plant growth regulators, alters different parameters involved in the growth and development of plants under environmental stress13. Uniconazole applied exogenously under stress increases resistance to freezing injury14, water deficit stress15, drought stress16, waterlogging-induced damage17, and salt stress18. Previous studies14–18 revealed that uniconazole improved net photosynthesis and transpiration rate, suggesting an increase in leaf water potential and chlorophyll (Chl) content. Uniconazole also has various functions in promoting proline and soluble sugar accumulation, strengthening the antioxidant defence system, stimulating protein levels and nitrate reductase (NR) activity, and altering endogenous hormone levels. These experiments were conducted for different plants, such as soybean, bean, winter rape, and datura. However, the response of hemp to uniconazole is unclear. Research in other plants has revealed that uniconazole produces a marked effect on plant hormone metabolism. It induced a reduction in the content of gibberellin-like substances, mainly by inhibiting the cytochrome P450 enzyme ent-kaurene oxidase, which catalyses the three oxidation steps of transformation from ent-kaurene to ent-kaurenoic acid19. Uniconazole also suppressed a major abscisic acid (ABA) catabolic enzyme, ABA 8’-hydroxylase20. In addition, Liu et al.21,22 performed transcriptome analysis to dissect uniconazole-regulated expression of genes involved in the starch metabolism pathway. Furthermore, transcriptome analysis was applied to identify the complex molecular mechanisms underlying root development after treatment with uniconazole23. The effects of uniconazole on key enzymes involved in plant hormone signal transduction, carbohydrate metabolism, and starch metabolism pathways had been investigated. However, the genes and metabolic pathways involved in the response of hemp to uniconazole treatment under drought stress are unknown.
In the present study, we investigated the influences of uniconazole on the physiological and genome-wide gene expression of Cannabis sativa L. under drought stress. Our findings potentially lead to a thorough understanding of the biological function and the molecular mechanisms underlying uniconazole activity in hemp under drought stress.
Results
Effect of uniconazole on Chl content and net photosynthetic rate under drought stress
Drought stress decreased Chl content, but uniconazole mitigated this downward trend prominently (Fig. 1). Chl a content in plants treated with uniconazole (DS) was 13.9%, 22.3%, and 27.6% higher than that in drought-treated (D) plants following exposure to drought stress for 4, 6, and 8 days, respectively (Fig. 1a). There was a significant increase in Chl b content of 26.5% in DS plants compared to D plants under drought stress for 8 days (Fig. 1b). The total Chl content of DS plants was 4.1–27.2% higher than that of D plants (Fig. 1c). At 4 and 8 days, carotenoid (Car) content rapidly improved by 11.8% and 39.3% in DS plants, respectively, compared with that in D plants (Fig. 1d).
Net photosynthetic rate (PN) was reduced in the leaves of industrial hemp subjected prolonged drought stress (Fig. 1e). Uniconazole treatment sharply increased PN of plants by 1.1–6.0-fold relative to D plants under drought stress (Fig. 1e).
Effect of uniconazole on endogenous hormone levels under drought stress
ABA content increased with the prolonged drought stress, but the difference was not significant (Fig. 4a). Under drought stress for 4 days, ABA content of DS plants apparently increased by 2.2% compared with that of D plants (Fig. 4a). The content of indole-3-acetic acid (IAA), gibberellic acid (GA3), and zeatin (ZT) markedly decreased by 60.7%, 7.2%, and 30.6%, respectively, under drought stress for 4 days compared with those in plants under drought stress for 2 days (Fig. 4b–d). Compared with untreated plants, exogenous uniconazole application apparently increased IAA content by 2.2-fold and 3.9-fold in plants under drought stress for 2 and 4 days, respectively (Fig. 4b). GA3 content in DS plants was 47.4% and 51.1% lower than that in D plants under drought stress for 2 and 4 days, respectively (Fig. 4c). Compared to plants without supplementation, uniconazole significantly enhanced ZT content by 33.2% and 74.1% under drought stress for 2 and 4 days, respectively (Fig. 4d).
Effect of uniconazole on the expression of differentially expressed genes (DEGs) under drought stress
Twelve samples were collected for transcriptome sequencing, the results from Supplementary Table S1 revealed that the selected reference genome assembly could meet the needs of information analysis. Identification of DEGs revealed 5736 DEGs between D2_vs_DS2 and D4_vs_DS4 (drought stress for 2 and 4 days/drought stress treated with uniconazole for 2 and 4 days). Among them, 1701 DEGs were commonly regulated across different treatment periods, accounting for 29.65%. Non-shared DEGs were 1071 and 2964 in number, accounting for 18.67% and 51.67%, respectively (Supplementary Fig. S1a). The results indicated that the exogenous application of uniconazole had already mobilized few genes to counteract the damage caused by drought stress for 2 days and promoted strong drought resistance ability after 4 days of drought stress.
The expression trends of 21 genes were identical between the results of RNA-Seq and qRT-PCR. Thus, the RNA-Seq data was reliable (Supplementary Fig. S1b, Supplementary Table S10).
Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichments at D_vs_DS was performed, with the top 20 metabolic pathways shown in Fig. 5. Pathways related to photosynthesis (KO: ko00195), starch and sucrose metabolism (KO: ko00500), photosynthesis-antenna proteins (KO: ko00196), porphyrin and Chl metabolism (KO: ko00860), nitrogen metabolism (KO: ko00910), plant hormone signal transduction (KO: ko04075), and carbon metabolism (KO: ko01200) were significantly enriched by uniconazole under drought stress.
DEGs involved in plant hormone signal transduction under drought stress
Plant hormone signal transduction plays a critical role in plant growth and abiotic stress tolerance. Among 65 DEGs associated with plant hormone signal transduction, 28 were regulated commonly to the two comparisons. These involved IAA, cytokinin (CTK), GA, ABA, brassinosteroid (BR), jasmonic acid (JA), and salicylic acid (SA) metabolic processes and signal transduction pathways. Auxin-related genes included one auxin transporter protein gene (AUX1; Cannabis_sativa_newGene_19036), two auxin-responsive protein genes (IAA; FN08651.1.g and FN27264.1.g), four SAUR family protein genes (SAUR; FN03063.1.g, FN03185.1.g, Cannabis_sativa_newGene_6221, and FN33877.1.g), and five auxin responsive GH3 gene family genes (GH3; FN00395.1.g, FN01004.1.g, FN15644.1.g, FN01628.1.g, and FN31094.1.g), all of which are involved in regulating the adaptive growth and development processes of plants, were induced by uniconazole in plants grown under drought stress. Four histidine kinase genes (AHK2_3_4; FN12857.1.g, FN18674.1.g, FN27301.1.g, and FN10311.1.g) and one two-component response regulator gene (ARR-A; FN19201.1.g), which are associated with biosynthesis and signal transduction of CKs, were induced by uniconazole in plants grown under drought stress. One DELLA protein (DELLA; Cannabis_sativa_newGene_9451), which is related to the signal transduction pathway of GA, was up-regulated by uniconazole in plants grown under drought stress. ABA-related genes including one protein phosphatase gene (PP2C; FN20151.1.g) and two serine/threonine-protein kinase genes (SnRK2; Cannabis_sativa_newGene_8004 and Cannabis_sativa_newGene_919) were down-regulated by uniconazole in plants grown under drought stress. Genes of BSK (FN03133.1.g), CYCD3 (FN12036.1.g), JAZ (Cannabis_sativa_newGene_14430 and Cannabis_sativa_newGene_18071), MYC2 (FN04245.1.g) were up-regulated, and two transcription factor TGA genes were down-regulated by uniconazole in plants grown under drought stress (Fig. 7c, Supplementary Table S8).
Discussion
When plants are subjected to drought stress, total Chl (e.g., Chla, Chl b, and total Chl) and Car levels may decrease24. Drought stress inhibits the synthesis of Chl a/b24, affects light energy absorption, and transmission and conversion of photosystem I (PSI) and photosystem II (PSII) as well as decreases the photosynthetic rate25. In our study, Chl content decreased under drought stress in industrial hemp (Fig. 1a–d). Our finding that PN was reduced in the leaves of industrial hemp seedlings under drought stress (Fig. 1e) was similar to that documented by Jhou et al.26. Reportedly, uniconazole application markedly improves contents of both Chl and Car and protects plants against drought27. Similar results were observed in this study (Fig. 1a–d). Many genes involved in porphyrin biosynthesis, such as those for glutamyl-tRNA reductase, protoporphyrinogen oxidase, and the Mg chelatase subunit, are regulated to provide better drought tolerance in plants28. Glutamyl-tRNA reductase, a key enzyme in plant porphyrin biosynthesis, affects Mg chelatase and enzymes involved in Chl synthesis29. Mg chelatase catalyses the first step of Chl synthesis30. Mg protoporphyrin IX monomethyl ester (oxidative) cyclase is a major enzyme involved in Chl biosynthesis, and induction of the corresponding genes improves Chl content and restores chloroplast structure31. Chlorophyll(ide) b reductase is a necessary enzyme for Chl b degradation, and inhibition of corresponding genes reduces Chl and Car degradation and promotes thylakoid grana retention32. In the present study, the expression of genes corresponding to glutamyl-tRNA reductase, Mg chelatase subunit, Mg protoporphyrin IX monomethyl ester (oxidative) cyclase, protochlorophyllide reductase, and Chlorophyll(ide) b reductase were regulated by uniconazole under drought stress (Fig. 6a, Supplementary Table S2). The differential expression of porphyrin and Chl metabolism pathway genes suggested that contents of Chl and Car were affected by uniconazole under drought stress (Fig. 1a–d). Some reports showed that under drought stress, plants treated with uniconazole had a higher PN compared to that in untreated plants15 same to our study (Fig. 1e). Uniconazole application under drought stress can increase the levels of compounds involved in osmoregulation, lessen lipid peroxidation16, and protect cell integrity. These positive effects can improve absorption of light energy and thereby enhance photosynthesis capacity. Light harvesting Chl a/b-binding proteins (LHCs), as the most abundant membrane proteins, are involved in ABA signal transduction by regulating ROS (reactive oxygen species) homeostasis to improve tolerance to drought stress and grain yield33. In the present study, photosynthesis-antenna proteins pathway genes encoding LHCs were up-regulated by uniconazole under drought stress (Fig. 6b, Supplementary Table S3). These results suggest that uniconazole may modulate the drought resistance mechanisms of plants by changes in LHC gene expression. Photosystem II is a unique photosynthesis pigment-protein complex and a site for photosynthetic water oxidation34. The three photosystem II oxygen-evolving enhancer proteins and two photosystem II reaction centre subunits encoded by psbQ, psbP, psb28, and psbW are required for oxidation of water to O2 and thus are essential enzymes (Fig. 6c, Supplementary Table S4). Cytochrome c6 and plastocyanin are soluble electron carriers that transfer electrons to PSI35. In our study, one cytochrome c6 gene and one PC gene were up-regulated by uniconazole under drought stress. Moreover, uniconazole treatment under drought stress regulated seven PSI reaction centre subunit genes coding for proteins that catalysed the transfer of electrons from PC to FD through a series of electron transporters (Fig. 6c, Supplementary Table S4). FD receives electrons from PSI and transfers electrons to FD-NADP reductase (FNR), FNR induces the transfer of electrons from reduced FD to NADP+, and generates reduced NADPH for CO2 fixation and other chloroplast metabolic processes in the Calvin cycle36. Our study showed that four FD genes and two FD-NADP reductase genes were induced by uniconazole in plants grown under drought stress (Fig. 6c, Supplementary Table S4). The differential expression of photosynthetic pathway genes indicated that photosynthesis plays a significant role in the response to drought stress by uniconazole.
Carbon and nitrogen metabolisms, as the two most basic metabolic processes in plant physiology, affect the formation, transformation, and transportation of photosynthetic products37. Plants grown under drought stress are characterized by reduced Chl contents and decreased photosynthesis, which inhibits starch synthesis. In addition, starch was converted into soluble sugar to maintain cell osmotic pressure38. Some studies showed that the content of soluble sugar and sucrose increases but of starch content decreases39, which may be related to the increase of carbohydrate metabolism activity under drought stress. Analogous consequences were observed in our study: soluble sugar and sucrose content improved, whereas starch content decreased with prolonged drought stress (Fig. 2a–c). SPS, SS, S-AI, and NI are key enzymes involved in starch biosynthesis. Zahoor et al.39 reported an increase in SPS and SS activities but a decrease in S-AI activity in response to drought stress. In other reports, SS activity was lowered40 and SPS activity was enhanced41 by stress. These results were replicated in our studies, indicating that S-AI, NI, and SS activities were significantly reduced but SPS activity was increased in drought-treated plants (Fig. 2d–g). SS can catalyse both the synthesis and decomposition of sucrose. Considering our results, SS might play a role in sucrose decomposition. Improved SPS activity changed the distribution of carbon assimilates, resulting in increased contents of soluble sugar and sucrose, which are important for regulating cell osmotic pressure to resist drought stress42. Uniconazole markedly raised the soluble sugar content during waterlogging43, a similar finding was observed here (Fig. 2a). It was reported that exogenous spraying with uniconazole increased starch accumulation, starch granule content, and key enzyme activities involved in starch synthesis 22, same effects were noted in our study (Fig. 2b). An increase in sucrose content and a decrease in starch content were induced by drought stress, indicating that drought promoted the distribution of photosynthetic products to sucrose. Exogenous uniconazole application improved starch content, whereas it reduced the sucrose content of industrial hemp leaves under drought stress (Fig. 2b, c). A reasonable explanation is that uniconazole could enhance the photosynthetic capacity of drought-stressed industrial hemp leaves, which possibly adjusted the synthesis of sucrose and starch influenced by drought stress. Our results revealed that in plants under drought stress, uniconazole supplementation increased S-AI, NI, and SS activities compared with no supplementation (Fig. 2d–f), suggesting that uniconazole promotes sucrose decomposition and regulates photo-assimilation and translocation. SPS affects the distribution of photosynthetic products between starch and sucrose, and its activity is negatively correlated with starch accumulation and positively with sucrose accumulation44. We observed similar results: uniconazole decreased SPS activity, increased starch content, and reduced sucrose content under drought stress (Fig. 2b–g). Drought stress has been reported to influence SPS transcripts levels, indicating that SPS enzymes play central roles in starch and sucrose metabolism pathways45. In the present study, five SPS genes were regulated by uniconazole under drought stress (Supplementary Table S5). Previous studies have shown that sucrose synthase genes are regulated by drought stress, improving hyperosmotic stress tolerance of plants in relation to the accumulation of starch and sucrose metabolism proteins46. The present study described that two SuSy genes were differentially expressed in drought stress induced by uniconazole (Fig. 7a, Supplementary Table S6). In addition, INV (1), scrK (4), bglX (8), GN4 (3), glgC (2), PYG (3), GBE1 (2), TPS (7), TREH (3), and beta-amylase (1) genes were differentially expressed (Fig. 7a, Supplementary Table S6). The differential expression of these genes suggested that the starch and sucrose metabolic pathways were involved in the responses of industrial hemp leaves to drought stress under uniconazole treatment.
Nitrogen metabolism has great significances in plant growth, ecosystem structure and function47. Plants absorb nitrate nitrogen (NO3−) and ammonium nitrogen (NH4+) from soil, which are then converted into amino acids, proteins, and nucleotides by the key enzymes of NR, nitrite reductase, GS, GOGAT, and GDH47 that are crucial to nitrogen assimilation. Reportedly, the activities of NR, GS, and GOGAT in the functional leaves of plants were significantly reduced by drought stress, which affected the assimilation of NO3− and NH4+, then inhibited nitrogen metabolism in the leaves of plants40. Our results were similar to these: NR, GS, and GOGAT activities of industrial hemp leaves were lower in drought stressed seedlings than in controls (Fig. 3a–c). Hessini et al.48 reported that GS levels decreased but GDH activity increased and proline levels improved, suggesting that GDH might be implicated in the synthesis of compatible solute under drought stress. In our results, the activity of GDH increased in drought stressed plants (Fig. 3d). Previous studies with uniconazole revealed that it has important roles in regulating levels of metabolites and promoting an increase in NR activity49. These results are similar to those seen in our study, where uniconazole showed increased NR activity under drought stress (Fig. 3a). Overexpression of the NR gene delayed the decrease in NR activity under drought stress, and the plants rapidly recovered N assimilation after rehydration50. In our study, one nitrate/nitrite transporter gene was up-regulated significantly by uniconazole under drought stress (Fig. 7b, Supplementary Table S7). Previous studies indicated that the regulation of a nitrate transporter gene implicated in N utilization decreased NO3− accumulation of guard cells and enhanced drought tolerance of plants51. In plants under drought stress, exogenous uniconazole application enhanced GS activity compared with no application (Fig. 3b), which is consistent with the results of a previous study52. The GS/GOGAT pathway is the major system for utilizing NH4+ in plants, in our study, exogenous uniconazole application increased the activities of GS and GOGAT to relieve the toxic effect of NH4+ over-accumulation on plant seedlings by converting NH4+ to glutamate53. Previous studies showed that the transcription of glutamine synthetase genes and glutamate synthase genes in different tissues was up-regulated or down-regulated with changes in nitrogen levels54. Overexpression of glutamine synthetase OsGS1;1 and OsGS2 genes could strengthen the tolerance of plants to adverse abiotic stresses at the seedling stage in rice55. Lu et al.56 revealed that the suppression of glutamate synthase genes affected nitrogen assimilation in plants and had significant roles in carbon and nitrogen metabolism. In our studies, one glutamine synthetase gene (glnA) and two glutamate synthase genes (GLT1) were regulated by uniconazole under drought stress (Fig. 7b, Supplementary Table S7). These results revealed that uniconazole promoted the GS/GOGAT cycle, which is important in nitrogen metabolism. Reportedly, GS activity was reduced and GDH activity was increased under conditions of adversity stress, suggesting that glutamate and proline were produced by GDH48,57. Our study showed that DS plants exhibited lower GDH activity than D plants (Fig. 3d). Lightfoot et al.58 found that gdhA transgenic plants had improved seed germination rate, growth level, and grain biomass production and showed enhanced plant resistance by promoting nitrogen uptake59. In our results, four glutamate dehydrogenase genes were induced by uniconazole under drought stress (Fig. 7b, Supplementary Table S7). In addition, four CA genes (cynT) were differentially expressed under drought stress treated with uniconazole (Fig. 7b, Supplementary Table S7), which may catalyse the conversion between CO2 and HCO3−, and play important functions in photosynthesizing cells of plants under environmental stress. The differential expression of nitrogen metabolism pathway genes revealed that nitrogen metabolism had an effect on plants under drought stress treated with uniconazole.
A series of changes in plant endogenous hormones under drought stress reflects the adaptability of plants to stress. Moreover, drought resistance in plants is achieved by the complex coordination of various hormones, rather than by a single hormone. ABA and ZT, as critical endogenous hormones, take part in regulating the abiotic stress response. Some studies showed that drought stress substantially enhanced endogenous hormone ABA level but reduced ZT content, indicating that ABA and ZT jointly controlled plant photosynthetic performance in response to drought stress60. Our results similarly showed that ABA content was increased but ZT level was decreased under drought stress (Fig. 4a, d). Some reports revealed that ABA and ZT levels significantly increased in plants grown under abiotic stress and treated with uniconazole14, promoting the accumulation of proline and abiotic stress-related amino acids13. We noticed that uniconazole application notably increased ABA and ZT contents under drought stress for 4 days (Fig. 4a, d). Previous studies showed that histidine kinases (AHKs) positively affected drought stress responses through the ABA signalling pathway61. Two-component response regulator genes are also employed by plants for stress adaptation62. A previous study revealed that ABA signalling key kinases interacted with the CTK signalling regulator ARR in response to drought stress63. Liu et al.21,22 reported that uniconazole regulated AHKs and type-B ARRs to affect Chl degradation. In our studies, histidine kinase genes (AHK2_3_4) and two-component response regulator gene (ARR-A) were induced by uniconazole under drought stress (Fig. 7c, Supplementary Table S8). Protein phosphatase genes have been isolated and cloned from plants and may play key roles in signal transduction in drought stress64. Serine/threonine-protein kinases are the key regulatory factors of ABA signalling pathways, which are involved in abiotic stress adaptation in plants65. In our studies, protein phosphatase (PP2C) and serine/threonine-protein kinase (SnRK2) were down-regulated by uniconazole as an adaptation to drought stress (Fig. 7c, Supplementary Table S8). Drought stress blocks IAA signal transduction, decreases IAA biosynthesis, and markedly reduces IAA content66. In our study, IAA content was decreased with prolonged drought stress (Fig. 4b). Exogenous uniconazole treatment increased IAA content compared with no treatment in seedlings under drought stress (Fig. 4b). A logical explanation is that uniconazole inhibited IAA oxidase activity under drought stress, thus increasing IAA concentration. Similar conditions have been identified wherein uniconazole-treated plants had increased IAA content under water deficit stress15. It has been reported that auxin transporter genes (AUX/LAX)67, auxin-responsive protein (IAA and SAUR)68, and indole-3-acetic acid-amido synthetase gene (GH3)69 are involved in the auxin signalling pathway and play important roles in drought tolerance. In the present study, auxin-related genes (AUX1, IAA, SAUR, and GH3) were induced by uniconazole in response to drought stress (Fig. 7c, Supplementary Table S8). Gibberellin plays a key role in different stages of plant growth, such as seed germination, internode elongation, and flower and fruit formation and is involved in regulating abiotic stress processes70. The GA3 level decreased under drought stress, restricting the growth of plants15. In our studies, the content of endogenous GA3 reduced during the final d of drought stress (Fig. 4c). Uniconazole is an inhibitor of gibberellin biosynthesis, reports also found that endogenous GA3 content decreased in uniconazole-treated plants21,22 under stress15. These results were similar to those obtained in our study, where supplemental uniconazole under drought stress decreased the GA3 content (Fig. 4c). DELLA proteins, members of the plant-specific GRAS family, negatively regulate plant growth71. In previous studies, abiotic stress promoted the accumulation of DELLA proteins and then reduced endogenous hormone GA levels to inhibit plant growth and enhance plant tolerance to stress70. However, a study had shown that the overexpression of DELLA proteins (S-della) increased stomatal sensitivity to ABA and reduced water loss72. Uniconazole affected not only the GA biosynthesis pathway but also related genes in this pathway23. A previous study showed that uniconazole regulated the expression of DELLA proteins to affect plant growth22. In our study, one DELLA protein was up-regulated by uniconazole under drought stress (Fig. 7c, Supplementary Table S8). In addition, BSK (1) and CYCD3 (1) genes involved in brassinosteroid biosynthesis; JAZ (2) and MYC2 (1) genes involved in jasmonic acid metabolism; and two transcription factor TGA genes involved in salicylic acid biosynthesis were differentially expressed (induced) by uniconazole (Fig. 7c, Supplementary Table S8). Numerous genes were differentially expressed, showing that plant hormone signal transduction pathways are involved in the response of industrial hemp leaves in plants grown under drought stress and treated with uniconazole.
Based on the above discussion, a hypothetical model of the effects of uniconazole on physiology and genome-wide gene expression in industrial hemp seedlings grown under drought stress is presented in Fig. 8.
Methods
Plant material and seedling treatments
‘Hanma No. 2’, an industrial hemp cultivar, was used in this study. The variety was bred and served at Daqing Branch of Heilongjiang Academy of Sciences (Daqing, China), and the seeds are available from the first author on reasonable request. Germinated seeds were sown in plastic pots with a height of 16 cm, mouth diameter of 20 cm, and bottom diameter of 12.5 cm filled with 2:1 (v/v) grass peat and sand. The properties of grass peat were as follows: maximum field capacity (FC), 46%; pH 6.3; organic matter content, 20.4%; available P, 13.7 mg/kg; NH4+, 315.3 mg/kg; and available K, 206 mg/kg. In total, 15 seedlings per pot were maintained after seed germination. Industrial hemp seedlings were grown using the weighing method, and an appropriate amount of water was added to maintain 70% of FC until the trifoliate stage in the mobile canopy. At the three-leaf stage, half of potted hemp plants (40 pots) were sprayed with 40 mg L−1 uniconazole solution73, and the others were sprayed with distilled water. At 48 h after spraying, the plants were subjected to drought stress for 0, 2, 4, 6, and 8 days. D and DS were used to represent drought stress after spraying with distilled water and drought stress after spraying with uniconazole, respectively. The soil moisture content at 0, 2, 4, 6, and 8 days was maintained at 32%, 28%, 23%, 18%, and 14%, respectively. For physiological experiments, the fully expanded third pair of leaves from the plant base was collected at each treatment; one biological repeat had twenty leaves, and four biological repeats were used. For endogenous hormone levels and RNA-seq, the completely expanded third pair of leaves from the base of plants subjected to drought stress for 2 and 4 days was harvested, and three biological repeats were used for D and DS treatments. The samples were frozen in liquid nitrogen immediately and stored at – 80 °C.
Chl content and net photosynthetic rate
Chl, including Chl a, Chl b, and total Chl, and Car contents were measured as reported by Arnon74 with a few modifications. Fresh leaves without veins (0.1 g) were soaked in ethanol and acetone (v/v, 1:1) until the leaves became colourless. Absorbances of the extracts were recorded at 470, 645, and 663 nm.
PN of the third pair of leaves was measured between 9:00 and 11:00 am using a portable photosynthesis system (LI-6400 XT, LI-COR, Inc. Lincoln, NE, USA) under light intensity of 1000 µmol m−2 s−1.
Endogenous levels of ABA, IAA, GA3, and ZT
Leaves (0.1 g) were ground in liquid nitrogen; then added into a 1 mL mixture of methanol, water, and acetic acid; and extracted overnight at 4 °C. The extracts were centrifuged at 8000×g for 10 min, and the supernatant was dried to the water phase with nitrogen. After adjusting the pH to 2–3 with citric acid, the solution was extracted three times with ethyl acetate. The ethyl acetate phases were combined and dried using a stream of nitrogen and then made up to a constant volume with methanol. The solution was filtered through a 0.22-μm organic filter and assessed using a liquid chromatography-tandem mass spectrometry system. The levels of the endogenous hormones ABA, IAA, GA3, and ZT in the leaves of hemp plants grown under D and DS for 2 and 4 days were analyzed using the method reported by Farrow and Emery80.
RNA-seq library preparation and sequencing
Total RNA was extracted from the fully developed third leaves of D2, D4, DS2, and DS4 using TRIzol reagent (Invitrogen, Carlsbad, CA, USA), and residual DNA was eliminated with DNase I (Fermentas, Vilnius, Lithuania). RNA concentration and quality were evaluated using a Nanophotometer Spectrophotometer (Implen, CA, United States) and an RNA Nano 6000 assay kit for the Agilent Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, CA, United States). The RNA sequencing library was constructed following the manufacture’s recommendation by using NEBNext Ultra RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). Equal amounts of total RNA from each sample were pooled together. Then the mRNAs were purified from total RNA by poly-T oligo-attached VAHTS mRNA Capture Beads (Vazyme BioTech, Nanjing, China). The mRNAs were fragmented to avoid priming bias when synthesizing cDNA by using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer. First-strand cDNAs were synthesized using random hexamer primer and M-MuLV Reverse Transcriptase (RNase H-). Second-strand cDNAs were subsequently synthesized with DNA Polymerase I and RNase H. Double-stranded cDNAs were purified by VATHS magnetic DNA clean beads (Vazyme BioTech). Then the End Repair Reaction Buffer and End Prep Enzyme Mix were added into the purified products for end repair and adenine (A) addition of 3'-end. The adaptor-ligated cDNAs were then added with USER Enzyme (NEB) for adaptor cleavage. Subsequently, library size selection was conducted using magnetic purification beads (Vazyme BioTech) to capture DNA fragments in a given size range, remove unwanted DNA fragments, and remove contaminants such as enzymes. Then PCR was performed with Phusion High-Fidelity DNA Polymerase, Universal PCR Primers and Index (X) Primer. At last, PCR products were purified with Vazyme VATHS DNA clean beads and library quality was assessed on the Agilent Bioanalyzer 2100 system. The RNA-seq library preparations were sequenced on an Illumina HiseqX-ten platform (San Diego, CA, USA) at Biomarker Technologies Co., Ltd. (Beijing, China), following the manufacturer’s instructions (Illumina).
Date processing and identification of DEGs
A large amount of raw data (raw reads) was screened with in-house perl scripts to obtain clean data (clean reads) by removing reads containing adapter and poly-N, and low quality reads. The GC content, Q30 and sequence duplication levels were calculated. High-quality clean reads were mapped to the C. sativa L. reference genome sequence with HISAT281.
Gene expression levels were represented using the fragments per kilobase of exon per million mapped read (FPKM) value of each transcript. The abundance of the same transcripts of D2_vs_DS2 and D4_vs_DS4 samples were compared to identify DEGs using DESeq282. This study compared the effects of uniconazole treatment on DEGs in the leaves of industrial hemp grown for 2 and 4 days under drought stress using fold change ≥ 2 and false discovery rate < 0.01 as the screening criteria.
KEGG pathway enrichment analysis
KEGG pathway enrichment analysis was performed using the KEGG pathway mapping tool (http://www.genome.jp/kegg)83–85 to analyze significantly enriched metabolic pathways or signal transduction pathways of DEGs.
qRT-PCR of mRNAs
Twenty-one genes showing different expression levels in RNA-Seq analyses were selected to confirm the results of RNA-seq with the method of quantitative real-time PCR (qRT-PCR) analyses. Total RNA was extracted and purified according to the above method. Subsequently, 1 μg of purified RNA was reverse-transcribed using the FastKing RT Kit (Tiangen BioTek, Beijing, China) to produce cDNA. The GAPDH gene was selected as the reference gene. The qRT-PCR system contained Power qPCR PreMix (Genecopoeia, Rockville, MD, USA), cDNA, and the forward and reverse gene-specific primers. The reaction procedure included 40 cycles of a denaturation step at 95 °C for 10 s, following by an annealing step at 60 °C for 40 s. The experiment consisted of three repetitions. Gene-specific primers were designed by Primer Premier 5 software (Premier Biosoft International, Palo Alto, CA, USA) and listed in Supplementary Table S9.
Data analysis
All experimental data were analyzed by analysis of variance using SPSS software version 20.0 (IBM Inc., Chicago, IL, USA) and presented as mean ± standard errors. Duncan’s multiple range test (P < 0.05) was applied to detect differences between means. All graphs were generated using OriginPro 9.1 software (OriginLab, Northampton, MA, USA).
The original RNA sequencing data were submitted to NCBI’s Gene Expression Omnibus (http://www.ncbi.nih.gov/geo/) under series number PRJNA635553.
All methods were performed in accordance with the relevant national or international guidelines/regulations/legislation.
Supplementary Information
Untitled section
Supplementary Information
The online version contains supplementary material available at 10.1038/s41598-021-93820-6.
Acknowledgements
The authors are grateful to Heilongjiang academy of sciences for supporting. This research was supported by Heilongjiang Academy of Sciences CXJQ2018DQ01 and National Key R&D Projects 2016YFC0501207-03.
Competing interests
The authors declare no competing interests.