Physiological responses, developmental programs, and cellular functions rely on complex networks of interactions at different levels and scales. Systems biology brings together high-throughput biochemical, genetic, and molecular approaches to generate omics data that can be analyzed and used in mathematical and computational models toward uncovering these networks on a global scale. Various approaches, including transcriptomics, proteomics, interactomics, and metabolomics, have been employed to obtain these data on the cellular, tissue, organ, and whole-plant level. We summarize progress on gene regulatory, cofunction, protein interaction, and metabolic networks. We also illustrate the main approaches that have been used to obtain these networks, with specific examples from Arabidopsis thaliana, and describe the pros and cons of each approach.

Brassinosteroids (BRs) regulate a wide range of developmental and physiological processes in plants through a receptor-kinase signaling pathway that controls the BZR transcription factors Here, we use transcript profiling and chromatin-immunoprecipitation microarray (ChIP-chip) experiments to identify 953 BR-regulated BZR1 target (BRBT) genes Functional studies of selected BRBTs further demonstrate roles in BR promotion of cell elongation The BRBT genes reveal numerous molecular links between the BR-signaling pathway and downstream components involved in developmental and physiological processes Furthermore, the results reveal extensive crosstalk between BR and other hormonal and light-signaling pathways at multiple levels For example, BZR1 not only controls the expression of many signaling components of other hormonal and light pathways but also coregulates common target genes with light-signaling transcription factors Our results provide a genomic map of steroid hormone actions in plants that reveals a regulatory network that integrates hormonal and light-signaling pathways for plant growth regulation

MetaCyc (MetaCyc.org) is a universal database of metabolic pathways and enzymes from all domains of life. The pathways in MetaCyc are curated from the primary scientific literature, and are experimentally determined small-molecule metabolic pathways. Each reaction in a MetaCyc pathway is annotated with one or more well-characterized enzymes. Because MetaCyc contains only experimentally elucidated knowledge, it provides a uniquely high-quality resource for metabolic pathways and enzymes. BioCyc (BioCyc.org) is a collection of more than 350 organism-specific Pathway/Genome Databases (PGDBs). Each BioCyc PGDB contains the predicted metabolic network of one organism, including metabolic pathways, enzymes, metabolites and reactions predicted by the Pathway Tools software using MetaCyc as a reference database. BioCyc PGDBs also contain predicted operons and predicted pathway hole fillerspredictions of which enzymes may catalyze pathway reactions that have not been assigned to an enzyme. The BioCyc website offers many tools for computational analysis of PGDBs, including comparative analysis and analysis of omics data in a pathway context. The BioCyc PGDBs generated by SRI are offered for adoption by any interested party for the ongoing integration of metabolic and genome-related information about an organism.

The great recent progress made in identifying the molecular parts lists of organisms revealed the paucity of our understanding of what most of the parts do. In this review, we introduce computational and statistical approaches and omics data used for inferring gene function in plants, with an emphasis on network-based inference. We also discuss caveats associated with network-based function predictions such as performance assessment, annotation propagation, the guilt-by-association concept, and the meaning of hubs. Finally, we note the current limitations and possible future directions such as the need for gold standard data from several species, unified access to data and tools, quantitative comparison of data and tool quality, and high-throughput experimental validation platforms for systematic gene function elucidation in plants.

Truncated transcription factor-like proteins called microProteins (miPs) can modulate transcription factor activities, thereby increasing transcriptional regulatory complexity. To understand their prevalence, evolution, and function, we predicted over 400 genes that encode putative miPs from Arabidopsis (Arabidopsis thaliana) using a bioinformatics pipeline and validated two novel miPs involved in flowering time and response to abiotic and biotic stress. We provide an evolutionary perspective for a class of miPs targeting homeodomain transcription factors in plants and metazoans. We identify domain loss as one mechanism of miP evolution and suggest the possible roles of miPs on the evolution of their target transcription factors. Overall, we reveal a prominent layer of transcriptional regulation by miPs, show pervasiveness of such proteins both within and across genomes, and provide a framework for studying their function and evolution.

All plants synthesize basic metabolites needed for survival (primary metabolism), but different taxa produce distinct metabolites that are specialized for specific environmental interactions (specialized metabolism). Because evolutionary pressures on primary and specialized metabolism differ, we investigated differences in the emergence and maintenance of these processes across 16 species encompassing major plant lineages from algae to angiosperms. We found that, relative to their primary metabolic counterparts, genes coding for specialized metabolic functions have proliferated to a much greater degree and by different mechanisms and display lineage-specific patterns of physical clustering within the genome and coexpression. These properties illustrate the differential evolution of specialized metabolism in plants, and collectively they provide unique signatures for the potential discovery of novel specialized metabolic processes.

Cellular membranes act as signaling platforms and control solute transport. Membrane receptors, transporters, and enzymes communicate with intracellular processes through protein-protein interactions. Using a split-ubiquitin yeast two-hybrid screen that covers a test-space of 6.4 x 10(6) pairs, we identified 12,102 membrane/signaling protein interactions from Arabidopsis. Besides confirmation of expected interactions such as heterotrimeric G protein subunit interactions and aquaporin oligomerization, >99% of the interactions were previously unknown. Interactions were confirmed at a rate of 32% in orthogonal in planta split-green flourescent protein interaction assays, which was statistically indistinguishable from the confirmation rate for known interactions collected from literature (38%). Regulatory associations in membrane protein trafficking, turnover, and phosphorylation include regulation of potassium channel activity through abscisic acid signaling, transporter activity by a WNK kinase, and a brassinolide receptor kinase by trafficking-related proteins. These examples underscore the utility of the membrane/signaling protein interaction network for gene discovery and hypothesis generation in plants and other organisms.

Plant biology is becoming a data-driven science. High-throughput technologies generate data quickly from molecular to ecosystem levels. Statistical and computational approaches enable describing and interpreting data quantitatively. We highlight the purpose, common problems, and general principles in data analysis. We use RNA sequencing (RNAseq) analysis to illustrate the rationale behind some of the choices made in statistical data analysis. Finally, we provide a list of free online resources that emphasize intuition behind quantitative data analysis.

Background: Gene Ontology (GO) has been used widely to study functional relationships between genes. The current semantic similarity measures rely only on GO annotations and GO structure. This limits the power of GO- based similarity because of the limited proportion of genes that are annotated to GO in most organisms. Results: We introduce a novel approach called NETSIM (network- based similarity measure) that incorporates information from gene co- function networks in addition to using the GO structure and annotations. Using metabolic reaction maps of yeast, Arabidopsis, and human, we demonstrate that NETSIM can improve the accuracy of GO term similarities. We also demonstrate that NETSIM works well even for genomes with sparser gene annotation data. We applied NETSIM on large Arabidopsis gene families such as cytochrome P450 monooxygenases to group the members functionally and show that this grouping could facilitate functional characterization of genes in these families. Conclusions: Using NETSIM as an example, we demonstrated that the performance of a semantic similarity measure could be significantly improved after incorporating genome- specific information. NETSIM incorporates both GO annotations and gene co- function network data as a priori knowledge in the model. Therefore, functional similarities of GO terms that are not explicitly encoded in GO but are relevant in a taxon- specific manner become measurable when GO annotations are limited. Supplementary information and software are available at http://www.msu.edu/similar to jinchen/NETSIM.