Skip to main content

Transcriptomic analysis identifies the key genes involved in stamen petaloid in lotus (Nelumbo nucifera)

Abstract

Background

Flower morphology, a phenomenon regulated by a complex network, is one of the vital ornamental features in Nelumbo nucifera. Stamen petaloid is very prevalent in lotus flowers. However, the mechanism underlying this phenomenon is still obscure.

Results

Here, the comparative transcriptomic analysis was performed among petal, stamen petaloid and stamen through RNA-seq. Using pairwise comparison analysis, a large number of genes involved in hormonal signal transduction pathways and transcription factors, especially the MADS-box genes, were identified as candidate genes for stamen petaloid in lotus.

Conclusions

Taken together, these results provide an insight into the molecular networks underlying lotus floral organ development and stamen petaloid.

Background

Lotus (Nelumbo), a basal eudicot in the family of Nelumbonaceae, is a perennial aquatic plant, comprising of two species Nelumbo nucifera Gaertn. and Nelumbo lutea Pers.. Its relative simple genome size (about 929 Mb) makes it ideal for cross breeding and genetic studies [1, 2]. The N. nucifera is mainly cultivated in Asia, especially in China as an important horticultural crop with very high economic value, and it has been cultivated for more than 2000 years [3, 4]. The sacred lotus can be divided into three groups according to their agronomic traits and primary utilization, namely, rhizome, seed, and flower lotus.

Flower lotus are widely cultivated for their aesthetic value. The color and shape of its flower are the two major ornamental features, which largely determine its ornamental values. The aspect of flower shape includes the number, size, and shape of the petals. There are four distinct organ types (sepals, petals, stamens, and carpels) that commonly constitute a flower and are arranged in different whorls. For the purpose of ornamentation, petaloid is one of the characteristics selected in ornamental plants breeding [5]. Based on previous studies on the ranuncukid genus Aquilegia, a model system for studying the evolution of petals [6], there are two major types of petaloid. The andropetaloidy, which originates from a modified stamen, and the other is derived from pre-existing sterile bracts variation [7]. In lotus, the flower morphologies are classified as a few-petalled group, double-petalled group, duplicate-petalled group and all-double-petalled group [3]. Several lotus cultivars have novel petaloid differentiation, such as the stamen petaloid and carpel petaloid. Lotus flowers exhibit numerous transition petal shapes with some having a stamen appendage, indicating that the petals may belong to andropetals. There are distinct stamen petaloid cultivars including ‘Fenhonglingxiao’, ‘Zhigaowushang’ and ‘Manjianghong’. However, the genetic pathways controlling this phenomenon of stamen petaloid remain unclear in lotus.

Due to its importance for both horticultural ornament and plant reproduction, floral patterning has been well studied, particularly in Arabidopsis and Antirrhinum. Based on previous studies, genetic pathways controlling petal organogenesis in Arabidopsis have been elucidated and the molecular mechanism of petal identity relies on specialized signal transduction [8, 9]. Majority of the identified genes involved in these processes encode hormone pathway or transcription factors [8, 10].

Numerous studies focusing on floral organ identity genes have been conducted, providing the well-known ‘ABCE model’ of flower development [11,12,13,14,15,16]. In Arabidopsis, A class genes include AP1 (APETALA 1) and AP2 (APETALA 2); B class genes includes AP3 (APETALA3) and PI (PISITTALA); C class contains only the AGAMOUS (AG); whereas, E class consist of four SEPALLATA (SEP1, 2, 3, 4) members. The above genes are MADS-box genes except for AP2 [11, 17]; with ‘quartet model’ suggesting that they interact with each other to form tetrameric complexes [15, 18]. Class A and E genes together with those in class B control petal development, with stamens being specifically controlled by class B genes in combination with those in class E and C [19]. Over the time, it has been widely believed for a long time that not all petals are homologous [20]. Loss of class C gene perhaps leads to the transition of the stamen to petal. Above-mentioned enlightenments indicate that the morphology of the stamen petaloid may be regulated by a complicated genetic pathway. However, it is still unknown whether this “ABCE” is also applicable in other plant species. Recently, several studies were performed using RNA-Seq to characterize the key genes controlling petal development in some non-model plant species, such as roses, Brassica rapa and chrysanthemum [21,22,23].

In lotus, the genome of an ancient lotus cultivar ‘China Antique’ has been sequenced and annotated [1], which facilitates studies of this plant on different aspects of molecular biology. Owing to the available genomic information, some transcriptomics studies have also been conducted on its rhizome growth and development, leaf development and flowering time [24,25,26,27]. However, there are still no studies conducted to explore its floral organ development. In the present study, to elucidate the mechanism of floral organ formation, especially the stamen petaloid, we conducted a comparative transcriptomic analysis among petal (P), stamen petaloid (Sp) and stamen (St) from an ornamental lotus cultivar ‘Fenhonglingxiao’ through RNA-seq technique. The results showed that several candidate genes might be involved in stamen petaloid in lotus.

Results

Stamen petaloid phenotype of N. nucifera

Flower shape is one of the important factors that determine the ornamental value of lotus. A normal flower contains four distinct organs, including sepals, petals, stamens, and carpels. The former three contribute greatly to the flower appearance and thus its ornamental value. As for lotus, the sequenced cultivar ‘China Antique’ contains normal petals, stamens, and carpels (Fig. 1a). For ornamentation, a series of flower lotus cultivars with different flower shape have been bred. Among them is the ‘Fenhonglingxiao’, a tropical lotus with floral aberration, in which stamens are converted to petal-like organs (Fig. 1b-f). Compared with ‘China Antique’, ‘Fenhonglingxiao’ has more petals and fewer stamens and the abnormal floral organ is known as stamen petaloid, which include stamen appendage being similar to petal (Fig. 1f). This feature provides a good model for the study of stamen petaloid.

Fig. 1
figure 1

Sacred lotus flower. a Flower of ‘China antique’. b and (c) Flower bud of ‘Fenhonglingxiao’. d Flower of ‘Fenhonglingxiao’ at the first day of blossom. e Flower of ‘Fenhonglingxiao’ at the fourth day of blossom. Irepresents Petal (P); IIrepresents Stamen (St); III represents carpel; IV represents Stamen petaloid (Sp). f Three floral organs, including St, Sp, and P. Bars are all 1 cm

