Manipulating large-scale Arabidopsis microarray expression data: identifying dominant expression patterns and biological process enrichment

Gene regulatory networks shape developmental plasticity of root cell types under water extremes in rice

Understanding how roots modulate development under varied irrigation or rainfall is crucial for development of climate-resilient crops. We established a toolbox of tagged rice lines to profile translating mRNAs and chromatin accessibility within specific cell populations. We used these to study roots in a range of environments: plates in the lab, controlled greenhouse stress and recovery conditions, and outdoors in a paddy. Integration of chromatin and mRNA data resolves regulatory networks of the following: cycle genes in proliferating cells that attenuate DNA synthesis under submergence; genes involved in auxin signaling, the circadian clock, and small RNA regulation in ground tissue; and suberin biosynthesis, iron transporters, and nitrogen assimilation in endodermal/exodermal cells modulated with water availability. By applying a systems approach, we identify known and candidate driver transcription factors of water-deficit responses and xylem development plasticity. Collectively, this resource will facilitate genetic improvements in root systems for optimal climate resilience.

Translational regulation contributes to the elevated CO2 response in two Solanum species

Understanding the impact of elevated CO2 (eCO2 ) in global agriculture is important given climate change projections. Breeding climate-resilient crops depends on genetic variation within naturally varying populations. The effect of genetic variation in response to eCO2 is poorly understood, especially in crop species. We describe the different ways in which Solanum lycopersicum and its wild relative S. pennellii respond to eCO2 , from cell anatomy, to the transcriptome, and metabolome. We further validate the importance of translational regulation as a potential mechanism for plants to adaptively respond to rising levels of atmospheric CO2 .

Dormancy cycling: translation-related transcripts are the main difference between dormant and non-dormant seeds in the field

Primary seed dormancy is a mechanism that orchestrates the timing of seed germination in order to prevent out-of-season germination. Secondary dormancy can be induced in imbibed seeds when they encounter prolonged unfavourable conditions. Secondary dormancy is not induced during dry storage, and therefore the mechanisms underlying this process have remained largely unexplored. Here, a 2-year seed burial experiment in which dormancy cycling was studied at the physiological and transcriptional level is presented. For these analyses six different Arabidopsis thaliana genotypes were used: Landsberg erecta (Ler) and the dormancy associated DELAY OF GERMINATION (DOG) near-isogenic lines 1, 2, 3, 6 and 22 (NILDOG1, 2, 3, 6 and 22). The germination potential of seeds exhumed from the field showed that these seeds go through dormancy cycling and that the dynamics of this cycling is genotype dependent. RNA-seq analysis revealed large transcriptional changes during dormancy cycling, especially at the time points preceding shifts in dormancy status. Dormancy cycling is driven by soil temperature and the endosperm is important in the perception of the environment. Genes that are upregulated in the low- to non-dormant stages are enriched for genes involved in translation, indicating that the non-dormant seeds are prepared for rapid seed germination.

Morphological and molecular responses of the sea urchin Paracentrotus lividus to highly contaminated marine sediments: The case study of Bagnoli-Coroglio brownfield (Mediterranean Sea)

Marine sediments store complex mixtures of compounds, including heavy metals, organotins and a large array of other contaminants. Sediment quality monitoring, characterization and management are priorities, due to potential impacts of the above compounds on coastal waters and their biota, especially in cases of pollutants released during dredging activities. Harbours and marinas, as well as estuaries and bays, where limited exchanges of water occurr, the accumulation of toxic compounds poses major concerns for human and environmental health. Here we report the effects of highly contaminated sediments from the site of national interest Bagnoli-Coroglio (Tyrrhenian Sea, Western Mediterranean) on the sea urchin Paracentrotus lividus, considered a good model for ecotoxicological studies. Adult sea urchins were reared one month in aquaria in the presence of contaminated sediment that was experimentally subject to different patterns of re-suspension events (mimicking the effect of natural storms occurring in the field), crossed with O₂ enrichment versus natural gas exchanges in the water. The development of embryos deriving from adult urchins exposed to such experimental conditions was followed until the pluteus stage, checking the power of contaminated sediment to induce morphological malformations and its eventual buffering by high oxygenation. Real-Time qPCR analysis revealed that the expression of several genes (among the fifty analyzed, involved in different functional processes) was targeted by contaminated sediments more than those exposed in oxygen-enriched condition. Our findings have biological and ecological relevance in terms of assessing the actual impact on local organisms of chronic environmental contamination by heavy metals and polycyclic aromatic hydrocarbons affecting the Bagnoli-Coroglio area, and of exploring enhanced sediment and water oxygenation as a promising tool to mitigate the effects of contamination in future environmental restoration actions.

ANAC032 regulates root growth through the MYB30 gene regulatory network

Reactive oxygen species (ROS) play important roles as root growth regulators. We previously reported a comprehensive transcriptomic atlas, which we named ROS-map, that revealed ROS-responsible genes in Arabidopsis root tips. By using ROS-map, we have characterised an early ROS response key transcription factor, MYB30, as a regulator of root cell elongation under ROS signals. However, there are other ROS-responsible transcription factors which have the potential to regulate root growth. In the present study, we characterised the function of another early ROS-responsible transcription factor, ANAC032, that was selected from ROS-map. Overexpression of ANAC032 fused with the transcriptional activation domain, VP16, inhibited root growth, especially decreasing cell elongation. By transcriptome analysis, we revealed that ANAC032 regulated many stress-responsible genes in the roots. Intriguingly, ANAC032 upregulated MYB30 and its target genes. The upregulation of MYB30 target genes was completely abolished in the ANAC032-VP16x2 OX and ANAC032 estradiol-inducible line in myb30-2 mutants. Moreover, root growth inhibition was alleviated in ANAC032-OX in myb30-2 mutants. Overall, we characterised an upstream transcription factor, ANAC032, of the MYB30 transcriptional cascade which is a key regulator for root cell elongation under ROS signalling.

