Evolutionary, molecular and genetic analyses of Tic22 homologues in Arabidopsis thaliana chloroplasts

The Tic22 protein was previously identified in pea as a putative component of the chloroplast protein import apparatus. It is a peripheral protein of the inner envelope membrane, residing in the intermembrane space. In Arabidopsis, there are two Tic22 homologues, termed atTic22-III and atTic22-IV, both of which are predicted to localize in chloroplasts. These two proteins defined clades that are conserved in all land plants, which appear to have evolved at a similar rates since their separation >400 million years ago, suggesting functional conservation. The atTIC22-IV gene was expressed several-fold more highly than atTIC22-III, but the genes exhibited similar expression profiles and were expressed throughout development. Knockout mutants lacking atTic22-IV were visibly normal, whereas those lacking atTic22-III exhibited moderate chlorosis. Double mutants lacking both isoforms were more strongly chlorotic, particularly during early development, but were viable and fertile. Double-mutant chloroplasts were small and under-developed relative to those in wild type, and displayed inefficient import of precursor proteins. The data indicate that the two Tic22 isoforms act redundantly in chloroplast protein import, and that their function is non-essential but nonetheless required for normal chloroplast biogenesis, particularly during early plant development.

An unusual dimeric small heat shock protein provides insight into the mechanism of this class of chaperones

Small heat shock proteins (sHSPs) are virtually ubiquitous stress proteins that are also found in many normal tissues and accumulate in diseases of protein folding. They generally act as ATP-independent chaperones to bind and stabilize denaturing proteins that can be later reactivated by ATP-dependent Hsp70/DnaK, but the mechanism of substrate capture by sHSPs remains poorly understood. A majority of sHSPs form large oligomers, a property that has been linked to their effective chaperone action. We describe AtHsp18.5 from Arabidopsis thaliana, demonstrating that it is dimeric and exhibits robust chaperone activity, which adds support to the model that suboligomeric sHSP forms are a substrate binding species. Notably, like oligomeric sHSPs, when bound to substrate, AtHsp18.5 assembles into large complexes, indicating that reformation of sHSP oligomeric contacts is not required for assembly of sHSP-substrate complexes. Monomers of AtHsp18.5 freely exchange between dimers but fail to coassemble in vitro with dodecameric plant cytosolic sHSPs, suggesting that AtHsp18.5 does not interact by coassembly with these other sHSPs in vivo. Data from controlled proteolysis and hydrogen-deuterium exchange coupled with mass spectrometry show that the N- and C-termini of AtHsp18.5 are highly accessible and lack stable secondary structure, most likely a requirement for substrate interaction. Chaperone activity of a series of AtHsp18.5 truncation mutants confirms that the N-terminal arm is required for substrate protection and that different substrates interact differently with the N-terminal arm. In total, these data imply that the core α-crystallin domain of the sHSPs is a platform for flexible arms that capture substrates to maintain their solubility.

A survey of dominant mutations in Arabidopsis thaliana

Following the recent publication of a comprehensive dataset of 2400 genes with a loss-of-function mutant phenotype in Arabidopsis (Arabidopsis thaliana), questions remain concerning the diversity of dominant mutations in Arabidopsis. Most of these dominant phenotypes are expected to result from inappropriate gene expression, novel protein function, or disrupted protein complexes. This review highlights the major classes of dominant mutations observed in model organisms and presents a collection of 200 Arabidopsis genes associated with a dominant or semidominant phenotype. Emphasis is placed on mutants identified through forward genetic screens of mutagenized or activation-tagged populations. These datasets illustrate the variety of genetic changes and protein functions that underlie dominance in Arabidopsis and may ultimately contribute to phenotypic variation in flowering plants.

Transcriptional and functional classification of the GOLVEN/ROOT GROWTH FACTOR/CLE-like signaling peptides reveals their role in lateral root and hair formation

