Interpreting omics data with pathway enrichment analysis.

Pathway enrichment analysis is indispensable for interpreting omics datasets and generating hypotheses. However, the foundations of enrichment analysis remain elusive to many biologists. Here, we discuss best practices in interpreting different types of omics data using pathway enrichment analysis and highlight the importance of considering intrinsic features of various types of omics data. We further explain major components that influence the outcomes of a pathway enrichment analysis, including defining background sets and choosing reference annotation databases. To improve reproducibility, we describe how to standardize reporting methodological details in publications. This article aims to serve as a primer for biologists to leverage the wealth of omics resources and motivate bioinformatics tool developers to enhance the power of pathway enrichment analysis.

A novel bivalent chromatin associates with rapid induction of camalexin biosynthesis genes in response to a pathogen signal in Arabidopsis.

Temporal dynamics of gene expression underpin responses to internal and environmental stimuli. In eukaryotes, regulation of gene induction includes changing chromatin states at target genes and recruiting the transcriptional machinery that includes transcription factors. As one of the most potent defense compounds in Arabidopsis thaliana, camalexin can be rapidly induced by bacterial and fungal infections. Though several transcription factors controlling camalexin biosynthesis genes have been characterized, how the rapid activation of genes in this pathway upon a pathogen signal is enabled remains unknown. By combining publicly available epigenomic data with in vivo chromatin modification mapping, we found that camalexin biosynthesis genes are marked with two epigenetic modifications with opposite effects on gene expression, trimethylation of lysine 27 of histone 3 (H3K27me3) (repression) and acetylation of lysine 18 of histone 3 (H3K18ac) (activation), to form a previously uncharacterized type of bivalent chromatin. Mutants with reduced H3K27me3 or H3K18ac suggested that both modifications were required to determine the timing of gene expression and metabolite accumulation at an early stage of the stress response. Our study indicates that the H3K27me3-H3K18ac bivalent chromatin, which we name as kairostat, plays an important role in controlling the timely induction of gene expression upon stress stimuli in plants.

In the fight against harmful fungi and bacteria, plants have an arsenal of chemicals at their disposal. For instance, species in the crucifer family ­ which includes mustard, cabbages and the model plant Arabidopsis thaliana ­ can defend themselves with camalexin, a compound produced soon after the plant receives signals from its attacker. What controls this precise timing, however, is still unclear. For the genes that rule the production of camalexin to be 'read', interpreted, and ultimately converted into proteins, their DNA sequences first need to be physically accessible to the cell. This availability is controlled, in part, by adding or removing chemical groups onto histones, the spool-like structures which DNA wraps around. These precisely controlled modifications ultimately help to activate or repress a gene. Sometimes, activating and inhibiting chemical groups can be present in the same location, creating what is known as a bivalent chromatin domain. Zhao et al. investigated whether histone modifications regulate when A. thaliana produces camalexin in response to an attack. A combination of bioinformatics and experimental approaches highlighted two chemical modifications (one repressive, the other activating) which were present on the same histone, forming a previously unknown bivalent chromatin domain. Mutant plants which did not carry these modifications could not produce camalexin at the right time. Further experiments showed that under normal conditions, both histone modifications were present. However, when the plant was under attack, the level of repressive and activating modifications respectively decreased and increased, leading to gene activation. Together, the results by Zhao et al. suggest that both histone modifications are required for camalexin genes to respond appropriately to signals from a harmful agent. A deeper understanding of this new mechanism could, in turn, allow scientists to engineer crops that are better at resisting disease.

Genome-enabled discovery of anthraquinone biosynthesis in Senna tora.

Senna tora is a widely used medicinal plant. Its health benefits have been attributed to the large quantity of anthraquinones, but how they are made in plants remains a mystery. To identify the genes responsible for plant anthraquinone biosynthesis, we reveal the genome sequence of S. tora at the chromosome level with 526 Mb (96%) assembled into 13 chromosomes. Comparison among related plant species shows that a chalcone synthase-like (CHS-L) gene family has lineage-specifically and rapidly expanded in S. tora. Combining genomics, transcriptomics, metabolomics, and biochemistry, we identify a CHS-L gene contributing to the biosynthesis of anthraquinones. The S. tora reference genome will accelerate the discovery of biologically active anthraquinone biosynthesis pathways in medicinal plants.

