Genome-wide identification and analysis of membrane-bound O-acyltransferase (MBOAT) gene family in plants

Membrane bound O-acyl transferase (MBOAT) family is composed of gene members encoding a variety of acyltransferase enzymes, which play important roles in plant acyl lipid metabolism. Here, we present the first genome-enabled identification and analysis of MBOAT gene models in plants. In total, we identified 136 plant MBOAT sequences from 14 plant species with complete genomes. Phylogenetic relationship analyses suggested the plant MBOAT gene models fell into four major groups, two of which likely encode enzymes of diacylglycerol acyltransferase 1 (DGAT1) and lysophospholipid acyltransferase (LPLAT), respectively, with one–three copies of paralogs present in each of the most plant species. A group of gene sequences, which are homologous to Saccharomyces cerevisiae glycerol uptake proteins (GUP), was identified in plants; copy numbers were conserved, with only one copy represented in each of the most plant species; analyses showed that residues essential for acyltransferases were more prone to be conserved than vertebrate orthologs. Among four groups, one was inferred to emerge in land plants and experience a rapid expansion in genomes of angiosperms, which suggested their important roles in adaptation of plants in lands. Sequence and phylogeny analyses indicated that genes in all four groups encode enzymes with acyltransferases. Comprehensive sequence identification of MBOAT family members and investigation into classification provide a complete picture of the MBOAT gene family in plants, and could shed light into enzymatic functions of different MBOAT genes in plants.     

Phylogenetic analyses provide the first insights into the evolution of OVATE family proteins in land plants

Background and Aims

The OVATE gene encodes a nuclear-localized regulatory protein belonging to a distinct family of plant-specific proteins known as the OVATE family proteins (OFPs). OVATE was first identified as a key regulator of fruit shape in tomato, with nonsense mutants displaying pear-shaped fruits. However, the role of OFPs in plant development has been poorly characterized.

Methods

Public databases were searched and a total of 265 putative OVATE protein sequences were identified from 13 sequenced plant genomes that represent the major evolutionary lineages of land plants. A phylogenetic analysis was conducted based on the alignment of the conserved OVATE domain from these 13 selected plant genomes. The expression patterns of tomato SlOFP genes were analysed via quantitative real-time PCR. The pattern of OVATE gene duplication resulting in the expansion of the gene family was determined in arabidopsis, rice and tomato.

Key Results

Genes for OFPs were found to be present in all the sampled land plant genomes, including the early-diverged lineages, mosses and lycophytes. Phylogenetic analysis based on the amino acid sequences of the conserved OVATE domain defined 11 sub-groups of OFPs in angiosperms. Different evolutionary mechanisms are proposed for OVATE family evolution, namely conserved evolution and divergent expansion. Characterization of the AtOFP family in arabidopsis, the OsOFP family in rice and the SlOFP family in tomato provided further details regarding the evolutionary framework and revealed a major contribution of tandem and segmental duplications towards expansion of the OVATE gene family.

Conclusions

This first genome-wide survey on OFPs provides new insights into the evolution of the OVATE protein family and establishes a solid base for future functional genomics studies on this important but poorly characterized regulatory protein family in plants.     

Computational identification and phylogenetic analysis of the oil-body structural proteins,oleosin and caleosin,in castor bean and flax

Oil bodies (OBs) are the intracellular particles derived from oilseeds. These OBs store lipids as a carbon resource, and have been exploited for a variety of industrial applications including biofuels. Oleosin and caleosin are the common OB structural proteins which are enabling biotechnological enhancement of oil content and OB-based pharmaceutical formations via stabilizing OBs. Although the draft whole genome sequence information for Ricinus communis L. (castor bean) and Linum usitatissimum L. (flax), important oil seed plants, is available in public database, OB-structural proteins in these plants are poorly indentified. Therefore, in this study, we performed a comprehensive bioinformatic analysis including analysis of the genome sequence, conserved domains and phylogenetic relationships to identify OB structural proteins in castor bean and flax genomes. Using comprehensive analysis, we have identified 6 and 15 OB-structural proteins from castor bean and flax, respectively. A complete overview of this gene family in castor bean and flax is presented, including the gene structures, phylogeny and conserved motifs, resulting in the presence of central hydrophobic regions with proline knot motif, providing an evolutionary proof that this central hydrophobic region had evolved from duplications in the primitive eukaryotes. In addition, expression analysis of L-oleosin and caleosin genes using quantitative real-time PCR demonstrated that seed contained their maximum expression, except that RcCLO-1 expressed maximum in cotyledon. Thus, our comparative genomics analysis of oleosin and caleosin genes and their putatively encoded proteins in two non-model plant species provides insights into the prospective usage of gene resources for improving OB-stability.     

