RESUMO
MOTIVATION: Identification of system-wide causal relationships can contribute to our understanding of long-distance, intercellular signalling in biological organisms. Dynamic transcriptome analysis holds great potential to uncover coordinated biological processes between organs. However, many existing dynamic transcriptome studies are characterized by sparse and often unevenly spaced time points that make the identification of causal relationships across organs analytically challenging. Application of existing statistical models, designed for regular time series with abundant time points, to sparse data may fail to reveal biologically significant, causal relationships. With increasing research interest in biological time series data, there is a need for new statistical methods that are able to determine causality within and between time series data sets. Here, a statistical framework was developed to identify (Granger) causal gene-gene relationships of unevenly spaced, multivariate time series data from two different tissues of Arabidopsis thaliana in response to a nitrogen signal. RESULTS: This work delivers a statistical approach for modelling irregularly sampled bivariate signals which embeds functions from the domain of engineering that allow to adapt the model's dependence structure to the specific sampling time. Using maximum-likelihood to estimate the parameters of this model for each bivariate time series, it is then possible to use bootstrap procedures for small samples (or asymptotics for large samples) in order to test for Granger-Causality. When applied to the A.thaliana data, the proposed approach produced 3078 significant interactions, in which 2012 interactions have root causal genes and 1066 interactions have shoot causal genes. Many of the predicted causal and target genes are known players in local and long-distance nitrogen signalling, including genes encoding transcription factors, hormones and signalling peptides. Of the 1007 total causal genes (either organ), 384 are either known or predicted mobile transcripts, suggesting that the identified causal genes may be directly involved in long-distance nitrogen signalling through intercellular interactions. The model predictions and subsequent network analysis identified nitrogen-responsive genes that can be further tested for their specific roles in long-distance nitrogen signalling. AVAILABILITY AND IMPLEMENTATION: The method was developed with the R statistical software and is made available through the R package 'irg' hosted on the GitHub repository https://github.com/SMAC-Group/irg where also a running example vignette can be found (https://smac-group.github.io/irg/articles/vignette.html). A few signals from the original data set are made available in the package as an example to apply the method and the complete A.thaliana data can be found at: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE97500. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
RESUMO
Plant responses to multiple environmental stimuli must be integrated to enable them to adapt their metabolism and development. Light and nitrogen (N) are two such stimuli whose downstream signaling pathways must be intimately connected to each other to control plant energy status. Here, we describe the functional role of the WRKY1 transcription factor in controlling genome-wide transcriptional reprogramming of Arabidopsis (Arabidopsis thaliana) leaves in response to individual and combined light and N signals. This includes a cross-regulatory network consisting of 724 genes regulated by WRKY1 and involved in both N and light signaling pathways. The loss of WRKY1 gene function has marked effects on the light and N response of genes involved in N uptake and assimilation (primary metabolism) as well as stress response pathways (secondary metabolism). Our results at the transcriptome and at the metabolite analysis level support a model in which WRKY1 enables plants to activate genes involved in the recycling of cellular carbon resources when light is limiting but N is abundant and upregulate amino acid metabolism when both light and N are limiting. In this potential energy conservation mechanism, WRKY1 integrates information about cellular N and light energy resources to trigger changes in plant metabolism.
