Thaliana Genome Research Articles

Cloning and Functional Analysis of the MADS-box CiMADS9 Gene from Carya illinoinensis

A MADS-box gene, CiMADS9, was cloned from the male flowers of Carya illinoinensis by rapid amplification of cDNA ends. The gene was 1 077bp with a 768bp open reading frame encoding 255 amino acids. Multiple sequence comparisons revealed that CiMADS9 is a typical MIKC-type MADS-box gene with a MADS-box domain and a K semi-conserved region. Phylogenetic analysis indicated that CiMADS9 belongs to the AGL15 group of the MADS-box gene family. Quantitative reverse transcription polymerase chain reaction analysis indicated that the expression levels in reproductive organs (i.e., flowers and young fruits) were considerably higher than in vegetative tissues (i.e., leaves and branches). The highest expression levels were observed in male flowers. An overexpression vector for CiMADS9 was constructed and the gene was inserted into the Arabidopsis thaliana genome. CiMADS9 expression was confirmed in all transgenic lines. Compared with wild-type plants, transgenic A. thaliana plants overexpressing CiMADS9 exhibited delayed flowering and an increased number of leaves.

Antagonistic control of flowering time by functionally specialized poly(A) polymerases in Arabidopsis thaliana.

Polyadenylation is a critical 3'-end processing step during maturation of pre-mRNAs, and the length of the poly(A) tail affects mRNA stability, nuclear export and translation efficiency. The Arabidopsis thaliana genome encodes three canonical nuclear poly(A) polymerase (PAPS) isoforms fulfilling specialized functions, as reflected by their different mutant phenotypes. While PAPS1 affects several processes, such as the immune response, organ growth and male gametophyte development, the roles of PAPS2 and PAPS4 are largely unknown. Here we demonstrate that PAPS2 and PAPS4 promote flowering in a partially redundant manner. The enzymes act antagonistically to PAPS1, which delays the transition to flowering. The opposite flowering-time phenotypes in paps1 and paps2 paps4 mutants are at least partly due to decreased or increased FLC activity, respectively. In contrast to paps2 paps4 mutants, plants with increased PAPS4 activity flower earlier than the wild-type, concomitant with reduced FLC expression. Double mutant analyses suggest that PAPS2 and PAPS4 act independently of the autonomous pathway components FCA, FY and CstF64. The direct polyadenylation targets of the three PAPS isoforms that mediate their effects on flowering time do not include FLC sense mRNA and remain to be identified. Thus, our results uncover a role for canonical PAPS isoforms in flowering-time control, raising the possibility that modulating the balance of the isoform activities could be used to fine tune the transition to flowering.

Arabidopsis STAY-GREEN, Mendel's Green Cotyledon Gene, Encodes Magnesium-Dechelatase.

Pheophytin a is an essential component of oxygenic photosynthetic organisms because the primary charge separation between chlorophyll a and pheophytin a is the first step in the conversion of light energy. In addition, conversion of chlorophyll a to pheophytin a is the first step of chlorophyll degradation. Pheophytin is synthesized by extracting magnesium (Mg) from chlorophyll; the enzyme Mg-dechelatase catalyzes this reaction. In this study, we report that Mendel's green cotyledon gene, STAY-GREEN (SGR), encodes Mg-dechelatase. The Arabidopsis thaliana genome has three SGR genes, SGR1, SGR2, and STAY-GREEN LIKE (SGRL). Recombinant SGR1/2 extracted Mg from chlorophyll a but had very low or no activity against chlorophyllide a; by contrast, SGRL had higher dechelating activity against chlorophyllide a compared with chlorophyll a All SGRs could not extract Mg from chlorophyll b Enzymatic experiments using the photosystem and light-harvesting complexes showed that SGR extracts Mg not only from free chlorophyll but also from chlorophyll in the chlorophyll-protein complexes. Furthermore, most of the chlorophyll and chlorophyll binding proteins disappeared when SGR was transiently expressed by a chemical induction system. Thus, SGR is not only involved in chlorophyll degradation but also contributes to photosystem degradation.

PanTools: representation, storage and exploration of pan-genomic data.

