This Article Abstract Full Text (PDF) Supplemental Data All Versions of this Article: 135/1/539    most recent pp.103.037036v1 Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Similar articles in Web of Science Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via CrossRef Citing Articles via Web of Science (12) Citing Articles via Google Scholar Google Scholar Articles by Qu, N. Articles by Baldwin, I. T. Search for Related Content PubMed PubMed Citation Articles by Qu, N. Articles by Baldwin, I. T. Agricola Articles by Qu, N. Articles by Baldwin, I. T.

PLANTS INTERACTING WITH OTHER ORGANISMS

Consistency of Nicotiana attenuata's Herbivore- and Jasmonate-Induced Transcriptional Responses in the Allotetraploid Species Nicotiana quadrivalvis and Nicotiana clevelandii^1,[w]

Department of Molecular Ecology, Max-Planck-Institute for Chemical Ecology, Beutenberg Campus, 07745 Jena, Germany

	ABSTRACT

TOP ABSTRACT RESULTS DISCUSSION MATERIALS AND METHODS LITERATURE CITED

 
We examined the consistency of the native diploid Nicotiana attenuata (Na)'s herbivore-induced transcriptional changes in the two allotetraploid natives, Nicotiana clevelandii (Nc) and Nicotiana quadrivalvis (Nq), which are thought to be derived from hybridizations with an ancestral Na. An analysis of nuclear-encoded chloroplast-expressed Gln synthetase gene (ncpGS) sequences found strong similarity between Nc and Na and between N. trigonophylla and the two allopolyploids. All species were elicited with methyl jasmonate (MeJA), or were wounded and treated with either water, Manduca sexta oral secretions and regurgitant (R), or the two most abundant fatty acid amino acid conjugates (F) in R to simulate herbivory. The induced transcriptional responses in all three species were compared with a cDNA microarray enriched in Na genes. Na had the fastest transcriptional responses followed by Nc and then Nq. Na's R- and F-elicited responses were more similar to those from Nq, while the MeJA- or wound-elicited responses were more consistent in Nc. Treatment of wounds with the full cocktail of elicitors found in R elicits more complex responses than does treatment with F. The species differ in their elicited JA responses, and these differences are mirrored in the expression of oxylipin genes (LOX, HPL, AOS, and -DOX) and downstream JA-elicited genes (TD). Elicitation decreases the expression of growth-related genes in all three species. We propose that this is a valuable system to examine the modification of complex, polygenic, adaptive responses during allopolyploid speciation.

The relationship between adaptation and polyploidy is complex. Some polyploids are more resistant to herbivore and pathogen attack than are their diploid relatives (Schoen et al., 1992), but the effects of polyploidy on herbivore resistance are not always observed in a single lineage (Thompson et al., 1997; Nuismer and Thompson, 2001; Janz and Thompson, 2002). For example, a 4-year study of diploid and tetraploid individuals of the host plant Heuchera grossulariifolia, with its three most abundant species of phytophagous insects (Greya politella, Greya piperella, and Eupithecia misturata), did not reveal consistent associations between ploidy levels and insect resistance (Nuismer and Thompson, 2001). Allopolyploidy, the combination of two or more distinct nuclear genomes, is a common speciation mechanism in many taxa. In these taxa, the merger paradox is particularly acute, as the species-specific adaptations of the ancestor species are likely disrupted by the genomic reorganization that accompanies this form of instantaneous speciation. The genus Nicotiana contains many examples of allopolyploid speciation as well as one of the best-documented examples of a complex polygenic adaptation: the induced responses to herbivore attack.