Spatial identification of transcripts and biological processes in laser micro-dissected sub-regions of waterlogged corn roots with altered expression of phytoglobin

Over-expression of the corn phytoglobin ZmPgb1.2 increases tolerance to waterlogging, while suppression of ZmPgb1.2 compromises plant growth. To unravel compartment-specific transcriptional changes evoked by ZmPgb1.2 during hypoxia, laser micro-dissected sub-regions from waterlogged roots of WT and ZmPgb1.2 overexpressing [ZmPgb1.2(S)] plants were probed for global transcriptional analysis using next generation RNA sequencing. These sub-regions included compartments within the meristematic, elongation, and maturation zone. Of the 149 genes differentially expressed by the up-regulation of ZmPgb1.2, 78 occurred within the meristematic region and included genes involved in jasmonic acid synthesis and response, ascorbic acid metabolism, and ethylene signalling. The ZmPgb1.2 regulation of these genes, discussed in the context of known functions of Pgbs, was further validated by monitoring their expression in meristematic cells of waterlogged roots suppressing ZmPgb1.2. Of the 27 genes differentially expressed by the over-expression of ZmPgb1.2 in the elongation zone, pyruvate kinase and alcohol dehydrogenase showed an expression pattern correlated to the level of ZmPgb1.2 in the tissue. The transcriptional induction of these two enzymes in hypoxic domains of the elongation zone over-expressing ZmPgb1.2 suggests the activation of the fermentation pathway which might be required to sustain metabolic flux and production of ATP in support of cell elongation.

The Lateral Root Cap Acts as an Auxin Sink that Controls Meristem Size

Plant developmental plasticity relies on the activities of meristems, regions where stem cells continuously produce new cells [1]. The lateral root cap (LRC) is the outermost tissue of the root meristem [1], and it is known to play an important role during root development [2-6]. In particular, it has been shown that mechanical or genetic ablation of LRC cells affect meristem size [7, 8]; however, the molecular mechanisms involved are unknown. Root meristem size and, consequently, root growth depend on the position of the transition zone (TZ), a boundary that separates dividing from differentiating cells [9, 10]. The interaction of two phytohormones, cytokinin and auxin, is fundamental in controlling the position of the TZ [9, 10]. Cytokinin via the ARABIDOPSIS RESPONSE REGULATOR 1 (ARR1) control auxin distribution within the meristem, generating an instructive auxin minimum that positions the TZ [10]. We identify a cytokinin-dependent molecular mechanism that acts in the LRC to control the position of the TZ and meristem size. We show that auxin levels within the LRC cells depends on PIN-FORMED 5 (PIN5), a cytokinin-activated intracellular transporter that pumps auxin from the cytoplasm into the endoplasmic reticulum, and on irreversible auxin conjugation mediated by the IAA-amino synthase GRETCHEN HAGEN 3.17 (GH3.17). By titrating auxin in the LRC, the PIN5 and the GH3.17 genes control auxin levels in the entire root meristem. Overall, our results indicate that the LRC serves as an auxin sink that, under the control of cytokinin, regulates meristem size and root growth.

Adapting INTACT to analyse cell-type-specific transcriptomes and nucleocytoplasmic mRNA dynamics in the Arabidopsis embryo

In the early embryo of vascular plants, the different cell types and stem cells of the seedling are specified as the embryo develops from a zygote towards maturity. How the key steps in cell and tissue specification are instructed by genome-wide transcriptional activity is poorly understood. Progress in defining transcriptional regulation at the genome-wide level in plant embryos has been hampered by difficulties associated with capturing cell-type-specific transcriptomes in this small and inaccessible structure. We recently adapted a two-component genetic nucleus labelling system called INTACT to isolate nuclei from distinct cell types at different stages of Arabidopsis thaliana embryogenesis. We have used these to generate a transcriptomic atlas of embryo development following microarray-based expression profiling. Here, we present a general description of the adapted INTACT procedure, including the two-component labelling system, seed isolation, nuclei preparation and purification, as well as transcriptomic profiling. We also compare nuclear and cellular transcriptomes from the early Arabidopsis embryo to assess nucleocytoplasmic differences and discuss how these differences can be used to infer regulation of gene activity.

Transcriptional regulation of nitrogen-associated metabolism and growth

Nitrogen is an essential macronutrient for plant growth and basic metabolic processes. The application of nitrogen-containing fertilizer increases yield, which has been a substantial factor in the green revolution¹. Ecologically, however, excessive application of fertilizer has disastrous effects such as eutrophication². A better understanding of how plants regulate nitrogen metabolism is critical to increase plant yield and reduce fertilizer overuse. Here we present a transcriptional regulatory network and twenty-one transcription factors that regulate the architecture of root and shoot systems in response to changes in nitrogen availability. Genetic perturbation of a subset of these transcription factors revealed coordinate transcriptional regulation of enzymes involved in nitrogen metabolism. Transcriptional regulators in the network are transcriptionally modified by feedback via genetic perturbation of nitrogen metabolism. The network, genes and gene-regulatory modules identified here will prove critical to increasing agricultural productivity.

