Formation and growth of tiller buds in Paphiopedilum
Mature Paphiopedilum plants typically have 4–5 distichous leaves that obscure the stem bases. Tiller buds could arise from the axial of each leaf on its abbreviated stem. Figure 1 shows the successive steps of tiller formation in two Paphiopedilum plants with different tillering abilities. Firstly, a tiller bud formed from the axillary meristem (Fig. 1a, h). Then the tiller bud developed the first, second, and third leaf primordium step by step (Fig. 1b-d, i-k). When the tiller bud developed the fourth leaf primordium (indicated by blue markers in Fig. 1e, l), it has matured and acquired the ability to rapidly elongate, break through the mother plant, and become a normal tiller. Each tiller bud develops sequentially as the leaves grow, with the tillering bud at the lowest stem node (designated as tillering bud A, marked in blue in Fig. 1) forming and maturing first, followed by tiller bud B (marked in yellow) and tiller bud C (marked in green) further up the stem. The period from formation to maturation of a tiller bud typically takes 3–4 months, and the ability of the mature tiller bud to outgrow is crucial for determining the number of tillers. The outgrowth of a tiller bud in Paphiopedilum involves two consecutive processes: breaking dormancy and sustained outgrowth. For P. callosum, the tiller ability is weak, and tiller buds after maturation enter a dormant period of 5–7 months. After the apical meristem (Fig. 1 marked in pink) changed from vegetative growth to reproductive growth and bolting, the mature tiller buds began to outgrowth, breaking through the mother plant (Fig. 1f). As a result, P. callosum usually exhibits a one-branch flowering phenotype. (Fig. 1g). In contrast, P. ‘SCBG Yingchun’, which has a stronger tillering ability, shows continuous elongation and growth of tiller buds after maturation, in sync with the growth of the main stem leaves, without dormancy (Fig. 1m). Before the formation of the flower organs, 2–3 tiller seedlings could be produced, ultimately resulting in a multi-branched flowering phenotype (Fig. 1n)
Exogenous BAP application could effectively promote both the tiller buds and roots development of P. callosum
Exogenous BAP treatment at concentrations of 100, 200, 400, and 1000 mg/L can promote the outgrowth of tiller buds in P. callosum. After 100 d of treatment, the formation of tiller buds can be observed in the mother plants (Fig. 2). Different concentrations of exogenous BAP have different promoting effects on the growth of tiller buds. At a BAP concentration of 400 mg/L, the number of tiller buds significantly increased to 3.41 after 100 d of treatment, with an average height of the first tiller bud at 4.32 cm. Increasing the BAP concentration to 1000 mg/L did not enhance the promotion effect on tiller buds compared to the 400 mg/L treatment (the number of tillering buds was 2.83, and the height of the first tiller bud was 3.30 cm at 100 days after treatment) (Table 1). Therefore, we selected the 400 mg/L BAP treatment as the optimal treatment for promoting tillering growth in P. callosum.
After 30 d of 400 mg/L BAP treatment, four protruding bud points can be seen on the basal stem of P. callosum after peeling off the outermost 1–2 leaves. This stage is referred to as the dormancy release stage (Fig. 2d-1), Subsequently, these tiller buds continue to expand and turn green, breaking through the mother stem at 60 days. This stage is denoted as the expansion stage (Fig. <link rid="fig2">2</link>d-2). Tiller buds continue their sustained outgrowth, following the development of new leaves, ultimately becoming new tiller shoots (Fig. 2d-3), which is referred to as the sustained outgrowth stage. Low concentrations (100 and 200 mg/L) of BAP treatment can promote the production of new roots, although the promotion effects on tiller buds are not as significant as with the high concentrations. Plants subjected to 100 and 200 mg/L BAP treatment at 100 d can produce new roots, the number of the new root per plant being 1.25 ± 0.16 and 1.46 ± 0.14, respectively (Fig. 2b-<link rid="fig3">3</link>, c-<link rid="fig3">3</link>).
Effects of 400mg/L BAP treatment on endogenous hormone content at dormancy release stage.
Morphological results of exogenous BAP treatment showed that tiller buds could break dormancy and start to grow after 30 days of 400mg/L BAP treatment. To determine whether the promotion of tiller growth by BAP is associated with changes in endogenous levels of key hormones, we quantified the content of various hormones at 30 days in both control and BAP-treated plants (Fig. 3).
The content of ABA, GA7, and isoprenoid CTKs (iP, tZ, and cZ) was significantly reduced (by 66%, 52%, 73%, 85%, 49%) due to BAP application compared to the control plants. However, the content of IAA and aromatic CTKs (BAP, oT, and pT) in the BAP-treated stems was significantly higher compared to that of the control plants, which were 10.48, 19.40 and 107.53-fold of those in control plants, respectively.
Transcriptomic analysis of tiller bud development in 400 mg/L BAP treatment
Based on the morphological characterization of tiller buds under 400mg/L BAP treatment, we assigned the stages of tiller buds treated with BAP at 0d, 30d, 60d and 100d as stages 0, 1, 2, 3, respectively. To further investigate the transcriptomic regulation mechanisms of tiller growth by BAP treatment, RNA-seq analysis was performed at these four developmental stages. Clean data (78.94 Gb in total) were obtained from 12 libraries (1 treatment × 3 biological replicates × 4 developmental stages). Detailed statistics of the clean reads are listed in Supplementary Data Table S1. By de novo assembly of the total clean reads, we obtained 54,728 unigenes with an average length and N50 length of 1866.14 and 2692 bp, respectively (Supplementary Data Table S2). Among these unigenes, 36,408 (66.53%) were annotated with at least one putative function in publicly available databases (Supplementary Data Table S3). Samples were separately clustered based on the developmental stages in the principal component analysis (PCA) (Supplementary Data Fig. S2a). Sequence comparisons showed a high similarity between P. callosum transcripts and those of Dendrobium catenatum (41.48%), Phalaenopsis equestris (11.18%), and Apostasia shenzhenica (9.66%) (Supplementary Data Fig. S2b).
Differentially expressed genes (DEGs) were detected using a screening criteria of False Discovery Rate (FDR) < 0.01 and log2 ratio ≥ 1.0 (Fold Change ≥ 2.0). Three comparison groups were designed between successive tiller growth stage. The S1 versus S2 group had the fewest DEGs, with 1227 DEGs, consisting of 809 up-regulated and 418 down-regulated genes. The S2 versus S3 group had the highest number of DEGs, with 2097 in total, including 1009 up-regulated and 1088 down-regulated genes (Fig. 4a). There were 80 overlapped DEGs shared by the three comparison groups (Fig. 4b), and 3797 DEGs were expressed in at least one comparison group. The expression patterns of these 3797 DEGs were distinct among the four stages (Fig. 4c). K-means clustering analysis revealed 15 distinct clusters, named K1-K15 (Fig. 4d, Supplementary Data Table S4), which provided insights into the potential roles of the DEGs. The highly expressed genes in S0 (represented by K1-K4), S1 (represented by K5-K8), S2 (represented by K9-K12), and S3 (represented by K13-K15) were primarily involved in DNA binding, ribosomes, starch and sucrose metabolism, and photosynthesis, respectively.
KEGG enrichment analysis of differentially expressed genes
KEGG enrichment analysis was conducted on the three comparison groups. In the S0 versus S1 group, the top-enriched KEGG pathways included ribosome, DNA replication, cutin, suberine and wax biosynthesis, and zeatin biosynthesis (Fig. 5a). Highly enriched pathways in the S1 versus S2 group were photosynthesis-antenna proteins, DNA replication, carotenoid biosynthesis, and circadian rhythm - plant (Fig. 5b). Highly enriched pathways in the S2 versus S3 group were cutin, suberine and wax biosynthesis, phenylpropanoid biosynthesis, DNA replication, and tyrosine metabolism (Fig. 5c).
These results indicate that after BAP treatment, the plant hormone signaling transduction, ribosome, and DNA replication pathways are activated in the early stage (S1), enhancing the plant's activity, which may promote the activation of axillary buds. In the mid-stage (S2), genes related to starch and sucrose metabolism, Nitrogen metabolism, and Fatty acid degradation were significantly enriched, suggesting that sugars, nitrogen, and fatty acid may be involved in the process of tiller bud breakthrough from the mother stem and further elongation in P. callosum. As the tiller buds enter the S3 stage, they require more energy for sustained growth, leading to the significant enrichment of pathways related to photosynthesis, photosynthesis-antenna proteins, and phenylpropanoid biosynthesis.
Expression profile of CTKs biosynthesis, degradation and signaling genes
CTKs are synthesized by Adenylate isopentenyltransferases (IPTs), cytochrome P450 monooxygenase (CYP735A) and riboside 5’-monophosphate phosphoribohydrolase (LONELY GUY, LOG) [19–21]. The deactivation of CTKs is accomplished by cytokinin oxidase/dehydrogenase (CKX) [22]. Cytokinin signaling is mediated by a two-component system. A- or B-type response regulators (A-type RRs or B-type RRs, respectively) negatively regulated CTKs signal transduction [23].
We found 4 homologs of IPT, 9 homologs of CKX, but only one homolog of CYP735A and LOG (Fig. 6). Three cytokinin-synthesis genes, IPT1, CYP735A1, and LOG3 shows differentially expressed patterns. Expression level of IPT1 was significantly down-regulated at S3, while CYP735A1 exhibited the opposite trend of expression. LOG3 was significantly up-regulated by BAP treatment at S1 and S3. Most members of CKX gene family exhibited up-regulating responses to BAP treatment.
P. callosum has nine paralogs of genes for HK, five for HP, seven for A-type RRs and nine for B-type RRs (Fig. 6, Supplementary Data Table S5). Our transcriptome profiling demonstrated that six genes coding for B-type RRs were significantly up-regulated by BAP treatment at S1, S2 and S3. However, expression level of two A-type RR genes significantly down-regulated at S1, S2 and S3. The decreased level of isoprenoid CTKs (Fig. 3) at S1, corresponded to the up-regulation of CTKs catabolism genes and CTKs signaling genes at S1.
Expression profile of DEGs related to biosynthesis, degradation and signaling of auxin and ABA.
Besides CTKs, auxin and ABA also play pivotal roles in tiller development [24, 25]. A total of 28 DEGs involved in auxin biosynthesis, degradation, transport and signaling transduction were annotated (Fig. 7a and Supplementary Data Table S6). The IPA pathway is a mainly conserved pathway for most auxin synthesized in plants [26]. Flavin monooxygenase (YUCCA) and tryptophan amino acid transferase-1 (TAA1) are two key enzymes that regulated via the IPA pathway in a two-step process [27]. DIOXYGENASE OF AUXIN OXIDATION (DAO) proteins is essential for keeping auxin homeostasis, which is responsible for catalyzing the conversion of active IAA into biologically inactive 2-oxoindole-3- acetic acid (OxIAA) [28]. PIN-FORMED (PIN) proteins, ATP-BINDING CASSETTE SUBFAMILY B (ABCB) proteins, AUXIN RESISTANT1/LIKE AUXIN RESISTANT1 (AUX1/LAX)proteins and PIN-LIKES (PILS) proteins are key auxin transporters involved in diverse transporter systems that control short-range and long-range auxin transport [29, 30]. Our studies showed that a single gene for YUCCA were up-regulated during S1-S3, While a single gene for DAO exhibited the opposite trend of expression. Four homologs of auxin transporter genes (PIN1C, ABCB19, PILS1, LAX3) were all up-regulated in response to BAP during S1-S3.
Auxin-responsive gene expression relies on the TRANSPORT INHIBITOR RESPONSE1/AUXIN SIGNALING F-BOX (TIR1/AFB) pathway to initiate auxin signaling [31]. The downstream of IAA signaling pathway includes the auxin response factor (ARF), Small Auxin Up RNA (SAUR), AUXIN/INDOLE-3ACETIC ACID (Aux/IAA), Gretchen Hagen3 (GH3) genes [32]. There are eight homologs for IAA/AUX, eight homologs for ARF, three homologs of GH3 genes and three homologs of SAUR genes. The expression of these genes displayed a complex expression profile. One gene for AFR(ARF19) was significantly down-regulated during S1-S3, while other homologs of AFR genes did not show synchronized expression profile. High level expression of IAA synthesis and transport genes and low level expression of degradation gene DAO during S1-S3 suggests that IAA functions actively in tiller outgrowth of P. callosum.
Zeaxanthin epoxidase (ZEP) and 9-cis-epoxycarotenoid dioxygenase (NCED) are two enzymes that regulate ABA biosynthesis. ABA levels in plants are not only regulated by synthesis but also by metabolism. CYP707A encoded ABA-8’-hydroxylase, which deactivated ABA [33]. PYR/PYL/RCARs, as ABA receptors, bind ABA in a complex with ABI1 or other phosphoprotein phosphatase 2Cs (PP2Cs). SNF1-related protein kinases (SnRKs) and ABRE-binding factors (ABFs) are core components that complete the downstream of ABA signaling [34]. In our analysis, there are three DEGs for ABA biosynthesis, two genes for ZEP up-regulated during S1-S3, however, one gene for NCED down-regulated during S1-S3. Four CYP707A family genes had different expression profile during S0-S3, in which CYP707A7 was significantly up-regulated during S1-S3. Six DEGs in the ABA signaling pathway showed complex expression patterns. The expression level of gene coding for PLY4 and SAPK7 showed synchronized expression patterns, which were up-regulated during S1-S3, while one gene for ABF2 exhibited the opposite trend of expression (Fig. 7b and Supplementary Data Table S6).
Comparison of expression patterns of tiller-related candidate genes in two Paphiopedilum species with different tillering abilities
Unlike P. callosum, P. ‘SCBG Yingchun’ exhibits stronger tillering ability, where the outgrowth of tiller buds is not inhibited by the main shoot. When the main shoot grows to 5 leaves, the tiller buds can protrude from the main shoot and develop into normal tillers without the need for the main shoot to bolt or exogenous BAP treatment (Fig. 1). Figure 8 shows the continuous processes of the growth of the primary tiller bud in P. ‘SCBG Yingchun’. Similar to P. callosum treated with 400 mg/L BAP, the growth of the tiller buds in P. ‘SCBG Yingchun’ can also be divided into four corresponding stages.
To investigate the expression of key genes during this process, we selected several DEGs that may be related to tiller bud outgrowth though in-depth analysis of transcriptome data. qRT-PCR was conducted in this study. Dormancy-associated gene-1 (DRM1) and dormancy-associated MADS-box/SHORT VEGETATIVE PHASE-like genes (DAM/SVP), known as negative regulators for axillary bud outgrowth and branching [35, 36]. As shown in Fig. 8, compared to the non-BAP-treated plants of P. callosum, BAP-treated plants of P. callosum showed decreased expression of these two dormancy-associated genes during bud outgrowth process. In P. ‘SCBG Yingchun’, the relative mRNA level of PcDRMH1 showed more than 7-fold lower than that in control plants of P. callosum during stage 0–3, and the expression level of PcSVP gene followed the same trend as that in BAP-treated plants of P. callosum. This indicated that BAP may induce the release of dormancy in tiller buds and promote tiller development. Furthermore, the expression of the gene involved in ABA catabolism (PcCYP707A) was below the initial level at S0 but rose dramatically after BAP treatment, while the expression level of PcCYP707A in control plants of P. callosum and P. ‘SCBG Yingchun’ did not show significant differences during stage 0–3. The auxin transporter gene PcPIN1 exhibited a significantly up-regulated expression pattern during the tiller outgrowth process in BAP-treated P. callosum and P. ‘SCBG Yingchun’, compared with the control plants of. This suggests that BAP affects the regulation of ABA catabolism and IAA transport during tiller bud outgrowth stage.
Validation of gene expression using qRT-PCR
We selected 12 DEGs related to endogenous hormone and sugar synthesis, metabolism and transport, as well as axillary bud formation and maintenance, and verified the accuracy of transcriptome data by qRT-PCR (Fig. 9). Consistent with the transcriptome data, qPCR showed that the expression levels of two genes involved in SLs biosynthesis and signaling (CCD8B, D14), and 3 genes involved in sugar biosynthesis and transporter (TPS7, TPP6, STP7) were the highest in S2 stage among the four tiller bud development stages. We also used qPCR to evaluate the expression profiles of 4 genes involved in CTKs biosynthesis and catabolism (LOG3, CKX10, CKX5, RR2) and 3 genes involved in auxin biosynthesis and transport (YUC2, PIN1C, ABCB19). Two CTKs catabolism genes (CKX5 and CKX10) and one gene RR2, which is a negative regulator of CTKs signaling pathway,showed significant up-regulation after 400mg/L BAP treatment. one IAA synthesis gene and two transport genes showed the same trend. The expression of LOG3 increased first, then decreased and then increased in the four stages of tiller growth under 400mg/L exogenous BAP treatment. SCL18 (BMK_Unigene_356084) is a gene related to axillary bud formation, and the relative expression of qRT-PCR showed no significant difference in the four stages of tiller bud. The relative qRT-PCR expression of the gene TB1 (BMK_Unigene_308669), which inhibits axillary bud elongation and growth, was significantly down-regulated at S1 stage. The relative expression levels of the target genes obtained by qRT-PCR were consistent with the FKPM of each gene.