The plant glutathione transferase gene family: genomic structure, functions, expression and evolution

Glutathione transferases (GSTs) are ubiquitous, multifunctional proteins encoded by large gene families. In different plant species this gene family is comprised of 25–60 members, that can be grouped into six classes on the basis of sequence identity, gene organization and active site residues in the protein. The Phi and Tau classes are the most represented and are plant specific, while Zeta and Theta GSTs are found also in animals. Despite pronounced sequence and functional diversification, GSTs have maintained a highly conserved three-dimensional structure through evolution. Most GSTs are cytosolic and active as dimers, performing diverse catalytic as well as non-catalytic roles in detoxification of xenobiotics, prevention of oxidative damage and endogenous metabolism. Among their catalytic activities are the conjugation of electrophilic substrates to glutathione, glutathione-dependent isomerizations and reductions of toxic organic hydroperoxides. Their main non-catalytic role is as hormone and flavonoid ligandins. GST genes are predominantly organized in clusters non-randomly distributed in the genome. Phylogenetic studies indicate that plant GSTs have mainly evolved after the divergence of plants, the two prevalent Phi and Tau classes being the result of recent, multiple duplication events.     

The Egg apparatus 1 gene from maize is a member of a large gene family found in both monocots and dicots

The maize ZmEA1 protein was recently postulated to be involved in short-range pollen tube guidance from the embryo sac. To date, EA1-like sequences had only been identified in monocot species. Using a more conserved C-terminal motif found in the monocot species, numerous ZmEA1-like sequences were retrieved in EST databases from dicot species, as well as from unannotated genomic sequences of Arabidopsis thaliana. RT-PCR analyses were produced for these unannotated genes and showed that these were indeed expressed genes. Further structural and phylogenetic analyses revealed that all members of the EA1-like (EAL) gene family shared a conserved 27–29 amino acid motif, termed the EA box near the C-terminal end, and appear to be secretory proteins. Therefore, the EA box proteins defines a new class of small secretory proteins, some of which being possibly involved in pollen tube guidance. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Phylogenomic Analysis of the PEBP Gene Family in Cereals

The TFL1 and FT genes, which are key genes in the control of flowering time in Arabidopsis thaliana, belong to a small multigene family characterized by a specific phosphatidylethanolamine-binding protein domain, termed the PEBP gene family. Several PEBP genes are found in dicots and monocots, and act on the control of flowering time. We investigated the evolution of the PEBP gene family in cereals. First, taking advantage of the complete rice genome sequence and EST databases, we found 19 PEBP genes in this species, 6 of which were not previously described. Ten genes correspond to five pairs of paralogs mapped on known duplicated regions of the rice genome. Phylogenetic analysis of Arabidopsis and rice genes indicates that the PEBP gene family consists of three main homology classes (the so-called TFL1-LIKE, MFT-LIKE, and FT-LIKE subfamilies), in which gene duplication and/or loss occurred independently in Arabidopsis and rice. Second, phylogenetic analyses of genomic and EST sequences from five cereal species indicate that the three subfamilies of PEBP genes have been conserved in cereals. The tree structure suggests that the ancestral grass genome had at least two MFT-like genes, two TFL1-like genes, and eight FT-like genes. A phylogenomic approach leads to some hypotheses about conservation of gene function within the subfamilies. [Reviewing Editor: Dr. Yves Van de Peer]     

Whole genome survey of the glutathione transferase family in the brown algal model Ectocarpus siliculosus

We report here an exhaustive analysis of the glutathione transferases (GSTs) in the model brown alga Ectocarpus siliculosus using available genomic resources. A genome survey revealed the presence of twelve cytosolic GSTs, belonging to the Sigma class, two pseudogenes, one GST of the Kappa class, and three microsomal GSTs of the MGST3 family of membrane associated protein involved in eicosanoid and glutathione metabolism. Gene structure and phylogenetic analyses demonstrated the partition of the Sigma GSTs into two clusters which have probably evolved by duplication events. Gene expression profiling was conducted after the addition of high concentrations of chemicals, such as H(2)O(2), herbicides, heavy metals, as well as fatty acid derivatives, in order to induce stress conditions and to monitor early response mechanisms. The results of these experiments suggested that E. siliculosus GST genes are recruited in different and specific conditions. In addition, heterologous expression in yeast of two E. siliculosus microsomal GST showed that these enzymes feature peroxidase rather than transferase activity. The potential involvement of E. siliculosus GST in the metabolism of oxygenated polyunsaturated fatty acids is discussed.     

