Natural genetic variation and negative density effects in plant-nematode interactions.

Arabidopsis thaliana is a suitable host for phytoparasitic nematodes of the genus Meloidogyne. Successful nematode infection leads to the formation of root galls. We tested for natural genetic variation and inoculation density effects on nematode reproductive success in the interaction between A. thaliana and Meloidogyne javanica. We inoculated different Arabidopsis genotypes with two sources of nematodes at two different doses, using a mild protocol for inoculum preparation. We counted root galls and egg masses 2 months after inoculation. We obtained a high number of successful nematode infections. Infection success differed among Arabidopsis genotypes in interaction with the nematode source. Overall, infection success and reproductive success of nematodes were lower at a higher inoculum dose of nematodes. Our results indicate that natural genetic variation in both host plants and nematodes, as well as short- and long-term negative density effects, shape nematode reproductive success.

Microevolution of Pieris butterfly genes involved in host plant adaptation along a host plant community cline.

Herbivorous insects have evolved counteradaptations to overcome the chemical defences of their host plants. Several of these counteradaptations have been elucidated at the molecular level, in particular for insects specialized on cruciferous host plants. While the importance of these counteradaptations for host plant colonization is well established, little is known about their microevolutionary dynamics in the field. In particular, it is not known whether and how host plant diversity shapes diversity in insect counteradaptations. In this study, we examine patterns of host plant use and insect counteradaptation in three Pieris butterfly species across Japan. The larvae of these butterflies express nitrile-specifier protein (NSP) and its paralogue major allergen (MA) in their gut to overcome the highly diversified glucosinolate-myrosinase defence system of their cruciferous host plants. Pieris napi and Pieris melete colonize wild Brassicaceae whereas Pieris rapae typically uses cultivated Brassica as a host, regardless of the local composition of wild crucifers. As expected, NSP and MA diversity was independent of the local composition of wild Brassicaceae in P. rapae. In contrast, NSP diversity correlated with local host plant diversity in both species that preferred wild Brassicaceae. Both P. melete and P. napi revealed two distinct major NSP alleles, which shaped diversity among local populations, albeit with different evolutionary trajectories. In comparison, MA showed no indication for local adaptation. Altogether, MA appeared to be evolutionary more conserved than NSP, suggesting that both genes play different roles in diverting host plant chemical defence.

The two copies of the zinc and cadmium ZIP6 transporter of Arabidopsis halleri have distinct effects on cadmium tolerance.

Plants have the ability to colonize highly diverse environments. The zinc and cadmium hyperaccumulator Arabidopsis halleri has adapted to establish populations on soils covering an extreme range of metal availabilities. The A. halleri ZIP6 gene presents several hallmarks of hyperaccumulation candidate genes: it is constitutively highly expressed in roots and shoots and is associated with a zinc accumulation quantitative trait locus. Here, we show that AhZIP6 is duplicated in the A. halleri genome. The two copies are expressed mainly in the vasculature in both A. halleri and Arabidopsis thaliana, indicative of conserved cis regulation, and acquired partial organ specialization. Yeast complementation assays determined that AhZIP6 is a zinc and cadmium transporter. AhZIP6 silencing in A. halleri or expression in A. thaliana alters cadmium tolerance, but has no impact on zinc and cadmium accumulation. AhZIP6-silenced plants display reduced cadmium uptake upon short-term exposure, adding AhZIP6 to the limited number of Cd transporters supported by in planta evidence. Altogether, our data suggest that AhZIP6 is key to fine-tune metal homeostasis in specific cell types. This study additionally highlights the distinct fates of duplicated genes in A. halleri.

[Biosynthesis of indole glucosinolates and ecological role of secondary modification pathways]. / Biosynthèse des glucosinolates indoliques et rôle écologique de leurs modifications secondaires.

