Maria Katherine Mejia-Guerra

AGRIS: the Arabidopsis Gene Regulatory Information Server, an update

The Arabidopsis Gene Regulatory Information Server (AGRIS; http://arabidopsis.med.ohio-state.edu/) provides a comprehensive resource for gene regulatory studies in the model plant Arabidopsis thaliana. Three interlinked databases, AtTFDB, AtcisDB and AtRegNet, furnish comprehensive and updated information on transcription factors (TFs), predicted and experimentally verified cis-regulatory elements (CREs) and their interactions, respectively. In addition to significant contributions in the identification of the entire set of TF–DNA interactions, which are the key to understand the gene regulatory networks that govern Arabidopsis gene expression, tools recently incorporated into AGRIS include the complete set of words length 5–15 present in the Arabidopsis genome and the integration of AtRegNet with visualization tools, such as the recently developed ReIN application. All the information in AGRIS is publicly available and downloadable upon registration.

Unraveling the KNOTTED1 regulatory network in maize meristems

KNOTTED1 (KN1)-like homeobox (KNOX) transcription factors function in plant meristems, self-renewing structures consisting of stem cells and their immediate daughters. We defined the KN1 cistrome in maize inflorescences and found that KN1 binds to several thousand loci, including 643 genes that are modulated in one or multiple tissues. These KN1 direct targets are strongly enriched for transcription factors (including other homeobox genes) and genes participating in hormonal pathways, most significantly auxin, demonstrating that KN1 plays a key role in orchestrating the upper levels of a hierarchical gene regulatory network that impacts plant meristem identity and function.

A Genome-Wide Regulatory Framework Identifies Maize Pericarp Color1 Controlled Genes

This study combines high-throughput RNA sequencing and chromatin immunoprecipitation and sequencing to examine the effect of Pericarp Color1 (P1) on global gene expression in maize pericarps and silks, finding that, in addition to regulating flavonoid biosynthesis genes, P1 modulates the expression of a much larger gene set involved in primary metabolism and production of other specialized compounds. Pericarp Color1 (P1) encodes an R2R3-MYB transcription factor responsible for the accumulation of insecticidal flavones in maize (Zea mays) silks and red phlobaphene pigments in pericarps and other floral tissues, which makes P1 an important visual marker. Using genome-wide expression analyses (RNA sequencing) in pericarps and silks of plants with contrasting P1 alleles combined with chromatin immunoprecipitation coupled with high-throughput sequencing, we show here that the regulatory functions of P1 are much broader than the activation of genes corresponding to enzymes in a branch of flavonoid biosynthesis. P1 modulates the expression of several thousand genes, and ∼1500 of them were identified as putative direct targets of P1. Among them, we identified F2H1, corresponding to a P450 enzyme that converts naringenin into 2-hydroxynaringenin, a key branch point in the P1-controlled pathway and the first step in the formation of insecticidal C-glycosyl flavones. Unexpectedly, the binding of P1 to gene regulatory regions can result in both gene activation and repression. Our results indicate that P1 is the major regulator for a set of genes involved in flavonoid biosynthesis and a minor modulator of the expression of a much larger gene set that includes genes involved in primary metabolism and production of other specialized compounds.

Emergence of Switch-Like Behavior in a Large Family of Simple Biochemical Networks

Bistability plays a central role in the gene regulatory networks (GRNs) controlling many essential biological functions, including cellular differentiation and cell cycle control. However, establishing the network topologies that can exhibit bistability remains a challenge, in part due to the exceedingly large variety of GRNs that exist for even a small number of components. We begin to address this problem by employing chemical reaction network theory in a comprehensive in silico survey to determine the capacity for bistability of more than 40,000 simple networks that can be formed by two transcription factor-coding genes and their associated proteins (assuming only the most elementary biochemical processes). We find that there exist reaction rate constants leading to bistability in ∼90% of these GRN models, including several circuits that do not contain any of the TF cooperativity commonly associated with bistable systems, and the majority of which could only be identified as bistable through an original subnetwork-based analysis. A topological sorting of the two-gene family of networks based on the presence or absence of biochemical reactions reveals eleven minimal bistable networks (i.e., bistable networks that do not contain within them a smaller bistable subnetwork). The large number of previously unknown bistable network topologies suggests that the capacity for switch-like behavior in GRNs arises with relative ease and is not easily lost through network evolution. To highlight the relevance of the systematic application of CRNT to bistable network identification in real biological systems, we integrated publicly available protein-protein interaction, protein-DNA interaction, and gene expression data from Saccharomyces cerevisiae, and identified several GRNs predicted to behave in a bistable fashion.

From plant gene regulatory grids to network dynamics

The regulation of gene expression is the most basic level at which genotypes encoded in DNA can manifest themselves into observable phenotypes. In eukaryotes, gene regulatory networks (GRNs) describe the regulatory web through which transcription factors and microRNAs tightly regulate the spatial and temporal expression of genes. In yeast, Escherichia coli, and animals the study of GRNs has uncovered many of the network properties responsible for creating complex regulatory behavior such as organism growth, development, and response to environmental stimuli. In plants, the study of GRNs is just starting to gain momentum thanks to new high quality genomes and the development of new tools for GRN mapping. Here, we review the latest advancements in the study of plant GRNs and describe the tools and techniques used to produce them. We also discuss the emerging field of network dynamics and the methods currently being developed to measure network dynamics and function in plants.

The Maize TFome – development of a transcription factor open reading frame collection for functional genomics

Establishing the architecture of the gene regulatory networks (GRNs) responsible for controlling the transcription of all genes in an organism is a natural development that follows elucidation of the genome sequence. Reconstruction of the GRN requires the availability of a series of molecular tools and resources that so far have been limited to a few model organisms. One such resource consists of collections of transcription factor (TF) open reading frames (ORFs) cloned into vectors that facilitate easy expression in plants or microorganisms. In this study, we describe the development of a publicly available maize TF ORF collection (TFome) of 2034 clones corresponding to 2017 unique gene models in recombination-ready vectors that make possible the facile mobilization of the TF sequences into a number of different expression vectors. The collection also includes several hundred co-regulators (CoREGs), which we classified into well-defined families, and for which we propose here a standard nomenclature, as we have previously done for TFs. We describe the strategies employed to overcome the limitations associated with cloning ORFs from a genome that remains incompletely annotated, with a partial full-length cDNA set available, and with many TF/CoREG genes lacking experimental support. In many instances this required the combination of genome-wide expression data with gene synthesis approaches. The strategies developed will be valuable for developing similar resources for other agriculturally important plants. Information on all the clones generated is available through the GRASSIUS knowledgebase (http://grassius.org/).

Core Promoter Plasticity Between Maize Tissues and Genotypes Contrasts with Predominance of Sharp Transcription Initiation Sites.

Genome-wide determination of maize transcription initiation sites reveals the prevalence of sharp promoters and underscores the role of alternate initiation as a determinant of protein function and localization. Core promoters are crucial for gene regulation, providing blueprints for the assembly of transcriptional machinery at transcription start sites (TSSs). Empirically, TSSs define the coordinates of core promoters and other regulatory sequences. Thus, experimental TSS identification provides an essential step in the characterization of promoters and their features. Here, we describe the application of CAGE (cap analysis of gene expression) to identify genome-wide TSSs used in root and shoot tissues of two maize (Zea mays) inbred lines (B73 and Mo17). Our studies indicate that most TSS clusters are sharp in maize, similar to mice, but distinct from Arabidopsis thaliana, Drosophila melanogaster, or zebra fish, in which a majority of genes have broad-shaped TSS clusters. We established that ∼38% of maize promoters are characterized by a broader TATA-motif consensus, and this motif is significantly enriched in genes with sharp TSSs. A noteworthy plasticity in TSS usage between tissues and inbreds was uncovered, with ∼1500 genes showing significantly different dominant TSSs, sometimes affecting protein sequence by providing alternate translation initiation codons. We experimentally characterized instances in which this differential TSS utilization results in protein isoforms with additional domains or targeted to distinct subcellular compartments. These results provide important insights into TSS selection and gene expression in an agronomically important crop.

Large-scale transcriptional profiling of lignified tissues in Tectona grandis

BackgroundCurrently, Tectona grandis is one of the most valuable trees in the world and no transcript dataset related to secondary xylem is available. Considering how important the secondary xylem and sapwood transition from young to mature trees is, little is known about the expression differences between those successional processes and which transcription factors could regulate lignin biosynthesis in this tropical tree. Although MYB transcription factors are one of the largest superfamilies in plants related to secondary metabolism, it has not yet been characterized in teak. These results will open new perspectives for studies of diversity, ecology, breeding and genomic programs aiming to understand deeply the biology of this species.ResultsWe present a widely expressed gene catalog for T. grandis using Illumina technology and the de novo assembly. A total of 462,260 transcripts were obtained, with 1,502 and 931 genes differentially expressed for stem and branch secondary xylem, respectively, during age transition. Analysis of stem and branch secondary xylem indicates substantial similarity in gene ontologies including carbohydrate enzymes, response to stress, protein binding, and allowed us to find transcription factors and heat-shock proteins differentially expressed. TgMYB1 displays a MYB domain and a predicted coiled-coil (CC) domain, while TgMYB2, TgMYB3 and TgMYB4 showed R2R3-MYB domain and grouped with MYBs from several gymnosperms and flowering plants. TgMYB1, TgMYB4 and TgCES presented higher expression in mature secondary xylem, in contrast with TgMYB2, TgHsp1, TgHsp2, TgHsp3, and TgBi whose expression is higher in young lignified tissues. TgMYB3 is expressed at lower level in secondary xylem.ConclusionsExpression patterns of MYB transcription factors and heat-shock proteins in lignified tissues are dissimilar when tree development was evaluated, obtaining more expression of TgMYB1 and TgMYB4 in lignified tissues of 60-year-old trees, and more expression in TgHsp1, TgHsp2, TgHsp3 and TgBi in stem secondary xylem of 12-year-old trees. We are opening a door for further functional characterization by reverse genetics and marker-assisted selection with those genes. Investigation of some of the key regulators of lignin biosynthesis in teak, however, could be a valuable step towards understanding how rigidity of teak wood and extractives content are different from most other woods. The obtained transcriptome data represents new sequences of T. grandis deposited in public databases, representing an unprecedented opportunity to discover several related-genes associated with secondary xylem such as transcription factors and stress-related genes in a tropical tree.

Identification and Characterization of Maize salmon silks Genes Involved in Insecticidal Maysin Biosynthesis

Characterization of the steps leading to maize natural resistance against insect pests. The century-old maize (Zea mays) salmon silks mutation has been linked to the absence of maysin. Maysin is a C-glycosyl flavone that, when present in silks, confers natural resistance to the maize earworm (Helicoverpa zea), which is one of the most damaging pests of maize in America. Previous genetic analyses predicted Pericarp Color1 (P1; R2R3-MYB transcription factor) to be epistatic to the sm mutation. Subsequent studies identified two loci as being capable of conferring salmon silks phenotypes, salmon silks1 (sm1) and sm2. Benefitting from available sm1 and sm2 mapping information and from knowledge of the genes regulated by P1, we describe here the molecular identification of the Sm1 and Sm2 gene products. Sm2 encodes a rhamnosyl transferase (UGT91L1) that uses isoorientin and UDP-rhamnose as substrates and converts them to rhamnosylisoorientin. Sm1 encodes a multidomain UDP-rhamnose synthase (RHS1) that converts UDP-glucose into UDP-l-rhamnose. Here, we demonstrate that RHS1 shows unexpected substrate plasticity in converting the glucose moiety in rhamnosylisoorientin to 4-keto-6-deoxy glucose, resulting in maysin. Both Sm1 and Sm2 are direct targets of P1, as demonstrated by chromatin immunoprecipitation experiments. The molecular characterization of Sm1 and Sm2 described here completes the maysin biosynthetic pathway, providing powerful tools for engineering tolerance to maize earworm in maize and other plants.

Important biological information uncovered in previously unaligned reads from chromatin immunoprecipitation experiments (ChIP-Seq)

Establishing the architecture of gene regulatory networks (GRNs) relies on chromatin immunoprecipitation followed by massively parallel sequencing (ChIP-Seq) methods that provide genome-wide transcription factor binding sites (TFBSs). ChIP-Seq furnishes millions of short reads that, after alignment, describe the genome-wide binding sites of a particular TF. However, in all organisms investigated an average of 40% of reads fail to align to the corresponding genome, with some datasets having as much as 80% of reads failing to align. We describe here the provenance of previously unaligned reads in ChIP-Seq experiments from animals and plants. We show that a substantial portion corresponds to sequences of bacterial and metazoan origin, irrespective of the ChIP-Seq chromatin source. Unforeseen was the finding that 30%–40% of unaligned reads were actually alignable. To validate these observations, we investigated the characteristics of the previously unaligned reads corresponding to TAL1, a human TF involved in lineage specification of hemopoietic cells. We show that, while unmapped ChIP-Seq read datasets contain foreign DNA sequences, additional TFBSs can be identified from the previously unaligned ChIP-Seq reads. Our results indicate that the re-evaluation of previously unaligned reads from ChIP-Seq experiments will significantly contribute to TF target identification and determination of emerging properties of GRNs.