Illumina HiSeq mRNA sequencing

To further investigate the molecular mechanisms underlying the stamen petaloid phenotype, comparative transcriptomic analysis among the petal (P), stamen petaloid (Sp), stamen (St) through RNA-Seq was conducted. The floral organs, including P, Sp, and St, were sampled when the flower bloomed on the first day (Fig. 1d), and were used to construct nine independent libraries: P1, P2, P3, Sp1, Sp2, Sp3, St1, St2, St3, with digital number representing the three biological replicates.

A total of 75.35 Gb clean data were generated from the nine libraries, with 73.11Gb high-quality (Q > 20) reads being selected, with each sample having over 6.51 Gb clean data used for further analysis. The number of transcripts identified in each sample, expressed in FPKMs (fragments per kilo base of transcript per million base pairs sequenced) are shown in Fig. 2a. Genes with normalized reads lower than 0.5 FPKM were removed from the analysis. In total, 27,053 transcripts were detected in the samples. Approximately 33% of the expressed genes were in the 0.5–10 FPKM range, while 36% were in > 10 FPKM (Fig. 2a). About 82.72–86.46% of the short clean reads were aligned against the ‘China Antique’ reference genome [1]. Among which, 73.63–76.29% of the reads were mapped to exon, 6.45–9.21% to intron, and 16.59–17.67% to intergenic (Table 1). For each of the three sampled floral organs, the Pearson correlation coefficient analysis showed they had good reproducibility among the three biological replicates (Additional file 1: Figure S1).

Fig. 2
figure 2

Information of the identified transcripts. a Abundance distribution of the expressed transcripts. b Venn diagram showing the overlap between the annotated genes from the lotus database (I) and the identified transcripts (II). c GO annotation of the novel genes

Table 1 Summary statistics of clean reads in the transcriptomes of lotus

Six databases (COG, GO, KEGG, KOG, SwissProt, Nr) were used to annotate all genes to provide comprehensive gene function information. Based on the annotation results (Table 2), 10,193 genes (36.66%) were annotated in COG; 21,434 genes (77.10%) in GO; 10,268 genes (36.93%) in KEGG; 14,775 genes (53.14%) in KOG; 20,768 genes (74.70%) in SwissProt and 27,784 genes (99.94%) in Nr. Altogether,27,802 genes were annotated. Moreover, 25,520 out of 27,802 genes were previously predicted in the lotus reference [1]. The rest 2282 genes could not be identified in the reference genes, representing novel genes (Fig. 2b). Of the 2282 novel transcripts, a total of 1836 genes were successfully annotated using all annotation databases (Additional file 1: Table S1). Also, 758 genes were annotated by GO assignments, being categorized into three major groups (cellular component, molecular function, and biological process), and 41 sub-categories (Fig. 2c).

Table 2 All genes annotation

Differentially expressed genes (DEGs) among floral organs

To identify the genes associated with the stamen petaloid, we conducted pairwise comparison among petal, stamen petaloid and stamen (Fig. 3). The FPKM was used to estimate the level of gene expression and the expression levels were compared in each floral organ to identify significant DEGs, with a fold change of no less than two and a False Discovery Rate (FDR) < 0.05. In total, 3900 genes were differentially expressed among these three organs (Additional file 2). We found that there were 1087, 2280 and 3114 DEGs between P vs Sp, St vs Sp, and St vs P, respectively (Fig. 3a). Of these DEGs, 364 genes were significantly differentially expressed among the three floral organs. The overlapped number of DEGs detected from P vs Sp and St vs Sp was 456, which might be specifically involved in the stamen petaloid formation. Additionally, comparing with Sp, 631 and 1824 genes were specifically regulated in petal and stamen, respectively. Moreover, comparing of P vs Sp, 982 and 105 transcripts were up- and down-regulated, respectively (Fig. 3b), while in St vs P and St vs Sp, their up-regulated genes were both less than the down-regulated genes. In the St vs P comparison, there were 776 and 2338 transcripts up- and down-regulated, respectively, and in St vs Sp, the up-regulated genes were 836 with the down-regulated genes being 1444.

Fig. 3
figure 3

DEGs among the three floral organs. a Venn diagram of the number of unique and common DEGs among the comparisons between P vs Sp, St vs P and St vs Sp. b The number of up-regulated and down-regulated DEGs in the comparisons between P vs Sp, St vs P and St vs Sp. The y-axis is the number of genes, and the x-axis is the different comparison

The clustered patterns of all DEGs expression during floral development were created using the Euclidean distance method associated with complete-linkage (Fig. 4a and b). K-means cluster analysis plotted the DEGs into six distinct expression patterns, K1-K6 (Fig. 4b). The K1 and K6 cluster contained 794 and 1185 genes respectively that had no obvious change in three different floral organs. The genes in cluster K2, K3 and K4 showed a similar expression pattern. K2 and K3 clusters having 394 and 484 genes, respectively, had an above the zero transcription level, and genes in cluster k2 were significantly changed among the three floral organs. Meanwhile, 748 genes in K4 cluster had almost all genes expression levels being below zero, except some in the stamen. The K5 cluster including 295 genes had a specific pattern being the only cluster with decreased expression in stamen.

Fig. 4
figure 4

Expressional patterns of the DEGs identified among the three floral organs. a The heatmap of DEGs in three floral organs. FPKM (fragments per kilo base of transcript per million base pairs sequenced) normalized log2 transformed counts were used to estimate the level of gene expression. The red (high expression) and green (low expression) represents the relative expression level. b The nine main clusters with different gene expression pattern types, namely K1-K6. The numbers of DEGs in each cluster are displayed