Indole glucosinolates are plant secondary metabolites derived from the amino acid tryptophan. They are part of a large group of sulfur-containing molecules almost exclusively found among Brassicales, which include the mustard family (Brassicaceae) with many edible plant species of major nutritional importance. These compounds mediate numerous interactions between these plants and their natural enemies and are therefore of major biological and economical interest. This literature review aims at taking stock of recent advances of our knowledge about the biosynthetic pathways of indole glucosinolates, but also about the defense strategies and ecological processes involving these metabolites.

An Insect Counteradaptation against Host Plant Defenses Evolved through Concerted Neofunctionalization.

Antagonistic chemical interactions between herbivorous insects and their host plants are often thought to coevolve in a stepwise process, with an evolutionary innovation on one side being countered by a corresponding advance on the other. Glucosinolate sulfatase (GSS) enzyme activity is essential for the Diamondback moth, Plutella xylostella, to overcome a highly diversified secondary metabolite-based host defense system in the Brassicales. GSS genes are located in an ancient cluster of arylsulfataselike genes, but the exact roles of gene copies and their evolutionary trajectories are unknown. Here, we combine a functional investigation of duplicated insect arylsulfatases with an analysis of associated nucleotide substitution patterns. We show that the Diamondback moth genome encodes three GSSs with distinct substrate spectra and distinct expression patterns in response to glucosinolates. Contrary to our expectations, early functional diversification of gene copies was not indicative of a coevolutionary arms race between host and herbivore. Instead, both copies of a duplicated arylsulfatase gene evolved concertedly in the context of an insect host shift to acquire novel detoxifying functions under positive selection, a pattern of duplicate gene retention that we call "concerted neofunctionalization."

Methyl Transfer in Glucosinolate Biosynthesis Mediated by Indole Glucosinolate O-Methyltransferase 5.

Indole glucosinolates (IGs) are plant secondary metabolites that are derived from the amino acid tryptophan. The product of Arabidopsis (Arabidopsis thaliana) IG core biosynthesis, indol-3-ylmethyl glucosinolate (I3M), can be modified by hydroxylation and subsequent methoxylation of the indole ring in position 1 (1-IG modification) or 4 (4-IG modification). Products of the 4-IG modification pathway mediate plant-enemy interactions and are particularly important for Arabidopsis innate immunity. While CYP81Fs encoding cytochrome P450 monooxygenases and IGMTs encoding indole glucosinolate O-methyltransferases have been identified as key genes for IG modification, our knowledge about the IG modification pathways is not complete. In particular, it is unknown which enzyme is responsible for methyl transfer in the 1-IG modification pathway and whether this pathway plays a role in defense, similar to 4-IG modification. Here, we analyze two Arabidopsis transfer DNA insertion lines with targeted metabolomics. We show that biosynthesis of 1-methoxyindol-3-ylmethyl glucosinolate (1MOI3M) from I3M involves the predicted unstable intermediate 1-hydroxyindol-3-ylmethyl glucosinolate (1OHI3M) and that IGMT5, a gene with moderate similarity to previously characterized IGMTs, encodes the methyltransferase that is responsible for the conversion of 1OHI3M to 1MOI3M. Disruption of IGMT5 function increases resistance against the root-knot nematode Meloidogyne javanica and suggests a potential role for the 1-IG modification pathway in Arabidopsis belowground defense.

The N-end rule pathway regulates pathogen responses in plants.

To efficiently counteract pathogens, plants rely on a complex set of immune responses that are tightly regulated to allow the timely activation, appropriate duration and adequate amplitude of defense programs. The coordination of the plant immune response is known to require the activity of the ubiquitin/proteasome system, which controls the stability of proteins in eukaryotes. Here, we demonstrate that the N-end rule pathway, a subset of the ubiquitin/proteasome system, regulates the defense against a wide range of bacterial and fungal pathogens in the model plant Arabidopsis thaliana. We show that this pathway positively regulates the biosynthesis of plant-defense metabolites such as glucosinolates, as well as the biosynthesis and response to the phytohormone jasmonic acid, which plays a key role in plant immunity. Our results also suggest that the arginylation branch of the N-end rule pathway regulates the timing and amplitude of the defense program against the model pathogen Pseudomonas syringae AvrRpm1.