A Tale of Plant Glutathione <Emphasis Type="Italic">S</Emphasis>-Transferases: Since 1970

Ubiquitously distributed multifunctional superfamily of Glutathione S-transferases (GST) generally constitute a dimeric enzymes and catalyse the conjugation of the thiol group of the glutathione (GSH) to diverse electrophilic centres on lipophilic molecules with the formation of rather less active end products. Besides their well investigated conjugation reaction for the detoxification of endogenous and xenobiotic compounds, they can also be involved in both GSH dependent peroxidation or isomerization reactions, and several other non-catalytic functions, like binding of non-substrate ligands, stress-induced signalling processes and preventing of apoptosis. Plant GSTs have been a focus of attention because of their roles in herbicide detoxification and today seven distinct classes of soluble (cytosolic) GSTs are presented as Phi, Tau, Theta, Zeta, Lambda, Dehydroascorbate reductases (DHARs) and Tetrachlorohydroquinone dehalogenase (TCHQD). While GSTs show overall sequence diversification within and between classes, they retain a high level of three-dimensional structure conservation over long evolutionary periods. In this review mainly the soluble plant GSTs will be considered by giving attention to their structures, subcellular localizations, genomic organizations, catalytic/noncatalytic functions, and comparisons given with respect to their mammalian counterparts where necessary.     

Proteomic Analysis of Glutathione S-Transferases of Arabidopsis thaliana Reveals Differential Salicylic Acid-Induced Expression of the Plant-Specific Phi and Tau Classes

Plant glutathione S -transferases (GSTs) are a large group of multifunctional proteins that are induced by diverse stimuli. Using proteomic approaches we identified 20 GSTs at the protein level in Arabidopsis cell culture with a combination of GST antibody detection, LC-MS/MS analysis of 23-30 kDa proteins and glutathione-affinity chromatography. GSTs identified were from phi, tau, theta, zeta and DHAR sub-sections of the GST superfamily of 53 members. We have uncovered preliminary evidence for post-translational modifications of plant GSTs and show that phosphorylation is unlikely to be responsible. Detailed analysis of GST expression in response to treatment with 0.01-1 mM of the plant defence signal salicylic acid (SA) uncovered some interesting features. Firstly, GSTs appear to display class-specific concentration-dependent SA induction profiles highlighting differences between the large, plant specific phi and tau classes. Secondly, different members of the same class, while sharing similar SA dose responses, may display differences in terms of magnitude and timing of induction, further highlighting the breadth of GST gene regulation. Thirdly, closely related members of the same class ( GSTF6 and GSTF7 ), arising via tandem duplication, may be regulated differently in terms of basal expression levels and also magnitude of induction raising questions about the role of subfunctionalisation within this family. Our results reveal that GSTs exhibit class specific responses to SA treatment suggesting that several mechanisms are acting to induce GSTs upon SA treatment and hinting at class-specific functions for this large and important, yet still relatively elusive gene family.     

The ankyrin repeat gene family in rice: genome-wide identification,classification and expression profiling

