Metabolic analysis of the interaction between plants and herbivores

Insect herbivores by necessity have to deal with a large arsenal of plant defence metabolites. The levels of defence compounds may be increased by insect damage. These induced plant responses may also affect the metabolism and performance of successive insect herbivores. As the chemical nature of induced responses is largely unknown, global metabolomic analyses are a valuable tool to gain more insight into the metabolites possibly involved in such interactions. This study analyzed the interaction between feral cabbage (Brassica oleracea) and small cabbage white caterpillars (Pieris rapae) and how previous attacks to the plant affect the caterpillar metabolism. Because plants may be induced by shoot and root herbivory, we compared shoot and root induction by treating the plants on either plant part with jasmonic acid. Extracts of the plants and the caterpillars were chemically analysed using Ultra Performance Liquid Chromatography/Time of Flight Mass Spectrometry (UPLCT/MS). The study revealed that the levels of three structurally related coumaroylquinic acids were elevated in plants treated on the shoot. The levels of these compounds in plants and caterpillars were highly correlated: these compounds were defined as the ‘metabolic interface’. The role of these metabolites could only be discovered using simultaneous analysis of the plant and caterpillar metabolomes. We conclude that a metabolomics approach is useful in discovering unexpected bioactive compounds involved in ecological interactions between plants and their herbivores and higher trophic levels.

The influence of metabolically engineered glucosinolates profiles in Arabidopsis thaliana on Plutella xylostella preference and performance

The oviposition preference and larval performance of the diamondback moth (DBM), Plutella xylostella, was studied using Arabidopsis thaliana plants with modified glucosinolate (GS) profiles containing novel GSs as a result of the introduction of individual CYP79 genes. The insect parameters were determined in a series of bioassays. The GS content of the plants as well as the number of trichomes were measured. Multivariate analysis was used to determine the possible relationships among insect and plant variables. The novel GSs in the tested lines did not appear to have any unequivocal effect on the DBM. Instead, the plant characteristics that affected larval performance and larval preference did not influence oviposition preference. Trichomes did not affect oviposition, but influenced larval parameters negatively. Although the tested A. thaliana lines had earlier been shown to influence disease resistance, in this study no clear results were found for P. xylostella

Antibacterial activity and mode of action of selected glucosinolate hydrolysis products against bacterial pathogens

Plants contain numerous components that are important sources of new bioactive molecules with antimicrobial properties. Isothiocyanates (ITCs) are plant secondary metabolites found in cruciferous vegetables that are arising as promising antimicrobial agents in food industry. The aim of this study was to assess the antibacterial activity of two isothiocyanates (ITCs), allylisothiocyanate (AITC) and 2-phenylethylisothiocyanate (PEITC) against Escherichia coli, Pseudomonas aeruginosa, Staphylococcus aureus and Listeria monocytogenes. The antibacterial mode of action was also characterized by the assessment of different physiological indices: membrane integrity, intracellular potassium release, physicochemical surface properties and surface charge. The minimum inhibitory concentration (MIC) of AITC and PEITC was 100 g/mL for all bacteria. The minimum bactericidal concentration (MBC) of the ITCs was at least 10 times higher than the MIC. Both AITC and PEITC changed the membrane properties of the bacteria decreasing their surface charge and compromising the integrity of the cytoplasmatic membrane with consequent potassium leakage and propidium iodide uptake. The surface hydrophobicity was also non-specifically altered (E. coli and L. monocytogenes become less hydrophilic; P. aeruginosa and S. aureus become more hydrophilic). This study shows that AITC and PEITC have strong antimicrobial potential against the bacteria tested, through the disruption of the bacterial cell membranes. Moreover, phytochemicals are highlighted as a valuable sustainable source of new bioactive products.This work was supported by the Operational Programme for Competitiveness Factors - COMPETE and by the Portuguese Foundation for Science and Technology through Project Phytodisinfectants - PTDC/DTP-SAP/1078/2012 (COMPETE: FCOMP-01-0124-FEDER-028765), the PhD grant awarded to Ana Abreu (SFRH/BD/84393/2012), and the post-doctoral grants awarded to Anabela Borges (SFRH/BPD/98684/2013) and Lucia C. Simoes (SFRH/BPD/81982/2011). Also, this work was undertaken as part of the European Research Project SUSCLEAN (Contract no FP7-KBBE-2011-5, project number: 287514) and the COST Action FA1202. The authors are solely responsible for this work. It does not represent the opinion of the European Community, and the Community is not responsible for any use that might be made of data appearing herein

Arabidopsis thaliana encodes a bacterial-type heterodimeric isopropylmalate isomerase involved in both Leu biosynthesis and the Met chain elongation pathway of glucosinolate formation

