Phylogenetic analysis, structural evolution and functional divergence of the 12-oxo-phytodienoate acid reductase gene family in plants

Oxylipins in moss development and defense

Oxylipins are oxygenated fatty acids that participate in plant development and defense against pathogen infection, insects, and wounding. Initial oxygenation of substrate fatty acids is mainly catalyzed by lipoxygenases (LOXs) and α-dioxygenases but can also take place non-enzymatically by autoxidation or singlet oxygen-dependent reactions. The resulting hydroperoxides are further metabolized by secondary enzymes to produce a large variety of compounds, including the hormone jasmonic acid (JA) and short-chain green leaf volatiles. In flowering plants, which lack arachidonic acid, oxylipins are produced mainly from oxidation of polyunsaturated C18 fatty acids, notably linolenic and linoleic acids. Algae and mosses in addition possess polyunsaturated C20 fatty acids including arachidonic and eicosapentaenoic acids, which can also be oxidized by LOXs and transformed into bioactive compounds. Mosses are phylogenetically placed between unicellular green algae and flowering plants, allowing evolutionary studies of the different oxylipin pathways. During the last years the moss Physcomitrella patens has become an attractive model plant for understanding oxylipin biosynthesis and diversity. In addition to the advantageous evolutionary position, functional studies of the different oxylipin-forming enzymes can be performed in this moss by targeted gene disruption or single point mutations by means of homologous recombination. Biochemical characterization of several oxylipin-producing enzymes and oxylipin profiling in P. patens reveal the presence of a wider range of oxylipins compared to flowering plants, including C18 as well as C20-derived oxylipins. Surprisingly, one of the most active oxylipins in plants, JA, is not synthesized in this moss. In this review, we present an overview of oxylipins produced in mosses and discuss the current knowledge related to the involvement of oxylipin-producing enzymes and their products in moss development and defense.

Cyclic nucleotide-gated ion channel gene family in rice, identification, characterization and experimental analysis of expression response to plant hormones, biotic and abiotic stresses

BACKGROUND

Cyclic nucleotide-gated channels (CNGCs) are Ca2+-permeable cation transport channels, which are present in both animal and plant systems. They have been implicated in the uptake of both essential and toxic cations, Ca2+ signaling, pathogen defense, and thermotolerance in plants. To date there has not been a genome-wide overview of the CNGC gene family in any economically important crop, including rice (Oryza sativa L.). There is an urgent need for a thorough genome-wide analysis and experimental verification of this gene family in rice.

RESULTS

In this study, a total of 16 full length rice CNGC genes distributed on chromosomes 1-6, 9 and 12, were identified by employing comprehensive bioinformatics analyses. Based on phylogeny, the family of OsCNGCs was classified into four major groups (I-IV) and two sub-groups (IV-A and IV- B). Likewise, the CNGCs from all plant lineages clustered into four groups (I-IV), where group II was conserved in all land plants. Gene duplication analysis revealed that both chromosomal segmentation (OsCNGC1 and 2, 10 and 11, 15 and 16) and tandem duplications (OsCNGC1 and 2) significantly contributed to the expansion of this gene family. Motif composition and protein sequence analysis revealed that the CNGC specific domain "cyclic nucleotide-binding domain (CNBD)" comprises a "phosphate binding cassette" (PBC) and a "hinge" region that is highly conserved among the OsCNGCs. In addition, OsCNGC proteins also contain various other functional motifs and post-translational modification sites. We successively built a stringent motif: (LI-X(2)-[GS]-X-[FV]-X-G-[1]-ELL-X-W-X(12,22)-SA-X(2)-T-X(7)-[EQ]-AF-X-L) that recognizes the rice CNGCs specifically. Prediction of cis-acting regulatory elements in 5' upstream sequences and expression analyses through quantitative qPCR demonstrated that OsCNGC genes were highly responsive to multiple stimuli including hormonal (abscisic acid, indoleacetic acid, kinetin and ethylene), biotic (Pseudomonas fuscovaginae and Xanthomonas oryzae pv. oryzae) and abiotic (cold) stress.

CONCLUSIONS

There are 16 CNGC genes in rice, which were probably expanded through chromosomal segmentation and tandem duplications and comprise a PBC and a "hinge" region in the CNBD domain, featured by a stringent motif. The various cis-acting regulatory elements in the upstream sequences may be responsible for responding to multiple stimuli, including hormonal, biotic and abiotic stresses.