Ankyrin repeat (ANK) containing proteins comprise a large protein family. Although many members of this family have been implicated in plant growth, development and signal transduction, only a few ANK genes have been reported in rice. In this study, we analyzed the structures, phylogenetic relationship, genome localizations and expression profiles of 175 ankyrin repeat genes identified in rice (OsANK). Domain composition analysis suggested OsANK proteins can be classified into ten subfamilies. Chromosomal localizations of OsANK genes indicated nine segmental duplication events involving 17 genes and 65 OsANK genes were involved in tandem duplications. The expression profiles of 158 OsANK genes were analyzed in 24 tissues covering the whole life cycle of two rice genotypes, Minghui 63 and Zhenshan 97. Sixteen genes showed preferential expression in given tissues compared to all the other tissues in Minghui 63 and Zhenshan 97. Nine genes were preferentially expressed in stamen of 1 day before flowering, suggesting that these genes may play important roles in pollination and fertilization. Expression data of OsANK genes were also obtained with tissues of seedlings subjected to three phytohormone (NAA, GA3 and KT) and light/dark treatments. Eighteen genes showed differential expression with at least one phytohormone treatment while under light/dark treatments, 13 OsANK genes showed differential expression. Our data provided a very useful reference for cloning and functional analysis of members of this gene family in rice. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Genome-wide analysis of the auxin response factor (ARF) gene family in maize (<Emphasis Type="Italic">Zea mays</Emphasis>)

Auxin response factors (ARFs) are an important family involved in auxin-mediated response through specific binding to auxin response elements (AuxREs). A few members of the ARF family have been functionally characterized in Arabidopsis, rice (Oryza sativa), Poplar (Populous trichocarpa). However, little is known about ARF genes in maize (Zea mays). We performed a comprehensive bioinformatics analysis of the maize ARF gene family including analysis of the genome sequence, conserved domains, chromosomal locations, phylogenetic relationships, gene duplication, and expression profiles. 35 ZmARF genes were identified and categorized into four groups (Class I, II, III, and IV). In addition, a segmental ZmARF duplication event was shown to play an important role in maize ARF gene expansion. 7 ZmARF genes had no expression in specific tissues we obtained, but presented in mixed tissues according to the NCBI EST database, respectively. These studies have laid the theoretical foundation for further functional verification of these ZmARF genes.     

Comparative genomic analysis of CIPK gene family in <Emphasis Type="Italic">Arabidopsis</Emphasis> and <Emphasis Type="Italic">Populus</Emphasis>

Calcium serves as a second messenger in various signal transduction pathways in plants. CBL-interacting protein kinases (CIPKs), which have a variety of functions, are involved in calcium signal transduction. Previous, the studies on CIPK family members focused on Arabidopsis and rice. Here, we present a comparative genomic analysis of the CIPK gene family in Arabidopsis and poplar, a model tree species. Twenty-seven potential CIPKs were identified from poplar using genome-wide analysis. Like the CIPK gene family from Arabidopsis, CIPK genes from poplar were also divided into intron-free and intron-harboring groups. In the intron-harboring group, the intron distribution of CIPKs is rather conserved during the genome evolutionary process. Many homologous gene pairs were found in the CIPK gene family, indicating duplication events might contribute to the amplification of this gene family. The phylogenetic comparison of CIPKs in combination with intron distribution analysis revealed that CIPK genes from both Arabidopsis and poplar might have an ancient origin, which formed earlier than the separation of these two eudicot species. Our genomic and bioinformatic analysis will provide an important foundation for further functional dissection of the CBL-CIPK signaling network in poplars. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The <Emphasis Type="Italic">Anopheles gambiae</Emphasis> glutathione transferase supergene family: annotation,phylogeny and expression profiles

Background  

Twenty-eight genes putatively encoding cytosolic glutathione transferases have been identified in the Anopheles gambiae genome. We manually annotated these genes and then confirmed the annotation by sequencing of A. gambiae cDNAs. Phylogenetic analysis with the 37 putative GST genes from Drosophila and representative GSTs from other taxa was undertaken to develop a nomenclature for insect GSTs. The epsilon class of insect GSTs has previously been implicated in conferring insecticide resistance in several insect species. We compared the expression level of all members of this GST class in two strains of A. gambiae to determine whether epsilon GST expression is correlated with insecticide resistance status.     

Marine Glutathione <Emphasis Type="Italic">S</Emphasis>-Transferases