A Framework for Discovering, Designing, and Testing MicroProteins to Regulate Synthetic Transcriptional Modules.

Transcription factors often form protein complexes and give rise to intricate transcriptional networks. The regulation of transcription factor multimerization plays a key role in the fine-tuning of the underlying transcriptional pathways and can be exploited to modulate synthetic transcriptional modules. A novel regulation of protein complex formation is emerging: microProteins-truncated transcription factors-engage in protein-protein interactions with transcriptional complexes and modulate their transcriptional activity. Here, we outline a strategy for the discovery, design, and test of putative miPs to fine-tune the activity of transcription factors regulating synthetic or natural transcriptional circuits.

microProtein Prediction Program (miP3): A Software for Predicting microProteins and Their Target Transcription Factors.

An emerging concept in transcriptional regulation is that a class of truncated transcription factors (TFs), called microProteins (miPs), engages in protein-protein interactions with TF complexes and provides feedback controls. A handful of miP examples have been described in the literature but the extent of their prevalence is unclear. Here we present an algorithm that predicts miPs and their target TFs from a sequenced genome. The algorithm is called miP prediction program (miP3), which is implemented in Python. The software will help shed light on the prevalence, biological roles, and evolution of miPs. Moreover, miP3 can be used to predict other types of miP-like proteins that may have evolved from other functional classes such as kinases and receptors. The program is freely available and can be applied to any sequenced genome.

Genetic basis and detection of unintended effects in genetically modified crop plants.

In January 2014, an international meeting sponsored by the International Life Sciences Institute/Health and Environmental Sciences Institute and the Canadian Food Inspection Agency titled "Genetic Basis of Unintended Effects in Modified Plants" was held in Ottawa, Canada, bringing together over 75 scientists from academia, government, and the agro-biotech industry. The objectives of the meeting were to explore current knowledge and identify areas requiring further study on unintended effects in plants and to discuss how this information can inform and improve genetically modified (GM) crop risk assessments. The meeting featured presentations on the molecular basis of plant genome variability in general, unintended changes at the molecular and phenotypic levels, and the development and use of hypothesis-driven evaluations of unintended effects in assessing conventional and GM crops. The development and role of emerging "omics" technologies in the assessment of unintended effects was also discussed. Several themes recurred in a number of talks; for example, a common observation was that no system for genetic modification, including conventional methods of plant breeding, is without unintended effects. Another common observation was that "unintended" does not necessarily mean "harmful". This paper summarizes key points from the information presented at the meeting to provide readers with current viewpoints on these topics.

Becoming data-savvy in a big-data world.

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.

A comprehensive analysis of microProteins reveals their potentially widespread mechanism of transcriptional regulation.

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.

Towards revealing the functions of all genes in plants.

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.

Towards understanding how molecular networks evolve in plants.

Residing beneath the phenotypic landscape of a plant are intricate and dynamic networks of genes and proteins. As evolution operates on phenotypes, we expect its forces to shape somehow these underlying molecular networks. In this review, we discuss progress being made to elucidate the nature of these forces and their impact on the composition and structure of molecular networks. We also outline current limitations and open questions facing the broader field of plant network analysis.

Creation of a genome-wide metabolic pathway database for Populus trichocarpa using a new approach for reconstruction and curation of metabolic pathways for plants.

