An in silico assessment of gene function and organization of the phenylpropanoid pathway metabolic networks in Arabidopsis thaliana and limitations thereof

The Arabidopsis genome sequencing in 2000 gave to science the first blueprint of a vascular plant. Its successful completion also prompted the US National Science Foundation to launch the Arabidopsis 2010 initiative, the goal of which is to identify the function of each gene by 2010. In this study, an exhaustive analysis of The Institute for Genomic Research (TIGR) and The Arabidopsis Information Resource (TAIR) databases, together with all currently compiled EST sequence data, was carried out in order to determine to what extent the various metabolic networks from phenylalanine ammonia lyase (PAL) to the monolignols were organized and/or could be predicted. In these databases, there are some 65 genes which have been annotated as encoding putative enzymatic steps in monolignol biosynthesis, although many of them have only very low homology to monolignol pathway genes of known function in other plant systems. Our detailed analysis revealed that presently only 13 genes (two PALs, a cinnamate-4-hydroxylase, a p-coumarate-3-hydroxylase, a ferulate-5-hydroxylase, three 4-coumarate-CoA ligases, a cinnamic acid O-methyl transferase, two cinnamoyl-CoA reductases) and two cinnamyl alcohol dehydrogenases can be classified as having a bona fide (definitive) function; the remaining 52 genes currently have undetermined physiological roles. The EST database entries for this particular set of genes also provided little new insight into how the monolignol pathway was organized in the different tissues and organs, this being perhaps a consequence of both limitations in how tissue samples were collected and in the incomplete nature of the EST collections. This analysis thus underscores the fact that even with genomic sequencing, presumed to provide the entire suite of putative genes in the monolignol-forming pathway, a very large effort needs to be conducted to establish actual catalytic roles (including enzyme versatility), as well as the physiological function(s) for each member of the (multi)gene families present and the metabolic networks that are operative. Additionally, one key to identifying physiological functions for many of these (and other) unknown genes, and their corresponding metabolic networks, awaits the development of technologies to comprehensively study molecular processes at the single cell level in particular tissues and organs, in order to establish the actual metabolic context.     

Clade classification of monolignol biosynthesis gene family members reveals target genes to decrease lignin in Lolium perenne

In monocots, lignin content has a strong impact on the digestibility of the cell wall fraction. Engineering lignin biosynthesis requires a profound knowledge of the role of paralogues in the multigene families that constitute the monolignol biosynthesis pathway. We applied a bioinformatics approach for genome‐wide identification of candidate genes in Lolium perenne that are likely to be involved in the biosynthesis of monolignols. More specifically, we performed functional subtyping of phylogenetic clades in four multigene families: 4CL, COMT, CAD and CCR. Essential residues were considered for functional clade delineation within these families. This classification was complemented with previously published experimental evidence on gene expression, gene function and enzymatic activity in closely related crops and model species. This allowed us to assign functions to novel identified L. perenne genes, and to assess functional redundancy among paralogues. We found that two 4CL paralogues, two COMT paralogues, three CCR paralogues and one CAD gene are prime targets for genetic studies to engineer developmentally regulated lignin in this species. Based on the delineation of sequence conservation between paralogues and a first analysis of allelic diversity, we discuss possibilities to further study the roles of these paralogues in lignin biosynthesis, including expression analysis, reverse genetics and forward genetics, such as association mapping. We propose criteria to prioritise paralogues within multigene families and certain SNPs within these genes for developing genotyping assays or increasing power in association mapping studies. Although L. perenne was the target of the analyses presented here, this functional subtyping of phylogenetic clades represents a valuable tool for studies investigating monolignol biosynthesis genes in other monocot species.     

Local coexpression domains of two to four genes in the genome of Arabidopsis

Expression of genes in eukaryotic genomes is known to cluster, but cluster size is generally loosely defined and highly variable. We have here taken a very strict definition of cluster as sets of physically adjacent genes that are highly coexpressed and form so-called local coexpression domains. The Arabidopsis (Arabidopsis thaliana) genome was analyzed for the presence of such local coexpression domains to elucidate its functional characteristics. We used expression data sets that cover different experimental conditions, organs, tissues, and cells from the Massively Parallel Signature Sequencing repository and microarray data (Affymetrix) from a detailed root analysis. With these expression data, we identified 689 and 1,481 local coexpression domains, respectively, consisting of two to four genes with a pairwise Pearson's correlation coefficient larger than 0.7. This number is approximately 1- to 5-fold higher than the numbers expected by chance. A small (5%-10%) yet significant fraction of genes in the Arabidopsis genome is therefore organized into local coexpression domains. These local coexpression domains were distributed over the genome. Genes in such local domains were for the major part not categorized in the same functional category (GOslim). Neither tandemly duplicated genes nor shared promoter sequence nor gene distance explained the occurrence of coexpression of genes in such chromosomal domains. This indicates that other parameters in genes or gene positions are important to establish coexpression in local domains of Arabidopsis chromosomes.     