Genome-wide analysis of the NADK gene family in plants

Background

NAD(H) kinase (NADK) is the key enzyme that catalyzes de novo synthesis of NADP(H) from NAD(H) for NADP(H)-based metabolic pathways. In plants, NADKs form functional subfamilies. Studies of these families in Arabidopsis thaliana indicate that they have undergone considerable evolutionary selection; however, the detailed evolutionary history and functions of the various NADKs in plants are not clearly understood.

Principal Findings

We performed a comparative genomic analysis that identified 74 NADK gene homologs from 24 species representing the eight major plant lineages within the supergroup Plantae: glaucophytes, rhodophytes, chlorophytes, bryophytes, lycophytes, gymnosperms, monocots and eudicots. Phylogenetic and structural analysis classified these NADK genes into four well-conserved subfamilies with considerable variety in the domain organization and gene structure among subfamily members. In addition to the typical NAD_kinase domain, additional domains, such as adenylate kinase, dual-specificity phosphatase, and protein tyrosine phosphatase catalytic domains, were found in subfamily II. Interestingly, NADKs in subfamily III exhibited low sequence similarity (∼30%) in the kinase domain within the subfamily and with the other subfamilies. These observations suggest that gene fusion and exon shuffling may have occurred after gene duplication, leading to specific domain organization seen in subfamilies II and III, respectively. Further analysis of the exon/intron structures showed that single intron loss and gain had occurred, yielding the diversified gene structures, during the process of structural evolution of NADK family genes. Finally, both available global microarray data analysis and qRT-RCR experiments revealed that the NADK genes in Arabidopsis and Oryza sativa show different expression patterns in different developmental stages and under several different abiotic/biotic stresses and hormone treatments, underscoring the functional diversity and functional divergence of the NADK family in plants.

Conclusions

These findings will facilitate further studies of the NADK family and provide valuable information for functional validation of this family in plants.

Cloning and molecular characterization of a flavin-dependent oxidoreductase gene from barley

Oxophytodienoate reductases (OPRs) are a small group of flavin-dependent oxidoreductases in plants. In this study, a new member of the OPR gene family (HvOPR2) was cloned from barley (Hordeum vulgare L.) using reverse transcription polymerase chain reaction (RT-PCR). The full-length cDNA of HvOPR2 was 1,206 bp with an open reading frame of 1,101 bp, encoding a 366 amino acids long polypeptide with a predicted molecular weight of 40.52 and a theoretical isoelectric point of 6.21. The corresponding genomic clone of HvOPR2 was isolated using the PCR amplification technique and was found to consist of five exons and four introns. Bioinformatic analysis revealed that the deduced HvOPR2 has a considerable homology with other plant OPRs and possessed the flavin oxidoreductase/NADH oxidase substrate-binding domain. Phylogenetic analysis showed that HvOPR2 codes for the OPR of subgroup I, which contains enzymes that are not required for jasmonic acid biosynthesis. Time-course transcriptional profiling of HvOPR2 was analyzed in response to a variety of abiotic stresses and hormonal treatments by semi-quantitative RT-PCR. The HvOPR2 gene was induced in response to drought, hydrogen peroxide, and wounding. Moreover, the corresponding mRNA transcripts were increased in response to jasmonic acid and salicylic acid, but not in response to abscisic acid. These results strongly suggested a role for HvOPR2 in barley defense/response to abiotic stresses and signaling molecules.

Carbon-carbon double-bond reductases in nature

Reduction of C = C bonds by reductases, found in a variety of microorganisms (e.g. yeasts, bacteria, and lower fungi), animals, and plants has applications in the production of metabolites that include pharmacologically active drugs and other chemicals. Therefore, the reductase enzymes that mediate this transformation have become important therapeutic targets and biotechnological tools. These reductases are broad-spectrum, in that, they can act on isolation/conjugation C = C-bond compounds, α,β-unsaturated carbonyl compounds, carboxylic acids, acid derivatives, and nitro compounds. In addition, several mutations in the reductase gene have been identified, some associated with diseases. Several of these reductases have been cloned and/or purified, and studies to further characterize them and determine their structure in order to identify potential industrial biocatalysts are still in progress. In this study, crucial reductases for bioreduction of C = C bonds have been reviewed with emphasis on their principal substrates and effective inhibitors, their distribution, genetic polymorphisms, and implications in human disease and treatment.

Adaptive functional divergence of the warm temperature acclimation-related protein (WAP65) in fishes and the ortholog hemopexin (HPX) in mammals

