Hormonal Regulation and Stimulation Response of Jatropha curcas L. Homolog Overexpression on Tobacco Leaf Growth by Transcriptome Analysis

The Flowering locus T (FT) gene encodes the florigen protein, which primarily regulates the flowering time in plants. Recent studies have shown that FT genes also significantly affect plant growth and development. The FT gene overexpression in plants promotes flowering and suppresses leaf and stem development. This study aimed to conduct a transcriptome analysis to investigate the multiple effects of Jatropha curcas L. homolog (JcFT) overexpression on leaf growth in tobacco plants. The findings revealed that JcFT overexpression affected various biological processes during leaf development, including plant hormone levels and signal transduction, lipid oxidation metabolism, terpenoid metabolism, and the jasmonic-acid-mediated signaling pathway. These results suggested that the effects of FT overexpression in plants were complex and multifaceted, and the combination of these factors might contribute to a reduction in the leaf size. This study comprehensively analyzed the effects of JcFT on leaf development at the transcriptome level and provided new insights into the function of FT and its homologous genes.


Introduction
The flowering locus T (FT) gene encodes a unique protein known as florigen, a protein hormone that plays a vital role in promoting flowering in plants [1]. FT genes are expressed in the phloem and can be transported to the shoot apices, where they promote the initiation of flowering [2]. The FT protein combines with the basic leucine zipper (bZIP) transcription factor flowering locus D (FD) to form FD/FT heterodimer complexes [3]. These complexes activate the suppression of co overexpression 1 (SOC1) and apetala 1 (AP1) proteins, which play a key role in triggering the transition from vegetative to reproductive growth [4]. Moreover, studies indicated that FT, FD, and 14-3-3 proteins form a flowering activation complex in rice that induces the transcription of downstream flowering-related genes, ultimately leading to the initiation of flowering [5]. So far, the molecular function of FT regulating flowering has been studied thoroughly [6,7].
FT-like proteins have been found to possess other important molecular functions besides their well-known role in regulating flowering. These functions include regulating the growth of shoot meristems and leaves, controlling the formation of callose plugs, modulating the ability of signals to reach the apical meristem, and regulating the transcriptional activity of genes controlling cell division [1,[8][9][10]. The overexpression of FT can also regulate H + -ATPase activity in guard cells and induce stomatal opening [11]. Florigen in

Overview of RNA-Seq Data and Biological Processes Affected by JcFT Overexpression
During the transcriptome analysis of the leaves of three JcFT OE samples and three control samples, 39.90 GB of clean data were obtained (National Center for Biotechnology Information (NCBI) Sequence Read Archive accession number: PRJNA943360). Each sample produced clean data with a size of 6.20 GB, and the Q30 base quality score was greater than 94.15%. Of the clean reads, 93.69-94.35% mapped to the Nicotiana tabacum "TN90"

Overview of RNA-Seq Data and Biological Processes Affected by JcFT Overexpression
During the transcriptome analysis of the leaves of three JcFT OE samples and three control samples, 39.90 GB of clean data were obtained (National Center for Biotechnology Information (NCBI) Sequence Read Archive accession number: PRJNA943360). Each sample produced clean data with a size of 6.20 GB, and the Q30 base quality score was greater than 94.15%. Of the clean reads, 93.69-94.35% mapped to the Nicotiana tabacum "TN90" genome, with 92. .12% mapped to exons, 3.02-4.41% mapped to introns, and 2.86-3.41% mapped to intergenic regions. The sequencing output data of each sample are presented in Table S1. The Pearson correlation coefficients (r) ranged from 0.996 to 0.997 between the three JcFT OE samples and from 0.989 to 0.998 in the three control samples. The r values of JcFT OE and control samples were relatively lower (0.958-0.991). The principal component analysis (PCA) results of the fragments per kilobase of exon model per million mapped fragments (FPKM) of each sample showed that the three JcFT OE samples and three control samples were each clustered together ( Figure S1), suggesting that sampling and sequencing data produced consistent results.
Upon JcFT overexpression, 1119 differentially expressed genes (DEGs) [log2 fold change (FC) ≥ 2 and false discovery rate (FDR) < 0.01] were identified in JcFT OE plants, with 729 upregulated and 390 downregulated ( Figure S2). The hierarchical clustering analysis of DEGs demonstrated that the expression pattern of related genes in both the three control samples and the three JcFT OE samples could be accurately distinguished (Figure 2a), indicating that these DEGs well reflected the effect of JcFT overexpression. The classification of the Eukaryotic Orthologous Groups (KOG) showed that the DEGs were involved in 25 biological processes ( Figure S3). The Gene Ontology (GO) enrichment analysis revealed that JcFT overexpression primarily influenced various processes. These processes included positive regulation of transfer from RNA polymerase II promoter, lipid oxidation, oxylipin biosynthetic process, regulation of defense response, response to wounding, and regulation of JA-mediated signaling pathway, among which genes involved in the last three processes were all upregulated ( Figure 2b). Additionally, according to gene set enrichment analysis (GSEA), genes involved in the positive regulation of transfer from RNA polymerase II promoter and lipid oxidation tended to be upregulated (Figure 2c,d).
With respect to the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses for DEGs, 11 pathways exhibited significant effects. These pathways included plant hormone signal transduction, flavonoid biosynthesis, zeatin biosynthesis, phenylpropanoid biosynthesis, fatty acid elongation, and alpha-linolenic acid metabolism, among others ( Figure 3a). The KEGG enrichment analysis for upregulated and downregulated DEGs revealed significant upregulation of the plant hormone signal transduction pathway, diterpenoid and monoterpenoid biosynthesis, and linoleic acid metabolism ( Figure 3b). Conversely, fatty acid metabolism, fatty acid elongation, and phenylpropanoid biosynthesis were enriched in downregulated DEGs (Figure 3c). Notably, both upregulated and downregulated DEGs were significantly enriched in zeatin biosynthesis and flavonoid biosynthesis (Figure 3b,c). With respect to the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses for DEGs, 11 pathways exhibited significant effects. These pathways included plant hormone signal transduction, flavonoid biosynthesis, zeatin biosynthesis, phenylpropanoid biosynthesis, fatty acid elongation, and alpha-linolenic acid metabolism, among others ( Figure 3a). The KEGG enrichment analysis for upregulated and downregulated DEGs revealed significant upregulation of the plant hormone signal transduction pathway, diterpenoid and monoterpenoid biosynthesis, and linoleic acid metabolism ( Figure 3b). Conversely, fatty acid metabolism, fatty acid elongation, and phenylpropanoid biosynthesis were enriched in downregulated DEGs (Figure 3c). Notably, both upregulated and downregulated DEGs were significantly enriched in zeatin biosynthesis and flavonoid biosynthesis (Figure 3b,c).
The DEGs were further analyzed using ClueGO in the Cytoscape (v3.8.2) software (showing only pathways with p ≤ 0.05). Twelve functional groups of GO terms or KEGG pathways were identified, and each group was represented by its most significant term (Table 1 and Figure S4), including the regulation of the JA-mediated signaling pathway, oxylipin metabolism, hormone level regulation, isoprenoid biosynthesis, carboxy-lyase activity, carbon-oxygen lyase activity, response to wounding and detoxification, among others. The number of DEGs and percentage of GO terms included in each group are presented in Table 1, and specific GO terms for each group are depicted in Table S2. Among these groups, the regulation of the JA-mediated signaling pathway had the most terms associated with it, with seven terms accounting for 33.33% and 84 DEGs. These identified    The DEGs were further analyzed using ClueGO in the Cytoscape (v3.8.2) software (showing only pathways with p ≤ 0.05). Twelve functional groups of GO terms or KEGG pathways were identified, and each group was represented by its most significant term (Table 1 and Figure S4), including the regulation of the JA-mediated signaling pathway, oxylipin metabolism, hormone level regulation, isoprenoid biosynthesis, carboxy-lyase activity, carbon-oxygen lyase activity, response to wounding and detoxification, among others. The number of DEGs and percentage of GO terms included in each group are presented in Table 1, and specific GO terms for each group are depicted in Table S2. Among these groups, the regulation of the JA-mediated signaling pathway had the most terms associated with it, with seven terms accounting for 33.33% and 84 DEGs. These identified functional groups suggested that JcFT overexpression significantly impacted the hormonerelated pathways in leaves, affecting the hormone levels, particularly the metabolism and response of JA.

Regulation of Cytokinin and Auxin Levels
A group related to the regulation of hormone levels was significantly annotated through the analysis of DEGs using ClueGO (Table 1). Further examination of these DEGs using BiNGO revealed that the regulation of hormone levels was primarily linked to cytokinin and auxin metabolic processes, where oxidation-reduction played a crucial role ( Figure 4a). The impact on cytokinin metabolic processes involved three significantly upregulated cytokinin dehydrogenase (CKX) and two lonely guy (LOG) gene families, while the auxin metabolic process involved three indole-3-acetic acid (IAA) amino acid hydrolase ILR1-like (ILL) gene families, and two of them were upregulated (Table S3). Moreover, the reverse transcription-quantitative polymerase chain reaction (RT-qPCR) analysis confirmed that JcFT overexpression led to changes in the expression of these eight genes, which was consistent with the findings of the transcriptome analysis (Figure 4b,c). CKX, one of the key enzymes regulating CK levels in plants, caused irreversible cytokinin degradation [21], indicating that the cytokinin levels of JcFT OE leaves were likely reduced.
which was consistent with the findings of the transcriptome analysis (Figure 4b,c). CKX, one of the key enzymes regulating CK levels in plants, caused irreversible cytokinin degradation [21], indicating that the cytokinin levels of JcFT OE leaves were likely reduced.

Regulation of Plant Hormone Signal Transduction
The FT protein is a proteinaceous hormone [1]. The KEGG pathway annotation showed that FT gene overexpression affected numerous biosynthesis pathways, metabolic pathways, and plant hormone signal transduction, which are relevant to eight plant hormone signal transduction pathways, including auxin, cytokinin, gibberellin, abscisic acid (ABA), ethylene, brassinosteroid, JA, and salicylic acid (SA). Transcriptome data indicated that all but 3 of the 17 pathways involved had upregulated DEGs associated with them ( Figure 5a). This suggested that the transcription of phytohormone signal transduction factors in JcFT OE plants was enhanced.
Of these pathways, the JA signal transduction pathway was found to be the most significantly affected, with a large number of DEGs involved, including 13 jasmonate zim domain (JAZ) and 4 basic helix-loop-helix (bHLH) transcription factor myelocytomatosis protein 2 (MYC2), all of which were significantly upregulated (Table S4). The expression patterns of these DEGs, as identified in the transcriptome analysis, were further validated by RT-qPCR analysis (Figure 5b), which confirmed that JcFT overexpression significantly impacted the expression of genes related to JA signal transduction.
The signal transduction pathway of auxin involved four types of DEGs, namely three auxin/indole-3-acetic acid (AUX/IAA) genes, one Gretchen Hagen 3 (GH3) gene, seven small auxin up RNA (SAUR) genes, and four auxin response factor (ARF) genes. Except for the

Regulation of Plant Hormone Signal Transduction
The FT protein is a proteinaceous hormone [1]. The KEGG pathway annotation showed that FT gene overexpression affected numerous biosynthesis pathways, metabolic pathways, and plant hormone signal transduction, which are relevant to eight plant hormone signal transduction pathways, including auxin, cytokinin, gibberellin, abscisic acid (ABA), ethylene, brassinosteroid, JA, and salicylic acid (SA). Transcriptome data indicated that all but 3 of the 17 pathways involved had upregulated DEGs associated with them ( Figure 5a). This suggested that the transcription of phytohormone signal transduction factors in JcFT OE plants was enhanced.
Of these pathways, the JA signal transduction pathway was found to be the most significantly affected, with a large number of DEGs involved, including 13 jasmonate zim domain (JAZ) and 4 basic helix-loop-helix (bHLH) transcription factor myelocytomatosis protein 2 (MYC2), all of which were significantly upregulated (Table S4). The expression patterns of these DEGs, as identified in the transcriptome analysis, were further validated by RT-qPCR analysis (Figure 5b), which confirmed that JcFT overexpression significantly impacted the expression of genes related to JA signal transduction.
The signal transduction pathway of auxin involved four types of DEGs, namely three auxin/indole-3-acetic acid (AUX/IAA) genes, one Gretchen Hagen 3 (GH3) gene, seven small auxin up RNA (SAUR) genes, and four auxin response factor (ARF) genes. Except for the downregulation of ARF expression, the expression of other factors was upregulated. RT-qPCR analysis further confirmed the expression of these DEGs, which was consistent with the results of transcriptome analysis (Figure 5b).

Effects of JcFT OE on Lipid Biosynthetic Process
In the ClueGO analysis of DEGs, the process groups of the oxylipin metabolic process and isoprenoid biosynthetic process were significantly enriched (Figure 6a), indicating that JcFT overexpression had a significant impact on the two interrelated processes of lipid biosynthesis in tobacco leaves. The DEGs involved in the oxylipin metabolic process were found to be the same as those in the oxylipin biosynthetic process. These DEGs were related to dioxygenase activity, oxidoreductase activity, and other functional groups. Similarly, DEGs in the enriched isoprenoid biosynthetic process were related to the functional group of carbon-oxygen lyase activity, acting on phosphates.
Oxylipins, which are derived from the oxidation of polyunsaturated fatty acids such as linoleic acid, are enzymatically converted by lipoxygenase (LOX) or α-dioxygenases (α-DOX) [27]. These enzymes play a role in plant signaling responses to both biotic and abiotic stresses [28]. In this study, the DEGs related to oxylipin metabolism included eight LOXs and one α-DOX, all of which were oxygenase genes involved in linoleic acid metabolism (Table S5). The LOXs comprised one LOX6, one LOX1.6, two LOX1.5, two LOX3.1, and two LOX2.1 genes. All but one, LOX2.1, were upregulated (the log2FC value increased from 1.17 to 6.10).
The GSEA suggested that the transcription of DEGs involved in linoleic acid metabolism tended to be upregulated (Figure 6b). Furthermore, RT-qPCR analysis of these upregulated DEGs confirmed the significant changes in their expression (Figure 6c). These results indicated that JcFT overexpression in tobacco leaves significantly increased the expression of genes related to linolenic acid oxidation, which was associated with increased The cytokinin signaling pathway in JcFT OE plants was found to be associated with one cytokinin receptor 1 (CRE1) and two type-A authentic response regulators (ARRs). CRE1 transduces cytokinin signals across the membrane, which is a class of plant hormones that plays an important role in regulating cell division and differentiation [22]. The proteins encoded by ARRs act as negative regulators of cytokinin response, and mutations in these genes have been shown to increase cytokinin sensitivity [23]. Therefore, the leaves in the JcFT OE group were less sensitive to cytokinins than those in the control group.
The signaling pathways related to gibberellin and ABA showed distinct responses in JcFT OE plants. The gibberellin signaling pathway was primarily associated with a downregulated gene, gibberellin-insensitive dwarf 1 (GID1), which is crucial for the positive regulation of gibberellin signals. On the contrary, the DEGs of the ABA signal transduction pathway were associated with two clade A protein phosphatases type-2C, three sucrose non-fermenting 1-related subfamily 2 (SnRK2), and two ABA-responsive element (ABRE) binding factors (ABFs). The SnRK2 protein kinase plays a key role in regulating the transcriptional response to ABA. The upregulation of SnRK2 suggested that the response to ABA was likely to increase in the leaves of JcFT OE plants.
The DEGs involved in the ethylene signal transduction pathway included one upregulated ethylene insensitive 4 (EIN4) and one upregulated constitutive triple response 1 (CTR1). EIN4 is an ethylene receptor that transmits signals to downstream CTR1, a protein kinase that negatively regulates ethylene signaling and suppresses ethylene responses in Arabidopsis [24,25]. The SA signal transduction pathway involves the nonexpresser of pr genes 1, which is a receptor of SA and plays a role in regulating SA-induced defense gene expression [26].
In summary, in plant signal transduction, JcFT OE affected the sensitivity of various phytohormones, including auxin, cytokinin, gibberellin, ABA, ethylene, JA, and SA. This was evidenced by the significant upregulation of DEGs associated with these pathways and the increased expression of related transcription factors. Additionally, the transmembrane transport capacity of cytokinin and brassinosteroid in leaves was also impacted by JcFT OE .

Effects of JcFT OE on Lipid Biosynthetic Process
In the ClueGO analysis of DEGs, the process groups of the oxylipin metabolic process and isoprenoid biosynthetic process were significantly enriched (Figure 6a), indicating that JcFT overexpression had a significant impact on the two interrelated processes of lipid biosynthesis in tobacco leaves. The DEGs involved in the oxylipin metabolic process were found to be the same as those in the oxylipin biosynthetic process. These DEGs were related to dioxygenase activity, oxidoreductase activity, and other functional groups. Similarly, DEGs in the enriched isoprenoid biosynthetic process were related to the functional group of carbon-oxygen lyase activity, acting on phosphates.
Int. J. Mol. Sci. 2023, 24, x FOR PEER REVIEW 9 of 16 oxylipin content. As linolenic acid is a precursor in JA biosynthesis [29], the upregulated LOXs and DOX1 may affect the JA-related pathways (Figure 7). The isoprenoid biosynthetic process involved 18 DEGs, of which 14 were upregulated. Six randomly selected upregulated genes were subjected to RT-qPCR analysis to validate the transcriptional analysis results, which showed consistent results (Figure 6d). These genes were involved in the biosynthesis of monoterpenoid, sesquiterpenoid, triterpenoid, terpenoid backbone, and carotenoid (Table S6). Notably, some genes, such as (−)-alpha terpineol synthase (LOCA107831250), viridiflorene synthase (LOC107769393), and alpha-farnesene synthase (LOC107800232), were significantly upregulated (all log2FC > 6), indicating that JcFT overexpression caused increased expression of genes related to oxylipins and isoprenoid synthesis in tobacco leaves.  Oxylipins, which are derived from the oxidation of polyunsaturated fatty acids such as linoleic acid, are enzymatically converted by lipoxygenase (LOX) or α-dioxygenases (α-DOX) [27]. These enzymes play a role in plant signaling responses to both biotic and abiotic stresses [28]. In this study, the DEGs related to oxylipin metabolism included eight LOXs and one α-DOX, all of which were oxygenase genes involved in linoleic acid metabolism (Table S5). The LOXs comprised one LOX6, one LOX1.6, two LOX1.5, two LOX3.1, and two LOX2.1 genes. All but one, LOX2.1, were upregulated (the log 2 FC value increased from 1.17 to 6.10).
The GSEA suggested that the transcription of DEGs involved in linoleic acid metabolism tended to be upregulated (Figure 6b). Furthermore, RT-qPCR analysis of these upregulated DEGs confirmed the significant changes in their expression (Figure 6c). These results indicated that JcFT overexpression in tobacco leaves significantly increased the expression of genes related to linolenic acid oxidation, which was associated with increased oxylipin content. As linolenic acid is a precursor in JA biosynthesis [29], the upregulated LOXs and DOX1 may affect the JA-related pathways (Figure 7).

Stimulus Response Induced by JcFT Overexpression
In the ClueGO analysis, several significantly enriched and interrelated groups were identified, including the regulation of the JA-mediated signaling pathway, response to wounding, cellular response to chemical stimulus, and detoxification (Table 1 and Figure  S4). The cellular response to the chemical stimulus was further classified into four aspects: cellular response to oxygen-containing compounds, cellular response to organic substances, cellular response to toxic substances, and cellular response to chemical stress (Figure 7, indicated by red arrows). Oxygen-containing compounds included fatty acids and superoxide, while organic substances mainly comprised lipids (primarily fatty acids) and hormones (primarily JAs). Detoxification of toxic substances and chemical stress involves the removal of superoxide radicals and superoxide dismutase activity. Therefore, JcFT overexpression in tobacco leaves resulted in a stimulated response to JA, fatty acids, and superoxide. This was primarily reflected in GO terms related to the JA-mediated signaling pathway, response to wounding, cellular response to chemical stimulus, and detoxification.
The DEGs involved in response to fatty acids were found to be the same as those involved in response to JA, thus indicating a response to JA stimulation (Figure 7, indi- The isoprenoid biosynthetic process involved 18 DEGs, of which 14 were upregulated. Six randomly selected upregulated genes were subjected to RT-qPCR analysis to validate the transcriptional analysis results, which showed consistent results (Figure 6d). These genes were involved in the biosynthesis of monoterpenoid, sesquiterpenoid, triterpenoid, terpenoid backbone, and carotenoid (Table S6). Notably, some genes, such as (−)-alpha terpineol synthase (LOCA107831250), viridiflorene synthase (LOC107769393), and alpha-farnesene synthase (LOC107800232), were significantly upregulated (all log 2 FC > 6), indicating that JcFT overexpression caused increased expression of genes related to oxylipins and isoprenoid synthesis in tobacco leaves.

Stimulus Response Induced by JcFT Overexpression
In the ClueGO analysis, several significantly enriched and interrelated groups were identified, including the regulation of the JA-mediated signaling pathway, response to wounding, cellular response to chemical stimulus, and detoxification (Table 1 and Figure S4). The cellular response to the chemical stimulus was further classified into four aspects: cellular response to oxygen-containing compounds, cellular response to organic substances, cellular response to toxic substances, and cellular response to chemical stress (Figure 7, indicated by red arrows). Oxygen-containing compounds included fatty acids and superoxide, while organic substances mainly comprised lipids (primarily fatty acids) and hormones (primarily JAs). Detoxification of toxic substances and chemical stress involves the removal of superoxide radicals and superoxide dismutase activity. Therefore, JcFT overexpression in tobacco leaves resulted in a stimulated response to JA, fatty acids, and superoxide. This was primarily reflected in GO terms related to the JAmediated signaling pathway, response to wounding, cellular response to chemical stimulus, and detoxification.
The DEGs involved in response to fatty acids were found to be the same as those involved in response to JA, thus indicating a response to JA stimulation (Figure 7, indicated by yellow arrows). Further analysis revealed that JcFT overexpression led to an increase in the expression of LOX-related genes, which caused the formation of oxylipins, an oxygencontaining compound, and eventually triggered a JA reaction. This response to JA involved 11 upregulated DEGs, including 9 TIFY (threonine, isoleucine, phenylalanine, tyrosine) genes (6B, 4B, 10B, and 10A) and 2 JAZ7, which play critical roles in plant growth and stress response as repressors of jasmonate responses [30,31]. Notably, 10 of the 11 JA-related DEGs were also JAZ-related DEGs of the plant hormone signal transduction pathway (Tables S4 and S7). The qPCR analysis further validated the upregulation of the DEGs involved in response to JA in the leaves of JcFT OE plants (Figure 5c). These results suggested that the response of tobacco leaves to JA-like stimulation, induced by JcFT overexpression, was mainly reflected in the signal transduction pathways.

Discussion
The mechanism underlying the FT-protein-induced flowering has been extensively studied [6,32]. Although FT overexpression usually leads to early flowering in plants, few FT genes have been used to produce early maturing varieties for seed harvesting, largely due to the negative impact of FT overexpression on plant growth, including stem thinning and reduced leaf number and size, which ultimately affects seed yield [33]. This study performed a transcriptome analysis of tobacco leaves induced by JcFT overexpression to elucidate the effects of FT on leaf development. The results indicated that JcFT overexpression significantly altered hormone signal transduction, lipid oxidation metabolism, and terpenoid metabolism during leaf development. Moreover, the transcription of genes related to the JA pathway and cell division also underwent changes, ultimately leading to a decrease in the leaf size. This result was consistent with previous studies on FT of Nicotiana tabacum (NtFT) in tobacco and provided further evidence for conserving FT-like genes in various plant species [34].
The leaf size of JcFT OE tobacco was significantly smaller than that of the control ( Figure S5a,b). A comparison of leaf epidermal cells confirmed that JcFT OE tobacco leaves had larger epidermal cells and fewer cell numbers than the control leaves ( Figure S5c,d). The expression of SAURs and CKXs in JcFT OE tobacco leaves was significantly altered (Figures 4b and 5a,b). SAUR proteins regulate protein phosphatases to control H + -ATPase activity and play a central role in auxin-induced acid growth [35]. SAUR overexpression increases the phosphorylation and activity of plasma membrane H + -ATPase, leading to cell wall acidification and promoting cell expansion [36,37]. In this study, a significant upregulation of SAUR expression in JcFT OE tobacco leaves suggested enhanced cell expansion. Kinoshita et al. demonstrated that FT could activate plasma membrane H + -ATPase in guard cells and is believed to be a general growth regulator by modulating H + -ATPase activity [11]. This study speculated that the increased expansion of leaf epidermal cells caused by FT overexpression might be due to the activation of H + -ATPase by increasing the SAUR expression. CKX is one of the key enzymes regulating the CK level in plants, which causes irreversible degradation of cytokinin [21]. The significant upregulation of CKX expression in JcFT OE tobacco leaves indicated a decrease in cytokinin levels in JcFT OE leaves, which might be an essential factor in reducing cell numbers. The JAZ proteins act as repressors of the JA signaling pathway. The transcription factor MYC2 plays a crucial role in the activation of JA response genes and is considered the "Master Switch" of the jasmonate signaling pathway [38,39]. In the absence of JA, MYC2 functions as a transcriptional repressor due to physical interaction with the members of the JAZ family [30]. However, when JA levels are elevated, the rate of JAZ turnover increases, leading to the release of MYC2 from repression. Bioactive JAs promote the interaction between JAZ protein and the coronatine insensitive 1 (COI1) component of the SCFCOI1 ubiquitin E3 ligase, resulting in the ubiquitination and degradation of JAZ by the 26S proteasome [40,41]. As a result, elevated levels of JA trigger MYC2 to release from the repression caused by JAZ interaction, resulting in the activation of JA response genes. The increased transcription of MYC2 and JAZ in JcFT OE plants might enhance the sensitivity of plants to jasmonate signaling.
Some researchers described florigen as a protein hormone [1], so its relationship with and influence on other plant hormones are essential topics. Combining the research on the effect of JcFT overexpression on tobacco stem growth [17] and the results of this study, it was found that JcFT overexpression had an impact on almost all types of hormone signaling pathways in tobacco stems and leaves, such as auxin, cytokinin, gibberellin, ABA, and SA. This suggested that FT and FT-like gene overexpression generally affected the signaling pathways of multiple hormones. Among these, the auxin-related signaling pathway involved more genes in both stems and leaves, most of which were significantly upregulated SAUR genes (7-30 genes). Additionally, the expression of response to JA genes involved in the JA signaling pathway was significantly increased in leaves but not in stems. Moreover, gene expression involving carbohydrates, such as hemicellulose and cellulose, was significantly affected in stems, while gene expression involving the lipid synthesis pathway was significantly affected in leaves. Wu et al. also showed that JcFT overexpression had a comprehensive impact on the tobacco stem at multiple levels of the cell wall, cytoplasm, and nucleus, including the core oscillators involving circular rhythm, cell cycle regulation, initiation of DNA replication, and expression of genes related to xylem development. The initiation of DNA replication and cell cycle regulation might be one of the main factors for the thinning of JcFT-overexpressed tobacco stems [17]. In this study, JcFT overexpression also affected multiple pathways in tobacco leaves but did not exhibit multidimensional effects as in stems. This might be related to the reactions of different tissues or to the period of material collection. This might be related to the reactions of different tissues or the time of material collection, with leaf collection occurring earlier.
In this study, JcFT overexpression in tobacco was achieved through genetic transformation, and the insertion of the JcFT gene into the tobacco genome may also have an impact on transcription, which is not the effect of JcFT. If the used transgenic independent lines would involve the insertion of the same functional gene, the result may reflect the functions of the inserted genes.
FT gene overexpression significantly affects the development of tobacco roots, stems, and leaves, besides their well-known function of promoting flowering [42][43][44][45][46]. Inhibition of leaf growth appears to be a fundamental effect of the FT gene, which has only evolved in seed plants. The FT protein is usually expressed in companion cells and then transported through the phloem to the shoot apical meristem to promote flowering [47]. However, when plants overexpressing FT are grafted onto an ft mutant or a wild-type plant, the ft or wild-type plants bloom earlier without undergoing visible changes such as those seen in FT-overexpressed plants [12,43,48,49]. These observations demonstrated that vascular bundle transport can bypass the effects of FT, resulting in the shrinking of leaves. When FT is naturally expressed, transported FT proteins may appear in cells other than companion cells, phloem, and stem apex cells. Thus, most cells in roots, stems, and leaves avoid the effects of FT, such as inhibiting cell division. As a result, although FT promotes flowering, plants have evolved unique expression and transportation mechanisms to avoid the influence of FT on most plant cells and tissues. Therefore, FT initiates reproduction, but plants prevent its negative effects on vegetative growth.

Vector Construction, Plant Transformation, and Cultivation
Li et al. [20] reported that the JcFT gene (accession no. NM_001308752) was overexpressed in N. tabacum "SR1." The JcFT sequence, isolated from J. curcas, was inserted into the SmaI and SacI sites of the binary vector pBI121 (DNA Cloning Service, Genewiz, Suzhou, China). The resulting vector was designated as CaMV 35S::JcFT, with an empty vector serving as the negative check (CK). The constructed vectors were introduced into the Agrobacterium tumefaciens strain EHA105 and then transformed into the wild-type N. tabacum "SR1" using a leaf disk method [50]. Transgenic plants were identified through PCR amplification of aminoglycoside 3 -phosphotransferase (NPTII) and JcFT genes, using leaf DNA as a template. The primers for NPTII and JcFT amplification are depicted in Table S8. The plants were grown under a 16 h light/8 h dark photoperiod at 26 • C. Transgenic seeds from the T5 generation with randomly selected four independent lines were harvested and used for this study. The JcFT OE plants were identified by morphological characteristics and one-step RT-PCR analysis using leaf total RNA as a template through the HiScript II One-Step RT-PCR Kit (Vazyme, Nanjing, China) ( Figure S6). The length and width of the fourth true leaf were recorded daily from 29 to 41 days of sowing.

RNA Isolation, Library Construction, and Sequencing
JcFT transgenic tobacco and control tobacco were simultaneously sown with 80 seeds each. Subsequently, transgenic plants were distinguished based on morphological characteristics and RT-PCR analysis. Select three out of four JcFT OE independent lines, each corresponding to a duplicate. A total of 20 JcFT transgenic plants, with each line consisting of 6-7 plants and an equal number of control plants, were prepared for further experimentation. On the 33rd and 35th days of sowing, 2-3 plants were randomly selected from each line to collect the fourth true leaf of each plant for transcriptome analysis. Three replicates each of JcFT OE and control were established, consisting of 2-3 mixed and ground leaves for RNA extraction.
Total RNA was extracted using TaKaRa MiniBEST Plant RNA Extraction reagents following the manufacturer's protocol (TaKaRa Bio, Inc., Dalian, China). The quality and concentration of RNA were measured using a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Dover, DE, USA) and an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA), respectively. A NEBNext Poly(A) mRNA Magnetic Isolation Module (E7490; NEB, Rowly, MA, USA) was used to isolate mRNA, which was then fragmented into approximately 200-nt RNA inserts. Next, double-stranded cDNA was synthesized from these fragments, and end-repair/dA-tail and adaptor ligation were performed. Suitable fragments were isolated using Agencourt AMPure XP beads (Beckman Coulter, Inc., Brea, CA, USA) and subsequently enriched through PCR amplification. The amplification procedure consisted of an initial denaturation at 98 • C for 30 s, 14 cycles of denaturation at 98 • C for 10 s, annealing at 65 • C for 30 s, and extension at 72 • C for 30 s, followed by a final extension at 72 • C for 5 min, and the process concluded with a 4 • C hold. Finally, the cDNA libraries of tobacco stems were sequenced by an Illumina HiSeq 2500 sequencing platform (Illumina, Inc., San Diego, CA, USA) using a flow cell.

Global and Differential Gene Expression Analysis of RNA-Seq Data
The reads were initially filtered to obtain clean reads by removing reads containing linkers and those of low quality, defined as reads with a proportion of undetermined bases exceeding 10% or with more than 50% of bases having a quality score ≤ 10, so as to ensure the accuracy of subsequent analyses [51]. The clean reads were then mapped to the tobacco cultivar TN90 genome (https://www.ncbi.nlm.nih.gov/bioproject/208209, accessed on 29 May 2014) using the HISAT2 program [52]. StringTie software (v1.3.3b) [53] was used to construct transcripts and evaluate gene expression. The Pearson correlation analysis was conducted on the expression levels of paired samples [54], and the coefficient of correlation (r) was used to evaluate the correlation strength. BLAST (v2.10.0) software [55] and various databases, including NR (non-redundant protein sequence database), Swiss-Prot, GO, COG, KOG, Pfam, and KEGG, were employed to annotate the genes. DEGseq [56] was used to identify DEGs between sample groups, according to the screening criteria of log 2 FC ≥ 2 and FDR < 0.01. The results of differential expression analysis and the interaction pairs included in the STRING database [57] were integrated to construct a DEG interaction network, which was then imported into Cytoscape version 3.8.1 [58] for visual analysis.
ClueGO was used to further analyze the functional grouping network of terms or pathways of the DEG sets, showing only pathways with p ≤ 0.05 [59]. The overrepresented GO terms in the network were identified and displayed as a network of significant GO terms using BiNGO with a significance level of 0.05 [60].

RT-qPCR Analysis
Total RNA was extracted from fresh plant tissue using a plant RNA extraction kit (Chengdu Biofit Biotechnologies Co., Ltd., Chengdu, China). First-strand cDNA was synthesized according to the manufacturer's protocols, using the HiScript III RT SuperMix for qPCR Kit with gDNA wiper (Vazyme). RT-qPCR analysis was performed using the AceQ qPCR SYBR green master mix (Vazyme) on a CFX Connect Real-Time System (Bio-Rad, Hercules, CA, USA). The primers used in RT-qPCR and RT-PCR are presented in Table S8. Two independent biological replicates and three technical replicates were used for each sample, and data were analyzed using the 2 −∆∆CT method [61]. The L25 ribosomal and Elongation factor 1a genes were used as internal reference genes [62] to ensure the accuracy of subsequent analyses.