Metabolic networks reconstructed from sequenced genomes or transcriptomes can help visualize and analyze large-scale experimental data, predict metabolic phenotypes, discover enzymes, engineer metabolic pathways, and study metabolic pathway evolution. We developed a general approach for reconstructing metabolic pathway complements of plant genomes. Two new reference databases were created and added to the core of the infrastructure: a comprehensive, all-plant reference pathway database, PlantCyc, and a reference enzyme sequence database, RESD, for annotating metabolic functions of protein sequences. PlantCyc (version 3.0) includes 714 metabolic pathways and 2,619 reactions from over 300 species. RESD (version 1.0) contains 14,187 literature-supported enzyme sequences from across all kingdoms. We used RESD, PlantCyc, and MetaCyc (an all-species reference metabolic pathway database), in conjunction with the pathway prediction software Pathway Tools, to reconstruct a metabolic pathway database, PoplarCyc, from the recently sequenced genome of Populus trichocarpa. PoplarCyc (version 1.0) contains 321 pathways with 1,807 assigned enzymes. Comparing PoplarCyc (version 1.0) with AraCyc (version 6.0, Arabidopsis [Arabidopsis thaliana]) showed comparable numbers of pathways distributed across all domains of metabolism in both databases, except for a higher number of AraCyc pathways in secondary metabolism and a 1.5-fold increase in carbohydrate metabolic enzymes in PoplarCyc. Here, we introduce these new resources and demonstrate the feasibility of using them to identify candidate enzymes for specific pathways and to analyze metabolite profiling data through concrete examples. These resources can be searched by text or BLAST, browsed, and downloaded from our project Web site (http://plantcyc.org).

Use and misuse of the gene ontology annotations.

The Gene Ontology (GO) project is a collaboration among model organism databases to describe gene products from all organisms using a consistent and computable language. GO produces sets of explicitly defined, structured vocabularies that describe biological processes, molecular functions and cellular components of gene products in both a computer- and human-readable manner. Here we describe key aspects of GO, which, when overlooked, can cause erroneous results, and address how these pitfalls can be avoided.

Bioinformatics and its applications in plant biology.

Bioinformatics plays an essential role in today's plant science. As the amount of data grows exponentially, there is a parallel growth in the demand for tools and methods in data management, visualization, integration, analysis, modeling, and prediction. At the same time, many researchers in biology are unfamiliar with available bioinformatics methods, tools, and databases, which could lead to missed opportunities or misinterpretation of the information. In this review, we describe some of the key concepts, methods, software packages, and databases used in bioinformatics, with an emphasis on those relevant to plant science. We also cover some fundamental issues related to biological sequence analyses, transcriptome analyses, computational proteomics, computational metabolomics, bio-ontologies, and biological databases. Finally, we explore a few emerging research topics in bioinformatics.

PubSearch and PubFetch: a simple management system for semiautomated retrieval and annotation of biological information from the literature.

For most systems in biology, a large body of literature exists that describes the complexity of the system based on experimental results. Manual review of this literature to extract targeted information into biological databases is difficult and time consuming. To address this problem, we developed PubSearch and PubFetch, which store literature, keyword, and gene information in a relational database, index the literature with keywords and gene names, and provide a Web user interface for annotating the genes from experimental data found in the associated literature. A set of protocols is provided in this unit for installing, populating, running, and using PubSearch and PubFetch. In addition, we provide support protocols for performing controlled vocabulary annotations. Intended users of PubSearch and PubFetch are database curators and biology researchers interested in tracking the literature and capturing information about genes of interest in a more effective way than with conventional spreadsheets and lab notebooks.

Using the Arabidopsis Information Resource (TAIR) to find information about Arabidopsis genes.

The Arabidopsis Information Resource (TAIR; http://www.arabidopsis.org) is a comprehensive Web resource of Arabidopsis biology for plant scientists. TAIR curates and integrates information about genes, proteins, gene expression, mutant phenotypes, biological materials such as DNA and seed stocks, genetic markers, genetic and physical maps, biochemical pathways, genome organization, images of mutant plants and protein sub-cellular localizations, publications, and the research community Data in TAIR are extensively interconnected and can be accessed through a variety of Web-based search and display tools. This unit primarily focuses on some basic methods for searching, browsing, visualizing, and analyzing information about Arabidopsis genes. Gene expression data from microarrays is a recent addition to the database and methods for accessing these data are also described. Two pattern identification programs are described for mining TAIR's unique Arabidopsis sequence data sets. We also describe how to use AraCyc for mining plant metabolic pathways.