RNAi-mediated silencing of the HD-Zip gene HD20 in Nicotiana attenuata affects benzyl acetone emission from corollas via ABA levels and the expression of metabolic genes
© Ré et al.; licensee BioMed Central Ltd. 2012
Received: 3 March 2012
Accepted: 1 May 2012
Published: 1 May 2012
The N. attenuata HD20 gene belongs to the homeodomain-leucine zipper (HD-Zip) type I family of transcription factors and it has been previously associated with the regulation of ABA accumulation in leaves and the emission of benzyl acetone (BA; 4-phenyl-2-butanone) from night flowers. In this study, N. attenuata plants stably reduced in the expression of HD20 (ir-hd20) were generated to investigate the mechanisms controlling the emission of BA from night flowers.
The expression of HD20 in corollas of ir-hd20 plants was reduced by 85 to 90% compared to wild-type plants (WT) without affecting flower morphology and development. Total BA emitted from flowers of ir-hd20 plants was reduced on average by 60%. This reduction occurred mainly at the late phase of BA emission and it was correlated with 2-fold higher levels of ABA in the corollas of ir-hd20 plants. When a 2-fold decline in ABA corolla levels of these plants was induced by salt stress, BA emissions recovered to WT levels. Supplying ABA to WT flowers either through the cuticle or by pedicle feeding reduced the total BA emissions by 25 to 50%; this reduction occurred primarily at the late phase of emission (similar to the reduction observed in corollas of ir-hd20 plants). Gene expression profiling of corollas collected at 12 pm (six hours before the start of BA emission) revealed that 274 genes changed expression levels significantly in ir-hd20 plants compared to WT. Among these genes, more than 35% were associated with metabolism and the most prominent group was associated with the metabolism of aromatic compounds and phenylpropanoid derivatives.
The results indicated that regulation of ABA levels in corollas is associated with the late phase of BA emission in N. attenuata plants and that HD20 affects this latter process by mediating changes in both ABA levels and metabolic gene expression.
Many plant species emit floral scents as long-distance attraction cues for pollinators, in particular moths, that search and visit flowers at night. The attraction of moths and other animal pollinators have important implications for plant reproductive success; floral scents are important in determining seed or fruit set in non-selfing plants and the frequency of outcross in selfing plants [1–4]. Thus, the biosynthesis and release of floral scents have important ecological and agricultural implications and the molecular mechanisms underlying the regulation of floral scent production are now beginning to be unraveled [5, 6].
Floral scents are usually complex mixtures of small volatile molecules and the most prevalent compounds in these mixtures are monoterpenoids, sesquiterpenoids, phenylpropanoids, benzenoid compounds and fatty acid derivatives . Volatile compounds can be either emitted from several parts of the flower or from a specific floral part . The expression of genes encoding floral scent biosynthetic enzymes is temporally and spatially regulated during flower development [1, 8]. In most cases studied so far, the expression of these genes correlates with the emission of the corresponding volatile compound indicating that these volatiles are de novo synthesized. Thus, scent production and emission are commonly regulated at the transcriptional level. Similar observations have been made in regard to the induction of genes involved in floral pigment production .
Nicotiana attenuata, a night-flowering tobacco that germinates after fires in the Southwestern United States, normally produces flowers that open at night and release benzyl acetone (BA; 4-phenyl-2-butanone) to attract night-active hawk moth pollinators (Manduca quinquemaculata and M. sexta). N. attenuata is a fully self-compatible species however it produces more than 30% of its seed from opportunistic out-crossing performed by pollinators [2, 11]. BA is the main constituent of the N. attenuata floral bouquet  and it begins to be emitted from the corolla limb in the evening [10, 12]. The emission of BA is synchronized with the development of the flower and it starts as the corolla limb opens . At present, the mechanisms that control the biosynthesis and release of BA from corollas are largely unknown. The BA biosynthesis pathway remains elusive, however, the chemical backbone may derived from the shikimate pathway  as N. attenuata plants with reduced expression of CHALCONE SYNTHASE 1 (CHAL1) are deficient in BA emission .
The homeodomain-leucine zipper (HD-Zip) family of transcription factors (TFs) is a plant-specific family of TFs in which the HD and Zip domains are combined in a single polypeptide . The HD is responsible for the specific binding to DNA and the Zip domain acts as a dimerization motif; dimerization is a prerequisite for binding to the phylogenetically conserved target sequence CAAT (A/T) ATTG [14, 15]. The HD-Zip family can be divided into four subfamilies (I to IV) according to the sequence similarity of the HD and Zip domains and to additional structural features outside these domains. The function of HD-Zips type-I has been associated to the regulation of development in response to changes in the environment [13, 16–26].
The N. attenuata HD20 gene was originally identified as an HD-Zip type-I whose expression is induced by multiple stress-associated stimuli including drought and wounding . In a previous study, we reduced the expression of HD20 by virus induced gene silencing (VIGS) and we demonstrated that HD20 plays not only a positive role in ABA accumulation in leaves during water stress but also in the emission of benzyl acetone (BA) from night flowers . However, the HD20-dependent mechanisms underlying the reduction in BA emission were not studied. Here we investigated these mechanisms by generating transgenic N. attenuata plants stably reduced in the expression of HD20. These plants were characterized by using a combination of volatile, phytohormone and gene expression profiling approaches to investigate in more detail the process of BA emission.
HD20 mRNA expression in corolla during development
Transformation of N. attenuata plants to reduce HD20 expression
HD20 positively regulates benzyl acetone emission from corollas
Changes in ABA levels correlate with changes in BA emission during corolla opening
To investigate whether the changes in ABA levels quantified in corollas of ir-hd20 plants were specific to this tissue and not the effect of altered ABA accumulation in the whole flower during development, phytohormone levels were quantified in whole flowers (corollas plus pistils, stamens, ovaries and nectaries) at progressive developmental stages. The results showed that accumulation of ABA in whole flowers was similar between WT and ir-hd20 plants (Additional file 4) and hence the changes in ABA accumulation in corollas of ir- hd20 plants do not result from a general alteration of whole flower ABA levels.
To further investigate this association, N. attenuata WT and ir-hd20 plants were grown under salt stress, a condition that maintained ABA in corollas at levels lower than 20 μg gFW-1 and almost completely abrogated the developmental fluctuation of this phytohormone in corollas of WT plants (Figure 4a). The salt stress did not affect BA emission from corollas of WT plants (compared to control growth conditions; Figure 4b). However, salt stress conditions recovered BA emission in ir-hd20 plants and the levels of emitted BA were similar to those in WT plants (Figure 4b). From these results, we hypothesized that high levels of ABA in corollas (i.e. 20 μg gFW-1 or more) have a negative effect on BA emission and that the developmentally controlled 2-fold reduction in ABA levels in corollas of WT plants contributes to the emission of this volatile. To test this hypothesis, exogenous ABA was supplied to flowers of WT plants and the level of BA emitted from the corollas was quantified. Two different approaches were used: (1) flowers attached to the plant were sprayed with 1 mL of a solution containing 1 μg mL-1 ABA at 8 pm, and (2) clipped flowers were placed in scintillation vials and continuously fed through the pedicle with 100 μL of a solution containing 1 μg mL-1 ABA. The levels of ABA used in these experiments were in the range of endogenous ABA levels quantified in corollas (for spraying, we estimated a 90% loss of solution by measuring the run off). In the first approach, the emitted BA was trapped from 3 pm to 9 am (18 h) and analyzed by GC-MS (Figure 4c). In the second approach, the emitted BA was analyzed every two hours starting at 6 pm with a zNOSETM (Figure 4d). Corollas from WT plants sprayed with ABA emitted 30% less BA than control-treated flowers during the night (Figure 4c). Consistently, continuously pedicle-fed flowers emitted 25 to 50% lower levels of BA at 12 am and 2 am compared to control-treated flowers (Figure 4d).
In summary, these experiments indicated that ABA levels in corollas contributed either directly or indirectly to the release of BA. The negative effect of ABA on BA emission was however only significant after 10 pm, suggesting that the developmentally controlled reduction in ABA levels in corollas of WT plants contributes to the late phase of BA release. This conclusion is consistent with the reduction in ABA levels after 8 pm in corollas of WT plants (Figure 4a).
Reduced levels of HD20 in corollas have strong effects on gene expression
List of selected genes involved in biochemical pathways
EC:188.8.131.52 - unspecific monooxygenase
EC:184.108.40.206 - L-amino-acid oxidase
EC:220.127.116.11 - aromatic-L-amino-acid decarboxylase
EC:18.104.22.168 - L-amino-acid oxidase
EC:22.214.171.124 - aromatic-L-amino-acid decarboxylase
EC:126.96.36.199 - tyrosine decarboxylase
EC:188.8.131.52 - L-amino-acid oxidase
EC:184.108.40.206 - aromatic-L-amino-acid decarboxylase
EC:220.127.116.11 - dihydroflavonol 4-reductase
EC:18.104.22.168 - chalcone synthase
EC:22.214.171.124 - leucocyanidin oxygenase
Starch and sucrose
EC:126.96.36.199 - starch synthase
EC:188.8.131.52 - fructokinase
EC:184.108.40.206 - alpha,alpha-trehalose-phosphate synthase (UDP-forming)
EC:220.127.116.11 - trehalose-phosphatase
Cell wall biosynthesis
EC:18.104.22.168 - cellulose synthase (UDP-forming)
EC:22.214.171.124 - cellulose synthase (UDP-forming)
EC:126.96.36.199 - phosphoethanolamine N-methyltransferase
EC:188.8.131.52 - glycerophosphodiester phosphodiesterase
Fatty acid biosynthesis
EC:184.108.40.206 - fatty-acyl-CoA synthase
List of selected genes involved in regulatory processes
GBF S pro-rich region-interacting factor 1
WRKY transcription factor
Fungal-specific transcription factor domain protein
Histone acetyltransferase complex component
RWP-RK domain-containing protein
Chromatin remodeling complex subunit
WRKY transcription factor 29
RNA polymerase sigma factor
Paired amphipathic helix protein SIN3-like 2
GRAS transcription factor
Protein kinases and phosphatases
Ser-Thr protein phosphatase
Protein phosphatase 2 C
TCTR2 protein kinase
Ethylene receptor 1 (ETR1)
S-locus lectin protein kinase
Coronatine insensitive 1 (COI1)
LRR receptor-like kinase
Other regulatory activities
circadian clock-associated FKF1
The expression of two regulatory genes was also analyzed, namely Casein Kinase (CK; the most strongly down-regulated gene; Additional file 5 and Table 2) and Protein Phosphatase 2 C (PP2C). The expression of CK was induced several fold during corolla opening in WT plants and reached a maximum at 8 pm (Figure 6e). In contrast, in corollas of ir-hd20 plants, the levels of this transcript remained largely non-induced (Figure 6e; univariate ANOVA, F 2,8 and P < 0.05 followed by a LSD post-hoc test P < 0.05 for both lines versus WT). Finally, the levels of PP2C mRNA were induced by 2-fold during corolla opening in WT plants and the levels remained approximately 2-fold lower in corollas of ir-hd20 plants (Figure 6f; univariate ANOVA, F 2,8 and P < 0.05 followed by a LSD post-hoc test P < 0.05 for both lines versus WT).
Changes in ABA levels in opening corollas are associated with BA emission
The expression of HD20 in flowers and corollas was consistent with the participation of this transcription factor in the process of BA emission. The levels of HD20 mRNA were highest in corollas and they increased during the time of corolla opening (Figure 1). The timing of induction correlated with the timing of BA emission which starts in the afternoon and peaks between 9 and 10 pm (Figure 3b). Several HD-Zip type-I transcription factors in different plants species are either highly expressed in flowers or have been associated with mechanisms regulating organ development. For example, LeHB-1 from tomato (Solanum lycopersicum) regulates floral organogenesis  and the Arabidopsis HB12 HB7 HB6 and HB1 are highly expressed in flowers [19, 20]. Based on phylogenetic analysis, the Arabidopsis HB12 and HB7 genes group together with HD20 in subgroup Ic [27, 33].
Consistent with our previous study in which the expression of HD20 was reduced by VIGS in N. attenuata plants , plants with stably reduced expression of HD20 were also reduced in the levels of BA released from corollas (Figure 3). Analysis of BA emission during the first and second nights showed that flowers from ir-hd20 plants produced less BA with no delay in its emission (Figure 3a). Moreover, analysis of the dynamics of BA emission showed that corollas of ir-hd20 plants released significantly smaller BA amounts after 8 pm, indicating that reduced HD20 expression affected primarily the late phase of BA emission (Figure 3b).
In WT plants, HD20 mRNA levels were detected throughout the time of corolla opening (i.e. from 12 pm to 12 am) and peaked at 8 pm (Figure 1b). ABA levels increased in corollas of WT plants from 12 pm to 4 pm and decrease from 8 pm to 9 am (Figure 4a). In corollas of ir-hd20 plants, this fluctuation in ABA levels was not observed (Figure 4a) however when ABA levels in corollas were reduced by salt stress treatment, the emission of BA was recovered to WT levels (Figure 4b). Moreover, feeding experiments with exogenous ABA decreased BA emission at the late phase of the process (i.e. after 10 pm; Figures 4c and d). Based on these observations, we concluded that there was a direct association between the decline of ABA levels in corolla and the late phase of BA emission. Thus, the results suggested that HD20 contributes partially to BA emission by directly or indirectly affecting ABA levels. How HD20 affects ABA levels in corollas is at present unknown. Genome-wide gene expression analysis performed with corolla tissue harvested at 12 pm (6 h before the start of BA emission) did not detect genes involved in ABA biosynthesis as differentially expressed. At present, however, it is unknown how the corolla controls ABA levels, this could be the result of de novo biosynthesis, import from adjacent tissue, metabolism (i.e. turnover) or a combination of these processes. Although at 12 pm, 2-fold higher levels of ABA were detected in corollas of ir-hd20 plants compared to WT (Figure 4a), only one gene homolog to Zea mays SnRK2 (known to be regulated by ABA ) was identified as differentially expressed (supplementary file 5). Little is known about the regulation of gene expression by ABA in corollas and ABA is known to affect multiple cellular processes (e.g., opening of ion channels, protein phosphorylation) in addition to gene expression . Moreover, the direct target genes of HD20 are thus far unknown in N. attenuata. The mechanisms underlying the reduction in BA emission mediated by HD20 and ABA are probably very complex and additional experimentation will be required to disentangle this complexity.
Changes in metabolic gene expression induced by HD20 are associated with BA emission
Chalcone Synthase 1 (CHAL1) is so far the only gene that has been involved in the biosynthesis of BA in N. attenuata plants . In other words, the biosynthesis pathway for this floral volatile is largely unknown. Chalcone synthases catalyze the first step of flavonoid biosynthesis using malonyl-coenzyme A (CoA) and 4-coumaroyl-CoA to produce naringenin chalcone . 4-coumaroyl-CoA is supplied by the phenylpropanoid pathway which utilizes phenylalanine produced by the shikimate pathway in the plastid. Tyrosine and tryptophan are also produced by this biosynthesis pathway. Naringenin chalcone can be further modified to form flavonoids and dihydroflavonol 4-reductase (DFR) is the first committed enzyme of anthocyanin biosynthesis in the flavonoid pathway . Benzoids are also derived from the phenylpropanoid or shikimate pathways and enzymes such as S-adenosyl-L-methionine:benzoic acid/Salicylic acid carboxyl methyltransferase (BAMT/SAMT) synthesize methyl esters (e.g., methyl benzoate and methyl salicylate) which are constituents of floral scents of many plant species . In Petunia hybrida, where benzenoids determine the floral scent, three transcription factors from the R2R3-type MYB family have been identified that regulate floral scent production. ODORANT1 (ODO1), EMISSION OF BENZENOIDS II (EOBII) and PhMYB4 affect the expression of several biosynthetic floral scent-related genes encoding enzymes from the phenylpropanoid and shikimate pathways [5, 6, 37].
The analysis of changes in gene expression in opening corollas of ir-hd20 plants showed that the most prevalent group of genes with altered expression corresponded to those involved in aromatic compound metabolism and derivatives of the phenylpropanoid biosynthesis pathway. Importantly, no candidate genes involved in BA biosynthesis could be identified, however, as the genes involved in this process are largely unknown and a large fraction of the genes differentially expressed in corollas of ir-hd20 are of undescribed function (supplementary file 5), it is plausible that among these genes some involved in BA biosynthesis are present.
The induction of expression of genes involved in aromatic compound metabolism and derivatives of the phenylpropanoid biosynthesis pathway during corolla opening (e.g., CHAL1, SAMT/BAMT, DFR) correlated positively with the emission of BA and with the turning of the corolla limb from green to white. Thus, these observations were consistent with the altered levels of BA production in ir-hd20 plants and suggested that reduced emission of this floral volatile could be the result of altered gene expression in BA biosynthesis genes. Paradoxically, however, in addition to CHAL1, several other genes involved in different steps of aromatic compound biosynthesis or metabolism were up-regulated in corollas of ir-hd20 plants (Table 1). Because several genes involved in carbohydrate metabolism were down-regulated in corollas of these plants (Table 1), one plausible explanation for the reduced levels of BA emission in ir-hd20 plants is that reduced carbon fluxes through the shikimate and phenylpropanoid pathways limit BA production. If this is the case, the increased levels of chalcone synthase and other phenylpropanoid biosynthesis genes observed in ir-hd20 plants may reflect a compensatory mechanism for reduced supply of carbon precursors. A detailed study of the metabolic fluxes of the shikimate and phenylpropanoid pathways in ir-hd20 plants will be necessary to test this hypothesis in the future.
In addition to metabolic genes, multiple regulatory factors also showed altered expression in opening corollas of ir-hd20 plants (Table 2). Among the genes in the group of protein kinases and phosphatases were N. attenuata homologues of casein kinase and PP2C (Table 2). Casein kinases (CKs) are ubiquitous Ser/Thr kinases that play critical roles in all higher organisms including plants. For example, the Arabidopsis CK2 plays important roles in light signaling, circadian rhythms, hormone responses, cell cycle control and flowering time [6, 38, 39]. The strong down-regulation of casein kinase suggests that this protein kinase may be important in regulating directly or indirectly BA emissions by HD20-dependent mechanisms. Interestingly, similar to HB12 and HB7 in Arabidopsis , the carboxyl terminus of HD20 presents three putative phosphorylation sites (data not shown), suggesting that phosphorylation/dephosphorylation mechanisms can be important for the regulation of its activity.
In connection with ABA signaling, a homologue of PP2C was down-regulated approximately two-fold in opening corollas of ir-hd20 plants (Table 2 and Figure 6f). PPC2s such as ABI1 and ABI2 (Abscisic Acid Insensitive 1 and 2, respectively) are rapidly inactivated upon the binding of ABA to RCARs/PYR1/PYLs [5, 40, 41]. Inactivation of PP2Cs triggers the activation of SNF1-type kinases that initiate ABA-dependent responses such as the activation of gene expression and the regulation of ion channels [5, 40, 41]. As their names indicate, mutations in ABI1 and ABI2 genes make plants insensitive to ABA . Functional ABI1 and ABI2 alleles are necessary for the induction of the Arabidopsis HB12 and HB7 by ABA . Thus, in addition to the de-regulation of ABA levels in corollas of ir-hd20 plants, sensitivity to this phytohormone could also be affected in corollas of these plants.
Other interesting genes with receptor activity that changed expression in opening corollas of ir-hd20 compared to WT plants were Ethylene receptor 1 (ETR1), Coronatine insensitive 1 (COI1) and the circadian clock-associated Flavin-binding Kelch repeat F-box 1 (FKF1) gene . Ethylene and JA-Ile are important phytohormones for flower development. For example, plants deficient in the expression of COI1 are sterile [31, 32, 44] and plants expressing a constitutively active allele of ETR1 have delayed corolla senescence [45, 46]. Although JA, JA-Ile and ethylene levels were not affected in corollas of ir-hd20 plants, changes in the sensitivity to these phytohormones may also participate in the regulation of BA emission.
In summary, from the results presented in this study it is clear that multiple metabolic pathways were affected when corollas were reduced in HD20 expression, pointing to the importance of this HD-Zip transcription factor in the regulation of biochemical processes in this tissue. Importantly, reduced expression of HD20 did not cause morphological and developmental changes in flowers, indicating that the main function of this transcription factor in this organ relates to the regulation of metabolism rather than development. The mechanisms underlying the regulation of BA levels by HD20 are clearly complex and they most likely involve a complex network of factors and ABA signaling. As the sequences of the promoter regions of N. attenuata genes become available with the sequencing of this plant’s genome, the analysis of the presence of the universally conserved HD-Zip type I binding element (CAAT(A/T)ATTG)[14, 15] in genes affected in their expression in ir-hd20 plants will facilitate the discovery of direct target genes of HD20.
Plant growth and treatments
Seeds of N. attenuata plants were germinated on agar plates containing Gamborg's B5 medium as previously described . Plates were maintained in a growth chamber (Snijders Scientific, Tilburg, Netherlands) at 26°C/16 h (155 μmol s-1 m-2 light), 24°C/8 h dark for 10 days. Ten-day old seedlings were transferred to TEKU pots (Pöppelmann GmbH & Co. KG, Lohne, Germany) with Klasmann plug soil (Klasmann-Deilmann GmbH, Geesten, Germany). After 10 days, seedlings were transferred to soil in 2 L pots and grown in the glasshouse under high-pressure sodium lamps (200-300 μmol s−1 m−2) with a day/night ratio of 16 h (26–28°C)/8 h (22–24°C) and 45–55% humidity. Plants were grown in pairs (one WT and one ir-hd20 per 2-L pot).
For salt stress treatments, elongated N. attenuata plants were watered daily under a regime of: 2 days with an aqueous solution of 75 mM NaCl, 3 days with 150 mM NaCl and 4 days with 300 mM NaCl. The plants were left to soak the saline water for 1 h and the excess of liquid was removed from the trays. The saline water volume was 0.6 L per tray and 5 plant pots were placed per tray.
For exogenous ABA application, a solution of 1 μg mL-1 ABA in 0.02% (v/v) Tween-20/water was sprayed (1 mL per flower) on corollas at 8 pm with a perfume dispenser. The control solution was 0.02% (v/v) Tween-20/water. Flowers were immediately enclosed in volatile traps for BA analysis by GC-MS (see below). For pedicle feeding of ABA, individual clipped flowers were placed inside 1.7 mL microcentrifuge tubes (with their caps removed) containing 100 μL of either water or ABA (1 μg mL-1) dissolved in water. The microcentrifuge tubes carrying the flowers were placed inside glass scintillation vials and the vials were screw-capped. BA levels were determined with a zNOSETM (Electronic Sensor Technology, Newbury Park, CA) every 2 h (see below), from 6 pm to 4 am. These experiments were performed with 6 biological replicates per genotype.
Generation of transgenic ir-hd20 lines
Transgenic N. attenuata (ir-hd20) plants reduced in the expression of HD20 were generated via elongated hypocotyl Agrobacterium-mediated transformation and seedling regeneration as previously described . The binary vector used for plant transformation was pSOL8  engineered to carry a 329 bp fragment of the HD20 mRNA corresponding to nucleotides 529 to 857 (GenBank accession: HM107874) and subcloned in inverted repeat orientation to generate pSOL8HD20. T1 transformed plants were analyzed by quantification of ABA levels (see below) and HD20 mRNA levels (see below) in leaves during water stress(i.e. by withholding irrigation for 5 days) and for T-DNA single insertion by Southern blot hybridization (see below). Segregation analysis of hygromycin resistance in T2 seedlings was performed on agar plates supplemented with hygromycin (0.025 mg mL-1). Two lines, ir-hd20 A-09-408 (ir-hd20 line 1) and A-09-411 (ir-hd20 line 2) had the lowest levels of ABA and HD20 mRNA accumulation at 5 days after withholding water and had a single T-DNA insertion in their genomes. These lines were used for all experiments. Southern blot analysis was performed as previously described .
Analysis of phytohormones
For analysis of JA, JA-Ile, SA and ABA, 0.2 g of frozen tissue were homogenized to a fine powder in the presence of liquid nitrogen. One mL of ethylacetate spiked with 200 ng [²H2JA and 40 ng 13 C6JA-Ile, 2 H6ABA, 2 H4SA was added to the samples and after vortexing the samples were centrifuged for 15 min at 13,200 g (4°C). The upper organic phase was transferred into a fresh tube and the leaf material was re-extracted with 0.5 mL ethylacetate. The organic phases were pooled and evaporated to dryness. The dry residue was reconstituted in 0.4 mL of 70/30 (v/v) methanol/water for analysis with an LC-ESI-MS/MS instrument (Varian 1200 Triple-Quadrupole-LC-MS system; Varian, Palo Alto, CA). Ten μL of the sample were injected in a ProntoSIL® column (C18-ace-EPS, 50 x 2 mm, 5 μm, 120 Å, Bischoff, Leonberg, Germany) connected to a pre-column (C18, 4 x 2 mm, Phenomenex, Torrance, CA). As mobile phases 0.05% formic acid in water (solvent A) and methanol (solvent B) were used in a gradient mode with the following conditions: time/concentration (min/%) for B: 0.0/15; 2.5/15; 4.5/98; 10.5/98; 12.0/15; 15.0/15; time/flow (min/mL): 0.0/0.4; 1.5/0.2; 1.5/0.2; 10.5/0.4; 15.0/0.4. Compounds were detected in the ESI negative mode and multiple reaction monitoring (MRM) according to published parameters .
Ethylene (ET) measurements
ET levels were analyzed by photoacoustic spectrometry (INVIVO; https://www.invivo-gmbh.de) as previously described . Corollas were excised from antherectomized flowers at different times (Additional file 3) and immediately placed in 100 mL sealed-glass vessels (four corollas per vessel) for 4.5 h before ET quantification. Three biological replicates per time and per genotype were performed. Antherectomization was performed in the morning before pollen release to avoid the production of ET induced by fertilization.
Quantification of benzyl acetone (BA) emission from flowers
Flower volatiles were collected by enclosing individual night flowers in plastic cups connected to Super-Q filter traps (ARS, Philadelphia, PA) under an air flow of 30 mL min-1. Volatiles were collected from 3 pm to 9 am. Filter traps were spiked with 400 ng of tetraline as an internal standard (IS) and eluted with 250 μL of dichloromethane. Volatile analysis was performed by GC-MS with a CP-4000 GC instrument (Varian, Palo Alto, CA) on a DB-5 column (Agilent, Waldbronn, Germany). BA was identified by comparing retention times and mass spectra with a commercial BA standard (Sigma, Taufkirchen, Germany). Peak areas were integrated and normalized to the area of the IS.
BA was also analyzed by the zNOSETM as previously described . In this case, individual flowers were placed in scintillation vials and the vials were left capped. Every 2 h (starting at 6 pm; Figure 4d) the needle of the zNOSETM was introduced into the vials to determine BA levels in the headspace. The vials were immediately capped again. The intensity of the peak area corresponding to BA was used to quantify BA levels per flower.
Analysis of gene expression by qPCR
Total RNA was extracted from different tissues by the TRIzol® reagent (Invitrogen) and DNase-I treated (Fermentas, St. Leon-Rot, Germany) according to commercial instructions. Five μg of total RNA were reverse transcribed using oligo(dT)18 and SuperScript reverse transcriptase II (Invitrogen) according to commercial instructions. Quantitative real-time PCR (qPCR) was performed with the Mx3005P Multiplex qPCR system (Stratagene, La Jolla, CA) and the qPCR Core kit for SYBR® Green I (Eurogentec, Liege, Belgium) using gene specific primers (Additional file 6). Quantification of HD20 mRNA levels was performed by normalization with the Elongation Factor 1A (EF1A) mRNA according to the ΔCt method . All the reactions were performed with at least three biological replicates.
Corollas from ir-hd20 and WT plants were collected at 12 pm (a total of three independent samples (biological replicates) were used per genotype and each sample was composed of three pooled corollas). Total RNA was extracted from corollas based on a previously described method  and its quality was checked by spectrophotometry (NanoDrop, Wilmington, DE). Genomic DNA was removed by DNAse treatment following commercial instructions (Turbo DNase; Ambion, Europe), RNA was cleaned up with RNeasy MinElute columns (Qiagen, Hilden, Germany) and the RNA quality was checked with the RNA 6000 n kit (Agilent, Santa Clara, CA) using an Agilent 2100 Bioanalyzer. Total RNA was used to generate labeled cRNA with the Quick Amp labeling kit (Agilent) and the yield was determined spectrophotometrically (NanoDrop). Labeled cRNA was hybridized using the gene expression hybridization kit (Agilent) following commercial instructions onto 44 K custom-designed 60-mer N. attenuata Agilent microarray containing 43,533 sequences (see accession numbers)[28, 29]. Microarrays were hybridized overnight at 65°C and slides were washed with the Gene Expression Wash Buffer kit (Agilent) as outlined in the One-Color Microarray-Based Gene Expression Analysis manual (Agilent). Three biological replicates were used per treatment with a total of six arrays. Arrays were scanned with an Agilent G2565BA scanner and image data was acquired with the Agilent Scan Control software (version A.7.0.1 for the B scanner). Data was extracted using the Agilent Feature Extraction software (version 9.5) and analyzed with the SAM (Significance Analysis of Microarrays) software . The q-values for each gene corresponded to a computed false discovery rate (FDR) of 5.2%. Significant changes in gene expression were considered when the fold change (FC; ir-hd20 versus WT) were equal or greater than 1.9 or equal or smaller than 0.6 (with q-values lower than 0.048; according to the FDR value calculated by SAM). Gene Ontology (GO) and enzyme code (EC) classification and sequence analysis was performed with the Blast2GO software .
Statistics were calculated using the SPSS software version 17.0. The data was subjected either to one-way analysis of variance (ANOVA; and means were compared by the lowest standard deviation (LSD) test) or to Student’s t-test as indicated in the text. The number of replicates (n) used in each experiment are detailed in the figure’s captions.
Sequence data from this article can be found under the following accession numbers: HD20 [GenBank:HM107874], CHAL1 [GenBank:EU503226], FDR [GenBank:JQ028693], PP2C [GenBank:JQ028694], SAMT/BAMT [GenBank: JQ028692], Agilent Chip platform [NCBI GEO:GPL13527], and microarray data [NCBI GEO: GSE33682].
The DAAD (fellowship to DR), the ERASMUS Program (fellowship to BR) and the Max Planck Society are acknowledged for funding.
- Dudareva N, Pichersky E: Biochemical and molecular genetic aspects of floral scents. Plant Physiol. 2000, 122 (3): 627-633. 10.1104/pp.122.3.627.PubMedPubMed CentralView ArticleGoogle Scholar
- Kessler D, Gase K, Baldwin IT: Field experiments with transformed plants reveal the sense of floral scents. Science. 2008, 321 (5893): 1200-1202. 10.1126/science.1160072.PubMedView ArticleGoogle Scholar
- Kessler D, Baldwin IT: Making sense of nectar scents: the effects of nectar secondary metabolites on floral visitors of Nicotiana attenuata. Plant J. 2007, 49 (5): 840-854. 10.1111/j.1365-313X.2006.02995.x.PubMedView ArticleGoogle Scholar
- Raguso RA: Flowers as sensory billboards: progress towards an integrated understanding of floral advertisement. Curr Opin Plant Biol. 2004, 7 (4): 434-440. 10.1016/j.pbi.2004.05.010.PubMedView ArticleGoogle Scholar
- Spitzer-Rimon B, Marhevka E, Barkai O, Marton I, Edelbaum O, Masci T, Prathapani NK, Shklarman E, Ovadis M, Vainstein A: EOBII, a gene encoding a flower-specific regulator of phenylpropanoid volatiles' biosynthesis in petunia. Plant Cell. 2010, 22 (6): 1961-1976. 10.1105/tpc.109.067280.PubMedPubMed CentralView ArticleGoogle Scholar
- Verdonk JC, Haring MA, van Tunen AJ, Schuurink RC: ODORANT1 regulates fragrance biosynthesis in petunia flowers. Plant Cell. 2005, 17 (5): 1612-1624. 10.1105/tpc.104.028837.PubMedPubMed CentralView ArticleGoogle Scholar
- Pichersky E, Dudareva N: Scent engineering: toward the goal of controlling how flowers smell. Trends Biotechnol. 2007, 25 (3): 105-110. 10.1016/j.tibtech.2007.01.002.PubMedView ArticleGoogle Scholar
- Dudareva N, Martin D, Kish CM, Kolosova N, Gorenstein N, Faldt J, Miller B, Bohlmann J: (E)-beta-ocimene and myrcene synthase genes of floral scent biosynthesis in snapdragon: function and expression of three terpene synthase genes of a new terpene synthase subfamily. Plant Cell. 2003, 15 (5): 1227-1241. 10.1105/tpc.011015.PubMedPubMed CentralView ArticleGoogle Scholar
- Martin C, Gerats T: Control of Pigment Biosynthesis Genes during Petal Development. Plant Cell. 1993, 5 (10): 1253-1264.PubMedPubMed CentralView ArticleGoogle Scholar
- Euler M, Baldwin IT: The chemistry of defense and apparency in the corollas of Nicotiana attenuata. Oecologia. 1996, 107 (1): 102-112. 10.1007/BF00582240.View ArticleGoogle Scholar
- Sime KR, Baldwin IT: Opportunistic out-crossing in Nicotiana attenuata (Solanaceae), a predominantly self-fertilizing native tobacco. BMC Ecol. 2003, 3: 6-10.1186/1472-6785-3-6.PubMedPubMed CentralView ArticleGoogle Scholar
- Baldwin IT, Preston C, Euler M, Gorham D: Patterns and consequences of benzyl acetone floral emissions from Nicotiana attenuata plants. . 1997, 23 (10): 2327-2343.Google Scholar
- Ariel FD, Manavella PA, Dezar CA, Chan RL: The true story of the HD-Zip family. Trends Plant Sci. 2007, 12 (9): 419-426. 10.1016/j.tplants.2007.08.003.PubMedView ArticleGoogle Scholar
- Sessa G, Morelli G, Ruberti I: DNA-binding specificity of the homeodomain-leucine zipper domain. J Mol Biol. 1997, 274 (3): 303-309. 10.1006/jmbi.1997.1408.PubMedView ArticleGoogle Scholar
- Palena CM, Gonzalez DH, Chan RL: A monomer-dimer equilibrium modulates the interaction of the sunflower homeodomain leucine-zipper protein Hahb-4 with DNA. Biochem J. 1999, 341 (Pt 1): 81-87.PubMedPubMed CentralGoogle Scholar
- Carabelli M, Sessa G, Baima S, Morelli G, Ruberti I: The Arabidopsis Athb-2 and −4 genes are strongly induced by far-red-rich light. Plant J. 1993, 4 (3): 469-479. 10.1046/j.1365-313X.1993.04030469.x.PubMedView ArticleGoogle Scholar
- Soderman E, Mattsson J, Engstrom P: The Arabidopsis homeobox gene ATHB-7 is induced by water deficit and by abscisic acid. Plant J. 1996, 10 (2): 375-381. 10.1046/j.1365-313X.1996.10020375.x.PubMedView ArticleGoogle Scholar
- Soderman E, Hjellstrom M, Fahleson J, Engstrom P: The HD-Zip gene ATHB6 in Arabidopsis is expressed in developing leaves, roots and carpels and up-regulated by water deficit conditions. Plant Mol Biol. 1999, 40 (6): 1073-1083. 10.1023/A:1006267013170.PubMedView ArticleGoogle Scholar
- Henriksson E, Olsson AS, Johannesson H, Johansson H, Hanson J, Engstrom P, Soderman E: Homeodomain leucine zipper class I genes in Arabidopsis. Expression patterns and phylogenetic relationships. Plant Physiol. 2005, 139 (1): 509-518. 10.1104/pp.105.063461.PubMedPubMed CentralView ArticleGoogle Scholar
- Olsson AS, Engstrom P, Soderman E: The homeobox genes ATHB12 and ATHB7 encode potential regulators of growth in response to water deficit in Arabidopsis. Plant Mol Biol. 2004, 55 (5): 663-677. 10.1007/s11103-004-1581-4.PubMedView ArticleGoogle Scholar
- Lee YH, Chun JY: A new homeodomain-leucine zipper gene from Arabidopsis thaliana induced by water stress and abscisic acid treatment. Plant Mol Biol. 1998, 37 (2): 377-384. 10.1023/A:1006084305012.PubMedView ArticleGoogle Scholar
- Johannesson H, Wang Y, Hanson J, Engstrom P: The Arabidopsis thaliana homeobox gene ATHB5 is a potential regulator of abscisic acid responsiveness in developing seedlings. Plant Mol Biol. 2003, 51 (5): 719-729. 10.1023/A:1022567625228.PubMedView ArticleGoogle Scholar
- Mayda E, Tornero P, Conejero V, Vera P: A tomato homeobox gene (HD-zip) is involved in limiting the spread of programmed cell death. Plant J. 1999, 20 (5): 591-600. 10.1046/j.1365-313X.1999.00633.x.PubMedView ArticleGoogle Scholar
- Lin Z, Hong Y, Yin M, Li C, Zhang K, Grierson D: A tomato HD-Zip homeobox protein, LeHB-1, plays an important role in floral organogenesis and ripening. Plant J. 2008, 55 (2): 301-310. 10.1111/j.1365-313X.2008.03505.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Deng X, Phillips J, Brautigam A, Engstrom P, Johannesson H, Ouwerkerk PB, Ruberti I, Salinas J, Vera P, Iannacone R, et al: A homeodomain leucine zipper gene from Craterostigma plantagineum regulates abscisic acid responsive gene expression and physiological responses. Plant Mol Biol. 2006, 61 (3): 469-489. 10.1007/s11103-006-0023-x.PubMedView ArticleGoogle Scholar
- Harris JC, Hrmova M, Lopato S, Langridge P: Modulation of plant growth by HD-Zip class I and II transcription factors in response to environmental stimuli. New Phytol. 2011, 190 (4): 823-837. 10.1111/j.1469-8137.2011.03733.x.PubMedView ArticleGoogle Scholar
- Re DA, Dezar CA, Chan RL, Baldwin IT, Bonaventure G: Nicotiana attenuata NaHD20 plays a role in leaf ABA accumulation during water stress, benzylacetone emission from flowers, and the timing of bolting and flower transitions. . 2011, 62 (1): 155-166. 101093/jxb/erm035Google Scholar
- Kallenbach M, Gilardoni PA, Allmann S, Baldwin IT, Bonaventure G: C(12) derivatives of the hydroperoxide lyase pathway are produced by product recycling through lipoxygenase-2 in Nicotiana attenuata leaves. New Phytol. 2011, 191 (4): 1054-1068. 10.1111/j.1469-8137.2011.03767.x.PubMedView ArticleGoogle Scholar
- Gilardoni PA, Hettenhausen C, Baldwin IT, Bonaventure G: Nicotiana attenuata LECTIN RECEPTOR KINASE1 suppresses the insect-mediated inhibition of induced defense responses during Manduca sexta herbivory. Plant Cell. 2011, 23 (9): 3512-3532. 10.1105/tpc.111.088229.PubMedPubMed CentralView ArticleGoogle Scholar
- Devoto A, Nieto-Rostro M, Xie DX, Ellis C, Harmston R, Patrick E, Davis J, Sherratt L, Coleman M, Turner JG: COI1 links jasmonate signalling and fertility to the SCF ubiquitin-ligase complex in Arabidopsis. Plant J. 2002, 32 (4): 457-466. 10.1046/j.1365-313X.2002.01432.x.PubMedView ArticleGoogle Scholar
- Paschold A, Halitschke R, Baldwin I: Co(i)-ordinating defenses: NaCOI1 mediates herbivore- induced resistance in Nicotiana attenuata and reveals the role of herbivore movement in avoiding defenses. Plant J. 2007, 51 (1): 79-91. 10.1111/j.1365-313X.2007.03119.x.PubMedView ArticleGoogle Scholar
- Li L, Zhao YF, McCaig BC, Wingerd BA, Wang JH, Whalon ME, Pichersky E, Howe GA: The tomato homolog of CORONATINE-INSENSITIVE1 is required for the maternal control of seed maturation, jasmonate-signaled defense responses, and glandular trichome development. Plant Cell. 2004, 16 (1): 126-143. 10.1105/tpc.017954.PubMedPubMed CentralView ArticleGoogle Scholar
- Arce AL, Raineri J, Capella M, Cabello JV, Chan RL: Uncharacterized conserved motifs outside the HD-Zip domain in HD-Zip subfamily I transcription factors; a potential source of functional diversity. BMC Plant Biol. 2011, 11: 42-10.1186/1471-2229-11-42.PubMedPubMed CentralView ArticleGoogle Scholar
- Boudsocq M, Barbier-Brygoo H, Lauriere C: Identification of nine sucrose nonfermenting 1-related protein kinases 2 activated by hyperosmotic and saline stresses in Arabidopsis thaliana. . 2004, 279 (40): 41758-41766. 101074/jbcM705780200Google Scholar
- Raghavendra AS, Gonugunta VK, Christmann A, Grill E: ABA perception and signalling. Trends Plant Sci. 2010, 15 (7): 395-401. 10.1016/j.tplants.2010.04.006.PubMedView ArticleGoogle Scholar
- Winkel-Shirley B: Flavonoid biosynthesis. A colorful model for genetics, biochemistry, cell biology, and biotechnology. Plant Physiol. 2001, 126 (2): 485-493. 10.1104/pp.126.2.485.PubMedPubMed CentralView ArticleGoogle Scholar
- Colquhoun TA, Kim JY, Wedde AE, Levin LA, Schmitt KC, Schuurink RC, Clark DG: PhMYB4 fine-tunes the floral volatile signature of Petunia x hybrida through PhC4H. . 2011, 62 (3): 1133-1143. 101093/jxb/erm035Google Scholar
- Cabello JV, Arce AL, Chan RL: The homologous HD-Zip I transcription factors HaHB1 and AtHB13 confer cold tolerance via the induction of pathogenesis-related and glucanase proteins. Plant J. 2011, 69 (1): 141-153.PubMedView ArticleGoogle Scholar
- Park CJ, Peng Y, Chen X, Dardick C, Ruan D, Bart R, Canlas PE, Ronald PC: Rice XB15, a protein phosphatase 2 C, negatively regulates cell death and XA21-mediated innate immunity. PLoS Biology. 2008, 6 (9): e231-10.1371/journal.pbio.0060231.PubMedPubMed CentralView ArticleGoogle Scholar
- Klingler JP, Batelli G, Zhu JK: ABA receptors: the START of a new paradigm in phytohormone signalling. . 2010, 61 (12): 3199-3210. 101093/jxb/erm035Google Scholar
- Hirayama T, Shinozaki K: Perception and transduction of abscisic acid signals: keys to the function of the versatile plant hormone ABA. Trends Plant Sci. 2007, 12 (8): 343-351. 10.1016/j.tplants.2007.06.013.PubMedView ArticleGoogle Scholar
- Finkelstein RR, Somerville CR: Three Classes of Abscisic Acid (ABA)-Insensitive Mutations of Arabidopsis Define Genes that Control Overlapping Subsets of ABA Responses. Plant Physiol. 1990, 94 (3): 1172-1179. 10.1104/pp.94.3.1172.PubMedPubMed CentralView ArticleGoogle Scholar
- Nelson DC, Lasswell J, Rogg LE, Cohen MA, Bartel B: FKF1, a clock-controlled gene that regulates the transition to flowering in Arabidopsis. Cell. 2000, 101 (3): 331-340. 10.1016/S0092-8674(00)80842-9.PubMedView ArticleGoogle Scholar
- Xie DX, Feys BF, James S, Nieto-Rostro M, Turner JG: COI1: An Arabidopsis gene required for jasmonate-regulated defense and fertility. Science. 1998, 280 (5366): 1091-1094. 10.1126/science.280.5366.1091.PubMedView ArticleGoogle Scholar
- von Dahl CC, Winz RA, Halitschke R, Kuhnemann F, Gase K, Baldwin IT: Tuning the herbivore-induced ethylene burst: the role of transcript accumulation and ethylene perception in Nicotiana attenuata. Plant J. 2007, 51 (2): 293-307. 10.1111/j.1365-313X.2007.03142.x.PubMedView ArticleGoogle Scholar
- von Dahl CC, Baldwin IT: Deciphering the role of ethylene in plant-herbivore interactions. Journal of Plant Growth Regulation. 2007, 26 (2): 201-209. 10.1007/s00344-007-0014-4.View ArticleGoogle Scholar
- Krügel T, Lim M, Gase K, Halitschke R, Baldwin IT: Agrobacterium-mediated transformation of Nicotiana attenuata, a model ecological expression system. Chemoecology. 2002, 12 (4): 177-183. 10.1007/PL00012666.View ArticleGoogle Scholar
- Gase K, Weinhold A, Bozorov T, Schuck S, Baldwin IT: Efficient screening of transgenic plant lines for ecological research. Molecular Ecology Resources. 2011, 11 (5): 890-902. 10.1111/j.1755-0998.2011.03017.x.PubMedView ArticleGoogle Scholar
- Bonaventure G, Schuck S, Baldwin IT: Revealing complexity and specificity in the activation of lipase-mediated oxylipin biosynthesis: a specific role of the Nicotiana attenuata GLA1 lipase in the activation of jasmonic acid biosynthesis in leaves and roots. Plant Cell Environ. 2011, 34 (9): 1507-1520. 10.1111/j.1365-3040.2011.02348.x.PubMedView ArticleGoogle Scholar
- Körner E, von Dahl C, Bonaventure G, Baldwin IT: Pectin methylesterase NaPME1 contributes to the emission of methanol during insect herbivory and to the elicitation of defence responses in Nicotiana. . 2009, 60: 2631-2640. 101093/jxb/erm035Google Scholar
- Kessler D, Diezel C, Baldwin IT: Changing pollinators as a means of escaping herbivores. Curr Biol. 2010, 20 (3): 237-242. 10.1016/j.cub.2009.11.071.PubMedView ArticleGoogle Scholar
- Gilardoni P, Schuck S, Jüngling R, Rotter B, Baldwin I, Bonaventure G: SuperSAGE analysis of the Nicotiana attenuata transcriptome after fatty acid-amino acid elicitation (FAC): identification of early mediators of insect responses. BMC Plant Biol. 2010, 10: 66-10.1186/1471-2229-10-66.PubMedPubMed CentralView ArticleGoogle Scholar
- Kistner C, Matamoros M: RNA isolation using phase extraction and LiClprecipitation. In Lotus japonicus Handbook. Edited by Márquez A.Netherlands: Springer; 2005:123–124.View ArticleGoogle Scholar
- Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498.PubMedPubMed CentralView ArticleGoogle Scholar
- Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.PubMedView ArticleGoogle Scholar