Next-generation sequencing technology is generating a wealth of highly similar genome sequences for many species, paving the way for a transition from single-genome to pan-genome analyses. Accordingly, genomics research is going to switch from reference-centric to pan-genomic approaches. We define the pan-genome as a comprehensive representation of multiple annotated genomes, facilitating analyses on the similarity and divergence of the constituent genomes at the nucleotide, gene and genome structure level. Current pan-genomic approaches do not thoroughly address scalability, functionality and usability. We introduce a generalized De Bruijn graph as a pan-genome representation, as well as an online algorithm to construct it. This representation is stored in a Neo4j graph database, which makes our approach scalable to large eukaryotic genomes. Besides the construction algorithm, our software package, called PanTools, currently provides functionality for annotating pan-genomes, adding sequences, grouping genes, retrieving gene sequences or genomic regions, reconstructing genomes and comparing and querying pan-genomes. We demonstrate the performance of the tool using datasets of 62 E. coli genomes, 93 yeast genomes and 19 Arabidopsis thaliana genomes. The Java implementation of PanTools is publicly available at http://www.bif.wur.nl sandra.smit@wur.nl.

Variation of 45S rDNA intergenic spacers in Arabidopsis thaliana

Approximately seven hundred 45S rRNA genes (rDNA) in the Arabidopsis thaliana genome are organised in two 4 Mbp-long arrays of tandem repeats arranged in head-to-tail fashion separated by an intergenic spacer (IGS). These arrays make up 5 % of the A. thaliana genome. IGS are rapidly evolving sequences and frequent rearrangements inside the rDNA loci have generated considerable interspecific and even intra-individual variability which allows to distinguish among otherwise highly conserved rRNA genes. The IGS has not been comprehensively described despite its potential importance in regulation of rDNA transcription and replication. Here we describe the detailed sequence variation in the complete IGS of A. thaliana WT plants and provide the reference/consensus IGS sequence, as well as genomic DNA analysis. We further investigate mutants dysfunctional in chromatin assembly factor-1 (CAF-1) (fas1 and fas2 mutants), which are known to have a reduced number of rDNA copies, and plant lines with restored CAF-1 function (segregated from a fas1xfas2 genetic background) showing major rDNA rearrangements. The systematic rDNA loss in CAF-1 mutants leads to the decreased variability of the IGS and to the occurrence of distinct IGS variants. We present for the first time a comprehensive and representative set of complete IGS sequences, obtained by conventional cloning and by Pacific Biosciences sequencing. Our data expands the knowledge of the A. thaliana IGS sequence arrangement and variability, which has not been available in full and in detail until now. This is also the first study combining IGS sequencing data with RFLP analysis of genomic DNA.

Genetic architecture and pleiotropy shape costs of Rps2-mediated resistance in Arabidopsis thaliana.

The mounting evidence that R-genes incur large fitness costs raises a question: how can there be a 5-10% fitness reduction for all 149 R-genes in the Arabidopsis thaliana genome? The R-genes tested to date segregate for insertion-deletion (indel) polymorphisms where susceptible alleles are complete deletions. Since costs of resistance are measured as the differential fitness of isolines carrying resistant and susceptible alleles, indels reveal costs that may be masked when susceptible alleles are expressed. Rps2 segregates for two expressed clades of alleles, one resistant and one susceptible. Plants with resistant Rps2 are not less fit than those with a susceptible Rps2 allele in the absence of disease. Instead, all alleles provide a fitness benefit relative to an artificial deletion, due to the role of RPS2 as a negative regulator of defense. Our results highlight the interplay between genomic architecture and the magnitude of costs of resistance.

Identification of a novel MLPK homologous gene MLPKn1 and its expression analysis in Brassica oleracea.