Metabolism, excretion and avoidance of cyanogenic glucosides in insects with different feeding specialisations.

Cyanogenic glucosides (CNglcs) are widespread plant defence compounds releasing toxic hydrogen cyanide when hydrolysed by specific ß-glucosidases after plant tissue damage. In contrast to specialist herbivores that have mechanisms to avoid toxicity from CNglcs, it is generally assumed that non-adapted herbivores are negatively affected by CNglcs. Recent evidence, however, implies that the defence potential of CNglcs towards herbivores may not be as effective as previously anticipated. Here, performance, metabolism and excretion products of insects not adapted to CNglcs were analysed, including species with different degrees of dietary specialisation (generalists, specialists) and different feeding modes (leaf-snipping lepidopterans, piercing-sucking aphids). Insects were reared either on cyanogenic or acyanogenic plants or on an artificial cyanogenic diet. Lepidopteran generalists (Spodoptera littoralis, Spodoptera exigua, Mamestra brassicae) were compared to lepidopteran glucosinolate-specialists (Pieris rapae, Pieris brassicae, Plutella xylostella), and a generalist aphid (Myzus persicae) was compared to an aphid glucosinolate-specialist (Lipaphis erysimi). All insects were tolerant to cyanogenic plants; in lepidopterans tolerance was mainly due to excretion of intact CNglcs. The two Pieris species furthermore metabolized aromatic CNglcs to amino acid conjugates (Cys, Gly, Ser) and derivatives of these, which is similar to the metabolism of benzylglucosinolates in these species. Aphid species avoided uptake of CNglcs during feeding. Our results imply that non-adapted insects tolerate plant CNglcs either by keeping them intact for excretion, metabolizing them, or avoiding uptake.

Hard selective sweep and ectopic gene conversion in a gene cluster affording environmental adaptation.

Among the rare colonizers of heavy-metal rich toxic soils, Arabidopsis halleri is a compelling model extremophile, physiologically distinct from its sister species A. lyrata, and A. thaliana. Naturally selected metal hypertolerance and extraordinarily high leaf metal accumulation in A. halleri both require Heavy Metal ATPase4 (HMA4) encoding a PIB-type ATPase that pumps Zn(2+) and Cd(2+) out of specific cell types. Strongly enhanced HMA4 expression results from a combination of gene copy number expansion and cis-regulatory modifications, when compared to A. thaliana. These findings were based on a single accession of A. halleri. Few studies have addressed nucleotide sequence polymorphism at loci known to govern adaptations. We thus sequenced 13 DNA segments across the HMA4 genomic region of multiple A. halleri individuals from diverse habitats. Compared to control loci flanking the three tandem HMA4 gene copies, a gradual depletion of nucleotide sequence diversity and an excess of low-frequency polymorphisms are hallmarks of positive selection in HMA4 promoter regions, culminating at HMA4-3. The accompanying hard selective sweep is segmentally eclipsed as a consequence of recurrent ectopic gene conversion among HMA4 protein-coding sequences, resulting in their concerted evolution. Thus, HMA4 coding sequences exhibit a network-like genealogy and locally enhanced nucleotide sequence diversity within each copy, accompanied by lowered sequence divergence between paralogs in any given individual. Quantitative PCR corroborated that, across A. halleri, three genomic HMA4 copies generate overall 20- to 130-fold higher transcript levels than in A. thaliana. Together, our observations constitute an unexpectedly complex profile of polymorphism resulting from natural selection for increased gene product dosage. We propose that these findings are paradigmatic of a category of multi-copy genes from a broad range of organisms. Our results emphasize that enhanced gene product dosage, in addition to neo- and sub-functionalization, can account for the genomic maintenance of gene duplicates underlying environmental adaptation.

