Comparative transcriptome analysis reveals the regulatory networks of cytokinin in promoting the floral feminization in the oil plant Sapium sebiferum
BMC Plant Biology volume 18, Article number: 96 (2018)
Sapium sebiferum, whose seeds contain high level of fatty acids, has been considered as one of the most important oil plants. However, the high male to female flower ratio limited the seed yield improvement and its industrial potentials. Thus, the study of the sex determination in S. sebiferum is of significant importance in increasing the seed yield.
In this study, we demonstrated that in S. sebiferum, cytokinin (CK) had strong feminization effects on the floral development. Exogenous application with 6-benzylaminopurine (6-BA) or thidiazuron (TDZ) significantly induced the development of female flowers and increased the fruit number. Interestingly, the feminization effects of cytokinin were also detected on the androecious genotype of S. sebiferum which only produce male flowers. To further investigate the mechanism underlying the role of cytokinin in the flower development and sex differentiation, we performed the comparative transcriptome analysis of the floral buds of the androecious plants subjected to 6-BA. The results showed that there were separately 129, 352 and 642 genes differentially expressed at 6 h, 12 h and 24 h after 6-BA treatment. Functional analysis of the differentially expressed genes (DEGs) showed that many genes are related to the hormonal biosynthesis and signaling, nutrients translocation and cell cycle. Moreover, there were twenty one flowering-related genes identified to be differentially regulated by 6-BA treatment. Specifically, the gynoecium development-related genes SPATULA (SPT), KANADI 2 (KAN2), JAGGED (JAG) and Cytochrome P450 78A9 (CYP79A9) were significantly up-regulated, whereas the expression of PISTILLATA (PI), TATA Box Associated Factor II 59 (TAFII59) and MYB Domain Protein 108 (MYB108) that were important for male organ development was down-regulated in response to 6-BA treatment, demonstrating that cytokinin could directly target the floral organ identity genes to regulate the flower sex.
Our work demonstrated that cytokinin is a potential regulator in female flower development in S. sebiferum. The transcriptome analysis of the floral sex transition from androecious to monoecious in response to cytokinin treatment on the androecious S. sebiferum provided valuable information related to the mechanism of sex determination in the perennial woody plants.
Sapium sebiferum, a member of the Euphorbiaceae family, is one of the most important perennial oil plants. It has attracted great attention due to its high oil production, pharmaceutical values, high ornamental values as landscape plant and its viability to grow in the marginal land (Peng et al. 2008; Yang et al., 2007). S. sebiferum has been cultivated as a valuable oil plant for more than 14 centuries, due to the high oil content of the seed coat and seed kernel, which can be used for the manufacture of industrial products, such as lubricant, soap and candles . Aside from being a source for industrial purposes, the oil extracted from the S. sebiferum is also considered as a kind of renewable oil resources, which can be further exploited for the production of diesel fuel and fatty acid alkyl esters [2, 3]. Nevertheless, low female to male flower ratio significantly limited the seed yield improvement and the industrial potentials of S. sebiferum. Most of the S. sebiferum plants are monoecious with separate male and female flowers on the same inflorescence, while some are androecious plants that only produce male flowers. There are over one hundred male flowers arrange in the narrow raceme-like inflorescence, while only with few female flowers growing at the base of the inflorescence. Thus, the manipulation of the flower sex differentiation by using the biochemical and genetic strategies was of significant importance in improving the seed yield in S. sebiferum.
Phytohormones are essential factors involved in the regulation of vegetative and reproductive growth in the plant kingdom. Phytohormones, such as cytokinin (CK), auxin, gibberellin (GA) and ethylene (ETH), were demonstrated to be involved in the regulation of flowering and sex determination in many species . Generally, cytokinin, auxin and ethylene were considered as positive regulators in gynoecium development, whereas GA (gibberellin) promoted the development of androecious organs [5,6,7,8]. In Jatropha curcas and Plukenetia volubilis, exogenous application with cytokinin can increase the female flower number [9, 10]. Nevertheless, the mechanisms how cytokinin regulated the sex determination remains largely unknown. In this study, we demonstrated that CK was a key regulator in sex determination in the oil plant S. sebiferum. CK treatment on the floral bud significantly induced the floral feminization and promoted the fruiting. As the seed yield improvement of S. sebiferum was greatly limited by the low female to male flower ratio, thus our findings will be of significant importance in the cultivation of S. sebiferum. Further, the feminization effect was also detected on the androecious genotype of S. sebiferum. To further investigate how CK regulated the transition of the floral bud from androecious to monecious, RNA-seq analysis was carried out on the floral buds of the androecious S. sebiferum at 6 h, 12 h and 24 h after 6-BA treatment. The present study not only provided an effective method in improving the seed yield during S. sebiferum cultivation, the transcriptome data also provided insights in understanding the mechanism of sex determination regulated by cytokinin in the woody plants.
Plant materials and growth conditions
In this experiment, three-year-old S. sebiferum trees, growing in the experimental field in Hefei Institutes of Physical Science (N 31° 49′, E117° 13′), Chinese Academy of Sciences, were used for the hormonal treatment. Trees were planted with 2.5 m × 2.5 m spacing, and fertilized every year before bud sprouting.
Phytohormone preparation and application
To make stock solutions (20 mM) of phytohormones, 6-BA, TDZ, trans-zeatin (tZ) or kinetin was dissolved in 0.2 M NaOH solution. The stock solutions were diluted with water to make different concentrations of working solutions containing 0.05% Tween-20 . All working solutions for each treatment had the same solvent. The floral buds (one week after its formation) were treated with working solutions using a 100 ml plastic sprayer.
Sample collection and RNA isolation for transcriptome sequencing and quantitative real time PCR
Floral buds of the androecious plants (Additional file 1: Figure S1) were treated with 6-BA (200 μM) and mock solution (containing the same solvent with 6-BA solution). The floral buds after removing the leaves (approximately 3 mm in length, Additional file 1: Figure S1) were collected at 6 h, 12 h and 24 h separately, and then immediately frozen in liquid nitrogen. Each treatment at each time point was prepared with six independent biological replicates (three for RNA seq, and three for quantitative real time PCR). The total RNA was extracted by using a Plant Total RNA Isolation Kit (Omega, China). RNA purity was checked by gel electrophoresis and NanoDrop 2000c spectrophotometer (Thermo Fisher, MA, USA). RNA concentration and integrity were separately determined using a Qubit 2.0 Flurometer (Life Technologies, Carlsbad, CA, USA) and Agilent Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, CA, USA).
Transcriptome sequencing and de novo transcriptome assembly
Sequencing libraries (each sample with three biological replicates) were prepared following the instructions of the NEBNext Ultra RNA Library Prep Kit for Illumina (NEB, USA). The literary quality was determined using the Agilent Bioanalyzer 2100 system. The clustering of the index-coded samples was performed according to the instructions of cBot Cluster Generation System. Then the libraries were sequenced on the Illumina HiSeq 2000 platform at the Beijing Genomics Institute (BGI, Shenzhen, China) and the raw reads were generated in 90-bp paired-end format. The low-quality reads (< 20 Phred scores) were removed by Fastq_clean  and further assessed by FastQC . The clean reads were assembled by using Trinity (version 2.0.6) . The paired-end reads were then mapped to the de novo assemblies by using Bowtie (version 1.1.1) . The abundance estimation was performed using Corset (version 1.03) .
Transcript annotation and gene expression analysis
All the assembled transcripts were searched against the NCBI non-redundant protein database (Nr), Swiss-Prot protein database and NCBI non-redundant nucleotide sequence (Nt) database. The GO functional classification was then carried out using Blast2GO program based on the Nr annotation . The GO analysis was carried out by using the WEGO software . The transcripts were aligned to the Clusters of Orthologous Groups of proteins (COG) database and Kyoto Encyclopedia of Genes and Genomes (KEGG) database [19, 20]. To quantify the differentially expressed unigenes, the read counts were calculated by using the fragments per kilobase of transcript per million fragments mapped method (FPKM) .
Quantitative real time PCR (qPCR) analysis
cDNA was synthesized according to the instructions of EasyScript One-Step gDNA removal and cDNA synthesis SuperMix (Transgen Biotechnology, Beijing, China). qPCR was performed using TaKaRa SYBR Premix Ex Taq™ II (TaKaRa Biotechnology, Dalian, China) on the Roche Light Cycler 96 System (Roche, Swiss). The PCR procedures were as follows: 95 °C for 5 min; 45 cycles of 95 °C for 10 s, 60 °C for 10 s, and 72 °C for 15 s; the melting curve analysis was from 65 °C to 95 °C. The information of the primers was listed as Additional file 2: Table S1.
Data were analyzed using the Statistical Product and Service Solution (SPSS) version 13.0 software. The significance between the hormone treated and control group was determined using Student’s test.
6-BA and TDZ had strong feminization effects on the floral development and promoted fruiting in S. sebiferum
Based on the floral structure, S. sebiferum can be mainly divided into two main genotypes, S. sebiferum var. conferticarpa and S. sebiferum var. laxiarpa . Most of S. sebiferum var. conferticarpa plants are monoecious, with few female flowers growing at the floral base, while some are androecious plants that only produce male flowers (Additional file 1: Figure S1). The results showed that exogenous 6-BA or TDZ treatment significantly induced the female flowers, which appeared at the position where the male flowers located (Fig. 1). The feminization effect was increased with increasing concentrations of 6-BA (Fig. 1a, b), or TDZ (Fig. 1c, d), reaching an average of 57.7 and 78.5 female flowers per inflorescence, separately by 6-BA and TDZ treatment. TDZ treatment had better effects in promoting the floral development, resulting in larger floral size and higher female flower number than that by 6-BA treatment (Fig. 1). However, over-dosed 6-BA or TDZ (e.g. 1 mM) was actually toxic to the floral bud (Fig. 1e). Approximately 20–30% of the 6-BA- or TDZ-induced female flowers could further develop into fruits (Fig. 2). It was worthy to mention here that although trans-zeatin (tZ) or kinetin were two other types of cytokinins. Exogenous application with these tZ or kinetin had no effects on the flower development or fruiting in S. sebiferum (Figs. 1f and 2d, e).
6-BA or TDZ treatment feminized the androecious genotype of S. sebiferum
In S. sebiferum, there are a small population of androecious plants that only produces male flowers and bares no fruits (Additional file 3: Figure S2). The effects of cytokinin on the floral development on the androecious genotype of S. sebiferum was also investigated in this study. The results showed that 6-BA or TDZ treatment can induce the female flower development on the male inflorescence of the androecious S. sebiferum (Fig. 3a, b). The induced female flowers, which had similar structures with that of the normal flowers (Fig. 3c), were located at the position of the male flowers all over the whole inflorescence, and can further develop into fruits (Fig. 3d).
Transcriptome sequencing of the floral buds of the androecious S. sebiferum in response to 6-BA treatment
The androecious genotype provided great advantages in investigating the mechanism of sex differentiation regulated by cytokinin in S. sebiferum. After 6-BA or TDZ treatment, abundant female flowers were induced in the “male inflorescence”, thus we presumed that cytokinin was not only the key factor in regulating the sex differentiation in S. sebiferum, but also possibly participated in the formation of the androecious plants. The transcriptome analysis of the floral buds at the early developing stages of the androecious plants would be helpful for identification of key genes involved in the female flower formation and development in response to cytokinin in S. sebiferum.
In S. sebiferum, the floral buds were located at the terminal of each branchlet. One week after the floral bud formation (approximately 3 mm in length, Additional file 1: Figure S1), 6-BA (200 μM) and mock solutions were directly applied to the floral buds. In total of six groups of samples (6-BA- and mock-treated floral buds at 6 h, 12 h and 24 h, each group contains six biological replications) were collected for the following transcriptome sequencing and quantitative real-time PCR, aiming at identification of the early responsive genes that were involved in the sex determination regulated by cytokinin. Eighteen cDNA libraries were constructed and sequenced on the Illumina high-throughput sequencing platform. The detailed information of the transcriptome assembly can be found in Table 1 and Additional file 4: Tables S2-S4. The differences of gene expression levels were plotted with colored dots and there were separately 129, 352 and 642 genes differentially expressed at 6 h, 12 h and 24 h after 6-BA treatment (Fig. 4, Additional file 5: Tables S5-S7). To validate the accuracy of the transcriptome sequencing results, the expression level of fourteen differentially expressed genes (DEGs) was determined by quantitative real-time PCR (qPCR). The expression profiles of these genes by qPCR were consistent with the values derived from the RNA-seq (Table 2). Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analysis were further carried out to generally describe the biological functions of these differentially expressed genes (DEGs) (Fig. 5 and Additional file 3: Figure S2). The results showed that the DEGs were mainly classified into 19 categories by KEGG analysis and the predominant categories were mostly related to “Metabolism” (58%) and “Genetic information processing” (21%) (Fig. 5). Further, the predominant categories by GO analysis were “metabolic process”, “cellular process”, “cell”, and “catalytic activity” (Additional file 6: Figure S3), which was closed correlated to the basic biological functions of cytokinin in promoting cell growth, replication and metabolism.
Differentially expressed genes related to cell cycle and nutrient translocation
The KEGG analysis revealed that abundant DEGs were classified into “Metabolism” and “Genetic information processing” categories (Fig. 5). Many cell cycle related DEGs, including four cyclin genes, CYCD3 (CL6841.Contig3_All), CYCD6 (CL3342.Contig16_All and CL3342.Contig6_All), CYCH1 (CL384.Contig6_All), and some other DEGs, such as E2F Transcription factor 3 (E2F3, CL1239.Contig4_All), ANAC068 (CL12342.Contig2_All), Block of Cell Proliferation 1 (BOP1, Unigene9027_All), NEDD1 (CL659.Contig13_All) and Origin of Replication Complex 1B (ORC1B, CL2740.Contig8_All) which encodes key transcription factors that involved in the regulation of cell division, were up-regulated after 6-BA treatment (Fig. 6a). Several DEGs that encoded DNA and RNA polymerases (Unigene7546_All, CL14256.Contig6_All, CL14312.Contig3_All and CL8049.Contig5_All) were also up-regulated in response to 6-BA treatment (Fig. 6a). In accompany with the enhanced cell division, the nutrients homeostasis in the floral bud was also affected by 6-BA treatment. The expression level of the nutrient translocation-related genes including Early Response to Dehydration 6 (ERD6, CL1189.Contig43_All), Glucose-6-Phosphate Translocator 2 (GPT2, CL17101.Contig1_All), Nucleotide/Sugar Transporter (NST, CL17241.Contig33_All) and UDP-Galactose Transporter 7 (UTR7, CL648.Contig66_All) that involved in sugar transportation, Oligopeptide Transporter 4 (OPT4, Unigene14375_All) in oligopeptide transportation, Cationic Amino Acid Transporter 7 (CAT7, Unigene35653_All) in amino acid transportation, ENOR3L5 (CL15097.Contig4_All) in magnesium transportation, Ammonium Transporter 2 (AMT2, CL16110.Contig3_All) in ammonium transportation, Embryo Defective 104 (EMB104, CL11518.Contig3_All) in nucleotide transportation, Mitochondrial Phosphate Transporter 3 (MPT3, CL17537.Contig2_All) in phosphate transportation, two Calcium Transporting ATPases ACA10 (CL1699.Contig5_All) and ECA3 (CL1459.Contig14_All) in calcium transportation, and Zinc-Induced Facilitator 2 (ZIF2, CL1189.Contiga42_All and CL1189.Contig17_All) in zinc transportation, were mostly up-regulated in response to 6-BA treatment (Fig. 6b). Recent studies demonstrated that trehalose-6-phosphate (T6P) was directly involved in the regulation of flowering, the biosynthesis of which was tightly regulated by the Thehalose-6-Phosphate Synthase (TPS) . The transcriptome data also showed that four TPS genes, TPS1 (CL7084.Contig2_All and CL7084.Contig4_All), TPS9 (CL9994.Contig3_All) and TPS10 (CL9994.Contig2_All), were differentially regulated by 6-BA treatment (Fig. 6b).
Differentially expressed genes related to cytokinin metabolism and signaling
After gene annotation, many cytokinin metabolism- and signaling-related genes were identified to be differentially expressed in response to 6-BA treatment (Fig. 7). Adenine Phophoribosyltransferases (APTs) and Cytokinin Oxidases (CKXs) are two key regulators in maintaining the endogenous cytokinin homeostasis. APTs participated in the interconversion of cytokinin bases into nucleotides, while CKXs can directly catalyze the degradation of cytokinins [24, 25]. In this study, three APT genes ATP1 (CL5503.Contig22_All), APT5 (CL5503.Contig21_All and Cl5503.Contig27_All), and four CKX genes CKX3 (CL13408.Contig1_All), CKX5 (CL13173.Contig2_All and CL13173.Contig1_All), CKX7 (CL12919.Contig2_All), were significantly up-regulated at 6 h, 12 h and 24 h after 6-BA treatment (Fig. 7), suggesting a strong feedback regulation of CK homeostasis after 6-BA treatment.
Previous studies suggested that cytokinin signal transduction occurs though a two-component phosphorelay system. Arabidopsis response regulators (ARRs) can be classified into two groups (type A and type B), both of which are important components in cytokinin signal transduction . Our results showed that the expression of four type A ARR genes, ARR3 (CL7852.Contig3_All), ARR4 (CL7852.Contig4_All), ARR5 (CL1788.Contig4_All) and ARR9 (Unigene12447), was up-regulated at 6 h, 12 h or 24 h after 6-BA treatment (Fig. 7). As the expression of type B ARRs was not significantly affected after 6-BA treatment, it suggested that in S. sebiferum, cytokinin signal transduction in the floral bud might be dependent on the regulation of type A ARRs.
Differentially expressed genes related to hormonal biosynthesis, metabolism and signaling
Phytohormones, such as CK, auxin, ABA, ETH and GA are the key regulators in controlling floral initiation, development, and sex differentiation. We further investigated whether the biosynthesis and signaling of the other hormones were affected by 6-BA treatment. The results showed that three auxin metabolism- and six auxin signaling-related genes were differentially expressed at 6 h, 12 h or 24 h after 6-BA treatment (Fig. 7). YUCCA2 (YUC2) and YUCCA2 (YUC4) which encoded two YUC flavin monooxygenases, played essential roles in auxin biosynthesis . In S. sebiferum, both YUC2 (CL10701.Contig1_All) and YUC4 (CL10984.Contig2_All) was up-regulated at 12 h and reached a higher expression level at 24 h after 6-BA treatment (Fig. 7). Meanwhile, the expression of Dioxygenase for Auxin Oxidation 1 (DAO1, CL2608.Contig20_All), which encoded a key auxin oxidase and was also a major contributor to IAA degradation , was firstly decreased at 12 h, then increased at 24 h after 6-BA treatment (Fig. 7), suggesting there could be a transient increase of endogenous auxin level after 6-BA treatment. Further, three ABA biosynthesis- and five ABA signaling-related genes were also regulated by 6-BA treatment. The expression of ABA biosynthesis-related genes XERICO (Unigene42548_All and CL8997.Contig4_All) and NCED3 (CL9818.Contig6_All) was significantly down-regulated by 6-BA treatment (Fig. 7). Four ABA signaling-related genes Regulatory Components of ABA Receptor 10 (RCAR10, Unigene31612_All), Ring Domain Ligase 2 (RGLG2, CL5656.Contig10_All), Growth-regulating Factor 3 (GRF3, CL16961.Contig9_All) and Abscisic Acid Responsive Elements-binding Factor 3 (ABF3, CL9171.Contig5_All) were negatively regulated in response to 6-BA treatment. Additionally, the expression of several unigenes that involved in the regulation of hormonal signaling of ethylene (CL1635.Contig5_All, CL8680.Contig2_All, CL15911.Contig2_All and CL2328.Contig15_All), gibberellin (CL3574.Contig42_All, CL13937.Contig2_All and CL13937.Contig3_All) and jasmonic acid (CL2992.Contig3_All, CL1493.Contig4_All and CL1493.Contig2_All) was also altered in response to 6-BA treatment (Fig. 7). These results suggested that cytokinin might regulate the floral development via interactions with the other phytohormones in S. sebiferum.
Identification of the differentially expressed genes related to the floral development
The floral organ identity genes that were differentially regulated after 6-BA treatment could directly contribute to the altered sex differentiation in S. sebiferum. The transcriptome data showed that a total of 21 unigenes that closely related to flowering were found to be differentially expressed in response to 6-BA treatment (Fig. 8). Five genes that related to the early floral meristem development, such as PUCHI (CL1635.Contig5_All), Target of Early Activation Tagged 1 (TOE1, CL1615.Contig68_All and CL46.Contig48_All), Ethylene Response DNA Binding Factor 4 (EDF4, CL4670.Contig3_All) and WIGGUM (WIG, CL3863.Contig4_All), were differentially regulated by 6-BA treatment (Fig. 8). Further, several key genes that directly regulated the gynoecium development, such as SPATULA (SPT, Unigene42678_All), KANADI2 (KAN2, CL1205.Contig3_All), JAGGED (JAG, CL11939.Contig4_All), Cytochrome P450 78A9 (CYP78A9, CL15451.Contig2_All), and APETALA2 (AP2, CL8973.Contig2_All) were differentially regulated at 6 h, 12 h or 24 h after 6-BA treatment (Fig. 8). Moreover, a series of unigenes such as MYB Domain Protein 108 (MYB108, CL847.Contig2_All), Cadmium Induced (CDI, CL12892.Contig3_All), NPK1-Activating Kinesin 1 (NACK1, CL2667.Contig7_All), Shaker Pollen Inward K+ Channel (SPIK, CL3798.Contig10_All), Vacuolar Protein Sorting 34 (VPS34, Unigene6428_All), TATA Box Associated Factor II 59 (TAFII59, CL361.Contig25_All), FERONIA (FER, CL6663.Contig3-All), Embryo Defective 1674 (EMB1674, Unigene33664_All), PISTILLATA (PI, CL12404.Contig2_All), which were related to the development of the male organs, were also differentially expressed after 6-BA treatment (Fig. 8). Conclusively, the identification and characterization of the cytokinin-responsive DEGs related to the floral organ identities provided insights in understanding the sex determination in S. sebiferum.
Cytokinin has strong feminization effects on the flower development in both monoecious and androecious genotypes of S. sebiferum
The regulation of the female to male flower ratio is critical for improvement of the seed yield in the monoecious plants. Various phytohormones were demonstrated to be important regulators not only in vegetative growth, but also in the flower initiation and sex differentiation. In plants, phytohormones such as ethylene , gibberellins  and auxins [5, 31], have feminization effects in various species. Cytokinins are essential for numerous decisions throughout the plant developmental processes and adaptation to the biotic and abiotic environment . In Jatropha curcas, Plukenetia volubilis and Luffa cylindrica, exogenous CK treatment can increase the number of female flowers [9, 10, 33]. Our present study in S. sebiferum were consistent with these findings that 6-BA or TDZ treatment has strong feminization effects on the floral development (Fig. 1). Although tZ and kinetin were two other types of cytokinins that can induce the division of plant cells , exogenous application with tZ and kinetin had no effects on the floral development in S. sebiferum (Fig. 1f). We postulated that in S. sebiferum, the regulation of floral sex might be controlled by specific cytokinin types.
Our results further demonstrated that 6-BA or TDZ application can also effectively induce the female flowers on the androecious S. sebiferum (Fig. 3). These findings suggested that cytokinin was not only the key regulator in controlling the female flower development, but also possibly was the potential determinant that involved in the formation of the androecious plants in S. sebiferum. The transcriptome sequencing of the male floral buds subjected to 6-BA treatment was carried out and abundant cytokinin-responsive genes that were possibly involved in female flower induction were further identified. The floral buds of the androecious S. sebiferum seemed to be less sensitive to 6-BA treatment in the transcriptomic level, since the number of the differentially expressed genes in response to 6-BA treatment (especially at 6 h) was fewer (Fig. 4), compared with the transcriptome results from the other plant species [7, 35,36,37]. We postulated that the formation of the androecious plants in S. sebiferum was possibly due to the genetic mutations that caused the less-sensitivity to cytokinins. Nevertheless, exogenous 6-BA treatment can still effectively induce the female flower development on the androcious S. sebiferum (Fig. 3). Transcriptome analysis of the transition of the inflorescences from androecious to monoecious on the androecious plants will provide valuable information of the regulative gene network of the cytokinin-regulated sex differentiation in S. sebiferum.
Cytokinin regulated the cell cycle and nutrient translocation
The originally identified function of cytokinin was to promote the cell division . In this study, ten cell cycle-related genes were identified in response to 6-BA treatment in the transcriptome database (Fig. 6a). Three D-type CYCD6 (CL3342.Contig16_All and CL3342.Contig11_All), CYCD3 (CL6841.Contig3_All) and one H-type CYCH1 (CL384.Contig6_All) cyclin genes that played important roles in the regulation of cell division , were up-regulated by 6-BA treatment (Fig. 6a). E2F Transcription Factor 3 (E2F3) was a key component in the cyclin D/retinoblastoma/E2F pathway, was also up-regulated at 24 h after 6-BA treatment . Origin of Replication Complex 1B (ORC1B) that regulates the initiation of DNA replication , was up-regulated at 24 h after 6-BA treatment (Fig. 6a). NEDD1 encoding a WD40 repeat protein that plays a critical role in the cell mitosis in interactions with tubulin complex , was significantly up-regulated by 6-BA treatment (Fig. 6a). Two other transcription factors Block of Cell Proliferation 1 (BOP1) and NAC Domain Containing Protein 68 (NAC068) that involved in the cell cycle regulation [43, 44], were also up-regulated in response to 6-BA (Fig. 6a). Additionally, several genes encoding the DNA/RNA polymerases, such as Nuclear RNA Polymerase C1 (NRPC1, CL2915.Contig1_All), Poly(ADP-Ribose) Polymerase 1 (PARP1, CL8049.Contig5_All), NRPB4 (CL4312.Contig3_All) and Poly(A) Polymerase 2 (PAPS2, CL14256.Contig6_All), were also up-regulated by 6-BA treatment (Fig. 6a). These results suggested that exogenous application with cytokinin could cause immediate cell replication and organ enlargement via regulating the cell-cycle related genes in the floral buds in S. sebiferum.
Further, a few genes that involved in the translocation of the nutrients, such as sugar, peptide, magnesium, ammonium, nucleotide, calcium, zinc and phosphate, were also differentially expressed in response to 6-BA treatment (Fig. 6b). The altered nutrients transportation to the floral bud after 6-BA treatment could be associated with the vigorous female flower induction in the floral bud. It has been demonstrated that sugars, such as sucrose, glucose and trehalose-6-phosphate were the initiation signals involved in the flower initiation in Arabidopsis . Trehalose-6-Phosphate Synthases (TPS) were the key enzymes involved in trehalose-6-phosphate biosynthesis. The homologs of TPS1 in S. sebiferum were up-regulated, whereas TPS9 and TPS10 were down-regulated in response to 6-BA treatment (Fig. 6b). Several genes such as Early Response to Dehydration 6 (ERD6), Glucose-6-Phosphate Translocator 2 (GPT2), Nucleotide/Sugar Transporter (NST) that involved in the sugar translocation [45,46,47], were also differentially regulated by 6-BA treatment (Fig. 6b). These results suggested that the cytokinin-responsive cell cycle and cell metabolism-related genes could be involved in the regulation of flowering.
Cytokinin regulated the biosynthesis, transportation and signaling of other phytohormones in the regulation of sex determination in S. sebiferum
Phytohormones, such as auxin, cytokinin (CK), abscisic acid (ABA), ethylene (ETH), brassinosteroid (BR), gibberellin (GA), jasmonic acid (JA) and salicylic acid (SA), coordinate plant growth and development by modulating various cellular changes in response to the environmental or intrinsic signals . The transcriptome data showed that many genes that were directly related to the biosynthesis, transportation and signaling of other phytohormones, were differentially regulated by 6-BA treatment in the floral buds of the androecious S. sebiferum. Specifically, ten genes were identified in auxin biosynthesis and signaling, eight genes in ABA biosynthesis and signaling, four genes in ethylene signaling, three gene in GA signaling, and three genes in JA signaling (Fig. 7). Previous studies demonstrated that cytokinin and auxin mutually regulated their biosynthesis and signaling in the regulation of plant growth and development . Exogenous application with cytokinin in young shoot up-regulated the expression of auxin biosynthesis genes, such as YUCCA (YUC) . In this study, two YUC genes YUC2 and YUC4 were found to be significantly up-regulated at 6 h, 12 h and 24 h in response to 6-BA treatment (Fig. 7). Ring Domain Ligase 2 (RGLG2) was previously reported to be involved in the regulation of auxin metabolism. Mutation of RGLG2 could cause significant decrease of auxin content in both leaf and root in Arabidopsis . The expression of RGLG2 was increased in response to 6-BA treatment in S. sebiferum (Fig. 7). The results also showed that the expression of Dioxygenase for Auxin Oxidation 1 (DAO1), which encoded a major IAA oxidase in plants , was found to be down-regulated at 12 h, whereas promoted at 24 h (Fig. 7). These results suggested that the auxin biosynthesis was transiently promoted, at least in a short time after 6-BA treatment at the floral buds. Additionally, several auxin signaling and transportation-related genes, such as Auxin Response Factor 2 (ARF2), Indole-3-Acetic Acid Inducible 13 (IAA13), Indole-3-Acetic Acid Inducible 26 (IAA26), GH3.9, Auxin Resistant 3 (AXR3) and Protein Phosphatase 2A (PP2A) were found to be mostly up-regulated in response to 6-BA treatment (Fig. 7). It has been demonstrated that auxin was important for the female reproductive organ development in Arabidopsis , and also exhibited a strong feminization effect on Cannabis sativus and Opuntia stenopetala [5, 31]. The expression of auxin signaling genes was also positively correlated with the transition to female flowers on the monoecious Ricinus communis . Our results were in accord with these findings that auxin could also play important role in the female flower induction regulated by cytokinin in S. sebiferum.
In cucumber, ABA did not change the sex expression, nevertheless it promoted the female tendency of gynoecious plants . XERICO and Nine-cis-epoxycarotenoid Dioxygenase 3 (NCED3) are two important regulators in ABA biosynthesis [54, 55]. The transcriptome data showed that these two genes were significantly down-regulated by 6-BA treatment (Fig. 7). The ABA signaling-related genes, such as Ring Domain Ligase 2 (RGLG2), ATMYB44 (MYB44), Growth Regulating Factor 3 (GRF3) and Abscisic Acid Responsive Elements-binding Factor 3 (ABF3) were down-regulated by 6-BA treatment (Fig. 7). Previous studies showed that ABA and cytokinin antagonistically regulated many aspects of plant growth and development [56, 57]. In S. sebiferum, the cytokinin responsiveness of the ABA biosynthesis and signaling-related genes suggested an antagonistic interactions between CK and ABA in the regulation of the floral development.
Cytokinin also interacted with ETH, GA and JA in the regulation of flowering . PUCHI, a key component in ethylene-activated signaling pathway, was required for morphogenesis in the early floral meristem growth . The transcript abundance of PUCHI in the floral bud was up-regulated in response to 6-BA treatment (Fig. 7). The S. sebiferum orthologs of three additional ethylene signaling-related genes, Ethylene Response Factor 74 (ERF74), Ethylene and Salt Inducible 3 (ESE3) and Transcription Factor IID (TFIID) were also differentially regulated by 6-BA treatment (Fig. 7), suggesting a potential interactions between CK and ETH in the regulation of sex determination in S. sebiferum. Among these differentially expressed genes after 6-BA treatment, three JA signaling genes ATMYC2 (MYC2), Jasmonate Insensitive 1 (JAI1), BHLH14 and three GA signaling genes GA-Stimulated Arabidopsis 1 (GASA1–1 and GASA1–2) and Glabrous Inflorescence Stems 2 (GIS2) were differentially regulated in response to 6-BA treatment (Fig. 7). Conclusively, these results implied that cytokinin interacted with auxin, ABA, ethylene, JA and GA in sex determination in S. sebiferum (Fig. 9).
Cytokinin regulated the expression of gynoecium and androecium identity-related genes
6-BA treatment profoundly induced the formation of female flowers in the androecious S. sebiferum (Fig. 3). The differential expression of ABCE floral organ-identity genes has been implicated in the regulation of male and female flower development [60, 61]. Thus it suggested that cytokinin could directly target the floral organ-identity genes in the regulation of sex determination. A total of 21 flowering-related genes were identified in response to 6-BA treatment at the floral bud of the androecious plants (Fig. 8). Several early floral meristem growth-related genes, such as PUCHI, Target of Early Activation Tagged 1 (TOE1), Ethylene Response DNA Binding Factor 4 (EDF4) and WIGGUM (WIG), were differentially regulated in response to 6-BA treatment (Fig. 8), suggesting cytokinin was directly involved in the regulation of the early floral meristem development in S. sebiferum. Among the DEGs, a few genes that were closely related to the gynoecium and androecium development, were found to be differentially regulated in response to 6-BA treatment. SPATULA (SPT), a bHLH transcription factor, was involved in gynoecium development , was up-regulated at 12 h after 6-BA treatment (Fig. 8). KANADI 2 (KAN2), which encodes a nuclear-localized protein involved in the development of the carpel and the outer integument of the ovule , was also up-regulated at 12 h (Fig. 8). The transcript abundance of two other genes JAGGED (JAG) and Cytochrome P450 78A9 (CYP79A9), which were positively related to gynoecium development [64, 65], were up-regulated at 12 and 24 h (Fig. 8). APETALA 2 (AP2), a member of the AP2/EREBP class of transcription factors, was involved in the specification of floral organ identity. In Arabidopsis, mutation of AP2 would cause the transition of perianth organs to carpels . The transcriptome results showed that the expression of AP2 in S. sebiferum was significantly decreased at 24 h after 6-BA treatment (Fig. 8), suggesting a feminization effect of cytokinin on the floral sex in S. sebiferum. Additionally, several genes, such as MYB Domain Protein 108 (MYB108), NPK1-Activating Kinesin 1 (NACK1), TATA Box Associated Factor II 59 (TAFII59) that were important in the regulation of androecium development [67,68,69], were down-regulated in response to 6-BA treatment (Fig. 8). PISTILLATA (PI) which encodes a MADS domain transcription factor in regulating the specification of the petal and stamen identities , was also down-regulated by 6-BA treatment (Fig. 8), suggesting cytokinin had inhibitory effects on the androecium development in S. sebiferum. These results demonstrated that cytokinin regulated the sex differentiation in S. sebiferum by directly targeting the genes related to the gynoecium and androecium development (Fig. 9).
Conclusively, this study revealed that exogenous application with cytokinin has strong feminization effects on the floral development in S. sebiferum. The transcriptome analysis of the transition from the androecious to monoecious inflorescences after 6-BA treatment on the androecious S. sebiferum contributed to the identification of the candidate genes that were involved in the regulation of sex determination in response to cytokinin. Our study provided a basis for further understanding the molecular mechanisms of the hormonal regulation of sex determination in S. sebiferum.
Differentially expressed genes
Kyoto Encyclopedia of Genes and Genomes
Polymerase chain reaction
Quantitative real time polymerase chain reaction
Reads per kilobase per million mapped reads
Xu I, Chikashige T, Meguro S, Kawachi S. Effective utilization of Stillingia or Chinese tallow tree (Sapium sebiferum) fruits. Mokuzai Gakkaishi. 1991;37(5):494–8.
Peng D, Zhang L, Tan XF, Yuan DY, Liu XM, Zhou B. Increasing seed oil content and altering oil quality of Brassica napus L. by over-expression of diacylglycerol acyltransferase 1 (SsDGAT1) from Sapium sebiferum (L.) Roxb. Mol Breeding. 2016;36(10):136.
Gao Y, Chen W, Lei H, Liu Y, Lin X, Ruan R. Optimization of transesterification conditions for the production of fatty acid methyl ester (FAME) from Chinese tallow kernel oil with surfactant-coated lipase. Biomass Bioenergy. 2009;33(2):277–82.
Golenberg EM, West NW. Hormonal interactions and gene regulation can link monoecy and environmental plasticity to the evolution of dioecy in plants. Am J Bot. 2013;100(6):1022–37.
Orozco-Arroyo G, Vazquez-Santana S, Camacho A, Dubrovsky JG, Cruz-Garcia F. Inception of maleness: auxin contribution to flower masculinization in the dioecious cactus Opuntia stenopetala. Planta. 2012;236(1):225–38.
Chen MS, Pan BZ, Fu Q, Tao YB, Martínezherrera J, Niu L, Ni J, Dong Y, Zhao ML, Xu ZF. Comparative transcriptome analysis between gynoecious and monoecious plants identifies regulatory networks controlling sex determination in Jatropha curcas. Front Plant Sci. 1953;7:2016.
Pan BZ, Chen MS, Ni J, Xu ZF. Transcriptome of the inflorescence meristems of the biofuel plant Jatropha curcas treated with cytokinin. BMC Genomics. 2014;15(1):974.
Yin TJ, Quinn JA. Tests of a mechanistic model of one hormone regulating both sexes in Cucumis sativus (Cucurbitaceae). Am J Bot. 1995;82(12):1537–46.
Pan BZ, Xu ZF. Benzyladenine treatment significantly increases the seed yield of the biofuel plant Jatropha curcas. J Plant Growth Regul. 2011;30(2):166–74.
Fu QT, Niu LJ, Zhang QF, Pan BZ, He HY, Xu ZF. Benzyladenine treatment promotes floral feminization and fruiting in a promising oilseed crop Plukenetia volubilis. Ind Crop Prod. 2014;59(x):295–8.
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(8):1655–66.
Zhang M, Sun H, Fei Z, Zhan F, Gong X, Gao S. Fastq_clean: An optimized pipeline to clean the Illumina sequencing data with quality control. In: IEEE International Conference on Bioinformatics and Biomedicine; 2014. p. 44–8.
Andrews S: FastQC a quality control tool for high throughput sequence data. 2010.
Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, Couger MB, Eccles D, Li B, Lieber M. De novo transcript sequence reconstruction from RNA-seq using the trinity platform for reference generation and analysis. Nature Protocol. 2013;8(8):1494–512.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.
Davidson NM, Oshlack A. Corset: enabling differential gene expression analysis for de novo assembled transcriptomes. Genome Biol. 2014;15:410.
Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.
Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L, et al. WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006;34:293–7.
Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40(1):109–14.
Tatusov RL, Natale DA, Garkavtsev IV, Tatusova TA, Shankavaram UT, Rao BS, Kiryutin B, Galperin MY, Fedorova ND, Koonin EV. The COG database: new developments in phylogenetic classification of proteins from complete genomes. Nucleic Acids Res. 2001;29(1):22–8.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature Biotechnol. 2010;28(5):511–174.
Gao R, Su Z, Yin Y, Sun L, Li S. Germplasm, chemical constituents, biological activities, utilization, and control of Chinese tallow ( Triadica sebifera (L.) small). Biol Invas. 2016;18(3):1–21.
Wahl V, Ponnu J, Schlereth A, Arrivault S, Langenecker T, Franke A, Feil R, Lunn JE, Stitt M, Schmid M. Regulation of flowering by trehalose-6-phosphate signaling in Arabidopsis thaliana. Science. 2013;339(6120):704–7.
Allen M, Qin WS, Moreau F, Moffatt B. Adenine phosphoribosyltransferase isoforms of Arabidopsis and their potential contributions to adenine and cytokinin metabolism. Physiol Plant. 2002;115(1):56–68.
Smehilova M, Dobruskova J, Novak O, Takac T, Galuszka P. Cytokinin-specific glycosyltransferases possess different roles in cytokinin homeostasis maintenance. Front Plant Sci. 2016;7:1264.
To JP, Haberer G, Ferreira FJ, Deruere J, Mason MG, Schaller GE, Alonso JM, Ecker JR, Kieber JJ. Type-a Arabidopsis response regulators are partially redundant negative regulators of cytokinin signaling. Plant Cell. 2004;16(3):658–71.
Cheng Y, Dai X, Zhao Y. Auxin synthesized by the YUCCA flavin monooxygenases is essential for embryogenesis and leaf formation in Arabidopsis. Plant Cell. 2007;19(8):2430–9.
Zhang J, Lin JE, Harris C, Campos Mastrotti Pereira F, Wu F, Blakeslee JJ, Peer WA. DAO1 catalyzes temporal and tissue-specific oxidative inactivation of auxin in Arabidopsis thaliana. Proc Natl Acad Sci. 2016;113(39):11010–5.
Ando S, Sato Y, Kamachi S, Sakai S. Isolation of a MADS-box gene (ERAF17) and correlation of its expression with the induction of formation of female flowers by ethylene in cucumber plants (Cucumis sativus L.). Planta. 2001;213(6):943–52.
Dellaporta SL, Calderonurrea A. The sex determination process in maize. Science. 1994;266(5190):1501–5.
Heslop-Harrison J. Auxin and Sexuality in Cannabis sativa. Physiol Plant. 2010;9(4):588–97.
Zuercher E, Mueller B. Cytokinin synthesis, signaling, and function-advances and new insights. Int Rev Cell Mol Biol. 2016;324:1–38.
Takahashi H, Suge H, Saito T. Sex expression as affected by N 6-benzylaminopurine in staminate inflorescence of Luffa cylindrica. Plant Cell Physiol. 1980;21(4):525–36.
Cohen I, Raz S, Malah D. Kinetin-a multiactive molecule. Int J Biol Macromol. 2007;40(3):182–92.
Fu Q, Niu L, Chen MS, Tao YB, Wang X, He H, Pan BZ, Xu ZF. De novo transcriptome assembly and comparative analysis between male and benzyladenine-induced female inflorescence buds of Plukenetia volubilis. J Plant Physiol. 2018;221:107–18.
Chen MS, Pan BZ, Wang GJ, Ni J, Niu L, Xu ZF. Analysis of the transcriptional responses in inflorescence buds of Jatropha curcas exposed to cytokinin treatment. BMC Plant Biol. 2014;14(1):318.
Raines T, Blakley IC, Tsai Y-C, Worthen JM, Franco-Zorrilla JM, Solano R, Schaller GE, Loraine AE, Kieber JJ. Characterization of the cytokinin-responsive transcriptome in rice. BMC Plant Biol. 2016;16(1):260.
Skoog FMC. Chemical regularion of growth and organ formation in plant fissue cultured. Symp Soc Exp Biol. 1957;11:118–31.
Rioukhamlichi C, Huntley R, Jacqmard A, Murray JA. Cytokinin activation of Arabidopsis cell division through a D-type cyclin. Science. 1999;283(5407):1541–4.
Kosugi S, Ohashi Y. Constitutive E2F expression in tobacco plants exhibits altered cell cycle control and morphological change in a cell type-specific manner. Plant Physiol. 2003;132(4):2012–22.
Paz Sanchez MDL, Gutierrez C. Arabidopsis ORC1 is a PHD-containing H3K4me3 effector that regulates transcription. Proc Natl Acad Sci U S A. 2009;106(6):2065–70.
Zeng CJT, Lee YRJ, Liu B. The WD40 repeat protein NEDD1 functions in microtubule organization during cell division in Arabidopsis thaliana. Plant Cell. 2009;21(4):1129–40.
Ahn CS, Cho HK, Lee DH, Sim HJ, Kim SG, Pai HS. Functional characterization of the ribosome biogenesis factors PES, BOP1, and WDR12 (PeBoW), and mechanisms of defective cell growth and proliferation caused by PeBoW deficiency in Arabidopsis. J Exp Bot. 2016;67(17):5217–32.
Kim YS, Kim SG, Park JE, Park HY, Lim MH, Chua NH, Park CM. A membrane-bound NAC transcription factor regulates cell division in Arabidopsis. Plant Cell. 2006;18(11):3132–44.
Kiyosue T, Abe H, Yamaguchi-Shinozaki K, Shinozaki K. ERD6, a cDNA clone for an early dehydration-induced gene of Arabidopsis, encodes a putative sugar transporter. Bioch Biophysica Act. 1998;1370(2):187–91.
Dyson BC, Allwood JW, Feil R, Yun XU, Miller M, Bowsher CG, Goodacre R, Lunn JE, Johnson GN. Acclimation of metabolism to light in Arabidopsis thaliana: the glucose 6-phosphate/phosphate translocator GPT2 directs metabolic acclimation. Plant Cell Environ. 2015;38(7):1404–17.
Handford M, Rodríguez-Furlán C, Marchant L, Segura M, Gómez D. Arabidopsis thaliana AtUTr7 encodes a golgi-localized UDP-glucose/UDP-galactose transporter that affects lateral root emergence. Mol Plant. 2012;5(6):1263–80.
Nordström A, Tarkowski P, Tarkowska D, Norbaek R, Åstot C, Dolezal K, Sandberg G. Auxin regulation of cytokinin biosynthesis in Arabidopsis thaliana: a factor of potential importance for auxin–cytokinin-regulated development. Proc Natl Acad Sci U S A. 2004;101(21):8039–44.
Jones B, Gunnerås SA, Petersson SV, Tarkowski P, Graham N, May S, Dolezal K, Sandberg G, Ljung K. Cytokinin regulation of auxin synthesis in Arabidopsis involves a homeostatic feedback loop regulated via auxin and cytokinin signal transduction. Plant Cell. 2010;22(9):2956–69.
Yin X-J, Volk S, Ljung K, Mehlmer N, Dolezal K, Ditengou F, Hanano S, Davis SJ, Schmelzer E, Sandberg G, et al. Ubiquitin lysine 63 chain–forming ligases regulate apical dominance in Arabidopsis. Plant Cell. 2007;19(6):1898–911.
Cucinotta M, Manrique S, Guazzotti A, Quadrelli NE, Mendes MA, Benkova E, Colombo L. Cytokinin response factors integrate auxin and cytokinin pathways for female reproductive organ development. Development. 2016;143(23):dev.143545.
Tan M, Xue J, Wang L, Huang J, Fu C, Yan X. Transcriptomic analysis for different sex types of Ricinus communis L. during development from apical buds to inflorescences by digital gene expression profiling. Front Plant Sci. 2016;6:1208.
Rudich J, Halevy AH. Involvement of abscisic acid in the regulation of sex expression in the cucumber. Plant Cell Physiol. 1974;15(4):635–42.
Ko JH, Yang SH, Han KH. Upregulation of an Arabidopsis RING-H2 gene, XERICO, confers drought tolerance through increased abscisic acid biosynthesis. Plant J. 2006;47(3):343–55.
Sussmilch FC, Brodribb TJ, Sam MA. Up-regulation of NCED3 and ABA biosynthesis occur within minutes of a decrease in leaf turgor but AHK1 is not required. J Exp Bot. 2017;68(11):2913–8.
Cowan AK, Cairns ALP, Bartelsrahm B. Regulation of abscisic acid metabolism: towards a metabolic basis for abscisic acid-cytokinin antagonism. J Exp Bot. 1999;50(334):595–603.
Cowan AK, Railton ID. Cytokinins and ancymidol inhibit abscisic acid biosynthesis in Persea gratissima. J Plant Physiol. 1987;130(2–3):273–7.
Elshowk S, Ruonala R, Helariutta Y. Crossing paths: cytokinin signalling and crosstalk. Development. 2013;140(7):1373–83.
Chandler JW, Werr W. DORNRÖSCHEN, DORNRÖSCHEN-LIKE, and PUCHI redundantly control floral meristem identity and organ initiation in Arabidopsis. J Exp Bot. 2017;68(13):3457–72.
Pfent C, Pobursky KJ, Sather DN, Golenberg E. Characterization of SpAPETALA3 and SpPISTILLATA, B class floral identity genes in Spinacia oleracea, and their relationship to sexual dimorphism. Dev Gen Evolut. 2005;215(3):132–42.
Li S-F, Zhang G-J, Zhang X-J, Yuan J-H, Deng C-L, Gao W-J. Comparative transcriptome analysis reveals differentially expressed genes associated with sex expression in garden asparagus (Asparagus officinalis). BMC Plant Biol. 2017;17(1):143.
Schuster C, Gaillochet C, Lohmann JU. Arabidopsis HECATE genes function in phytohormone control during gynoecium development. Development. 2015;142(19):3343–50.
Pires HR, Monfared MM, Shemyakina EA, Fletcher JC. ULTRAPETALA trxG genes interact with KANADI transcription factor genes to regulate Arabidopsis gynoecium patterning. Plant Cell. 2014;26(11):4345–61.
Dinneny JR, Yadegari R, Fischer RL, Yanofsky MF, Weigel D. The role of JAGGED in shaping lateral organs. Development. 2004;131(5):1101–10.
Sotelosilveira M, Cucinotta M, Chauvin AL, Colombo L, De FS. Cytochrome P450 CYP78A9 is involved in Arabidopsis reproductive development. Plant Physiol. 2013;162(2):779–99.
Kunst L, Klenz JE, Martinez-Zapater J, Haughn GW. AP2 gene determines the identity of perianth organs in flowers of Arabidopsis thaliana. Plant Cell. 1989;1(12):1195–208.
Mandaokar A, Browse J. MYB108 acts together with MYB24 to regulate jasmonate-mediated stamen maturation in Arabidopsis. Plant Physiol. 2009;149(2):851–62.
Oh SA, Bourdon V, Pal MD, Dickinson H, Twell D. Arabidopsis kinesins HINKEL and TETRASPORE act redundantly to control cell plate expansion during cytokinesis in the male gametophyte. Mol Plant. 2008;1(5):794–9.
Lago C, Clerici E, Dreni L, Horlow C, Caporali E, Colombo L, Kater MM. The Arabidopsis TFIID factor AtTAF6 controls pollen tube growth. Dev Biol. 2005;285(1):91–100.
Wuest SE, Wellmer F. Molecular basis for the specification of floral organs by APETALA3 and PISTILLATA. Proc Natl Acad Sci U S A. 2012;109(33):13452.
The authors are grateful to Mr. Li for the assistance of the field work. We also thank Dr. Kaiqin Ye for the suggestions for the manuscript.
This work was funded by National Natural Science Foundation of China (31500531), the Science and Technology Service program of Chinese Academy of Sciences (KFJ-STS-ZDTP-002), the key program of 13th five-year plan, CASHIPS (No. kp-2017-21), the Grant of the President Foundation of Hefei Institutes of Physical Science of Chinese Academy of Sciences (YZJJ201619), the major special project of Anhui Province (16030701103) and Natural Science Foundation of Anhui province (1708085QC70).
Availability of data and materials
All the data in the study has been included in the tables and figures in this manuscript. RNA-seq data of this project have been deposited in Genbank under the SRA accession number PRJNA438002.
Ethics approval and consent to participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Photographs of the floral buds used for hormonal treatment and RNA-sequencing. (TIF 5697 kb)
Table S1. Primers for quantitative real-time PCR (DOCX 16 kb)
Figure S2. Photograph of the inflorescences of different genotypes of S. sebiferum. (TIF 1797 kb)
Tables S2-S4. Information of the assembled genes at 6 h, 12 h and 24 h. (XLSX 48555 kb)
Tables S5-S7. Differentially expressed genes at 6 h, 12 h and 24 h after 6-BA treatment. (XLSX 363 kb)
Figure S3. GO analysis of the differentially expressed genes at 6 h (A), 12 h (B) and 24 h (C) after 6-BA and mock treatment. (TIF 1842 kb)
About this article
Cite this article
Ni, J., Shah, F.A., Liu, W. et al. Comparative transcriptome analysis reveals the regulatory networks of cytokinin in promoting the floral feminization in the oil plant Sapium sebiferum. BMC Plant Biol 18, 96 (2018). https://doi.org/10.1186/s12870-018-1314-5
- Floral feminization
- Bioenergy plant
- Sapium sebiferum