M locus protein kinase, one of the SRK-interacting proteins, is a necessary positive regulator for the self-incompatibility response in Brassica. In B. rapa, MLPK is expressed as two different transcripts, MLPKf1 and MLPKf2, and either isoform can complement the mlpk/mlpk mutation. The AtAPK1B gene has been considered to be the ortholog of BrMLPK, and AtAPK1B has no role in self-incompatibility (SI) response in A. thaliana SRK-SCR plants. Until now, what causes the MLPK and APK1B function difference during SI response in Brassica and A. thaliana SRKb-SCRb plants has remained unknown. Here, in addition to the reported MLPKf1/2, we identified the new MLPKf1 homologous gene MLPKn1 from B. oleracea. BoMLPKn1 and BoMLPKf1 shared nucleotide sequence identity as high as 84.3%, and the most striking difference consisted in two fragment insertions in BoMLPKn1. BoMLPKn1 and BoMLPKf1 had a similar gene structure; both their deduced amino acid sequences contained a typical plant myristoylation consensus sequence and a Ser/Thr protein kinase domain. BoMLPKn1 was widely expressed in petal, sepal, anther, stigma and leaf. Genome-wide survey revealed that the B. oleracea genome contained three MLPK homologous genes: BoMLPKf1/2, BoMLPKn1 and Bol008343n. The B. rapa genome also contained three MLPK homologous genes, BrMLPKf1/2, BraMLPKn1 and Bra040929. Phylogenetic analysis revealed that BoMLPKf1/2 and BrMLPKf1/2 were phylogenetically more distant from AtAPK1A than Bol008343n, Bra040929, BraMLPKn1 and BoMLPKn1, Synteny analysis revealed that the B. oleracea chromosomal region containing BoMLPKn1 displayed high synteny with the A. thaliana chromosomal region containing APK1B, whereas the B. rapa chromosomal region containing BraMLPKn1 showed high synteny with the A. thaliana chromosomal region containing APK1B. Together, these results revealed that BoMLPKn1/BraMLPKn1, and not the formerly reported BoMLPKf1/2 (BrMLPKf1/2), was the orthologous genes of AtAPK1B, and no ortholog of BoMLPKf1/2 (BrMLPKf1/2) was found in the A. thaliana genome. We speculated that Brassica MLPKf1/2 might have emerged after speciation of Brassica and A. thailiana, and that it was recruited to the SRK-triggered SI signaling cascade in Brassica.

Structure of a Berberine Bridge Enzyme-Like Enzyme with an Active Site Specific to the Plant Family Brassicaceae.

Berberine bridge enzyme-like (BBE-like) proteins form a multigene family (pfam 08031), which is present in plants, fungi and bacteria. They adopt the vanillyl alcohol-oxidase fold and predominantly show bi-covalent tethering of the FAD cofactor to a cysteine and histidine residue, respectively. The Arabidopsis thaliana genome was recently shown to contain genes coding for 28 BBE-like proteins, while featuring four distinct active site compositions. We determined the structure of a member of the AtBBE-like protein family (termed AtBBE-like 28), which has an active site composition that has not been structurally and biochemically characterized thus far. The most salient and distinguishing features of the active site found in AtBBE-like 28 are a mono-covalent linkage of a histidine to the 8α-position of the flavin-isoalloxazine ring and the lack of a second covalent linkage to the 6-position, owing to the replacement of a cysteine with a histidine. In addition, the structure reveals the interaction of a glutamic acid (Glu426) with an aspartic acid (Asp369) at the active site, which appear to share a proton. This arrangement leads to the delocalization of a negative charge at the active site that may be exploited for catalysis. The structure also indicates a shift of the position of the isoalloxazine ring in comparison to other members of the BBE-like family. The dioxygen surrogate chloride was found near the C(4a) position of the isoalloxazine ring in the oxygen pocket, pointing to a rapid reoxidation of reduced enzyme by dioxygen. A T-DNA insertional mutant line for AtBBE-like 28 results in a phenotype, that is characterized by reduced biomass and lower salt stress tolerance. Multiple sequence analysis showed that the active site composition found in AtBBE-like 28 is only present in the Brassicaceae, suggesting that it plays a specific role in the metabolism of this plant family.

ENDOSOMAL RAB EFFECTOR WITH PX-DOMAIN, an Interacting Partner of RAB5 GTPases, Regulates Membrane Trafficking to Protein Storage Vacuoles in Arabidopsis.