We used GO assignments to classify the functions of DEGs in pairwise comparisons of cDNA libraries between different floral organs (Additional file 1: Figure S2A). The GO enrichment revealed 20 categories of biological processes, 16 categories of cellular component and 17 categories of molecular function which had similar pattern in the DEGs and all genes among P vs Sp, St vs Sp and St vs P, except ‘structural molecular activity’ category in P vs Sp, whose percentage was not corresponding to the pattern (Additional file 1: Figure S2A). To further explore the functional classification and pathway assignment of DEGs, we performed KEGG analysis. Dozens of genes were enriched in phenylpropanoid biosynthesis, phenylalanine metabolism, flavonoid biosynthesis, pentose and glucoronateinter conversion, nitrogen metabolism, amino sugar and nucleotide sugar metabolism, ascorbate and aldarate metabolism, starch and sucrose metabolism, alpha-linolenic acid metabolism and photosynthesis (Additional file 1: Figure S2B).

DEGs involved in stamen petaloid

To elucidate the genetic regulation of the stamen petaloid formation phenomenon, the differentially expressed genes were filtered for those believed to be involved in stamen petaloid. As stated above, 364 DEGs, could be the specific genes through which stamen petaloid may be regulated. There were numerous genes assigned to phytohormone pathways, including the auxin (AUX), abscisic acid (ABA), jasmonic acid (JA), cytokinine (CTK), gibberellin (GA), brassinosteroid (BR), ethylene (ETH), and salicylic acid (SA) pathways. The transcriptome annotated as being involved in known phytohormone biosynthesis and signaling pathways are displayed in Additional file 3 with DEGs in the different floral organ presented in Additional file 1: Figure S3. Eleven out of 364 DEGs were annotated as being related to phytohormone (Fig. 5a and Additional file 1: Figure S4). Four of them were involved in AUX, namely, YUCCA (NNU_01284), AUX (NNU_25947), and GH3 (NNU_16281 and NNU_05194). The protein phosphatase 2C (PP2C, NNU_01507) and one short-chain dehydrogenase/reductase (SDR, NNU_07187) were shown to be involved in the ABA pathway. Three genes encoding the jasmonate ZIM domain-containing protein (JAZ, NNU_01965 and NNU_11760), the lipoxygenase (LOX, NNU_05041) involved in the JA pathway, GA biosynthesis gene encoding for gibberellin 2-beta-dioxygenases (GA2ox; NNU_16063), and protein EIN4-like (ETR, NNU_20013) involved in ETH pathway, considerably varied in expression in the three floral organs and peaked in stamen. The expression patterns of ten out of eleven genes in three organs were similar, while NNU_05194 had an opposite pattern (Fig. 5a and Additional file 1: Figure S4). However, all of them exhibited a medium expression level in stamen petaloid. The expression of BR, SA, and CTK pathway genes (not listed above) also showed varied expression among the three samples.

Fig. 5
figure 5

Heat maps of phytohormone-related genes and transcription factors (TFs) in comparison of three floral organs. a Eleven phytohormone-related genes. b 22 transcription factors

A large number of DEGs were transcription factors (TFs) encoding genes. In our data, there were 89 TFs in P vs Sp, 172 TFs in St vs P, and 102 TFs in St vs Sp with the total number of TFs identified being 225 (Additional file 1: Figure S5). In the pair comparison of the three floral organs, there were 22 TFs amongst the 364 DEGs. 21 TFs, including eight MYBs like (NNU_02286, NNU_02357, NNU_05105, NNU_09095, NNU_14378, NNU_15849, NNU_17837 and NNU_24843), three putative homologs of LIM family members (NNU_08737, NNU_19193, and NNU_21034), two transcripts encoding WRKY (NNU_01367 and NNU_22208), two homologs of BTB/POZ and TAZ domain-containing proteins (BT: NNU_13743 and NNU_22020), one transcript encoding bHLH (NNU_25984), AP2/ERF and B3 domain–containing TEM1 (NNU_00569), Homeobox-leucine zipper protein GLABRA 2 (NNU_21,434), trihelix transcription factor GT-3b-like (NNU_12600), one homolog of NAC (NNU_09288), and AGAMOUS-like MADS-box protein AGL15 (NNU_22830) showed a down-regulation change in St compared with P or Sp. Nevertheless, only one TFs (GATA22, NNU_11982) was up-regulated in St vs P (Fig. 5b and Additional file 1: Figure S6). Similar to the hormone pathway genes, they also had a medium expression level in stamen petaloid. In addition, the homolog of NAC (NNU_09288) and AGL15 (NNU_22830) were down-regulated with more than a five-fold threshold. The MADS-box family members are key regulators of floral organ identity [28]. We regard the functions of MADS-box with extremely high interest.

In total, there were 12 members of floral homeotic protein, including 11 MADS-box TFs and one APETALA2 (NNU_17043) observed in our data (Additional file 1: Figure S7 and Additional file 4). The MADS-box family members have been reported to be involved in the floral development [15]. Seven floral homeotic protein encoding genes exhibited the substantial changes in each pairwise comparison with a medium expression level in stamen petaloid (Additional file 1: Figure S7 and Additional file 4). Similar to NNU_22830, the AGL15 homolog gene (NNU_01653), two homologs of floral homeotic protein AGAMOUS (NNU_10192 and NNU_26656), and the homolog of AGL80 (NNU_09213) were down-regulated in St vs P. Nonetheless, the expression levels of two AGL6 like (NNU_08847 and NNU_26581) were up-regulated in the same pairwise comparison.

Confirmation of gene expression through qRT-PCR

To assess the quality of the RNA-seq data, a total of 29 genes, including 12 homeotic floral genes (11 MADS-box genes and one AP2), were selected as representatives of specific DEGs for qRT-PCR (Additional file 1: Figure S8). The specific primers are listed in Additional file 1: Table S2. We showed that 29 genes had significant correlations between qRT-PCR result and RNA-Seq data and thus providing guaranteed data reliability (Additional file 1: Figure S9).