MYB30 links ROS signaling, root cell elongation, and plant immune responses

Reactive oxygen species (ROS) are known to be important signal molecules that are involved in biotic and abiotic stress responses as well as in growth regulation. However, the molecular mechanisms by which ROS act as a growth regulator, as well as how ROS-dependent growth regulation relates to its roles in stress responses, are not well understood. We performed a time-course microarray analysis of Arabidopsis root tips upon treatment with hydrogen peroxide, which we named "ROS-map." Using the ROS-map, we identified an MYB transcription factor, MYB30, which showed a strong response to ROS treatment and is the key regulator of a gene network that leads to the hydrogen peroxide-dependent inhibition of root cell elongation. Intriguingly, this network contained multiple genes involved in very-long-chain fatty acid (VLCFA) transport. Finally, we showed that MYB30 is necessary for root growth regulation during defense responses, thus providing a molecular link between these two ROS-associated processes.

Transcriptome dynamics revealed by a gene expression atlas of the early Arabidopsis embryo

During early plant embryogenesis, precursors for all major tissues and stem cells are formed. While several components of the regulatory framework are known, how cell fates are instructed by genome-wide transcriptional activity remains unanswered-in part because of difficulties in capturing transcriptome changes at cellular resolution. Here, we have adapted a two-component transgenic labelling system to purify cell-type-specific nuclear RNA and generate a transcriptome atlas of early Arabidopsis embryo development, with a focus on root stem cell niche formation. We validated the dataset through gene expression analysis, and show that gene activity shifts in a spatio-temporal manner, probably signifying transcriptional reprogramming, to induce developmental processes reflecting cell states and state transitions. This atlas provides the most comprehensive tissue- and cell-specific description of genome-wide gene activity in the early plant embryo, and serves as a valuable resource for understanding the genetic control of early plant development.

Metabolite profiling and associated gene expression reveal two metabolic shifts during the seed-to-seedling transition in Arabidopsis thaliana

Metabolic and transcriptomic correlation analysis identified two distinctive profiles involved in the metabolic preparation for seed germination and seedling establishment, respectively. Transcripts were identified that may control metabolic fluxes. The transition from a quiescent metabolic state (dry seed) to the active state of a vigorous seedling is crucial in the plant's life cycle. We analysed this complex physiological trait by measuring the changes in primary metabolism that occur during the transition in order to determine which metabolic networks are operational. The transition involves several developmental stages from seed germination to seedling establishment, i.e. between imbibition of the mature dry seed and opening of the cotyledons, the final stage of seedling establishment. We hypothesized that the advancement of growth is associated with certain signature metabolite profiles. Metabolite-metabolite correlation analysis underlined two specific profiles which appear to be involved in the metabolic preparation for seed germination and efficient seedling establishment, respectively. Metabolite profiles were also compared to transcript profiles and although transcriptional changes did not always equate to a proportional metabolic response, in depth correlation analysis identified several transcripts that may directly influence the flux through metabolic pathways during the seed-to-seedling transition. This correlation analysis also pinpointed metabolic pathways which are significant for the seed-to-seedling transition, and metabolite contents that appeared to be controlled directly by transcript abundance. This global view of the transcriptional and metabolic changes during the seed-to-seedling transition in Arabidopsis opens up new perspectives for understanding the complex regulatory mechanism underlying this transition.

Global transcriptome and coexpression network analyses reveal cultivar-specific molecular signatures associated with seed development and seed size/weight determination in chickpea

Seed development is an intricate process regulated via a complex transcriptional regulatory network. To understand the molecular mechanisms governing seed development and seed size/weight in chickpea, we performed a comprehensive analysis of transcriptome dynamics during seed development in two cultivars with contrasting seed size/weight (small-seeded, Himchana 1 and large-seeded, JGK 3). Our analysis identified stage-specific expression for a significant proportion (>13%) of the genes in each cultivar. About one half of the total genes exhibited significant differential expression in JGK 3 as compared with Himchana 1. We found that different seed development stages can be delineated by modules of coexpressed genes. A comparative analysis revealed differential developmental stage specificity of some modules between the two cultivars. Furthermore, we constructed transcriptional regulatory networks and identified key components determining seed size/weight. The results suggested that extended period of cell division during embryogenesis and higher level of endoreduplication along with more accumulation of storage compounds during maturation determine large seed size/weight. Further, we identified quantitative trait loci-associated candidate genes harboring single nucleotide polymorphisms in the promoter sequences that differentiate small- and large-seeded chickpea cultivars. The results provide a valuable resource to dissect the role of candidate genes governing seed development and seed size/weight in chickpea.

New inter-correlated genes targeted by diatom-derived polyunsaturated aldehydes in the sea urchin Paracentrotus lividus

