Next Article in Journal
Lipid Alterations and Metabolism Disturbances in Selected Inflammatory Skin Diseases
Previous Article in Journal
Recent Advances in Genome-Editing Technology with CRISPR/Cas9 Variants and Stimuli-Responsive Targeting Approaches within Tumor Cells: A Future Perspective of Cancer Management
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 8
10.3390/ijms24087051
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Transcriptome-Based Construction of the Gibberellin Metabolism and Signaling Pathways in Eucalyptus grandis × E. urophylla, and Functional Characterization of GA20ox and GA2ox in Regulating Plant Development and Abiotic Stress Adaptations

by
Wenfei Wu
1,2,
Linhui Zhu
1,2,
Pan Wang
1,2,
Yuwu Liao
1,2,
Lanjuan Duan
1,2,
Kai Lin
1,2,
Xin Chen
1,2,
Lijie Li
1,2,
Jiajing Xu
1,2,
Hao Hu
1,2,
Zeng-Fu Xu
1,2,* and
Jun Ni
1,2,*
1
State Key Laboratory for Conservation and Utilization of Subtropical Agro-Bioresources, College of Forestry, Guangxi University, Nanning 530004, China
2
Guangxi Key Laboratory of Forest Ecology and Conservation, College of Forestry, Guangxi University, Nanning 530004, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(8), 7051; https://doi.org/10.3390/ijms24087051
Submission received: 10 March 2023 / Revised: 5 April 2023 / Accepted: 9 April 2023 / Published: 11 April 2023
(This article belongs to the Section Molecular Plant Sciences)
Browse Figures
Versions Notes

Abstract

:
Gibberellins (GAs) are the key regulators controlling plant growth, wood production and the stress responses in perennial woody plants. The role of GA in regulating the above-mentioned processes in Eucalyptus remain largely unclear. There is still a lack of systematic identification and functional characterization of GA-related genes in Eucalyptus. In this study, a total of 59,948 expressed genes were identified from the major vegetative tissues of the E. grandis × E. urophylla using transcriptome sequencing. Then, the key gene families in each step of GA biosynthesis, degradation and signaling were investigated and compared with those of Arabidopsis, rice, and Populus. The expression profile generated using Real-time quantitative PCR showed that most of these genes exhibited diverse expression patterns in different vegetative organs and in response to abiotic stresses. Furthermore, we selectively overexpressed EguGA20ox1, EguGA20ox2 and EguGA2ox1 in both Arabidopsis and Eucalyptus via Agrobacterium tumefaciens or A. rhizogenes-mediated transformation. Though both Arabidopsis EguGA20ox1- and EguGA20ox2-overexpressing (OE) lines exhibited better vegetative growth performance, they were more sensitive to abiotic stress, unlike EguGA2ox1-OE plants, which exhibited enhanced stress resistance. Moreover, overexpression of EguGA20ox in Eucalyptus roots caused significantly accelerated hairy root initiation and elongation and improved root xylem differentiation. Our study provided a comprehensive and systematic study of the genes of the GA metabolism and signaling and identified the role of GA20ox and GA2ox in regulating plant growth, stress tolerance, and xylem development in Eucalyptus; this could benefit molecular breeding for obtaining high-yield and stress-resistant Eucalyptus cultivars.

1. Introduction

Eucalyptus is one of the fastest growing tree species. It is widely planted in many tropical and subtropical regions for pulp and paper production. Eucalyptus hybrids (e.g., E. grandis × E. urophylla) combine fast growth with excellent wood quality and relatively high stress resistance; thus, they have been extensively used in commercial plantations, mostly in clonal forestry worldwide [1,2,3]. The excellent performance of Eucalyptus hybrids is the consequence of interactions between good genotype and environment [2]. Thus, ‘elite’ clones are the prerequisites for the Eucalyptus industry. Recent advancements in genomic and transcriptomic technologies have led to the sequencing of the E. grandis genome and the transcriptomes of several other Eucalyptus species [4,5,6,7,8,9]. Theoretically, the important economic characteristics of Eucalyptus could be “on-target” increased by using molecular breeding strategies. Nevertheless, genetic studies and molecular breeding in Eucalyptus are still lacking. The transcriptomic or genomic investigation of E. grandis × E. urophylla would help to clarify the formation of heterosis in the hybrid at the molecular level.
Gibberellins (GAs), a class of tetracyclic diterpenoid phytohormones, play an important role in many aspects of plant growth and development, including seed germination, wood formation, flowering, and stress adaptation [10,11,12]. In the past few decades, the metabolism and signaling of GAs have been well studied in many plant species, especially in Arabidopsis and rice. GA biosynthesis is a complicated pathway that is mainly catalyzed by seven key steps: ent-copalyl diphosphate synthase (CPS), ent-kaurene synthase (KS), ent-kaurene oxidase (KO), ent-kaurenoic acid oxidase (KAO), GA13-oxidase (GA13ox), GA20-oxidase (GA20ox), and GA3-oxidase (GA3ox) [10]. During this synthesis route, abundant GA intermediates (GA12, GA15, GA24, GA9, GA53, GA44, GA19, and GA20) are produced and then catalyzed into the bioactive forms (GA1, GA3, and GA4) at the final catalytic step by GA3ox. Although many gene families are involved in late GA synthesis, GA20ox seems to be the limiting step for GA synthesis. GA20ox2 was previously regarded as the “green revolution” gene, the mutation of which caused a relative dwarf phenotype in rice, which led to yield improvement [13,14]. In hybrid aspen and tobacco, overexpression of GA3ox did not cause a significant increase in bioactive GAs [15,16]. In pea, overexpression of GA3ox influenced GA biosynthesis, but it was timing- or localization-dependent [17]. Previous studies also showed that overexpression of the early-step genes of GA biosynthesis, such as CPS and KS, did not affect the bioactive GAs but rather the GA intermediates [18,19]. Thus, it is apparent that GA20ox, but not GA3ox or other catalytic enzymes, are the rate-limiting enzymes in GA synthesis. Therefore, GA20oxs are the key targets for genetic engineering aimed at manipulating the important economic traits that are governed by GAs.
The deactivation of bioactive GAs is regulated by several distinct groups of enzymes, including GA2ox, CYP714D1, methyl transferase, and GA glucosyl transferase [10]. Among these, GA2ox-mediated catalysis of active GAs into nonactive forms has been the most universal and well-studied step in recent years [20,21]. GA2oxs belong to the 2OG-FE (II) oxygenase superfamily. Similar to the GA20ox and GA3ox families, GA2oxs are also encoded by a small group of genes, which can be divided into two major families, the members of which act primarily on C19-GA and C20-GA substrates [10,22]. To date, nine GA2ox genes have been identified in Arabidopsis [21], four in rice [23], eleven in Populus [24], five in tomato [25], seven in peach [26], and eight in Liriodendrom chinense [27]. In L. chinense and peach, different members of the GA2ox family exhibited diverse expression pattern in different plant tissues [26,27]. The expression of several GA2oxs, such as LcGA2ox1, LcGA2ox4 and LcGA2ox7, was specially up-regulated by salt, osmotic, or chilling stresses [27]. Overexpression of GA2ox in tobacco reduced the bioactive GA content and thus induced the dwarfism phenotype [26]. As GA2oxs have essential functions in counterbalancing GA concentration, the flexible regulation of GA levels, achieved by the temporal and spatial expression of multiple GA2ox family members, is very important for maintaining regular plant growth and development and in responses to environmental clues [21].
To date, abscisic acid (ABA), jasmonic acid (JA), salicylic acid (SA), strigolactone (SL), and melatonin are the most well-studied phytohormones in the regulation of abiotic and biotic stress [28,29,30,31,32,33]. Although a few studies have reported that the biosynthesis and metabolism of GA are differentially expressed in response to biotic and abiotic stresses [34,35,36,37], the molecular mechanisms remain elusive. In Populus, several GA2oxs are specifically highly expressed in the roots, and overexpression of GA2ox modulates biomass accumulation and lateral root formation via interactions with other phytohormones [38,39]. In Arabidopsis, overexpression of cucumber GA20ox1 promoted the primary root growth and increased the lateral root number [39]. Root architecture, including root length, number of lateral roots, and root hairs, is an important trait for plants to adapt to abiotic stresses [40]. Thus, GA-modulated plant root growth might also contribute to stress adaptation. Moreover, under abiotic stresses, the level of bioactive GAs could be immediately downregulated by inducing the expression of different GA2oxs [34,37,41]. Recent reports have shown that stress responsive AtGA2ox9 and AtGA2ox10 directly regulate plant adaptation to cold stress in Arabidopsis [42]. Thus, the stress-responsive modulation of endogenous bioactive GA levels and GA-regulated root architecture might also play an important role in regulating plant adaptation to environmental changes.
In several perennial woody plant species, GAs are important regulators that control multiple developmental aspects, including wood formation and stress adaptations. However, in Eucalyptus, the role of GA in regulating the above-mentioned processes largely remains elusive. In addition, there is still a lack of systemic identification and characterization of the key genes involved in GA synthesis, degradation, signaling, and of investigation of their molecular functions in the regulation of diverse aspects of plant growth and development. Based on the third-generation transcriptome sequencing strategy by PacBio Iso-Seq, the genome-wide expressed genes in the main vegetative organs of E. grandis × E. urophylla were identified. Based on this, we comprehensively identified the key genes involved in GA synthesis, degradation, and early signaling components and constructed GA-related networks. In addition, we also investigated the expression pattern of the main GA synthesis and degradation genes in response to abiotic stresses (salt and osmotic). To further investigate the molecular functions of the identified EguGA20ox and EguGA2ox genes in E. grandis × E. urophylla, we selectively cloned and overexpressed EguGA20ox1, EguGA20ox2, and EguGA2ox1 in Arabidopsis and Eucalyptus and analyzed their role in controlling plant growth, root initiation and growth, xylem development, and abiotic stress responses. This study not only systematically identified the key gene families involved in the regulation of GA synthesis, degradation, and signaling for the first time, based on third-generation transcriptome sequencing, but also provides evidence that GA is stress-responsive and might directly regulate stress adaptation via regulation of GA homeostasis in Eucalyptus. Therefore, this study aimed to identify GA-related gene families and their expression patterns in response to environmental factors, and to characterize the functions of key GA metabolism-related genes; this identification and characterization could greatly improve our understanding of the role of GA in the regulation of plant growth and stress-adaptation in Eucalyptus, and facilitate the selection of high-yield and stress-resistant cultivars though genetic modification strategies.