Comparative genomic and transcriptomic characterization of the toxigenic marine dinoflagellate Alexandrium ostenfeldii.

Many dinoflagellate species are notorious for the toxins they produce and ecological and human health consequences associated with harmful algal blooms (HABs). Dinoflagellates are particularly refractory to genomic analysis due to the enormous genome size, lack of knowledge about their DNA composition and structure, and peculiarities of gene regulation, such as spliced leader (SL) trans-splicing and mRNA transposition mechanisms. Alexandrium ostenfeldii is known to produce macrocyclic imine toxins, described as spirolides. We characterized the genome of A. ostenfeldii using a combination of transcriptomic data and random genomic clones for comparison with other dinoflagellates, particularly Alexandrium species. Examination of SL sequences revealed similar features as in other dinoflagellates, including Alexandrium species. SL sequences in decay indicate frequent retro-transposition of mRNA species. This probably contributes to overall genome complexity by generating additional gene copies. Sequencing of several thousand fosmid and bacterial artificial chromosome (BAC) ends yielded a wealth of simple repeats and tandemly repeated longer sequence stretches which we estimated to comprise more than half of the whole genome. Surprisingly, the repeats comprise a very limited set of 79-97 bp sequences; in part the genome is thus a relatively uniform sequence space interrupted by coding sequences. Our genomic sequence survey (GSS) represents the largest genomic data set of a dinoflagellate to date. Alexandrium ostenfeldii is a typical dinoflagellate with respect to its transcriptome and mRNA transposition but demonstrates Alexandrium-like stop codon usage. The large portion of repetitive sequences and the organization within the genome is in agreement with several other studies on dinoflagellates using different approaches. It remains to be determined whether this unusual composition is directly correlated to the exceptionally genome organization of dinoflagellates with a low amount of histones and histone-like proteins.

Natural diversity and adaptation in plant secondary metabolism.

Technological advances in metabolomics, transcriptomics and genomics have facilitated the detection of genes that contribute to diversification in plant secondary metabolism. Statistical tools from molecular population genetics may help in evaluating whether the corresponding genes or genomic regions carry a signature of selection and answering the question of whether novel compounds are 'adaptive'. Gene duplication fuels diversification in plant secondary metabolism and the evolutionary mechanism for adaptation may follow a path of neofunctionalization subsequent to gene duplication, or gene duplication may occur subsequent to--and resolve--an adaptive conflict present in a single ancestral gene sequence.

Metabolic engineering in Nicotiana benthamiana reveals key enzyme functions in Arabidopsis indole glucosinolate modification.

Indole glucosinolates, derived from the amino acid Trp, are plant secondary metabolites that mediate numerous biological interactions between cruciferous plants and their natural enemies, such as herbivorous insects, pathogens, and other pests. While the genes and enzymes involved in the Arabidopsis thaliana core biosynthetic pathway, leading to indol-3-yl-methyl glucosinolate (I3M), have been identified and characterized, the genes and gene products responsible for modification reactions of the indole ring are largely unknown. Here, we combine the analysis of Arabidopsis mutant lines with a bioengineering approach to clarify which genes are involved in the remaining biosynthetic steps in indole glucosinolate modification. We engineered the indole glucosinolate biosynthesis pathway into Nicotiana benthamiana, showing that it is possible to produce indole glucosinolates in a noncruciferous plant. Building upon this setup, we demonstrate that all members of a small gene subfamily of cytochrome P450 monooxygenases, CYP81Fs, are capable of carrying out hydroxylation reactions of the glucosinolate indole ring, leading from I3M to 4-hydroxy-indol-3-yl-methyl and/or 1-hydroxy-indol-3-yl-methyl glucosinolate intermediates, and that these hydroxy intermediates are converted to 4-methoxy-indol-3-yl-methyl and 1-methoxy-indol-3-yl-methyl glucosinolates by either of two family 2 O-methyltransferases, termed indole glucosinolate methyltransferase 1 (IGMT1) and IGMT2.