The marine environment is continually subjected to the action of stressors (including natural toxins), which represent a constant danger for benthic communities. In the present work using network analysis we identified ten genes on the basis of associated functions (FOXA, FoxG, GFI-1, nodal, JNK, OneCut/Hnf6, TAK1, tcf4, TCF7, VEGF) in the sea urchin Paracentrotus lividus, having key roles in different processes, such as embryonic development and asymmetry, cell fate specification, cell differentiation and morphogenesis, and skeletogenesis. These genes are correlated with three HUB genes, Foxo, Jun and HIF1A. Real Time qPCR revealed that during sea urchin embryonic development the expression levels of these genes were modulated by three diatom-derived polyunsaturated aldehydes (PUAs), decadienal, heptadienal and octadienal. Our findings show how changes in gene expression levels may be used as an early indicator of stressful conditions in the marine environment. The identification of key genes and the molecular pathways in which they are involved represents a fundamental tool in understanding how marine organisms try to afford protection against toxicants, to avoid deleterious consequences and irreversible damages. The genes identified in this work as targets for PUAs can be considered as possible biomarkers to detect exposure to different environmental pollutants.

SeqEnrich: A tool to predict transcription factor networks from co-expressed Arabidopsis and Brassica napus gene sets

Transcription factors and their associated DNA binding sites are key regulatory elements of cellular differentiation, development, and environmental response. New tools that predict transcriptional regulation of biological processes are valuable to researchers studying both model and emerging-model plant systems. SeqEnrich predicts transcription factor networks from co-expressed Arabidopsis or Brassica napus gene sets. The networks produced by SeqEnrich are supported by existing literature and predicted transcription factor–DNA interactions that can be functionally validated at the laboratory bench. The program functions with gene sets of varying sizes and derived from diverse tissues and environmental treatments. SeqEnrich presents as a powerful predictive framework for the analysis of Arabidopsis and Brassica napus co-expression data, and is designed so that researchers at all levels can easily access and interpret predicted transcriptional circuits. The program outperformed its ancestral program ChipEnrich, and produced detailed transcription factor networks from Arabidopsis and Brassica napus gene expression data. The SeqEnrich program is ideal for generating new hypotheses and distilling biological information from large-scale expression data.

Transcriptome analysis of the Brassica napus-Leptosphaeria maculans pathosystem identifies receptor, signaling and structural genes underlying plant resistance

The hemibiotrophic fungal pathogen Leptosphaeria maculans is the causal agent of blackleg disease in Brassica napus (canola, oilseed rape) and causes significant loss of yield worldwide. While genetic resistance has been used to mitigate the disease by means of traditional breeding strategies, there is little knowledge about the genes that contribute to blackleg resistance. RNA sequencing and a streamlined bioinformatics pipeline identified unique genes and plant defense pathways specific to plant resistance in the B. napus-L. maculans LepR1-AvrLepR1 interaction over time. We complemented our temporal analyses by monitoring gene activity directly at the infection site using laser microdissection coupled to quantitative PCR. Finally, we characterized genes involved in plant resistance to blackleg in the Arabidopsis-L. maculans model pathosystem. Data reveal an accelerated activation of the plant transcriptome in resistant host cotyledons associated with transcripts coding for extracellular receptors and phytohormone signaling molecules. Functional characterization provides direct support for transcriptome data and positively identifies resistance regulators in the Brassicaceae. Spatial gradients of gene activity were identified in response to L. maculans proximal to the site of infection. This dataset provides unprecedented spatial and temporal resolution of the genes required for blackleg resistance and serves as a valuable resource for those interested in host-pathogen interactions.

Tissue-specific laser microdissection of the Brassica napus funiculus improves gene discovery and spatial identification of biological processes

The three primary tissue systems of the funiculus each undergo unique developmental programs to support the growth and development of the filial seed. To understand the underlying transcriptional mechanisms that orchestrate development of the funiculus at the globular embryonic stage of seed development, we used laser microdissection coupled with RNA-sequencing to produce a high-resolution dataset of the mRNAs present in the epidermis, cortex, and vasculature of the Brassica napus (canola) funiculus. We identified 7761 additional genes in these tissues compared with the whole funiculus organ alone using this technology. Differential expression and enrichment analyses were used to identify several biological processes associated with each tissue system. Our data show that cell wall modification and lipid metabolism are prominent in the epidermis, cell growth and modification occur in the cortex, and vascular tissue proliferation and differentiation occur in the central vascular strand. We provide further evidence that each of the three tissue systems of the globular stage funiculus are involved in specific biological processes that all co-ordinate to support seed development. The identification of genes and gene regulators responsible for tissue-specific developmental processes of the canola funiculus now serves as a valuable resource for seed improvement research.

A Predictive Coexpression Network Identifies Novel Genes Controlling the Seed-to-Seedling Phase Transition in Arabidopsis thaliana

The transition from a quiescent dry seed to an actively growing photoautotrophic seedling is a complex and crucial trait for plant propagation. This study provides a detailed description of global gene expression in seven successive developmental stages of seedling establishment in Arabidopsis (Arabidopsis thaliana). Using the transcriptome signature from these developmental stages, we obtained a coexpression gene network that highlights interactions between known regulators of the seed-to-seedling transition and predicts the functions of uncharacterized genes in seedling establishment. The coexpressed gene data sets together with the transcriptional module indicate biological functions related to seedling establishment. Characterization of the homeodomain leucine zipper I transcription factor AtHB13, which is expressed during the seed-to-seedling transition, demonstrated that this gene regulates some of the network nodes and affects late seedling establishment. Knockout mutants for athb13 showed increased primary root length as compared with wild-type (Columbia-0) seedlings, suggesting that this transcription factor is a negative regulator of early root growth, possibly repressing cell division and/or cell elongation or the length of time that cells elongate. The signal transduction pathways present during the early phases of the seed-to-seedling transition anticipate the control of important events for a vigorous seedling, such as root growth. This study demonstrates that a gene coexpression network together with transcriptional modules can provide insights that are not derived from comparative transcript profiling alone.

