Skip to main content

Comparative transcriptomic analysis of the flower induction and development of the Lei bamboo (Phyllostachys violascens)



Bamboo is a very important forest resource. However, the prolonged vegetative stages and uncertainty of flowering brings difficulties in bamboo flowers sampling. Until now, the flowering mechanism of bamboo is still unclear.


In this study, three successive stages of flowering buds and the corresponding vegetative buds (non-flowering stage) from Lei bamboo (Phyllostachys violascens) were collected for transcriptome analysis using Illumina RNA-Seq method. We generated about 442 million clean reads from the above samples, and 132,678 unigenes were acquired with N50 of 1080 bp. A total of 7266 differentially expressed genes (DEGs) were determined. According to expression profile and gene function analysis, some environmental stress responsive and plant hormone-related DEGs were highly expressed in the inflorescence meristem formation stage (TF_1) while some floral organ development related genes were up-regulated significantly in floral organs determination stage (TF_2) and floral organs maturation (TF_3) stage, implying the essential roles of these DEGs in flower induction and maturation of Lei bamboo. Additionally, a total of 25 MADS-box unigenes were identified. Based on the expression profile, B, C/D and E clade genes were more related to floral organs development compared with A clade genes in Lei bamboo.


This transcriptome data presents fundamental information about the genes and pathways involved in flower induction and development of Lei bamboo. Moreover, a critical sampling method is provided which could be benefit for bamboo flowering mechanism study.


In plant kingdom, flowering is necessary for the transition from vegetative stage to reproductive stage as one of the most important process. Investigation into Arabidopsis thaliana reveals that the phenomenon of the flowering is controlled by diverse environmental and endogenous factors, such as temperature, light signals, day length and plant hormones [1,2,3,4]. For example, plant endogenous hormones Gibberellin (GA) could promote flower development and low temperature will cause flowering delay in Arabidopsis [5, 6]. Additionally, the ABCDE model genes are widely used for understanding floral development [7], and the MADS-box genes are well known as ABC model factors to control the floral organ identity during flower development [8]. According to ABC model, B, C and E genes control petal and stamen identity; C and E genes control carpel formation; D and E genes are involved in ovule identity [9,10,11].

Bamboo is one of the most important non-timber forest resources in the world. Its fast growth and strong nitrogen fixation capacity have received much attention in economy and ecology. However, the prolong vegetative period and the death of bamboo after flowering, lead to great economic and ecological losses [12]. The uncertainty of bamboo flowering time increased the difficulty of sampling, further inhibiting the study of flower induction to bamboo. Fortunately, the completion of moso bamboo (Phyllostachys edulis) draft genome sequencing is beneficial to screen the flowering related genes. For example, the transcription factors play an important role in gene expression, thereby influencing the floral transition, the PheMADS14 and PheMADS4 obtained from moso bamboo may play vital regulatory roles in flower development [13,14,15,16]. Lei bamboo (Phyllostachys violascens), widely distributed over south of China, has high economic value because of its bamboo shoots edibility. The income for intensively managed Lei bamboo is about 20 times higher than that of rice [17]. During flowering, however, shoot production decreases sharply. Different from moso bamboo flowering gregariously with 60–120 years period [18], Lei bamboo flowers sporadically [19]. It suggests a similar but not identical flowering mechanism between moso bamboo and Lei bamboo. Recently, some flowering genes in Lei bamboo have also been identified, such as PvFRI-l, PvSOC1 and PvMADS56 [20,21,22]. However, the flowering mechanism in Lei bamboo is still poorly understood.

Next-generation sequencing [23] technologies such as the Illumina Solexa, Roche 454, and ABI SOLID platforms have revolutionized biological research by providing genomic and transcriptome data rapidly and inexpensively [24]. In this study, we used RNA-sequencing (RNA-seq) to investigate the flowering mechanism of Lei bamboo by transcriptome analysis [25]. The goal of this study was to obtain a complete set of assembled unigenes and transcripts for Lei bamboo and to identify flowering-related genes. Our data will be useful for Lei bamboo flowering study and provides a critical sampling method to study flowering mechanism in bamboos.


Plant materials

After observing flowering characteristics of Lei bamboo for more than 10 years, we began to learn and pre-judge before flowering occurs with a certain probability [26]. Samples of Lei bamboo were collected from bamboo forests cultivated in Lin’an, Zhejiang province, China. Flowering buds and vegetative buds were collected at about 2 pm from expected flowering and non-flowering plants from the same rhizome, at different stages: March 8th and 29th and April 12th, 2013. The expected bud samples confirmed after bamboo flowering were used for transcriptome sequencing. Flowering buds from a flowering bamboo plant were sampled at three stages, formation of inflorescence meristem (TF_1), formation and identity determination of floral organs (TF_2), and maturation of floral organs (TF_3) stages. Vegetative buds (TV) from a non-flowering bamboo plant growing on the same rhizome of the flowering bamboo plant were both sampled at the same time (Fig. 1g). The first two stages (TF_1 and TV_1) were sampled on March 8th, 2013, when the inflorescence meristem of bamboo began to form, but the meristem morphology of flower buds was similar to that of vegetative buds (Fig. 1a, Fig. 1d). The second two stages (TF_2 and TV_2) were sampled on March 29th, 2013, when the inflorescence appeared, and the floral organs started to take shape (Fig. 1b, Fig. 1e). The third two stages (TF_3 and TV_3) were sampled on April 12th, 2013, when the floral organs continued to develop and gradually become mature (Fig. 1c, Fig. 1f). The six samples were used for cDNA library construction and RNA-seq separately. All samples were transferred to liquid nitrogen immediately and stored in − 80 °C.

