Deciphering herbivory-induced gene-to-metabolite dynamics in Nicotiana attenuata tissues using a multifactorial...

In response to biotic stresses, such as herbivore attack, plants reorganize their transcriptomes and reconfigure their physiologies not only in attacked tissues but throughout the plant. These whole-organismic reconfigurations are coordinated by a poorly understood network of signal transduction cascades. To explore tissue-based interdependencies in the resistance of Nicotiana attenuata to insect attack, we conducted time-series transcriptome and metabolome profiling of herbivory-elicited source leaves and unelicited sink leaves and roots. To probe the multidimensionality of these molecular responses, we designed a novel approach of combining an extended self-organizing maps-based dimensionality reduction method with bootstrap-based nonparametric analysis of variance models to identify the onset and context of signaling and metabolic pathway activations. We illustrate the value of this analysis by revisiting dynamic changes in the expression of regulatory and structural genes of the oxylipin pathway and by studying nonlinearities in gene-metabolite associations involved in the acyclic diterpene glucoside pathway after selectively extracting modules based on their dynamic response patterns. This novel dimensionality reduction approach is broadly applicable to capture the dynamic rewiring of gene and metabolite networks in experimental design with multiple factors.