Identification of a novel cis-acting element conferring sulfur deficiency response in Arabidopsis roots

SULTR1;1 high-affinity sulfate transporter is highly regulated in the epidermis and cortex of Arabidopsis roots responding to sulfur deficiency (-S). We identified a novel cis-acting element involved in the -S-inducible expression of sulfur-responsive genes in Arabidopsis. The promoter region of SULTR1;1 was dissected for deletion and gain-of-function analysis using luciferase (LUC) reporter gene in transgenic Arabidopsis. The 16-bp sulfur-responsive element (SURE) from -2777 to -2762 of SULTR1;1 promoter was sufficient and necessary for the -S-responsive expression, which was reversed when supplied with cysteine and glutathione (GSH). The SURE sequence contained an auxin response factor (ARF) binding sequence (GAGACA). However, SURE was not responsive to naphthalene acetic acid, indicating its specific function in the sulfur response. The base substitution analysis indicated the significance of a 5-bp sequence (GAGAC) within the conserved ARF binding site as a core element for the -S response. Microarray analysis of early -S response in Arabidopsis roots indicated the presence of SURE core sequences in the promoter regions of -S-inducible genes on a full genome GeneChip array. It is suggested that SURE core sequences may commonly regulate the expression of a gene set required for adaptation to the -S environment.     

Promoter analysis of iron-deficiency-inducible barley IDS3 gene in Arabidopsis and tobacco plants.

Under conditions of iron deficiency, graminaceous plants induce the expression of genes involved in the biosynthesis of mugineic acid family phytosiderophores. We previously identified the novel cis-acting elements IDE1 and IDE2 (iron-deficiency-responsive element 1 and 2) through promoter analysis of the barley (Hordeum vulgare L.) iron-deficiency-inducible IDS2 gene in tobacco (Nicotiana tabacum L.). To gain further insight into plant gene regulation under iron deficiency, we analyzed the barley iron-deficiency-inducible IDS3 gene, which encodes mugineic acid synthase. IDS3 promoter fragments were fused to the beta-glucuronidase (GUS) gene, and this construct was introduced into Arabidopsis thaliana L. and tobacco plants. In both Arabidopsis and tobacco, GUS activity driven by the IDS3 promoter showed strongly iron-deficiency-inducible and root-specific expression. Expression occurred mainly in the epidermis of Arabidopsis roots, whereas expression was dominant in the pericycle, endodermis, and cortex of tobacco roots, resembling the expression pattern conferred by IDE1 and IDE2. Deletion analysis revealed that a sequence within -305 nucleotides from the translation start site was sufficient for specific expression in both Arabidopsis and tobacco roots. Gain-of-function analysis revealed functional regions at -305/-169 and -168/-93, whose coexistence was required for the induction activity in Arabidopsis roots. Multiple IDE-like sequences were distributed in the IDS3 promoter and were especially abundant within the functional region at -305/-169. A sequence moderately homologous to that of IDE1 was also present within the -168/-93 region. These IDE-like sequences would be the first candidates for the functional iron-deficiency-responsive elements in the IDS3 promoter.     

Analysis of sequence, map position, and gene expression reveals conserved essential genes for iron uptake in Arabidopsis and tomato

Arabidopsis (Arabidopsis thaliana) and tomato (Lycopersicon esculentum) show similar physiological responses to iron deficiency, suggesting that homologous genes are involved. Essential gene functions are generally considered to be carried out by orthologs that have remained conserved in sequence and map position in evolutionarily related species. This assumption has not yet been proven for plant genomes that underwent large genome rearrangements. We addressed this question in an attempt to deduce functional gene pairs for iron reduction, iron transport, and iron regulation between Arabidopsis and tomato. Iron uptake processes are essential for plant growth. We investigated iron uptake gene pairs from tomato and Arabidopsis, namely sequence, conserved gene content of the regions containing iron uptake homologs based on conserved orthologous set marker analysis, gene expression patterns, and, in two cases, genetic data. Compared to tomato, the Arabidopsis genome revealed more and larger gene families coding for the iron uptake functions. The number of possible homologous pairs was reduced if functional expression data were taken into account in addition to sequence and map position. We predict novel homologous as well as partially redundant functions of ferric reductase-like and iron-regulated transporter-like genes in Arabidopsis and tomato. Arabidopsis nicotianamine synthase genes encode a partially redundant family. In this study, Arabidopsis gene redundancy generally reflected the presumed genome duplication structure. In some cases, statistical analysis of conserved gene regions between tomato and Arabidopsis suggested a common evolutionary origin. Although involvement of conserved genes in iron uptake was found, these essential genes seem to be of paralogous rather than orthologous origin in tomato and Arabidopsis.     