Gene duplication is an important mechanism that leads to genetic novelty. Different, nonexclusive processes are likely involved, and many adaptive and nonadaptive events may contribute to the maintenance of duplicated genes. In some teleosts, a duplicate copy of the mammalian ortholog Hemopexin (HPX) is present, known as the warm temperature acclimation-related protein (WAP65). Both WAP65 and HPX have been associated with iron homeostasis due to the affinity to bind the toxic-free heme circulating in the blood stream. We have assessed the evolutionary dynamics of WAP65 and HPX genes to understand the adaptive role of positive selection at both nucleotide and amino acid level. Our results showed an asymmetrical evolution between the paralogs WAP65-1 and WAP65-2 after duplication with a slight acceleration of the evolutionary rate in WAP65-1, but not in WAP65-2, and few sites contributing to the functional distinction between the paralogs, whereas the majority of the protein remained under negative selection or relaxed negative selection. WAP65-1 is functionally more distinct from the ancestral protein function than WAP65-2. HPX is phylogenetically closer to WAP65-2 but even so functional divergence was detected between both proteins. In addition, HPX showed a fast rate of evolution when compared with both WAP65-1 and WAP65-2 genes. The assessed 3-dimensional (3-D) structure of WAP65-1 and WAP65-2 suggests that the functional differences detected are not causing noticeable structural changes in these proteins. However, such subtle changes between WAP65 paralogs may be important to understand the differential gene retention of both copies in 20 out of 30 teleosts species studied.

Genome-wide identification, evolutionary and expression analysis of the aspartic protease gene superfamily in grape

Background

Aspartic proteases (APs) are a large family of proteolytic enzymes found in almost all organisms. In plants, they are involved in many biological processes, such as senescence, stress responses, programmed cell death, and reproduction. Prior to the present study, no grape AP gene(s) had been reported, and their research on woody species was very limited.

Results

In this study, a total of 50 AP genes (VvAP) were identified in the grape genome, among which 30 contained the complete ASP domain. Synteny analysis within grape indicated that segmental and tandem duplication events contributed to the expansion of the grape AP family. Additional analysis between grape and Arabidopsis demonstrated that several grape AP genes were found in the corresponding syntenic blocks of Arabidopsis, suggesting that these genes arose before the divergence of grape and Arabidopsis. Phylogenetic relationships of the 30 VvAPs with the complete ASP domain and their Arabidopsis orthologs, as well as their gene and protein features were analyzed and their cellular localization was predicted. Moreover, expression profiles of VvAP genes in six different tissues were determined, and their transcript abundance under various stresses and hormone treatments were measured. Twenty-seven VvAP genes were expressed in at least one of the six tissues examined; nineteen VvAPs responded to at least one abiotic stress, 12 VvAPs responded to powdery mildew infection, and most of the VvAPs responded to SA and ABA treatments. Furthermore, integrated synteny and phylogenetic analysis identified orthologous AP genes between grape and Arabidopsis, providing a unique starting point for investigating the function of grape AP genes.

Conclusions

The genome-wide identification, evolutionary and expression analyses of grape AP genes provide a framework for future analysis of AP genes in defining their roles during stress response. Integrated synteny and phylogenetic analyses provide novel insight into the functions of less well-studied genes using information from their better understood orthologs.

Molecular evolution and functional divergence of soluble starch synthase genes in cassava (manihot esculenta crantz)

Soluble starch synthases (SSs) are major enzymes involved in starch biosynthesis in plants. Cassava starch has many remarkable characteristics, which should be influenced by the evolution of SS genes in this starchy root crop. In this work, we performed a comprehensive phylogenetic and evolutionary analysis of the soluble starch synthases in cassava. Genome-wide identification showed that there are 9 genes encoding soluble starch synthases in cassava. All of the soluble starch synthases encoded by these genes contain both Glyco_transf_5 and Glycos_transf_1 domains, and a correlation analysis showed evidence of coevolution between these 2 domains in cassava SS genes. The SS genes in land plants can be divided into 6 subfamilies that were formed before the origin of seed plants, and species-specific expansion has contributed to the evolution of this family in cassava. A functional divergence analysis for this family provided statistical evidence for shifted evolutionary rates between the subfamilies of land plant soluble starch synthases. Although the main selective pressure acting on land plant SS genes was purifying selection, our results also revealed that point mutation with positive selection contributed to the evolution of 2 SS genes in cassava. The remarkable cassava starch characteristics might be the result of both the duplication and adaptive selection of SS genes.