To validate the involvement of these gene in petaloid formation, three DEGs encoding AUX pathway, 12 transcription factors included seven MADS-box genes, were chosen to examine their gene expression profiles in eight tissues (curly leaf, unfolded leaf, rhizome, flower bud, petal, stamen petaloid, stamen, and carpel) (Fig. 6). Systematic analysis of the expression patterns of 15 candidate genes may reveal their morphologic roles in different tissues. Six genes, including BT4 (NNU_13743), YUUCCA10 (NNU_01284), MYB (NNU_02357), NAC (NNU_09288), AGL80 (NNU_09213), and AGL15 (NNU_22830), showed similar uniform expression patterns. Their expression levels showed higher enrichment in floral organs and peaked in stamen but were absent in vegetative tissues. GH3 (NNU_16281) also showed similar expression profiles as these six genes, except for the transcript from rhizome. This indicates that NNU_16281 could also be essential for rhizome formation. However, GH3 (NNU_05194) had a higher expression level in flower bud and petal and weakly expressed in other organs. MYB (NNU_24843) was mostly detected in stamen and lowly expressed in petal. In our result, GATA (NNU_11982) was the exclusive gene, which was expressed at high level in vegetative tissues (curly leaf and unfolded leaf but rhizome). Meanwhile, it also had enrichment in flower bud, petal, and stamen petaloid, and a lower level in stamen and carpel. Similarly, two AGL6 like (NNU_08847 and NNU_26581) had high abundant in petal, with their transcript levels being higher in flower bud. Two AG (NNU_10192 and NNU_26656) had relatively higher expression levels in flower organs except for petal, while low expressions were detected in other organs. In this study, the spatial expression of all seven MADS-box genes significantly higher in floral organs compared to the low expression in vegetative organs (Fig. 6).

Fig. 6
figure 6

The expression pattern of fifteen candidate genes in lotus. CL: Curly Leaf; UL: Unfolded Leaf; Rhizome: R; FB: Flower Bud; Petal: P; Stamen petaloid: Sp; Stamen: St. Relative gene expressions were normalized by comparison with the expression of lotus β-actin (NNU_24,864), and using the 2-ΔΔCT method. Each gene used three biological replicates, with three technical replicates per experiment in all qRT-PCRs. The Error bars represented the SD for three biological replicates

Discussion

Morphological and transcriptomic analyses showing the petal features of petaloid stamen

Floral morphology, formed by specific organs, specifically the petal, contributes to the ornamental values in flowering plants. Stamen petaloid has been frequently selected as one of the paramount ornamental features in flower lotus breeding. Here, we selected a lotus cultivar ‘Fenhonglingxiao’ showing stamen petaloid (Fig. 1). Based on the petal, stamen petaloid, and stamen morphology, it could be easily concluded that petaloid stamen lies in the middle of transition from stamen to petal. To investigate the mechanism underlying stamen petaloid in this cultivar, a comprehensive transcriptome analysis was conducted among the petal, stamen petaloid, and stamen. The number of DEGs detected in P vs Sp was much less than those in St vs P and St vs Sp (Fig. 3), suggesting that the petaloid stamen is more similar to petal than stamen, which was further supported by the Pearson correlation coefficient analysis (Additional file 1: Figure S1). The total number of identified genes was 27,802, a high value than the annotated genes in the lotus database. Although 1065 annotated genes were not identified, there were 2282 novel genes. Together with our previous study [26], it’s clear that more transcriptomic studies are necessary to optimize the lotus genome annotation.

Involvement of phytohormones and transcription factors in stamen petaloid in lotus

Petals and stamens, being among the most important floral organs are located in the second and third whorls, respectively. The stamen petaloid grows between petal and stamen, and its morphology contains both of their features (Fig. 1a). This phenomenon may be regulated at the gene level as the petal spill out into the whorl of stamen or else the failure in expression of stamen identified genes. From the commonly expressed genes in St vs P, St vs Sp, and P vs Sp, genes controlling specific regulation in Sp were obtained.

Phytohormones are very crucial factors that affect organ genesis and development. Numerous genes involved in different phytohormone metabolism or signaling were differentially expressed among the three tissues. In this study, we specifically focused on the auxin pathway that controls floral organ formation. Auxin has been reported to play a crucial role in determining the identity of floral organs [29, 30]. In this work, auxin biosynthesis related gene YUCCA10 (NNU_01284) had a high expression in floral tissues (Fig. 5 and Additional file 3), especially in stamen, indicating its critical role in the formation and development of floral organs in lotus. Previous studies reported that anthers are the major sites that produce high concentrations of auxin in Arabidopsis [31]. Auxin might be synthesized more in stamen than in petal, but then possibly transported into petal for its action. Several auxin signaling related genes, showed higher expression levels in petal and petaloid stamen than stamen (Additional file 1: Figure S3 and Additional file 3), which is consistent with previous study [23], indicating that auxin may promote the formation of petaloid stamen. However, one auxin signaling related gene GH3 (NNU_16281) exhibited a predominant expression in stamen of lotus (Fig. 5 and Additional file 3). This result indicated that auxin signaling related genes have complicated function in regulation of floral organs formation.

The expression of numerous TFs changed drastically in lotus floral organs (Fig. 5 and Additional file 5). 225 TFs were identified as part of the DEGs and for further analysis, we focused on 22 TFs from the 364 DEGs obtained from the three organ comparison. Previously, it was reported that MYB is essential for floral organ development and plays a key role in pollen development in various plant species, such as rice, cotton and Arabidopsis [32,33,34]. In lotus, we have reported that there are potential 116 MYB genes in the lotus genome, with a number of them being related to flavonoid biosynthesis [35]. However, the mechanisms underlying floral morphologies are still obscure. Here, two genes encoding MYB (NNU_02357 and NNU_24843), showed a higher expression level in stamen and necessitates further investigation. It has been demonstrated that NAC transcription factor participates in floral-boundary morphogenesis in tomato [36]. NNU-09288, which was annotated as NAC-like TFs, probably participated in stamen development. However, NAC was also reported to control ethylene-regulated cell expansion in petals [37]. It has been shown that GATA transcription factors are negatively regulated by the floral homeotic genes APETALA3 and PISILLATA in petals and stamens [38]. Similarly, a GATA3 family member GNC participated in stamen and carpel development [39]. HAN encoding a GATA transcription factor transcribes at the boundaries between different floral whorls, whose mutant diminished boundary formation, suggesting that it is involved in controlling cell differentiation [40]. The key role of HAN in determining petal number also was demonstrated [41]. In our study, a GATA-like gene (NNU_11982) had a higher expression in the petal compared to stamens and stamen petaloid, indicating its potential roles in breaking down the boundary of stamen and petal and hence promoting stamen conversion to petal-like. The DEGs probably involved in the conversion of stamen to stamen petaloid is summarized in Fig. 7.