Transcriptional programs regulated by both LEAFY and APETALA1 at the time of flower formation

Two key regulators of the switch to flower formation and of flower patterning in Arabidopsis are the plant-specific helix-turn-helix transcription factor LEAFY (LFY) and the MADS box transcription factor APETALA1 (AP1). The interactions between these two transcriptional regulators are complex. AP1 is both a direct target of LFY and can act in parallel with LFY. Available genetic and molecular evidence suggests that LFY and AP1 together orchestrate the switch to flower formation and early events during flower morphogenesis by altering transcriptional programs. However, very little is known about target genes regulated by both transcription factors. Here, we performed a meta-analysis of public datasets to identify genes that are likely to be regulated by both LFY and AP1. Our analyses uncovered known and novel direct LFY and AP1 targets with a role in the control of onset of flower formation. It also identified additional families of proteins and regulatory pathways that may be under transcriptional control by both transcription factors. In particular, several of these genes are linked to response to hormones, to transport and to development. Finally, we show that the gibberellin catabolism enzyme ELA1, which was recently shown to be important for the timing of the switch to flower formation, is positively feedback-regulated by AP1. Our study contributes to the elucidation of the regulatory network that leads to formation of a vital plant organ system, the flower.

RNA-seq analysis of the response of the halophyte, Mesembryanthemum crystallinum (ice plant) to high salinity

Understanding the molecular mechanisms that convey salt tolerance in plants is a crucial issue for increasing crop yield. The ice plant (Mesembryanthemum crystallinum) is a halophyte that is capable of growing under high salt conditions. For example, the roots of ice plant seedlings continue to grow in 140 mM NaCl, a salt concentration that completely inhibits Arabidopsis thaliana root growth. Identifying the molecular mechanisms responsible for this high level of salt tolerance in a halophyte has the potential of revealing tolerance mechanisms that have been evolutionarily successful. In the present study, deep sequencing (RNAseq) was used to examine gene expression in ice plant roots treated with various concentrations of NaCl. Sequencing resulted in the identification of 53,516 contigs, 10,818 of which were orthologs of Arabidopsis genes. In addition to the expression analysis, a web-based ice plant database was constructed that allows broad public access to the data. The results obtained from an analysis of the RNAseq data were confirmed by RT-qPCR. Novel patterns of gene expression in response to high salinity within 24 hours were identified in the ice plant when the RNAseq data from the ice plant was compared to gene expression data obtained from Arabidopsis plants exposed to high salt. Although ABA responsive genes and a sodium transporter protein (HKT1), are up-regulated and down-regulated respectively in both Arabidopsis and the ice plant; peroxidase genes exhibit opposite responses. The results of this study provide an important first step towards analyzing environmental tolerance mechanisms in a non-model organism and provide a useful dataset for predicting novel gene functions.

Application of transcriptomics in Chinese herbal medicine studies

Transcriptomics using DNA microarray has become a practical and popular tool for herbal medicine study because of high throughput, sensitivity, accuracy, specificity, and reproducibility. Therefore, this article focuses on the overview of DNA microarray technology and the application of DNA microarray in Chinese herbal medicine study. To understand the number and the objectives of articles utilizing DNA microarray for herbal medicine study, we surveyed 297 frequently used Chinese medicinal herbs listed in Pharmacopoeia Commission of People's Republic of China. We classified these medicinal herbs into 109 families and then applied PudMed search using “microarray” and individual herbal family as keywords. Although thousands of papers applying DNA microarray in Chinese herbal studies have been published since 1998, most of the articles focus on the elucidation of mechanisms of certain biological effects of herbs. Construction of the bioactivity database containing large-scaled gene expression profiles of quality control herbs can be applied in the future to analyze the biological events induced by herbs, predict the therapeutic potential of herbs, evaluate the safety of herbs, and identify the drug candidate of herbs. Moreover, the linkage of systems biology tools, such as functional genomics, transcriptomics, proteomics, metabolomics, pharmacogenomics and toxicogenomics, will become a new translational platform between Western medicine and Chinese herbal medicine.

Discovering functional modules across diverse maize transcriptomes using COB, the Co-expression Browser

Tools that provide improved ability to relate genotype to phenotype have the potential to accelerate breeding for desired traits and to improve our understanding of the molecular variants that underlie phenotypes. The availability of large-scale gene expression profiles in maize provides an opportunity to advance our understanding of complex traits in this agronomically important species. We built co-expression networks based on genome-wide expression data from a variety of maize accessions as well as an atlas of different tissues and developmental stages. We demonstrate that these networks reveal clusters of genes that are enriched for known biological function and contain extensive structure which has yet to be characterized. Furthermore, we found that co-expression networks derived from developmental or tissue atlases as compared to expression variation across diverse accessions capture unique functions. To provide convenient access to these networks, we developed a public, web-based Co-expression Browser (COB), which enables interactive queries of the genome-wide networks. We illustrate the utility of this system through two specific use cases: one in which gene-centric queries are used to provide functional context for previously characterized metabolic pathways, and a second where lists of genes produced by mapping studies are further resolved and validated using co-expression networks.