Fig. 1
figure 1

The reproductive and vegetative bud samples of Lei bamboo. a-c The flowering buds sampled on March 8th, March 29th and April 12th are named as TF_1, TF_2 and TF_3 respectively. d-f Vegetative buds sampled on March 8th, March 29th and April 12th are named TV_1, TV_2 and TV_3 respectively. g The same rhizome of flowering and non-flowering plants

RNA-seq sequencing and transcriptome assembly

We extracted total RNA from samples by using TRIzol reagent (Invitrogen) and treated with RNase-free DNase I (New England Biolabs) to remove any contaminating genomic DNA. Then cDNA libraries were constructed by using NEBNext® Ultra™ Directional RNA Library Prep Kit for Illumina® (NEB, USA) instruction. Raw data (raw reads) in the fastq format were first processed by using in-house perl scripts. To ensure high data quality, we acquired the clean data by removing reads containing an adapter or poly-N and low-quality reads, then the clean-data quality including Q20, Q30, GC-content and error rate were calculated (Table 1). Finally, the clean data was assembled by using Trinity software with min_kmer_cor set to 2 by default; all other parameters were set to the default [27].

Table 1 The analysis of data output quality. Sample: Sample name_1, left reads; Sample name_2, right reads. The total number of clean reads is left + right. Q20: the percentage of bases with a phred value > 20. Q30: the percentage of bases with a phred value > 30. GC content: the GC ratio of the total base number

Analysis of differential expressed genes

Expression level of genes is calculated by using RPKM method [28, 29]. For each sequenced library, the read counts were adjusted by using the edgeR program package with one scaling normalized factor. The DEGseq R package was used for differential expression analyses. Qvalue < 0.05 & |log2 (fold change)| > 1 were set as the thresholds for significantly differential expression. Gene function was annotated by using the following databases: Nr (NCBI non-redundant protein sequences), Nt (NCBI non-redundant nucleotide sequences), Pfam (Protein family), COG (Clusters of Orthologous Groups of proteins), Swiss-Port (A manually annotated and reviewed protein sequence database), KO (KEGG Ortholog database) and GO (Gene Ontology) [30,31,32]. We used NCBI blast 2.2.28+to search genes in the Nr, Nt and Swiss-Port databases with E-value 1e-5, and the COG database with E-value 1e-3 [33]. We predicted the protein domain from Pfam using HMMER 3.0 software and hmmscan with E-value 1e-5 [34]. Blast2GO v2.5 was used for GO function annotation based on the annotation of the NR and Pfam databases [35]. GO enrichment analysis of differentially expressed genes (DEGs) involved in the GOseq R packages based on Wallenius non-central hyper-geometric distribution, when the corrected P-value < 0.05, the function is regarded as an enrichment item [36]. We used KOBAS software to test the statistical enrichment of DEGs in KEGG pathways with E-value < 10− 5 and rank ≤5 [37].

Identification and phylogenetic analysis of the MADS-box gene family