RAB5 GTPases act as molecular switches that regulate various endosomal functions in animal cells, including homotypic fusion of early endosomes, endosomal motility, endosomal signaling, and subcompartmentalization of the endosomal membrane. RAB5 proteins fulfill these diverse functions through interactions with downstream effector molecules. Two canonical RAB5 members, ARA7 and RAB HOMOLOG1 (RHA1), are encoded in the Arabidopsis thaliana genome. ARA7 and RHA1 play crucial roles in endocytic and vacuolar trafficking pathways. Plant RAB5 GTPases function via interactions with effector molecules, whose identities and functions are currently unclear. In this study, we searched for canonical RAB5 effector molecules of Arabidopsis and identified a candidate, which we called ENDOSOMAL RAB EFFECTOR WITH PX-DOMAIN (EREX). The intimate genetic interaction between EREX and RAB5 members, the results from subcellular colocalization experiments, and the direct interaction observed in an in vitro pull-down assay strongly suggest that EREX is a genuine effector of canonical RAB5s in Arabidopsis. We further found that close homologs of EREX play partially redundant functions with EREX in the transport of seed storage proteins. Our results indicate that canonical plant RAB5s acquired distinct effector molecules from those of non-plant systems to fulfill their functions.

Quantitative Trait Loci for Morphological Traits and their Association with Functional Genes in Raphanus sativus.

Identification of quantitative trait loci (QTLs) governing morphologically important traits enables to comprehend their potential genetic mechanisms in the genetic breeding program. In this study, we used 210 F2 populations derived from a cross between two radish inbred lines (Raphanus sativus) “835” and “B2,” including 258 SSR markers were used to detect QTLs for 11 morphological traits that related to whole plant, leaf, and root yield in 3 years of replicated field test. Total 55 QTLs were detected which were distributed on each linkage group of the Raphanus genome. Individual QTLs accounted for 2.69–12.6 of the LOD value, and 0.82–16.25% of phenotypic variation. Several genomic regions have multiple traits that clustered together, suggested the existence of pleiotropy linkage. Synteny analysis of the QTL regions with A. thaliana genome selected orthologous genes in radish. InDels and SNPs in the parental lines were detected in those regions by Illumina genome sequence. Five identified candidate gene-based markers were validated by co-mapping with underlying QTLs affecting different traits. Semi-quantitative reverse transcriptase PCR analysis showed the different expression levels of these five genes in parental lines. In addition, comparative QTL analysis with B. rapa revealed six common QTL regions and four key major evolutionarily conserved crucifer blocks (J, U, R, and W) harboring QTL for morphological traits. The QTL positions identified in this study will provide a valuable resource for identifying more functional genes when whole radish genome sequence is released. Candidate genes identified in this study that co-localized in QTL regions are expected to facilitate in radish breeding programs.

Tyrosine Phosphorylation Based Homo-dimerization of Arabidopsis RACK1A Proteins Regulates Oxidative Stress Signaling Pathways in Yeast.

Scaffold proteins are known as important cellular regulators that can interact with multiple proteins to modulate diverse signal transduction pathways. RACK1 (Receptor for Activated C Kinase 1) is a WD-40 type scaffold protein, conserved in eukaryotes, from Chlamydymonas to plants and humans, plays regulatory roles in diverse signal transduction and stress response pathways. RACK1 in humans has been implicated in myriads of neuropathological diseases including Alzheimer and alcohol addictions. Model plant Arabidopsis thaliana genome maintains three different RACK1 genes termed RACK1A, RACK1B, and RACK1C with a very high (85–93%) sequence identity among them. Loss of function mutation in Arabidopsis indicates that RACK1 proteins regulate diverse environmental stress signaling pathways including drought and salt stress resistance pathway. Recently deduced crystal structure of Arabidopsis RACK1A- very first among all of the RACK1 proteins, indicates that it can potentially be regulated by post-translational modifications, like tyrosine phosphorylations and sumoylation at key residues. Here we show evidence that RACK1A proteins, depending on diverse environmental stresses, are tyrosine phosphorylated. Utilizing site-directed mutagenesis of key tyrosine residues, it is found that tyrosine phosphorylation can potentially dictate the homo-dimerization of RACK1A proteins. The homo-dimerized RACK1A proteins play a role in providing UV-B induced oxidative stress resistance. It is proposed that RACK1A proteins ability to function as scaffold protein may potentially be regulated by the homo-dimerized RACK1A proteins to mediate diverse stress signaling pathways.