Predicting transcriptional circuitry underlying seed coat development

Filling, protection, and dispersal of angiosperm seeds are largely dependent on the development of the maternally derived seed coat. The development of the seed coat in plants such as Arabidopsis thaliana and Glycine max (soybean) is regulated by a complex network of genes and gene products responsible for the establishment and identity of this multicellular structure. Recent studies support the hypothesis that the structure, development, and function of the seed coat are under the control of transcriptional regulators that are specified in space and time. Furthermore, these transcriptional regulators can act in combination to orchestrate the expression of large gene sets. We discuss the underlying transcriptional circuits of the seed coat sub-regions through the interrogation of large-scale datasets, and also provide some ideas on how the identification and analysis of these datasets can be further improved in these two model oilseed systems.

Gibberellin acts positively then negatively to control onset of flower formation in Arabidopsis

The switch to reproductive development is biphasic in many plants, a feature important for optimal pollination and yield. We show that dual opposite roles of the phytohormone gibberellin underpin this phenomenon in Arabidopsis. Although gibberellin promotes termination of vegetative development, it inhibits flower formation. To overcome this effect, the transcription factor LEAFY induces expression of a gibberellin catabolism gene; consequently, increased LEAFY activity causes reduced gibberellin levels. This allows accumulation of gibberellin-sensitive DELLA proteins. The DELLA proteins are recruited by SQUAMOSA PROMOTER BINDING PROTEIN-LIKE transcription factors to regulatory regions of the floral commitment gene APETALA1 and promote APETALA1 up-regulation and floral fate synergistically with LEAFY. The two opposing functions of gibberellin may facilitate evolutionary and environmental modulation of plant inflorescence architecture.

Comprehensive developmental profiles of gene activity in regions and subregions of the Arabidopsis seed

Seeds are complex structures that consist of the embryo, endosperm, and seed-coat regions that are of different ontogenetic origins, and each region can be further divided into morphologically distinct subregions. Despite the importance of seeds for food, fiber, and fuel globally, little is known of the cellular processes that characterize each subregion or how these processes are integrated to permit the coordinated development of the seed. We profiled gene activity genome-wide in every organ, tissue, and cell type of Arabidopsis seeds from fertilization through maturity. The resulting mRNA datasets offer the most comprehensive description of gene activity in seeds with high spatial and temporal resolution,providing unique insights into the function of understudied seed regions. Global comparisons of mRNA populations reveal unexpected overlaps in the functional identities of seed subregions. Analyses of coexpressed gene sets suggest that processes that regulate seed size and filling are coordinated across several subregions. Predictions of gene regulatory networks based on the association of transcription factors with enriched DNA sequence motifs upstream of coexpressed genes identify regulators of seed development. These studies emphasize the utility of these data sets as an essential resource for the study of seed biology.

Microgenomics: genome-scale, cell-specific monitoring of multiple gene regulation tiers

The expression of nuclear protein-coding genes is controlled by dynamic mechanisms ranging from DNA methylation, chromatin modification, and gene transcription to mRNA maturation, turnover, and translation and the posttranslational control of protein function. A genome-scale assessment of the spatiotemporal regulation of gene expression is essential for a comprehensive understanding of gene regulatory networks. However, there are major obstacles to the precise evaluation of gene regulation in multicellular plant organs; these include the monitoring of regulatory processes at levels other than steady-state transcript abundance, resolution of gene regulation in individual cells or cell types, and effective assessment of transient gene activity manifested during development or in response to external cues. This review surveys the advantages and applications of microgenomics technologies that enable panoramic quantitation of cell-type-specific expression in plants, focusing on the importance of querying gene activity at multiple steps in the continuum, from histone modification to selective translation.

A map of cell type-specific auxin responses

In plants, changes in local auxin concentrations can trigger a range of developmental processes as distinct tissues respond differently to the same auxin stimulus. However, little is known about how auxin is interpreted by individual cell types. We performed a transcriptomic analysis of responses to auxin within four distinct tissues of the Arabidopsis thaliana root and demonstrate that different cell types show competence for discrete responses. The majority of auxin-responsive genes displayed a spatial bias in their induction or repression. The novel data set was used to examine how auxin influences tissue-specific transcriptional regulation of cell-identity markers. Additionally, the data were used in combination with spatial expression maps of the root to plot a transcriptomic auxin-response gradient across the apical and basal meristem. The readout revealed a strong correlation for thousands of genes between the relative response to auxin and expression along the longitudinal axis of the root. This data set and comparative analysis provide a transcriptome-level spatial breakdown of the response to auxin within an organ where this hormone mediates many aspects of development.

Isolation and analysis of mRNAs from specific cell types of plants by ribosome immunopurification

Multiple ribosomes assemble onto an individual mRNA to form a polyribosome (polysome) complex. The epitope tagging of specific ribosomal proteins can enable the immunopurification of polysomes from crude cell extracts derived from cryopreserved tissue samples. Through expression of the epitope-tagged ribosomal protein in cell-type and regional specific domains of Arabidopsis thaliana and other organisms it is feasible to quantitatively assess the mRNAs that are associated with ribosomes with cell-specific resolution. Here we present detailed methods for development of transgenics that express a FLAG-tagged version of ribosomal protein L18 (RPL18) under the direction of individual promoters with specific domains of expression, the immunopurification of ribosomes, and bioinformatic analyses of the resultant datasets obtained by microarray profiling. This methodology provides researchers with the opportunity to assess rapid changes at the organ, tissue, regional or cell-type specific level of mRNAs that are associated with ribosomes and therefore engaged in translation.