Global Analysis of Ankyrin Repeat Domain C3HC4-Type RING Finger Gene Family in Plants

Ankyrin repeat (ANK) C3HC4-type RING finger (RF) genes comprise a large family in plants and play important roles in various physiological processes of plant life. In this study, we identified 187 ANK C3HC4-type RF proteins from 29 species with complete genomes and named the ANK C3HC4-type RF proteins the XB3-like proteins because they are structurally related to the rice (Oryza sativa) XB3. A phylogenetic relationship analysis suggested that the XB3-like genes originated from ferns, and the encoded proteins fell into 3 major groups. Among these groups, we found that the spacing between the metal ligand position 6 and 7, and the conserved residues, which was in addition to the metal ligand amino acids, in the C3HC4-type RF were different. Using a wide range of protein structural analyses, protein models were established, and all XB3-like proteins were found to contain two to seven ANKs and a C3HC4-type RF. The microarray data for the XB3-like genes of Arabidopsis, Oryza sative, Zea mays and Glycine max revealed that the expression of XB3-like genes was in different tissues and during different life stages. The preferential expression of XB3-like genes in specified tissues and the response to phytohormone and abiotic stress treatments of Arabidopsis and Zea mays not only confirmed the microarray analysis data but also demonstrated that the XB3-like proteins play roles in plant growth and development as well as in stress responses. Our data provide a very useful reference for the identification and functional analysis of members of this gene family and also provide a new method for the genome-wide analysis of gene families.     

Comparative genomic analysis of the <Emphasis Type="Italic">Lipase3</Emphasis> gene family in five plant species reveals distinct evolutionary origins

Lipases are physiologically important and ubiquitous enzymes that share a conserved domain and are classified into eight different families based on their amino acid sequences and fundamental biological properties. The Lipase3 family of lipases was reported to possess a canonical fold typical of α/β hydrolases and a typical catalytic triad, suggesting a distinct evolutionary origin for this family. Genes in the Lipase3 family do not have the same functions, but maintain the conserved Lipase3 domain. There have been extensive studies of Lipase3 structures and functions, but little is known about their evolutionary histories. In this study, all lipases within five plant species were identified, and their phylogenetic relationships and genetic properties were analyzed and used to group them into distinct evolutionary families. Each identified lipase family contained at least one dicot and monocot Lipase3 protein, indicating that the gene family was established before the split of dicots and monocots. Similar intron/exon numbers and predicted protein sequence lengths were found within individual groups. Twenty-four tandem Lipase3 gene duplications were identified, implying that the distinctive function of Lipase3 genes appears to be a consequence of translocation and neofunctionalization after gene duplication. The functional genes EDS1, PAD4, and SAG101 that are reportedly involved in pathogen response were all located in the same group. The nucleotide diversity (Dxy) and the ratio of nonsynonymous to synonymous nucleotide substitutions rates (Ka/Ks) of the three genes were significantly greater than the average across the genomes. We further observed evidence for selection maintaining diversity on three genes in the Toll-Interleukin-1 receptor type of nucleotide binding/leucine-rich repeat immune receptor (TIR-NBS LRR) immunity-response signaling pathway, indicating that they could be vulnerable to pathogen effectors.     