The last steps of the Leu biosynthetic pathway and the Met chain elongation cycle for glucosinolate formation share identical reaction types suggesting a close evolutionary relationship of these pathways. Both pathways involve the condensation of acetyl-CoA and a 2-oxo acid, isomerization of the resulting 2-malate derivative to form a 3-malate derivative, the oxidation-decarboxylation of the 3-malate derivative to give an elongated 2-oxo acid, and transamination to generate the corresponding amino acid. We have now analyzed the genes encoding the isomerization reaction, the second step of this sequence, in Arabidopsis thaliana. One gene encodes the large subunit and three encode small subunits of this enzyme, referred to as isopropylmalate isomerase (IPMI) with respect to the Leu pathway. Metabolic profiling of large subunit mutants revealed accumulation of intermediates of both Leu biosynthesis and Met chain elongation, and an altered composition of aliphatic glucosinolates demonstrating the function of this gene in both pathways. In contrast, the small subunits appear to be specialized to either Leu biosynthesis or Met chain elongation. Green fluorescent protein tagging experiments confirms the import of one of the IPMI small subunits into the chloroplast, the localization of the Met chain elongation pathway in these organelles. These results suggest the presence of different heterodimeric IPMIs in Arabidopsis chloroplasts with distinct substrate specificities for Leu or glucosinolate metabolism determined by the nature of the different small subunit

A Systems Biology Approach Identifies a R2R3 MYB Gene Subfamily with Distinct and Overlapping Functions in Regulation of Aliphatic Glucosinolates

BACKGROUND: Glucosinolates are natural metabolites in the order Brassicales that defend plants against both herbivores and pathogens and can attract specialized insects. Knowledge about the genes controlling glucosinolate regulation is limited. Here, we identify three R2R3 MYB transcription factors regulating aliphatic glucosinolate biosynthesis in Arabidopsis by combining several systems biology tools. METHODOLOGY/PRINCIPAL FINDINGS: MYB28 was identified as a candidate regulator of aliphatic glucosinolates based on its co-localization within a genomic region controlling variation both in aliphatic glucosinolate content (metabolite QTL) and in transcript level for genes involved in the biosynthesis of aliphatic glucosinolates (expression QTL), as well as its co-expression with genes in aliphatic glucosinolate biosynthesis. A phylogenetic analysis with the R2R3 motif of MYB28 showed that it and two homologues, MYB29 and MYB76, were members of an Arabidopsis-specific clade that included three characterized regulators of indole glucosinolates. Over-expression of the individual MYB genes showed that they all had the capacity to increase the production of aliphatic glucosinolates in leaves and seeds and induce gene expression of aliphatic biosynthetic genes within leaves. Analysis of leaves and seeds of single knockout mutants showed that mutants of MYB29 and MYB76 have reductions in only short-chained aliphatic glucosinolates whereas a mutant in MYB28 has reductions in both short- and long-chained aliphatic glucosinolates. Furthermore, analysis of a double knockout in MYB28 and MYB29 identified an emergent property of the system since the absence of aliphatic glucosinolates in these plants could not be predicted by the chemotype of the single knockouts. CONCLUSIONS/SIGNIFICANCE: It seems that these cruciferous-specific MYB regulatory genes have evolved both overlapping and specific regulatory capacities. This provides a unique system within which to study the evolution of MYB regulatory factors and their downstream targets

Tri-trophic effects of inter- and intra-population variation in defence chemistry of wild cabbage (Brassica oleracea)

The effect of direct chemical defences in plants on the performance of insect herbivores and their natural enemies has received increasing attention over the past 10 years. However, much less is known about the scale at which this variation is generated and maintained, both within and across populations of the same plant species. This study compares growth and development of the large cabbage butterfly, Pieris brassicae, and its gregarious pupal parasitoid, Pteromalus puparum, on three wild populations [Kimmeridge (KIM), Old Harry (OH) and Winspit (WIN)] and two cultivars [Stonehead (ST), and Cyrus (CYR)] of cabbage, Brassica oleracea. The wild populations originate from the coast of Dorset, UK, but grow in close proximity with one another. Insect performance and chemical profiles were made from every plant used in the experiment. Foliar glucosinolates (GS) concentrations were highest in the wild plants in rank order WIN > OH > KIM, with lower levels found in the cultivars. Caterpillar-damaged leaves in the wild cabbages also had higher GS levels than undamaged leaves. Pupal mass in P. brassicae varied significantly among populations of B. oleracea. Moreover, development time in the host and parasitoid were correlated, even though these stages are temporally separated. Parasitoid adult dry mass closely approximated the development of its host. Multivariate statistics revealed a correlation between pupal mass and development time of P. brassicae and foliar GS chemistry, of which levels of neoglucobrassicin appeared to be the most important. Our results show that there is considerable variation in quantitative aspects of defensive chemistry in wild cabbage plants that is maintained at very small spatial scales in nature. Moreover, the performance of the herbivore and its parasitoid were both affected by differences in plant quality