The Arabidopsis MADS protein sequences were acquired from TAIR ( database. We used Arabidopsis MADS domain sequences as query to blast against the Lei bamboo transcriptome database. Multiple sequence alignment was performed using the full-length protein sequences in the program Clustal W. The un-rooted neighbor-joining (NJ) [23] and maximum likelihood (ML) [38] trees were constructed by using the MEGA6 package [39] with bootstrap support of 1000 replicates.

qRT-PCR analysis

Several DEGs involved in flower development of Lei bamboo were chosen for validation with quantitative real-time PCR (qRT-PCR). qRT-PCR reaction was analyzed with the CFX96TM real-time PCR detection system (Bio-Rad, Hercules, USA) with the SYBR Green Master Mix (Takara, Dalian, China), and amplified with 1 ul of cDNA template, 10 ul of 2 x SYBR Green Master Mix, and 1ul of each primer (10 umol/ul), to a final volume of 20 ul by adding water. The amplification program consisted of 1 cycle of 95 °C for 30s, followed by 39 cycles of 95 °C for 5 s, TM for 20s, 72 °C for 20s. The Actin was used as reference gene. All the qPCR assays were performed with three biological and four technical replicates, and the quantitative analysis used the 2-ΔΔCT method.


Identification of total unigenes during Lei bamboo flower development

A total of 468,502,728 raw reads were generated by Illumina sequencing of TF and TV, resulting in 442,717,036 clean reads after moving low quality sequences (Table 1). Assembly of the clean reads produced 894,499,502 unigenes with sequence length range of 201–16,212 bp. The N50 length of unigenes was 1080 bp and the N90 length was 269 bp (Table 2).

Table 2 The length distribution analysis

In this study, we analyzed the function of the total unigenes by annotating the unigenes with seven databases (Additional file 1). A total of 2931(2.2%) unigenes were annotated in all databases and 58,628 (44.18%) unigenes were annotated in at least one database. According to the GO analysis (Fig. 2a), the genes associated with “cellular process”, “metabolic process”, “single-organism process”, “biological regulation” were enriched. To identify biochemical pathways, we mapped the annotated sequences in the KEGG database (Fig. 2b). Most unigenes were involved in “translation”, “energy metabolism” and “carbohydrate metabolism”, with the least in “membrane transport”, “signaling molecules” and “degradation”. These pathway assignments provide valuable information for investigating specific biochemical and development processes.

Fig. 2
figure 2

The functional analysis of total unigenes in Lei bamboo. a GO functional classification (b) KEGG classification

The phytohormone related and stress responsive DEGs are involved in Lei bamboo flower induction

A total of 7266 unigenes were identified as DEGs from 16 comparisons (TF_1 vs TF_2, TF_1 vs TF_3, TF_2 vs TF_3, TF_23 vs TF_1, TF13 vs TF_2, TF12 vs TF_3, TV_2 vs TV_1, TV_3 vs TV_1, TV_3 vs TV_2, TV23 vs TV_1, TV_13 vs TV_2, TV_12 vs TV_3, TF_1 vs TV_1, TF_2 vs TV_2, TF_3 vs TV_3, TF vs TV). According to the expression pattern (Additional file 2), more than half of the DEGs were highly expressed at TF_3. In the rest DEGs, most were highly expressed in TV stages, but only a few genes have higher expressions in former stages of flower development (TF_1 and TF_2). In order to better to analyze the flowering mechanism of Lei bamboo, we categorized the samples according to their development stages and analyzed the DEGs involved in TV and TF stages (Fig. 3a and b). For TV stages, 388, 837 and 160 DEGs were identified in the comparisons of TV_3 vs TV_1, TV_2 vs TV_1 and TV_3 vs TV_2, while for TF stages, 127, 1856 and 1757 DEGs were found in the comparisons of TF_2 vs TF_1,TF_3 vs TF_1 and TF_2 vs TF_3. After removing the duplicated ones separately, a total of 958 DEGs in vegetative buds and 2224 DEGs in flowering buds were finally identified.

Fig. 3
figure 3

The venn diagram of differentially expressed genes (DEGs) in Lei bamboo. a the DEGs in TV stages (b) the DEGs in TF stages (c) the DEGs between TF and TV stages

As described before, the samples of TF_1 and TV_1 were collected at the same time when the inflorescence meristem began to form, but there was no flower meristem morphology in TF_1 (Fig. 1). It indicates the TF_1 might be flower induction stage in Lei bamboo. Although the two periods of TF_1 and TV_1 differed in flower induction, they experienced the same biological processes, such as vegetative buds and stems development. Based on KEGG analysis (Additional file 3), the DEGs in TF_1 (compared with TF_2 and TF_3) were enriched in a number of same biological pathways as well as in TV_1 (compared with TV_2 and TV_3) such as the “photosynthesis”, “cell cycle”, “metabolic” and “DNA replication” pathways. To better to obtain the flower induction related genes, we removed the coincident genes in TF_1 compared to TV_1. Finally, a total of 395 genes (28 from the comparison of TF_1 vs TF_2, and 367 from the comparison of TF_1 vs TF_3) were differential expressed significantly specific in TF_1. According to KEGG analysis, the “plant hormone signal transduction” pathway was abundant in TF_1 compared with corresponding TV_1 (Additional file 3) and a serious of highly-expressed related genes (mainly includes ARF, CRE1, BRI1, BSK, NPR1, AUX1, TGA, TIR1, B-ARR and CTR1) were identified. The expression profile suggests these plant hormones related DEGs should be involved in the regulating of Lei bamboo flower induction (Fig. 4). Additionally, we found the “detection of abiotic stimulus”, “detection of external stimulus”, “detection of light stimulus” related genes are up-regulated but “long-day photoperiodism, flowering” related genes are down-regulated in TF_1 (compared with TV_1), respectively (Fig. 4b). It reflects that the flower induction of Lei bamboo possibly affected by environmental stress factors. Therefore, the genes involved in stress responsive and phytohormone should be essential for the process of Lei bamboo flower induction.

Fig. 4
figure 4

The heatmap of the selected DEGs in Lei bamboo flowering (a). the genes involved in plant hormone signaling pathway (b) the genes involved in photoperiod and stress resistance

Floral organ development related DEGs are involved in Lei bamboo flower maturation

With the development of Lei bamboo flowers (TF_2), the “lipid metabolic and biosynthetic process”, “oxidoreductase activity” related genes were up-regulated (compared with TF_1). They were involved in the “flavonoid biosynthesis”, “circadian rhythm”, “biosynthesis of secondary metabolites” and “metabolic” pathways (Additional file 3 and Additional file 4), indicating their essential roles in formation of floral organ primordia. Moreover, a total of 1149 genes were up-regulated in TF_3 compared with TF_1 and TF_2. In this stage, the floral organ began to mature. GO analysis showed that the “single-organism metabolic process”, “cytoplasmic part”, “catalytic activity” related genes were abundant. They were enriched in a large number of metabolic processes of “Ribosome”, “beta-Alanine metabolism”, “Ascorbate and aldarate metabolism” (Additional file 3 and Additional file 4). It reflects the maturation of flower organs, especially the reproductive organs of carpels and stamens, requires a great deal of material and energy exchange in the TF_3 stage of Lei bamboo.

To better determining the critical genes during Lei bamboo flower organ development, a comparative analysis was made between two later reproductive stages (TF_2 and TF_3) with the corresponding vegetative stages (TV_2 and TV_3), and a total of 293 and 2419 DEGs were get. Among the DEGs, 1377 genes were found to be up-regulated during floral organs identity determination (TF_2) and floral organs maturation (TF_3) stages, that were enriched in the functions of “oxidation-reduction process”, “single-organism metabolic process”, “cytoplasmic part” and “catalytic activity”. They were abundant in “beta-Alanine metabolism”, “Ribosome”, “Ascorbate and aldarate metabolism”, “Amino sugar and nucleotide sugar metabolism” pathways. Additionally, some genes involved in “female gamete generation”, “flower development”, “pigment biosynthetic process”, “anther morphogenesis” were found to differential expressed significantly in the latest stage (TF_3), such as comp135901_c0, comp128343_c0, comp136560_c0, comp136560_c0, comp136560_c0 (Fig. 5), reflecting their critical roles in flower organ formation and maturation.

Fig. 5
figure 5

The heatmap of the selected DEGs involved in floral organ development of Lei bamboo

The B, C/D and E MADS genes are involved in flower organ development of Lei bamboo

From the conserved 103 Arabidopsis MADS proteins in the hidden Markov model [40], we acquired 25 Lei bamboo transcriptome ORF sequences designed as MADS genes. The phylogenetic trees were constructed by using the NJ (Fig. 6a) and ML (Fig. 6b) methods based the MADS protein classification of Arabidopsis [41]. In the NJ and ML trees, a bit interior branches exist some minor differences, but in general, the classification results are very similar. It suggests that most proteins have similar origin in both trees. Here, we used NJ phylogenetic tree for further analysis. Accordingly, 11 unigenes were found to belong to the Mα, Mβ and Mγ, two were grouped in the Mδ type and other 12 were in the MIKC type.

Fig. 6
figure 6

Phylogenetic analysis of MADS-box proteins in Arabidopsis and Lei bamboo. a A total of 25 open reading frame sequences in Lei bamboo and 103 proteins in Arabidopsis was used to construct the neighbor-joining Tree. b A total of 25 open reading frame sequences in Lei bamboo and 103 proteins in Arabidopsis was used to construct the maximum likelihood tree

According to the phylogenetic tree (Fig. 6a), we found eight unigenes belonging to ABCDE model genes. A class genes represented by APETALA1(AP1) and APETALA2(AP2) in Arabidopsis [42], and three unigenes (comp128538_c2_orf1, comp133968_c1_orf1, comp127325_c0_orf1) were in a same group with them. PI(PISTILLATA) and AP3(APETALA3) belong to B class genes in Arabidopsis [43], while there are two unigenes (comp132107_c0_orf1, comp99121_c0_orf1) belonging to B class genes. In Arabidopsis, AG(AGAMOUS) is recognized as C genes, STK is recognized as D genes, they have similar expression patterns [44,45,46], but we identified only one gene (comp124536_c0) belonging to the C/D genes family. The representative gene of the E class gene is SEP(SEPALLATA)1/2/3/4 [47], and two unigenes (comp122881_c0_orf1, comp131728_c1_orf1) were found to belong to E class genes. According to RNA-seq results, eight genes showed different patterns of expression in different samples (Fig. 7). In these, comp99121_c0, comp128538_c2 and comp133968_c1 were highly expressed in TV stages while the rest four genes (comp132107_c0, comp124536_c0, comp122881_c0, comp131728_c1) that belong to B, C/D and E clades were up-regulated significantly in TF stages, suggesting the essentials of B, C/D and E genes in flower organ development of Lei bamboo. Moreover, the comp119979_c0_orf1 and comp127580_co_orf1 were identified to belong to SVP (SHORT VEGATATIVE PHASE), while comp131813_co_orf1 was belonged to SOC1 (SUPPRESSOR OF OVEREXPRESSION OF CO1), indicating their functional roles as flower regulators in plants [48,49,50].

Fig. 7
figure 7

The heatmap of the MADs genes belong to ABCDE model in Lei bamboo

Verification of the gene expression by qRT-PCR

Transcriptional regulation revealed by RNA-seq was confirmed in a biologically independent experiment using RT-PCR, a total of 12 genes were chosen to design gene-specific primers (Additional file 5). Finally, 10 genes (comp104979_c0, comp139179_c0, comp111894_c0, comp130563_c0, comp139255_c0, comp136544_c2, comp139166_c1, comp139643_c0, comp132107_c0 and comp128073_c0) showed significant correlations (P = 0.05 & P = 0.01) between the qRT-PCR data and the RNA-seq results (Fig. 8), suggesting a good reproducibility between transcript abundance assayed by RNA-seq and the expression profile revealed by qRT-PCR.

Fig. 8
figure 8

Quantitative real-time RT-PCR confirmation of 12 candidate genes at the six stages

Discussion and conclusion

As the randomness and uncertainty of bamboo flowering, how to obtain the accurate flower samples is the most crucial problem for bamboo flowering mechanism study. Here, the sampling time of bamboo flowering materials was earlier than that in most usual researches. For example, the moso bamboo flower buds were took when the flower morphology has appeared [51], but we pre-judged and began to collect the samples in earlier stages (TF_1 and TV_1) when the floral or leaf buds have no obvious different phenotype. This method makes it better to identify the flower induction related genes in bamboo. Moreover, our sampling method is “one-to-one” pattern. As the development of vegetative buds and stems may accompany the flowering process of bamboo, the flower development samples and corresponding vegetative samples were “co-sampled” and “co-sequenced” in our study to better ensure the same genetic background (Fig. 1g). It is not only beneficial to avoid interference from vegetative growth related genes in bamboo flowering mechanism study, but also helpful to identify the related genes involved in vegetative periods in Lei bamboo. Above all, our research will provide a new sampling method for studying the flowering mechanism of bamboo.

We obtained 132,678 assembled unigenes using RNA-seq, with 7266 DEGs identified between flowering and vegetative buds at the three stages of Lei bamboo. Although the RNA-seq databases of certain bamboo species, such as moso bamboo, Ma bamboo (Dendrocalamus latiflorus) and Bambusa edulis [52,53,54], have been established, we provided transcriptome data in Lei bamboo for the first time. Finally, a number of DEGs were determined according to the expression profiles, indicating their functional roles in early flowering initiation and flower development of Lei bamboo. For example, in the TF_1 stage, the “detection of abiotic stimulus”, “detection of external stimulus”, “detection of light stimulus” genes were found to be enriched. But the “photoperiod” genes were downregulated compared with TV_1. Interestingly, the stress-responsive genes were also identified to be highly-expressed in moso bamboo, but the genes employed in typical flowering promotion pathways (such as those in the photoperiod, gibberellins, ambient-temperature or light-quality pathways) and floral pathway integrator (FPI) genes were expressed lowly in these floral tissues [52]. It is highly likely that the flowering of bamboo associated more with the environmental stress factors than photoperiod or other traditional pathways. Plant hormones play necessary role during flower development [47–51]. In Lei bamboo, the “plant hormone signaling transduction” related genes were abundant in TF_1 compared with TV_1. In previous study, the endogenous hormone concentrations were found to be changed with the development of flower bud morphological differentiation [55]. It seems that the auxin (IAA) and Cytokinine (CTK) are act as promoters during flower bud differentiation [56]. As we known, AUX/LAX genes encode auxin influx transporters that have distinct developmental functions and regulatory mechanisms in Arabidopsis, while CTK regulates the activity of reproductive meristems, flower organ size, ovule formation, and thus seed yield [57]. Here, the DEGs CRE1 (comp138523_c0) and ARR-B (comp139166_c1 and comp134657_c0) enriched in the CTK pathway and AUX1/LAX (comp138485_c1, comp137757_c1), TIR1(comp136544_c2) and ARF (comp139738_c1) enriched in the auxin pathway were all significantly upregulated in TF_1 stage. It implies the essential roles of these CTK and IAA related genes in Lei bamboo flower induction.

In TF_3, some “female gamete generation”, “flower development”, “pigment biosynthetic process”, “anther morphogenesis” related genes were up-regulated significantly compared with vegetative stage (TV_3). For higher plants, the anther and carpel developments were recognized as the most important stages of organogenesis regarding sexual reproduction processes [46]. Therefore, these DEGs might have essential roles in the sextual reproduction process of Lei bamboo. The family of MADS-box genes play important roles in plant flower development. They are expressed regionally specific in particular flower tissues, during differentiation growth of the flower development. According to previous reports, MADS-box genes participate in floral organ development and flowering time regulation [58]. The MIKC-type proteins represent a class of MADS-box transcription factors, and the MIKC genes could control flowering time and flower organ identity, especially ABC model genes [59, 60]. In rice, 34 MADS-box genes have been identified, and 15 of which regulate flower development [61]. In moso bamboo, 34 MADS-box genes were identified and the expression patterns of most genes were found to similar to Arabidopsis and rice [23]. In our study, we found eight MADS box genes belonging to ABC model genes. Combing with the expression profiles, the B, C/D and E genes should be associated more with flower development than A genes. In Arabidopsis, B genes(PI and AP3) could control petal identity with A class genes, and together with C class genes form the stamen [47]. We found comp132107_c0 have high expression especially in maturation process of floral organs, hence we infer that comp132107_c0 might mainly control petal and stamen identity in later floral stages. Additionally, it is generally accepted that the C class genes (AG) control the development of stamen and pistil [45], and the D class genes are involved in ovule identity [11]. But in petunia, D genes are not essential to confer ovule identity, instead, this function is redundantly shared among all AG members, and D genes might participate in floral determinacy [62]. In our study, comp124536_c0 had high expression in maturation of floral organs, so it might be a key gene involved in floral determinacy in later development stages. The E genes were known to participate in regulation of flower organ characteristics [63]. SEP 1/2/3/4 are key genes in floral organ development, and their functions are redundant in floral identity. These genes play fundamental roles by interacting with other MADS-box genes products [47, 64, 65]. Two genes comp122881_c0 and comp131728_c1 had high expressions in flower plants and low expressions in non-flower plants, so they might exhibit similar functions as the SEP1/2/3/4 gene in Lei bamboo. The ABCDE class MADS-domain transcription factors can assemble into an organ-specific quaternary protein complex to control downstream genes [58]. Whether this interaction is suitable for Lei bamboo needs to be further explored.

Availability of data and materials

All data generated or analyzed during this study are included in this published article。.



Differentially expressed genes


Gene ontology


Kyoto encyclopedia of genes and genomes


KEGG ortholog database


Clusters of orthologous groups of proteins


NCBI non-redundant protein sequences


NCBI non-redundant nucleotide sequences


Protein family


Reads Per Kilobase per Million mapped reads


quantitative real-time PCR


A manually annotated and reviewed protein sequence database


Flowering buds


Flowering buds of formation of inflorescence meristem


Flowering buds of formation and identity determination of floral organs


Flowering buds of maturation of floral organs


Vegetative buds


Vegetative buds from a non-flowering bamboo plant growing on the same rhizome of the flowering bamboo plant were sampled at the same time with TF_1


Vegetative buds from a non-flowering bamboo plant growing on the same rhizome of the flowering bamboo plant were sampled at the same time with TF_2


Vegetative buds from a non-flowering bamboo plant growing on the same rhizome of the flowering bamboo plant were sampled at the same time with TF_3


  1. Balasubramanian S, Sureshkumar S, Lempe J, et al. Potent induction of Arabidopsis thaliana flowering by elevated growth temperature. PLoS Genet. 2006;2(7):e106.

  2. Sawa MA, Nusinow DA, Kay SA, et al. FKF1 and GIGANTEA complex formation is required for day-length measurement in Arabidopsis. Science. 2007;318(5848):261.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Blázquez MA, Weigel D. Integration of floral inductive signals in Arabidopsis. Nature. 2000;404(6780):889–92.

    Article  PubMed  Google Scholar 

  4. Morris K, Thomber S, Codrai L, et al. DAY NEUTRAL FLOWERING represses CONSTANS to prevent Arabidopsis flowering early in short days. Plant Cell. 2010;22(4):1118–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Balasubramanian S, Weigel D. Temperature induced flowering in Arabidopsis thaliana. Plant Signal Behav. 2006;1(5):227–8.

  6. Mouradov A, Cremer F, Coupland G. Control of flowering time: interacting pathways as a basis for diversity. Plant Cell. 2002;14(Suppl):S111.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Bowman JL, Smyth DR, Meyerowitz EM. The ABC model of flower development: then and now. Development. 2012;139(22):4095.

    Article  CAS  PubMed  Google Scholar 

  8. Schwarzsommer Z, Huijser P, Nacken W, et al. Genetic control of flower development by homeotic genes in Antirrhinum majus. Science. 1990;250(4983):931–6.

    Article  CAS  Google Scholar 

  9. Honma T, Goto K. Complexes of MADS-box proteins are sufficient to convert leaves into floral organs. Nature. 2001;409(6819):525–9.

    Article  CAS  PubMed  Google Scholar 

  10. Theißen G, Saedler H. Floral quartets. Nature. 2001;409:469-71.

  11. Kramer EM, Jaramillo MA, Di SV. Patterns of gene duplication and functional evolution during the diversification of the AGAMOUS subfamily of MADS box genes in angiosperms. Genetics. 2004;166(2):1011–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Zheng Z, Yang X, Fu Y, et al. Overexpression ofPvPin1, a Bamboo Homolog ofPIN1-Type Parvulin 1, Delays Flowering Time in TransgenicArabidopsisand Rice. Front Plant Sci. 2017;8:1526.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Ge W, Zhang Y, Cheng Z, et al. Main regulatory pathways, key genes and microRNAs involved in flower formation and development of moso bamboo (Phyllostachys edulis). Plant Biotechnol J. 2017;15(1):82–96.

    Article  CAS  PubMed  Google Scholar 

  14. Shen Z, Bao W, Huang D-S. Recurrent neural network for predicting transcription factor binding sites. Sci Rep. 2018;8(1):15270.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Guo WL, Huang DS. An efficient method to transcription factor binding sites imputation via simultaneous completion of multiple matrices with positional consistency. Mol BioSyst. 2017;13(9):1827–37.

    Article  CAS  PubMed  Google Scholar 

  16. Zhang Y, Tang D, Lin X, et al. Genome-wide identification of MADS-box family genes in moso bamboo (Phyllostachys edulis) and a functional analysis of PeMADS5 in flowering. BMC Plant Biol. 2018;18(1):176.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Song X, Zhou G, Jiang H, et al. Carbon sequestration by Chinese bamboo forests and their ecological benefits: assessment of potential, problems, and future challenges. Environ Rev. 2011;19:418–28.

    Article  CAS  Google Scholar 

  18. John CK, Nadgauda RS. Bamboo flowering and famine. Curr Sci. 2002;82(3):261–2.

    Google Scholar 

  19. Lin XC, Yuan XL, Lin R, et al. Studies on floral biology of Phyllostachys violascens (in Chinese). J Fujian Coll Forestry. 2010;4(4):333–7.

    CAS  Google Scholar 

  20. Liu SN, Zhu LF, Lin XC, et al. Overexpression of the repressor gene PvFRI-L from Phyllostachys violascens delays flowering time in transgenic Arabidopsis thaliana. Biol Plant. 2016;60(3):1–9.

  21. Liu S, Ma T, Ma L, et al. Ectopic expression of PvSOC1, a homolog of SOC1 from Phyllostachys violascens, promotes flowering in Arabidopsis and rice. Acta Physiol Plant. 2016;38(7):166.

    Article  CAS  Google Scholar 

  22. Liu S, Qi TT, Ma JJ, et al. Ectopic expression of a SOC1 homolog from Phyllostachys violascens alters flowering time and identity of floral organs in Arabidopsis thaliana. Trees. 2016;30(6):1–13.

  23. Cheng Z, Ge W, Li L, et al. Analysis of MADS-box gene family reveals conservation in floral organ ABCDE model of Moso bamboo (Phyllostachys edulis). Front Plant Sci. 2017;8:656.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Metzker ML. Sequencing technologies - the next generation. Nat Rev Genet. 2010;11(1):31–46.

    Article  CAS  PubMed  Google Scholar 

  25. Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Lin X, Yuan X, Lin R, et al. Morphogenesis of indefinite inflorescence of Phyllostachys violascens (Carr.) Riv. J Fujian Coll Forestry. 2012;32(2):141-45.

  27. Grabherr M, Haas BJ, Yassour M, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12(1):323.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Mortazavi A, Williams BA, Mccue K, et al. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8.

    Article  CAS  PubMed  Google Scholar 

  30. Stefan GT, Juan Miguel GG, Javier T, et al. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008;36(10):3420–35.

    Article  CAS  Google Scholar 

  31. Finn RD, Alex B, Jody C, et al. Pfam: the protein families database. Nucleic Acids Res. 2014;42(Database issue):222–30.

    Article  CAS  Google Scholar 

  32. Altschul SF, Madden TL, Schäffer AA, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. 1997.

    Google Scholar 

  33. Altschul SF, Gish W, Miller W, et al. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.

    Article  CAS  PubMed  Google Scholar 

  34. Finn RD, Clements J, Eddy SR. HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011;39:29–37.

    Article  CAS  Google Scholar 

  35. Conesa A, Gotz S, Garciagomez JM, et al. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.

    Article  CAS  PubMed  Google Scholar 

  36. Young MD, Wakefield MJ, Smyth GK, et al. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):1–12.

    Article  CAS  Google Scholar 

  37. Xie C, Mao X, Huang J, et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39:316–22.

    Article  CAS  Google Scholar 

  38. Chin AR, Fong MY, Somlo G, et al. Cross-kingdom inhibition of breast cancer growth by plant miR159. Cell Res. 2016;26(2):217.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Kumar S, Nei M, Dudley JT, et al. MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform. 2008;9(4):299–306.

    Article  CAS  PubMed  Google Scholar 

  40. Eddy SR. Profile hidden Markov models. Bioinformatics. 1998;14(9):755–63.

    Article  CAS  PubMed  Google Scholar 

  41. Parenicova L, De Folter S, Kieffer M, et al. Molecular and phylogenetic analyses of the complete MADS-box transcription factor family in Arabidopsis new openings to the MADS world. Plant Cell. 2003;15(7):1538–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Pařenicová L, Falasca G, Masiero S, et al. Functional characterization of OsMADS18, a member of the AP1/SQUA subfamily of MADS box genes. Plant Physiol. 2004;135(4):2207–19.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Melzer R, Theißen G. Reconstitution of ‘floral quartets’ in vitro involving class B and class E floral homeotic proteins. Nucleic Acids Res. 2009;37(8):2723–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Paolacci AR, Tanzarella OA, Porceddu E, et al. Molecular and phylogenetic analysis of MADS-box genes of MIKC type and chromosome location of SEP-like genes in wheat (Triticum aestivum L.). Mol Gen Genomics. 2007;278(6):689.

    Article  CAS  Google Scholar 

  45. Sablowski R. Flowering and determinacy in Arabidopsis. J Exp Bot. 2007;58(5):899–907.

    Article  CAS  PubMed  Google Scholar 

  46. Pinyopich A, Ditta GS, Savidge B, et al. Assessing the redundancy of MADS-box genes during carpel and ovule development. Nature. 2003;424(6944):85–8.

    Article  CAS  PubMed  Google Scholar 

  47. Pelaz S, Ditta GS, Baumann E, et al. B and C floral organ identity functions require SEPALLATA MADS-box genes. Nature. 2000;405(6783):200–3.

    Article  CAS  PubMed  Google Scholar 

  48. Immink RGH, Pose D, Ferrario S, et al. Characterization of SOC1's central role in flowering by the identification of its upstream and downstream regulators. Plant Physiol. 2012;160(1):433–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Lee J, Lee I. Regulation and function of SOC1, a flowering pathway integrator. J Exp Bot. 2010;61(9):2247–54.

    Article  CAS  PubMed  Google Scholar 

  50. Lee JH, Yoo SJ, Park SH, et al. Role of SVP in the control of flowering time by ambient temperature in Arabidopsis. Genes Dev. 2007;21(4):397–402.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Gao J, Zhang Y, Zhang C, et al. Characterization of the floral Transcriptome of Moso bamboo (Phyllostachys edulis) at different flowering developmental stages by Transcriptome sequencing and RNA-Seq analysis. PLoS One. 2014;9(6):e98910.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Peng Z, Lu Y, Li L, et al. The draft genome of the fast-growing non-timber forest species moso bamboo (Phyllostachys heterocycla). Nat Genet. 2013;45(4):456–61.

    Article  CAS  PubMed  Google Scholar 

  53. Zhang XM, Zhao L, Larsonrabin Z, et al. De novo sequencing and characterization of the floral transcriptome of Dendrocalamus latiflorus (Poaceae: Bambusoideae). PLoS One. 2012;7(8):e42082.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Shih MC, Chou ML, Yue JJ, et al. BeMADS1 is a key to delivery MADSs into nucleus in reproductive tissues-De novo characterization of Bambusa edulis transcriptome and study of MADS genes in bamboo floral development. BMC Plant Biol. 2014;14(1):179.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Ya-Ting LU, Yuan XL, Lin XC, et al. Endogenous hormone changes during floral bud morphological differentiation of Phyllostachys violascens. J Zhejiang A F Univ. 2012;29(2):161–5.

    Google Scholar 

  56. Ding X. Hormonal regulation mechanism of flowering of Phyllostachys violascens under mulching cultivation. Sci Silvae Sin. 2007;43(07):10–5.

    CAS  Google Scholar 

  57. Bartrina I, Otto E, Strnad M, et al. Cytokinin regulates the activity of reproductive meristems, flower organ size, ovule formation, and thus seed yield in Arabidopsis thaliana. Plant Cell. 2011;23(1):69–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Theissen G, Becker A, Di RA, et al. A short history of MADS-box genes in plants. Plant Mol Biol. 2000;42(1):115–49.

    Article  CAS  PubMed  Google Scholar 

  59. Becker A, Theissen G. The major clades of MADS-box genes and their role in the development and evolution of flowering plants. Mol Phylogenet Evol. 2003;29(3):464–89.

    Article  CAS  PubMed  Google Scholar 

  60. Coen E, Meyerowitz EM. The war of the whorls: genetic interactions controlling flower development. Nature. 1991;353(6339):31–7.

    Article  CAS  PubMed  Google Scholar 

  61. Lee S, Kim J, Son JS, et al. Systematic reverse genetic screening of T-DNA tagged genes in rice for functional genomic analyses: MADS-box genes as a test case. Plant Cell Physiol. 2003;44(12):1403.

    Article  CAS  PubMed  Google Scholar 

  62. Heijmans K, Ament K, Rijpkema AS, et al. Redefining C and D in the petunia ABC. Plant Cell. 2012;24(6):2305.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Theissen G. Development of floral organ identity: stories from the MADS house. Curr Opin Plant Biol. 2001;4(1):75–85.

    Article  CAS  PubMed  Google Scholar 

  64. Robles P. The SEP4 gene of Arabidopsis thaliana functions in floral organ and meristem identity. Curr Biol. 2004;14(21):1935–40.

    Article  PubMed  CAS  Google Scholar 

  65. Theißen G, Melzer R. Combinatorial Control of Floral Organ Identity by MADS-domain Transcription Factors: Wiley-Blackwell; 2007. p. 253–65.

    Google Scholar 