A comprehensive expression analysis of all members of a gene family encoding cell-wall enzymes allowed us to predict cis-regulatory regions involved in cell-wall construction in specific organs of Arabidopsis.

The Arabidopsis thaliana genome sequencing project has revealed that multigene families, such as those generated by genome duplications, are more abundant among plant genomes than among animal genomes. To gain insight into the evolutionary implications of the multigene families in higher plants, we examined the XTH gene family, a group of genes encoding xyloglucan endotransglucosylase/hydrolase, which are responsible for cell-wall construction in plants. Expression analysis of all members (33 genes) of this family, using quantitative real-time RT-PCR, revealed that most members exhibit distinct expression profiles in terms of tissue specificity and responses to hormonal signals, with some members exhibiting similar expression patterns. By comparing the flanking sequences of individual genes, we identified four sets of large-segment duplications and two sets of solitary gene duplications. In each set of gene duplicates, long nucleotide sequences, ranging from one to two hundred base pairs, are conserved. Furthermore, gene duplicates exhibit similar organ-specific expression profiles. These facts allowed us to predict putative cis-regulatory regions, particularly those responsible for cell-wall construction, and hence for morphogenesis, that are specific for certain organs or tissues in plants.     

Genetic diversity of NBS–LRR class disease-resistance gene analogs in cultivated and wild eggplants

The genes encoding the nucleotide-binding site (NBS) and leucine-rich repeat (LRR) motifs constitute a large gene family in plants and have attracted much interest, because most of the plant disease-resistance genes that have been cloned are from this gene family. In this study, degenerate oligonucleotide primers, designed on the basis of conserved regions of the NBS domains from known plant resistance genes, were used to isolate resistance gene analogs (RGAs) from cultivated and wild eggplants, i.e., S. melongena, S. aethiopicum gr. Gilo, S. linnaeanum, S. integrifolium, S. sisymbriifolium, and S. khasianum. Sequence analysis indicated that the cloned eggplant RGAs belong to the non-TIR–NBS–LRR type, which are very similar to the R genes or the RGAs identified in other plant species, especially Solanaceae plants, suggesting the existence of common ancestors. Wide genetic diversity of eggplant RGAs was observed both in interspecific and intraspecific sequences, and eight distinct families of eggplant RGAs were identified. Further studies revealed a high average ratio of synonymous to non-synonymous substitution and a low level of recombination. These results suggest that NBS-encoding sequences of RGAs in cultivated and wild eggplants are subject to gradual accumulation of mutations leading to purifying selection. This is the first report of NBS–LRR class RGAs in eggplants.     

Genome-wide Comparative Analysis of Annexin Superfamily in Plants

Most annexins are calcium-dependent, phospholipid-binding proteins with suggested functions in response to environmental stresses and signaling during plant growth and development. They have previously been identified and characterized in Arabidopsis and rice, and constitute a multigene family in plants. In this study, we performed a comparative analysis of annexin gene families in the sequenced genomes of Viridiplantae ranging from unicellular green algae to multicellular plants, and identified 149 genes. Phylogenetic studies of these deduced annexins classified them into nine different arbitrary groups. The occurrence and distribution of bona fide type II calcium binding sites within the four annexin domains were found to be different in each of these groups. Analysis of chromosomal distribution of annexin genes in rice, Arabidopsis and poplar revealed their localization on various chromosomes with some members also found on duplicated chromosomal segments leading to gene family expansion. Analysis of gene structure suggests sequential or differential loss of introns during the evolution of land plant annexin genes. Intron positions and phases are well conserved in annexin genes from representative genomes ranging from Physcomitrella to higher plants. The occurrence of alternative motifs such as K/R/HGD was found to be overlapping or at the mutated regions of the type II calcium binding sites indicating potential functional divergence in certain plant annexins. This study provides a basis for further functional analysis and characterization of annexin multigene families in the plant lineage.     

Isolation of a monocot 3-hydroxy-3-methylglutaryl coenzyme A reductase gene that is elicitor-inducible