Deciphering energy-associated gene networks operating in the response of Arabidopsis plants to stress and nutritional cues

Plants need to continuously adjust their transcriptome in response to various stresses that lead to inhibition of photosynthesis and the deprivation of cellular energy. This adjustment is triggered in part by a coordinated re-programming of the energy-associated transcriptome to slow down photosynthesis and activate other energy-promoting gene networks. Therefore, understanding the stress-related transcriptional networks of genes belonging to energy-associated pathways is of major importance for engineering stress tolerance. In a bioinformatics approach developed by our group, termed 'gene coordination', we previously divided genes encoding for enzymes and transcription factors in Arabidopsis thaliana into three clusters, displaying altered coordinated transcriptional behaviors in response to multiple biotic and abiotic stresses (Plant Cell, 23, 2011, 1264). Enrichment analysis indicated further that genes controlling energy-associated metabolism operate as a compound network in response to stress. In the present paper, we describe in detail the network association of genes belonging to six central energy-associated pathways in each of these three clusters described in our previous paper. Our results expose extensive stress-associated intra- and inter-pathway interactions between genes from these pathways, indicating that genes encoding proteins involved in energy-associated metabolism are expressed in a highly coordinated manner. We also provide examples showing that this approach can be further utilized to elucidate candidate genes for stress tolerance and functions of isozymes.

The protein expression landscape of the Arabidopsis root

Because proteins are the major functional components of cells, knowledge of their cellular localization is crucial to gaining an understanding of the biology of multicellular organisms. We have generated a protein expression map of the Arabidopsis root providing the identity and cell type-specific localization of nearly 2,000 proteins. Grouping proteins into functional categories revealed unique cellular functions and identified cell type-specific biomarkers. Cellular colocalization provided support for numerous protein-protein interactions. With a binary comparison, we found that RNA and protein expression profiles are weakly correlated. We then performed peak integration at cell type-specific resolution and found an improved correlation with transcriptome data using continuous values. We performed GeLC-MS/MS (in-gel tryptic digestion followed by liquid chromatography-tandem mass spectrometry) proteomic experiments on mutants with ectopic and no root hairs, providing complementary proteomic data. Finally, among our root hair-specific proteins we identified two unique regulators of root hair development.

The protein expression landscape of the Arabidopsis root

Because proteins are the major functional components of cells, knowledge of their cellular localization is crucial to gaining an understanding of the biology of multicellular organisms. We have generated a protein expression map of the Arabidopsis root providing the identity and cell type-specific localization of nearly 2,000 proteins. Grouping proteins into functional categories revealed unique cellular functions and identified cell type-specific biomarkers. Cellular colocalization provided support for numerous protein-protein interactions. With a binary comparison, we found that RNA and protein expression profiles are weakly correlated. We then performed peak integration at cell type-specific resolution and found an improved correlation with transcriptome data using continuous values. We performed GeLC-MS/MS (in-gel tryptic digestion followed by liquid chromatography-tandem mass spectrometry) proteomic experiments on mutants with ectopic and no root hairs, providing complementary proteomic data. Finally, among our root hair-specific proteins we identified two unique regulators of root hair development.

Cell identity regulators link development and stress responses in the Arabidopsis root

Stress responses in plants are tightly coordinated with developmental processes, but interaction of these pathways is poorly understood. We used genome-wide assays at high spatiotemporal resolution to understand the processes that link development and stress in the Arabidopsis root. Our meta-analysis finds little evidence for a universal stress response. However, common stress responses appear to exist with many showing cell type specificity. Common stress responses may be mediated by cell identity regulators because mutations in these genes resulted in altered responses to stress. Evidence for a direct role for cell identity regulators came from genome-wide binding profiling of the key regulator SCARECROW, which showed binding to regulatory regions of stress-responsive genes. Coexpression in response to stress was used to identify genes involved in specific developmental processes. These results reveal surprising linkages between stress and development at cellular resolution, and show the power of multiple genome-wide data sets to elucidate biological processes.

Transcriptomic characterization of a synergistic genetic interaction during carpel margin meristem development in Arabidopsis thaliana

In flowering plants the gynoecium is the female reproductive structure. In Arabidopsis thaliana ovules initiate within the developing gynoecium from meristematic tissue located along the margins of the floral carpels. When fertilized the ovules will develop into seeds. SEUSS (SEU) and AINTEGUMENTA (ANT) encode transcriptional regulators that are critical for the proper formation of ovules from the carpel margin meristem (CMM). The synergistic loss of ovule initiation observed in the seu ant double mutant suggests that SEU and ANT share overlapping functions during CMM development. However the molecular mechanism underlying this synergistic interaction is unknown. Using the ATH1 transcriptomics platform we identified transcripts that were differentially expressed in seu ant double mutant relative to wild type and single mutant gynoecia. In particular we sought to identify transcripts whose expression was dependent on the coordinated activities of the SEU and ANT gene products. Our analysis identifies a diverse set of transcripts that display altered expression in the seu ant double mutant tissues. The analysis of overrepresented Gene Ontology classifications suggests a preponderance of transcriptional regulators including multiple members of the REPRODUCTIVE MERISTEMS (REM) and GROWTH-REGULATING FACTOR (GRF) families are mis-regulated in the seu ant gynoecia. Our in situ hybridization analyses indicate that many of these genes are preferentially expressed within the developing CMM. This study is the first step toward a detailed description of the transcriptional regulatory hierarchies that control the development of the CMM and ovule initiation. Understanding the regulatory hierarchy controlled by SEU and ANT will clarify the molecular mechanism of the functional redundancy of these two genes and illuminate the developmental and molecular events required for CMM development and ovule initiation.