2. Results

2.1. Transcriptome Sequencing of Five Major Vegetative Organs of E. grandis × E. urophylla

The vegetative organs of Eucalyptus, such as wood and fibers, as well as cineole, are the primary resources for industrial purposes. To investigate the global gene expression in the vegetative organs of E. grandis × E. urophylla, tissues from the stem, root, mature leaf, young leaf, and shoot apex were collected for transcriptome sequencing. The third-generation transcriptome sequencing platform PacBio Iso-Seq was used for the analysis. A total of 71.59 Gb (approximately 100-fold the E. grandis genome size) of subreads were obtained, and a total of 26 M subreads were collected. After error correction, a total of 59,948 consensus transcripts were obtained, with a mean length of 2716 bp, N50 of 2883, and N90 of 1914 (Figure 1A). The read length range is shown in Figure 1A, in which the most abundant read length was from 2000 to 4000 bp (Figure 1B). As the genomes of E. grandis × E. urophylla and E. urophylla are lacking, the consensus transcripts were mapped to the genome of one of the hybrid’s parents, E. grandis. The results showed that 91.47% of the transcripts could be mapped to the E. grandis genome (sequence similarity >97%), 0.78% were multiple-mapped, and 7.75% were unmapped (Figure 1A). The genome of E. grandis × E. urophylla included both parental genomes; thus, we postulated that 7.75% of unmapped transcripts could be specifically transcribed from the E. urophylla genome. Nevertheless, the possibility that E. grandis and E. urophylla might share high sequence similarity of many genes in each other’s genomes cannot be ruled out. Moreover, the results showed that the density of the transcript length of E. grandis × E. urophylla was higher than that of E. grandis (Figure 1C). A total of 91.47% of the transcripts were finely mapped to the eleven E. grandis chromosomes, and some transcripts were mapped to the mitochondrial genome (Figure 1D).
Typically, with increasing chromosome size, the number of expressed genes would increase accordingly. Interestingly, NC_052619.1 and NC_052618.1 were the largest chromosomes, but the number of mapped reads of these two chromosomes were significantly lower than that of NC_052617.1. NC_052621.1 was the smallest chromosome, but the number of mapped transcripts was no less than that of NC_052614.1 (Figure 1D,E). These results suggested that the gene intensity might not always be in accordance with chromosome size. However, we cannot rule out the possibility that some high-similarity genes from the E. urophylla chromosomes were also mapped to their homologs in the E. grandis chromones, which might have caused the changed gene density at the chromosome level. Gene annotation of all the mapped and unmapped transcripts was then carried out using seven gene annotation databases: NR, NT, Pfam, KOG/COG, SwissProt, KEGG, and GO. A total of 18,957 mapped transcripts and 4643 unmapped transcripts were annotated (Figure 1F); this facilitated further characterization of the gene functions of E. grandis × E. urophylla.

2.2. Identification of Key Gene Families in Each Step of GA Metabolism and Signaling

The E. grandis × E. urophylla hybrid exhibits excellent growth rate, wood quality, and abiotic stress resistance. As GA is a multifaceted phytohormone that directly or in part regulates the abovementioned traits [10,11,12,43], investigation of the role of GA in Eucalyptus might help to clarify the molecular mechanism governing the formation of these economic traits and further facilitate genetic modifications. Although the genes of the main steps in GA biosynthesis, degradation, and signaling have been well studied in Arabidopsis and rice in recent decades [10], there is still a lack of integrative investigation of the whole pathway based on genomic and transcriptomic data for most species, including Eucalyptus. Here, using the protein sequences of Arabidopsis CPS, KS, KO, KAO, GA20ox, GA3ox, GA2ox, GID, and DELLAs as the blast query for BLAST analysis, together with conserved domain analysis using TBtools [44], we identified their potential homologs in the transcriptome of E. grandis × E. urophylla. Then, the synthesis, degradation, and signaling routes of Arabidopsis, rice, Populus trichocarpa, and E. grandis × E. urophylla were constructed and compared (Figure 2). It can be clearly seen that the number of family members in some GA-related gene families was conserved, including CPS, KS, and KAO (Figure 2). In this Eucalyptus hybrid, a total of 14 GA receptor GID genes were identified, whereas there was only one in rice, three in Arabidopsis, and four in Populus (Figure 2). Among all the families, GA2ox was one of the largest gene families in these species, with 15 members in Populus, 13 in E. grandis × E. urophylla, 10 in rice, and 9 in Arabidopsis (Figure 2). Moreover, phylogenetic analysis was carried out to investigate the phylogenetic relationships of the GA2ox and GA20ox families among these four plant species (Figure 3). The results showed that in E. grandis × E. urophylla, the 13 members of EguGA2oxs can be divided into two major types, C19-GA2ox and C20-GA2ox (Figure 3A), which is similar to that in peach [26]. Moreover, the phylogenetic tree also showed that C19-GA2ox can be further divided into two clades, C19-GA2ox-I and C19-GA2ox-II (Figure 3A). The results also showed that, with the exception EguGA20ox2, the other EguGA20oxs showed lower genetic homology with the GA20oxs from the other three species (Figure 3B). These phylogenetic trees suggested the diverse evolutionary formation processes of different GA2ox and GA20ox members and indicated the functional differences among different members.

2.3. Expression Profile of the Genes Involved in the Regulation of GA Metabolism in Different Vegetative Organs

Different members of a family might have similar functions in enzymatic activities or signal transduction; however, they might be temporally or spatially regulated in plants to achieve more precise and flexible regulation, as with the GA2ox family genes in peach and Liriodendron chinense [26,27]. In Eucalyptus, each step of GA metabolism is regulated by multiple gene members. To further investigate the expression pattern of these genes, Real-time quantitative PCR (RT qPCR) was used to survey the expression pattern of all the identified GA metabolism-related genes in five major vegetative organs of Eucalyptus (root, stem, mature leaf, young leaf, and shoot apex). The results showed that several early GA synthesis genes, such as EguCPS, EguKS1, EguKS2 and EguKO1, were highly expressed in the root, while EguKAO1 and EguKAO2 were mainly expressed in the stem and mature leaf in the E. grandis × E. urophylla seedlings (Figure 4). EguKAO3 showed relatively high expression in both mature leaves and roots (Figure 4). GA20ox is the key rate-limiting enzyme catalyzing the synthesis of GA intermediates. The results showed that four EguGA20oxs only showed high expression levels in the leaves, with EguGA20ox1 and EguGA20ox3 in the mature leaves and EguGA20ox2 and EguGA20ox4 in the young leaves. Thirteen EguGA2oxs also showed diverse expression patterns in the five different organs. EguGA2ox1, 4, and 8 were found to be highly expressed in the root, EguGA2ox2 in the young leaf, EguGA2ox3, 12, and 13 in the mature leaf, EguGA2ox5 and 6 in the stem, and EguGA2ox9 and 11 in the shoot apex (Figure 4).

2.4. Diverse Expression Patterns of GA Metabolism-Related Genes in Response to Abiotic Stress

Gene expression patterns are closely related to biological functions. We further investigated the expression of GA metabolism genes in response to abiotic stresses (salt and osmotic stress) in three different organs (leaf, stem, and root) of E. grandis × E. urophylla seedlings at 24 h after treatment with 200 mM NaCl and 300 mM mannitol, respectively. The results showed that the expression of early GA synthesis genes, including EguCPS, EguKO, EguKS, and EguKAO, was mostly upregulated in one or more tissues in response to either stress treatment, except that EguKO2 was downregulated in the leaves and roots after mannitol treatment (Figure 5). Four EguGA20ox genes exhibited diverse responses to both stresses, whereas the expression of EguGA3ox1 remained unaffected or was inhibited in response to the stresses (Figure 5). GA2ox genes play a pivotal role in counterbalancing bioactive GAs. Most of the EguGA2ox genes were found to be upregulated in response to the stresses (Figure 5), suggesting that GA might be negatively regulated in response to the stresses in Eucalyptus.

2.5. Overexpression of EguGA20ox Promoted Plant Growth and Biomass Accumulation in Arabidopsis

The expression of several GA metabolism-related genes was diversely regulated in response to abiotic stress treatment. Specifically, several EguGA2oxs, including EguGA2ox1, 3, 4, 7, and 11, were upregulated, whereas EguGA20ox2, 3, and 4 were downregulated after stress treatment (Figure 5). This stress-responsive expression pattern of the GA-metabolism genes in Eucalyptus is similar to that in Liriodendron chinense, where the expression of several GA2oxs was significantly up-regulated in response to different abiotic stresses [27]. Thus, we postulated that GA might be a negative regulator in controlling stress adaptations. To further investigate the molecular functions of the genes involved in GA synthesis (EguGA20ox) or catabolism (EguGA2ox), we selectively cloned and overexpressed EguGA20ox1, EguGA20ox2, and EguGA2ox1 in Arabidopsis. For each gene, over ten independent transgenic Arabidopsis lines with similar phenotypes were obtained. The gene expression of the transgenic Arabidopsis lines was analyzed by RT-qPCR, as shown in Figure S1. The results showed that overexpression of either of the EguGA20oxs significantly increased plant height and caused an early flowering phenotype, while EguGA2ox1-overexpressing (OE) lines exhibited a delayed flowering phenotype (Figure 6A,B), suggesting a classic role of GA in the promotion of flowering in Arabidopsis. Surprisingly, we also detected that overexpression of either EguGA20ox1 or EguGA20ox2 led to significantly enhanced vegetative growth, which resulted in enlarged leaf size and increased biomass accumulation (Figure 6C–G). The GA signal is primarily transduced via the degradation of DELLA proteins. In this study, we discovered that the della mutant did not exhibit similar phenotypes to the EguGA20ox1- or EguGA20ox2-OE lines, except for the early flowering phenotype (Figure 6A). Due to early flowering, the final number of rosette leaves of the EguGA20ox lines was lower than that of the Col-0 plant (Figure S2). Nevertheless, before the plant started bolting, the number of leaves of either EguGA20ox line was significantly greater than that of Col-0 or della mutant (Figure S2).