The gene controlling the indole glucosinolate modifier1 quantitative trait locus alters indole glucosinolate structures and aphid resistance in Arabidopsis.

Glucosinolates are defensive secondary compounds that display large structural diversity in Arabidopsis thaliana and related plants. Much attention has been paid to variation in the biosynthesis of Met-derived aliphatic glucosinolates and its ecological consequences, but little is known about the genes that cause qualitative and quantitative differences in Trp-derived indole glucosinolates. We use a combination of quantitative trait locus (QTL) fine-mapping and microarray-based transcript profiling to identify CYP81F2 (At5g57220), encoding a cytochrome P450 monooxygenase, as the gene underlying Indole Glucosinolate Modifier1 (IGM1), a metabolic QTL for the accumulation of two modified indole glucosinolates, 4-hydroxy-indole-3-yl-methyl and 4-methoxy-indole-3-yl-methyl glucosinolate. We verify CYP81F2 function with two SALK T-DNA insertion lines and show that CYP81F2 catalyzes the conversion of indole-3-yl-methyl to 4-hydroxy-indole-3-yl-methyl glucosinolate. We further show that the IGM1 QTL is largely caused by differences in CYP81F2 expression, which results from a combination of cis- and trans-acting expression QTL different from known regulators of indole glucosinolate biosynthesis. Finally, we elucidate a potential function of CYP81F2 in plant-insect interactions and find that CYP81F2 contributes to defense against the green peach aphid (Myzus persicae) but not to resistance against herbivory by larvae from four lepidopteran species.

Evolution of metal hyperaccumulation required cis-regulatory changes and triplication of HMA4.

Little is known about the types of mutations underlying the evolution of species-specific traits. The metal hyperaccumulator Arabidopsis halleri has the rare ability to colonize heavy-metal-polluted soils, and, as an extremophile sister species of Arabidopsis thaliana, it is a powerful model for research on adaptation. A. halleri naturally accumulates and tolerates leaf concentrations as high as 2.2% zinc and 0.28% cadmium in dry biomass. On the basis of transcriptomics studies, metal hyperaccumulation in A. halleri has been associated with more than 30 candidate genes that are expressed at higher levels in A. halleri than in A. thaliana. Some of these genes have been genetically mapped to broad chromosomal segments of between 4 and 24 cM co-segregating with Zn and Cd hypertolerance. However, the in planta loss-of-function approaches required to demonstrate the contribution of a given candidate gene to metal hyperaccumulation or hypertolerance have not been pursued to date. Using RNA interference to downregulate HMA4 (HEAVY METAL ATPASE 4) expression, we show here that Zn hyperaccumulation and full hypertolerance to Cd and Zn in A. halleri depend on the metal pump HMA4. Contrary to a postulated global trans regulatory factor governing high expression of numerous metal hyperaccumulation genes, we demonstrate that enhanced expression of HMA4 in A. halleri is attributable to a combination of modified cis-regulatory sequences and copy number expansion, in comparison to A. thaliana. Transfer of an A. halleri HMA4 gene to A. thaliana recapitulates Zn partitioning into xylem vessels and the constitutive transcriptional upregulation of Zn deficiency response genes characteristic of Zn hyperaccumulators. Our results demonstrate the importance of cis-regulatory mutations and gene copy number expansion in the evolution of a complex naturally selected extreme trait. The elucidation of a natural strategy for metal hyperaccumulation enables the rational design of technologies for the clean-up of metal-contaminated soils and for bio-fortification.

Teaching (and learning from) metabolomics: the 2006 PlantMetaNet ETNA Metabolomics Research School.