The aquatic environment is generally affected by the presence of environmental xenobiotic compounds. One of the major xenobiotic detoxifying enzymes is glutathione S-transferase (GST), which belongs to a family of multifunctional enzymes involved in catalyzing nucleophilic attack of the sulfur atom of glutathione (γ-glutamyl-cysteinylglycine) to an electrophilic group on metabolic products or xenobiotic compounds. Because of the unique nature of the aquatic environment and the possible pollution therein, the biochemical evolution in terms of the nature of GSTs could by uniquely expressed. The full complement of GSTs has not been studied in marine organisms, as very few aquatic GSTs have been fully characterized. The focus of this article is to present an overview of the GST superfamily and their critical role in the survival of organisms in the marine environment, emphasizing the critical roles of GSTs in the detoxification of marine organisms and the unique characteristics of their GSTs compared to those from non-marine organisms.     

Genome-wide analysis of defense-responsive genes in bacterial blight resistance of rice mediated by the recessive<Emphasis Type="Italic"> R</Emphasis> gene<Emphasis Type="Italic"> xa13</Emphasis>

Defense responses triggered by dominant and recessive disease resistance ( R) genes are presumed to be regulated by different molecular mechanisms. In order to characterize the genes activated in defense responses against bacterial blight mediated by the recessive R gene xa13, two pathogen-induced subtraction cDNA libraries were constructed using the resistant rice line IRBB13—which carries xa13 —and its susceptible, near-isogenic, parental line IR24. Clustering analysis of expressed sequence tags (ESTs) identified 702 unique expressed sequences as being involved in the defense responses triggered by xa13; 16% of these are new rice ESTs. These sequences define 702 genes, putatively encoding a wide range of products, including defense-responsive genes commonly involved in different host-pathogen interactions, genes that have not previously been reported to be associated with pathogen-induced defense responses, and genes (38%) with no homology to previously described functional genes. In addition, R -like genes putatively encoding nucleotide-binding site/leucine rich repeat (NBS-LRR) and LRR receptor kinase proteins were observed to be induced in the disease resistance activated by xa13. A total of 568 defense-responsive ESTs were mapped to 588 loci on the rice molecular linkage map through bioinformatic analysis. About 48% of the mapped ESTs co-localized with quantitative trait loci (QTLs) for resistance to various rice diseases, including bacterial blight, rice blast, sheath blight and yellow mottle virus. Furthermore, some defense-responsive sequences were conserved at similar locations on different chromosomes. These results reveal the complexity of xa13 -mediated resistance. The information obtained in this study provides a large source of candidate genes for understanding the molecular bases of defense responses activated by recessive R genes and of quantitative disease resistance.Electronic Supplementary Material Supplementary material is available in the online version of this article at The first two authors contributed equally to this workCommunicated by R. Hagemann     

Plant glutathione transferases

The soluble glutathione transferases (GSTs, EC 2.5.1.18) are encoded by a large and diverse gene family in plants, which can be divided on the basis of sequence identity into the phi, tau, theta, zeta and lambda classes. The theta and zeta GSTs have counterparts in animals but the other classes are plant-specific and form the focus of this article. The genome of Arabidopsis thaliana contains 48 GST genes, with the tau and phi classes being the most numerous. The GST proteins have evolved by gene duplication to perform a range of functional roles using the tripeptide glutathione (GSH) as a cosubstrate or coenzyme. GSTs are predominantly expressed in the cytosol, where their GSH-dependent catalytic functions include the conjugation and resulting detoxification of herbicides, the reduction of organic hydroperoxides formed during oxidative stress and the isomerization of maleylacetoacetate to fumarylacetoacetate, a key step in the catabolism of tyrosine. GSTs also have non-catalytic roles, binding flavonoid natural products in the cytosol prior to their deposition in the vacuole. Recent studies have also implicated GSTs as components of ultraviolet-inducible cell signaling pathways and as potential regulators of apoptosis. Although sequence diversification has produced GSTs with multiple functions, the structure of these proteins has been highly conserved. The GSTs thus represent an excellent example of how protein families can diversify to fulfill multiple functions while conserving form and structure.