2.6. Overexpression of EguGA20ox and EguGA2ox in Arabidopsis Led to Diverse Responses to Salt Stress

Salt stress treatment was used to investigate the abiotic stress tolerance of EguGA20ox1, EguGA20ox2, and EguGA2ox1 transgenic lines. Seven days after subjecting the roots to treatment with 200 mM NaCl, the growth status of EguGA20ox1, EguGA20ox2, EguGA2ox1, and Col-0 plants was recorded. The results showed that both the EguGA20ox1- and EguGA20ox2-OE lines were hypersensitive to salt stress and had a more severe leaf decoloring and withering phenotype, while the EguGA2ox1-OE line exhibited higher stress tolerance one week after treatment (Figure 6C,D). Due to stress hypersensitivity, the biomass and leaf size of the EguGA20ox-OE lines were more likely to be reduced in response to salt stress, while plant growth seemed to be unaffected in EguGA2ox1 plants (Figure 6E–G).
The photosynthesis parameters were analyzed to further clarify whether overexpression of EguGA20ox or EguGA2ox could alter the key photosynthetic parameters, which were also regarded as important markers for stress tolerance. The results showed that salt stress could impair the photosynthetic ability of EguGA20ox-OE plants; this was reflected by the decreased chlorophyll fluorescence in some leaves (Figure 7A). Intriguingly, the quantum efficiency of photosystem II (Fv/Fm) of EguGA20ox was not significantly affected but was only slightly lower than that of EguGA2ox1 (Figure 7B). However, it can be clearly seen that the maximum electron transport rate (ETRmax) of EguGA20ox1-OE, EguGA20ox1-OE, and della plants was approximately two-fold lower than that of Col-0 plants, while EguGA2ox1-OE exhibited the highest ETRmax. These photosynthetic parameters reflected the higher tolerance to salt stress of the EguGA2ox1-OE lines. Moreover, we also found that the EguGA2ox1 plants exhibited much higher chlorophyll (chla and chlb) and anthocyanin levels (Figure 7D,E). As anthocyanins are the key antioxidants [45], the significantly higher anthocyanin level in EguGA2ox1 plants greatly facilitates adaptations to severe environmental stresses (Figure 6).

2.7. GA20ox Promotes Root Initiation and Lateral Root Growth in Eucalyptus

In addition to the significant promotive role of EguGA20ox in shoot growth, we also detected that the expression of EguGA20ox and EguGA2ox affected root development in Arabidopsis. Previous reports showed that GA could directly affect lateral root growth but exhibited diverse effects in different species [38,39]. Here, our results showed that the number of lateral roots of the EguGA20ox1- and EguGA20ox2-OE Arabidopsis lines was significantly increased, while the EguGA2ox1-OE line exhibited the lowest number of lateral roots among all tested Arabidopsis lines (Figure 8A,B). GA signal transduction was mostly dependent on the degradation of DELLA proteins. The results also revealed that the della mutant also showed an increase in the number of lateral roots. These results suggested that GA might be positively correlated with lateral root development in Arabidopsis.
The gene expression profile of the GA biosynthesis and metabolism genes in different organs of Eucalyptus showed that several early GA synthesis genes (EguCPS, two EguKSs, and EguKO1) and GA degradation genes (EguGA2ox1, EguGA2ox4, and EguGA2ox8) exhibited high expression levels in the root (Figure 4), suggesting that GA might be important for root growth and development in Eucalyptus. Here, we used Agrobacterium rhizogenes-mediated Eucalyptus transformation to investigate whether overexpression of EguGA20ox had any effects on root development in E. grandis seedlings. In comparison with two control groups that were separately inoculated with A. rhizogenes K599 or K599 harboring the 35S::GFP plasmid, the seedlings inoculated with EguGA20ox1 (K599) showed significantly accelerated hairy root induction at one week after inoculation at the seedling hypocotyl (Figure 8C); it took almost one month for the control groups to start root formation. Within 3 to 4 days after inoculation with EguGA20ox1 (K599), obvious callus formation can be clearly seen at the inoculation site (Stage 1). One week later, some inoculated seedlings showed obvious lateral root growth (Stage 2). Frequently, more than one lateral root could be induced at the wound site by EguGA20ox1 (indicated in Figure 8C, Stage 3). Two weeks after inoculation, the hairy root induction ratio reached 25%, and a more significant increase (to 47%) was observed at four weeks (Figure 8E). These results suggested that overexpression of EguGA20ox1 conferred fast A. rhizogenes-induced hairy root formation.
For both EguGA20ox1 and mock-inoculated seedlings, the plants were then transplanted into liquid medium one week after hairy root induction for further analysis of root growth and development. Two weeks after transplantation, the hairy roots overexpressing EguGA20ox1 showed a significantly increased number of lateral roots (Figure 8F), suggesting the potential pivotal role of GA in the regulation of lateral root growth in Eucalyptus. These results suggest that GA might play a pivotal role in the regulation of root growth and development in Eucalyptus. Nevertheless, the molecular mechanisms of the GA-regulated signaling pathway governing root growth in Eucalyptus require further investigation.

2.8. Overexpression of GA20ox1 Significantly Promoted Xylem Development in Eucalyptus

In Arabidopsis and some other perennial woody plants, such as Populus, GA is one of the most important regulators involved in xylem development [10]. Nevertheless, the role of GA in controlling woody tissue formation has scarcely been studied in Eucalyptus. After two weeks in liquid medium, the hairy roots of mock- and EguGA20ox1-OE plants, as shown in Figure 8D, were used for histochemical microscopic examination. The results showed that overexpression of EguGA20ox1 significantly improved xylem differentiation in the transgenic Eucalyptus hairy roots (Figure 9A); this is consistent with previous results in another tree species hybrid (aspen), which showed that increased GA levels could lead to enhanced xylem differentiation and cell elongation [46]. As a comparison, we investigated the xylem development of the hypocotyl of the Arabidopsis lines separately overexpressing EguGA20ox1, EguGA20ox2, and EguGA2ox1 (Figure 9B). The results showed, similar to those in Eucalyptus (Figure 9C), that the EguGA20ox1- and EguGA20ox2-OE Arabidopsis lines exhibited a high level of xylem differentiation; this was reflected by the significantly increased xylem-to-stele ratio (Figure 9D). The della mutant also showed significantly improved xylem development, whereas the EguGA2ox1-OE line showed no difference from Col-0 (Figure 9D). Conclusively, these results suggested that overexpression of EguGA20ox1 could improve xylem differentiation in Eucalyptus.

3. Discussion