The Arabidopsis Information Resource (TAIR): a model organism database providing a centralized,curated gateway to Arabidopsis biology,research materials and community

Arabidopsis thaliana is the most widely-studied plant today. The concerted efforts of over 11 000 researchers and 4000 organizations around the world are generating a rich diversity and quantity of information and materials. This information is made available through a comprehensive on-line resource called the Arabidopsis Information Resource (TAIR) (http://arabidopsis.org), which is accessible via commonly used web browsers and can be searched and downloaded in a number of ways. In the last two years, efforts have been focused on increasing data content and diversity, functionally annotating genes and gene products with controlled vocabularies, and improving data retrieval, analysis and visualization tools. New information include sequence polymorphisms including alleles, germplasms and phenotypes, Gene Ontology annotations, gene families, protein information, metabolic pathways, gene expression data from microarray experiments and seed and DNA stocks. New data visualization and analysis tools include SeqViewer, which interactively displays the genome from the whole chromosome down to 10 kb of nucleotide sequence and AraCyc, a metabolic pathway database and map tool that allows overlaying expression data onto the pathway diagrams. Finally, we have recently incorporated seed and DNA stock information from the Arabidopsis Biological Resource Center (ABRC) and implemented a shopping-cart style on-line ordering system.     

Small Glycosylated Lignin Oligomers Are Stored in Arabidopsis Leaf Vacuoles

Lignin is an aromatic polymer derived from the combinatorial coupling of monolignol radicals in the cell wall. Recently, various glycosylated lignin oligomers have been revealed in Arabidopsis thaliana. Given that monolignol oxidation and monolignol radical coupling are known to occur in the apoplast, and glycosylation in the cytoplasm, it raises questions about the subcellular localization of glycosylated lignin oligomer biosynthesis and their storage. By metabolite profiling of Arabidopsis leaf vacuoles, we show that the leaf vacuole stores a large number of these small glycosylated lignin oligomers. Their structural variety and the incorporation of alternative monomers, as observed in Arabidopsis mutants with altered monolignol biosynthesis, indicate that they are all formed by combinatorial radical coupling. In contrast to the common believe that combinatorial coupling is restricted to the apoplast, we hypothesized that the aglycones of these compounds are made within the cell. To investigate this, leaf protoplast cultures were cofed with ¹³C₆-labeled coniferyl alcohol and a ¹³C₄-labeled dimer of coniferyl alcohol. Metabolite profiling of the cofed protoplasts provided strong support for the occurrence of intracellular monolignol coupling. We therefore propose a metabolic pathway involving intracellular combinatorial coupling of monolignol radicals, followed by oligomer glycosylation and vacuolar import, which shares characteristics with both lignin and lignan biosynthesis.     

Genomic analysis of the terpenoid synthase (<Emphasis Type="Italic">AtTPS</Emphasis>) gene family of<Emphasis Type="Italic"> Arabidopsis thaliana</Emphasis>

A family of 40 terpenoid synthase genes ( AtTPS) was discovered by genome sequence analysis in Arabidopsis thaliana. This is the largest and most diverse group of TPS genes currently known for any species. AtTPS genes cluster into five phylogenetic subfamilies of the plant TPS superfamily. Surprisingly, thirty AtTPS closely resemble, in all aspects of gene architecture, sequence relatedness and phylogenetic placement, the genes for plant monoterpene synthases, sesquiterpene synthases or diterpene synthases of secondary metabolism. Rapid evolution of these AtTPS resulted from repeated gene duplication and sequence divergence with minor changes in gene architecture. In contrast, only two AtTPS genes have known functions in basic (primary) metabolism, namely gibberellin biosynthesis. This striking difference in rates of gene diversification in primary and secondary metabolism is relevant for an understanding of the evolution of terpenoid natural product diversity. Eight AtTPS genes are interrupted and are likely to be inactive pseudogenes. The localization of AtTPS genes on all five chromosomes reflects the dynamics of the Arabidopsis genome; however, several AtTPS genes are clustered and organized in tandem repeats. Furthermore, some AtTPS genes are localized with prenyltransferase genes ( AtGGPPS, geranylgeranyl diphosphate synthase) in contiguous genomic clusters encoding consecutive steps in terpenoid biosynthesis. The clustered organization may have implications for TPS gene evolution and the evolution of pathway segments for the synthesis of terpenoid natural products. Phylogenetic analyses highlight events in the divergence of the TPS paralogs and suggest orthologous genes and a model for the evolution of the TPS gene family.