DNA Gyrase Is the Target for the Quinolone Drug Ciprofloxacin in Arabidopsis thaliana

The Arabidopsis thaliana genome contains four genes that were originally annotated as potentially encoding DNA gyrase: ATGYRA, ATGYRB1, ATGYRB2, and ATGYRB3. Although we subsequently showed that ATGYRB3 does not encode a gyrase subunit, the other three genes potentially encode subunits of a plant gyrase. We also showed evidence for the existence of supercoiling activity in A. thaliana and that the plant is sensitive to quinolone and aminocoumarin antibiotics, compounds that target DNA gyrase in bacteria. However, it was not possible at that time to show whether the A. thaliana genes encoded an active gyrase enzyme, nor whether that enzyme is indeed the target for the quinolone and aminocoumarin antibiotics. Here we show that an A. thaliana mutant resistant to the quinolone drug ciprofloxacin has a point mutation in ATGYRA. Moreover we show that, as in bacteria, the quinolone-sensitive (wild-type) allele is dominant to the resistant gene. Further we have heterologously expressed ATGYRA and ATGYRB2 in a baculovirus expression system and shown supercoiling activity of the partially purified enzyme. Expression/purification of the quinolone-resistant A. thaliana gyrase yields active enzyme that is resistant to ciprofloxacin. Taken together these experiments now show unequivocally that A. thaliana encodes an organelle-targeted DNA gyrase that is the target of the quinolone drug ciprofloxacin; this has important consequences for plant physiology and the development of herbicides.

Elucidation of an Array of Nitrate Transporter Paralogs in Arabidopsis Thaliana Genome

Plants derive gene paralogs after both small-scale duplication (SSD) and whole genome duplication (WGD). Plants are known to retain duplicated genes and enhance their gene repertoire which intensively contributes to the evolution of their genomes. The general view is that such paralogs are both selected and preserved in evolution because they express variable levels of proteins in different temporal and spatial patterns. In Arabidopsis thaliana, paralogs involved in signalling and transcriptional regulation are very often retained after gene and genome duplication events. We analysed duplication events of a primitive nitrate transporter gene NTR across the genome of Arabidopsis thaliana. Majority of the paralogs have acquired the function of both high, low and also dual affinity nitrate transporters, some being regulated by nitrate levels themselves and involved in a variety of plant functions such as embryonic development, repression of lateral root initiation, stomatal opening, xylem and phloem transport and wound response etc. For the present analysis, 19 nitrate transporter protein sequences were retrieved from TAIR database and analysed by a web service tool PINDA as well as manually which gave analogous results. Several transporters were found to be acting as paralogs as close sequence homology was observed amongst them. The results point towards the idea that sub functionalization of nitrate transporter proteins in higher plants happened as a result of a series of gene duplication events.

Arabidopsis Rab Geranylgeranyltransferases Demonstrate Redundancy and Broad Substrate Specificity in Vitro

Posttranslational lipid modifications mediate the membrane attachment of Rab GTPases, facilitating their function in regulating intracellular vesicular trafficking. In Arabidopsis, most Rab GTPases have two C-terminal cysteines and potentially can be double-geranylgeranylated by heterodimeric Rab geranylgeranyltransferases (Rab-GGTs). Genes encoding two putative α subunits and two putative β subunits of Rab-GGTs have been annotated in the Arabidopsis thaliana genome, but little is known about Rab-GGT activity in Arabidopsis. In this study, we demonstrate that four different heterodimers can be formed between putative Arabidopsis Rab-GGT α subunits RGTA1/RGTA2 and β subunits RGTB1/RGTB2, but only RGTA1·RGTB1 and RGTA1·RGTB2 exhibit bona fide Rab-GGT activity, and they are biochemically redundant in vitro. We hypothesize that RGTA2 function might be disrupted by a 12-amino acid insertion in a conserved motif. We present evidence that Arabidopsis Rab-GGTs may have preference for prenylation of C-terminal cysteines in particular positions. We also demonstrate that Arabidopsis Rab-GGTs can not only prenylate a great variety of Rab GTPases in the presence of Rab escort protein but, unlike Rab-GGT in yeast and mammals, can also prenylate certain non-Rab GTPases independently of Rab escort protein. Our findings may help to explain some of the phenotypes of Arabidopsis protein prenyltransferase mutants.