Under the auspices of the European Training and Networking Activity programme of the European Union, a 'Metabolic Profiling and Data Analysis' Plant Genomics and Bioinformatics Summer School was hosted in Potsdam, Germany between 20 and 29 September 2006. Sixteen early career researchers were invited from the European Union partner nations and the so-called developing nations (Appendix). Lectures from invited leading European researchers provided an overview of the state of the art of these fields and seeded discussion regarding major challenges for their future advancement. Hands-on experience was provided by an example experiment - that of defining the metabolic response of Arabidopsis to treatment of a commercial herbicide of defined mode of action. This experiment was performed throughout the duration of the course in order to teach the concepts underlying extraction and machine handling as well as to provide a rich data set with which the required computation and statistical skills could be illustrated. Here we review the state of the field by describing both key lectures given at and practical aspects taught at the summer school. In addition, we disclose results that were obtained using the four distinct technical platforms at the different participating institutes. While the effects of the chosen herbicide are well documented, this study looks at a broader number of metabolites than in previous investigations. This allowed, on the one hand, not only to characterise further effects of the herbicide than previously observed but also to detect molecules other than the herbicide that were obviously present in the commercial formulation. These data and the workshop in general are all discussed in the context of the teaching of metabolomics.

Different transcript patterns in response to specialist and generalist herbivores in the wild Arabidopsis relative Boechera divaricarpa.

BACKGROUND: Plants defend themselves against herbivorous insects, utilizing both constitutive and inducible defenses. Induced defenses are controlled by several phytohormone-mediated signaling pathways. Here, we analyze transcriptional changes in the North American Arabidopsis relative Boechera divaricarpa in response to larval herbivory by the crucifer specialist lepidopteran Plutella xylostella (diamondback moth) and by the generalist lepidopteran Trichoplusia ni (cabbage semilooper), and compare them to wounding and exogenous phytohormone application. METHODOLOGY/PRINCIPAL FINDINGS: We use a custom macroarray constructed from B. divaricarpa herbivory-regulated cDNAs identified by suppression subtractive hybridization and from known stress-responsive A. thaliana genes for transcript profiling after insect herbivory, wounding and in response to jasmonate, salicylate and ethylene. In addition, we introduce path analysis as a novel approach to analyze transcript profiles. Path analyses reveal that transcriptional responses to the crucifer specialist P. xylostella are primarily determined by direct effects of the ethylene and salicylate pathways, whereas responses to the generalist T. ni are influenced by the ethylene and jasmonate pathways. Wound-induced transcriptional changes are influenced by all three pathways, with jasmonate having the strongest effect. CONCLUSIONS/SIGNIFICANCE: Our results show that insect herbivory is distinct from simple mechanical plant damage, and that different lepidopteran herbivores elicit different transcriptional responses.

Mapping of QTL for resistance against the crucifer specialist herbivore Pieris brassicae in a new Arabidopsis inbred line population, Da(1)-12 x Ei-2.

BACKGROUND: In Arabidopsis thaliana and other crucifers, the glucosinolate-myrosinase system contributes to resistance against herbivory by generalist insects. As yet, it is unclear how crucifers defend themselves against crucifer-specialist insect herbivores. METHODOLOGY/PRINCIPAL FINDINGS: We analyzed natural variation for resistance against two crucifer specialist lepidopteran herbivores, Pieris brassicae and Plutella xylostella, among Arabidopsis thaliana accessions and in a new Arabidopsis recombinant inbred line (RIL) population generated from the parental accessions Da(1)-12 and Ei-2. This RIL population consists of 201 individual F(8) lines genotyped with 84 PCR-based markers. We identified six QTL for resistance against Pieris herbivory, but found only one weak QTL for Plutella resistance. To elucidate potential factors causing these resistance QTL, we investigated leaf hair (trichome) density, glucosinolates and myrosinase activity, traits known to influence herbivory by generalist insects. We identified several previously unknown QTL for these traits, some of which display a complex pattern of epistatic interactions. CONCLUSIONS/SIGNIFICANCE: Although some trichome, glucosinolate or myrosinase QTL co-localize with Pieris QTL, none of these traits explained the resistance QTL convincingly, indicating that resistance against specialist insect herbivores is influenced by other traits than resistance against generalists.