Fig. 7
figure 7

Summary of transcription-level regulation of stamen petaloid. Differentially expressed genes (DEGs) that may regulate stamen petaloid (those with more than a one-fold change in expression, FDR < 0.05). DEGs were either up-regulated (orange) or down-regulated (blue) as petaloid formation

Interestingly, we found fewer up-regulated DEGs than down-regulated DEGs in St vs Sp. After our critical selection, only two up-regulated DEGs were considered as the most potential candidates that cause the stamen petaloid formation, including auxin signaling related gene GH3 (NNU_05194) and GATA (NNU_11982). A recent study reports that GH3.5 may be a direct target of HAN (a member of GATA family) and the functions of HAN as a key repressor that regulates floral development, along with genes involved in hormone action [39]. Some GATA family transcript factors integrate auxin, GA, and CTK signaling pathway to control plant development [41, 42]. GATA transcription factors play roles in multiple developmental courses. Thus, the possible interaction between the GH3 (NNU_05194) and GATA (NNU_11982) could promote the change of stamen to stamen petaloid.

Regulation of stamen petaloid by floral homeotic genes

Members of the MADS-box family are known to be involved in floral organ specification [43]. Transcription factors including 40 MADS-box genes in sacred lotus were identified [2]. In this study, we detected 12 floral homeotic genes, including 11 MADS-Box genes and one AP2 (Additional file 4). Among them, seven MADS-box genes further analyzed through RT-qPCR on their expression pattern. We found that all of them had abundance expression in floral organs (Fig. 6). Previous studies have demonstrated that AG plays several key roles in floral organ determination, such as stamen identity and carpel specify [44,45,46]. In our study, two homologs of AG, NNU_10192 and NNU_26656 were both highly expressed in stamen petaloid and stamen (Fig. 6 and Additional file 1: Figure S8). The MADS-domain factors AGL15 like, NNU_22830 and NNU_01653 were most abundant in stamen; especially, the accumulation of NNU_22830 was down-regulated over five-fold in St vs P (Fig. 6 and Additional file 1: Figure S8). NNU_22830 was the only MADS-box gene common in St vs P, St vs Sp, and P vs Sp (Fig. 5). That indicated that NNU_22830 may play a vital function in stamen petaloid. Previous studies demonstrated that AGL15 and AGL18, along with SVP and AGL24, are inevitable and restrict floral gene expression [47, 48]. In rice, a AGL6 like gene (OsMADS6) determines floral meristem and floral organ identity [49]. Auxin-Responsive OsMGH3/OsGH3–8 was modulated by OsMADS6 and affected stamen differentiation [50]. Silencing SlAGL6 resulted in pale green petals in plants suggesting that SlAGL6 is involved in petal development [51]. In the present study, we found that two AGL6 genes (NNU_08847and NNU_26581) that were mainly expressed in petal, with their expression confirmed by qRT-PCR (Fig. 6 and Additional file 1: Figure S8). Thus, these two genes were speculated to be part of the genes that regulate petal development and promote stamen conversion to petal-like organs in lotus. One type I MADS-box gene, AGL80 (NNU_09213) expression was detected in the stamen within high-level. These findings are consistent with a prior report on this gene in Arabidopsis [52]. Meanwhile, it sharply declined in stamen petaloid indicating that losing NNU_09213 expression may lead to lose of stamen identity. These observations suggest that NNU_09213 is a crucial requirement for stamen determination.

As described above, we found that 15 candidate genes had medium expression in stamen petaloid compared with petal and stamen. And they are all involved in petal or stamen formation, which displayed obscure expression domains crossing the borders of petal and stamen causing the stamen petaloid phenomena.

Conclusion

This study covers the analysis of the comprehensive transcriptomic data to provide an overview of the transcriptional changes in different floral organs in lotus cultivar ‘Fenhonglingxiao’. In addition, through pairwise comparison analyses between three floral organs, we identified 3900 DEGs likely involved in floral organ identity. The DEGs were assigned to hormone pathways and transcription factors, which were reported previously as related to flower formation. Outstandingly, the main 11 hormone related genes and 22 TFs are involved in the conversion of stamen to stamen-petaloid. We report that 29 DEGs were detected in floral organs and confirmed by qRT-PCR, and had a significant accordance with RNA-seq result. Additionally, the expression pattern of 15 candidate genes was evaluated in eight different tissues by qRT-PCR. Two of candidate genes GH3.5 and GATA22, were specifically up-regulated in St vs Sp. We hypothesise that they function together in converting stamen to stamen petaloid. We identified some MADS-box genes being potential regulators of ambiguous boundaries of petal and stamen. Overall, further analysis on the mechanisms of floral formation in lotus will provide an insight and enable a deeper understanding of genetic involvement in the regulation of floral organ differentiation.

Methods

Plant materials

Nelumbo necifera cultivar ‘Fenhonglingxiao’ was used for transcriptome analysis. The cultivar has been conserved via rhizome in Wuhan Botanical Garden, Chinese Academy of Sciences (WBGCAS), Hubei Province, China for many years. The rhizomes of ‘Fenhonglingxiao’ were planted in three separate pools (3 m × 3 m) with three rhizomes per pool in the trial plot of WBGCAS on April 4th, 2016. When the flowers of ‘Fenhonglingxiao’ bloomed, petal, petaloid stamen and stamen in the same flower were sampled, respectively (Fig. 1). Meanwhile, curly leaf, unfolded leaf, flower bud, rhizome and carpel were also collected. All samples were collected at 10 A.M, snap frozen in liquid nitrogen, and stored at − 80 °C until used for RNA and DNA extraction. The petal, stamen petaloid, and stamen samples were used to construct nine libraries for RNA-seq, which were named as P1, P2, P3, Sp1, Sp2, Sp3, St1, St2 and St3, respectively.

RNA extraction, library construction, and sequencing