Gibberellins play a pivotal role in regulating diverse aspects of plant growth, development, and responses to environmental stimuli, including seed germination, hypocotyl and stem elongation, wood formation, flowering, and stress adaptations [47,48]. The homeostasis of GA signaling is strictly regulated by various steps in GA biosynthesis, degradation, or signaling, most of which are controlled by a group of gene families, (GA20ox, GA3ox, and especially GA2ox) [43]. The temporal and spatial expression of these family genes in different plant organs, and in response to diverse environmental stimuli, lead to a more precise manipulation of GA homeostasis; this is very important for maintaining the adequate balance between fast growth and stress adaptions in plants. In Eucalyptus, one of the fastest growing woody plants, the roles of GA in regulating the abovementioned biological processes are poorly understood. Thus, identification and functional characterization of GA synthesis, degradation, and signaling pathways could further facilitate the molecular breeding of Eucalyptus. In this study, based on the third-generation transcriptome sequencing strategy PacBio Iso-Seq, we first identified the potential genes involved in GA synthesis, degradation, and signaling in the Eucalyptus hybrid E. grandis × E. urophylla, investigated the expression pattern of GA biosynthesis and degradation genes in response to abiotic stresses, and functionally characterized a few genes related to GA synthesis and degradation in Arabidopsis.
Gibberellin is one of the major signals that both regulates vegetative growth and controls the vegetative to reproductive transition in plants [43]. It positively regulates flowering in Arabidopsis [49], though it inhibits flowering in several woody plants, such as Jatropha curcas [50] and Prunus persica [51]. Overexpression of either EguGA20ox1 or EguGA20ox2 significantly shortened the flowering time, while EguGA2ox1 delayed flowering in Arabidopsis (Figure 6), suggesting a conserved role of these genes from Eucalyptus in the regulation of flowering. Countless studies have demonstrated that early flowering is frequently accompanied by decreased vegetative growth. Early flowering caused a significant reduction in biomass accumulation by overexpressing CO and FT in Arabidopsis [52,53]. Our results showed that although the EguGA20ox1- and EguGA20ox2-OE lines both exhibited early flowering in Arabidopsis, they also promoted vegetative growth, resulting in significantly larger rosette leaves and increased numbers of leaves and biomass accumulation under normal growth conditions (Figure 6). The promotive effects of overexpressing GA20ox in vegetative growth were also reported in a few studies in tobacco and Arabidopsis [54,55,56,57]. In contrast, overexpression of EguGA2ox1 caused a significant dwarf phenotype and was accompanied by decreased biomass and a dark-leaf phenotype, which is similar to that in rice [58,59,60]. As a comparison, we also analyzed the phenotype of the quintuple DELLA Arabidopsis mutant (della). The della mutant exhibited a significantly early flowering phenotype but had no growth-promoting phenotypes, as seen in the EguGA20ox-OE plants (Figure 6). In trees, high doses of GA might be toxic and can cause abnormal growth and cell death, as found in our previous study that showed different doses of GA in controlling shoot branching in J. curcas [61]. The diverse vegetative growth status of GA20ox OE and della mutants suggested that GA-regulated plant growth might be dosage-dependent, and that a moderate increase in bioactive GA or GA signaling in plants could improve vegetative growth.
Plants have developed elaborate strategies by spatially and temporally regulating the endogenous hormonal homeostasis to facilitate plant growth or adaptation to complex environments [27]. Phytohormones, such as abscisic acid (ABA), ethylene, jasmonic acid (JA), and salicylic acid (SA), are crucial factors in controlling plant adaptations to environmental stresses [62,63]. In recent years, an increasing number of studies have demonstrated that GA synthesis and degradation are tightly regulated in response to abiotic stresses [27,35,64,65,66,67]. In this study, more than half of the GA metabolism-related genes in Eucalyptus were significantly affected by salt or drought treatment (Figure 5). In Arabidopsis, stresses positively regulate the expression of GA degradation-related genes, while downregulating GA synthesis-related genes [21,37]. Moreover, overexpression of GA2ox confers abiotic stress tolerance in Arabidopsis, rice, and switchgrass [35,65,68,69]. However, the expression of GA20ox2 can be induced by flooding stress in rice; this could promote stem elongation and enhance the flooding stress resistance [64]. In Eucalyptus, our results not only showed that the expression pattern of the key GA metabolism-related genes is closely correlated with stress responses, but further evidence in Arabidopsis showed that the transgenic EguGA2ox1-OE plants exhibited enhanced tolerance to abiotic stress (Figure 5 and Figure 6). On the contrary, EguGA20ox-OE plants were hypersensitive to salt stress and exhibited a significant decrease of the biomass accumulation by salt treatment. The significantly increased anthocyanin content in the leaf as a result of overexpressing EguGA2ox might contribute to the enhanced abiotic stress tolerance (Figure 7E). The GA-DELLA signaling pathway is central in regulating the anthocyanin synthesis in plants [70,71]. Thus, GA depletion by overexpressing GA degradation genes might enhance DELLA-mediated anthocyanin production, which not only promotes the abiotic stress tolerance, but also prevents chlorophyll degradation (Figure 7D) and the maintenance of photosynthetic capacity (Figure 7C), as previously reported [45]. Conclusively, our results agree with those in found in Arabidopsis and rice [37,65], which might suggest a common negative role of GAs in the regulation of abiotic stress adaptation. However, the molecular mechanisms of GA in the regulation of abiotic stress adaptations need further investigation in Eucalyptus.
Root morphogenesis and architecture are mainly controlled by auxin, cytokinin, and ethylene [40]. In recent years, more discoveries have shown that strigolactones also directly control root morphogenesis in response to environmental stimuli, mainly phosphate and nitrogen levels [30]. However, the role of GA in the regulation of root architecture is poorly understood. Here, we discovered that overexpression of EguGA20ox promoted the lateral root growth in both Arabidopsis and Eucalyptus (Figure 8). In Populus, overexpression of GA2ox1 seems to increase the lateral root (LR) number and LR length [38,55], suggesting a negative role of GA in the regulation of root growth. Nevertheless, a recent report examining tobacco and Melia azedarach showed that ectopic expression of the GA20ox gene driven by a root preferential promoter significantly promoted root growth and biomass accumulation [55]. Here, we discovered that in both Arabidopsis and Eucalyptus, overexpression of EguGA20ox1 could increase LR growth (Figure 8). These results might suggest the diverse roles of GA in regulating LR growth and development across different species. Intriguingly, we also discovered that in Eucalyptus seedlings inoculated with K599 (EguGA20ox1), transgenic hairy roots can be efficiently induced within days; hairy root induction typically takes a much longer time (at least three-four weeks), suggesting that GA might be a positive regulator in controlling root formation and development. Moreover, many GA metabolism-related genes exhibited high expression levels in Eucalyptus roots, especially early synthesis genes (Figure 4), which suggests that Eucalyptus root could be a potential place for the synthesis of GA or GA intermediates. Most previous studies of GA have focused on plant aerial organ growth and development, such as the elongation of the hypocotyl, stem, and leaves, and the regulation of seed germination and flowering. As we showed in this study, GA also directly participates in root growth in Eucalyptus, which should be highlighted in future research on perennial woody plants.
The role of GA and its interactions with other phytohormones in the regulation of xylem development have received considerable attention in recent years. However, wood formation in perennial woody plants seems to be more complex and lacks sufficient molecular and genetic evidence due to the lack of a genetic transformation strategy in most woody plants. In our study, we first obtained transgenic Eucalyptus roots overexpressing EguGA20ox1 based on a modified transformation strategy, as reported previously [72], and investigated its role in the regulation of Eucalyptus root growth, including xylem development. The results showed a conserved role of GA synthesis genes in the regulation of xylem development, as reported in Populus and Arabidopsis [73,74]. Importantly, this system in Eucalyptus might be a convenient strategy for further investigation of the role of GA in the regulation of diverse aspects of root development and may facilitate the high-throughput selection and functional characterization of candidate genes in Eucalyptus.

4. Materials and Methods

4.1. Plant Materials and Growth Conditions

E. grandis × E. urophylla hybrid “Guanglin 9” clonal seedlings (three months after transplanting in the soil) were provided by the Guangxi Academy of Forestry (Nanning, China). E. grandis seeds were provided by the Guangxi Dongmen Forestry Farm (Chongzuo, China). The hybrid clonal seedlings and E. grandis seedlings were all planted in a growth chamber (25 °C, 100 μmol·m−2·s−1, 16-h light/8-h dark photoperiod) and fertilized every week with 1/2 Hoagland solution. Wild-type Arabidopsis Col-0 was used for genetic transformation. Arabidopsis plants were planted in a growth chamber (22 °C, 16-h light/8-h dark photoperiod). The Arabidopsis della quinary mutant (gai; rga; rgl1; rgl2; rgl3, CS16298) was ordered from the Nottingham Arabidopsis Stock Centre.

4.2. RNA Extraction, cDNA Synthesis and Quantitative Real-Time PCR

Total RNA of Arabidopsis and Eucalyptus was extracted using the Plant Total RNA Extraction Kit (Omega, Shanghai, China). The quantity and quality of RNA was verified by using gel electrophoresis and a Nanodrop 2000c ultraviolet spectrophotometer (Eppendorf, Framingham, MA, USA). The purified RNA was reverse transcribed into cDNA using the HiScript III 1st Strand cDNA Synthesis Kit (+gDNA wiper) (Vazyme, Nanjing, China). The quality of the cDNA was verified by amplification of the reference genes ACT2 and PP2A from Arabidopsis and Eucalyptus using PCR. The relative gene expression level was quantified by qPCR using SYBR Master Mix (Vazyme, Nanjing, China) according to the manufacturer’s protocol on a LightCycler 96 system (Roche, Basel, Switzerland). PCR conditions were set as follows: 95 °C for 5 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. After amplification, a melting curve was obtained to confirm the presence of a simple amplicon. The relative gene expression level to the reference genes Arabidopsis ACT2 or Eucalyptus PP2A-1 and PP2A-3 [75] was calculated based on the 2−ΔΔCt method. The qPCR experiments were conducted in triplicate. Information on the primers is listed in Table S1.

4.3. Gene Cloning, Overexpression Vector Construction and Arabidopsis Transformation

The full-length sequences of EguGA20ox1, EguGA20ox2, and EguGA2ox1 were amplified from the cDNA of E. grandis × E. urophylla using primers containing restriction sites and inserted into the pOCA30 expression vector after sequencing. The overexpression vectors EguGA20ox1, EguGA20ox2, and EguGA2ox1 were introduced into the Agrobacterium tumefaciens strain EHA105 and used for Arabidopsis transformation. Four-week-old Arabidopsis thaliana Columbia (Col-0) plants were transformed using A. tumefaciens and the floral dip method as previously reported [76]. More than ten independent transformants for each gene were obtained. The T2 homozygous Arabidopsis transgenic lines were used for further physiological and biochemical analysis.

4.4. Eucalyptus Transformation

Eucalyptus transformation to induce transgenic hairy roots was performed as previously described [72,77]. Briefly, E. grandis seeds were germinated in sterilized peat soil. One-week-old Eucalyptus seedlings were used for transformation. The seedling hypocotyls were infected with a needle swabbed with Agrobacterium rhizogenes (K599, OD600 = 2.0) by stabbing four holes along the hypocotyl. The infected seedlings were then transferred to a plastic box (20 cm × 10 cm × 5 cm) containing sterilized peat soil that was saturated with 1/2 Hoagland solution before use. Then, the plastic box was covered with a preservative film to maintain high humidity and kept under dim light conditions. One week after inoculation, the induction efficiency and the root phenotype were recorded. To calculate the transformation efficiency, at least three independent transformation experiments were performed.

4.5. Sample Collection, Iso-Seq Library Construction and Transcriptome Sequencing

Five different plant tissues from the three-month-old E. grandis × E. urophylla clonal seedlings were collected and frozen at −80 °C. The total RNA of the mixture of the five tissues was isolated following the aforementioned method. The quality of the RNA was further verified using Qubit and Agilent 2100. The Iso-Seq library was prepared following the isoform sequencing protocol of the Clontech SMARTer PCR cDNA Synthesis Kit and BluePippin Size Selection System, as described by Pacific Biosciences. Sequence data were processed using SMRTlink 5.0 software. Additional nucleotide errors in consensus reads were corrected using the Illumina RNA-seq data with LoRDEC v0.9 software. The consensus reads were then aligned to the reference using GMAP with the default parameters. The functional annotation of the unmapped transcripts and novel gene transcripts was based on the following databases: NR (NCBI nonredundant protein sequences), NT (NCBI nonredundant nucleotide sequences), Pfam (Protein family), KOG/COG (Clusters of Orthologous Groups of proteins), Swiss-Prot (a manually annotated and reviewed protein sequence database), KO (KEGG Ortholog database), and GO (Gene Ontology).

4.6. Sequence Analysis of GA Synthesis, Degradation and Signaling Genes

To identify the GA synthesis, degradation, and signaling genes in the E. grandis × E. urophylla, A. thaliana, Oryza sativa, and Populus trichocarpa genome, BLASTP with default parameters of the data files of different species was performed based on the conserved amino acid sequences using the amino acid sequences of Arabidopsis GA synthesis, degradation, and signaling-related genes as queries on the TBtools v1.115 software [44]. The search results were aligned by the ClustalW algorithm with default parameters to reduce duplicate and redundant sequences. The amino acid sequences of all GA-related proteins of Arabidopsis, rice, and Populus used in this study were downloaded from the NCBI database. The phylogenetic trees were constructed using the neighbor-joining method in MEGA (version 5.1). Unless specifically mentioned, the default parameters were used for all bioinformatic analyses. All the detailed information of the used sequences are listed as the Table S2.

4.7. Salt and Osmotic Stress Treatment of Eucalyptus and Arabidopsis