Coordinated gene networks regulating Arabidopsis plant metabolism in response to various stresses and nutritional cues

The expression pattern of any pair of genes may be negatively correlated, positively correlated, or not correlated at all in response to different stresses and even different progression stages of the stress. This makes it difficult to identify such relationships by classical statistical tools such as the Pearson correlation coefficient. Hence, dedicated bioinformatics approaches that are able to identify groups of cues in which there is a positive or negative expression correlation between pairs or groups of genes are called for. We herein introduce and discuss a bioinformatics approach, termed Gene Coordination, that is devoted to the identification of specific or multiple cues in which there is a positive or negative coordination between pairs of genes and can further incorporate additional coordinated genes to form large coordinated gene networks. We demonstrate the utility of this approach by providing a case study in which we were able to discover distinct expression behavior of the energy-associated gene network in response to distinct biotic and abiotic stresses. This bioinformatics approach is suitable to a broad range of studies that compare treatments versus controls, such as effects of various cues, or expression changes between a mutant and the control wild-type genotype.

Differentiation of Arabidopsis guard cells: analysis of the networks incorporating the basic helix-loop-helix transcription factor, FAMA

Nearly all extant land plants possess stomata, the epidermal structures that mediate gas exchange between the plant and the environment. The developmental pathways, cell division patterns, and molecules employed in the generation of these structures are simple examples of processes used in many developmental contexts. One specific module is a set of "master regulator" basic helix-loop-helix transcription factors that regulate individual consecutive steps in stomatal development. Here, we profile transcriptional changes in response to inducible expression of Arabidopsis (Arabidopsis thaliana) FAMA, a basic helix-loop-helix protein whose actions during the final stage in stomatal development regulate both cell division and cell fate. Genes identified by microarray and candidate approaches were then further analyzed to test specific hypothesis about the activity of FAMA, the shape of its regulatory network, and to create a new set of stomata-specific or stomata-enriched reporters.

Spatiotemporal regulation of cell-cycle genes by SHORTROOT links patterning and growth

The development of multicellular organisms relies on the coordinated control of cell divisions leading to proper patterning and growth. The molecular mechanisms underlying pattern formation, particularly the regulation of formative cell divisions, remain poorly understood. In Arabidopsis, formative divisions generating the root ground tissue are controlled by SHORTROOT (SHR) and SCARECROW (SCR). Here we show, using cell-type-specific transcriptional effects of SHR and SCR combined with data from chromatin immunoprecipitation-based microarray experiments, that SHR regulates the spatiotemporal activation of specific genes involved in cell division. Coincident with the onset of a specific formative division, SHR and SCR directly activate a D-type cyclin; furthermore, altering the expression of this cyclin resulted in formative division defects. Our results indicate that proper pattern formation is achieved through transcriptional regulation of specific cell-cycle genes in a cell-type- and developmental-stage-specific context. Taken together, we provide evidence for a direct link between developmental regulators, specific components of the cell-cycle machinery and organ patterning.

Detecting separate time scales in genetic expression data

Background

Biological processes occur on a vast range of time scales, and many of them occur concurrently. As a result, system-wide measurements of gene expression have the potential to capture many of these processes simultaneously. The challenge however, is to separate these processes and time scales in the data. In many cases the number of processes and their time scales is unknown. This issue is particularly relevant to developmental biologists, who are interested in processes such as growth, segmentation and differentiation, which can all take place simultaneously, but on different time scales.

Results

We introduce a flexible and statistically rigorous method for detecting different time scales in time-series gene expression data, by identifying expression patterns that are temporally shifted between replicate datasets. We apply our approach to a Saccharomyces cerevisiae cell-cycle dataset and an Arabidopsis thaliana root developmental dataset. In both datasets our method successfully detects processes operating on several different time scales. Furthermore we show that many of these time scales can be associated with particular biological functions.

Conclusions

The spatiotemporal modules identified by our method suggest the presence of multiple biological processes, acting at distinct time scales in both the Arabidopsis root and yeast. Using similar large-scale expression datasets, the identification of biological processes acting at multiple time scales in many organisms is now possible.

Omics meet networks - using systems approaches to infer regulatory networks in plants

Many genomic-scale datasets in plants have been generated over the last few years. This substantial achievement has led to impressive progress, including some of the most detailed molecular maps in any multicellular organism. Networks and pathways have been reconstructed using transcriptome, genome-wide transcription factor binding, proteome and metabolome data, and subsequently used to infer functional interactions among genes, proteins, and metabolites. However, more sophisticated systems biology approaches are needed to integrate different omics datasets. Ultimately, the integration of diverse and massive datasets into coherent models will improve our understanding of the molecular networks that underlie biological processes.