Activation of Defense Mechanisms against Pathogens in Mosses and Flowering Plants

During evolution, plants have developed mechanisms to cope with and adapt to different types of stress, including microbial infection. Once the stress is sensed, signaling pathways are activated, leading to the induced expression of genes with different roles in defense. Mosses (Bryophytes) are non-vascular plants that diverged from flowering plants more than 450 million years ago, allowing comparative studies of the evolution of defense-related genes and defensive metabolites produced after microbial infection. The ancestral position among land plants, the sequenced genome and the feasibility of generating targeted knock-out mutants by homologous recombination has made the moss Physcomitrella patens an attractive model to perform functional studies of plant genes involved in stress responses. This paper reviews the current knowledge of inducible defense mechanisms in P. patens and compares them to those activated in flowering plants after pathogen assault, including the reinforcement of the cell wall, ROS production, programmed cell death, activation of defense genes and synthesis of secondary metabolites and defense hormones. The knowledge generated in P. patens together with comparative studies in flowering plants will help to identify key components in plant defense responses and to design novel strategies to enhance resistance to biotic stress.

Role of plant peroxisomes in the production of jasmonic acid-based signals

Jasmonates are a family of oxylipins derived from linolenic acid that control plant responses to biotic and abiotic stress factors and also regulate plant growth and development. Jasmonic acid (JA) is synthesized through the octadecanoid pathway that involves the translocation of lipid intermediates from the chloroplast membranes to the cytoplasm and later on into peroxisomes. The peroxisomal steps of the pathway involve the reduction of cis-(+)-12-oxophytodienoic acid (12-OPDA) and dinor-OPDA, which are the final products of the choroplastic phase of the biosynthetic pathway acting on 18:3 and 16:3 fatty acids, respectively. Further shortening of the carbon side-chain by successive rounds of β-oxidation reactions are required to complete JA biosynthesis. After peroxisomal reactions are completed, (+)-7-iso-JA is synthesized and then transported to the cytoplasm where is conjugated to the amino acid isoleucine to form the bioactive form of the hormone (+)-7-iso-JA-Ile (JA-Ile). Further regulatory activity of JA-Ile triggering gene activation in the jasmonate-dependent signaling cascades is exerted through a process mediated by the perception via the E3 ubiquitin ligase COI1 and further ligand-activated interaction with the family of JAZ repressor proteins. Upon interaction, JAZ are ubiquitinated and degraded by the proteasome, thus releasing transcription factors such as MYC2 from repression and allowing the activation of JA-responsive genes.

Jasmonate: A decision maker between cell death and acclimation in the response of plants to singlet oxygen

Under stress conditions that bring about excessive absorption of light energy in the chloroplasts, the formation of singlet oxygen ( (1)O2) can be strongly enhanced, triggering programmed cell death. However, the (1)O2 signaling pathway can also lead to acclimation to photooxidative stress, when (1)O2 is produced in relatively low amounts. This acclimatory response is associated with a strong downregulation of the jasmonate biosynthesis pathway and the maintenance of low jasmonate levels, even under high light stress conditions that normally induce jasmonate synthesis. These findings suggest a central role for this phytohormone in the orientation of the (1)O2 signaling pathway toward cell death or acclimation. This conclusion is confirmed here in an Arabidopsis double mutant obtained by crossing the (1)O2-overproducing mutant ch1 and the jasmonate-deficient mutant dde2. This double mutant was found to be constitutively resistant to (1)O2 stress and to display a strongly stimulated growth rate compared with the single ch1 mutant. However, the involvement of other phytohormones, such as ethylene, cannot be excluded.

PGDD: a database of gene and genome duplication in plants