Studies on nickel uptake in transgenicArabidopsis thalianaintroduced with TgMTP1 gene encoding metal tolerance protein

To enhance phytoremediation, which removes heavy metal from soil, transgenic plants were applied to contaminated soil. We constructed a transformation vector expressing both (T. goesingense metal tolerance protein):HA and TgMTP:GFP genes. Transgenic plants were generated using an Agrobacterium-mediated transformation system that expressed the two vectors. Screening and analysis confirmed the incorporation of foreign genes into the Arabidopsis thaliana genome. Callus was induced in the 116 T3 line. These transgenic plants and calli were used for further analyses on the accumulation of Ni. The 116 T3-line plants and calli from selected lines were resistant to heavy metals and accumulated Ni in their leaves. The expression level of TgMTP RNA was equal in all leaves, but protein stability increased in the leaves with Ni treatment. According to these results, we suggest that -overexpressing plants may be useful for phytoremediation of soil.

Towards understanding peroxisomal phosphoregulation in Arabidopsis thaliana.

This work identifies new protein phosphatases and phosphatase-related proteins targeting peroxisomes, and raises the question of a novel protein import pathway from ER to peroxisomes involving peroxisomal targeting signal type 1 (PTS1) Plant peroxisomes are essential for several processes, for example lipid metabolism, free radical detoxification, development, and stress-related functions. Although research on peroxisomes has been intensified, reversible phosphorylation as a control mechanism in peroxisomes is barely studied. Therefore, it is crucial to identify all peroxisomal proteins involved in phosphoregulation. We here started with protein phosphatases, and searched the Arabidopsis thaliana genome for phosphatase-related proteins with putative peroxisomal targeting signals (PTS). Five potential peroxisomal candidates were detected, from which four were confirmed to target peroxisomes or have a functional PTS. The highly conserved Ser-Ser-Met> was validated for two protein phosphatase 2C (PP2C) family members (POL like phosphatases, PLL2 and PLL3) as a functional peroxisomal targeting signal type 1 (PTS1). Full-length PLL2 and PLL3 fused with a reporter protein targeted peroxisomes in two plant expression systems. A putative protein phosphatase, purple acid phosphatase 7 (PAP7), was found to be dually targeted to ER and peroxisomes and experiments indicated a possible trafficking to peroxisomes via the ER depending on peroxisomal PTS1. In addition, a protein phosphatase 2A regulator (TIP41) was validated to harbor a functional PTS1 (Ser-Lys-Val>), but the full-length protein targeted cytosol and nucleus. Reverse genetics indicated a role for TIP41 in senescence signaling. Mass spectrometry of whole seedlings and isolated peroxisomes was employed, and identified new putative phosphorylated peroxisomal proteins. Previously, only one protein phosphatase, belonging to the phospho-protein phosphatase (PPP) family, was identified as a peroxisomal protein. The present work implies that members of two other main protein phosphatase families, i.e. PP2C and PAP, are also targeting peroxisomes.

Boosting AthaMap Database Content with Data from Protein Binding Microarrays

The AthaMap database generates a map of predicted transcription factor binding sites (TFBS) and small RNA target sites for the whole Arabidopsis thaliana genome. With the advent of protein binding microarrays (PBM), the number of known TFBS for A. thaliana transcription factors (TFs) has increased dramatically. Using 113 positional weight matrices (PWMs) generated from a single PBM study and representing a total number of 68 different TFs, the number of predicted TFBS in AthaMap was now increased by about 3.8 × 10(7) to 4.9 × 10(7). The number of TFs with PWM-predicted TFBS annotated in AthaMap has increased to 126, representing a total of 29 TF families and newly including ARF, AT-Hook, YABBY, LOB/AS2 and SRS. Furthermore, links from all Arabidopsis TFs and genes to the newly established Arabidopsis Information Portal (AIP) have been implemented. With this qualitative and quantitative update, the improved AthaMap increases its value as a resource for the analysis of A. thaliana gene expression regulation at www.athamap.de.