Download references


We are thankful to Dr. Dan Hou in College of Forestry and Biotechnology, Zhejiang Agricultural and Forestry University for article revising.

About this supplement

This article has been published as part of BMC Bioinformatics Volume 20 Supplement 25, 2019: Proceedings of the 2018 International Conference on Intelligent Computing (ICIC 2018) and Intelligent Computing and Biomedical Informatics (ICBI) 2018 conference: bioinformatics. The full contents of the supplement are available online at


Publication costs are funded by National Natural Science Foundation of China (31971735), the Natural Science Foundation of Zhejiang Province (LZ20C160002) and the State Key Laboratory of Subtropical Silviculture (ZY20180203).

Author information

Authors and Affiliations



XLin and XLi designed the experiments. YJ and QH analyzed data and wrote the manuscript. YZ, LZ, TM, HZ, QZ made substantial contributions to the materials collection and data analysis. XLin revised and edited the manuscript. All authors read and approved the paper.

Corresponding authors

Correspondence to Xuan Li or Xinchun Lin.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1.

Summary of the annotations for the assembled unigenes in public databases.

Additional file 2.

Cluster analysis of differentially expressed genes. (a) Heat map of DEGs expression. TF_1: the expression value in TF_1; TF_2: the expression value in TF_2; TF_3: the expression value in TF_3; TV_1: the expression value in TV_1; TV_2: the expression value in TV_2; the expression value in TV_3; TF_12: average of the expression value in TF_1 and TF_2; TF_13: average of the expression value in TF_1 and TF_3; TF_23: average of the expression value in TF_2 and TF_3; TV_12: average of the expression value in TV_1 and TV_2; TV_13: average of the expression value in TV_1 and TV_3; TV_23: average of the expression value in TV_2 and TV_3; TF: average of the expression value in TF_1, TF_2 and TF_3; TV: average of the expression value in TV_1, TV_2 and TV_3. Expression values for 6 libraries are presented as RPKM normalized by transformed counts. (b) Six main clusters are shown based on the heat map of DEGs.

Additional file 3.

KEGG pathways and the corresponding unigene numbers and IDs in the transcriptome of Lei bamboo.

Additional file 4.

Summary of the GO enrichment analysis of the DEGs in each comparison.

Additional file 5.

Primer sequences for qRT-PCR analysis.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Jiao, Y., Hu, Q., Zhu, Y. et al. Comparative transcriptomic analysis of the flower induction and development of the Lei bamboo (Phyllostachys violascens). BMC Bioinformatics 20 (Suppl 25), 687 (2019).

Download citation

  • Published:

  • DOI: