- Research article
- Open Access
Molecular diversity of phospholipase D in angiosperms
BMC Genomics volume 3, Article number: 2 (2002)
The phospholipase D (PLD) family has been identified in plants by recent molecular studies, fostered by the emerging importance of plant PLDs in stress physiology and signal transduction. However, the presence of multiple isoforms limits the power of conventional biochemical and pharmacological approaches, and calls for a wider application of genetic methodology.
Taking advantage of sequence data available in public databases, we attempted to provide a prerequisite for such an approach. We made a complete inventory of the Arabidopsis thaliana PLD family, which was found to comprise 12 distinct genes. The current nomenclature of Arabidopsis PLDs was refined and expanded to include five newly described genes. To assess the degree of plant PLD diversity beyond Arabidopsis we explored data from rice (including the genome draft by Monsanto) as well as cDNA and EST sequences from several other plants. Our analysis revealed two major PLD subfamilies in plants. The first, designated C2-PLD, is characterised by presence of the C2 domain and comprises previously known plant PLDs as well as new isoforms with possibly unusual features-catalytically inactive or independent on Ca2+. The second subfamily (denoted PXPH-PLD) is novel in plants but is related to animal and fungal enzymes possessing the PX and PH domains.
The evolutionary dynamics, and inter-specific diversity, of plant PLDs inferred from our phylogenetic analysis, call for more plant species to be employed in PLD research. This will enable us to obtain generally valid conclusions.
Phospholipase D (PLD, EC 18.104.22.168.) is a ubiquitous eukaryotic enzyme participating in various cellular processes (for a review see [1, 2]). Biochemically distinct types of PLDs have been described, but only two, the mammalian glycosylphosphatidylinositol-specific PLD (GPI-PLD) and a family usually referred to as phosphatidylcholine-specific PLD (PC-PLD), have been characterised also on the molecular level. Two distinct PC-PLD genes have been identified in mammals; they seem to be involved in signal transduction and vesicular trafficking. The yeast Saccharomyces cerevisiae contains only one gene from the PC-PLD family and its function in sporulation has been recognised.
Plants are a traditional model for PLD research. Indeed, PLD activity was first described from a plant source , and the first cloned eukaryotic cDNA coding for a PLD was isolated from the castor bean, Ricinus communis. Using mainly biochemical and pharmacological approaches, plant PLD has been implicated in many cellular processes (reviewed in ). Beside its roles in membrane degradation and turnover during senescence, seed germination and under stress conditions, plant PLD is emerging as an important component of signal transduction cascades, e.g. in response to wounding, abscisic acid  or Nod factors . Earlier pharmacological evidence for the involvement of heterotrimeric G-proteins in plant PLD regulation  has been recently strengthened by a report on direct interaction of an alpha subunit of a G-protein with PLDα in tobacco . Also products of PLD action, i.e. phosphatidic acid (PA), diacylglycerol and N-acylethanolamine, are potential signalling molecules in plants (reviewed in [8, 9]).
Up to now about 20 PLDs have been cloned from plants. Multiple isoforms have been found in some species, complicating the study of plant PLD. Application of a reverse genetic approach, combining the knowledge of genomic sequences and molecular genetic techniques, holds greatest promises here. This can be documented e.g. by successful inactivation of the AtPLDα1 gene in Arabidopsis thaliana by antisense strategy, which allowed identification of a novel PLD activity in plants [10, 11].
Thorough characterisation of the gene family concerned is an obvious prerequisite for productive application of the reverse-genetic approach. Here we present the results of a detailed comparative analysis of the Arabidopsis and rice PLD families, combining data from the complete Arabidopsis genome sequence , publicly available rice (Oryza sativa) genomic and cDNA sequences, and the draft rice genome data made available by Monsanto . Extensive EST collections from three plant species, tomato (Lycopersicon esculentum), Medicago truncatula and Sorghum bicolor, have been included into the analysis to provide the insight into the inter-specific variability of plant PLDs. Our results indicate that the angiosperm PLD family, comprising two major subfamilies (C2- and PXPH-PLDs), is evolutionarily very dynamic, and conclusions based on a single species (such as Arabidopsis) might not therefore be simply applicable to others.
Results and Discussion
A dozen Arabidopsis PLDs
Up to now several cDNAs representing six distinct PLD-encoding genes have been reported from Arabidopsis (Table 1). Using the cloned PLDs from Arabidopsis and other organisms we conducted exhaustive BLAST searches of the Arabidopsis sequences available from GenBank and found 12 genes from the eukaryotic PC-PLD family, five of them not yet recorded in the literature. All genes found code for proteins containing all the conserved sequence motifs characteristic of eukaryotic PLDs, including two copies of the invariant catalytic HxKxxxxD motif , suggesting that all probably posses the genuine PLD enzymatic activity (thougt this must be proven experimentally). The catalytic HxKxxxxD motif is shared also by other proteins put together with the eukaryotic PLDs into the PLD superfamily , we however did not identify any other members of the superfamily in Arabidopsis besides the 12 PC-PLDs.
Before attempting a detailed phylogenetic analysis of the Arabidopsis PLDs, we used a combination of computational tools, comparison with cDNAs/ESTs and information from protein alignments to verify the exon-intron structures proposed by AGI annotators (see Materials and Methods). In several cases, prediction ambiguities and cloning or sequencing errors have been uncovered, and refined gene models and protein sequence predictions have been obtained and used in further analysis (see the discussion below and Additional file 1.
All Arabidopsis PLDs can be classified into two subfamilies (Fig. 1A). Since these subfamilies differ by the presence of distinct N-terminal phospholipid-binding domains (C2 vs. PX-PH), we propose denoting them as C2-PLDs and PXPH-PLDs. The C2-PLD subfamily contains 10 Arabidopsis isoforms harbouring a phospholipid/Ca2+-binding fold called the C2 (or CalB) domain (Fig. 2, 5). Seven genes have been already found in the genome, including AtPLDα1 (formerly PLDα), AtPLDβ1 (previously PLDβ) and a tandem triplication of AtPLDγ1, AtPLDγ2 and AtPLDγ3. A gene tentatively designated as PLDδ1  should in our view be more suitably labelled AtPLDβ2 (see below), while we keep the term AtPLDδ for a gene recently reported with this designation [15, 16]. (see Table 1).
We suggest the term AtPLDα2 for a newly identified homolog closely related to AtPLDα1 (88 % sequence identity at the protein level). Interestingly, AtPLDα1 and AtPLDα2 genes reside within one of several large-scale intragenomic duplications believed to be remnants of a tetraploidisation event dated 112 Myr ago , pointing toward probable evolutionary origin of these two paralogs. The remaining Arabidopsis C2-PLD genes do not correspond to any of the previously established group, so we propose terming them AtPLDε and AtPLDζ. (although a PLDε has already been mentioned in a recent review , it is not clear from the text to which of the Arabidopsis PLD genes it corresponds to.)
Members of the PXPH-PLD subfamily typically bear different two phospholipid-binding domains, the PX (phox) and the PH (pleckstrin-homology) domain, in the N-terminal region (Fig. 2). PXPH-PLDs have previously been known only from animals and fungi, so identification of two Arabidopsis genes belonging to this subfamily adds an important new dimension to the picture of plant PLD. We propose terming these genes AtPLDp1 and AtPLDp2 (the "p" from PX and PH) to underline the principal difference between the C2-PLD branch (Greek letters) and the PXPH-PLD subfamily. Very recently, a cDNA has appeared in GenBank corresponding to AtPLDp1 gene but unfortunately annotated as PLD zeta1 (Table 1). As we believe that our nomenclature better reflects structural and phylogenetic aspects of the plant PLD family, in the rest of this text, the gene names proposed by us are used.
Functionality of those Arabidopsis PLD genes, for which full-length cDNAs have been cloned, is undisputed. Moreover, proteins encoded by three of these genes, PLDβ1, PLDγ1 and PLDδ, have been characterised biochemically (see [2, 16]). Expression of several other isoforms is documented by ESTs in GenBank, but there are currently no ESTs available cognate for AtPLDα2, AtPLDβ2 and AtPLDζ genes (Table 1). However, absence of cognate ESTs is not exceptional, since in general only about 60% of predicted Arabidopsis genes are recorded in available EST collections . It is therefore very likely that expression of most of the genes without ESTs is very low or limited only to some special developmental stages or conditions.
Exon-intron organisation of Arabidopsis PLD genes
Limitations of theoretical prediction of exon-intron structures are well known and cDNA sequencing is often necessary for building accurate gene models. This proves true also for many of the Arabidopsis PLDs. Unfortunately, four reported cDNAs, i.e. AtPLDα1, AtPLDβ1, AtPLDγ1 and AtPLDγ2 (Table 1), contain mismatches compared to the highly accurate genomic sequences (reported to contain less than 1 error per 104-105 bp; ). While some of the discrepancies may represent a natural polymorphism, others, particularly those associated with frame shifts, are most likely due to sequencing errors or cloning artefacts. This suspicion is also supported by available EST sequences, which nearly always match the genomic sequences, not the cDNAs.
For example, within the coding portion of the AtPLDα1 cDNA there are four regions with the reading frame shifted relative to the genomic sequence. As a result, the protein sequence derived from the cDNA (AAC49274.1) is highly divergent from other PLDs in these four regions, while that predicted from the genome data (NP_188194.1) matches well the PLD consensus. We found similar discrepancies also for AtPLDβ1, AtPLDγ1 and AtPLDγ2 cDNAs. Moreover, published AtPLDγ1 and AtPLDγ2 cDNAs appear to be chimeric, perhaps due to cloning artefacts. The last ~180 nucleotides of AtPLDγ1 cDNA apparently originate from a gene encoding a pseudo-response regulator (AB046955, chromosome 5). Similarly, the 3' third of the cDNA reported as AtPLDγ2 is actually derived from the AtPLDγ3 gene. We therefore believe that the cDNA sequences have to be interpreted very cautiously, and we base our conclusions mainly on the genome project data. In several cases, however, we proposed corrections of the AGI annotation of PLD genes. Details and refined coding sequences can be found in the Additional files, most important aspects are also discussed below.
Despite sequencing errors, the AtPLDα1 cDNA is in good agreement with the previously suggested gene structure, the gene contains three coding exons and a 5' non-coding one (Fig. 3), similarly to other characterised PLD genes from the alpha subgroup . Coding portion of the AtPLDα2 gene appears to be arranged in the same way, but whether there is also a non-coding exon in the 5' UTR that could be proven only by cloning of the respective cDNA. Two other predicted Arabidopsis genes, AtPLDε and AtPLDζ appear also to exhibit structures related to alpha-type PLDs, but differ by the presence of an additional intron. This intron resides at different positions in both genes, suggesting that it has been acquired independently by AtPLDε and AtPLDζ.
AtPLDβ1 gene was found to consist of 10 exons . Current database annotation should, however, be corrected in some points (see Table 1 and Additional file for details). Most importantly, there is a long region devoid of in-frame STOP-codons upstream from the first predicted exon, and the ORF could be thus extended in the 5' direction (Fig. 3). This would add an unusual, ~250 aa long N-terminal projection to the PLDβ1 protein. However, the potential initiation codon upstream of the predicted one is just out of the region covered by the respective cDNA, suggesting that the cDNA might be 5' truncated. There are several reasons arguing for inclusion of the N-terminal extension. First, such a long sequence would have most likely accumulated STOP-codons if not translated. Second, we found such an N-terminal extension (though loosely conserved) in other PLDs of the beta group, including also Arabidopsis PLDβ2, tomato PLDβ2 (; AY013256) and two PLDβ isoforms from rice (see Table 2). Third, a cotton cDNA coding for a PLDβ (AF159139) contains a very long 5' leader sequence (more than 1300 bp) harbouring a potential ORF coding for a peptide ~300 aa long and similar to the N-terminal extensions of both beta PLDs from Arabidopsis. The 3' end of the leader ORF overlaps with the beginning of the PLD ORF but resides in a different reading frame. It is tempting to speculate that the frame shift is due to a cloning or sequencing error and that the cotton PLDβ possesses an N-terminal extension similar to other beta PLDs. The extensions in beta PLDs generally do not appear to resemble a sorting signal for any cellular compartment, so we suggest that they might e.g. mediate interaction with regulatory factors.
The AtPLDβ2 gene was originally described as PLDδ1, and 11 exons predicted by AGI were proposed as a unique feature . However, this hypothetical gene structure was not supported by the results of gene-finding programmes that we employed (see Materials and Methods). Inclusion of the 1st originally predicted intron into the ORF, which is supported by the programmes, introduces a conserved portion of the C2 domain and adjusts the splicing pattern to the 10-exonic scheme exhibited by several other PLD genes. The resulting predicted protein sequence belongs clearly to a beta PLD type (Fig. 1A,B).
The three very similar paralogs of PLDγ reside in a tandem triplication (arranged AtPLDγ1 – AtPLDγ3 – AtPLDγ2) on the Arabidopsis chromosome 4 , indicating a relatively recent origin of the triplet. The predicted gene structure of all three genes fits the 10-exonic scheme typical for some other PLD types. However, there appears to be a probable obscure intron of 96 nucleotides in the AtPLDγ2 gene delimited by GT-AG borders and supported by a matching cDNA sequence (see Fig. 3). Exclusion of this intron (possibly as a result of alternative splicing) deletes a region from within the C2 domain that might substantially affect function and/or regulation of AtPLDγ2. A similar potential intron may be present also in AtPLDγ3 gene. The AtPLDγ2 gene further contains an additional non-coding exon at the 5' end, but it is currently unknown whether the 5' UTR of the other two genes is organised in the same way (the available AtPLDγ1 cDNA may be 5' truncated).
The AtPLDδ gene had been predicted by AGI annotators as consisting of 16 exons, but, as revealed by EST sequences and cDNAs, AtPLDδ possesses only 10 conserved exons shared with beta and gamma PLDs. Interestingly, there is an evidence for alternative splicing of the AtPLDδ gene, because one of the independently cloned cDNAs (AB031047, ) differs from the others by extension of the second exon at the 3' boundary by 33 nucleotides (Table 1, Fig. 3). Interestingly, the apparently more abundant shorter alternative (as defined also by three ESTs covering the respective region) uses an unconventional 5'GC intron boundary instead of GT (see http://www.Arabidopsis.org/splice_site_excep.html). The longer and the shorter putative splice variants have recently been denoted PLDδa and PLDδb, respectively . The PLDδa protein bears an insertion of 11 aa in an otherwise relatively conserved region, which could have profound functional consequences.
The two genes classified into the PXPH-PLD subfamily appear to exhibit the most complex exon-intron structure of all Arabidopsis PLDs. A corresponding full-length cDNA has been reported only for AtPLDp1 (Table 1), so the prediction of the AtPLDp2 gene remains tentative. A minor correction should perhaps be introduced into the current database prediction of the AtPLDp2 gene to restore a highly conserved region (see Table 1 and the Additional files). Despite a difference in the number of exons (20 and 16, respectively), the structures of AtPLDp1 and AtPLDp2 genes are clearly related, as the difference is due to 4 introns probably lost from AtPLDp2 (20-exonic structure seems to be primordial in the plant PXPH-PLD subfamily, see the rice homologs below and Fig. 1).
PLDs in rice: an alternative view
Arabidopsis is presently the only plant for which the complete PLD set can be catalogued. Nonetheless, other species are emerging as important models for genome-wide studies. Rice genome sequencing is highly advanced, with a substantial portion (more than 230 Mbp up to now, see http://rgp.dna.affrc.go.jp for updates) already sequenced by the International Rice Genome Sequencing Project (IRGSP). Even greater portion of the genome (about 250 Mbp) has been sequenced by the Monsanto company, who have made their sequences publicly available . With redundancy between the two data resources taken into account, we could analyse about three quarters of the whole rice genome.
We identified at least 16 complete or partial sequences of putative rice PLD genes (Table 2). Five of them have been cloned individually, 13 genes or their portions have been already sequenced by the IRGSP, sequences coming from 13 genes could be found in the Monsanto genome draft and fragments of at least one PLD gene are available only as EST or GSS sequences. Since a systematic nomenclature of rice PLD genes has not been established, we propose a terminology that would reflect phylogenetic and structural relationship with PLDs from other species (Table 2, Fig. 1A,B). Several rice PLDs appear to be orthologous to Arabidopsis genes, including two alpha-type PLDs and two beta-like PLDs. At least one presumed ortholog of AtPLDδ is available only as GSS and EST sequences too fragmentary to be included directly into phylogenetic analysis. Their assignment to the delta type is however supported by an analysis employing a higly similar barley PLD assembled from EST sequences and used as a placeholder (Fig. 1B). Two genes from the PXPH-PLD subfamily were also found and denoted OsPLDp1 and OsPLDp2. Remaining genes cannot be assigned to any of the classes established for Arabidopsis. Although OsPLDλ, tends to cluster together with AtPLDε, the two proteins share only about 42% identical amino acids and the genes appear to differ in the number of introns, supporting classification of OsPLDλ as a novel subtype. Similarly, according to the phylogenetic analysis OsPLDθ and OsPLDμ appear to be related, but different intron numbers and the second catalytic motif missing from OsPLDθ (see below) justify a separate classification.
Only the five individually cloned PLD genes have been annotated. Complete cDNA has been reported for OsPLDα1; the gene has the exon-intron structure closely related to other alpha-type PLDs [2, 21]. Delimitation of the coding region of the OsPLDη1 gene has also been verified experimentally and found to have a similar organisation . Structures of the other genes could be predicted only theoretically, but comparison with EST sequences and other PLD genes proved helpful, as exon-intron junctions appear to be highly conserved within individual subgroups of the PLD family (see below; predicted or corrected coding sequences available in Additional files). Thus, we introduced a minor correction into the previously proposed OsPLDν1 gene structure (see Table 2 and Additional files). Annotated OsPLDη2 and OsPLDη3 appear to have a similar splicing pattern as OsPLDη1. Interestingly, we found the three OsPLDη genes residing in the genome adjacently in a series OsPLDη2-OsPLDη3-OsPLDη1, but, in contrast to the AtPLDγ cluster in Arabidopsis, OsPLDη2 is inverted with respect to the remaining two genes. Exon-intron structures proposed by us for other rice PLD genes reflect phylogenetic affinity to Arabidopsis orthologs (compare Table 1 and Table 2). The novel OsPLDλ, and OsPLDθ genes probably have 3 coding exons with introns occupying conserved positions shared with the PLDα and PLDη prototypes. The novel OsPLDμ gene also resembles PLDα and PLDη, although comparison with a highly similar barley EST revealed 4 coding exons. The second exon is very short and encodes a part of the first non-conserved loop of the C2 domain (Fig. 5). This exon might be completely novel or, more probably, has been insulated from the downstream coding region by acquisition of an intron. Gene structure of OsPLDκ, comprising six coding exons, appears to be derived from the 10-exonic scheme, retained by related PLDs (i.e. β, γ, δ and ν), by loss of four introns. Two identified rice PXPH-PLDs are very similar to each other including exon-intron organisation, which is obviously shared with AtPLDp1 from Arabidopsis.
In summary, comparison of Arabidopsis and rice PLD genes revealed that they exhibit generally non-conserved exon-intron structures (Table 1, Table 2; Fig. 3), and positions of introns do not reflect boundaries between functional domains. However, three clusters of plant PLD genes can be recognised, which differ completely in exon-intron organisation from each other, but the organisation appear conserved among genes within each cluster. Independent acquisition of introns seems to be a plausible explanation (Fig. 1A,B; Fig. 3). Described plant PXPH-PLD genes appear to have primarily 20 conserved exons, with the exception of AtPLDp2 lacking four introns. The C2-PLD subfamily comprises a clade characterised primarily with 3 coding exons, and a group of originally 10-exonic genes. Independent acquisition of additional introns (in AtPLDε, AtPLDζ and probably also OsPLDμ) or intron losses (suggested for OsPLDκ and AtPLDp2) may be common during the evolution of plant PLDs.
PLD diversity and expression as recorded in the EST collections
Databases of expressed sequence tags (ESTs) are available for a number of plant species and represent invaluable resource for both functional and evolutionary studies, providing information on both genetic diversity and expression profiles. To assess these aspects of the angiosperm PLD family, we identified a number of PLD-derived ESTs from Arabidopsis, rice, tomato, Medicago truncatula and Sorghum bicolor (see Table 1 and Table 2, and Additional files).
For exploration of PLD diversity beyond Arabidopsis and rice, tomato is a suitable starting point, with more than 140,000 ESTs available and five full-length PLD cDNAs cloned representing three PLDα and two PLDβ genes [19, 22, 23]. With the exception of LePLDα2, all cloned tomato PLDs are recorded among ESTs, but expression of additional isoforms is documented, too, including a PLD similar to alpha types, at least two putative delta isoforms, a PLD most similar to AtPLDε and at least one gene from the PXPH-PLD group. The second species analysed was Medicago truncatula with more than 137,000 ESTs in GenBank. No PLD has yet been reported from this plant, but ESTs again indicate the presence of a complex PLD family, comprising at least two indisputable PLDα homologs highly similar to each other, at least two additional genes less similar to alpha types, potentially three PLDβ isoforms, at least one delta ortholog, a PLD most similar to AtPLDε and two members of the PXPH-PLD subfamily. As a monocotyledonous model for EST analysis we chose Sorghum bicolor, for which more that 84,000 ESTs had been sequenced. Multiple homologs could again be found among the ESTs, including at least two obvious alpha PLDs, a gene related to the rice PLDη1, one PLDβ and a PLD most similar to the rice PLDμ.
In summary, our EST analysis revealed that PLD types identified in Arabidopsis and rice are widespread in angiosperms, but there might be additional types not yet characterised. With the help of EST clones, full-length genes/cDNAs can be easily isolated and characterised, so deeper insight into PLD diversity in plants can become soon available (the list of ESTs analysed is available in Additional files).
The relative abundance of ESTs can provide information on expression of individual genes . Unfortunately, for most PLD genes there are too few ESTs for statistically significant estimation of their expression in specific tissues, developmental stages or conditions, and only general level of expression can be inferred. According to the total number of cognate ESTs in the GenBank, the most highly expressed PLD gene in Arabidopsis is AtPLDα1 (40 EST entries) followed with AtPLDδ (29 entries), while other genes seem to be expressed at a considerably lower level or not recorded at all (see Table 1). In rice, expression of OsPLDα1 predominates to a similar extent as in Arabidopsis (~40% of all ESTs from PLD genes), and expression of the other genes is markedly lower as well (Table 2). The EST collection from Medicago provides even more strongly substantiated evidence for an expression bias with 63% (of 74 PLD-derived ESTs in total) matching one of multiple PLDα paralogs.
Among tomato ESTs that can be assigned to the cloned cDNAs, 14 come from LePLDα1, 2 from LePLDα3, 4 from LePLDβ1 and 1 from LePLDβ2, suggesting that expression of LePLDα1 might again be prevalent. However, LePLDα2 and LePLDα3 have higher expression levels than LePLDα1 when measured by Northern blots , so caution must be paid when few ESTs are used for conclusions on expression profiles. In Sorghum one of two alpha-type genes also accounts for ~40% of PLD-derived ESTs, but, in contrast to previous collections, a PLDβ-like isoform appear to be sampled to a similar extent. Interestingly, all the EST corresponding to the latter gene are derived from a cDNA library prepared from a pathogen-infected plants. It is tempting to speculate that the expression of this PLDβ gene might be induced by a pathogen-derived signal, similarly to LePLDβ1 reported to be induced upon treatment with an elicitor xylanase .
Predominant expression of alpha-type PLDs inferred from our EST analysis fits with biochemical experience, since the enzymatic activity usually ascribed to alpha-type PLDs is much more abundant in plant tissues compared to the activity of the beta and gamma types . Interestingly, in most plants studied (except for tomato), two PLDα genes could be found, but only one of them was highly expressed (see Table 1, Table 2 and Additional files). Similarly, two PLDα paralogs have been cloned from the resurrection plant Craterostigma plantagineum, one of them expressed constitutively and the other one induced only upon desiccation stress . Differential expression mode for two very similar PLDs has been observed also in tomato, where the elicitor xylanase stimulated expression of LePLDβ1 but not of LePLDβ2. Henceforth, if the differences in expression did relate to differences in physiological function, it could be concluded that there is only little functional redundancy within plant PLD family, even among highly similar isoforms.
Functional aspects of the primary structure of plant PLDs
As already noted, all known eukaryotic PC-PLDs belong to two subfamilies differing in their N-terminal portion (Fig. 2), but within the core of the enzyme several highly conserved regions shared by all PLDs have been recognised [1, 6]. Among these, the most important are two copies of HxKxxxxD (or HKD) motif (Fig. 3, 4A) with the three residues absolutely conserved. The two motifs form together a catalytic site and their mutation abolishes enzyme activity. Our inspection of protein sequences of both characterised and predicted plant PLDs revealed that both HKD motifs are present in all isoforms, with the exception of rice PLDθ, which, unexpectedly, has the second copy mutated (Fig. 4A). Also the tomato PLDα3 has been reported to lack the aspartate residue from the second HKD motif, and its functionality has been questioned . However, we found that LePLDα3 actually does posses the second catalytic motif conserved, as the conclusion of Laxalt et al.  is based on a protein sequence deduced from an un-spliced transcript. Indeed, a spliced cDNA corresponding to the same gene has been cloned independently and codes for a functional enzyme . On the contrary, the replacement of histidine and lysine residues in the second motif of OsPLDθ does not seem to be an artefact (the region has been sequenced independently by Monsanto and IRGSP). It is tempting to speculate that the catalytically inactive OsPLDθ protein, if expressed, may fulfil regulatory roles via competition with functional PLDs and/or other proteins in binding of phospholipids.
C2- and PXPH-PLDs are believed to differ in their dependence on Ca2+. Animal and fungal PLDs are not directly dependent on Ca2+, and the same is likely also for plant PXPH-PLDs. On the other hand, most (but perhaps not all, see below) C2-PLDs will exhibit dependence on and regulation by Ca2+, as the C2 domains usually bind phospholipids in a Ca2+-dependent manner . Structural characterisation of several C2 domains revealed that three Ca2+-coordinating sites occupied by aspartate or asparagine residues are used generally, while other ligands are specific for individual domains (see  and Fig. 5). Different concentrations of Ca2+ required for optimal stimulation of alpha, beta and gamma PLD types have been explained in terms of the number of expected Ca2+-coordinating residues present within the C2 domains of individual PLDs. It was noted that beta and gamma PLDs in Arabidopsis posses all the binding residues conserved in characterised C2 domains, whereas PLDα lacks two of them because of substitution [27, 28]. However, this hypothesis relies on a protein derived from the non-accurate sequence of AtPLDα1 cDNA, with an asparagine residue accidentally substituted for by an isoleucine at the position of the third conserved Ca2+-coordinating site (Fig. 1 in ). The asparagine residue, occurring actually in all sequenced alpha-type PLDs (Fig. 5), is a Ca2+-binding ligand in the C2 domain of phospholipase A2. Nevertheless, only some PLDα isoforms posses an aspartate residue at the site corresponding to the second generally shared Ca2+-binding ligand, and the first position is mutated in all alpha-type PLDs (Fig. 5), in agreement with the higher [Ca2+] requirement of PLDα compared to PLDβ and γ , and with lower affinity to Ca2+ of the C2 from AtPLDα1 compared to the domain from AtPLDβ1 . Indeed, all PLDs belonging to the 10-exonic cluster (i.e. PLDβ, γ, δ, κ and v) do posses all the three ligands conserved and most closely resemble the C2 domain from PLCδ1. On the other hand, the remaining PLD types, i.e. PLDε, ζ, η, θ and λ, lack all the three Ca2+-binding residues due to substitution or deletion, while OsPLDμ retains only one (Fig. 5). Some C2 domains have been described that do not bind Ca2+, but nonetheless do target respective proteins to the membranes . It is therefore tempting to speculate that the C2 domains of the latter PLD types might function similarly and that these PLDs are perhaps Ca2+-independent. However, still different modes of Ca2+ binding employing different residues may be discovered when additional C2 domains are structurally characterised, and threfore no definite functional predictions can be made from sequence only.
All characterised PC-PLDs from both major subfamilies are stimulated by or even dependent on PIP2 under physiological or near-physiological conditions [1, 2, 29]. Mammalian and yeast PLDs appear to interact with PIP2 by the PH domain  and via a novel highly conserved motif located between the two copies of the catalytic HKD motif (Fig 2; ). Interestingly, we found that three arginine residues in this motif involved in binding of PIP2 are shared also by plant members of the PXPH-PLD subfamily (Fig. 4B). In contrast, two of these residues are replaced by non-conserved amino acids in C2-PLDs, suggesting that this subfamily might adopt different mechanism of interaction with PIP2.
Two motifs rich in basic residues and allegedly similar to a polyphosphoinositide-binding motif from gelsolin or phospholipase C ([KR]X3-4KX [KR] [KR]) have been found in plant PLDs flanking the second catalytic HKD domain and proposed to mediate PIP2 binding by PLDs . It was claimed that all the basic residues are conserved only in Arabidopsis PLDβ1, whereas some are replaced with non-polar or acidic residues in AtPLDα1 and AtPLDγ1. However, our inspection of revised AtPLDβ1 sequence shows that the first "motif" of AtPLDβ1 has actually also only three basic residues, since the original motif definition was based on the inaccurate cDNA sequence. A genuine gelsolin/PLC consensus motif is found only in OsPLDβ1, OsPLDv1 and OsPLDv2 (Fig. 4A). The second proposed motif does not correspond exactly to the consensus of gelsolin family and PLC, but as if it was inverted (RKXRX4R). Not only that merely AtPLDβ1 out of all C2-PLDs retains all four basic residues (Fig. 4A), but also there is no experimental evidence that such a motif superficially resembling gelsolin/PLC consensus really binds phosphoinositides. Available experimental data [27–29] do not exclude the possibility that the C2 domain instead is responsible for regulatory effects of PIP2 on C2-PLDs.
Our analysis of Arabidopsis and rice genomic data complemented by searches of EST sequences revealed that plant PLDs are unexpectedly structurally diverse in two aspects.
First, individual plant genomes harbour various PLD types from both main PLD subfamilies. This is in sharp contrast to other large eukaryotic lineages. Species with completely or almost completely sequenced genomes, i.e. Saccharomyces cervisiae, Schizosaccharomyces pombe, Caenorhabditis elegans and Drosophila melanogaster, all posses only one gene from the PC-PLD family, and mammalian diversity is perhaps limited to two thoroughly characterised isoforms (our findings and ). All characterised animal and fungal PLDs belong to the PXPH-PLD subfamily. The occurrence of C2-PLDs beyond plants is unsure, they have been described only from angiosperms, and mosses are the most remote group for which C2-PLD sequences can be reliably found in databases (at least four distinct genes in Physcomitrella patens, see Additional files). It can be inferred from our phylogenetic analysis (see Fig. 1A) that the last common ancestor of animal, fungal and plant clades did harbour at least one gene from the PXPH-PLD subfamily as well as a gene that gave rise to the plant C2-PLDs. The latter PLD did not have to posses the C2 domain, this could be acquired later during evolution. In any case, the true C2-PLDs or their predecessors have been lost from the lineage leading to animals and fungi.
The second aspect of plant PLD diversity relates to inter-specific differences in the repertoire of distinct PLD types. For instance, there are no Arabidopsis orthologs of rice OsPLDη, OsPLDθ or OsPLDκ, while rice may lack counterparts of AtPLDγ or AtPLDζ from Arabidopsis. Similarly, only LePLDα1 from tomato is a true ortholog of other dicotyledonous alpha PLDs, while LePLDα2 and LePLDα3 form together a separate lineage within the PLDα cluster (Fig. 1B). Multiple independent losses of distinct PLD types must have occurred in individual lineages of angiosperms (e.g. PLDη or PLDκ lost in the lineage leading to Arabidopsis, PLDγ disappeared from the lineage toward rice, see Fig. 1A,B). On the other hand, independent multiplication within individual genomes seems to be common as well, exemplified by pairs of PLDα or PLDβ in Arabidopsis and rice (Fig. 1A,B; Table 1 and Table 2). Two main mechanisms for gene multiplication have apparently contributed to the diversity of plant PLDs. Origin of the two PLDα isoforms in Arabidopsis can be accounted for by polyploidisation, while the PLDγ triplet in Arabidopsis and the rice PLDη cluster have probably arisen by a non-reciprocal crossing-over (see above).
Diversity of plant PLDs raises the question of functional specificities of individual isoforms. Although only limited functional predictions can be made solely on the basis of sequence data, the principal difference in domain structure between C2- and PXPH-PLDs suggests that their cellular functions will also differ. PXPH-PLDs in animals and yeasts appear to be involved in regulation of vesicular and membrane trafficking (reviewed in ), and plant orthologs could be used in a similar context . Ca2+-independent PLD activity, which is probably exhibited by all PXPH-PLDs, has not been reported from plant tissues, but this is perhaps due to overabundant activity of C2-PLDs (especially PLDα) and to the notably low abundance of regulatory enzymes in general. Moreover, some stimulatory factors might be necessary for measurable activity of plant PXPH-PLDs, as is the case for mammalian PLD1 .
On the other hand, C2-PLDs may fulfil plant-specific tasks. Evolutionary dynamics of this subfamily in angiosperms indicates that environmental factors might exert big influence on these enzymes. Recognised role for C2-PLDs in processes such as response to wounding, pathogen attack and multiple abiotic stresses seems to fit this view, but other processes including membrane degradation during senescence also have to be considered . Functioning in signalling cascades may be common to both C2- and PXPH-PLDs, although the distinction between signalling function and the previously suggested roles does not have to be unambiguous.
Directions for future research on the plant PLD are straightforward. Besides the routinely used biochemical or pharmacological approaches, methods of reverse genetics (including anti-sense silencing and screening for insertional mutants) have to be employed. Partial functional redundancy, which can be expected for some plant PLD isoforms, could be coped with by generation of multiple mutants, accompanied by monitoring of expression of individual genes upon various circumstances and by experimental analysis of promoters. For deeper understanding of PLD regulation and interconnections within cellular context, attention must be focused on possible posttranslation modifications and interacting partners. Coordination of all these approaches has the potential to answer the question why plants farm so many PLDs.
Materials and Methods
For searches of public data we used BLAST toolkit at the National Centre for Biotechnology Information (http://www.ncbi.nlm.nih.gov/BLAST; [33, 34]). Searches were done in parallel with low complexity filter on and off, respectively, other parameters were kept default. All sequence databases containing plant data were exploited, including non-redundant nucleotide database, HTGS, GSS and EST databases and the non-redundant protein database. The final check of these databases was done between January 24, 2002, and January 26, 2002. Rice sequence data generated by Monsanto were searched using BLASTN and TBLAST facility at the rice-research.org web-page http://www.rice-research.org/. Hits from all BLAST searches with E-value above 0.1 were not considered for further analysis.
Multiple alignments were constructed by CLUSTALW (version 1.8) at the BCM Search Launcher (http://searchlauncher.bcm.tmc.edu/multi-align/multi-align.html; ), with default parameters. Manual editing of the alignments was done with the assistance of GENEDOC (Free Software Foundation, Inc.). Alternatively, multiple alignments were constructed using MACAW , with PAM 120 matrix used for protein sequences.
Exon-intron structures of Arabidopsis PLD genes were predicted employing GENESCAN (http://genes.mit.edu/GENSCAN.html; ), GRAIL http://grail.lsd.ornl.gov/Grail-1.3/, NetGene2 http://www.cbs.dtu.dk/services/NetGene2/, FGENEP http://dot.imgen.bcm.tmc.edu:9331/gene-finder/gf.html and SplicePredictor http://bioinformatics.iastate.edu/cgi-bin/sp.cgi. Models proposed by each programme were compared and final structures were proposed relating to the information from cognate cDNAs, ESTs and multiple alignments. Rice genes were predicted manually with the assistance of GENESCAN (with the options set up for maize). These predictions were again checked by comparison with ESTs, cDNAs and protein sequences of PLDs.
Specific domains in PLD protein sequences were searched by SMART (http://smart.embl-heidelberg.de/; ) and the Search Pfam tool (http://www.sanger.ac.uk/Software/Pfam/; ). Searches for targeting signals were performed using the TargetP programme http://www.cbs.dtu.dk/services/TargetP/.
Phylogenetic trees were inferred from multiple alignments of protein sequences using appropriate programmes from the PHYLIP package, version 3.57c . Neighbour-joining trees were constructed as described previously , PROTPARS programme was employed for maximum parsimony methods and confidence of the tree topology was estimated from 500 bootstrap replications. In the case that multiple alignments generated by CLUSTALW were used for phylogenetic inference, regions that could not be aligned unambiguously or containing deletions/insertions had been removed prior. For phylogenetic inference from alignments generated by MACAW only the most conserved boxes were used.
Levels of sequence identity/similarity occasionally noted through the text refer to values calculated by the BLAST 2 Sequences programme  with the low complexity filter off.
Note added in proof
A reannotation of the Arabidopsis genome released into GenBank after submission of the manuscript removes some inaccuracies in predictions of exon-intron structures of PLD genes independently uncovered also by our analysis. An updated list of Arabidopsis PLD genes has been deposited into the TAIR gene families database (http://www.arabidopsis.org/info/genefamily.html).
Liscovitch M, Czarny M, Fiucci G, Tang X: Phospholipase D: molecular and cell biology of a novel gene family. Biochem. 2000, 345: 401-415. 10.1042/0264-6021:3450401.
Wang X: Multiple forms of phospholipase D in plants: the gene family, catalytic and regulatory properties, and cellular functions. Prog Lipid Res. 2000, 39: 109-149. 10.1016/S0163-7827(00)00002-3.
Hanahan DJ, Chaikoff IL: A new phospholipide-splitting enzyme specific for the ester linkage between the nitrogenous base and the phosphoric acid grouping. J Biol Chem. 1947, 169: 699-705.
Wang X, Xu L, Zheng L: Cloning and expression of phosphatidylcholine-hydrolyzing phospholipase D from Ricinus communis L. J Biol Chem. 1994, 269: 20312-20317.
den Hartog M, Musgrave A, Munnik T: Nod factor-induced phosphatidic acid and diacylglycerol pyrophosphate formation: a role for phospholipase C and D in root hair deformation. Plant J. 2001, 25: 55-65. 10.1046/j.1365-313X.2001.00931.x.
Munnik T, Irvine RF, Musgrave A: Phospholipid signalling in plants. Biochim Biophys Acta. 1998, 1389: 222-272. 10.1016/S0005-2760(97)00158-6.
Lein W, Saalbach G: Cloning and direct G-protein regulation of phospholipase D from tobacco. Biochim Biophys Acta. 2001, 1530: 172-183. 10.1016/S1388-1981(00)00182-7.
Chapman KD: Emerging physiological roles for N-acylphosphatidylethanolamine metabolism in plants: signal transduction and membrane protection. Chem Phys Lipids. 2000, 108: 221-229. 10.1016/S0009-3084(00)00198-5.
Munnik T: Phosphatidic acid: an emerging plant lipid second messenger. Trends Plant Sci. 2001, 6: 227-233. 10.1016/S1360-1385(01)01918-5.
Fan L, Zheng S, Wang X: Antisense suppression of phospholipase D alpha retards abscisic acid- and ethylene-promoted senescence of postharvest Arabidopsis leaves. Plant Cell. 1997, 9: 2183-2196. 10.1105/tpc.9.12.2183.
Pappan K, Qin W, Dyer JH, Zheng L, Wang X: Molecular cloning and functional analysis of polyphosphoinositide-dependent phospholipase D, PLDbeta, from Arabidopsis. J Biol Chem. 1997, 272: 7055-7061. 10.1074/jbc.272.11.7055.
AGI: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature. 2000, 408: 796-815. 10.1038/35048692.
Barry GF: The use of the Monsanto draft rice genome sequence in research. Plant Physiol. 2001, 125: 1164-1165. 10.1104/pp.125.3.1164.
Ponting CP, Kerr ID: A novel family of phospholipase D homologues that includes phospholipid synthases and putative endonucleases: Identification of duplicated repeats and potential active site residues. Protein Sci. 1996, 5: 914-922.
Katagiri T, Takahashi S, Shinozaki S: Involvement of a novel Arabidopsis phospholipase D, AtPLDδ, in dehydration-inducible accumulation of phosphatidic acid in stress signalling. Plant J. 2001, 26: 595-606. 10.1046/j.1365-313X.2001.01060.x.
Wang C, Wang X: A Novel Phospholipase D of Arabidopsis That Is Activated by Oleic Acid and Associated with the Plasma Membrane. Plant Physiol. 2001, 127: 1102-12. 10.1104/pp.127.3.1102.
Wang X: Plant phospholipases. Annu Rev Plant Physiol Plant Mol Biol. 2001, 52: 211-231. 10.1146/annurev.arplant.52.1.211.
Qin W, Dyer JH, Zheng Z, Wang X: Isolation and nucleotide sequence of the fourth phospholipase D (accession no. Af138281), PLD-gamma 2, from Arabidopsis thaliana. (PGR99-084). Plant Physiol. 1999, 120: 635-
Laxalt AM, ter Riet B, Verdonk JC, Parigi L, Tameling WIL, Vossen J, Haring M, Musgrave A, Munnik T: Characterization of five tomato phospholipase D cDNAs: rapid and specific expression of LePLDbeta1 on elicitation with xylanase. Plant. 2001, 26: 237-247. 10.1046/j.1365-313X.2001.01023.x.
Ueki J, Morioka S, Komari T, Kumashiro T: Purification and characterization of phospholipase-D (PLD) from rice (Oryza-sativa L) and cloning of cDNA for PLD from rice and maize (Zea-mays L). Plant Cell Physiol. 1995, 36: 903-914.
Morioka S, Ueki J, Komari T: Characterization of two distinctive genomic clones (accession nos. AB001919 and AB001920) for phospholipase D from rice (PGR97-076). Plant Physiol. 1997, 114: 396-
Almquist KC, Paliyath G: Cloning and sequencing of a full-length cDNA coding for phospholipase D alpha (Accession Number 201661) from Lycopersicon esculentum Mill. (PGR00-013). Plant Physiol. 2000, 122: 292-
Whitaker BD, Smith DL, Green KC: Cloning, characterization and functional expression of a phospholipase D alpha cDNA from tomato fruit. Physiol Plant. 2001, 112: 87-94. 10.1034/j.1399-3054.2001.1120112.x.
Ewing RM, Ben Kahla A, Poirot O, Lopez F, Audic S, Claverie JM: Large-scale statistical analyses of rice ESTs reveal correlated patterns of gene expression. Genome Res. 1999, 9: 950-959. 10.1101/gr.9.10.950.
Frank W, Munnik T, Kerkmann K, Salamini F, Bartels D: Water deficit triggers phospholipase D activity in the resurrection plant Craterostigma plantagineum. Plant Cell. 2000, 12: 111-124. 10.1105/tpc.12.1.111.
Rizo J, Sudhof TC: C2-domains, structure and function of a universal Ca2+-binding domain. J Biol Chem. 1998, 273: 15879-15882. 10.1074/jbc.273.26.15879.
Qin W, Pappan K, Wang X: Molecular heterogeneity of phospholipase D (PLD). Cloning of PLDgamma and regulation of plant PLDgamma, -beta, and -alpha by polyphosphoinositides and calcium. J Biol Chem. 1997, 272: 28267-28273. 10.1074/jbc.272.45.28267.
Zheng L, Krishnamoorthi R, Zolkiewski M, Wang X: Distinct Ca2+ binding properties of novel C2 domains of plant phospholipase dalpha and beta. J Biol Chem. 2000, 275: 19700-19706. 10.1074/jbc.M001945200.
Pappan K, Wang X: Plant phospholipase Dalpha is an acidic phospholipase active at near-physiological Ca(2+) concentrations. Arch Biochem Biophys. 1999, 368: 347-353. 10.1006/abbi.1999.1325.
Hodgkin MN, Masson MR, Powner D, Saqib KM, Ponting CP, Wakelam MJ: Phospholipase D regulation and localisation is dependent upon a phosphatidylinositol 4,5-biphosphate-specific PH domain. Curr Biol. 2000, 10: 43-46. 10.1016/S0960-9822(99)00264-X.
Sciorra VA, Rudge SA, Prestwich GD, Frohman MA, Engebrecht J, Morris AJ: Identification of a phosphoinositide binding motif that mediates activation of mammalian and yeast phospholipase D isoenzymes. EMBO J. 1999, 18: 5911-5921. 10.1093/emboj/18.21.5911.
Munnik T, Musgrave A: Phospholipid signaling in plants: holding on to phospholipase D. Sci STKE. 2001, 111: PE42-
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410. 10.1006/jmbi.1990.9999.
Altschul SF, Madden TL, Schaffer AA, Zhang JH, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
Thompson JD, Higgins DG, Gibson TJ: Clustal-W – improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680.
Schuler GD, Altschul SF, Lipman DJ: A workbench for multiple alignment construction and analysis. Prot Struct Fund Genet. 1991, 9: 180-190.
Burge C, Karlin S: Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997, 268: 78-94. 10.1006/jmbi.1997.0951.
Schultz J, Copley RR, Doerks T, Ponting CP, Bork P: SMART: a web-based tool for the study of genetically mobile domains. Nucleic Acids Res. 2000, 28: 231-234. 10.1093/nar/28.1.231.
Bateman A, Birney E, Durbin R, Eddy SR, Howe KL, Sonnhammer ELL: The Pfam protein families database. Nucleic Acids Res. 2000, 28: 263-266. 10.1093/nar/28.1.263.
Felsenstein J: Inferring phylogenies from protein sequences by parsimony, distance, and likelihood methods. Methods Enzymol. 1996, 266: 418-427.
Cvrckova F: Are plant formins integral membrane proteins?. Genome Biology. 2000, 1: 1-001.
Tatusova TA, Madden TL: Blast 2 Sequences, a new tool for comparing protein and nucleotide sequences. FEMS Microbiology Letters. 1999, 174: 247-250. 10.1016/S0378-1097(99)00149-4.
Dyer JH, Zheng L, Wang X: Cloning and nucleotide sequence of a cDNA (Accession No. U36381) encoding phospholipase D from Arabidopsis. (PGR95-096). Plant Physiol. 1995, 109: 1497-10.1104/pp.109.4.1497.
Gardiner JC, Harper JD, Weerakoon ND, Collings DA, Ritchie S, Gilroy S, Cyr RJ, Marc J: A 90-kd phospholipase d from tobacco binds to microtubules and the plasma membrane. Plant Cell. 2001, 13: 2143-58. 10.1105/tpc.13.9.2143.
Nalefski AE, Falke JJ: The C2 domain calcium-binding motif: Structural and functional diversity. Protein Sci. 1996, 12: 2375-2390.
This work was supported by MSMT CR LN00A081, J13/98:113100003, and GACR 206/99/1138 grants.
We are highly indebted to Monsanto for being allowed to search their rice genome data and to make the results publicly available. We thank Marta Čadyová for perfect technical assistance.
Electronic supplementary material
About this article
Cite this article
Eliáš, M., Potocký, M., Cvrčková, F. et al. Molecular diversity of phospholipase D in angiosperms. BMC Genomics 3, 2 (2002). https://doi.org/10.1186/1471-2164-3-2
- Medicago Truncatula
- International Rice Genome Sequencing Project
- Resurrection Plant Craterostigma Plantagineum
- Plant PLDs
- Tandem Triplication