Population-Related Variation in Plant Defense more Strongly Affects Survival of an Herbivore than Its Solitary Parasitoid Wasp

The performance of natural enemies, such as parasitoid wasps, is affected by differences in the quality of the host’s diet, frequently mediated by species or population-related differences in plant allelochemistry. Here, we compared survival, development time, and body mass in a generalist herbivore, the cabbage moth, Mamestra brassicae, and its solitary endoparasitoid, Microplitis mediator, when reared on two cultivated (CYR and STH) and three wild (KIM, OH, and WIN) populations of cabbage, Brassica oleracea. Plants either were undamaged or induced by feeding of larvae of the cabbage butterfly, Pieris rapae. Development and biomass of M. brassicae and Mi. mediator were similar on both cultivated and one wild cabbage population (KIM), intermediate on the OH population, and significantly lower on the WIN population. Moreover, development was prolonged and biomass was reduced on herbivore-induced plants. However, only the survival of parasitized hosts (and not that of healthy larvae) was affected by induction. Analysis of glucosinolates in leaves of the cabbages revealed higher levels in the wild populations than cultivars, with the highest concentrations in WIN plants. Multivariate statistics revealed a negative correlation between insect performance and total levels of glucosinolates (GS) and levels of 3-butenyl GS. However, GS chemistry could not explain the reduced performance on induced plants since only indole GS concentrations increased in response to herbivory, which did not affect insect performance based on multivariate statistics. This result suggests that, in addition to aliphatic GS, other non-GS chemicals are responsible for the decline in insect performance, and that these chemicals affect the parasitoid more strongly than the host. Remarkably, when developing on WIN plants, the survival of Mi. mediator to adult eclosion was much higher than in its host, M. brassicae. This may be due to the fact that hosts parasitized by Mi. mediator pass through fewer instars, and host growth is arrested when they are only a fraction of the size of healthy caterpillars. Certain aspects of the biology and life-history of the host and parasitoid may determine their response to chemical challenges imposed by the food plant

Molecular and life-history effects of a natural toxin on herbivorous and non-target soil arthropods

Natural toxins, such as isothiocyanate (ITC), are harmful secondary metabolites produced by plants. Many natural toxins occur in commercial crops, yet their possible negative repercussions on especially non-target soil organisms are largely unknown. This study examined life-history and gene transcriptional responses to 2-phenylethyl ITC on two soil arthropod species: Folsomia candida and Protaphorura fimata. To that end the standardized ISO guideline for ecotoxicological tests and a microarray for F. candida were used. The dissipation of 2-phenylethyl ITC in natural soil was investigated using GC-MS/MS for quantification. Half-lives, tested at four concentration levels in natural soil, were on average 16 h with biodegradation as the plausible main removal process. Regardless, toxic effects on reproduction were shown for F. candida and P. fimata, with EC50 values of around 11.5 nmol/g soil illustrating the toxic character of this compound. Gene expression profiles revealed the importance of fatty acid metabolism at low exposure concentrations (EC10), which is associated with the lipophilic nature of 2-phenylethyl ITC. At higher concentrations (EC50) gene expression became more ubiquitous with over-expression of especially stress-related genes and sugar metabolism. The regulation of a gene encoding a precursor of follistatin, furthermore, implied the inhibition of reproduction and may be an important molecular target that can be linked to the observed adverse effect of life-history traits

What Can Causal Networks Tell Us about Metabolic Pathways?

Graphical models describe the linear correlation structure of data and have been used to establish causal relationships among phenotypes in genetic mapping populations. Data are typically collected at a single point in time. Biological processes on the other hand are often non-linear and display time varying dynamics. The extent to which graphical models can recapitulate the architecture of an underlying biological processes is not well understood. We consider metabolic networks with known stoichiometry to address the fundamental question: “What can causal networks tell us about metabolic pathways?”. Using data from an Arabidopsis BaySha population and simulated data from dynamic models of pathway motifs, we assess our ability to reconstruct metabolic pathways using graphical models. Our results highlight the necessity of non-genetic residual biological variation for reliable inference. Recovery of the ordering within a pathway is possible, but should not be expected. Causal inference is sensitive to subtle patterns in the correlation structure that may be driven by a variety of factors, which may not emphasize the substrate-product relationship. We illustrate the effects of metabolic pathway architecture, epistasis and stochastic variation on correlation structure and graphical model-derived networks. We conclude that graphical models should be interpreted cautiously, especially if the implied causal relationships are to be used in the design of intervention strategies

Mapping of QTL for Resistance against the Crucifer Specialist Herbivore Pieris brassicae in a New Arabidopsis Inbred Line Population, Da(1)-12×Ei-2

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.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.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