The rice (Oryza sativa) phytoalexins, momilactones and oryzalexins, are synthesized by the isoprenoid pathway. An early step in this pathway, one that is rate-limiting in mammalian systems, is catalyzed by the enzyme 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR). A gene that encodes this enzyme has been isolated from rice, and found to contain an open reading frame of 1527 bases. The encoded protein sequence of the rice HMGR appears to be conserved with respect to other HMGR proteins, and 1 or 2 membrane-spanning domains characteristic of plant HMGRs are predicted by a hydropathy plot of the amino acid sequence. The protein is truncated at its 5 end, and shows reduced sequence conservation in this region as compared to other plant sequences. The rice genome contains a small family of HMGR genes. The isolated gene, HMGR I, is expressed at low levels in both vegetative and floral organs of rice plants. It is not induced in plants by wounding, but is strongly and rapidly induced in suspension cells by a fungal cell wall elicitor from the pathogenMagnaporthe grisea, causal agent of rice blast disease. This suggests that HMGR I may be important in the induction of rice phytoalexin biosynthesis in response to pathogen attack, and therefore may play a key role as a component of the inducible defense mechanism in rice.     

The F-box gene family is expanded in herbaceous annual plants relative to woody perennial plants

F-box proteins are generally responsible for substrate recognition in the Skp1-Cullin-F-box complexes that are involved in protein degradation via the ubiquitin-26S proteasome pathway. In plants, F-box genes influence a variety of biological processes, such as leaf senescence, branching, self-incompatibility, and responses to biotic and abiotic stresses. The number of F-box genes in Populus (Populus trichocarpa; approximately 320) is less than half that found in Arabidopsis (Arabidopsis thaliana; approximately 660) or Oryza (Oryza sativa; approximately 680), even though the total number of genes in Populus is equivalent to that in Oryza and 1.5 times that in Arabidopsis. We performed comparative genomics analysis between the woody perennial plant Populus and the herbaceous annual plants Arabidopsis and Oryza in order to explicate the functional implications of this large gene family. Our analyses reveal interspecific differences in genomic distribution, orthologous relationship, intron evolution, protein domain structure, and gene expression. The set of F-box genes shared by these species appear to be involved in core biological processes essential for plant growth and development; lineage-specific differences primarily occurred because of an expansion of the F-box genes via tandem duplications in Arabidopsis and Oryza. The number of F-box genes in the newly sequenced woody species Vitis (Vitis vinifera; 156) and Carica (Carica papaya; 139) is similar to that in Populus, supporting the hypothesis that the F-box gene family is expanded in herbaceous annual plants relative to woody perennial plants. This study provides insights into the relationship between the structure and composition of the F-box gene family in herbaceous and woody species and their associated developmental and physiological features.     

Complete mitochondrial genome sequences of the three pelagic chaetognaths Sagitta nagae,Sagitta decipiens and Sagitta enflata

The complete nucleotide sequences of the mitochondrial genomes were determined for the three pelagic chaetognaths, Sagitta nagae, Sagitta decipiens, and Sagitta enflata. The mitochondrial genomes of these species which were 11,459, 11,121, and 12,631 bp in length, respectively, contained 14 genes (11 protein-coding genes, one transfer RNA gene, and two ribosomal RNA genes), and were found to have lost 23 genes that are present in the typical metazoan mitochondrial genome. The same mitochondrial genome contents have been reported from the benthic chaetognaths belonging to the family Spadellidae, Paraspadella gotoi and Spadella cephaloptera. Within the phylum Chaetognatha, Sagitta and Spadellidae are distantly related, suggesting that the gene loss occurred in the ancestral species of the phylum. The gene orders of the three Sagitta species are markedly different from those of the other non-Chaetognatha metazoans. In contrast to the region with frequent gene rearrangements, no gene rearrangements were observed in the gene cluster encoding COII–III, ND1–3, srRNA, and tRNA^met. Within this conserved gene cluster, gene rearrangements were not observed in the three Sagitta species or between the Sagitta and Spadellidae species. The gene order of this cluster was also assumed to be the ancestral state of the phylum.