Differential transcriptome analysis reveals insight into monosymmetric corolla development of the crucifer Iberis amara
© Busch et al.; licensee BioMed Central Ltd. 2014
Received: 3 June 2014
Accepted: 14 October 2014
Published: 19 November 2014
In the co-evolution between insects and plants, the establishment of floral monosymmetry was an important step in angiosperm development as it facilitated the interaction with insect pollinators and, by that, likely enhanced angiosperm diversification. In Antirrhinum majus, the TCP transcription factor CYCLOIDEA is the molecular key regulator driving the formation of floral monosymmetry. Although most Brassicaceae form a polysymmetric corolla, six genera develop monosymmetric flowers with two petal pairs of unequal size. In the monosymmetric crucifer Iberis amara, formation of the different petal pairs coincides with a stronger expression of the CYC-homolog IaTCP1 in the small, adaxial petals.
In this study, RNA-Seq was employed to reconstruct the petal transcriptome of the non-model species Iberis amara. About 9 Gb of sequence data was generated, processed and re-assembled into 18,139 likely Iberis unigenes, from which 15,983 showed high sequence homology to Arabidopsis proteins. The transcriptome gives detailed insight into the molecular mechanisms governing late petal development. In addition, it was used as a scaffold to detect genes differentially expressed between the small, adaxial and the large, abaxial petals in order to understand the molecular mechanisms driving unequal petal growth. Far more genes are expressed in adaxial compared to abaxial petals implying that IaTCP1 activates more genes than it represses. Amongst all genes upregulated in adaxial petals, a significantly enhanced proportion is associated with cell wall modification and cell-cell signalling processes. Furthermore, microarrays were used to detect and compare quantitative differences in TCP target genes in transgenic Arabidopsis plants ectopically expressing different TCP transcription factors.
The increased occurrences of genes implicated in cell wall modification and signalling implies that unequal petal growth is achieved through an earlier stop of the cell proliferation phase in the small, adaxial petals, followed by the onset of cell expansion. This process, which forms the monosymmetric corolla of Iberis amara, is likely driven by the enhanced activity of IaTCP1 in adaxial petals.
KeywordsBrassicaceae Monosymmetry CYC TCP1 RNA-Seq Microarray
In the co-evolution between plants and insects, the development of floral monosymmetry was an important step, facilitating angiosperm speciation as a response to the adaptation to the visual senses of insect pollinators . Monosymmetry, thus, likely functioned as a morphological key innovation and evolved independently in different angiosperm lineages ,.
Groundbreaking research almost two decades ago identified the TCP transcription factor CYCLOIDEA (CYC) as the molecular key regulator of monosymmetry development in Antirrhinum . CYC and its paralog DICHOTOMA are expressed in the adaxial part of developing Antirrhinum flowers, where they guide the acquisition of adaxial identities of second and third whorl organs ,. CYC belongs to the CYC2 clade of the TCP transcription factor family  and in all core eudicot species analysed so far, monosymmetry development is controlled by CYC2 clade genes (e.g. -).
The majority of crucifers (Brassicaceae) develop a polysymmetric corolla and only six genera form flowers with two petal pairs of different sizes . In Iberis amara, unequal petal pair formation correlates with a stronger expression of the CYC2 clade gene IaTCP1 in the smaller, adaxial petals. Comparison of adaxial and abaxial epidermal cell sizes revealed that petal size differences are due to a differential rate of cell proliferation . In a peloric flower variant, forming only large, abaxialized petals, the IaTCP1 expression is dramatically decreased. Transgenic Arabidopsis plants overexpressing the cruciferous CYC2 transcription factors IaTCP1 from Iberis amara and TCP1 from Arabidopsis, both, produce similar flowers with smaller petals. For plants overexpressing IaTCP1 this was shown to be due to a reduction in cell number . Contrarily, ectopic expression of CYC from Antirrhinum results in transgenic plants forming flowers with larger petals, a consequence of an increase in cell size . This demonstrates that the function of the two crucifer proteins is mainly conserved, whereas that of CYC from the more distantly related species Antirrhinum likely diverged . Iberis petals are initiated simultaneously as little bulges and the onset of an unequal size development can be detected around the start of stamen differentiation. From this point on, adaxial and abaxial petals develop differentially throughout flower development. The major difference in petal size, however, is acquired during late flower development, when a size difference of 1.6-fold, just after anthesis (stage A1) increases to 3.7-fold in fully mature flowers (stage A2) .
This raises the question about the molecular network that realises differential petal growth. Comprehensive research has been conducted analysing the genetic basis of general floral organ size determination, which is regulated through several independent pathways (reviewed in ,). Initial petal growth is achieved through cell proliferation that is later maintained only in restricted regions . Growth via cell division ceases and petals acquire their final size through cell elongation, a transition that seems to occur during later stages of flower development, after the maturation of microspores -. The switch to cell elongation goes along with an increased expression of cell wall synthesis and cell wall metabolization genes .
Thus, important determinators of final organ size are factors that control the timing of the cell proliferation arrest and the onset of cell expansion. The control of plant growth depends, in most cases, on an interplay between different hormones, which can affect either cell division (e.g. cytokinin), elongation (e.g. brassinosteroids, gibberellin) or, as auxin does, both processes . Growth regulators, like the transcription factors AINTEGUMENTA (ANT) and JAGGED, the cytochrome P540 monooxygenase KLUH or the auxin-inducible gene ARGOS (through ANT) have a positive effect on petal cell proliferation, mainly via affecting its duration -. The restriction of the cell division period during petal development, on the other hand, is accomplished e.g. through protein degradation by the E3 ubiquitin ligase BIG BROTHER and by the predicted ubiquitin receptor DA1 ,. MED25, a component of the mediator complex, has a dual function in restricting the periods of cell proliferation as well as cell expansion . Similarly, also AUXIN-RESPONSE FACTOR8 is a negative regulator of cell division early in petal development and represses late cell expansion together with the transcription factor BIG PETAL ,. In contrast, ARGOS LIKE promotes cell expansion .
Employing RNA-Seq, a de novo re-assembled petal transcriptome was generated for global analysis of Iberis amara petal development and to identify genes differentially expressed in the two petal pairs. This allowed detailed insight into the molecular network realizing unequal petal growth in the non-model organism Iberis amara. A significantly enhanced occurrence of genes involved in cell wall organization and modification processes in adaxial petals implies that the formation of smaller, adaxial petals is a result of an earlier onset of cell expansion, compared to the large, abaxial petals. In addition, microarray analyses of transgenic Arabidopsis plants ectopically expressing different CYC2 transcription factors were performed in order to detect differences in target gene regulation.
Analysis of spatial and temporal cell division patterns during Iberis amarapetal development
The temporal and spatial cell division pattern during petal development was analysed by conducting in situ hybridisation studies in younger Iberis flower buds before anthesis with the cell cycle marker gene Histone4 from Iberis amara (IaH4). Initially, IaH4 is expressed uniformly in petal primordia (Figure 1D). In older flowers, after the start of petal differentiation, IaH4 expression declines in the proximal region and accumulates in the distal petal part (Figure 1E). This indicates that cell divisions cease first in the proximal region where cells expand into elongated claw cells, while cell divisions in the blade continue. Moreover, there are more cells expressing IaH4 in the abaxial petal blade compared to adaxial petals (Figure 1E and F), as shown for an earlier stage by . Around anthesis, when petals completely cover reproductive structures, IaH4 expression concentrates predominantly at the distal petal tips and only few IaH4 expressing cells are detectable in the proximal claw region (Figure 1F).
These data show that Iberis petals mature along a proximo-distal gradient with a decline in cell division moving from the base to the petal tip. Furthermore, cell division still persists around anthesis in the blade tips, promoting the establishment of the final 3.7-fold size difference between mature adaxial and abaxial petals.
Reconstruction of the Iberis amarapetal transcriptome
Prior to the detection of genes, which are differentially expressed in the Iberis amara corolla by RNA-Seq, a reference petal transcriptome was reconstructed. Two replicates were harvested from adaxial and abaxial petals of stage A1 flowers as IaTCP1 shows at this stage the highest expression in adaxial petals and the largest expression difference between adaxial and abaxial petals . In order to avoid genetic variance due to a lack of isogenic lines from this non-model organism that barely self-pollinates, all petal material was harvested successively from one individual plant.
De novo assembly statistics
Read counts & nucleotides (Gb)
Total reads (4 libraries)
Statistics of contigs
Max sequence size (nt)
Min sequence size (nt)
Average sequence size (nt)
Contig length (nt)
Number of contigs
> = 10,000
Total contig number
About 59% of all contigs are smaller than 500 nucleotides, indicating that likely a larger portion of transcripts is not covered completely by their corresponding contigs. About 22% of all contigs have a size of up to 1,000 nt, 14% are between 1,000-2,000 nt long, followed by about 5% between 2,000-5,000 nt and 0.2% between 5,000-10,000 nt. Only four contigs are larger than 10,000 nucleotides (Table 1C). With the established petal transcriptome re-assembly a global analysis of petal development could next be accomplished.
Annotation of Iberisgenes involved in petal development
All 52,081 contigs were annotated based on BLASTX searches against an Arabidopsis protein database in order to allow a functional categorization of Iberis petal transcripts and to facilitate further comparison with Arabidopsis expression data. 45,001 Iberis contigs matched to a corresponding Arabidopsis Genome Initiative (AGI) code. 59.6% of matching contigs retrieved an AGI code that was found also as a match to other contigs. This is likely due to the fact that a large proportion of Iberis transcripts are not covered completely by their respective contigs and therefore several contigs were assigned to the same AGI code. Similarly, detection of an unexpected large number of genes encoding for transcription factors in Marchantia by RNA-Seq was explained as being caused by a "fragmentary" dataset, likely exceeding the number of actual transcription factors present . For this reason, here, contigs with the same AGI code were treated as partial transcripts of the same gene and subtracted, yielding 18,139 putative unigenes. Next, applying an e-value cutoff of e−4 resulted in 15,983 (88%) likely unigenes with a putatively conserved function in Arabidopsis. Contigs without a hit might be either assembly artefacts or are unique to Iberis.
Based on their respective AGI codes, the 15,983 Iberis unigenes were grouped into functional categories provided by MapMan using the classification superviewer tool ,. For an inter-species transcriptome comparison, microarray data from 13,492 genes expressed in Arabidopsis petals at stage 15 was taken. Stage 15 starts approximately one day after flower anthesis  and is, thus, closest in development to the Iberis stage A1, defined as the stage just after anthesis .
Characterization of Iberispetal development
RNA-Seq reads and mapping statistics
Total clean reads
Percent mapped uniquely
Functional categories of most abundant stage A1 Iberis petal transcripts
Number of contigs in
Identification of transcripts differentially expressed in adaxial and abaxial petals
Classification and analysis of genes differentially expressed between adaxial and abaxial petals
↑ Adaxial petals
↑ Abaxial petals
1 - 1.9
2 - 4.9
5 - 9.9
Fold change ≥ 2
Functional category (contig nr)
Not assigned (81)
Cell wall (21)
aD petal genes
Cell wall organization
Cell wall modification
Clathrin coat assembly
phs, light reaction
phs, light harvesting
Overall, the data indicate that a significantly higher gene number governs the development of small adaxial petals compared to the control of large abaxial petal formation.
Determination of gene categories in the adaxial petal transcriptome
To get an overview of developmental processes participating in the formation of the two petal types, differentially expressed genes with a minimal fold-change of two were classified into functional MapMan categories. No categories that consisted of more than three members, one of the applied threshold criteria, were found for the 21 transcripts upregulated in abaxial petals. Amongst the 313 transcripts upregulated in adaxial petals, about 25% (81) were not assigned into any category, encoding for proteins with unknown functions (Table 4B and Additional file 3). A large number of genes (26) is grouped together as encoding for various enzymes in the category misc and 25 transcripts are involved in protein metabolism, such as protein modification, targeting and degradation. The category transport comprises 24 genes, encoding for proteins transporting inorganic compounds, primary metabolites, ions, metal or auxin. Categories cell wall, signalling and photosynthesis were also detected as being significantly enriched when applying gene ontology (GO) terms of AGI codes linked to the adaxial-specific transcripts and are described below. Amongst the category RNA (23 members), 20 encode for various transcription factors, including IaTCP1. In addition, 16 adaxial upregulated genes perform a role in the adaptation to biotic and abiotic stress.
To gain a more detailed understanding of the developmental mechanisms and the respective genes involved in generating smaller adaxial petals, significantly enriched GO terms amongst the 313 adaxial-specific transcripts were revealed using FatiGO .
Individual genes are repeatedly allocated to related GO terms. GO terms associated to cell wall-related processes, such as cell wall organization and modification are significantly enriched in adaxial petals (Table 4C). The in total 15 genes exerting cell wall-related functions encode for pectate lyases, pectin methyl esterases/inhibitors implicated in pectin degradation  and other enzymes that function in cell wall remodelling and loosening. Even more genes were assigned to the MapMan category cell wall (21) (Table 4B) and include, in addition to cell wall degrading/modifying enzymes, COBRA like 10 and an arabinogalactan protein. Both proteins constitute a putative link between the cell wall and the cytoplasm .
The GO term cell-cell signalling contains five transcripts encoding for small signal peptides that play a role in plant cell growth (Table 4C). These peptides were first discovered in tobacco leaf extracts due to their ability to induce a rapid pH increase in cell suspension cultures and were, hence, termed RALFs (rapid alkalinization factors) ,.
The GO term flower morphogenesis includes three transcription factors involved in petal development, namely PETAL LOSS, BLADE ON PETIOLE 1 and 2. The categories endocytosis, membrane budding and clathrin coat assembly include genes that encode for ENTH/ANTH/VHS superfamily proteins, key regulators of the endosomal system involved in membrane trafficking .
Unexpectedly, 17 genes were found to be associated with photosynthesis, including the GO terms photosynthesis, light harvesting, light reaction and carbon fixation. In order to analyse whether photosynthesis-related processes contribute to flower monosymmetry formation in Iberis, this finding was scrutinized by detailed qPCR analyses. Therefore, a phosphoribulokinase homolog from Iberis was analysed as a representative for the category photosynthesis and confirmed the RNA-Seq data showing a stronger expression in adaxial petals (Additional file 4A). Next, qPCR was performed on separately harvested petal blade and claw material of the two petal types to test if the greenish claw cells, which are present in both petal types (Figure 1B), contribute to the detection of photosynthesis-related genes. The phosphoribulokinase transcript is not differentially expressed between adaxial and abaxial petal blades or claws (Additional file 4B), implying that photosynthesis genes do not contribute to the formation of smaller adaxial petals. In both petal types, the greenish claw region, where the phosphoribulokinase is stronger expressed (Additional file 4C), is similar in size and shape but blade sizes differ significantly. Therefore, the differential expression of photosynthesis-related genes detected with RNA-Seq might be due to the fact that harvesting whole petals diluted out claw-expressed genes stronger in large abaxial petals compared to small adaxial petals. As a consequence, photosynthesis-related genes were excluded from further analyses.
Expression analysis of genes upregulated in adaxial Iberispetals via qPCR
From these data, a general picture emerges, where adaxial petals are characterized by a pool of transcripts enriched in genes involved in cell wall organization and signalling. This implies that smaller petals in comparison to larger petals are the result of an earlier stop of cell division. Concomitantly, up-regulation of genes participating in cell elongation and differentiation occurs earlier in the small, adaxial petals compared to the larger, abaxial petals, the latter remaining longer in the cell proliferation phase. Further genes were selected from transcripts upregulated in adaxial petals (Additional file 3) based on their annotation and a putative relation to such processes.
A first interesting candidate for a participation in cell expansion processes is COBRA-like 10 (COBL10), a member of the COB multigene family. COBRA, the first analysed member, is strongly expressed in the root elongation zone and was proposed to be a key regulator of oriented cell expansion . Similarly, COBL10 is involved in the control of the directional and likely female-guided pollen tube growth . In both biological samples, the Iberis COBL10 is upregulated in adaxial petals compared to abaxial petals (Figure 4C), as shown with RNA-Seq. Another gene upregulated in adaxial Iberis petals is P-glycoprotein (PGP) 13 (Figure 4F), which belongs to the plant ABC transporter superfamily from the category transport (Table 4B). PGPs have been shown to fulfil a role in plant growth and development as they mediate the cellular and long-distance transport of auxin . Cytokinin oxidase 3 (CKX 3) catalyzes the breakdown of the plant hormone cytokinin and by that exerts a negative influence on cell division . Therefore, CKX 3 with its RNA-Seq expression fold change of 1.9 being slightly below the two-fold threshold was included into further qPCR analysis. The stronger expression of the Iberis CKX 3 homolog in adaxial petals (Figure 4G) might be involved in reducing the rate of cell division, prior to the switch into cell elongation. Also Rcd1-like, a protein involved in cell differentiation from the category RNA (Table 4B) shows a higher expression in adaxial petals (Figure 4H). The mammalian Rcd1 gene is a transcriptional cofactor with a role in the transition of cells to differentiate . The stronger expression in adaxial petals of a gene that shows highest homology to the Arabidopsis gene At5G42680 encoding for a protein of unknown function (Figure 4I), confirms the differential expression of a member of this highly abundant class of transcripts upregulated in adaxial petals.
In summary, qPCR analysis strengthens the observation that differential petal size development is accomplished via an earlier onset of cell elongation and differentiation in smaller adaxial petals, as indicated by an upregulation of genes that are involved in cell wall remodelling, repression of cell division, cell differentiation and signalling.
Interspecies comparison of CYC2 target genes via ectopic CYC2 expression in the heterologous system Arabidopsis thaliana
Transgenic Arabidopsis plants overexpressing the two crucifer CYC2 transcription factors IaTCP1 and TCP1 produce flowers with smaller petals, whereas ectopic expression of CYC from Antirrhinum results in the formation of larger petals ,. To compare the effect of the CYC, IaTCP1 and TCP1 proteins on downstream regulatory networks in Arabidopsis, a microarray analysis was conducted. Material for RNA extraction was harvested from young inflorescences, comprising flower bud stages before anthesis, of transgenic T2 plants overexpressing these genes.
The CYC2 transcription factors CYC and IaTCP1, which govern monosymmetry development, address both a larger number of target genes in the heterologous Arabidopsis system than TCP1. This tendency is more obvious regarding only the number of target genes, which are exclusively addressed by a single CYC2 protein, namely 348 by CYC, 249 by IaTCP1 and only 42 by TCP1. The capability of the CYC and IaTCP1 genes to address a larger number of downstream targets might reflect their function in the control of monosymmetry formation, even when expressed in a heterologous species forming a polysymmetric corolla such as Arabidopsis. 117 target genes are commonly addressed by all three transcription factors and might thus not be linked to specific functions in symmetry regulation. The highest number of overlappingly regulated genes is found between the two crucifer transcription factors (152), whereas the smallest number of targets is mutually addressed by CYC and TCP1 (19).
Enriched GO terms in the pool of IaTCP1 target genes in Arabidopsis inflorescences
Genes ↓ by IaTCP1
Arabidopsis infl. genes
Cell wall organization
Cell wall modification
Plant type cell wall org.
Unidimensional cell growth
Plant type cell wall loosening
sec. metabolic process
H2O2 catabolic process
Response to red light
RNA-Seq and microarray analyses, both, detected a differential regulation of genes involved in cell wall processes. However, in young transgenic Arabidopsis inflorescences all cell wall-related transcripts are downregulated following an ectopic IaTCP1 activity, whereas they are upregulated in older Iberis petals after anthesis, in correlation with and presumably also in response to a high endogenous IaTCP1-expression. This indicates that the spatial and temporal context of CYC2 activity is crucial for its proper effect on cell growth and differentiation.
Monosymmetry in Iberisis achieved through differential cell proliferation during late flower development
Iberis petals are initiated simultaneously and unequal growth begins around the onset of stamen differentiation, from where it proceeds and gains its maximal size difference after anthesis, likely driven by differential cell proliferation . This requires cell division to continue throughout later phases of flower development, which was demonstrated by monitoring the cell division pattern via Histone4 expression from Iberis in different petal stages. In older petals, cell division still takes place close to anthesis in both petal types, where it is restricted mainly to the petal tips. Similar to Iberis, cell division in Arabidopsis petals increases throughout flower development until stage 11 and can still be visualized in stage 12, just before anthesis .
After anthesis, adaxial Iberis petals increase 2.8-fold, whereas abaxial petals increase 6.7-fold to reach their final size. In contrast, cells of both petal types only double in size from anthesis to maturity . These data indicate that cell division – at least in abaxial petals – proceeds beyond anthesis.
Transcriptome analysis of petal development of Iberis amara
RNA-Seq is being used more and more often as a convenient tool for comparative transcriptome analyses of species with inavailable genome sequences -. Applying an RNA-Seq approach, we investigated the Iberis transcriptome at the late petal developmental stage A1, just after anthesis of flowers.
Comparison of the Arabidopsis and Iberis petal transcript distribution into functional categories indicates that molecular processes governing late petal development are similar in the two crucifers. Interestingly, for the largest group of the Iberis petal transcriptome, comprising 31% of the transcripts, no function could be assigned yet. A similar situation with 28% of these genes is observed for Arabidopsis. Interestingly, comparative analysis of different petal stages in Antirrhinum, where flowers form an even more elaborate and complex shaped corolla, revealed that the frequency of the unclassified protein group accounts for over 50% amongst all transcripts expressed in petals . This suggests that a large proportion of genes exert an intriguing, yet uncharacterized, function during petal morphogenesis.
Analysis of the most abundant transcripts in Iberis petal stage A1 revealed that various enzymes, subsumed in the category misc, are highly expressed in the two petal types, amongst them GDSL-lipases and lipid transfer proteins. Both enzyme types have a function in lipid metabolism and have been described as putative target genes of the petal and stamen-specifying B-class MADS box transcription factor during petal development in Antirrhinum and Arabidopsis ,. The cytochrome P450 gene CYP76C3, which is amongst the most highly expressed Iberis petal genes, is also strongly expressed in various Arabidopsis floral organs, where it participates in the linalool metabolism. Linalool, a well-studied attractant for pollinators, is converted further into its oxigenated derivates by CYP76C3 in autogamous Arabidopsis flowers , which likely function in plant defence against pathogens and insect herbivores . Two lipoxygenases, LOX1 and LOX2, which exert functions in hormone metabolism, are also highly expressed in Iberis petals. The LOX2 protein catalyzes the first dedicated step in jasmonic acid (JA) biosynthesis and is activated by the transcription factor TCP4 in leaves . TCP4 is also expressed in both Iberis petal types (Additional file 1), and since LOX2 is localized to chloroplasts , LOX2-dependent induction of JA likely takes place in the green petal claws, where JA might exert a negative effect on cell proliferation .
The monosymmetric corolla formation of Iberisis realized through differential timing of the switch from cell division to cell elongation
With almost 15-times more genes upregulated in the small adaxial petals compared to the large abaxial petals, it seems that more genes are required to restrict adaxial petal size during Iberis corolla development. The fact that the monosymmetry key regulator IaTCP1 is highly expressed in adaxial petals and represents one of the genes with the largest expression difference implies that IaTCP1 activates a larger number of target genes compared to its inhibitory transcription effects.
Differential Iberis petal development is likely driven by unequal cell division processes. Consistently, genes exerting negative effects on cell division are amongst the pool of transcripts upregulated in adaxial petals. In Arabidopsis, cytokinin oxidase (CKX) 3 catalyses the metabolic breakdown of the plant hormone cytokinin and plants with defective ckx genes form flowers with enlarged petals due to an increased cell number . The trihelix transcription factor PTL is expressed in regions of developing Arabidopsis flowers, where growth is repressed. PTL expression in the inter-sepal zone is assumed to keep cell proliferation under control to allow the outgrowth of petal primordia. The loss of PTL function results in an increased cell proliferation in that region . Petals are initiated in response to localized auxin accumulation, and Lampugnani and colleagues  speculate this could be evoked via a supportive role of PTL in non-polar auxin efflux through PGPs. PGPs belong to the superfamily of ABC-transporters and function in mediating auxin transport . PGP13, which is also transcriptionally enriched in adaxial Iberis petals, belongs to the same clade as the auxin exporters PGP1 and PGP19 , for which double mutant plants show severe growth defects .
In addition to transcripts negatively affecting cell proliferation, a large number of genes related to cell wall modifying processes are upregulated in adaxial petals. After an initial phase of organ growth by cell division, final organ size is realised by cell elongation ,. The primary cell wall of land plants is composed of a strong cellulose microfibril framework cross-linked through xyloglucan (XG) chains and embedded in a stabilizing matrix of pectin polysaccharides and structural proteins. Plant cell walls provide stability by withstanding tensile stress, and cell growth results from changes in wall extensibility by relenting to this stress ,. Amongst the group of cell wall modifying proteins in adaxial Iberis petals are several pectinases, enzymes that function in the de-/methylesterification or degradation of pectin  and, by that, might contribute to cell growth via precisely tuning the proportion of pectate in the cell wall . A Xyloglucan Endotransglucosylase (XET), which cleaves and re-attaches xyloglucan chains in vitro was found to be upregulated in adaxial petals. XET activity correlates with regions undergoing cell expansion, indicating that XETs have a role in cell expansion, e.g. through the biogenesis of newly formed cell walls (reviewed in ). In addition, several Iberis glycosyl hydrolases (GHs) were upregulated in adaxial petals, which form a diverse group of enzymes that hydrolyze glycosidic bonds . One Arabidopsis homolog of the GHs belongs to the subfamily of the endo-1,4-ß-glucanases KORRIGAN, which is essential for cell elongation and cellulose synthesis in Arabidopsis hypocotyls . Several transcripts encoding for a COBRA-like protein, two proline rich extension-like receptor kinases (PERKs), an arabinogalactan protein and a formin are upregulated in adaxial petals and might act as a potential connection between the cytoplasma and cell wall. These proteins are, thus, part of the molecular continuum that connects the cytoskeletal network, the cytoplasm and the cell wall, and which contributes to cell expansion control in Arabidopsis (reviewed in ).
Arabidopsis homologs of several of the floral cell wall-related genes, such as a polygalacturonase (At1G02790), a XET (At1G65310) and a pectin lyase like member (At5G07430) are also predominately expressed in the cell elongation zone of root tips (visualized with the Arabidopsis eFP browser at bar.utoronto.ca) and only weakly in the preceeding meristematic region, where growth is realized by cell divisions. This supports their general function in cell wall modification, both in a floral and root context, contributing to tissue growth by mediating cell elongation processes.
The determination of final organ size requires the precise coordination of both, cell division and cell elongation processes, and likely depends on a tight cellular communication. Amongst all adaxial-petal specific genes, a group of transcripts involved in cell-cell signalling is significantly enriched. These RALFs are small peptides, which are important for many developmental processes and regarded as novel plant hormones . In grass leaves, RALFs were found to be mainly expressed in young elongating cells beyond the meristematic zone, and expression is almost absent in mature cells . As RALFs are thought to exert a negative effect on cell growth, their expression strength and location is likely precisely adjusted during organ development to achieve an accurate guidance of cell expansion without completely blocking it ,.
Summing up, a significantly increased number of transcripts involved in cell wall modification and cell-cell signalling, together with genes regulating the restriction of cell division and auxin transport were detected in the pool of 313 transcripts enriched in small, adaxial Iberis petals. Hence, unequal petal growth is likely achieved via an earlier stop of cell division in adaxial petals with a subsequent earlier onset of cell expansion, whilst abaxial petal cells continue dividing. In Iberis, cell division ceases in an acropetal gradient, where cells at the tip keep on dividing longest, similar as in other plant species ,. Thus, an unequal stop of cell division might be realised through a faster acropetal progression of a cell division stop front in adaxial petals compared to abaxial petals.
Ectopic expression of IaTCPcauses expression changes of cell wall related genes in Arabidopsis
Microarray analyses aimed at comparing the molecular networks that are addressed in inflorescences by IaTCP1 and TCP1 causing smaller petal formation, and by CYC, giving rise to the development of larger petals.
A considerable number of genes (117 out of 600, 564 and 328, respectively) is addressed by all three CYC2 proteins together. These might be transcripts that constitute a common basal set of TCP/CYC2 transcription factor targets genes in Arabidopsis thaliana. Interestingly, in the heterologous Arabidopsis system, more genes are exclusively addressed by IaTCP1 from Iberis (249) and CYC from Antirrhinum (348), developing a monosymmetric flower, than by TCP1 from Arabidopsis (42), forming a polysymmetric corolla. These quantitative differences indicate that IaTCP1 and CYC might regulate more target genes, implying that a larger downstream regulatory network realizes monosymmetric compared to polysymmetric corolla development. Therefore, in addition to cis-regulatory changes that led to a heterochronic CYC2 gene expression shift in the Brassicaceae, promoting corolla monosymmetry evolution ,, also a divergence of trans-acting CYC2 proteins likely contributed to the establishment of this morphological novelty.
The largest overlap in target gene number is shared between IaTCP1 and TCP1. These 152 genes comprise several transcripts involved in cell wall-related and signalling processes (Additional files 6 and 7) and the activity of these exclusive crucifer targets could be essential for the petal size restriction in Arabidopsis, as only transgenic plants ectopically expressing a crucifer CYC2 protein develop flowers with smaller petals . Amongst the population of target genes downregulated in response to heterologous IaTCP1 expression in Arabidopsis, again a high proportion is involved in cell wall modification processes. This indicates that IaTCP1 addresses these genes in Arabidopsis and that the high abundance of this gene class detected with RNA-Seq in adaxial Iberis petals is the effect of the strong IaTCP1 expression there. While cell wall genes are positively affected in older adaxial Iberis petals, where IaTCP1 expression is strongly enhanced, IaTCP1 activity has a negative effect on cell wall genes in young Arabidopsis inflorescences overexpressing IaTCP1. TCP transcription factors are known to exert either an activating or repressive effect on organ growth, depending on their tissue specific expression ,,. As TCP proteins act as dimers , the outcome of their activity might differ in dependence of available interaction partners in different organs or developmental stages. Similar development-dependent regulatory effects were shown previously for the KNOX/BEL transcription factor BELLRINGER  (BLR), which interacts with other BEL and KNOX proteins ,. One of its targets, the pectin methylesterase PME5 is repressed by BLR in the shoot meristem, to allow coordinated formation of primordia, but is activated by BLR in internodes, where it promotes cell elongation .
An RNA-Seq expression analysis on adaxial and abaxial stage A1 Iberis petals reveals a significantly increased number of cell wall-related transcripts in adaxial petals. This indicates that differential petal growth, a characteristic of the monosymmetric Iberis corolla, is achieved through an earlier stop of the cell proliferation phase followed by the onset of cell expansion and differentiation in the smaller adaxial petals, while abaxial petal cells likely still continue growth via cell division. As ectopic expression of IaTCP1 in Arabidopsis addresses cell wall-related genes, the enhanced occurrence of this gene class in adaxial Iberis petals is very likely also a direct consequence of the strong IaTCP1 expression there. Moreover, crucifer monosymmetry evolution likely relied on the combination of cis-regulatory changes, leading to a late adaxial IaTCP1 expression and also on trans-regulatory changes mediated by differences in the CYC2 proteins themselves, altering the set of potential target genes that contribute to the large petal size differences in Iberis amara.
Plant material and growth conditions
Iberis amara seeds were purchased from Saatgut Kiepenkerl (Everswinkel, Germany) and have been cultivated in the greenhouses for several generations. Plants were grown under standard greenhouse conditions at the University of Osnabrück. Transgenic Arabidopsis thaliana (Col.) plants were grown at 20°C with 13-hrs light/11-hrs dark cycles in growth chambers.
Illumina sequencing, de novoassembly and RNA-Seq analysis
For RNA-Seq, two technical replicates of stage A1 adaxial (samples aD_1 & aD_2, Table 2) and abaxial (samples aB_1 & aB_2, Table 2) petals were harvested from one single plant, in order to avoid genetic variation, that would hinder the subsequent transcriptome de novo assembly. Plant material was flash frozen in liquid nitrogen and total RNA was isolated with the RNeasy Plant Mini kit (Qiagen, Valencia, CA), including an on-column DNA digestion step.
The library preparation, Illumina sequencing, de novo assembly, read mapping and significance statistics, which are described below, were performed by BaseClear BV, the Netherlands (www.baseclear.com). From the total RNA, mRNA was purified, fragmented and converted to double stranded cDNA using the Illumina TruSeq RNA-Seq library preparation kit. Paired-end reads were obtained after 50 cycles on the TruSeq Illumina sequencer. Sequence reads were generated using the Illumina Casava pipeline version 1.8.2 and after passing quality assessments, filtered and trimmed reads, excluding duplicates, were assembled using the CLC Genomics Workbench version 5.1.0. Clean reads, including duplicates, from each sample were mapped individually to the consensus transcriptome and expression values for contigs from each sample were calculated using the CLC Genomics Workbench. A read consists of intact paired-end read pairs and single reads which, when mapping uniquely, are both counted as one hit. Broken pairs, where one partner is missing, were omitted from the quantification. Expression of each contig is displayed through an RPKM value for each petal type, which is the mean expression from the two technical replicates. Significant expression of each contig was evaluated performing a Baggerley's test , followed by a Bonferroni correction for multiplicative errors.
Data analysis and global gene expression comparisons
All contig sequences were blasted against an Arabidopsis protein database (TAIR10_pep_20110103_representative_gene_model) using BLASTX. Multiple contigs that retrieved the same Arabidopsis AGI-code were treated as partial transcripts of the same gene. Consequently, respective expression values from contigs with the same AGI code were totalled. Differentially expressed genes were defined as such when they had a p-value for statistical significance below 0.01 and had an expression fold change equal to or above two.
Functional annotation of Iberis sequences is based on the assignment of each contig to its respective Arabidopsis ortholog via a BLAST-search. For all downstream analyses only contig/AGI code matches that are supported with an e-value equal to or below e-4 were considered. Functional classification was performed using the classification superviewer tool at bar.utoronto.ca , based on MapMan annotations (file Ath_AGI_LOCUS_TAIR10_Aug2012.txt). Only categories with a p-value below 0.05 and with more than three members are listed. Note that in Figure 2 several MapMan categories are combined.
A statistical analysis of GO term enrichment (p-value <0.05) was done using FatiGO  against the 15,983 unigenes of the Iberis transcriptome (Table 4C) or against the Arabidopsis inflorescence transcriptome downloaded from TAIR (Table 5).
The Arabidopsis stage 15 petal and inflorescence transcriptomic data originate from a microarray analysis (TAIR-ME00319; AtGenExpress: Expression Atlas of Arabidopsis Development) and were downloaded from TAIR (accessions numbers: HybData:1007126238–40 and HybData:1007126271–73). Only elements that were present in all three replicates were considered as being expressed, AGI codes appearing more than once were removed.
Quantitative real time PCR validation
Total RNA was isolated using the RNeasy Plant Mini kit (Qiagen, Valencia, CA), including an on-column DNA digestion step. Quantitative reverse transcriptase PCR was performed as described in . Primers were designed online with the program Primer3  and primer sequences are given in Additional file 8.
Specificity of each PCR was assessed by gel analysis and examination of the melt curve, generated after each reaction. The mean normalized expression (MNE) for each sample (calculated as in ) is based on three technical repetitions of each reaction, except for contigs 5537, 12149, 16003, 6640 and 7844 (Additional file 1 and 3), where for one single sample only two technical replicates were used. Error bars are standard error from technical replicates. The Iberis Ran3 GTPase [GenBank: EU145777] was used as a housekeeping gene for normalization.
Real time PCR was performed on adaxial and abaxial petals from stage A1 Iberis flowers on three biological replicates (replicates 5–7 in Figure 4D-I and Additional file 4A). The expression of Arabidopsis orthologs of Iberis genes to be tested was analysed using the Arabidopsis eFP browser at bar.utoronto.ca . In cases, where orthologs were strongly expressed in pollen, two biological replicates of pollen-free adaxial and abaxial petals from stage A1 flowers were used (replicates 3, 4 in Figure 4A-C). Expression differences of Arabidopsis orthologs of Iberis genes assigned to the functional GO-categories connected to cell wall and photosynthesis were verified on two biological replicates of (pollen-free) adaxial and abaxial petals of stage A1 Iberis flowers, which were dissected into the upper petal blade and the lower petal claw (replicates 1, 2 in Figure 3 and Additional file 4B and C).
For the validation of gene expression differences detected with the RNA-Seq approach, qPCR was performed on stage A1 petal pairs of at least two biological replicates. A pair wise expression comparison of a given gene was performed and the expression value of abaxial petals (or petal claws, in cases of blade/claw comparison) was set to 1. An expression increase in adaxial petals of 1.5-times or above in at least two biological replicates was considered as a differential expression and consequently as a further corroboration of RNA-Seq data.
Two technical samples of total RNA were prepared from each transgenic T2 Line expressing IaTCP1, TCP1 and CYC, respectively, under the control of the CaMV35S-promoter  or the empty vector pBAR35S [GenBank: AJ251014] as a reference. Total RNA was extracted from plants showing clear flower phenotypes in the main inflorescence, harvesting secondary inflorescences that carried only young, closed flower buds. RNA isolation was conducted using the RNA isolation kit NucleoSpin RNA XS (Macherey and Nagel). Probe preparation, hybridization to ATH1 Arabidopsis Genome Arrays (Affymetrix, Santa Clara, California) and statistical data analysis were carried out at the Integrated Functional Genomic service unit from the University of Münster (Germany). Only genes showing a minimal expression fold change of two and a p-value ≤0.05 were considered as being differentially expressed. Microarray expression data were determined in comparison to transgenic plants harbouring the empty vector. Data were corroborated for selected genes by reverse transcriptase PCR performed as in  (Additional file 9), using the same plant material as for microarray analysis. Values for transcript accumulation derive from three PCR repetitions with an annealing temperature of 58°C, normalized to the expression strength of Ran3 GTPase (At5g55190). Primer sequences are given in Additional file 8, microarray data is given in Additional files 6, 10, 7.
For in situ hybridization experiments an IaH4 [GenBank: EU145778] antisense probe was hybridized to longitudinal sections of I. amara flower buds before anthesis, carried out as described in .
For scanning electronic microscopy, petals were frozen in liquid nitrogen, sputter coated with palladium in a cryo-preparation system (K1250X, Emitech) and examined with a digital scanning electron microscope (Zeiss Auriga). Stage A1 and A2 petals pictures were taken with a binocular with integrated camera (Leica M165 FC, Wetzlar).
Availability of supporting data
The sequence datasets supporting the results of this article are available from the Short Read Archive (SRA) database at NCBI under the accession SRP048782. Assembled contigs are provided in FASTA format in Additional file 11. The microarray data have been deposited in the Gene Expression Omnibus (GEO) database at NCBI and are accessible under the accession GSE62213.
We are grateful to Claudia Gieshoidt for help with harvesting plant material, qPCR and plant cultivation. We are also much obliged to Dr. Kurt Stüber f nt Breeding Research in Cologne for his kind help with the BLAST analysis of Iberis contigs against an Arabidopsis protein database. This research was supported by the Deutsche Forschungsgemeinschaft.
- Dilcher D: Toward a new synthesis: major evolutionary trends in the angiosperm fossil record. Proc Natl Acad Sci U S A. 2000, 97 (13): 7030-7036. 10.1073/pnas.97.13.7030.PubMed CentralView ArticlePubMedGoogle Scholar
- Sargent RD: Floral symmetry affects speciation rates in angiosperms. Proc R Soc Land B Biol Sci. 2004, 271: 603-608. 10.1098/rspb.2003.2644.View ArticleGoogle Scholar
- Stebbins G: Flowering plants: evolution above the species level. 1974, Harvard University Press, Cambridge (MA)View ArticleGoogle Scholar
- Luo D, Carpenter R, Vincent C, Copsey L, Coen E: Origin of floral asymmetry in Antirrhinum . Nature. 1996, 383: 794-799. 10.1038/383794a0.View ArticlePubMedGoogle Scholar
- Luo D, Carpenter R, Copsey L, Vincent C, Clark J, Coen E: Control of organ asymmetry in flowers of Antirrhinum . Cell. 1999, 99: 367-376. 10.1016/S0092-8674(00)81523-8.View ArticlePubMedGoogle Scholar
- Howarth DG, Donoghue MJ: Phylogenetic analysis of the „ECE"(CYC/TB1) clade reveals duplications predating the core eudicots. Proc Natl Acad Sci. 2006, 103: 9101-9106. 10.1073/pnas.0602827103.PubMed CentralView ArticlePubMedGoogle Scholar
- Broholm SK, Tähtiharju S, Laitinen RAE, Albert VA, Teeri TH, Elomaa P: TCP domain transcription factor controls flower type specification along the radial axis of the Gerbera (Asteraceae) inflorescence. Proc Natl Acad Sci. 2008, 105 (26): 9117-9122. 10.1073/pnas.0801359105.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim M, Cui ML, Cubas P, Gillies A, Lee K, Chapman MA, Abbott RJ, Coen E: Regulatory genes control a key morphological and ecological trait transferred between species. Science. 2008, 322: 1116-1119. 10.1126/science.1164371.View ArticlePubMedGoogle Scholar
- Wang Z, Luo Y, Li X, Wang L, Xu S, Yang J, Weng L, Sato S, Tabata S, Ambrose M, Rameau C, Feng X, Hu X, Luo D: Genetic control of floral zygomorphy in pea (Pisum sativum L.). Proc Natl Acad Sci. 2008, 105: 10414-10419. 10.1073/pnas.0803291105.PubMed CentralView ArticlePubMedGoogle Scholar
- Busch A, Zachgo S: Control of corolla monosymmetry in the Brassicaceae Iberis amara . Proc Natl Acad Sci U S A. 2007, 104: 16714-16719. 10.1073/pnas.0705338104.PubMed CentralView ArticlePubMedGoogle Scholar
- Feng X, Zhao Z, Tian Z, Xu S, Luo Y, Cai Z, Wang Y, Yang J, Wang Z, Weng L, Chen J, Zheng L, Guo X, Luo J, Sato S, Tabata S, Ma W, Cao X, Hu X, Sun C, Luo D: Control of petal shape and floral zygomorphy in Lotus japonicas . Proc Natl Acad Sci. 2006, 103: 4970-4975. 10.1073/pnas.0600681103.PubMed CentralView ArticlePubMedGoogle Scholar
- Busch A, Horn S, Mühlhausen A, Mummenhoff K, Zachgo S: Corolla monosymmetry: evolution of a morphological novelty in the Brassicaceae Family. Mol Biol Evol. 2012, 29 (4): 1241-1254. 10.1093/molbev/msr297.View ArticlePubMedGoogle Scholar
- Costa MM, Fox S, Hanna AI, Baxter C, Coen C: Evolution of regulatory interactions controlling floral asymmetry. Development. 2005, 132: 5093-5101. 10.1242/dev.02085.View ArticlePubMedGoogle Scholar
- Krizek BA, Anderson JT: Control of flower size. J Exp Bot. 2013, 64: 1427-1437. 10.1093/jxb/ert025.View ArticlePubMedGoogle Scholar
- Weiss J, Degado-Benarroch L, Egea-Cortines M: Genetic control of floral size and proportions. Int J Dev Biol. 2005, 49: 513-525. 10.1387/ijdb.051998jw.View ArticlePubMedGoogle Scholar
- Bey M, Stüber K, Fellenberg K, Schwarz-Sommer Z, Sommer H, Saedler H, Zachgo S: Characterization of antirrhinum petal development and identification of target genes of the class B MADS box gene DEFICIENS . Plant Cell. 2004, 16: 3197-3215. 10.1105/tpc.104.026724.PubMed CentralView ArticlePubMedGoogle Scholar
- Vincent CA, Coen ES: A temporal and morphological framework for flower development in Antirrhinum majus . Can J Bot. 2004, 82: 681-690. 10.1139/b04-042.View ArticleGoogle Scholar
- Izhaki A, Borochov A, Zamski E, Weiss D: Gibberellin regulates post-microsporogenesis processes in petunia anthers. Physiol Plant. 2002, 115: 442-447. 10.1034/j.1399-3054.2002.1150314.x.View ArticlePubMedGoogle Scholar
- Martin C, Gerats T: Control of pigment biosynthesis genes during petal development. Plant Cell. 1993, 5: 1253-1264. 10.1105/tpc.5.10.1253.PubMed CentralView ArticlePubMedGoogle Scholar
- Cosgrove DJ: Expansive growth of plant cell walls. Plant Physiol Biochem. 2000, 38 (1/2): 109-124. 10.1016/S0981-9428(00)00164-9.View ArticlePubMedGoogle Scholar
- Ramirez-Parra E, Desvoyes B, Guiterrez C: Balance between cell division and differentiation during plant development. Int J Dev Biol. 2005, 49: 467-477. 10.1387/ijdb.052001er.View ArticlePubMedGoogle Scholar
- Anastasiou E, Kenz S, Gerstung M, MacLean D, Timmer J, Fleck C, Lenhard M: Control of plant organ size by KLUH/CYP78A5-dependant intercellular signaling. Dev Cell. 2007, 13: 843-856. 10.1016/j.devcel.2007.10.001.View ArticlePubMedGoogle Scholar
- Dinneny JR, Yadegari R, Fischer RL, Yanofsky MF, Weigel D: The role of JAGGED in shaping lateral organs. Development. 2004, 131: 1101-1110. 10.1242/dev.00949.View ArticlePubMedGoogle Scholar
- Hu Y, Xie Q, Chua NH: The Arabidopsis Auxin-inducible gene ARGOS controls lateral organ size. Plant Cell. 2003, 15: 1951-1961. 10.1105/tpc.013557.PubMed CentralView ArticlePubMedGoogle Scholar
- Mizukami Y, Fischer RL: Plant organ size control: AINTEGUMENTA regulates growth and cell numbers during organogenesis. Proc Natl Acad Sci U S A. 2000, 97: 942-947. 10.1073/pnas.97.2.942.PubMed CentralView ArticlePubMedGoogle Scholar
- Li Y, Zheng L, Corke F, Smith C, Bevan MW: Control of final seed and organ size by the DA1 gene family in Arabidopsis thaliana . Genes Dev. 2008, 22: 1331-1336. 10.1101/gad.463608.PubMed CentralView ArticlePubMedGoogle Scholar
- Disch S, Anastasiou E, Sharma VK, Laux T, Fletcher JC, Lenhard M: The E3 ubiquitin ligase BIG BROTHER controls Arabidopsis organ size in a dosage-dependant manner. Curr Biol. 2006, 16: 272-279. 10.1016/j.cub.2005.12.026.View ArticlePubMedGoogle Scholar
- Xu R, Li Y: Control of final organ size by mediator complex subunit 25 in Arabidopsis thaliana . Development. 2011, 138: 4545-4554. 10.1242/dev.071423.View ArticlePubMedGoogle Scholar
- Varaud E, Brioiudes F, Sz°Csi J, Leroux J, Brown S, Perrot-Rechenmann C, Bendahmanea M: Auxin response factor8 regulates Arabidopsis petal growth by interacting with the bHLH transcription factor BIGPETALp. Plant Cell. 2011, 23: 973-983. 10.1105/tpc.110.081653.PubMed CentralView ArticlePubMedGoogle Scholar
- Sz°Csi J, Joly C, Bordji K, Varaud E, Cock M, Dumas C, Bendahmane M: BIGPETALp, a bHLH transcription factor is involved in the control of Arabidopsis petal size. EMBO J. 2006, 25: 3912-3920. 10.1038/sj.emboj.7601270.View ArticleGoogle Scholar
- Hu Y, Poh HM, Chua N-H: The Arabidopsis ARGOS-LIKE gene regulates cell expansion during organ growth. Plant Jnl. 2006, 47: 1-9. 10.1111/j.1365-313X.2006.02750.x.View ArticleGoogle Scholar
- Guan S, Lu Y: Dissecting organ-specific transcriptomes through RNA-sequencing. Plant Methods 2013, 9:42.,Google Scholar
- Bräutigam A, Kajala K, Wullenweber J, Sommer M, Gagneul D, Weber KL, Carr KM, Gowik U, Mass J, Lercher MJ, Westhoff P, Hibberd JM, Weber APM: An mRNA blueprint for C4 photosynthesis derived from comparative transcriptomics of closely related C3 and C4 species. Plant Physiol. 2011, 155: 142-156. 10.1104/pp.110.159442.PubMed CentralView ArticlePubMedGoogle Scholar
- Sharma N, Bhalla PL, Singh MB: Transcriptome-wide profiling and expression analysis of transcription factor families in a liverwort, Marchantia polymorpha . BMC Genomics 2013, 14:915.,Google Scholar
- Thimm O, Bläsing O, Gibon Y, Nagel A, Meyer S, Krüger P, Selbig J, Müller LA, Rhee SY, Stitt M: MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004, 37: 914-939. 10.1111/j.1365-313X.2004.02016.x.View ArticlePubMedGoogle Scholar
- Provart N, Zhu T: A browser-based functional classification SuperViewer for Arabidopsis genomics. Curr Comput Mol Biol 2003, 271-272.,Google Scholar
- Smyth D, Bowman J, Meyerowitz E: Early flower development in Arabidopsis . Plant Cell. 1990, 2: 755-767. 10.1105/tpc.2.8.755.PubMed CentralView ArticlePubMedGoogle Scholar
- Al-Shahrour F, Dfaz-Uriarte R, Dopazo J: FatiGO: a web tool for finding significant associations of gene ontology terms with groups of genes. Bioinformatics. 2004, 20: 578-580. 10.1093/bioinformatics/btg455.View ArticlePubMedGoogle Scholar
- Cao J: The pectine lyases in Arabidopsis thaliana: evolution, selection and expression profiles. PLoS One. 2012, 7 (10): e46944-10.1371/journal.pone.0046944.PubMed CentralView ArticlePubMedGoogle Scholar
- Kanneganti V, Gupta AK: Wall associated kinases form plants – an overview. Physiol Mol Biol Plants. 2008, 14 (1&2): 109-118. 10.1007/s12298-008-0010-6.PubMed CentralView ArticlePubMedGoogle Scholar
- Cao J, Shi F: Evolution of the RALF gene family in plants: gene duplication and selection patterns. Evol Bioinformatics. 2012, 8: 271-292. 10.4137/EBO.S9652.View ArticleGoogle Scholar
- Pearce G, Moura DS, Stratmann J, Ryan CA: RALF, a 5-kDa ubiquitous polypeptide in plants arrests root growth and development. Proc Natl Acad Sci U S A. 2001, 98: 12843-12847. 10.1073/pnas.201416998.PubMed CentralView ArticlePubMedGoogle Scholar
- De Craene J-O, Ripp R, Lecompte O, Thompson JD: Evolutionary analysis of the ENTH/ANTH/VHS protein superfamily reveals a coevolution between membrane trafficking and metabolism. BMC Genomics 2012, 13:297.,Google Scholar
- Schindelmann G, Morikami A, Jung J, Baskin TI, Carpita NC, Derbyshire P, McCann MC, Benfey PN: COBRA encodes a putative GPI-anchored protein, which is polarly localized and necessary for oriented cell expansion in Arabidopsis . Genes Dev. 2001, 15: 1115-1127. 10.1101/gad.879101.View ArticleGoogle Scholar
- Li S, Ge F-R, Xu M, Zhao X-Y, Huang G-Q, Zhou L-Z, Wang J-G, Kombrink A, McCormick S, Zhang XS, Zhang Y: Arabidopsis COBRA-LIKE 10, a GPI-anchored protein, mediates directional growth of pollen tubes. Plant J. 2013, 74: 486-497. 10.1111/tpj.12139.View ArticlePubMedGoogle Scholar
- Geisler M, Murphy AS: The ABC of auxin transport: the role of p-glycoproteins in plant development. FEBS Lett. 2006, 580: 1094-1102. 10.1016/j.febslet.2005.11.054.View ArticlePubMedGoogle Scholar
- Bartina I, Otto E, Strnad M, Werner T, Schmülling T: Cytokinin regulates the activity of reproductive meristems, flower organ size, ovule formation and thus seed yield in Arabidopsis thaliana . Plant Cell. 2011, 23: 69-80. 10.1105/tpc.110.079079.View ArticleGoogle Scholar
- Hiroi N, Ito T, Yamamoto H, Ochiya T, Jinno S, Okayama H: Rcd1 is a novel transcriptional cofactor that mediates retinoic acid-induced cell differentiation. EMBO. 2002, 21 (19): 5235-5244. 10.1093/emboj/cdf521.View ArticleGoogle Scholar
- Li X, Luo J, Yan T, Xiang L, Feng J, Qin D, Sun C, Xie M: Deep sequencing-based analysis of the Cymbidium ensifolium floral transcriptome. PLoS One. 2013, 8 (12): e85480-10.1371/journal.pone.0085480.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu Z, Ma L, Nan Z, Wang Y: Comparative transcriptional profiling provides insights into the evolution and development of the zygomorphic flower of Vicia sativa (Papilionoideae). PLoS One. 2013, 8 (2): e57338-10.1371/journal.pone.0057338.PubMed CentralView ArticlePubMedGoogle Scholar
- Soeraert SSA, Neste CMF, Vandewoestyne ML, Head SRm Goossens A, van Nieuwerburgh CW, Degorce DLD: Differential transcriptome analysis of glandular and filamentous trichomes in Artemisia annua . BMC Plant Biol 2013, 13:220.,Google Scholar
- Zik M, Irish VF: Global identification of target genes regulated by APETALA3 and PISTILLATA floral homeotic gene action. Plant Cell. 2003, 15: 207-222. 10.1105/tpc.006353.PubMed CentralView ArticlePubMedGoogle Scholar
- Ginglinger J-F, Boachon B, Höfer R, Paetz C, Köllner TG, Miesch L, Lugan R, Baltenweck R, Mutterer J, Ullmann P, Beran F, Claudel P, Verstappen F, Fischer MJC, Karst F, Bouwmeester H, Miesch M, Schneider B, Gershenzon J, Ehlting J, Werck-Reichhart D: Gene coexpression analysis reveals complex metabolism of the monoterpene alcohol linalool in Arabidopsis flowers. Plant Cell. 2013, 25: 4640-4657. 10.1105/tpc.113.117382.PubMed CentralView ArticlePubMedGoogle Scholar
- Herken EN, Celik A, Aslan M, Aydinlik N: The constituents of essential oils: antimicrobial and antioxidant activity of Micromeria congesta Boiss. & Hausskn. Ex Boiss. From East Anatolia. J Med Food. 2012, 15 (9): 835-839. 10.1089/jmf.2011.0315.View ArticlePubMedGoogle Scholar
- Schommer C, Palatnik JF, Aggarwal P, Chételat A, Cubas P, Farmer EE, Nath U, Weigel D: Control of jasmonate biosynthesis and senescence by miR319 targets. PLoS Biol. 2008, 6 (9): e230-10.1371/journal.pbio.0060230.PubMed CentralView ArticlePubMedGoogle Scholar
- Bell E, Creelmann RA, Mullet JE: A chloroplast lipoxygenase is required for wound-induced jasmonic acid accumulation in Arabidopsis . Proc Natl Acad Sci. 1995, 92: 8675-8679. 10.1073/pnas.92.19.8675.PubMed CentralView ArticlePubMedGoogle Scholar
- Pauwels L, Morreel K, De Witte E, Lammertyn F, van Montagu M, Boerjan W, Inzé D, Goossens A: Mapping methyl jasmonate-mediated transcriptional reprogramming of metabolites and cell cycle progression in cultured Arabidopsis cells. Proc Natl Acad Sci U S A. 2008, 105: 1380-1385. 10.1073/pnas.0711203105.PubMed CentralView ArticlePubMedGoogle Scholar
- Lampugnani ER, Kilinc A, Smyth DR: PETAL LOSS is a boundary gene that inhibits growth between developing sepals in Arabidopsis thaliana . Plant J. 2012, 71: 724-735. 10.1111/j.1365-313X.2012.05023.x.View ArticlePubMedGoogle Scholar
- Lampugnani ER, Kilinc A, Smyth DR: Auxin controls petal initiation in Arabidopsis . Development. 2013, 140: 185-194. 10.1242/dev.084582.View ArticlePubMedGoogle Scholar
- Noh B, Murphy AS, Spalding EP, Cosgrove DJ: Multidrug resistance-like genes of Arabidopsis required for auxin transport and auxin-mediated development. Growth Plant Cell Wall Nat Rev. 2005, 6: 850-861.Google Scholar
- Cosgrove DJ: Growth of the plant cell wall. Nat Rev. 2005, 6: 850-861. 10.1038/nrm1746.View ArticleGoogle Scholar
- Peaucell A, Braybrook S, Höfte H: Cell wall mechanics and growth control in plants: The role of pectins revisited. 2012, Frontiers in, Plant SciGoogle Scholar
- Campbell P, Braam J: Xyloglucan endotransglycosylases diversity of genes, enzymes and potential wall-modifying functions. Trends Plant Sci. 1999, 4 (9): 361-366. 10.1016/S1360-1385(99)01468-5.View ArticlePubMedGoogle Scholar
- Urbanowicz BR, Bennett AB, del Campillo E, Catalá C, Hayashi T, Henrissat B, Höfte H, McQueen-Mason SJ, Patterson SE, Shoseyov O, Teeri TT, Rose JKC: Structural organization and a standardized nomenclature for plant endo-1-4-ß-glucanases (cellulases) of glycosyl hydrolase family 9. Plant Physiol. 2007, 144: 1693-1696. 10.1104/pp.107.102574.PubMed CentralView ArticlePubMedGoogle Scholar
- Sato S, Kato T, Kakegawa K, Ishii T, Liu Y-G, Awano T, Takabe K, Nishiyama Y, Kuga S, Sato S, Nakamura Y, Tabata S, Shibata D: Role of the putative membrane-bound endo-1,4-glucanase KORRIGAN in cell elongation and cellulose synthesis in Arabidopsis thaliana . Plant Cell Physiol. 2001, 42 (3): 251-263. 10.1093/pcp/pce045.View ArticlePubMedGoogle Scholar
- Mingossi FB, Matos JL, Rizzato AP, Medeiros AH, Falco MC, Silva-Filho MC, MourA DS: SacRALF1, a peptide signal from the grass sugarcane (Saccharum spp.) is potentially involved in the regulation of tissue expansion. Plant Mol Biol. 2010, 73: 271-281. 10.1007/s11103-010-9613-8.View ArticlePubMedGoogle Scholar
- Covey PA, Subbaiah CC, Parsons RL, Pearce G, Lay FT, Anderson MA, Ryan CA, Bedinger PA: A pollen-specific RALF from tomato that regulates pollen tube elongation. Plant Physiol. 2010, 153: 703-715. 10.1104/pp.110.155457.PubMed CentralView ArticlePubMedGoogle Scholar
- Reale L, Porcessu A, Lanfaloni L, Moretti C, Zenoni S, Pezzotti M, Romano B, Ferranti F: Patterns of cell division and expansion in developing petals of Petunia hybrida . Sex Plant Reprod. 2002, 15: 123-132. 10.1007/s00497-002-0150-8.View ArticleGoogle Scholar
- Crawford BCW, Nath U, Carpenter R, Coen ES: CINCINNATA controls both cell differentiation and growth in petal lobes and leaves of Antirrhinum. Plant Physiol. 2004, 135: 244-253. 10.1104/pp.103.036368.PubMed CentralView ArticlePubMedGoogle Scholar
- Kosugi S, Ohashi Y: DNA binding and dimerization specificity and potential targets for the TCP protein family. Plant J. 2002, 30 (3): 337-348. 10.1046/j.1365-313X.2002.01294.x.View ArticlePubMedGoogle Scholar
- Peaucelle A, Louvet R, Johansen JN, Salsac F, Morin H, Fournet F, Belcram K, Gillet F, Höfte H, Laufs P, Mouille G, Pelloux J: The transcription factor BELLRINGER modulates phyllotaxis by regulating the expression of a pectin methylesterase in Arabidopsis . Development. 2011, 138: 4733-4741. 10.1242/dev.072496.View ArticlePubMedGoogle Scholar
- Rutjens B, Bao D, van Eck-Stouten E, Brand M, Smeekens S, Proveniers M: Shoot apical meristem function in Arabidopsis requires the combined activities of the three BEL1-like homeodomain proteins. Plant J. 2009, 58: 641-654. 10.1111/j.1365-313X.2009.03809.x.View ArticlePubMedGoogle Scholar
- Smith HMS, Hake S: The interaction of two homeobox genes BREVIPEDICELLUS and PENNYWISE, regulates internode patterning in the Arabidopsis inflorescence. Plant Cell. 2003, 15: 1717-1727. 10.1105/tpc.012856.PubMed CentralView ArticlePubMedGoogle Scholar
- Baggerly KA, Deng L, Morris JS, Aldaz CM: Differential expression in SAGE: Accounting for normal between-library variation. Bioinformatics. 2003, 19 (12): 1477-1483. 10.1093/bioinformatics/btg173.View ArticlePubMedGoogle Scholar
- Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2000, Humana Press, Totowa (NJ), 365-386.Google Scholar
- Winter D, Vinegar B, Nahal H, Ammar R, Wilson GV, Provart NJ: An „Electronic Fluorescent Pictograph"Browser for exploring and analyzing large-scale biological data sets. PLoS One. 2007, 2 (8): e718-10.1371/journal.pone.0000718.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.