Genome duplication (GD) has permanently shaped the architecture and function of many higher eukaryotic genomes. The angiosperms (flowering plants) are outstanding models in which to elucidate consequences of GD for higher eukaryotes, owing to their propensity for chromosomal duplication or even triplication in a few cases. Duplicated genome structures often require both intra- and inter-genome alignments to unravel their evolutionary history, also providing the means to deduce both obvious and otherwise-cryptic orthology, paralogy and other relationships among genes. The burgeoning sets of angiosperm genome sequences provide the foundation for a host of investigations into the functional and evolutionary consequences of gene and GD. To provide genome alignments from a single resource based on uniform standards that have been validated by empirical studies, we built the Plant Genome Duplication Database (PGDD; freely available at http://chibba.agtec.uga.edu/duplication/), a web service providing synteny information in terms of colinearity between chromosomes. At present, PGDD contains data for 26 plants including bryophytes and chlorophyta, as well as angiosperms with draft genome sequences. In addition to the inclusion of new genomes as they become available, we are preparing new functions to enhance PGDD.

Functional and evolutionary analysis of DXL1, a non-essential gene encoding a 1-deoxy-D-xylulose 5-phosphate synthase like protein in Arabidopsis thaliana

The synthesis of 1-deoxy-D-xylulose 5-phosphate (DXP), catalyzed by the enzyme DXP synthase (DXS), represents a key regulatory step of the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway for isoprenoid biosynthesis. In plants DXS is encoded by small multigene families that can be classified into, at least, three specialized subfamilies. Arabidopsis thaliana contains three genes encoding proteins with similarity to DXS, including the well-known DXS1/CLA1 gene, which clusters within subfamily I. The remaining proteins, initially named DXS2 and DXS3, have not yet been characterized. Here we report the expression and functional analysis of A. thaliana DXS2. Unexpectedly, the expression of DXS2 failed to rescue Escherichia coli and A. thaliana mutants defective in DXS activity. Coherently, we found that DXS activity was negligible in vitro, being renamed as DXL1 following recent nomenclature recommendation. DXL1 is targeted to plastids as DXS1, but shows a distinct expression pattern. The phenotypic analysis of a DXL1 defective mutant revealed that the function of the encoded protein is not essential for growth and development. Evolutionary analyses indicated that DXL1 emerged from DXS1 through a recent duplication apparently specific of the Brassicaceae lineage. Divergent selective constraints would have affected a significant fraction of sites after diversification of the paralogues. Furthermore, amino acids subjected to divergent selection and likely critical for functional divergence through the acquisition of a novel, although not yet known, biochemical function, were identified. Our results provide with the first evidences of functional specialization at both the regulatory and biochemical level within the plant DXS family.

Physcomitrella patens activates reinforcement of the cell wall, programmed cell death and accumulation of evolutionary conserved defence signals, such as salicylic acid and 12-oxo-phytodienoic acid, but not jasmonic acid, upon Botrytis cinerea infection

The moss Physcomitrella patens is an evolutionarily basal model system suitable for the analysis of plant defence responses activated after pathogen assault. Upon infection with the necrotroph Botrytis cinerea, several defence mechanisms are induced in P. patens, including the fortification of the plant cell wall by the incorporation of phenolic compounds and the induced expression of related genes. Botrytis cinerea infection also activates the accumulation of reactive oxygen species and cell death with hallmarks of programmed cell death in moss tissues. Salicylic acid (SA) levels also increase after fungal infection, and treatment with SA enhances transcript accumulation of the defence gene phenylalanine ammonia-lyase (PAL) in P. patens colonies. The expression levels of the genes involved in 12-oxo-phytodienoic acid (OPDA) synthesis, including lipoxygenase (LOX) and allene oxide synthase (AOS), increase in P. patens gametophytes after pathogen assault, together with a rise in free linolenic acid and OPDA concentrations. However, jasmonic acid (JA) could not be detected in healthy or infected tissues of this plant. Our results suggest that, although conserved defence signals, such as SA and OPDA, are synthesized and are probably involved in the defence response of P. patens against B. cinerea infection, JA production appears to be missing. Interestingly, P. patens responds to OPDA and methyl jasmonate by reducing moss colony growth and rhizoid length, suggesting that jasmonate perception is present in mosses. Thus, P. patens can provide clues with regard to the evolution of different defence pathways in plants, including signalling and perception of OPDA and jasmonates in nonflowering and flowering plants.

Physcomitrella patens activates reinforcement of the cell wall, programmed cell death and accumulation of evolutionary conserved defence signals, such as salicylic acid and 12-oxo-phytodienoic acid, but not jasmonic acid, upon Botrytis cinerea infection

The moss Physcomitrella patens is an evolutionarily basal model system suitable for the analysis of plant defence responses activated after pathogen assault. Upon infection with the necrotroph Botrytis cinerea, several defence mechanisms are induced in P. patens, including the fortification of the plant cell wall by the incorporation of phenolic compounds and the induced expression of related genes. Botrytis cinerea infection also activates the accumulation of reactive oxygen species and cell death with hallmarks of programmed cell death in moss tissues. Salicylic acid (SA) levels also increase after fungal infection, and treatment with SA enhances transcript accumulation of the defence gene phenylalanine ammonia-lyase (PAL) in P. patens colonies. The expression levels of the genes involved in 12-oxo-phytodienoic acid (OPDA) synthesis, including lipoxygenase (LOX) and allene oxide synthase (AOS), increase in P. patens gametophytes after pathogen assault, together with a rise in free linolenic acid and OPDA concentrations. However, jasmonic acid (JA) could not be detected in healthy or infected tissues of this plant. Our results suggest that, although conserved defence signals, such as SA and OPDA, are synthesized and are probably involved in the defence response of P. patens against B. cinerea infection, JA production appears to be missing. Interestingly, P. patens responds to OPDA and methyl jasmonate by reducing moss colony growth and rhizoid length, suggesting that jasmonate perception is present in mosses. Thus, P. patens can provide clues with regard to the evolution of different defence pathways in plants, including signalling and perception of OPDA and jasmonates in nonflowering and flowering plants.

Genome-wide identification and analysis of the TIFY gene family in grape

BACKGROUND

The TIFY gene family constitutes a plant-specific group of genes with a broad range of functions. This family encodes four subfamilies of proteins, including ZML, TIFY, PPD and JASMONATE ZIM-Domain (JAZ) proteins. JAZ proteins are targets of the SCF(COI1) complex, and function as negative regulators in the JA signaling pathway. Recently, it has been reported in both Arabidopsis and rice that TIFY genes, and especially JAZ genes, may be involved in plant defense against insect feeding, wounding, pathogens and abiotic stresses. Nonetheless, knowledge concerning the specific expression patterns and evolutionary history of plant TIFY family members is limited, especially in a woody species such as grape.

METHODOLOGY/PRINCIPAL FINDINGS

A total of two TIFY, four ZML, two PPD and 11 JAZ genes were identified in the Vitis vinifera genome. Phylogenetic analysis of TIFY protein sequences from grape, Arabidopsis and rice indicated that the grape TIFY proteins are more closely related to those of Arabidopsis than those of rice. Both segmental and tandem duplication events have been major contributors to the expansion of the grape TIFY family. In addition, synteny analysis between grape and Arabidopsis demonstrated that homologues of several grape TIFY genes were found in the corresponding syntenic blocks of Arabidopsis, suggesting that these genes arose before the divergence of lineages that led to grape and Arabidopsis. Analyses of microarray and quantitative real-time RT-PCR expression data revealed that grape TIFY genes are not a major player in the defense against biotrophic pathogens or viruses. However, many of these genes were responsive to JA and ABA, but not SA or ET.

CONCLUSION

The genome-wide identification, evolutionary and expression analyses of grape TIFY genes should facilitate further research of this gene family and provide new insights regarding their evolutionary history and regulatory control.

The plant proteome folding project: structure and positive selection in plant protein families

Despite its importance, relatively little is known about the relationship between the structure, function, and evolution of proteins, particularly in land plant species. We have developed a database with predicted protein domains for five plant proteomes (http://pfp.bio.nyu.edu) and used both protein structural fold recognition and de novo Rosetta-based protein structure prediction to predict protein structure for Arabidopsis and rice proteins. Based on sequence similarity, we have identified ∼15,000 orthologous/paralogous protein family clusters among these species and used codon-based models to predict positive selection in protein evolution within 175 of these sequence clusters. Our results show that codons that display positive selection appear to be less frequent in helical and strand regions and are overrepresented in amino acid residues that are associated with a change in protein secondary structure. Like in other organisms, disordered protein regions also appear to have more selected sites. Structural information provides new functional insights into specific plant proteins and allows us to map positively selected amino acid sites onto protein structures and view these sites in a structural and functional context.

MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity

MCScan is an algorithm able to scan multiple genomes or subgenomes in order to identify putative homologous chromosomal regions, and align these regions using genes as anchors. The MCScanX toolkit implements an adjusted MCScan algorithm for detection of synteny and collinearity that extends the original software by incorporating 14 utility programs for visualization of results and additional downstream analyses. Applications of MCScanX to several sequenced plant genomes and gene families are shown as examples. MCScanX can be used to effectively analyze chromosome structural changes, and reveal the history of gene family expansions that might contribute to the adaptation of lineages and taxa. An integrated view of various modes of gene duplication can supplement the traditional gene tree analysis in specific families. The source code and documentation of MCScanX are freely available at http://chibba.pgml.uga.edu/mcscan2/.

Analyses of the oligopeptide transporter gene family in poplar and grape

Background

Oligopeptide transporters (OPTs) are a group of membrane-localized proteins that have a broad range of substrate transport capabilities and that are thought to contribute to many biological processes. The OPT proteins belong to a small gene family in plants, which includes about 25 members in Arabidopsis and rice. However, no comprehensive study incorporating phylogeny, chromosomal location, gene structure, expression profiling, functional divergence and selective pressure analysis has been reported thus far for Populus and Vitis.

Results

In the present study, a comprehensive analysis of the OPT gene family in Populus (P. trichocarpa) and Vitis (V. vinifera) was performed. A total of 20 and 18 full-length OPT genes have been identified in Populus and Vitis, respectively. Phylogenetic analyses indicate that these OPT genes consist of two classes that can be further subdivided into 11 groups. Gene structures are considerably conserved among the groups. The distribution of OPT genes was found to be non-random across chromosomes. A high proportion of the genes are preferentially clustered, indicating that tandem duplications may have contributed significantly to the expansion of the OPT gene family. Expression patterns based on our analyses of microarray data suggest that many OPT genes may be important in stress response and functional development of plants. Further analyses of functional divergence and adaptive evolution show that, while purifying selection may have been the main force driving the evolution of the OPTs, some of critical sites responsible for the functional divergence may have been under positive selection.

Conclusions

Overall, the data obtained from our investigation contribute to a better understanding of the complexity of the Populus and Vitis OPT gene family and of the function and evolution of the OPT gene family in higher plants.

Comparative characterization, expression pattern and function analysis of the 12-oxo-phytodienoic acid reductase gene family in rice

The 12-oxo-phytodienoic acid reductases (OPRs) belong to the old yellow enzyme family of flavoenzymes and form multiple subfamilies in angiosperm plants. In our previous study, a comparative genomic analysis showed that five OPR subfamilies (subs. I-V) occur in monocots, and two subfamilies (subs. I and II) in dicots. Here, a comparative study of five OsOPR genes, representing five subfamilies (I-V) in rice, was performed to provide insights into OPR biochemical properties and physiological importance. Comparative analysis of the three-dimensional structure by homology modeling indicated all five OsOPR proteins contained a highly conserved backbone with (α/β)(8)-barrels, while two middle variable regions (MVR i and ii) were also detected and defined. Analysis of enzymatic characteristics revealed that all five OsOPR fusion proteins exhibit distinct substrate specificity. Different catalytic activity was observed using racemic OPDA and trans-2-hexen-1-al as substrates, suggesting OsOPR family genes participate in two main branches of the octadecanoid pathway, including the allene oxide synthase and hydroperoxide lyase pathways which regulate various developmental processes and/or defense responses. The transcript profiles of five OsOPR genes exhibited strong tissue-specific and inducible expression patterns under abiotic stress, hormones and plant wounding treatments. Furthermore, the transcriptions of OsOPR04-1 (OsOPR11) and OsOPR08-1 (OsOPR7), representing subs. I and II, respectively, were observed in all six selected tissues and with all above-stress treatments. This suggests that these two subfamilies play an important role during different developmental stages and in response to stresses; while the expressions of OsOPR06-1 (OsOPR6), OsOPR01-1 (OsOPR10) and OsOPR02-1 (OsOPR8), representing subs. III, IV and V respectively, were strongly up-regulated with abscisic acid (ABA) and indoleacetic acid (IAA) treatments in roots, suggesting these three subfamilies play an important role in responding to hormones especially ABA and IAA signals in roots.

Origin and evolutionary analysis of the plant-specific TIFY transcription factor family

A substantial number of transcription factor families have been investigated from all kingdoms of life, but a particular class of plant-specific TIFY transcription factors, characterized by a highly conserved TIFY domain, lacks a systemic analysis of its origin and evolutionary relationships among different plant species. After exhaustive genome-wide searches against 14 genomes, TIFY transcription factors were identified and classified into four subfamilies TIFY, PPD, JAZ and ZML according to their different domain architectures. Results show that the TIFY domain of the ZML subfamily possesses a core "TLS[F/Y]XG" motif rather than the "TIFYXG" motif that is dominant in the other three subfamilies. A comprehensive survey of the TIFY family allowed us to discover a new group within the JAZ subfamily and to identify several novel conserved motifs via phylogenetic analysis. Evolutional analysis indicates that whole genome duplication and tandem duplication contributed to the expansion of the TIFY family in plants.

Phylogenetic analysis and evolution of aromatic amino acid hydroxylase

A study was performed to investigate the phylogenetic relationship among AAAH members and to statistically evaluate sequence conservation and functional divergence. In total, 161 genes were identified from 103 species. Phylogenetic analysis showed that well-conserved subfamilies exist. Exon-intron structure analysis showed that the gene structures of AAAH were highly conserved across some different lineage species, while some species-specific introns were also found. The dynamic distribution of ACT domain suggested one gene fusion event has occurred in eukaryota. Significant functional divergence was found between some subgroups. Analysis of the site-specific profiles revealed critical amino acid residues for functional divergence. This study highlights the molecular evolution of this family and may provide a starting point for further experimental verifications.

The moss Physcomitrella patens contains cyclopentenones but no jasmonates: mutations in allene oxide cyclase lead to reduced fertility and altered sporophyte morphology

• Two cDNAs encoding allene oxide cyclases (PpAOC1, PpAOC2), key enzymes in the formation of jasmonic acid (JA) and its precursor (9S,13S)-12-oxo-phytodienoic acid (cis-(+)-OPDA), were isolated from the moss Physcomitrella patens. • Recombinant PpAOC1 and PpAOC2 show substrate specificity against the allene oxide derived from 13-hydroperoxy linolenic acid (13-HPOTE); PpAOC2 also shows substrate specificity against the allene oxide derived from 12-hydroperoxy arachidonic acid (12-HPETE). • In protonema and gametophores the occurrence of cis-(+)-OPDA, but neither JA nor the isoleucine conjugate of JA nor that of cis-(+)-OPDA was detected. • Targeted knockout mutants for PpAOC1 and for PpAOC2 were generated, while double mutants could not be obtained. The ΔPpAOC1 and ΔPpAOC2 mutants showed reduced fertility, aberrant sporophyte morphology and interrupted sporogenesis.

Cloning, functional expression and phylogenetic analysis of plant sterol 24C-methyltransferases involved in sitosterol biosynthesis

Sterol 24C-methyltransferases (SMTs) constitute a group of sequence-related proteins that catalyze the distinct patterns of 24-alkyl sterols that occur throughout nature. Two SMT cDNAs (SMT2-1 and SMT2-2) were cloned by homology based PCR methods from young leaves of Glycine max (soybean) and the corresponding enzymes were expressed functionally in Escherichia coli. The full-length cDNA for SMT2-1 and SMT2-2 have open reading frames of 1086 bp and 1092 bp, respectively, and encode proteins of 361 and 363 residues with a calculated molecular mass of 40.3 and 40.4 kDa, respectively. The substrate preference of the two isoforms was similar yet they differed from SMT1; kinetically SMT2-1 and SMT2-2 generated k(cat) values for the optimal substrate 24(28)methylene lophenol of 0.8 min(-1) and 1.34 min(-1), respectively, compared to the activity of SMT1 that generated a k(cat) for the optimal substrate cycloartenol of 0.6 min(-1). SMT2-2 was purified to homogeneity and the subunit organization shown to be tetrameric in similar fashion to other cloned SMTs. Analysis of the accumulated products catalyzed by the recombinant enzymes demonstrated that soybean SMT2-1 and SMT2-2 operate transalkylation activities analogous to the soybean plant SMT1. Metabolite analyses correlated with transcript profiling of the three SMT isoforms during soybean maturation clearly demonstrated that SMT isoform expression determines specific C24-methyl to C24-ethyl ratios to flowering whereas with seed development there is a disconnection such that the SMT transcript levels decrease against an increase in sterol content; generally SMT2-2 is expressed more than SMT2-1 or SMT1. These observations suggest that the genes that encode SMT1 and SMT2 in sitosterol biosynthesis may have undergone divergent evolution. In support of this proposition, the genomic organization for SMT1 of fungi and protozoa align very closely with one another and to those of the plant SMT2; both sets of SMTs lack introns. Unexpectedly, the SMT1 from Glycine max and other embryophytes of diverse origin possess disparate intron-exon characteristics that can be shown relates back to the algae. Our results suggest that the order of SMT1 appearing before SMT2 in phytosterol synthesis arose recently in plant evolution in response to duplication of a more primitive SMT gene likely to have been bifunctional and catalytically promiscuous.