Total RNAs were extracted using a RNA reagent (OminiPlant RNA Kit, CWBIO, China), and treated with RNase-free DNaseI (Thermo, Shanghai, China) to remove genomic DNA contamination. RNA integrity was evaluated, with RNA integrity number (RIN) being above 6.5 for all samples. More than 1.5 μg of total RNA per sample was used as input material for the RNA sample preparations. For each floral organ of the cultivar, these nine RNA samples were subsequently used to create cDNA libraries and Illumina sequencing by Beijing Novogene Bioinformatics Technology Co., Ltd. The cDNA libraries were sequenced using the Illumina HiSeq 2500 high throughput sequencing platform. Transcriptome assembly was accomplished using Trinity (version r20140413p1) [53] with min_kmer_cov set to 2 by default, and all other parameters set default. All clean reads data were deposited in NCBI Sequence Read Archive (SRA, https://trace.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?) with accession number PRJNA417869.

Data processing and differential expression

Sequenced reads containing adapter, poly N or low-quality sequences (Q < 20) were removed, and the remainder were termed as clean reads. All the downstream analyses were based on clean data with high quality. The clean data was uploaded in the BMKCloud (https://www.biocloud.net/) for analysis. The data was mapped to the reference genome of Nelumbo nucifera (China lotus 1.1) [1] using TopHat2 Software v2.0.9 [54]. The default settings were considered as tolerance parameter, enabling mismatches of less than two bases. The transcripts were annotated using Cufflinks v2.1.1 [55]. For annotations of novel genes, they were aligned using BLAST (E value < 10− 5) [56] to National Center for Biotechnology Information (NCBI) non-redundant protein (Nr) database, NCBI non-redundant nucleotide sequence (Nt) database, and Swiss-Prot database and the sequences with the highest similitudes were retrieved. The differentially expressed genes (DEGs) between the two samples were carried out by using DESeq package (http://bioconductor.org/packages/release/bioc/html/DESeq.html). The transcript abundance (read counts) was estimated by RSEM v1.2.15 [57]. The expression level of genes was measured by FPKM method [55]. The DEGs were identified by false discovery rate (FDR) < 0.05 and a fold change ≥2. Any genes with an adjusted P-value < 0.05 were assigned as differentially expressed. Expression pattern analysis using Multiple Experiment Viewer software (MeV 4.9.0, https://sourceforge.net/projects/mev-tm4/files/mev-tm4/MeV%204.9.0/) was applied. Euclidean distance was used in the Distance method and K-means for hierarchical clustering was performed.

Gene ontology and KEGG Ortholog enrichment analysis

Based on Wallenius non-central hyper-geometric distribution [58], GO enrichment analysis of the differentially expressed genes (DEGs) was implemented by the GOseq R packages 1.10.0, which can be adjusted for gene length bias in DEGs. KOBAS (KEGG Orthology-Based Annotation System) software was used to test the statistical enrichment of DEGs in KEGG pathways (FDR ≤ 0.05).

Quantitative real-time PCR validation of RNA-seq data

A total of 29 DEGs (log2 FC > 1 or < − 1 and FDR < 0.05) involved in floral development were chosen for validation using quantitative real-time PCR (qRT-PCR) in petal, stamen petaloid and stamen. Among the 29 DEGs, 15 genes were chosen for analysis expression pattern in eight different tissues. Primers for qRT-PCR were listed in Additional file 1: Table S2. The qRT-PCR reactions were performed on the Bio-Rad using the SYBR (BioRad, http://www.bio-rad.com/). The reaction was initiated at 95 °C for 10 s, followed by 40 cycles of 95 °C for 15 s, 60 °C for 15 s and 72 °C for 30 s. Fluorescent products were detected in the last step of each cycle. Melting curve analysis was performed at the end of 40 cycles to ensure proper amplification of target fragments. The data presented for qRT-PCR experiments are the average relative quantities from three biological replicates, where every biological replicate is the mean of three technical repeats. Relative gene expressions were normalized by comparison with the expression of lotus β-actin (NNU-24864), and analyzed using the 2-ΔΔCT method. The data was indicated as mean ± SD.

Abbreviations

AG:

AGAMOUS

AP:

APETALA

DEG:

Differential expression gene

FDR:

False discovery rate

FPKM:

Fragments per kilo base of transcript per million base pairs sequenced

MADS:

MCM1, AG, DEFA and SRF

P:

Petal

PI:

PISITTALA

qRT-PCR:

Quantitative real-time PCR

Sp:

Stamen petaloid

St:

Stamen

TF:

Transcription factor

References

  1. Ming R, Vanburen R, Liu Y, Yang M, Han Y, Li LT, et al. Genome of the long-living sacred lotus (Nelumbo nucifera Gaertn). Genome Biol. 2013;14(5):R41.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  2. Wang Y, Fan G, Liu Y, Sun F, Shi C, Liu X, et al. The sacred lotus genome provides insights into the evolution of flowering plants. Plant J. 2013;76(4):557–67.

    Article  PubMed  CAS  Google Scholar 

  3. Wang Q, Zhang X. Colored illustration of lotus cultivars in China. Beijing: China Forestry Publishing House; 2005.

    Google Scholar 

  4. Guo H. Cultivation of lotus (Nelumbo nucifera Gaertn. Ssp. nucifera) and its utilization in China. Genet Resour Crop Evol. 2009;56(3):323–30.

    Article  Google Scholar 

  5. Tian X, Yu Q, Liu H, Liao J. Temporal-spatial transcriptome analyses provide insights into the development of petaloid androecium in Canna indica. Front Plant Sci. 2016;7:1194.

    PubMed  PubMed Central  Google Scholar 

  6. Kramer EM, Hodges SA. Aquilegia as a model system for the evolution and ecology of petals. Philos Trans R Soc B. 2010;365(1539):477–90.

    Article  Google Scholar 

  7. Takhtajan A. Evolutionary trends in flowering plants. New York: Columbia University Press; 1991.

    Google Scholar 

  8. Huang T, Irish VF. Gene networks controlling petal organogenesis. J Exp Bot. 2016;67(1):61.

    Article  PubMed  CAS  Google Scholar 

  9. Brioudes F, Joly C, Szécsi J, Varaud E, Leroux J, Bellvert F, et al. Jasmonate controls late development stages of petal growth in Arabidopsis thaliana. Plant J. 2009;60(6):1070–80.

    Article  PubMed  CAS  Google Scholar 

  10. Reeves PH, Ellis CM, Ploense SE, Wu MF, Yadav V, Tholl D, et al. A regulatory network for coordinated flower maturation. PLoS Genet. 2012;8(2):e1002506.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

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

    Article  PubMed  CAS  Google Scholar 

  12. Weigel D, Meyerowitz EM. The ABCs of floral homeotic genes. Cell. 1994;78(2):203.

    Article  PubMed  CAS  Google Scholar 

  13. Pelaz S, Ditta GS, Baumann E, Wisman E, Yanofsky MF. B and C floral organ identity functions require SEPALLATA MADS-box genes. Nature. 2000;405(6783):200.

    Article  PubMed  CAS  Google Scholar 

  14. Causier B, Schwarz-Sommer Z, Davies B. Floral organ identity: 20 years of ABCs. Semin Cell Dev Biol. 2010;21(1):73–9.

    Article  PubMed  CAS  Google Scholar 

  15. Theißen G, Melzer R, Rümpler F. MADS-domain transcription factors and the floral quartet model of flower development: linking plant development and evolution. Development. 2016;143(18):3259.

    Article  PubMed  CAS  Google Scholar 

  16. Irish V. The ABC model of floral development. Curr Biol. 2017;27(17):R887.

    Article  PubMed  CAS  Google Scholar 

  17. Bowman JL, Smyth DR, Meyerowitz EM. Genetic interactions among floral homeotic genes of Arabidopsis. Development. 1991;112:1–20.

    PubMed  CAS  Google Scholar 

  18. Thei G, Saedler H. Plant biology: floral quartets. Nature. 2001;409(6819):469–72.

    Article  Google Scholar 

  19. Wellmer F, Graciet E, Riechmann JL. Specification of floral organs in Arabidopsis. J Exp Bot. 2014;65(1):2047–8.

    Article  CAS  Google Scholar 

  20. Litt A, Kramer EM. The ABC model and the diversification of floral organ identity. Semin Cell Dev Biol. 2010;21(1):129–37.

    Article  PubMed  CAS  Google Scholar 

  21. Han Y, Wan H, Cheng T, Wang J, Yang W, Pan H, et al. Comparative RNA-seq analysis of transcriptome dynamics during petal development in Rosa chinensis. Sci Rep. 2017;7:43382.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Chang L, Zhiyong L, Xueling Y, Yiheng W, Qian M, Hui F. Transcriptome analysis of the Chinese cabbage (Brassica rapa ssp. pekinensis) petal using RNA-Seq. J Hortic. 2014;1(3):114.

    Google Scholar 

  23. Wang J, Wang H, Ding L, Song A, Shen F, Jiang J, et al. Transcriptomic and hormone analyses reveal mechanisms underlying petal elongation in Chrysanthemum morifolium ‘Jinba’. Plant Mol Biol. 2017;93(6):593–606.

    Article  PubMed  CAS  Google Scholar 

  24. Yang M, Zhu L, Xu L, Pan C, Liu Y. Comparative transcriptomic analysis of the regulation of flowering in temperate and tropical lotus (Nelumbo nucifera) by RNA-Seq. Ann Appl Biol. 2014;165(1):73–95.

    Article  CAS  Google Scholar 

  25. Yang M, Xu L, Liu Y, Yang P. RNA-Seq uncovers SNPs and alternative splicing events in Asian lotus (Nelumbo nucifera). PLoS One. 2015;10(4):e0125702.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Yang M, Zhu L, Pan C, Xu L, Liu Y, Ke W, et al. Transcriptomic analysis of the regulation of rhizome formation in temperate and tropical lotus (Nelumbo nucifera). Sci Rep. 2015;5(4):13059.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Yang M, Zhu L, Li L, Li J, Xu L, Feng J, et al. Digital gene expression analysis provides insight into the transcript profile of the genes involved in aporphine alkaloid biosynthesis in lotus (Nelumbo nucifera). Front Plant Sci. 2017;8:80.

    PubMed  PubMed Central  Google Scholar 

  28. Sablowski R. Control of patterning, growth, and differentiation by floral organ identity genes. J Exp Bot. 2015;66(4):1065.

    Article  PubMed  CAS  Google Scholar 

  29. Lampugnani ER, Kilinc A, Smyth DR. Auxin controls petal initiation in Arabidopsis. Development. 2013;140(1):185–94.

    Article  PubMed  CAS  Google Scholar 

  30. Robert SNLH, Crhak LK, Mroue S, Benkov E. The importance of localized auxin production for morphogenesis of reproductive organs and embryos in Arabidopsis. J Exp Bot. 2015;66(16):5029.

    Article  PubMed  CAS  Google Scholar 

  31. Aloni R, Aloni E, Langhans M, Ullrich CI. Role of auxin in regulating Arabidopsis flower development. Planta. 2006;223(2):315–28.

    Article  PubMed  CAS  Google Scholar 

  32. Kaneko M, Inukai Y, Ueguchi-Tanaka M, Itoh H, Izawa T, Kobayashi Y, et al. Loss-of-function mutations of the rice GAMYB gene impair alpha-amylase expression in aleurone and flower development. Plant Cell. 2004;16(1):33–44.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Leydon AR, Beale KM, Woroniecka K, Castner E, Chen J, Horgan C, et al. Three MYB transcription factors control pollen tube differentiation required for sperm release. Curr Biol. 2013;23(13):1209–14.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Xu Y, Iacuone S, Li SF, Parish RW. MYB80 homologues in Arabidopsis, cotton and Brassica: regulation and functional conservation in tapetal and pollen development. BMC Plant Biol. 2014;14(1):278.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Deng J, Li M, Huang L, Yang M, Yang P. Genome-wide analysis of the R2R3 MYB subfamily genes in lotus (Nelumbo nucifera). Plant Mol Biol Report. 2016;34(5):1–11.

    Article  CAS  Google Scholar 

  36. Hendelman A, Stav R, Zemach H, Arazi T. The tomato NAC transcription factor SlNAM2 is involved in flower-boundary morphogenesis. J Exp Bot. 2013;64(18):5497.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Pei H, Ma N, Tian J, Luo J, Chen J, Li J, et al. An NAC transcription factor controls ethylene-regulated cell expansion in flower petals. Plant Physiol. 2013;163(2):775–91.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Mara CD, Irish VF. Two GATA transcription factors are downstream effectors of floral homeotic gene action in Arabidopsis. Plant Physiol. 2008;147(2):707.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Zhang X, Zhou Y, Ding L, Wu Z, Liu R, Meyerowitz EM. Transcription repressor HANABA TARANU controls flower development by integrating the actions of multiple hormones, floral organ specification genes, and GATA3 family genes in Arabidopsis. Plant Cell. 2013;25(1):83–101.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Zhao Y, Medrano L, Ohashi K, Fletcher JC, Yu H, Sakai H, et al. HANABA TARANU is a GATA transcription factor that regulates shoot apical meristem and flower development in Arabidopsis. Plant Cell. 2004;16(10):2586–600.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. Ding L, Yan S, Jiang L, Zhao W, Ning K, Zhao J, et al. HANABA TARANU (HAN) bridges meristem and organ primordia boundaries through PINHEAD, JAGGED, BLADE-ON-PETIOLE2 and CYTOKININ OXIDASE 3 during flower development in Arabidopsis. PLoS Genet. 2015;11(9):e1005479.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. Richter R, Behringer C, Zourelidou M, Schwechheimer C. Convergence of auxin and gibberellin signaling on the regulation of the GATA transcription factors GNC and GNL in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2013;110(32):13192–7.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Smaczniak C, Immink RGH, Angenent GC, Kaufmann K. Developmental and evolutionary diversity of plant MADS-domain factors: insights from recent studies. Development. 2012;139(17):3081.

    Article  PubMed  CAS  Google Scholar 

  44. Ó'Maoiléidigh DS, Wuest SE, Rae L, Raganelli A, Ryan PT, Kwaśniewska K, et al. Control of reproductive floral organ identity specification in Arabidopsis by the C function regulator AGAMOUS. Plant Cell. 2013;25(7):2482.

    Article  PubMed Central  CAS  Google Scholar 

  45. Dreni L, Kater MM. MADS reloaded: evolution of the AGAMOUS subfamily genes. New Phytol. 2014;201(3):717–32.

    Article  PubMed  CAS  Google Scholar 

  46. Huang Z, Shi T, Zheng B, Yumul RE, Liu X, You C, et al. APETALA2 antagonizes the transcriptional activity of AGAMOUS in regulating floral stem cells in Arabidopsis thaliana. New Phytol. 2016;215(3):1197–209.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. Gregis V, Sessa A, Dorcafornell C, Kater MM. The Arabidopsis floral meristem identity genes AP1, AGL24 and SVP directly repress class B and C floral homeotic genes. Plant J. 2009;60(4):626–37.

    Article  PubMed  CAS  Google Scholar 

  48. Fernandez DE, Wang CT, Zheng Y, Adamczyk B, Singhal R, Hall PK, et al. The MADS-domain factors AGL15 and AGL18, along with SVP and AGL24, are necessary to block floral gene expression during the vegetative phase. Plant Physiol. 2014;165(4):1591–603.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  49. Li H, Liang W, Jia R, Yin C, Zong J, Kong H, et al. The AGL6-like gene OsMADS6 regulates floral organ and meristem identities in rice. Cell Res. 2010;20(3):299.

    Article  PubMed  CAS  Google Scholar 

  50. Yadav SR, Khanday I, Majhi BB, Veluthambi K, Vijayraghavan U. Auxin-responsive OsMGH3, a common downstream target of OsMADS1 and OsMADS6, controls Rice floret fertility. Plant Cell Physiol. 2011;52(12):2123–35.

    Article  PubMed  CAS  Google Scholar 

  51. Yu X, Chen G, Guo X, Lu Y, Zhang J, Hu J, et al. Silencing SlAGL6 , a tomato AGAMOUS-LIKE6 lineage gene, generates fused sepal and green petal. Plant Cell Rep. 2017;36(6):1–11.

    Article  CAS  Google Scholar 

  52. Portereiko MF, Lloyd A, Steffen JG, Punwani JA, Otsuga D, Drews GN. AGL80 is required for central cell and endosperm development in Arabidopsis. Plant Cell. 2006;18(8):1862–72.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  53. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  54. Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  55. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  56. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

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

    Article  CAS  Google Scholar 

  58. Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):R14.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

We would like to thank the reviewers for their helpful comments on the original manuscript.

Funding

This work was supported by National Natural Science Foundation of China (No. 31700197).

Availability of data and materials

All the data supporting our findings are contained within the manuscript. The clean sequence reads dataset supporting the conclusions of this article are available in the NCBI Sequence Read Archive with accession number PRJNA417869.

Author information

Authors and Affiliations

Authors

Contributions

Conceived and design the experiment: PY, ST and ZL; Performed the experiments: ZL, RND and JL; Analyzed the data: ST and ZL; Writing of the manuscript: ZL, RND and PY. All authors read and approved the final version of the manuscript.

Corresponding author

Correspondence to Pingfang Yang.

Ethics declarations

Ethics approval and consent to participate

In this study, the samples of N. nucifera ‘Fenhonglingxiao’ used in the experiment were planted in Wuhan botanical garden. This is an introduced species from Bangkok and a registered cultivar in China. No specific permission is required for use of this cultivar for experimental purpose.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1:

Additional figures and tables. Tables and figures in this file are referred to as Table S1-S2 and Figure S1-S9 in the main text. (PDF 1672 kb)

Additional file 2:

3900 DEGs. (XLSX 749 kb)

Additional file 3:

Hormone related genes. (XLSX 65 kb)

Additional file 4:

Floral organ homeotic genes. (XLSX 12 kb)

Additional file 5:

225 TFs. (XLSX 77 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) 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

Lin, Z., Damaris, R.N., Shi, T. et al. Transcriptomic analysis identifies the key genes involved in stamen petaloid in lotus (Nelumbo nucifera). BMC Genomics 19, 554 (2018). https://doi.org/10.1186/s12864-018-4950-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-018-4950-0

Keywords