The GOLVEN (GLV)/ROOT GROWTH FACTORS/CLE-Like small signaling peptide family is encoded by 11 genes in Arabidopsis (Arabidopsis thaliana). Some of them have already been shown to control root meristem maintenance, auxin fluxes, and gravitropic responses. As a basis for the detailed analysis of their function, we determined the expression domains for each of the 11 GLV genes with promoter-reporter lines. Although they are collectively active in all examined plant parts, GLV genes have highly specific transcription patterns, generally restricted to very few cells or cell types in the root and shoot and in vegetative and reproductive tissues. GLV functions were further investigated with the comparative analysis of root phenotypes induced by gain- and loss-of-function mutants or in treatments with GLV-derived synthetic peptides. We identified functional classes that relate to the gene expression domains in the primary root and suggest that different GLV signals trigger distinct downstream pathways. Interestingly, GLV genes transcribed at the early stages of lateral root development strongly inhibited root branching when overexpressed. Furthermore, transcription patterns together with mutant phenotypes pointed to the involvement of GLV4 and GLV8 in root hair formation. Overall, our data suggest that nine GLV genes form three subgroups according to their expression and function within the root and offer a comprehensive framework to study the role of the GLV signaling peptides in plant development.

Recruitment of Arf1-GDP to Golgi by Glo3p-type ArfGAPs is crucial for golgi maintenance and plant growth

ADP-ribosylation factor1 (Arf1), a member of the small GTP-binding proteins, plays a pivotal role in protein trafficking to multiple organelles. In its GDP-bound form, Arf1 is recruited from the cytosol to organelle membranes, where it functions in vesicle-mediated protein trafficking. However, the mechanism of Arf1-GDP recruitment remains unknown. Here, we provide evidence that two Glo3p-type Arf GTPase-activating proteins (ArfGAPs), ArfGAP domain8 (AGD8) and AGD9, are involved in the recruitment of Arf1-GDP to the Golgi apparatus in Arabidopsis (Arabidopsis thaliana). RNA interference plants expressing low levels of AGD8 and AGD9 exhibited abnormal Golgi morphology, inhibition of protein trafficking, and arrest of plant growth and development. In RNA interference plants, Arf1 was poorly recruited to the Golgi apparatus. Conversely, high levels of AGD8 and AGD9 induced Arf1 accumulation at the Golgi and suppressed Golgi disruption and inhibition of vacuolar trafficking that was caused by overexpression of AGD7. Based on these results, we propose that the Glo3p-type ArfGAPs AGD8 and AGD9 recruit Arf1-GDP from the cytosol to the Golgi for Arf1-mediated protein trafficking, which is essential for plant development and growth.

The association of a La module with the PABP-interacting motif PAM2 is a recurrent evolutionary process that led to the neofunctionalization of La-related proteins

La-related proteins (LARPs) are largely uncharacterized factors, well conserved throughout evolution. Recent reports on the function of human LARP4 and LARP6 suggest that these proteins fulfill key functions in mRNA metabolism and/or translation. We report here a detailed evolutionary history of the LARP4 and 6 families in eukaryotes. Genes coding for LARP4 and 6 were duplicated in the common ancestor of the vertebrate lineage, but one LARP6 gene was subsequently lost in the common ancestor of the eutherian lineage. The LARP6 gene was also independently duplicated several times in the vascular plant lineage. We observed that vertebrate LARP4 and plant LARP6 duplication events were correlated with the acquisition of a PABP-interacting motif 2 (PAM2) and with a significant reorganization of their RNA-binding modules. Using isothermal titration calorimetry (ITC) and immunoprecipitation methods, we show that the two plant PAM2-containing LARP6s (LARP6b and c) can, indeed, interact with the major plant poly(A)-binding protein (PAB2), while the third plant LARP6 (LARP6a) is unable to do so. We also analyzed the RNA-binding properties and the subcellular localizations of the two types of plant LARP6 proteins and found that they display nonredundant characteristics. As a whole, our results support a model in which the acquisition by LARP4 and LARP6 of a PAM2 allowed their targeting to mRNA 3' UTRs and led to their neofunctionalization.

Functional divergence in the Arabidopsis LOB-domain gene family

The Arabidopsis LOB-domain (LBD) gene family is composed by 43 members divided in two classes based on amino acid conservation within the LOB-domain. The LOB domain is known to be responsible for DNA binding and protein-protein interactions. There is very little functional information available for most genes in the LBD family and many lbd single mutants do not exhibit conspicuous phenotypes. One plausible explanation for the limited loss-of-function phenotypes observed in this family is that LBD genes exhibit significant functional redundancy. Here we discuss an example of one phylogenetic subgroup of the LBD family, in which genes that are closely related based on phylogeny exhibit distinctly different expression patterns and do not have overlapping functions. We discuss the challenges of using phylogenetic analyses to predict redundancy in gene families.

Genome-wide classification and expression analysis of MYB transcription factor families in rice and Arabidopsis