Three-month-old E. grandis × E. urophylla clonal seedlings with similar growth status were used separately for salt and drought treatments using 200 mM NaCl and 300 mM mannitol (to mimic drought stress), respectively. The growth pots were directly submerged in different solutions until the soil was completely saturated. After that, the plants were kept under normal growth conditions. Then, different tissues of the seedlings were collected for gene expression analysis. For the analysis of the stress response of different Arabidopsis lines, the same approaches were carried out on the three-week-old Arabidopsis seedlings. The survival rate, growth status, and other physiological parameters were recorded after stress treatment. For all the treatments, the experiments were conducted in at least three replicates.

4.8. Analysis of the Photosynthesis Parameters of Arabidopsis Seedlings after Stress Treatment

Light responses were monitored by simultaneously obtaining the chlorophyll fluorescence and other photosynthetic parameters in the leaves by using the Imaging-PAM M-Series Chlorophyll Fluorometer (Heinz Walz, Effeltrich, Germany). To investigate whether abiotic stresses could have diverse effects on photosynthetic activities, three-week-old transgenic Arabidopsis seedlings (EguGA20ox1-, EguGA20ox2-, and EguGA2ox1-OE lines, Col-0, and della mutant) subjected to salt stress for one week were used for analysis. Three mature leaves per plant (>3 plants per line) were selected for the parameter detection.

4.9. Paraffin Section and Microscopic Observation

Hypocotyl and root cross-sections of Arabidopsis and Eucalyptus plants were collected and histologically analyzed. Briefly, the plant samples were fixed in FAA solution and then dehydrated in an increasing concentration of ethanol from 30% to 100%. After hyalinization in xylene, the samples were embedded in paraffin. Then, the samples were cut into slices using a microtome. After dewaxing, the samples were double stained with 1% sarranine and 0.5% fast green. Images of stem sections were used to measure the xylem diameter using ImageJ v1.54d software.

4.10. Statistical Analysis

Multiple comparisons between different groups were carried out using R (version 4.2.2). The data were analyzed using one-way ANOVA and Duncan’s test, with significance set at p < 0.05. Student’s t test was also used to analyze the significant differences between the indicated groups and the control. The heatmaps were constructed using the pheatmap package (1.0.12).

5. Conclusions

In this study, the regulatory gene network of gibberellin biosynthesis, degradation, and signaling was comprehensively analyzed, based on the transcriptome of five major vegetative organs of E. grandis × E.urophylla. Our results showed that the genes involved in GA synthesis and degradation were diversely regulated in different vegetative organs and in response to the biotic stresses in Eucalyptus. the expression of most GA2ox genes was positively regulated in response to stress treatment. Importantly, the functional characterization of GA synthesis genes EguGA20ox1 and EguGA20ox2, and GA degradation gene EguGA2ox1 in Eucalyptus and Arabidopsis, demonstrated that they participated in the regulation of plant growth, root initiation and growth, abiotic responses, and xylem development. Our study provides the first comprehensive understanding of regulatory pathways of GA metabolism and signaling, and demonstrated the key role of GA in regulating the balance between plant growth and stress responses. This could further benefit molecular breeding for obtaining high-yield and stress-resistant Eucalyptus cultivars.

Supplementary Materials

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

Author Contributions