CressInt: A user-friendly web resource for genome-scale exploration of gene regulation in Arabidopsis thaliana

The thale cress Arabidopsis thaliana is a powerful model organism for studying a wide variety of biological processes. Recent advances in sequencing technology have resulted in a wealth of information describing numerous aspects of A. thaliana genome function. However, there is a relative paucity of computational systems for efficiently and effectively using these data to create testable hypotheses. We present CressInt, a user-friendly web resource for exploring gene regulatory mechanisms in A. thaliana on a genomic scale. The CressInt system incorporates a variety of genome-wide data types relevant to gene regulation, including transcription factor (TF) binding site models, ChIP-seq, DNase-seq, eQTLs, and GWAS. We demonstrate the utility of CressInt by showing how the system can be used to (1) Identify TFs binding to the promoter of a gene of interest; (2) identify genetic variants that are likely to impact TF binding based on a ChIP-seq dataset; and (3) identify specific TFs whose binding might be impacted by phenotype-associated variants. CressInt is freely available at http://cressint.cchmc.org.

Genome-Wide Prediction and Validation of Intergenic Enhancers in Arabidopsis Using Open Chromatin Signatures.

Enhancers are important regulators of gene expression in eukaryotes. Enhancers function independently of their distance and orientation to the promoters of target genes. Thus, enhancers have been difficult to identify. Only a few enhancers, especially distant intergenic enhancers, have been identified in plants. We developed an enhancer prediction system based exclusively on the DNase I hypersensitive sites (DHSs) in the Arabidopsis thaliana genome. A set of 10,044 DHSs located in intergenic regions, which are away from any gene promoters, were predicted to be putative enhancers. We examined the functions of 14 predicted enhancers using the β-glucuronidase gene reporter. Ten of the 14 (71%) candidates were validated by the reporter assay. We also designed 10 constructs using intergenic sequences that are not associated with DHSs, and none of these constructs showed enhancer activities in reporter assays. In addition, the tissue specificity of the putative enhancers can be precisely predicted based on DNase I hypersensitivity data sets developed from different plant tissues. These results suggest that the open chromatin signature-based enhancer prediction system developed in Arabidopsis may serve as a universal system for enhancer identification in plants.

Genome-Wide Analysis of PAPS1-Dependent Polyadenylation Identifies Novel Roles for Functionally Specialized Poly(A) Polymerases in Arabidopsis thaliana.

The poly(A) tail at 3’ ends of eukaryotic mRNAs promotes their nuclear export, stability and translational efficiency, and changes in its length can strongly impact gene expression. The Arabidopsis thaliana genome encodes three canonical nuclear poly(A) polymerases, PAPS1, PAPS2 and PAPS4. As shown by their different mutant phenotypes, these three isoforms are functionally specialized, with PAPS1 modifying organ growth and suppressing a constitutive immune response. However, the molecular basis of this specialization is largely unknown. Here, we have estimated poly(A)-tail lengths on a transcriptome-wide scale in wild-type and paps1 mutants. This identified categories of genes as particularly strongly affected in paps1 mutants, including genes encoding ribosomal proteins, cell-division factors and major carbohydrate-metabolic proteins. We experimentally verified two novel functions of PAPS1 in ribosome biogenesis and redox homoeostasis that were predicted based on the analysis of poly(A)-tail length changes in paps1 mutants. When overlaying the PAPS1-dependent effects observed here with coexpression analysis based on independent microarray data, the two clusters of transcripts that are most closely coexpressed with PAPS1 show the strongest change in poly(A)-tail length and transcript abundance in paps1 mutants in our analysis. This suggests that their coexpression reflects at least partly the preferential polyadenylation of these transcripts by PAPS1 versus the other two poly(A)-polymerase isoforms. Thus, transcriptome-wide analysis of poly(A)-tail lengths identifies novel biological functions and likely target transcripts for polyadenylation by PAPS1. Data integration with large-scale co-expression data suggests that changes in the relative activities of the isoforms are used as an endogenous mechanism to co-ordinately modulate plant gene expression.