Positive selection driving diversification in plant secondary metabolism.

In Arabidopsis thaliana and related plants, glucosinolates are a major component in the blend of secondary metabolites and contribute to resistance against herbivorous insects. Methylthioalkylmalate synthases (MAM) encoded at the MAM gene cluster control an early step in the biosynthesis of glucosinolates and, therefore, are central to the diversification of glucosinolate metabolism. We sequenced bacterial artificial chromosomes containing the MAM cluster from several Arabidopsis relatives, conducted enzyme assays with heterologously expressed MAM genes, and analyzed MAM nucleotide variation patterns. Our results show that gene duplication, neofunctionalization, and positive selection provide the mechanism for biochemical adaptation in plant defense. These processes occur repeatedly in the history of the MAM gene family, indicating their fundamental importance for the evolution of plant metabolic diversity both within and among species.

Natural variation in MAM within and between populations of Arabidopsis lyrata determines glucosinolate phenotype.

The genetic variation that underlies the glucosinolate phenotype of Arabidopsis lyrata ssp. petraea was investigated between and within populations. A candidate glucosinolate biosynthetic locus (MAM, containing methylthioalkylmalate synthase genes) was mapped in A. lyrata to a location on linkage group 6 corresponding to the homologous location for MAM in A. thaliana. In A. thaliana MAM is responsible for side chain elongation in aliphatic glucosinolates, and the MAM phenotype can be characterized by the ratios of long- to short-chain glucosinolates. A quantitative trait loci (QTL) analysis of glucosinolate ratios in an A. lyrata interpopulation cross found one QTL at MAM. Additional QTL were identified for total indolic glucosinolates and for the ratio of aliphatic to indolic glucosinolates. MAM was then used as the candidate gene for a within-population cosegregation analysis in a natural A. lyrata population from Germany. Extensive variation in microsatellite markers at MAM was found and this variation cosegregated with the same glucosinolate ratios as in the QTL study. The combined results indicate that both between- and within-population genetic variation in the MAM region determines phenotypic variation in glucosinolate side chains in A. lyrata.

Expression profiling and local adaptation of Boechera holboellii populations for water use efficiency across a naturally occurring water stress gradient.

We studied the physiological basis of local adaptation to drought in Boechera holboellii, a perennial relative of Arabidopsis thaliana, and used cDNA-AFLPs to identify candidate genes showing differential expression in these populations. We compared two populations of B. holboellii from contrasting water environments in a reciprocal transplant experiment, as well as in a laboratory dry-down experiment. We continuously measured the water content of soils using time domain reflectometery (TDR). We compared populations for their water use efficiency (WUE), root/shoot ratios (R:S) and leaf mass per unit area (LMA) in the field and in the laboratory, and identified candidate genes that (i) responded plastically to water stress and (ii) were differentially expressed between the two populations. Genotypes from the drier site had higher WUE, which was attributable to a large reduction in transpirational water loss. The xeric-adapted population also had increased investment in root biomass and greater leaf mass per unit area. Reciprocal transplants in the field had significantly greater survival in their native habitat. In total, 450 cDNA-AFLP fragments showed significant changes between drought and control treatments. Furthermore, some genes showed genotype (population)-specific patterns of up- or down-regulation in response to drought. Three hundred cDNA-AFLP bands were sequenced leading to the identification of cDNAs coding for proteins involved in signal transduction, transcriptional regulation, redox regulation, oxidative stress and pathways involved in stress adaptation. Some of these proteins could contribute a physiological advantage under drought, making them potential targets for natural selection.