Conceptualization, J.N.; Investigation, W.W., J.N., L.Z., P.W., K.L., Y.L., H.H. and L.D.; Data analysis, W.W., J.N., L.L., J.X., H.H. and X.C.; Supervision, J.N.; Writing and editing, J.N., W.W. and Z.-F.X. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Guangxi Natural Science Foundation (2021GXNSFBA220062 and AD22035175), and National Natural Science Foundation of China (32001345), Guangxi Key Research and Development Program (JB22035001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data of the transcriptome is available in Genbank with the accession numberPRJNA944269.

Acknowledgments

We thank Kaiqin Ye for giving critical advice in improving the manuscript. We thank Jianzhong Wang and Bingfa Qiu for their help on the experiments!

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gominho, J.; Figueira, J.; Rodrigues, J.C.; Pereira, H. Within-tree variation of heartwood, extractives and wood density in the eucalypt hybrid urograndis (Eucalyptus grandis × E. urophylla). Wood Fiber SCI 2001, 33, 3–8. [Google Scholar]
  2. Gomes da Cunha, T.Q.; Santos, A.C.; Novaes, E.; Santiago Hansted, A.L.; Yamaji, F.M.; Sette, C.R., Jr. Eucalyptus expansion in Brazil: Energy yield in new forest frontiers. Biomass Bioenerg. 2021, 144, 105900. [Google Scholar] [CrossRef]
  3. Hakamada, R.E.; Hubbard, R.M.; Moreira, G.G.; Stape, J.L.; Campoe, O.; de Barros Ferraz, S.F. Influence of stand density on growth and water use efficiency in Eucalyptus clones. Forest Ecol. Manag. 2020, 466, 118125. [Google Scholar] [CrossRef]
  4. Myburg, A.A.; Grattapaglia, D.; Tuskan, G.A.; Hellsten, U.; Hayes, R.D.; Grimwood, J.; Jenkins, J.; Lindquist, E.; Tice, H.; Bauer, D.; et al. The genome of Eucalyptus grandis. Nature 2014, 510, 356–362. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Vining, K.J.; Romanel, E.; Jones, R.C.; Klocko, A.; Alves-Ferreira, M.; Hefer, C.A.; Amarasinghe, V.; Dharmawardhana, P.; Naithani, S.; Ranik, M.; et al. The floral transcriptome of Eucalyptus grandis. New Phytol. 2015, 206, 1406–1422. [Google Scholar] [CrossRef] [PubMed]
  6. Ayala, P.G.; Acevedo, R.M.; Luna, C.V.; Rivarola, M.; Acuna, C.; Poltri, S.M.; Gonzalez, A.M.; Sansberro, P.A. Transcriptome dynamics of rooting zone and leaves during In vitro adventitious root formation in Eucalyptus nitens. Plants 2022, 11, 3301. [Google Scholar] [CrossRef] [PubMed]
  7. Shang, X.; Zhang, P.; Liu, G.; Zhan, N.; Wu, Z. Comparative transcriptomics analysis of contrasting varieties of Eucalyptus camaldulensis reveals wind resistance genes. PEERJ 2022, 10, e12954. [Google Scholar] [CrossRef]
  8. Oates, C.N.; Carsten, K.; Myburg, A.A.; Bernard, S.; Sanushka, N. The transcriptome and terpene profile of Eucalyptus grandis reveals mechanisms of defense against the insect pest, Leptocybe invasa. Plant Cell Physiol. 2015, 7, 1418–1428. [Google Scholar] [CrossRef] [Green Version]
  9. Zhang, Y.; Li, J.J.; Li, C.R.; Chen, S.K.; Tang, Q.L.; Xiao, Y.F.; Zhong, L.X.; Chen, Y.Y.; Chen, B. Gene expression programs during callus development in tissue culture of two Eucalyptus species. BMC Plant Biol. 2022, 22, 1. [Google Scholar] [CrossRef]
  10. Hedden, P. The current status of research on gibberellin biosynthesis. Plant Cell Physiol. 2020, 61, 1832–1849. [Google Scholar] [CrossRef]
  11. Gao, S.; Chu, C. Gibberellin metabolism and signaling: Targets for improving agronomic performance of crops. Plant Cell Physiol. 2020, 61, 1902–1911. [Google Scholar] [CrossRef] [PubMed]
  12. Leopoldina Ojeda-Barrios, D.; Raquel Orozco-Melendez, L.; Cano-Medrano, R.; Sanchez-Chavez, E.; Angel Parra-Quezada, R.; Calderon-Jurado, M.; Jacbo-Cuellar, J.; Hernandez-Ordonez, E.; Cruz-Alvarez, O. Non-structural carbohydrates, foliar nutrients, yield components and oxidative metabolism in pecan trees in response to foliar applications of growth regulators. Agriculture 2022, 12, 688. [Google Scholar] [CrossRef]
  13. Hedden, P. The genes of the green revolution. Trends Genet. 2003, 19, 5–9. [Google Scholar] [CrossRef]
  14. Lopez-Cristoffanini, C.; Serrat, X.; Jauregui, O.; Nogues, S.; Lopez-Carbonell, M. Phytohormone profiling method for rice: Effects of GA20ox mutation on the gibberellin content of Japonica rice varieties. Front. Plant Sci. 2019, 10, 733. [Google Scholar] [CrossRef]
  15. Gallego-Giraldo, L.; Ubeda-Tomas, S.; Gisbert, C.; Garcia-Martinez, J.L.; Moritz, T.; Lopez-Diaz, I. Gibberellin homeostasis in tobacco is regulated by gibberellin metabolism genes with different gibberellin sensitivity. Plant Cell Physiol. 2008, 49, 679–690. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Israelsson, M.; Mellerowicz, E.; Chono, M.; Gullberg, J.; Moritz, T. Cloning and overproduction of gibberellin 3-oxidase in hybrid aspen trees. Effects of gibberellin homeostasis and development. Plant Physiol. 2004, 135, 221–230. [Google Scholar] [CrossRef] [Green Version]
  17. Reinecke, D.M.; Wickramarathna, A.D.; Ozga, J.A.; Kurepin, L.V.; Jin, A.L.; Good, A.G.; Pharis, R.P. Gibberellin 3-oxidase gene expression patterns influence gibberellin biosynthesis, growth, and development in pea. Plant Physiol. 2013, 163, 929–945. [Google Scholar] [CrossRef] [Green Version]
  18. Oikawa, T.; Koshioka, M.; Kojima, K.; Yoshida, H.; Kawata, M. A role of OsGA20ox1, encoding an isoform of gibberellin 20-oxidase, for regulation of plant stature in rice. Plant Mol. Biol. 2004, 55, 687–700. [Google Scholar] [CrossRef]
  19. Fleet, C.M.; Yamaguchi, S.; Hanada, A.; Kawaide, H.; David, C.J.; Kamiya, Y.; Sun, T.P. Overexpression of AtCPS and AtKS in Arabidopsis confers increased ent-kaurene production but no increase in bioactive gibberellins. Plant Physiol. 2003, 132, 830–839. [Google Scholar] [CrossRef] [Green Version]
  20. Hernandez-Garcia, J.; Briones-Moreno, A.; Blazquez, M.A. Origin and evolution of gibberellin signaling and metabolism in plants. Semin. Cell Dev. Biol. 2021, 109, 46–54. [Google Scholar] [CrossRef]
  21. Li, C.; Zheng, L.L.; Wang, X.N.; Hu, Z.B.; Zheng, Y.; Chen, Q.H.; Hao, X.C.; Xiao, X.; Wang, X.B.; Wang, G.D.; et al. Comprehensive expression analysis of Arabidopsis GA2-oxidase genes and their functional insights. Plant Sci. 2019, 285, 1–13. [Google Scholar] [CrossRef]
  22. Lange, M.J.P.; Liebrandt, A.; Arnold, L.; Chmielewska, S.-M.; Felsberger, A.; Freier, E.; Heuer, M.; Zur, D.; Lange, T. Functional characterization of gibberellin oxidases from cucumber, Cucumis sativus L. Phytochem 2013, 90, 62–69. [Google Scholar] [CrossRef]
  23. Sakamoto, T.; Miura, K.; Itoh, H.; Tatsumi, T.; Ueguchi-Tanaka, M.; Ishiyama, K.; Kobayashi, M.; Agrawal, G.K.; Takeda, S.; Abe, K.; et al. An overview of gibberellin metabolism enzyme genes and their related mutants in rice. Plant Physiol. 2004, 134, 1642–1653. [Google Scholar] [CrossRef] [Green Version]
  24. Gou, J.Q.; Ma, C.; Kadmiel, M.; Gai, Y.; Strauss, S.; Jiang, X.N.; Busov, V. Tissue-specific expression of Populus C-19 GA 2-oxidases differentially regulate above- and below-ground biomass growth through control of bioactive GA concentrations. New Phytol. 2011, 192, 626–639. [Google Scholar] [CrossRef]
  25. Martinez-Bello, L.; Moritz, T.; Lopez-Diaz, I. Silencing C-19-GA 2-oxidases induces parthenocarpic development and inhibits lateral branching in tomato plants. J. Exp. Bot. 2015, 66, 5897–5910. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Cheng, J.; Ma, J.; Zheng, X.; Lv, H.; Zhang, M.; Tan, B.; Ye, X.; Wang, W.; Zhang, L.L.; Li, Z.Q.; et al. Functional analysis of the Gibberellin 2-oxidase gene family in peach. Front. Plant Sci. 2021, 12, 619158. [Google Scholar] [CrossRef]
  27. Hu, L.F.; Wang, P.K.; Hao, Z.D.; Lu, Y.; Xue, G.X.; Cao, Z.J.; Qu, H.X.; Cheng, T.L.; Shi, J.S.; Chen, J. Gibberellin oxidase gene family in L. chinense: Genome-wide identification and gene expression analysis. Int. J. Mol. Sci. 2021, 22, 7167. [Google Scholar] [CrossRef] [PubMed]
  28. Aerts, N.; Mendes, M.P.; Van Wees, S.C.M. Multiple levels of crosstalk in hormone networks regulating plant defense. Plant J. 2021, 105, 489–504. [Google Scholar] [CrossRef] [PubMed]
  29. Sun, C.; Liu, L.; Wang, L.; Li, B.; Jin, C.; Lin, X. Melatonin: A master regulator of plant development and stress responses. J. Integ. Plant Biol. 2021, 63, 126–145. [Google Scholar] [CrossRef] [PubMed]
  30. Yu, Z.; Duan, X.; Luo, L.; Dai, S.; Ding, Z.; Xia, G. How plant hormones mediate salt stress responses. Trends Plant Sci. 2020, 25, 1117–1130. [Google Scholar] [CrossRef] [PubMed]
  31. Zhang, H.; Zhu, J.; Gong, Z.; Zhu, J.-K. Abiotic stress responses in plants. Nat. Rev. Genet. 2022, 23, 104–119. [Google Scholar] [CrossRef] [PubMed]
  32. Cardinale, F.; Krukowski, P.K.; Schubert, A.; Visentin, I. Strigolactones: Mediators of osmotic stress responses with a potential for agrochemical manipulation of crop resilience. J. Exp. Bot. 2018, 69, 2291–2303. [Google Scholar] [CrossRef] [PubMed]
  33. Wang, Q.; Smith, S.M.; Huang, J. Origins of strigolactone and karrikin signaling in plants. Trends Plant Sci. 2022, 27, 450–459. [Google Scholar] [CrossRef]
  34. Colebrook, E.H.; Thomas, S.G.; Phillips, A.L.; Hedden, P. The role of gibberellin signalling in plant responses to abiotic stress. J. Exp. Biol. 2014, 217, 67–75. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Magome, H.; Yamaguchi, S.; Hanada, A.; Kamiya, Y.; Oda, K. The DDF1 transcriptional activator upregulates expression of a gibberellin-deactivating gene, GA2ox7, under high-salinity stress in Arabidopsis. Plant 2008, 56, 613–626. [Google Scholar] [CrossRef] [PubMed]
  36. He, H.H.; Liang, G.P.; Lu, S.X.; Wang, P.P.; Liu, T.; Ma, Z.H.; Zuo, C.W.; Sun, X.M.; Chen, B.H.; Mao, J. Genome-wide identification and expression analysis of GA2ox, GA3ox, and GA20ox are related to gibberellin oxidase genes in grape (Vitis vinifera L.). Genes 2019, 10, 680. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Lange, M.J.P.; Lange, T. Touch-induced changes in Arabidopsis morphology dependent on gibberellin breakdown. Nat. Plants 2015, 1, 14025. [Google Scholar] [CrossRef]
  38. Gou, J.Q.; Strauss, S.H.; Tsai, C.J.; Fang, K.; Chen, Y.R.; Jiang, X.N.; Busov, V.B. Gibberellins regulate lateral root formation in Populus through interactions with auxin and other hormones. Plant Cell 2010, 22, 623–639. [Google Scholar] [CrossRef] [Green Version]
  39. Sun, H.; Pang, B.; Yan, J.; Wang, T.; Wang, L.; Chen, C.H.; Li, Q.; Ren, Z.H. Comprehensive analysis of cucumber gibberellin oxidase family genes and functional characterization of CsGA20ox1 in root development in Arabidopsis. Int. J. Mol. Sci. 2018, 19, 3135. [Google Scholar] [CrossRef] [Green Version]
  40. Karlova, R.; Boer, D.; Hayes, S.; Testerink, C. Root plasticity under abiotic stress. Plant Physiol. 2021, 187, 1057–1070. [Google Scholar] [CrossRef]
  41. Ptoskova, K.; Szecowka, M.; Jaworek, P.; Tarkowska, D.; Petrik, I.; Pavlovic, I.; Novak, O.; Thomas, S.G.; Phillips, A.L.; Hedden, P. Changes in the concentrations and transcripts for gibberellins and other hormones in a growing leaf and roots of wheat seedlings in response to water restriction. BMC Plant Biol. 2022, 22, 284. [Google Scholar] [CrossRef]
  42. Lange, T.; Kramer, C.; Pimenta Lange, M.J. The Class III Gibberellin 2-Oxidases AtGA2ox9 and AtGA2ox10 contribute to cold stress tolerance and fertility. Plant Physiol. 2020, 184, 478–486. [Google Scholar] [CrossRef]
  43. Claeys, H.; De Bodt, S.; Inze, D. Gibberellins and DELLAs: Central nodes in growth regulatory networks. Trends Plant Sci. 2014, 19, 231–239. [Google Scholar] [CrossRef]
  44. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  45. Xu, Z.; Rothstein, S.J. ROS-Induced anthocyanin production provides feedback protection by scavenging ROS and maintaining photosynthetic capacity in Arabidopsis. Plant Signal. Behav. 2018, 13, e1451708. [Google Scholar] [CrossRef] [Green Version]
  46. Eriksson, M.E.; Israelsson, M.; Olsson, O.; Moritz, T. Increased gibberellin biosynthesis in transgenic trees promotes growth, biomass production and xylem fiber length. Nat. Biotechnol. 2000, 18, 784–788. [Google Scholar] [CrossRef] [PubMed]
  47. Yamauchi, Y.; Takeda-Kamiya, N.; Hanada, A.; Ogawa, M.; Kuwahara, A.; Seo, M.; Kamiya, Y.; Yamaguchi, S. Contribution of gibberellin deactivation by AtGA2ox2 to the suppression of germination of dark-imbibed Arabidopsis thaliana seeds. Plant Cell Physiol. 2007, 48, 555–561. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Kijidani, Y.; Tsuyama, T.; Takata, K. Seasonal variations of auxin and gibberellin A4 levels in cambial-region tissues of three conifers (Pinus elliottii, Chamaecyparis obtusa, and Cryptomeria japonica) with inherently different wood densities. J. Wood Sci. 2021, 67, 46. [Google Scholar] [CrossRef]
  49. Bao, S.; Hua, C.; Shen, L.; Yu, H. New insights into gibberellin signaling in regulating flowering in Arabidopsis. J. Integr. Plant Biol. 2020, 62, 118–131. [Google Scholar] [CrossRef] [Green Version]
  50. Li, J.L.; Pan, B.-Z.; Niu, L.J.; Chen, M.-S.; Tang, M.; Xu, Z.-F. Gibberellin inhibits floral initiation in the perennial woody plant Jatropha curcas. J. Plant Growth Regul. 2018, 37, 999–1006. [Google Scholar] [CrossRef] [Green Version]
  51. Garcia-Pallas, I.; Val, J.; Blanco, A. The inhibition of flower bud differentiation in ‘Crimson Gold’ nectarine with GA(3) as an alternative to hand thinning. Sci. Hortic. 2001, 90, 265–278. [Google Scholar] [CrossRef] [Green Version]
  52. Yoo, S.K.; Chung, K.S.; Kim, J.; Lee, J.H.; Hong, S.M.; Yoo, S.J.; Lee, J.S.; Ahn, J.H. CONSTANS activates SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 through FLOWERING LOCUS T to promote flowering in Arabidopsis. Plant Physiol. 2005, 139, 770–778. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Teper-Bamnolker, P.; Samach, A. The flowering integrator FT regulates SEPALLATA3 and FRUITFULL accumulation in Arabidopsis leaves. Plant Cell 2005, 17, 2661–2675. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Garcia-Hurtado, N.; Carrera, E.; Ruiz-Rivero, O.; Pilar Lopez-Gresa, M.; Hedden, P.; Gong, F.; Luis Garcia-Martinez, J. The characterization of transgenic tomato overexpressing gibberellin 20-oxidase reveals induction of parthenocarpic fruit growth, higher yield, and alteration of the gibberellin biosynthetic pathway. J. Exp. Bot. 2012, 63, 5803–5813. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Ly, L.K.; Bui, T.P.; Le, A.V.T.; Nguyen, P.V.; Ong, P.X.; Pham, N.B.; Zhang, Z.J.; Do, P.T.; Chu, H. Enhancing plant growth and biomass production by overexpression of GA20ox gene under control of a root preferential promoter. Transg. Res. 2022, 31, 73–85. [Google Scholar] [CrossRef]
  56. Rieu, I.; Ruiz-Rivero, O.; Fernandez-Garcia, N.; Griffiths, J.; Powers, S.J.; Gong, F.; Linhartova, T.; Eriksson, S.; Nilsson, O.; Thomas, S.G.; et al. The gibberellin biosynthetic genes AtGA20ox1 and AtGA20ox2 act, partially redundantly, to promote growth and development throughout the Arabidopsis life cycle. Plant J. 2008, 53, 488–504. [Google Scholar] [CrossRef]
  57. Wang, J.; Wu, B.; Li, J.; Yin, H.; Fan, Z.; Li, X.; He, L. Overexpression and silent expression of CrGA20ox1 from Camellia reticulata ‘Hentiangao’ and its effect on morphological alterations in transgenic tobacco plants. Plant Breeding 2018, 137, 903–911. [Google Scholar] [CrossRef]
  58. Huang, J.; Tang, D.; Shen, Y.; Qin, B.; Hong, L.; You, A.; Li, M.; Wang, X.; Yu, H.K.; Gu, M.H.; et al. Activation of gibberellin 2-oxidase 6 decreases active gibberellin levels and creates a dominant semi-dwarf phenotype in rice (Oryza sativa L.). J. Genet. Genomics 2010, 37, 23–36. [Google Scholar] [CrossRef]
  59. Lo, S.-F.; Yang, S.-Y.; Chen, K.-T.; Hsing, Y.-l.; Zeevaart, J.A.D.; Chen, L.-J.; Yu, S.-M. A novel class of Gibberellin 2-Oxidases control semidwarfism, tillering, and root development in rice. Plant Cell 2008, 20, 2603–2618. [Google Scholar] [CrossRef] [Green Version]
  60. Schomburg, F.M.; Bizzell, C.M.; Lee, D.J.; Zeevaart, J.A.D.; Amasino, R.M. Overexpression of a novel class of gibberellin 2-oxidases decreases gibberellin levels and creates dwarf plants. Plant Cell 2003, 15, 151–163. [Google Scholar] [CrossRef] [Green Version]
  61. Ni, J.; Gao, C.; Chen, M.-S.; Pan, B.-Z.; Ye, K.; Xu, Z.-F. Gibberellin promotes shoot branching in the perennial woody plant Jatropha curcas. Plant Cell Physiol. 2015, 56, 1655–1666. [Google Scholar] [CrossRef] [Green Version]
  62. Verma, V.; Ravindran, P.; Kumar, P.P. Plant hormone-mediated regulation of stress responses. BMC Plant Biol. 2016, 16, 86. [Google Scholar] [CrossRef] [Green Version]
  63. Ku, Y.S.; Sintaha, M.; Cheung, M.Y.; Lam, H.M. Plant hormone signaling crosstalks between biotic and abiotic stress responses. Int. J. Mol. Sci. 2018, 19, 3206. [Google Scholar] [CrossRef] [Green Version]
  64. Kuroha, T.; Nagai, K.; Gamuyao, R.; Wang, D.R.; Furuta, T.; Nakamori, M.; Kitaoka, T.; Adachi, K.; Minami, A.; Mori, Y.; et al. Ethylene-gibberellin signaling underlies adaptation of rice to periodic flooding. Science 2018, 361, 181–186. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Lo, S.F.; Ho, T.H.D.; Liu, Y.L.; Jiang, M.J.; Hsieh, K.T.; Chen, K.T.; Yu, L.C.; Lee, M.H.; Chen, C.Y.; Huang, T.P.; et al. Ectopic expression of specific GA2 oxidase mutants promotes yield and stress tolerance in rice. Plant Biotechnol. J. 2017, 15, 850–864. [Google Scholar] [CrossRef] [Green Version]
  66. Achard, P.; Gong, F.; Cheminant, S.; Alioua, M.; Hedden, P.; Genschik, P. The cold-inducible CBF1 factor-dependent signaling pathway modulates the accumulation of the growth-repressing DELLA proteins via its effect on gibberellin metabolism. Plant Cell 2008, 20, 2117–2129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Zhang, Y.; Liu, Z.; Liu, R.; Hao, H.; Bi, Y. Gibberellins negatively regulate low temperature-induced anthocyanin accumulation in a HY5/HYH-dependent manner. Plant Signal. Behav. 2011, 6, 632–634. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Wuddineh, W.A.; Mazarei, M.; Zhang, J.; Poovaiah, C.R.; Mann, D.G.J.; Ziebell, A.; Sykes, R.W.; Davis, M.F.; Udvardi, M.K.; Stewart, C.N. Identification and overexpression of gibberellin 2-oxidase (GA2ox) in switchgrass (Panicum virgatum L.) for improved plant architecture and reduced biomass recalcitrance. Plant Biotechnol. J. 2015, 13, 636–647. [Google Scholar] [CrossRef] [Green Version]
  69. Shan, C.; Mei, Z.; Duan, J.; Chen, H.; Feng, H.; Cai, W. OsGA2ox5, a gibberellin metabolism enzyme, is involved in plant growth, the root gravity response and salt stress. PLoS ONE 2014, 9, e87110. [Google Scholar] [CrossRef] [Green Version]
  70. Xie, Y.; Tan, H.; Ma, Z.; Huang, J. DELLA proteins promote anthocyanin biosynthesis via sequestering MYBL2 and JAZ suppressors of the MYB/bHLH/WD40 complex in Arabidopsis thaliana. Mol. Plant 2016, 9, 711–721. [Google Scholar] [CrossRef] [Green Version]
  71. Zhang, Y.; Liu, Z.; Liu, J.; Lin, S.; Wang, J.; Lin, W.; Xu, W. GA-DELLA pathway is involved in regulation of nitrogen deficiency-induced anthocyanin accumulation. Plant Cell Rep. 2017, 36, 557–569. [Google Scholar] [CrossRef]
  72. Plasencia, A.; Soler, M.; Dupas, A.; Ladouce, N.; Silva-Martins, G.; Martinez, Y.; Lapierre, C.; Franche, C.; Truchet, I.; Grima-Pettenati, J. Eucalyptus hairy roots, a fast, efficient and versatile tool to explore function and expression of genes involved in wood formation. Plant Biotechnol. J. 2016, 14, 1381–1393. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Park, E.J.; Kim, H.T.; Choi, Y.I.; Lee, C.; Nguyen, V.P.; Jeon, H.W.; Cho, J.S.; Funada, R.; Pharis, R.P.; Kurepin, L.V.; et al. Overexpression of gibberellin 20-oxidase1 from Pinus densiflora results in enhanced wood formation with gelatinous fiber development in a transgenic hybrid poplar. Tree Physiol. 2015, 35, 1264–1277. [Google Scholar] [PubMed] [Green Version]
  74. Jeon, H.-W.; Cho, J.-S.; Park, E.-J.; Han, K.-H.; Choi, Y.-I.; Ko, J.-H. Developing xylem-preferential expression of PdGA20ox1, a gibberellin 20-oxidase 1 from Pinus densiflora, improves woody biomass production in a hybrid poplar. Plant Biotechnol. J. 2016, 14, 1161–1170. [Google Scholar] [CrossRef] [PubMed]
  75. Cassan-Wang, H.; Soler, M.; Yu, H.; Camargo, E.L.O.; Carocha, V.; Ladouce, N.; Savei, B.; Paiva, J.A.P.; Leple, J.C.; Grima-Pettenati, J. Reference genes for high-throughput quantitative reverse transcription-PCR analysis of gene expression in organs and tissues of Eucalyptus grown in various environmental conditions. Plant Cell Physiol. 2012, 53, 2101–2116. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Bent, A. Arabidopsis thaliana floral dip transformation method. Methods Mol. Biol. 2006, 343, 87–103. [Google Scholar]
  77. Soler, M.; Plasencia, A.; Larbat, R.; Pouzet, C.; Jauneau, A.; Rivas, S.; Pesquet, E.; Lapierre, C.; Truchet, I.; Grima-Pettenati, J. The Eucalyptus linker histone variant EgH1.3 cooperates with the transcription factor EgMYB1 to control lignin biosynthesis during wood formation. New Phytol. 2017, 213, 287–299. [Google Scholar] [CrossRef] [PubMed]
Ijms 24 07051 g001 550
Figure 1. Overview of the transcriptome sequencing of the five major vegetative organs of Eucalyptus grandis × E. urophylla. (A) Statistics of GMAP mapping. (B) FLNC length distribution of the transcripts. (C) Comparison of the transcript length distribution against the reference genome of E. grandis. (D) Number of the mapped reads in different chromosomes of E. grandis. (E) Reads density in different chromosomes. (F) Gene annotation in different databases.
Figure 1. Overview of the transcriptome sequencing of the five major vegetative organs of Eucalyptus grandis × E. urophylla. (A) Statistics of GMAP mapping. (B) FLNC length distribution of the transcripts. (C) Comparison of the transcript length distribution against the reference genome of E. grandis. (D) Number of the mapped reads in different chromosomes of E. grandis. (E) Reads density in different chromosomes. (F) Gene annotation in different databases.
Ijms 24 07051 g001
Ijms 24 07051 g002 550
Figure 2. Overview of the gene families in the synthesis, degradation, and signaling of gibberellin. The number of the family members in E. grandis × E. urophylla, Populus, Oryza sativa, and Arabidopsis was showed in the heatmaps.
Figure 2. Overview of the gene families in the synthesis, degradation, and signaling of gibberellin. The number of the family members in E. grandis × E. urophylla, Populus, Oryza sativa, and Arabidopsis was showed in the heatmaps.
Ijms 24 07051 g002
Ijms 24 07051 g003 550
Figure 3. Phylogenetic analysis of GA2ox (A) and GA20 (B) families of different plant species (E. grandis × E. urophylla, Arabidopsis thaliana, Oryza sativa, and Populus trichocarpa). The sequences were aligned using the Mega 5.0 software. A maximum likelihood tree was calculated based on the alignment and evaluated using 1000 bootstrap repetitions.
Figure 3. Phylogenetic analysis of GA2ox (A) and GA20 (B) families of different plant species (E. grandis × E. urophylla, Arabidopsis thaliana, Oryza sativa, and Populus trichocarpa). The sequences were aligned using the Mega 5.0 software. A maximum likelihood tree was calculated based on the alignment and evaluated using 1000 bootstrap repetitions.
Ijms 24 07051 g003
Ijms 24 07051 g004 550
Figure 4. Gene expression profile of the GA synthesis and degradation genes in five different vegetative organs of E. grandis × E. urophylla. The expression level of each gene was normalized to the reference genes PP2A-1 and PP2A-3. The heatmaps were plotted using data obtained by RT-qPCR with three biological replicates. SA, shoot apex; YL, young leaf; ML, mature leaf; S, stem; R, root.
Figure 4. Gene expression profile of the GA synthesis and degradation genes in five different vegetative organs of E. grandis × E. urophylla. The expression level of each gene was normalized to the reference genes PP2A-1 and PP2A-3. The heatmaps were plotted using data obtained by RT-qPCR with three biological replicates. SA, shoot apex; YL, young leaf; ML, mature leaf; S, stem; R, root.
Ijms 24 07051 g004
Ijms 24 07051 g005 550
Figure 5. Gene expression profile of the GA synthesis and degradation genes in the leaves, stems, and roots of the three-month-old E. grandis × E. urophylla seedlings at 24 h after being subjected to salt (200 mM NaCl) or mannitol (300 mM) treatment. The expression level of each gene was normalized to the reference genes PP2A-1 and PP2A-3. The heatmaps were plotted using data obtained by RT-qPCR with three biological replicates. L, leaf; S, stem; R, root.
Figure 5. Gene expression profile of the GA synthesis and degradation genes in the leaves, stems, and roots of the three-month-old E. grandis × E. urophylla seedlings at 24 h after being subjected to salt (200 mM NaCl) or mannitol (300 mM) treatment. The expression level of each gene was normalized to the reference genes PP2A-1 and PP2A-3. The heatmaps were plotted using data obtained by RT-qPCR with three biological replicates. L, leaf; S, stem; R, root.
Ijms 24 07051 g005
Ijms 24 07051 g006 550
Figure 6. Overexpression of either EguGA20ox1 or EguGA20ox2 promoted vegetative growth but decreased salt stress tolerance in Arabidopsis. (A) Plant growth status of different Arabidopsis lines at 45 d. (B) Plant height at 45 d (n = 20). (C) Plant status of different Arabidopsis lines before (upper layer), and one week after salt stress treatment (under layer). (D) Plant growth status two weeks after stress treatment. Fresh weight (E), leaf length (F), and width (G) of the mature rosette leaves of different Arabidopsis lines were calculated at one week after mock or stress treatment (n > 10). Data represents as mean ± sd. In each group, different lowercase letters represent the degree of significant difference (p < 0.05). Bars = 1 cm.
Figure 6. Overexpression of either EguGA20ox1 or EguGA20ox2 promoted vegetative growth but decreased salt stress tolerance in Arabidopsis. (A) Plant growth status of different Arabidopsis lines at 45 d. (B) Plant height at 45 d (n = 20). (C) Plant status of different Arabidopsis lines before (upper layer), and one week after salt stress treatment (under layer). (D) Plant growth status two weeks after stress treatment. Fresh weight (E), leaf length (F), and width (G) of the mature rosette leaves of different Arabidopsis lines were calculated at one week after mock or stress treatment (n > 10). Data represents as mean ± sd. In each group, different lowercase letters represent the degree of significant difference (p < 0.05). Bars = 1 cm.
Ijms 24 07051 g006
Ijms 24 07051 g007 550
Figure 7. EguGA2ox1 overexpression increased the maximum electron transport ratio and the content of chlorophyll and anthocyanin in the leaf. (A) Picture of the photosynthetic inflorescence. Fv/Fm (B) and ETRmax (C) of different Arabidopsis lines one week after salt stress treatment (n = 4). (D,E) Chlorophyll and anthocyanin content of three-week-old Arabidopsis seedlings (n = 3). Data represents as mean ± sd. Different lowercase letters represent the degree of significant difference (p < 0.05). Bars = 1 cm.
Figure 7. EguGA2ox1 overexpression increased the maximum electron transport ratio and the content of chlorophyll and anthocyanin in the leaf. (A) Picture of the photosynthetic inflorescence. Fv/Fm (B) and ETRmax (C) of different Arabidopsis lines one week after salt stress treatment (n = 4). (D,E) Chlorophyll and anthocyanin content of three-week-old Arabidopsis seedlings (n = 3). Data represents as mean ± sd. Different lowercase letters represent the degree of significant difference (p < 0.05). Bars = 1 cm.
Ijms 24 07051 g007
Ijms 24 07051 g008 550
Figure 8. Overexpression of EguGA20ox in Arabidopsis and Eucalyptus improved the root growth. (A,B) Lateral root number of different Arabidopsis lines was counted two weeks after germination on the MS plates (n = 12). (C) Four stages of the hairy root induction in Eucalyptus seedlings. (D) Root growth status after growing in the liquid culture for two weeks. (E) Hairy root induction ratio after transformation with mock, GFP, or EguGA20ox1 (n = 3, 90 seedlings for each replicate). (F) Lateral root number of the transgenic hairy roots of E. grandis overexpressing EguGA20ox1 (n = 10). Data represents as mean ± sd. **, p < 0.01. Different lowercase letters represent the degree of significant difference (p < 0.05).
Figure 8. Overexpression of EguGA20ox in Arabidopsis and Eucalyptus improved the root growth. (A,B) Lateral root number of different Arabidopsis lines was counted two weeks after germination on the MS plates (n = 12). (C) Four stages of the hairy root induction in Eucalyptus seedlings. (D) Root growth status after growing in the liquid culture for two weeks. (E) Hairy root induction ratio after transformation with mock, GFP, or EguGA20ox1 (n = 3, 90 seedlings for each replicate). (F) Lateral root number of the transgenic hairy roots of E. grandis overexpressing EguGA20ox1 (n = 10). Data represents as mean ± sd. **, p < 0.01. Different lowercase letters represent the degree of significant difference (p < 0.05).
Ijms 24 07051 g008
Ijms 24 07051 g009 550
Figure 9. Overexpression of EguGA20ox1 promotes xylem development in Eucalyptus and Arabidopsis. (A) Cross-section of the hairy root of control and EguGA20ox1 OE lines (upper layer, 50-fold enlarged, under layer, 100-fold enlarged). (B) Cross-section of the hypocotyl of three-week-old Arabidopsis seedlings. Ratio of xylem to stele area of Eucalyptus (C) and Arabidopsis (D) (n > 4). Data represents as mean ± sd. **, p < 0.01. Different lowercase letters represent the degree of significant difference (p < 0.05). Xy, xylem; Ph, phloem; Ct, cortex.
Figure 9. Overexpression of EguGA20ox1 promotes xylem development in Eucalyptus and Arabidopsis. (A) Cross-section of the hairy root of control and EguGA20ox1 OE lines (upper layer, 50-fold enlarged, under layer, 100-fold enlarged). (B) Cross-section of the hypocotyl of three-week-old Arabidopsis seedlings. Ratio of xylem to stele area of Eucalyptus (C) and Arabidopsis (D) (n > 4). Data represents as mean ± sd. **, p < 0.01. Different lowercase letters represent the degree of significant difference (p < 0.05). Xy, xylem; Ph, phloem; Ct, cortex.
Ijms 24 07051 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wu, W.; Zhu, L.; Wang, P.; Liao, Y.; Duan, L.; Lin, K.; Chen, X.; Li, L.; Xu, J.; Hu, H.; et al. Transcriptome-Based Construction of the Gibberellin Metabolism and Signaling Pathways in Eucalyptus grandis × E. urophylla, and Functional Characterization of GA20ox and GA2ox in Regulating Plant Development and Abiotic Stress Adaptations. Int. J. Mol. Sci. 2023, 24, 7051. https://doi.org/10.3390/ijms24087051

AMA Style

Wu W, Zhu L, Wang P, Liao Y, Duan L, Lin K, Chen X, Li L, Xu J, Hu H, et al. Transcriptome-Based Construction of the Gibberellin Metabolism and Signaling Pathways in Eucalyptus grandis × E. urophylla, and Functional Characterization of GA20ox and GA2ox in Regulating Plant Development and Abiotic Stress Adaptations. International Journal of Molecular Sciences. 2023; 24(8):7051. https://doi.org/10.3390/ijms24087051

Chicago/Turabian Style

Wu, Wenfei, Linhui Zhu, Pan Wang, Yuwu Liao, Lanjuan Duan, Kai Lin, Xin Chen, Lijie Li, Jiajing Xu, Hao Hu, and et al. 2023. "Transcriptome-Based Construction of the Gibberellin Metabolism and Signaling Pathways in Eucalyptus grandis × E. urophylla, and Functional Characterization of GA20ox and GA2ox in Regulating Plant Development and Abiotic Stress Adaptations" International Journal of Molecular Sciences 24, no. 8: 7051. https://doi.org/10.3390/ijms24087051

APA Style

Wu, W., Zhu, L., Wang, P., Liao, Y., Duan, L., Lin, K., Chen, X., Li, L., Xu, J., Hu, H., Xu, Z. -F., & Ni, J. (2023). Transcriptome-Based Construction of the Gibberellin Metabolism and Signaling Pathways in Eucalyptus grandis × E. urophylla, and Functional Characterization of GA20ox and GA2ox in Regulating Plant Development and Abiotic Stress Adaptations. International Journal of Molecular Sciences, 24(8), 7051. https://doi.org/10.3390/ijms24087051

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

Article Metrics

Back to TopTop