BACKGROUND

The MYB gene family comprises one of the richest groups of transcription factors in plants. Plant MYB proteins are characterized by a highly conserved MYB DNA-binding domain. MYB proteins are classified into four major groups namely, 1R-MYB, 2R-MYB, 3R-MYB and 4R-MYB based on the number and position of MYB repeats. MYB transcription factors are involved in plant development, secondary metabolism, hormone signal transduction, disease resistance and abiotic stress tolerance. A comparative analysis of MYB family genes in rice and Arabidopsis will help reveal the evolution and function of MYB genes in plants.

RESULTS

A genome-wide analysis identified at least 155 and 197 MYB genes in rice and Arabidopsis, respectively. Gene structure analysis revealed that MYB family genes possess relatively more number of introns in the middle as compared with C- and N-terminal regions of the predicted genes. Intronless MYB-genes are highly conserved both in rice and Arabidopsis. MYB genes encoding R2R3 repeat MYB proteins retained conserved gene structure with three exons and two introns, whereas genes encoding R1R2R3 repeat containing proteins consist of six exons and five introns. The splicing pattern is similar among R1R2R3 MYB genes in Arabidopsis. In contrast, variation in splicing pattern was observed among R1R2R3 MYB members of rice. Consensus motif analysis of 1kb upstream region (5' to translation initiation codon) of MYB gene ORFs led to the identification of conserved and over-represented cis-motifs in both rice and Arabidopsis. Real-time quantitative RT-PCR analysis showed that several members of MYBs are up-regulated by various abiotic stresses both in rice and Arabidopsis.

CONCLUSION

A comprehensive genome-wide analysis of chromosomal distribution, tandem repeats and phylogenetic relationship of MYB family genes in rice and Arabidopsis suggested their evolution via duplication. Genome-wide comparative analysis of MYB genes and their expression analysis identified several MYBs with potential role in development and stress response of plants.

Genome-wide classification and expression analysis of MYB transcription factor families in rice and Arabidopsis

BACKGROUND

The MYB gene family comprises one of the richest groups of transcription factors in plants. Plant MYB proteins are characterized by a highly conserved MYB DNA-binding domain. MYB proteins are classified into four major groups namely, 1R-MYB, 2R-MYB, 3R-MYB and 4R-MYB based on the number and position of MYB repeats. MYB transcription factors are involved in plant development, secondary metabolism, hormone signal transduction, disease resistance and abiotic stress tolerance. A comparative analysis of MYB family genes in rice and Arabidopsis will help reveal the evolution and function of MYB genes in plants.

RESULTS

A genome-wide analysis identified at least 155 and 197 MYB genes in rice and Arabidopsis, respectively. Gene structure analysis revealed that MYB family genes possess relatively more number of introns in the middle as compared with C- and N-terminal regions of the predicted genes. Intronless MYB-genes are highly conserved both in rice and Arabidopsis. MYB genes encoding R2R3 repeat MYB proteins retained conserved gene structure with three exons and two introns, whereas genes encoding R1R2R3 repeat containing proteins consist of six exons and five introns. The splicing pattern is similar among R1R2R3 MYB genes in Arabidopsis. In contrast, variation in splicing pattern was observed among R1R2R3 MYB members of rice. Consensus motif analysis of 1kb upstream region (5' to translation initiation codon) of MYB gene ORFs led to the identification of conserved and over-represented cis-motifs in both rice and Arabidopsis. Real-time quantitative RT-PCR analysis showed that several members of MYBs are up-regulated by various abiotic stresses both in rice and Arabidopsis.

CONCLUSION

A comprehensive genome-wide analysis of chromosomal distribution, tandem repeats and phylogenetic relationship of MYB family genes in rice and Arabidopsis suggested their evolution via duplication. Genome-wide comparative analysis of MYB genes and their expression analysis identified several MYBs with potential role in development and stress response of plants.

Arabidopsis PIZZA has the capacity to acylate brassinosteroids

Brassinosteroids (BRs) affect a wide range of developmental processes in plants and compromised production or signalling of BRs causes severe growth defects. To identify new regulators of plant organ growth, we searched the Arabidopsis FOX (Full-length cDNA Over-eXpressor gene) collection for mutants with altered organ size and isolated two overexpression lines that display typical BR deficient dwarf phenotypes. The phenotype of these lines, caused by an overexpression of a putative acyltransferase gene PIZZA (PIZ), was partly rescued by supplying exogenous brassinolide (BL) and castasterone (CS), indicating that endogenous BR levels are rate-limiting for the growth of PIZ overexpression lines. Our transcript analysis further showed that PIZ overexpression leads to an elevated expression of genes involved in BR biosynthesis and a reduced expression of BR inactivating hydroxylases, a transcriptional response typical to low BR levels. Taking the advantage of relatively high endogenous BR accumulation in a mild bri1-301 background, we found that overexpression of PIZ results in moderately reduced levels of BL and CS and a strong reduction of typhasterol (TY) and 6-deoxocastasterone (6-deoxoCS), suggesting a role of PIZ in BR metabolism. We tested a set of potential substrates in vitro for heterologously expressed PIZ and confirmed its acyltransferase activity with BL, CS and TY. The PIZ gene is expressed in various tissues but as reported for other genes involved in BR metabolism, the loss-of-function mutants did not display obvious growth phenotypes under standard growth conditions. Together, our data suggest that PIZ can modify BRs by acylation and that these properties might help modulating endogenous BR levels in Arabidopsis.

A land-plant-specific glycerol-3-phosphate acyltransferase family in Arabidopsis: substrate specificity, sn-2 preference, and evolution

Arabidopsis (Arabidopsis thaliana) has eight glycerol-3-phosphate acyltransferase (GPAT) genes that are members of a plant-specific family with three distinct clades. Several of these GPATs are required for the synthesis of cutin or suberin. Unlike GPATs with sn-1 regiospecificity involved in membrane or storage lipid synthesis, GPAT4 and -6 are unique bifunctional enzymes with both sn-2 acyltransferase and phosphatase activity resulting in 2-monoacylglycerol products. We present enzymology, pathway organization, and evolutionary analysis of this GPAT family. Within the cutin-associated clade, GPAT8 is demonstrated as a bifunctional sn-2 acyltransferase/phosphatase. GPAT4, -6, and -8 strongly prefer C16:0 and C18:1 ω-oxidized acyl-coenzyme As (CoAs) over unmodified or longer acyl chain substrates. In contrast, suberin-associated GPAT5 can accommodate a broad chain length range of ω-oxidized and unsubstituted acyl-CoAs. These substrate specificities (1) strongly support polyester biosynthetic pathways in which acyl transfer to glycerol occurs after oxidation of the acyl group, (2) implicate GPAT specificities as one major determinant of cutin and suberin composition, and (3) argue against a role of sn-2-GPATs (Enzyme Commission 2.3.1.198) in membrane/storage lipid synthesis. Evidence is presented that GPAT7 is induced by wounding, produces suberin-like monomers when overexpressed, and likely functions in suberin biosynthesis. Within the third clade, we demonstrate that GPAT1 possesses sn-2 acyltransferase but not phosphatase activity and can utilize dicarboxylic acyl-CoA substrates. Thus, sn-2 acyltransferase activity extends to all subbranches of the Arabidopsis GPAT family. Phylogenetic analyses of this family indicate that GPAT4/6/8 arose early in land-plant evolution (bryophytes), whereas the phosphatase-minus GPAT1 to -3 and GPAT5/7 clades diverged later with the appearance of tracheophytes.

WNP: a novel algorithm for gene products annotation from weighted functional networks

Predicting the biological function of all the genes of an organism is one of the fundamental goals of computational system biology. In the last decade, high-throughput experimental methods for studying the functional interactions between gene products (GPs) have been combined with computational approaches based on Bayesian networks for data integration. The result of these computational approaches is an interaction network with weighted links representing connectivity likelihood between two functionally related GPs. The weighted network generated by these computational approaches can be used to predict annotations for functionally uncharacterized GPs. Here we introduce Weighted Network Predictor (WNP), a novel algorithm for function prediction of biologically uncharacterized GPs. Tests conducted on simulated data show that WNP outperforms other 5 state-of-the-art methods in terms of both specificity and sensitivity and that it is able to better exploit and propagate the functional and topological information of the network. We apply our method to Saccharomyces cerevisiae yeast and Arabidopsis thaliana networks and we predict Gene Ontology function for about 500 and 10000 uncharacterized GPs respectively.

Plant vegetative and animal cytoplasmic actins share functional competence for spatial development with protists

Actin is an essential multifunctional protein encoded by two distinct ancient classes of genes in animals (cytoplasmic and muscle) and plants (vegetative and reproductive). The prevailing view is that each class of actin variants is functionally distinct. However, we propose that the vegetative plant and cytoplasmic animal variants have conserved functional competence for spatial development inherited from an ancestral protist actin sequence. To test this idea, we ectopically expressed animal and protist actins in Arabidopsis thaliana double vegetative actin mutants that are dramatically altered in cell and organ morphologies. We found that expression of cytoplasmic actins from humans and even a highly divergent invertebrate Ciona intestinalis qualitatively and quantitatively suppressed the root cell polarity and organ defects of act8 act7 mutants and moderately suppressed the root-hairless phenotype of act2 act8 mutants. By contrast, human muscle actins were unable to support prominently any aspect of plant development. Furthermore, actins from three protists representing Choanozoa, Archamoeba, and green algae efficiently suppressed all the phenotypes of both the plant mutants. Remarkably, these data imply that actin's competence to carry out a complex suite of processes essential for multicellular development was already fully developed in single-celled protists and evolved nonprogressively from protists to plants and animals.

Fatty acid phytyl ester synthesis in chloroplasts of Arabidopsis

During stress or senescence, thylakoid membranes in chloroplasts are disintegrated, and chlorophyll and galactolipid are broken down, resulting in the accumulation of toxic intermediates, i.e., tetrapyrroles, free phytol, and free fatty acids. Chlorophyll degradation has been studied in detail, but the catabolic pathways for phytol and fatty acids remain unclear. A large proportion of phytol and fatty acids is converted into fatty acid phytyl esters and triacylglycerol during stress or senescence in chloroplasts. We isolated two genes (PHYTYL ESTER SYNTHASE1 [PES1] and PES2) of the esterase/lipase/thioesterase family of acyltransferases from Arabidopsis thaliana that are involved in fatty acid phytyl ester synthesis in chloroplasts. The two proteins are highly expressed during senescence and nitrogen deprivation. Heterologous expression in yeast revealed that PES1 and PES2 have phytyl ester synthesis and diacylglycerol acyltransferase activities. The enzymes show broad substrate specificities and can employ acyl-CoAs, acyl carrier proteins, and galactolipids as acyl donors. Double mutant plants (pes1 pes2) grow normally but show reduced phytyl ester and triacylglycerol accumulation. These results demonstrate that PES1 and PES2 are involved in the deposition of free phytol and free fatty acids in the form of phytyl esters in chloroplasts, a process involved in maintaining the integrity of the photosynthetic membrane during abiotic stress and senescence.

TOPLESS co-repressor interactions and their evolutionary conservation in plants

Large-scale protein-protein interaction studies recently demonstrated that the Arabidopsis TPL/TPR family of transcriptional co-repressors is involved in a broad range of developmental processes. TPL/TPRs predominantly interact with transcription factors that contain repression domain (RD) sequences. Interestingly, RDs reported in the literature are quite diverse in sequence, yet TPL/TPRs interact with proteins containing all of the known motifs. These data lead us to conclude that the TPL/TPRs act as general repressors of gene transcription in plants. To investigate this further, we examined interactions between TPL/TPR proteins encoded by the moss Physcomitrella patens genome and components of the auxin signaling pathway. As in Arabidopsis, moss TPL proteins interact with AUX/IAA and ARF proteins, suggesting that they act in both forms of ARF-mediated transcriptional repression. These data suggest that the involvement of TPL in auxin signaling has been conserved across evolution, since mosses and angiosperms diverged approximately 450 million years ago.

The Arabidopsis tetratricopeptide thioredoxin-like gene family is required for osmotic stress tolerance and male sporogenesis

TETRATRICOPEPTIDE THIOREDOXIN-LIKE (TTL) proteins are characterized by the presence of six tetratricopeptide repeats in conserved positions and a carboxyl-terminal region known as the thioredoxin-like domain with homology to thioredoxins. In Arabidopsis (Arabidopsis thaliana), the TTL gene family is composed by four members, and the founder member, TTL1, is required for osmotic stress tolerance. Analysis of sequenced genomes indicates that TTL genes are specific to land plants. In this study, we report the expression profiles of Arabidopsis TTL genes using data mining and promoter-reporter β-glucuronidase fusions. Our results show that TTL1, TTL3, and TTL4 display ubiquitous expression in normal growing conditions but differential expression patterns in response to osmotic and NaCl stresses. TTL2 shows a very different expression pattern, being specific to pollen grains. Consistent with the expression data, ttl1, ttl3, and ttl4 mutants show reduced root growth under osmotic stress, and the analysis of double and triple mutants indicates that TTL1, TTL3, and TTL4 have partially overlapping yet specific functions in abiotic stress tolerance while TTL2 is involved in male gametophytic transmission.

Comparison of plasma membrane proteomic changes of Arabidopsis suspension-cultured cells (T87 Line) after cold and ABA treatment in association with freezing tolerance development

The plasma membrane (PM) is the primary site of freezing injury in plants. To determine global changes in PM protein profiles in association with freezing tolerance development, proteome analysis of the purified PM of Arabidopsis suspension-cultured cells (T87 line) was conducted with label-free protein quantification technology. Freezing tolerance of Arabidopsis cells at the lag growth phase (8 d old) increased after cold acclimation (CA) or ABA treatment. Proteome analysis assigned 658 proteins in the PM in total, of which 45.3% (298 proteins) were predicted to have transmembrane domains. They were classified into several functional categories, with the primary categories being proteins in transporters, signal transduction, protein destination and storage, and cell structure. After CA, 271 proteins increased and 111 proteins decreased. ABA treatment resulted in 185 increased and 56 decreased proteins. Of these, 139 increased and 49 decreased proteins were identified in common after both CA and ABA treatment. In addition, there were proteins specifically expressed in cold- (132 increased and 62 decreased) or ABA- (46 increased and 7 decreased) treated cells. Collectively, our results clearly show that (i) responses of the PM proteome to CA and ABA treatment overlap substantially but, at the same time, some proteins exhibited different response patterns in each treatment; and (ii) the majority of ABA-responsive proteins are CA-responsive proteins but not vice versa, suggesting complex interactions of CA and ABA signaling pathways in the PM proteome responses.

The phytocyanin gene family in rice (Oryza sativa L.): genome-wide identification, classification and transcriptional analysis

BACKGROUND

Phytocyanins (PCs) are plant-specific blue copper proteins involved in electron transport, and a large number of known PCs are considered to be chimeric arabinogalactan proteins (AGPs). To date there has not been a genome-wide overview of the OsPC gene family. Therefore, as the first step and a useful strategy to elucidate the functions of OsPCs, there is an urgent need for a thorough genome-wide analysis of this gene family.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, a total of 62 OsPC genes were identified through a comprehensive bioinformatics analysis of the rice (Oryza sativa L.) genome. Based on phylogeny and motif constitution, the family of OsPCs was classified into three subclasses: uclacyanin-like proteins (OsUCLs), stellacyanin-like proteins (OsSCLs) and early nodulin-like proteins (OsENODLs). Structure and glycosylation prediction indicated that 46 OsPCs were glycosylphosphatigylinositol-anchored proteins and 38 OsPCs were chimeric AGPs. Gene duplication analysis revealed that chromosomal segment and tandem duplications contributed almost equally to the expansion of this gene family, and duplication events were mostly happened in the OsUCL subfamily. The expression profiles of OsPC genes were analyzed at different stages of vegetative and reproductive development and under abiotic stresses. It revealed that a large number of OsPC genes were abundantly expressed in the various stages of development. Moreover, 17 genes were regulated under the treatments of abiotic stresses.

CONCLUSIONS/SIGNIFICANCE

The genome-wide identification and expression analysis of OsPC genes should facilitate research in this gene family and give new insights toward elucidating their functions in higher plants.

Formation of plant metabolic gene clusters within dynamic chromosomal regions

In bacteria, genes with related functions often are grouped together in operons and are cotranscribed as a single polycistronic mRNA. In eukaryotes, functionally related genes generally are scattered across the genome. Notable exceptions include gene clusters for catabolic pathways in yeast, synthesis of secondary metabolites in filamentous fungi, and the major histocompatibility complex in animals. Until quite recently it was thought that gene clusters in plants were restricted to tandem duplicates (for example, arrays of leucine-rich repeat disease-resistance genes). However, operon-like clusters of coregulated nonhomologous genes are an emerging theme in plant biology, where they may be involved in the synthesis of certain defense compounds. These clusters are unlikely to have arisen by horizontal gene transfer, and the mechanisms behind their formation are poorly understood. Previously in thale cress (Arabidopsis thaliana) we identified an operon-like gene cluster that is required for the synthesis and modification of the triterpene thalianol. Here we characterize a second operon-like triterpene cluster (the marneral cluster) from A. thaliana, compare the features of these two clusters, and investigate the evolutionary events that have led to cluster formation. We conclude that common mechanisms are likely to underlie the assembly and control of operon-like gene clusters in plants.

A 5-formyltetrahydrofolate cycloligase paralog from all domains of life: comparative genomic and experimental evidence for a cryptic role in thiamin metabolism

A paralog (here termed COG0212) of the ATP-dependent folate salvage enzyme 5-formyltetrahydrofolate cycloligase (5-FCL) occurs in all domains of life and, although typically annotated as 5-FCL in pro- and eukaryotic genomes, is of unknown function. COG0212 is similar in overall structure to 5-FCL, particularly in the substrate binding region, and has distant similarity to other kinases. The Arabidopsis thaliana COG0212 protein was shown to be targeted to chloroplasts and to be required for embryo viability. Comparative genomic analysis revealed that a high proportion (19%) of archaeal and bacterial COG0212 genes are clustered on the chromosome with various genes implicated in thiamin metabolism or transport but showed no such association between COG0212 and folate metabolism. Consistent with the bioinformatic evidence for a role in thiamin metabolism, ablating COG0212 in the archaeon Haloferax volcanii caused accumulation of thiamin monophosphate. Biochemical and functional complementation tests of several known and hypothetical thiamin-related activities (involving thiamin, its breakdown products, and their phosphates) were, however, negative. Also consistent with the bioinformatic evidence, the COG0212 proteins from A. thaliana and prokaryote sources lacked 5-FCL activity in vitro and did not complement the growth defect or the characteristic 5-formyltetrahydrofolate accumulation of a 5-FCL-deficient (ΔygfA) Escherichia coli strain. We therefore propose (a) that COG0212 has an unrecognized yet sometimes crucial role in thiamin metabolism, most probably in salvage or detoxification, and (b) that is not a 5-FCL and should no longer be so annotated.

Evolutionary history of the GH3 family of acyl adenylases in rosids

GH3 amino acid conjugases have been identified in many plant and bacterial species. The evolution of GH3 genes in plant species is explored using the sequenced rosids Arabidopsis, papaya, poplar, and grape. Analysis of the sequenced non-rosid eudicots monkey flower and columbine, the monocots maize and rice, as well as spikemoss and moss is included to provide further insight into the origin of GH3 clades. Comparison of co-linear genes in regions surrounding GH3 genes between species helps reconstruct the evolutionary history of the family. Combining analysis of synteny with phylogenetics, gene expression and functional data redefines the Group III GH3 genes, of which AtGH3.12/PBS3, a regulator of stress-induced salicylic acid metabolism and plant defense, is a member. Contrary to previous reports that restrict PBS3 to Arabidopsis and its close relatives, PBS3 syntelogs are identified in poplar, grape, columbine, maize and rice suggesting descent from a common ancestral chromosome dating to before the eudicot/monocot split. In addition, the clade containing PBS3 has undergone a unique expansion in Arabidopsis, with expression patterns for these genes consistent with specialized and evolving stress-responsive functions.

Recruitment of a ribosomal release factor for light- and stress-dependent regulation of petB transcript stability in Arabidopsis chloroplasts

Land plant genomes encode four functional ribosomal peptide chain release factors (Prf) of eubacterial origin, two (PrfA and PrfB homologs) for each endosymbiotic organelle. Formerly, we have shown that the Arabidopsis thaliana chloroplast-localized PrfB homolog, PrfB1, is required not only for termination of translation but also for stabilization of UGA stop codon-containing chloroplast transcripts. A previously undiscovered PrfB-like protein, PrfB3, is localized to the chloroplast stroma in a petB RNA-containing complex and found only in vascular plants. Highly conserved positions of introns unequivocally indicate that PrfB3 arose from a duplication of PrfB1. Notably, PrfB3 is lacking the two most important tripeptide motifs characteristic for all eubacterial and organellar PrfB homologs described so far: the stop codon recognition motif SPF and the catalytic center GGQ for peptidyl-tRNA hydrolysis. Complementation studies, as well as functional and molecular analyses of two allelic mutations in Arabidopsis, both of which lead to a specific deficiency of the cytochrome b₆f complex, revealed that PrfB3 is essentially required for photoautotrophic growth. Plastid transcript, polysome, and translation analyses indicate that PrfB3 has been recruited in vascular plants for light- and stress-dependent regulation of stability of 3' processed petB transcripts to adjust cytochrome b₆ levels.

Speeding cis-trans regulation discovery by phylogenomic analyses coupled with screenings of an arrayed library of Arabidopsis transcription factors

Transcriptional regulation is an important mechanism underlying gene expression and has played a crucial role in evolution. The number, position and interactions between cis-elements and transcription factors (TFs) determine the expression pattern of a gene. To identify functionally relevant cis-elements in gene promoters, a phylogenetic shadowing approach with a lipase gene (LIP1) was used. As a proof of concept, in silico analyses of several Brassicaceae LIP1 promoters identified a highly conserved sequence (LIP1 element) that is sufficient to drive strong expression of a reporter gene in planta. A collection of ca. 1,200 Arabidopsis thaliana TF open reading frames (ORFs) was arrayed in a 96-well format (RR library) and a convenient mating based yeast one hybrid (Y1H) screening procedure was established. We constructed an episomal plasmid (pTUY1H) to clone the LIP1 element and used it as bait for Y1H screenings. A novel interaction with an HD-ZIP (AtML1) TF was identified and abolished by a 2 bp mutation in the LIP1 element. A role of this interaction in transcriptional regulation was confirmed in planta. In addition, we validated our strategy by reproducing the previously reported interaction between a MYB-CC (PHR1) TF, a central regulator of phosphate starvation responses, with a conserved promoter fragment (IPS1 element) containing its cognate binding sequence. Finally, we established that the LIP1 and IPS1 elements were differentially bound by HD-ZIP and MYB-CC family members in agreement with their genetic redundancy in planta. In conclusion, combining in silico analyses of orthologous gene promoters with Y1H screening of the RR library represents a powerful approach to decipher cis- and trans-regulatory codes.

Subunit stoichiometry, evolution, and functional implications of an asymmetric plant plastid ClpP/R protease complex in Arabidopsis

The caseinolytic protease (Clp) protease system has been expanded in plant plastids compared with its prokaryotic progenitors. The plastid Clp core protease consists of five different proteolytic ClpP proteins and four different noncatalytic ClpR proteins, with each present in one or more copies and organized in two heptameric rings. We determined the exact subunit composition and stoichiometry for the intact core and each ring. The chloroplast ClpP/R protease was affinity purified from clpr4 and clpp3 Arabidopsis thaliana null mutants complemented with C-terminal StrepII-tagged versions of CLPR4 and CLPP3, respectively. The subunit stoichiometry was determined by mass spectrometry-based absolute quantification using stable isotope-labeled proteotypic peptides generated from a synthetic gene. One heptameric ring contained ClpP3,4,5,6 in a 1:2:3:1 ratio. The other ring contained ClpP1 and ClpR1,2,3,4 in a 3:1:1:1:1 ratio, resulting in only three catalytic sites. These ClpP1/R1-4 proteins are most closely related to the two subunits of the cyanobacterial P3/R complex and the identical P:R ratio suggests conserved adaptation. Furthermore, the plant-specific C-terminal extensions of the ClpP/R subunits were not proteolytically removed upon assembly, suggesting a regulatory role in Clp chaperone interaction. These results will now allow testing ClpP/R structure-function relationships using rationale design. The quantification workflow we have designed is applicable to other protein complexes.

The plant-specific actin binding protein SCAB1 stabilizes actin filaments and regulates stomatal movement in Arabidopsis

Microfilament dynamics play a critical role in regulating stomatal movement; however, the molecular mechanism underlying this process is not well understood. We report here the identification and characterization of STOMATAL CLOSURE-RELATED ACTIN BINDING PROTEIN1 (SCAB1), an Arabidopsis thaliana actin binding protein. Plants lacking SCAB1 were hypersensitive to drought stress and exhibited reduced abscisic acid-, H(2)O(2)-, and CaCl(2)-regulated stomatal movement. In vitro and in vivo analyses revealed that SCAB1 binds, stabilizes, and bundles actin filaments. SCAB1 shares sequence similarity only with plant proteins and contains a previously undiscovered actin binding domain. During stomatal closure, actin filaments switched from a radial orientation in open stomata to a longitudinal orientation in closed stomata. This switch took longer in scab1 plants than in wild-type plants and was correlated with the delay in stomatal closure seen in scab1 mutants in response to drought stress. Our results suggest that SCAB1 is required for the precise regulation of actin filament reorganization during stomatal closure.