Cytologically, Nicotiana attenuata (Na; synonymous with Nicotiana torreyana Nelson and Macbr.) is a 12-paired diploid species and is thought to be the common ancestor of two other North American tetraploid species, Nicotiana quadrivalvis (Nq; synonymous with Nicotiana Bigelovii; Chase et al., 2003) and Nicotiana clevelandii (Nc), according to their similarity in caryotypes and leaf, inflorescence, and trichome morphology (Goodspeed, 1954). Goodspeed also hypothesized that the other ancestor of the two tetraploids was an alatoid species. In agreement with this hypothesis, recently published genomic in situ hybridization (GISH) experiments used DNA probes from Na as well as Nicotiana sylvestris hybridized to the genomes of Nc and Nq, but the specificity of the probes is unclear (Chase et al., 2003). Phylogenies of the genus Nicotiana based on the internal transcribed spacer (ITS) regions and maturase K (matK) sequences (an analysis which did not include Nc) suggest that an ancestor of the species in N. section Trigonphyllae might have contributed to the genomes of the two tetraploids (Komarnyts'kyi et al., 1998; Aoki and Ito, 2000; Chase et al., 2003). Chloroplast-expressed Gln synthetase is a nuclear-encoded gene (ncpGS) that occurs as a single copy in most diploid taxa examined and is useful for phylogenetic studies (Emshwiller and Doyle, 1999). To examine Goodspeed's 1954 hypothesis, we cloned and compared the ncpGS chromosomal sequences from Na, Nq, and Nc, as well as from eight other wild species from the subgenus Petunioides.

Plant-herbivore interactions have been intensively investigated using ecological, chemical, and molecular approaches in the Na-Manduca sexta system (Baldwin and Preston, 1999). Extensive analysis reveals that Na can recognize attack from the larvae of its specialist Sphingid herbivore, M. sexta, as evidenced by an endogenous jasmonic acid (JA) burst (Baldwin et al., 1994). Following the JA burst, proteinase inhibitor (PI) accumulation and the release of volatile organic compounds (VOCs) increased dramatically (Halitschke et al., 2001). PIs act as direct defenses, slowing the growth of the larvae and keeping them at a size which is vulnerable to generalist predators, while VOCs, as an indirect defense, attract generalist predators to the feeding larvae (Pohnert et al., 1999; Kessler and Baldwin, 2001, 2002; Glawe et al., 2003). In addition, a large-scale transcriptional reorganization (Schittko et al., 2000, 2001; Hermsmeier et al., 2001; Voelckel et al., 2001; Halitschke et al., 2003; Hui et al., 2003; Voelckel and Baldwin, 2003) accompanies the up-regulation of direct and indirect defenses elicited by M. sexta attack. The available evidence suggests that many of the responses are adaptive, increasing the fitness of plants attacked by M. sexta larvae. Moreover, because the responses are clearly polygenic, they provide an ideal system in which to study how complex polygenic adaptations are conserved or modified during allopolyploid speciation. To conduct such a comparative analysis, the responses must be elicited in a consistent and reproducible manner.

When herbivores attack plants, they cause wounding and introduce herbivore-specific elicitors into the wounds, which results in complex transcriptional and metabolic changes that are different from those caused by mechanical wounding or simple application of JA. Several types of elicitors have been found in the oral secretions and regurgitant (R) of the herbivores, such as -glycosidase (Mattiacci et al., 1995) and Glc oxidase (Musser et al., 2002) as well as fatty acid-amino acid conjugates (FACs), which can alter plants' responses to wounds (Halitschke et al., 2001) and mimic herbivore-induced changes. Eight FACs in the R from M. sexta account for most of the herbivore-specific changes in Na defense-responses, including JA and PI bursts and the release of VOCs (Halitschke et al., 2001). Moreover, the treatment of wounds with only two (of the eight) FACs found in Manduca R at the R-equivalent concentrations elicits more than 55% of the R-specific transcriptional changes. Hence, Manduca-specific changes can be elicited by applying R to the wounds, and parts of the responses can be elicited by applying water or FACs to wounds, or by treating plants with the methyl ester of JA (MeJA).

M. sexta larvae, one of the three most abundant and damaging herbivores for Na, has also been found on the two allotetraploid descendants, Nq and Nc, in their native habitats (I.T. Baldwin, unpublished data). Comparisons of M. sexta recognition (JA burst) and resistance responses (PI, nicotine, VOCs, and other secondary metabolites) revealed that most of the M. sexta-induced responses are retained with modifications in Nq but lost in Nc (Lou and Baldwin, 2003). Here we provide a transcriptional analysis of the Manduca-specific responses in Na, Nc, and Nq using a cDNA-microarray containing 246 Manduca-R-induced Na genes.

To examine the Manduca-specific responses in the three species, the source-sink transition leaves from a plant at the rosette stage of growth were either (1) wounded and immediately treated with 20 µL water, 1:10 diluted R, or the two most abundant FACs (F) in the R; or (2) directly treated with MeJA in lanolin paste. The transcriptional responses were examined 24 h after elicitation in all three species. Since transcriptional elicitation by M. sexta R is a rapidly induced and waning response in Na (Halitschke et al., 2003), while R-elicited metabolic responses are slower in Nq and Nc than in Na (Lou and Baldwin, 2003), the M. sexta R- and F-specific transcriptional responses from all three species were also compared 2 h after elicitation. Since the induced responses occur in both local and systemic tissues, albeit with a delayed kinetic, we examined the transcriptional responses in both directly treated and systemic leaves 24 h after elicitation.

	RESULTS

TOP ABSTRACT RESULTS DISCUSSION MATERIALS AND METHODS LITERATURE CITED

 

Molecular Support for Goodspeed's Phylogeny

With the primer pairs GS57, GS2, NQ40, and NQ41, we cloned and sequenced ncpGS genes from Na, Nq, Nc, and 8 other species closely related to Na. The sequence analysis of 8 clones from Nc contained two ncpGS isoforms, but the 12 clones from Nq all proved to be of the same isoform. Consensus lengths of approximately 650 nucleotides of ncpGS genes covering one exon and intron completely and another exon and intron partially are used for multiple alignments with ClustalW. To root the phylogenetic tree, the cDNA sequence from Solanum tuberosum was also included in which the introns were counted as indels. The distances were calculated according to the Jukes-Cantor one parameter evolution model. The phylogram obtained from the neighbor-joining algorithm was drawn by TreeCon. Two clades are recognizable in the neighbor-joining tree of Nicotiana: (1) one consisting of N. trigonophylla and two species from the N. section of Bigelovianae, Nq and Nc. In this clade, Nq and Nc are placed in the same subclade as sister taxa; and (2) one consisting of N. sylvestris, Nc, and all the examined species of N. section Acuminatae, which include Na. This clade is split again into two subclades, one consisting of only two species from the Acuminatae section, N. linearis and N. spegazzinii, and the rest of the Acuminatae section, with the addition of Nc. One ncpGS gene from Nc grouped closely with Na and had a bootstrap value of 100% (Fig. 1). This phylogeny is consistent with Goodspeed's hypothesis that all three species (Na, Nq, and Nc) share a recent common ancestor and suggests that a progenitor of N. trigonophylla is involved in the formation of Nq and Nc.

View larger version (13K): [in this window] [in a new window]   Figure 1. Neighbor-joining tree of the ncpGS sequences from selected North American Nicotiana species to test Goodspeed's hypothesis for the relationships among the diploid N. attenuata and the two allopolyploids, N. clevelandii and N. quadrivalvis. The chloroplast Gln synthetase gene from S. tuberosum was included to root the phylogenic tree. The branch lengths (bar) and bootstrap (bold, as percentages of 100 replicates) and the chromosome number (n value) of each species (from Goodspeed, 1954) are given. GenBank accession numbers: AF302113 (S. tuberosum); AY566989 (N. sylvestris); AY568077 (N. corymbosa); AY568079 (N. pauciflora); AY568078 (N. acuminate); AY568082 (N. linearis); AY568081 (N. spegazzinii); AY568084 (N. quadrivalvis); AY568083 (N. trigonophylla); AY568080 (N. miersii); AY566990 (N. clevelandii 1); AY566991 (N. clevelandii 2); and AY183657 (N. attenuata).

The 24-h transcriptional responses elicited by either MeJA treatment or wounds (W) treated with water (W + W), R (W + R) or F (W + F) are examined separately in the treated leaves and the untreated systemic leaves from all three species (Table I; Supplemental Table I, available at www.plantphysiol.org). The analysis of the 24 hybridized arrays from the 24 h treatments revealed that 217 (88% of 246) genes were differentially expressed with 91.2% (198), 65% (141), and 55.3% (120) being differentially regulated in Na, Nq, and Nc, respectively. Of these 217 genes, 150 genes were regulated by the different treatments in both diploid and allotetraploid species, and 67 genes were elicited only in the allotetraploid (19) or diploid species (48). The MeJA treatment allowed the 150-µg dose of MeJA to diffuse from the lanolin paste into the leaf over a protracted period of time (Zhang and Baldwin, 1997); hence we expected this treatment to produce the largest transcriptional response at the 24 h harvest. MeJA elicitation differentially regulated 133, 81, and 70 genes in the treated leaves and 33, 63, and 24 genes in the untreated systemic leaves of Na, Nq, and Nc, respectively. Of the MeJA-elicited responses, 51.1% and 45.1% (local) or 36% and 30% (systemic) from Na were also observed in Nq and Nc. Most of those elicited in the treated leaves of allotetraploid (84%, Nq and 86%, Nc) were also elicited in the treated leaves of Na. However, fewer of the MeJA-elicited genes from allotetraploid species (19%, Nq and 41%, Nc) in systemic leaves were elicited in Na (Table I).

To provide a global analysis of the Na-specific induced transcriptional responses in Nq and Nc, we log-10 transformed the expression ratios and conducted a principal component analysis (PCA). A two-dimensional presentation of the PCA accounted for about one-half of the variation in the data; axes 1 and 2 explain 55.1% and 59.2% of the total variance observed in treated (Fig. 2A) and systemic (Fig. 2B) leaves, respectively. This analysis revealed that the responses in treated leaves are more distinct than those observed in systemic leaves. The local responses of the three species were largely resolved along axis 2 (accounting for 17.7% of the variance; Fig. 2A), while no clear species-specific pattern was observed in the systemic (and younger) leaves. MeJA elicitation, as expected, produced transcriptional imprints (transient changes in the transcriptome) that are distinct in both local and systemic leaves, but again the differences are clearer in the elicited leaves. The transcriptional imprints resulting from the treatment of wounds with R or F modifies the wound-induced response in all three species, but the differences are more apparent in the systemic tissues.

View larger version (17K): [in this window] [in a new window]   Figure 2. A PCA of the logarithmic expression ratios from 24 cDNA microarrays hybridized with fluorescently labeled cDNA isolated from Na, Nq, and Nc. At the rosette stage of growth, node 0 leaves from each species were wounded and immediately treated with: (1) 20 µL water, 1:10 diluted R, or F in R; or (2) directly treated with MeJA in lanolin paste. The treated leaf (A- local) at node 0 as well as the next younger leaf at position –1 (B- systemic) were harvested for mRNA extraction 24 h after treatment. Identical leaves were also harvested from noninduced control plants. Microarrays were hybridized with cDNA reverse-transcribed from mRNA pooled from 20 independently treated plants and Cy3-labeled (elicited plants) or Cy5-labeled (control plants). The MeJA treatment does not include wound responses, while water, R, and F treatments include wound responses. Together, axes 1 and 2 of the PCA account for 55.1% (local) and 59.2% (systemic) the variance. White circles, squares, and black squares represent Na, Nq, and Nc, respectively.

In summary, MeJA elicitation leads to more distinctive transcriptional changes compared to R or F elicitation in all three species. MeJA-elicited local responses from two allotetraploids are more similar to those from Na than are MeJA-elicited systemic responses. In comparison to F-elicited responses, the R-elicited responses are less similar among three species. Local recognition of R and F elicitation was faster in comparison to the systemic responses in all three species and, furthermore, Na showed the fastest R- and F-elicited responses (followed by Nc and then Nq).

Two-Hour Transcriptional Responses

From the analysis conducted at 24 h, it is clear that by then elicited recognition responses had largely waned to control levels, particularly in the fast-responding diploid Na. To remove the effect of wounding from the transcriptional signatures and to capture the earlier responses, we hybridized arrays with samples taken from locally elicited leaves 2 h after the elicitation in which mechanical wounds had been treated with either R or F, with samples taken at the same time from plants whose wounds had been treated only with water (for R), or triton-containing water (for F).

This analysis from the three different species revealed that the expression of 138 transcripts (56.3% of 246) was differentially regulated: 84% (116), 24.6% (34), and 50.7% (70) in Na, Nq, and Nc, respectively. Of these 138 regulated transcripts, 58 genes were regulated by R or F treatments in all 3 species, while 80 genes were elicited only in the allotetrapolyploids (23) or in Na (57). The comparison of R- and F-elicited transcriptional changes in all 3 species showed that R-induced responses (25%, Nq and 40.6%, Nc) from Na are more similar to those from the 2 allotetraploids than are F-induced responses (16.1%, Nq and 12.9%, Nc); moreover, the R- and F-induced responses from Nq (76.2%, R and 62.5%, F) are more similar than those from Na compared to Nc (47.3%, R and 41.4%, F). Treatment of F accounted for 64% (41) of the R-specific elicitation observed in Na. Similar analyses of F- and R-induced genes in Nc and Nq revealed that F treatment accounted for 21.8% (12 of 55) and 52.4% (11 of 21) of the total R- induced genes, respectively (Table I).

The PCA analysis indicates that the R- or F-induced transcriptional changes from Na differ distinctly from those belonging to the two allotetraploids as is readily seen along ordination axis 1, which accounts for 52.9% of the variance in the data (Fig. 3). While the R-treated plants of all three species cluster similarly along axis 1, the F treatments in Na are displaced along the vector in a positive direction, while in both Nc and Nq, the F-treated samples are displaced in a negative direction. The second axis (accounting for 16.1% of the variance) clearly distinguishes the transcriptional profiles of the two allopolyploids.

View larger version (14K): [in this window] [in a new window]   Figure 3. A PCA of the logarithmic expression ratios from six cDNA microarrays hybridized with fluorescently labeled cDNA isolated from Na, Nq, and Nc, respectively. At the rosette stage of growth, node 0 leaves from each species were wounded and immediately treated with 20 µL of water, 1:10 R, F in the R, or 0.0005% triton solution, respectively. The treated leaves were harvested for mRNA extraction 2 h after treatment. Microarrays were hybridized with cDNA reverse transcribed from mRNA pooled from 20 independently treated plants and Cy3-labeled (plants elicited with R or F) or Cy5-labeled (plants elicited with water or triton). In these hybridizations, the wound-induced responses were removed. Axes 1 and 2 of the PCA account for 79% of the variance. White circles, squares, and black squares represent Na, Nq, and Nc, respectively.

Individual Genes

Genes involved in oxylipin signaling, antimicrobial defense, remodeling of metabolism, and transcriptional machinery are shown by this study to be regulated. Continuous release of MeJA in the plants leads to significant up-regulation of genes that encode for enzymes involved in the oxylipin signaling cascade (13-lipoxygenase, LOX; allene oxide synthase, AOS; hydroperoxide lyase, HPL; alpha-dioxygenase, -DIOX) locally in all three species, while the significant up-regulation of all three genes in systemic leaves was observed only in Nq (Fig. 4). MeJA elicted a 2-fold increase of the expressions of LOX and AOS in Nq compared to those from Na, while the expressions of HPL and LOX from Nc were lower than those from Na. Transcripts involved in the oxylipin-signaling cascade were also regulated by R or F treatments in all three species, but the regulation differed in timing among the three species as the expression of HPL, LOX, and AOS in Na and Nc shows; these genes were already up-regulated by R or F at 2 h, while the up-regulation of these genes in Nq was observed only 24 h after elicitation. The expression of LOX in Nc was up-regulated only by R or F at 2 h, while its expression in Nq was up-regulated by F at 24 h. HPL showed strong up-regulation by R or F in Nc at 2 and 24 h, whereas this pattern was observed in Nq only 24 h after elicitation. In addition, the R- or F-elicited expression of dehydration-induced protein and unknown genes such as RE166 and RE172 also differed in timing among the three species, showing the fastest response in Na, followed by Nc and then Nq. In addition, a 2-fold F-induced increase of the expression of HPL, LOX, and AOS in Nq compared to that of Na and Nc was also observed.

View larger version (28K): [in this window] [in a new window]   Figure 4. Mean (± SD) of the logarithmic expression ratios from microarrays with eight replicate cDNA spots of Na sequences of: genes mediating oxylipin cascade (LOX, HPL, AOS, and -DOX); an example of JA-elicited genes (TD); and an example of a growth-related gene (RUBISCO SSU: RbcsS). Microarrays were hybridized with fluorescently labeled probes derived from Na, Nq, and Nc. At the rosette stage of growth, node 0 leaves (local) from each species were wounded and immediately treated with either 20 µL water (wound + H20), 1:10 diluted R (wound + R) or F in the R (wound + F), or directly treated with MeJA in lanolin paste (MeJA). The wound + H20, wound + R, wound + F, or MeJA samples were hybridized against the untreated control samples, while R or F samples were hybridized against the wound + water or wound + triton treatments.

Genes involved in the remodeling of transcriptional machinery and metabolism are also regulated in response to different elicitors in different species. The transcripts encoding photosynthetic enzymes (N. sylvestris mRNA for the small subunit of ribulose bisphosphate carboxylase; tobacco ribulose-1,5-bisphosphate carboxylase small subunit pseudogene) tend to be down-regulated in all three species and the expression patterns are retained up to 24 h after elicitation. RB493, which shares the sequence similarity with Zymomonas mobilis rrnB operon and 23S ribosomal RNA gene, were locally up-regulated by R/F in both Na and Nq, and systemically regulated in Nc (Fig. 4).

	DISCUSSION

TOP ABSTRACT RESULTS DISCUSSION MATERIALS AND METHODS LITERATURE CITED

 
The species Na, Nq, and Nc are morphologically similar annuals whose distributions in North America differ. Based on the caryotypes and morphology of leaf, inflorescence, and trichomes, Goodspeed deduced that Nq and Nc are members of Nicotiana section Bigelovianae which originated from amphidiploidy involving Na and an ‘alatoid’ progenitor such as N. sylvestris (Goodspeed, 1954). Phylogenetic relationships in the genus Nicotiana were investigated by parsimony analyses of plastid matK and internal transcribed spacer (ITS) regions of nuclear ribosomal DNA (nrDNA). The possible progenitor of N. section Bigelovianae is still not clear, although the patterns of relationship in Nicotiana are largely congruent (Goodspeed, 1954; Komarnyts'kyi et al., 1998; Aoki and Ito, 2000). In accord with the published data, the phylogenic analysis of ncpGS splits the N. section Acuminatae into two clades. N. linearis is placed together with N. spegazzinii and the ancestor of Na is supposed to be the common origin of N. corymbosa, N. pauciflora, N. acuminate, and N. miersii. In light of all the published data, as well as of the two ncpGS genes from Nc that belong to two clades, we can conclude that neither Nq nor Nc is autotetraploid and that an ancestor of Na is involved in the formation of both allopolyploids, Nq and Nc. This conclusion further supports Goodspeed's phylogenetic hypothesis.

N. sylvestris, according to Goodspeed's hypothesis, represents the other line of parents for Nq and Nc. Although N. sylvestris is placed in the N. section Alatae (Goodspeed, 1954) mostly based on flower morphology and distribution, its chromosomes physically resemble those of N. section Acuminatae. Both ITS and ncpGS analyses suggest that N. sylvestris is closely related to Na. Furthermore, the genomic DNA from N. sylvestris and Na hybridize to the same chromosome set in Nc and Nq. These results all suggest that species other than Na and N. sylvestris must be involved in the speciation of N. section Bigelovianae and that an ancestor of N. trigonophylla may be involved in the formation of N. section Bigelovianae. Polyploid speciation is frequently associated with adaptive radiations and polyploids can exhibit an enhanced ability to survive under unfavorable conditions (Song et al., 1995). How this occurs is not well understood, but polyploidy and its associated duplication of genes are thought to offer opportunities for the evolution of a new function in the redundant copy of the gene. To compensate for the dosage effects, recent studies have demonstrated that extra copies are frequently silenced or lost due to mutations, or otherwise epigenetically regulated (Pichersky et al., 1990; Comai et al., 2000; Soltis and Soltis, 2000; Wendel, 2000; Kashkusha et al., 2002). Unlike autotetraploidy (in which the same genome is duplicated), the allopolyploid caryotype is fixed when the ancestral genomes that are combined in the new species are sufficiently different to prevent tetrasomic recombination. Our examination of the induced transcriptional profiles from two allotetraploids, Nq and Nc, in comparison to their diploid ancestor Na provides an evolutionarily-motivated analysis which reveals three insights into ecological adaptation on the transcriptional level.

First, the examination of their transcriptional changes in response to herbivore-specific elicitor R shows that all three plants can recognize M. sexta attack, but the recognition of herbivore attack on the transcription level differs in timing among the three species. The quickest activation of gene expression is found in Na; 64 and 56 transcripts from local leaves are significantly regulated by R at 2 and 24 h. Nq, on the other hand, shows the slowest responses and these may last longer, as evidenced by the fact that 21 and 45 transcripts from local leaves are regulated by R at 2 and 24 h. In addition, herbivore-specific recognition is much faster in leaves that are locally treated as compared to the responses in systemic untreated leaves. The largest difference among transcripts occurs in Na followed by Nc, then Nq in systemic leaves, a pattern which is exactly the opposite of that which characterizes the locally-treated leaves 24 h after the elicitation.

Second, the transcriptional regulation of Na is more similar to that in Nq than in Nc, although Nq is slower to recognize herbivore attack. The comparison of R-induced transcriptional changes among the three species at 2 h revealed that 76.2% and 47.2% of R-induced transcripts from Nq and Nc are regulated similarly to those from Na. The difference in the consistency of Na-specific R recognition between two allotetraploids is not due to the plant's inability to react to R. About 47% of MeJA-elicited transcriptional changes from Na are observed in Nq and Nc, while 84% of MeJA-elicited transcripts from both tetraploids are regulated similarly to those from Na. Furthermore, the expression profiles of LOX, AOS, and HPL involved in the synthesis of jasmonic acid, as well as -DIOX, are coordinated in response to different elicitors from all three species. This coordinated regulation is correlated with induced changes in the jasmonate pool in the different species (Lou and Baldwin, 2003). Jasmonic acid, a fatty-acid-derived signaling molecule, is known to be involved in several aspects of plant biology, including pollen and seed development as well as defense against attack from microbial pathogens and insect pests (Dong, 1998; Reymond and Farmer, 1998; Bostock, 1999; Kunkel and Books, 2002; Thaler et al., 2002).

Third, the recognition of R- and F-elicitation is different among the three species. The comparisons of R- and F-elicited transcriptional changes in the three species indicate that F treatments account for a majority of R-induced transcriptional changes in Na (2 h: 64%; 24 h: 77%), but not in Nq (2 h, 52%; 24 h, 33%) or in Nc (2 h, 22%; 24 h, 52%). With its full complement of elicitors, R brings out more complex transcriptional responses, which are both positive and negative. R treatments differentially regulate F-elicited transcripts by 56% (2 h) and 47% (24 h), 54% (2 h) and 74% (24 h), and 59% (2 h) and 44% (24 h) in Na, Nq, and Nc, respectively. The proportion of R-induced transcripts accounted for by F treatments increases with the timing of elicitation in Na (2 h, 64%; 24 h, 76.8%) and Nc (2 h, 22%; 24 h, 52%) but not in Nq (2 h, 52%; 24 h, 33%), indicating that F elicits longer-lasting or slower responses in Na and Nc than does R, while in Nq, F treatment leads to fast responses.

In addition, the coordinated regulation between the up-regulated genes involved in the remodeling of transcriptional machinery (DH017, DH219, and RB493) and down-regulated photosynthesis genes is also observed in the three species. This observation is consistent with a more general response to stress in which the remobilization of limiting resources can minimize the fitness consequences of losing tissues to herbivores, or activate metabolically demanding defense responses in all three species (Karban and Baldwin, 1997; Voelckel et al., 2001; Heil and Baldwin, 2002; Strauss et al., 2002; Zangerl et al., 2002; Hui et al., 2003; Zavala et al., 2004).

The ecological significance of the retention of MeJA and larval R recognition abilities in the two allopolyploid species remains unknown, but it likely provides important insights into the maintenance or modification of polygenic adaptive traits during allopolyploid speciation. Nq and Nc descended from an ancestor of modern Na due to polyploidy that happened several million years ago. The differences among the levels of consistency of Na-specific MeJA-, R- or F-induced transcriptional responses in the two allopolyploids may be due to the chromosomal rearrangements that occurred during speciation (allopolyploidy) as has been shown in wheat (Ozkan et al., 2001; He et al., 2003). Alternatively or perhaps in addition, the differences may be due to the fact that Na, Nq, and Nc are distributed in different habitats and therefore are likely subjected to different selection pressures. Nq is found in coastal California; Nc has a more southern distribution in both coastal and interior Baja California and southern California; Na is found throughout the Great Basin Desert and adjacent areas as an ephemeral member of the post-fire community in sagebrush and pinyon-juniper habitats. Hence, differences between the levels of consistency of Na-specific induced transcription profiles in two allopolyploids may also be due to the effects of natural selection during the interval since the speciation occurred. It is important to note that the Na genome has likely undergone evolutionary change since the hybridization events. Since the genes on the array are derived from a modern Na genome, some of the differences in transcriptional signatures among the species may reflect differences between the modern Na cDNAs and the Na genome that were incorporated into the other species at the time of speciation, which has not likely been subjected to the same selection pressures. Moreover, the differences in the speed of responses among the species are inferred from analyzing responses at a few time points. It is quite possible that, in nature, all the genes that respond to these treatments in Na are also capable of responding in the other species in the context of their own natural habitats. The experiments were conducted under the same growth conditions and with plants in the same developmental stage, but these growth conditions have been optimized for Na and environmental signals important for the regulation of Nc and Nq genes may have been lacking.

These evolutionary hypotheses can be falsified by determining whether the M. sexta-induced responses in Na, Nq, and Nc are indeed defensive (e.g. increase plant fitness in environments with herbivores) and by comparing the responses in newly created allopolyploid species (synthetic hybrids from Na and N. trigonophylla). Furthermore, the fact that larger-scale transcriptional changes can be elicited by either R or F and the coordinated regulation of genes involved in the signaling and remodeling of the metabolism in all three species point to the existence of unknown trans-activating factors and corresponding cis-elements. These elements may be involved in the reorganization of the complex herbivore-specific transcription changes in plants and may have been under strong selection in the allotetraploids.

	MATERIALS AND METHODS

TOP ABSTRACT RESULTS DISCUSSION MATERIALS AND METHODS LITERATURE CITED

 

Plant Growth

Nicotiana sylvestris, Nicotiana corymbosa, Nicotiana pauciflora, Nicotiana acuminata, Nicotiana linearis, Nicotiana spegazzinii, Nq, Nicotiana trigonophylla, Nicotiana miersii, and Nc seeds were kindly supplied by Dr. Verne A. Sisson (Oxford Tobacco Research Station, Oxford, NC) and originated from collections made by Dr. H. Goodspeed (Goodspeed, 1954). Na Torr. Ex Watts (an inbred line collected from the DI Ranch, Utah (Baldwin et al., 1998) seeds were germinated in smoke-treated soil and plants were grown in individual 1-L pots as previously described (Kruegel et al., 2002). All plants were grown in the glasshouse at 26°C to 28°C under 16 h light supplemented by Philips Sun-T Agro 400 or 600 W Na lights (Eindhoven, The Netherlands).

PCR Amplification, Cloning, and Sequencing of ncpGS

Using young leaves, total genomic DNA was extracted according to the cetyl-trimethyl-ammonium bromide method (Richards, 1997). Two pairs of primers GS57 (5'-AGA C/T CC TTT CCG TGG A/G/C/T GG-3'), GS2 (5'-CAG CTC CGG TTC CAC AGT AGT-3') and NQ 40 (5'-TAA AGG CCA CTT CAC ATT TT-3'), NQ 41 (5'-TGA ATG TTT GAT CTC ATT GG-3') were designed to amplify Gln synthetase isoforms. Amplification was performed in the total volume of 50 µL including 10 ng of template genomic DNA, 1x PCR buffer, 1.5 mM MgCl₂, 0.2 µM of each primer, 200 µM of each dNTP, and 1.25 units of Taq DNA polymerase. The PCR conditions were as follows: an initial denaturation step at 94°C for 5 min, 30 cycles at 94°C for 1 min, 57°C for 1 min (50°C during the first two cycles for the primer pair GS57/GS2) and 72°C for 2 min, and a final extension step at 72°C for 10 min. The PCR fragments were gel-purified and cloned into pCR2.1 TOPO (in Escherichia coli TOP10F') or pGEM-7Zf(+) vector. Sequencing was conducted using an ABI PRISM 377 automated sequencer (Sunnyvale, CA). Each clone was sequenced in both sense and antisense directions and at least four clones were sequenced for each fragment.

Treatments

For microarray analysis, the node 0 leaves (source-sink transition) of 20 individual, developmentally synchronized plants per treatment in the rosette stage of growth from Na, Nq, and Nc, were wounded by creating three rows of puncture wounds on each leaf half (Ohnmeiss and Baldwin, 1994). Wounds were immediately treated with either 20 µL of deionized water, M. sexta R (diluted in 1:10 deionized water) from fourth to fifth instar larvae, or a synthetic mixture of the two most abundant FACs in M. sexta R: N-linolenoyl-L-Gln and N-linolenoyl-L-Glu (0.012 mM N-linolenoyl-L-Gln and 0.034 mM N-linolenoyl-L-Glu). Since the FACs were dissolved in R-equivalent concentrations in an aqueous solution of 0.005% Triton X-100, we applied 0.0005% (1:10 dilution) triton solution to the wounds to control the possible effects of triton. Nonwounded plants were treated with 150 µg MeJA in 20 µL of lanolin paste, which was applied to adaxial lamina. The treated or next younger position leaves (position –1) were harvested at 2 or 24 h after the treatment and frozen in liquid nitrogen. The identical leaves from uninduced plants were harvested at 24 h as control.

cDNA-Microarray Analysis

Two PCR-amplified fragments from each of 246 herbivore-regulated genes isolated by either differential-display reverse transcription-PCR or cDNA-amplified fragment-length polymorphism from Na were amino-linked on either strand and spotted on epoxy-coated slides as previously described (Hermsmeier et al., 2001; Halitschke et al., 2003; Hui et al., 2003). Leaf samples were homogenized in liquid nitrogen and the total RNA was isolated according to the methods described in Hermsmeier et al. (2001). Twenty-four cDNA microarrays were competitively hybridized with the Cy3-labeled cDNA from a 24 h harvest of treated node 0 leaves or –1 leaves and Cy5-labeled cDNA generated from identical leaves from uninduced plants. To monitor the herbivore-specific recognition responses without the potentially confounding effects of the wound response, we hybridized an additional six microarrays with Cy3-labeled cDNA from leaves from the three species elicited by treating wounds immediately with either R or FACs for 2 h and Cy5-labeled cDNA from leaves elicited by treating wounds immediately with water or triton-containing solution. Hence, in these six arrays, the wound-response was removed from the analysis.

The cDNA microarrays were scanned and evaluated with methods described and verified in previous publications (Halitschke et al., 2003; Hui et al., 2003). Briefly, an array-specific normalization factor was calculated from the middle 75% of the distribution of expression ratios (ERs). The ratios of normalized ERs of each individual spot and the mean of the four replicate spots for each cDNA (two for each gene) were calculated. A transcript was defined as being differentially regulated only if the following three criteria were met: the sum of normalized background corrected Cy3 and Cy5 signals was higher than 1,000; both individual ERs were equal to or exceeded the arbitrary thresholds for differential expression (1.25 and 0.75), respectively; and both individual expression ratios were significantly different from one as determined by t test analysis of the eight replicate spots. The expressions of genes are summarized in Supplemental Table I in which the expression patterns are represented as up, down, or -, depending on whether or not the three criteria are met.

To compare the full transcriptional responses of three species to different treatments at different harvest times, the logarithmic expression ratios of all transcripts from all 30 cDNA-microarrays were analyzed by detrended correspondence analysis. Since the longest gradient from DCA analysis was shorter than 3, a principal components analysis, a linear type of ordination model, was performed to examine the interarrays distances (Canoco for Windows 4.5, Microcomputer Power, Ithaca, NY; Yeung and Ruzzo, 2001; ter Braak and Smilauer, 2002).

Sequence data from this article have been deposited with the EMBL/GenBank data libraries under accession numbers AF302113, AY566989, AY568077, AY568079, AY568078, AY568082, AY568081, AY568084, AY568083, AY568080, AY566990, AY566991, and AY183657.

	ACKNOWLEDGMENTS

 
We thank R. Ropte and W. Kröber for technical assistant, R. Halitschke for providing the FACs, S. Kutschbach, K. Gase, T. Hahn, M. Held, C. Voelckel, and C. McInerney for invaluable assistance in microarray hybridization, reading, and data analysis, E. Wheeler for editorial assistance.

Received November 30, 2003; returned for revision March 10, 2004; accepted March 10, 2004.

	FOOTNOTES

 
¹ This work was funded by Max Planck Gesellschaft and the Deutsche Forschungsgemeinschaft (SPP1152).

² Present address: PHARMAPLANT GmbH, Am Westbahnhof 4, 06556 Artern, Germany.

^[w] The online version of this article contains Web-only data.

Article, publication date, and citation information can be found at www.plantphysiol.org/cgi/doi/10.1104/pp.103.037036.

^* Corresponding author; e-mail address baldwin{at}ice.mpg.de; fax 49–3641–571102.

	LITERATURE CITED

TOP ABSTRACT RESULTS DISCUSSION MATERIALS AND METHODS LITERATURE CITED

 
Aoki S, Ito M (2000) Molecular phylogeny of Nicotiana (solanaceae) based on the nucleotid sequence of the matK gene. Plant Biol 2: 316–324[CrossRef]

Baldwin IT, Gorham D, Schmelz EA, Lewandowski C, Lynds GY (1998) Allocation of nitrogen to an inducible defense and seed production in Nicotiana attenuata. Oecologia 115: 541–552[CrossRef]

Baldwin IT, Preston CA (1999) The eco-physiological complexity of plant responses to insect. Planta 208: 137–145[CrossRef]

Baldwin IT, Schmelz EA, Ohnmeiss TE (1994) Wound-induced changes in root and shoot jasmonic acid pools correlate with induced nicotine synthesis in Nicotiana sylvestris Spegazzini and Comes. J Chem Ecol 20: 2139–2157[CrossRef]

Bostock RM (1999) Signal conflicts and synergies in induced resistance to multiple attackers. Physiol Mol Plant Pathol 55: 99–109[CrossRef]

Chase MW, Knapp S, Cox AV, Clarkson JJ, Butsko Y, Joseph J, Savolainen V, Parokonny AS (2003) Molecular systematics, GISH and the origin of hybrid taxa in Nicotiana (Solanaceae). Ann Bot (Lond) 92: 107–127[Abstract/Free Full Text]

Comai L, Tyagi AP, Winter K, Holmes-Davis R, Reynolds SH, Stevens Y, Byers B (2000) Phenotypic instability and rapid gene silencing in newly formed Arabidopsis allotetraploids. Plant Cell 12: 1551–1568[Abstract/Free Full Text]

Dong X (1998) SA, JA, ethylene, and disease resistance in plants. Curr Opin Plant Biol 1: 316–323[CrossRef][Web of Science][Medline]

Emshwiller E, Doyle JJ (1999) Chloroplast-expressed glutamine synthetase (ncpGS): potential utility for phylogenetic studies with an example from Oxalis (Oxalidaceae). Mol Phylogenet Evol 12: 310–319[CrossRef][Web of Science][Medline]

Glawe GA, Zavala JA, Kessler A, Van Dam NM, Baldwin IT (2003) Ecological costs and benefits correlated with trypsin protease inhibitor production in Nicotiana attenuata. Ecology 84: 79–90[CrossRef][Web of Science]

Goodspeed TH (1954) The Genus Nicotiana. Chronica Botanic, Waltham, MA

Halitschke R, Gase K, Hui D, Schmidt D, Baldwin IT (2003) Molecular interactions between the specialist herbivore Manduca sexta (Lepidoptera, Sphingidae) and its natural host Nicotiana attenuata. VI. Microarray analysis reveals that most herbivore-specific transcriptional changes are mediated by fatty acid-amino acid conjugates. Plant Physiol 131: 1894–1902[Abstract/Free Full Text]

Halitschke R, Schittko U, Pohnert G, Boland W, Baldwin IT (2001) Molecular interactions between the specialist herbivore Manduca sexta (Lepidoptera, Sphingidae) and its natural host Nicotiana attenuata. III. Fatty acid-amino acid conjugates in herbivore oral secretions are necessary and sufficient for herbivore-specific plant responses. Plant Physiol 125: 711–717[Abstract/Free Full Text]

He P, Friebe BR, Gill BS, Zhou JM (2003) Allopolyploidy alters gene expression in the highly stable hexaploid wheat. Plant Mol Biol 52: 401–414[CrossRef][Web of Science][Medline]

Heil M, Baldwin IT (2002) Fitness costs of induced resistance: emerging experimental support for a slippery concept. Trends Plant Sci 7: 61–67[CrossRef][Web of Science][Medline]

Hermsmeier D, Schittko U, Baldwin IT (2001) Molecular interactions between the specialist herbivore Manduca sexta (Lepidoptera, Sphingidae) and its natural host Nicotiana attenuata. I. Large-scale changes in the accumulation of growth- and defense-related plant mRNAs. Plant Physiol 125: 683–700[Abstract/Free Full Text]

Hieter J, Griffiths T (1999) Polyploidy-more is more or less. Science 285: 210–211[Free Full Text]

Hilu KW (1993) Polyploidy and the evolution of domesticated plants. Am J Bot 80: 1494–1499[CrossRef]

Hui D, Iqbal J, Lehmann K, Gase K, Saluz HP, Baldwin IT (2003) Molecular interactions between the specialist herbivore Manduca sexta (Lepidoptera, Sphingidae) and its natural host Nicotiana attenuata. V. Microarray analysis and further characterization of large-scale changes in herbivore-induced mRNAs. Plant Physiol 131: 1877–1893[Abstract/Free Full Text]

Janz N, Thompson JN (2002) Plant polyploidy and host expansion in an insect herbivore. Oecologia 130: 570–575[CrossRef]

Karban R, Baldwin IT (1997) Induced Response to Herbivory. Chicago University Press, Chicago

Kashkusha K, Feldmana M, Levya AA (2002) Gene loss, silencing and activation in a newly synthesized wheat allotetraploid. Genetics 160: 1651–1659[Abstract/Free Full Text]

Kellogg EA (1998) Relationships of cereal crops and other grasses. Proc Natl Acad Sci USA 95: 2005–2010[Abstract/Free Full Text]

Kessler A, Baldwin IT (2001) Defensive function of herbivore-induced plant volatile emissions in nature. Science 291: 2141–2144[Abstract/Free Full Text]

Kessler A, Baldwin IT (2002) Manduca quinquemaculata's optimization of intra-plant oviposition to predation, food quality, and thermal constraints. Ecology 83: 2346–2354[CrossRef][Web of Science]

Komarnyts'kyi SI, Komarnyts'kyi IK, Cox A, Parokonnyi AS (1998) The evolution of the sequences of the internal spacer of nuclear ribosomal DNA for American species in the genus Nicotiana. Tsitol Genet 32: 69–76[Medline]

Kroymann J, Donnerhacke S, Schnabelrauch D, Mitchell-Olds T (2003) Evolutionary dynamics of an Arabidopsis insect resistance quantitative trait locus. Proc Natl Acad Sci USA 100: 14587–14592[Abstract/Free Full Text]

Kruegel T, Lim M, Gase K, Halitschke R, Baldwin IT (2002) Agrobacterium-mediated transformation of Nicotiana attenuata, a model ecological expression system. Chemoecology 12: 177–183

Kunkel BN, Books DM (2002) Cross talk between signaling pathways in pathogen defense. Curr Opin Plant Biol 5: 325–331[CrossRef][Web of Science][Medline]

Lou YG, Baldwin IT (2003) Manduca sexta recognition and resistance among allopolyploid Nicotiana host plants. Proc Natl Acad Sci USA 100: 14513–14600[Free Full Text]

Masterson J (1994) Stomatal size in fossil plants: evidence for polyploidy in majority of angiosperms. Science 264: 421–424[Abstract/Free Full Text]

Mattiacci L, Dicke M, Posthumus MA (1995) Glucosidase: an elicitor of herbivore-induced plant odor that attracts host-searching parasitic wasps. Proc Natl Acad Sci USA 92: 2036–2040[Abstract/Free Full Text]

Musser RO, Hum-Musser SM, Eichenseer H, Peiffer M, Ervin G, Murphy JB, Felton GW (2002) Caterpillar saliva beats plant defences: a new weapon emerges in the evolutionary arms race between plants and herbivores. Nature 416: 599–600[CrossRef][Medline]

Nuismer SL, Thompson JN (2001) Plant polyploidy and non-uniform effects on insect herbivores. Proc R Soc Lond Ser B Biol Sci 268: 1937–1940[Medline]

Ohnmeiss TE, Baldwin IT (1994) The allometry of nitrogen allocation to growth and an inducible defense under nitrogen-limited growth. Ecology 75: 995–1002[CrossRef]

Ozkan H, Levy AA, Feldman M (2001) Allopolyploidy-induced rapid genome evolution in the wheat (Aegilops-Triticum) group. Plant Cell 13: 1735–1747[Abstract/Free Full Text]

Pichersky E, Soltis D, Soltis P (1990) Defective chlorophyll a/b-binding protein genes in the genome of a homosporous fern. Proc Natl Acad Sci USA 87: 195–199[Abstract/Free Full Text]

Pohnert GV, Jung E, Haukioja KL, Boland W (1999) New fatty acid amides from regurgitant of lepidopteran (Noctuidae, Geometridae) caterpillars. Tetrahedron Lett 55: 11275–11280

Reymond P, Farmer EE (1998) Jasmonate and salicylate as global signals for defense gene expression. Curr Opin Plant Biol 1: 404–411[CrossRef][Web of Science][Medline]

Richards EJ (1997) Preparation of plant DNA using CTAB. In F Ausuble, R Brent, RE Kingston, DD Moore, JG Seidman, JA Smith, K Struhl, eds, Short Protocols in Molecular Biology. John Wiley & Sons, New York, pp 2.10–2.11

Schittko U, Hermsmeier D, Baldwin IT (2001) Molecular interactions between the specialist herbivore Manduca sexta (Lepidoptera, Sphingidae) and its natural host Nicotiana attenuata. II. Accumulation of plant mRNAs in response to insect-derived cues. Plant Physiol 125: 701–710[Abstract/Free Full Text]

Schittko U, Preston CA, Baldwin IT (2000) Eating the evidence? Manduca sexta larvae can not disrupt specific jasmonate induction in Nicotiana attenuata by rapid consumption. Planta 210: 343–346[CrossRef][Web of Science][Medline]

Schoen DJ, Burdon JJ, Brown AHD (1992) Resistance of Glycine tomentella to soybean leaf rust Phakopsora pachyrhizi in relation to ploidy level and geographic distribution. Theor Appl Genet 83: 827–832

Shoemaker RC, Polzin K, Labate J, Specht J, Brummer EC, Olson T, Young N, Concibido V, Wilcox J, Tamulonis JP, et al. (1996) Genome duplication in soybean (Glycine subgenus soja). Genetics 144: 329–338[Abstract]

Soltis PS, Soltis DE (2000) The role of genetic and genomic attributes in the success of polyploids. Proc Natl Acad Sci USA 20: 7051–7057

Song K, Lu P, Tang K, Osborn TC (1995) Rapid genome change in synthetic polyploids of Brassica and its implications for polyploid evolution. Proc Natl Acad Sci USA 92: 7719–7723[Abstract/Free Full Text]

Strauss SY, Rudgers JA, Lau JA, Irwin RE (2002) Direct and ecological costs of resistance to herbivory. Trends Ecol Evol 17: 278–285[CrossRef]

ter Braak CJF, Smilauer P (2002) Canoco Reference Manual and CanoDraw for Windows User's Guide, Version 4.5. Microcomputer Power, Ithaca, NY

Thaler JS, Karban R, Ullman DE, Boege K, Bostock RM (2002) Cross-talk between jasmonate and salicylate plant defense pathways: effects on several plant parasites. Oecologia 131: 227–235

Thompson JN, Cunningham BM, Segraves KA, Althoff DM, Wagner D (1997) Plant polyploidy and insect/plant interaction. Am Nat 150: 731–743

Voelckel C, Baldwin IT (2003) Detecting herbivore-specific transcriptional responses in plants with multiple DDRT-PCR and subtractive library procedures. Physiol Plant 118: 240–252[CrossRef]

Voelckel C, Schittko U, Baldwin IT (2001) Herbivore-induced ethylene burst reduces fitness costs of jasmonate- and oral secretion-induced defenses in Nicotiana attenuata. Oecologia 127: 274–280[CrossRef]

Wendel JF (2000) Genome evolution in polyploids. Plant Mol Biol 42: 225–249[CrossRef][Web of Science][Medline]

Yeung KY, Ruzzo WL (2001) Principal component analysis for clustering gene expression data. Bioinformatics 17: 763–774[Abstract/Free Full Text]

Zangerl AR, Hamilton JG, Miller TJ, Crofts AR, Oxborough K, Berenbaum MR, de Lucia EH (2002) Impact of folivory on photosynthesis is greater than the sum of its holes. Proc Natl Acad Sci USA 99: 1088–1091[Abstract/Free Full Text]

Zavala JA, Patankar AG, Gase K, Baldwin IT (2004) Constitutive and inducible trypsin proteinase inhibitor production incurs large fitness costs in Nicotiana attenuata. Proc Natl Acad Sci USA 101: 1607–1612[Abstract/Free Full Text]

Zhang ZP, Baldwin IT (1997) Transport of [2-¹⁴C]jasmonic acid from leaves to roots mimics wound-induced changes in endogenous jasmonic acid pools in Nicotiana sylvestris. Planta 203: 436–441[CrossRef][Web of Science]

This article has been cited by other articles:

				S. Anssour, T. Krugel, T. F. Sharbel, H. P. Saluz, G. Bonaventure, and I. T. Baldwin Phenotypic, genetic and genomic consequences of natural and synthetic polyploidization of Nicotiana attenuata and Nicotiana obtusifolia Ann. Bot., June 1, 2009; 103(8): 1207 - 1217. [Abstract] [Full Text] [PDF]

				G. Vandenborre, O. Miersch, B. Hause, G. Smagghe, C. Wasternack, and E. J.M. Van Damme Spodoptera littoralis-Induced Lectin Expression in Tobacco Plant Cell Physiol., June 1, 2009; 50(6): 1142 - 1155. [Abstract] [Full Text] [PDF]

				M. J Hegarty, G. L Barker, A. C Brennan, K. J Edwards, R. J Abbott, and S. J Hiscock Changes to gene expression associated with hybrid speciation in plants: further insights from transcriptomic studies in Senecio Phil Trans R Soc B, September 27, 2008; 363(1506): 3055 - 3069. [Abstract] [Full Text] [PDF]

				M. Skibbe, N. Qu, I. Galis, and I. T. Baldwin Induced Plant Defenses in the Natural Environment: Nicotiana attenuata WRKY3 and WRKY6 Coordinate Responses to Herbivory PLANT CELL, July 1, 2008; 20(7): 1984 - 2000. [Abstract] [Full Text] [PDF]

				I. J. Leitch, L. Hanson, K. Y. Lim, A. Kovarik, M. W. Chase, J. J. Clarkson, and A. R. Leitch The Ups and Downs of Genome Size Evolution in Polyploid Species of Nicotiana (Solanaceae) Ann. Bot., April 1, 2008; 101(6): 805 - 814. [Abstract] [Full Text] [PDF]

This Article Abstract Full Text (PDF) Supplemental Data All Versions of this Article: 135/1/539    most recent pp.103.037036v1 Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Similar articles in Web of Science Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via CrossRef Citing Articles via Web of Science (12) Citing Articles via Google Scholar Google Scholar Articles by Qu, N. Articles by Baldwin, I. T. Search for Related Content PubMed PubMed Citation Articles by Qu, N. Articles by Baldwin, I. T. Agricola Articles by Qu, N. Articles by Baldwin, I. T.