Biosynthesis and function of 7-deazaguanine derivatives in bacteria and phages

SUMMARYDeazaguanine modifications play multifaceted roles in the molecular biology of DNA and tRNA, shaping diverse yet essential biological processes, including the nuanced fine-tuning of translation efficiency and the intricate modulation of codon-anticodon interactions. Beyond their roles in translation, deazaguanine modifications contribute to cellular stress resistance, self-nonself discrimination mechanisms, and host evasion defenses, directly modulating the adaptability of living organisms. Deazaguanine moieties extend beyond nucleic acid modifications, manifesting in the structural diversity of biologically active natural products. Their roles in fundamental cellular processes and their presence in biologically active natural products underscore their versatility and pivotal contributions to the intricate web of molecular interactions within living organisms. Here, we discuss the current understanding of the biosynthesis and multifaceted functions of deazaguanines, shedding light on their diverse and dynamic roles in the molecular landscape of life.

Epoxyqueuosine Reductase QueH in the Biosynthetic Pathway to tRNA Queuosine Is a Unique Metalloenzyme

Queuosine is a structurally unique and functionally important tRNA modification, widely distributed in eukaryotes and bacteria. The final step of queuosine biosynthesis is the reduction/deoxygenation of epoxyqueuosine to form the cyclopentene motif of the nucleobase. The chemistry is performed by the structurally and functionally characterized cobalamin-dependent QueG. However, the queG gene is absent from several bacteria that otherwise retain queuosine biosynthesis machinery. Members of the IPR003828 family (previously known as DUF208) have been recently identified as nonorthologous replacements of QueG, and this family was renamed QueH. Here, we present the structural characterization of QueH from Thermotoga maritima. The structure reveals an unusual active site architecture with a [4Fe-4S] metallocluster along with an adjacent coordinated iron metal. The juxtaposition of the cofactor and coordinated metal ion predicts a unique mechanism for a two-electron reduction/deoxygenation of epoxyqueuosine. To support the structural characterization, in vitro biochemical and genomic analyses are presented. Overall, this work reveals new diversity in the chemistry of iron/sulfur-dependent enzymes and novel insight into the last step of this widely conserved tRNA modification.

The EFI Web Resource for Genomic Enzymology Tools: Leveraging Protein, Genome, and Metagenome Databases to Discover Novel Enzymes and Metabolic Pathways

The assignment of functions to uncharacterized proteins discovered in genome projects requires easily accessible tools and computational resources for large-scale, user-friendly leveraging of the protein, genome, and metagenome databases by experimentalists. This article describes the web resource developed by the Enzyme Function Initiative (EFI; accessed at https://efi.igb.illinois.edu/ ) that provides "genomic enzymology" tools ("web tools") for (1) generating sequence similarity networks (SSNs) for protein families (EFI-EST); (2) analyzing and visualizing genome context of the proteins in clusters in SSNs (in genome neighborhood networks, GNNs, and genome neighborhood diagrams, GNDs) (EFI-GNT); and (3) prioritizing uncharacterized SSN clusters for functional assignment based on metagenome abundance (chemically guided functional profiling, CGFP) (EFI-CGFP). The SSNs generated by EFI-EST are used as the input for EFI-GNT and EFI-CGFP, enabling easy transfer of information among the tools. The networks are visualized and analyzed using Cytoscape, a widely used desktop application; GNDs and CGFP heatmaps summarizing metagenome abundance are viewed within the tools. We provide a detailed example of the integrated use of the tools with an analysis of glycyl radical enzyme superfamily (IPR004184) found in the human gut microbiome. This analysis demonstrates that (1) SwissProt annotations are not always correct, (2) large-scale genome context analyses allow the prediction of novel metabolic pathways, and (3) metagenome abundance can be used to identify/prioritize uncharacterized proteins for functional investigation.

Discovery of novel bacterial queuine salvage enzymes and pathways in human pathogens

Queuosine (Q) is a complex tRNA modification widespread in eukaryotes and bacteria that contributes to the efficiency and accuracy of protein synthesis. Eukaryotes are not capable of Q synthesis and rely on salvage of the queuine base (q) as a Q precursor. While many bacteria are capable of Q de novo synthesis, salvage of the prokaryotic Q precursors preQ0 and preQ1 also occurs. With the exception of Escherichia coli YhhQ, shown to transport preQ0 and preQ1, the enzymes and transporters involved in Q salvage and recycling have not been well described. We discovered and characterized 2 Q salvage pathways present in many pathogenic and commensal bacteria. The first, found in the intracellular pathogen Chlamydia trachomatis, uses YhhQ and tRNA guanine transglycosylase (TGT) homologs that have changed substrate specificities to directly salvage q, mimicking the eukaryotic pathway. The second, found in bacteria from the gut flora such as Clostridioides difficile, salvages preQ1 from q through an unprecedented reaction catalyzed by a newly defined subgroup of the radical-SAM enzyme family. The source of q can be external through transport by members of the energy-coupling factor (ECF) family or internal through hydrolysis of Q by a dedicated nucleosidase. This work reinforces the concept that hosts and members of their associated microbiota compete for the salvage of Q precursors micronutrients.

Gene Graphics: a genomic neighborhood data visualization web application

Summary

The examination of gene neighborhood is an integral part of comparative genomics but no tools to produce publication quality graphics of gene clusters are available. Gene Graphics is a straightforward web application for creating such visuals. Supported inputs include National Center for Biotechnology Information gene and protein identifiers with automatic fetching of neighboring information, GenBank files and data extracted from the SEED database. Gene representations can be customized for many parameters including gene and genome names, colors and sizes. Gene attributes can be copied and pasted for rapid and user-friendly customization of homologous genes between species. In addition to Portable Network Graphics and Scalable Vector Graphics, produced representations can be exported as Tagged Image File Format or Encapsulated PostScript, formats that are standard for publication. Hands-on tutorials with real life examples inspired from publications are available for training.

Availability and implementation

Gene Graphics is freely available at https://katlabs.cc/genegraphics/ and source code is hosted at https://github.com/katlabs/genegraphics.

Contact

katherinejh@ufl.edu or remizallot@ufl.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

'Democratized' genomic enzymology web tools for functional assignment

The protein databases contain an exponentially growing number of sequences as a result of the recent increase in ease and decrease in cost of genome sequencing. The rate of data accumulation far exceeds the rate of functional studies, producing an increase in genomic 'dark matter', sequences for which no precise and validated function is defined. Publicly accessible, that is 'democratized,' genomic enzymology web tools are essential to leverage the protein and genome databases for discovery of the in vitro activities and in vivo functions of novel enzymes and proteins belonging to the dark matter. In this review, we discuss the use of web tools that have proven successful for functional assignment. We also describe a mechanism for ensuring the capture of published functional data so that the quality of both curated and automated annotations transfer can be improved.

Identification of a Novel Epoxyqueuosine Reductase Family by Comparative Genomics

The reduction of epoxyqueuosine (oQ) is the last step in the synthesis of the tRNA modification queuosine (Q). While the epoxyqueuosine reductase (EC 1.17.99.6) enzymatic activity was first described 30 years ago, the encoding gene queG was only identified in Escherichia coli in 2011. Interestingly, queG is absent from a large number of sequenced genomes that harbor Q synthesis or salvage genes, suggesting the existence of an alternative epoxyqueuosine reductase in these organisms. By analyzing phylogenetic distributions, physical gene clustering, and fusions, members of the Domain of Unknown Function 208 (DUF208) family were predicted to encode for an alternative epoxyqueuosine reductase. This prediction was validated with genetic methods. The Q modification is present in Lactobacillus salivarius, an organism missing queG but harboring the duf208 gene. Acinetobacter baylyi ADP1 is one of the few organisms that harbor both QueG and DUF208, and deletion of both corresponding genes was required to observe the absence of Q and the accumulation of oQ in tRNA. Finally, the conversion oQ to Q was restored in an E. coli queG mutant by complementation with plasmids harboring duf208 genes from different bacteria. Members of the DUF208 family are not homologous to QueG enzymes, and thus, duf208 is a non-orthologous replacement of queG. We propose to name DUF208 encoding genes as queH. While QueH contains conserved cysteines that could be involved in the coordination of a Fe/S center in a similar fashion to what has been identified in QueG, no cobalamin was identified associated with recombinant QueH protein.

The Escherichia coli COG1738 Member YhhQ Is Involved in 7-Cyanodeazaguanine (preQ₀) Transport

Queuosine (Q) is a complex modification of the wobble base in tRNAs with GUN anticodons. The full Q biosynthesis pathway has been elucidated in Escherichia coli. FolE, QueD, QueE and QueC are involved in the conversion of guanosine triphosphate (GTP) to 7-cyano-7-deazaguanine (preQ₀), an intermediate of increasing interest for its central role in tRNA and DNA modification and secondary metabolism. QueF then reduces preQ₀ to 7-aminomethyl-7-deazaguanine (preQ₁). PreQ₁ is inserted into tRNAs by tRNA guanine(34) transglycosylase (TGT). The inserted base preQ₁ is finally matured to Q by two additional steps involving QueA and QueG or QueH. Most Eubacteria harbor the full set of Q synthesis genes and are predicted to synthesize Q de novo. However, some bacteria only encode enzymes involved in the second half of the pathway downstream of preQ₀ synthesis, including the signature enzyme TGT. Different patterns of distribution of the queF, tgt, queA and queG or queH genes are observed, suggesting preQ₀, preQ₁ or even the queuine base being salvaged in specific organisms. Such salvage pathways require the existence of specific 7-deazapurine transporters that have yet to be identified. The COG1738 family was identified as a candidate for a missing preQ₀/preQ₁ transporter in prokaryotes, by comparative genomics analyses. The existence of Q precursor salvage was confirmed for the first time in bacteria, in vivo, through an indirect assay. The involvement of the COG1738 in salvage of a Q precursor was experimentally validated in Escherichia coli, where it was shown that the COG1738 family member YhhQ is essential for preQ₀ transport.

Arabidopsis TH2 Encodes the Orphan Enzyme Thiamin Monophosphate Phosphatase

To synthesize the cofactor thiamin diphosphate (ThDP), plants must first hydrolyze thiamin monophosphate (ThMP) to thiamin, but dedicated enzymes for this hydrolysis step were unknown and widely doubted to exist. The classical thiamin-requiring th2-1 mutation in Arabidopsis thaliana was shown to reduce ThDP levels by half and to increase ThMP levels 5-fold, implying that the THIAMIN REQUIRING2 (TH2) gene product could be a dedicated ThMP phosphatase. Genomic and transcriptomic data indicated that TH2 corresponds to At5g32470, encoding a HAD (haloacid dehalogenase) family phosphatase fused to a TenA (thiamin salvage) family protein. Like the th2-1 mutant, an insertional mutant of At5g32470 accumulated ThMP, and the thiamin requirement of the th2-1 mutant was complemented by wild-type At5g32470 Complementation tests in Escherichia coli and enzyme assays with recombinant proteins confirmed that At5g32470 and its maize (Zea mays) orthologs GRMZM2G148896 and GRMZM2G078283 are ThMP-selective phosphatases whose activity resides in the HAD domain and that the At5g32470 TenA domain has the expected thiamin salvage activity. In vitro and in vivo experiments showed that alternative translation start sites direct the At5g32470 protein to the cytosol and potentially also to mitochondria. Our findings establish that plants have a dedicated ThMP phosphatase and indicate that modest (50%) ThDP depletion can produce severe deficiency symptoms.

Functional Annotations of Paralogs: A Blessing and a Curse

Gene duplication followed by mutation is a classic mechanism of neofunctionalization, producing gene families with functional diversity. In some cases, a single point mutation is sufficient to change the substrate specificity and/or the chemistry performed by an enzyme, making it difficult to accurately separate enzymes with identical functions from homologs with different functions. Because sequence similarity is often used as a basis for assigning functional annotations to genes, non-isofunctional gene families pose a great challenge for genome annotation pipelines. Here we describe how integrating evolutionary and functional information such as genome context, phylogeny, metabolic reconstruction and signature motifs may be required to correctly annotate multifunctional families. These integrative analyses can also lead to the discovery of novel gene functions, as hints from specific subgroups can guide the functional characterization of other members of the family. We demonstrate how careful manual curation processes using comparative genomics can disambiguate subgroups within large multifunctional families and discover their functions. We present the COG0720 protein family as a case study. We also discuss strategies to automate this process to improve the accuracy of genome functional annotation pipelines.

Systematic identification and analysis of frequent gene fusion events in metabolic pathways

Background

Gene fusions are the most powerful type of in silico-derived functional associations. However, many fusion compilations were made when <100 genomes were available, and algorithms for identifying fusions need updating to handle the current avalanche of sequenced genomes. The availability of a large fusion dataset would help probe functional associations and enable systematic analysis of where and why fusion events occur.

Results

Here we present a systematic analysis of fusions in prokaryotes. We manually generated two training sets: (i) 121 fusions in the model organism Escherichia coli; (ii) 131 fusions found in B vitamin metabolism. These sets were used to develop a fusion prediction algorithm that captured the training set fusions with only 7 % false negatives and 50 % false positives, a substantial improvement over existing approaches. This algorithm was then applied to identify 3.8 million potential fusions across 11,473 genomes. The results of the analysis are available in a searchable database at http://modelseed.org/projects/fusions/. A functional analysis identified 3,000 reactions associated with frequent fusion events and revealed areas of metabolism where fusions are particularly prevalent.

Conclusions

Customary definitions of fusions were shown to be ambiguous, and a stricter one was proposed. Exploring the genes participating in fusion events showed that they most commonly encode transporters, regulators, and metabolic enzymes. The major rationales for fusions between metabolic genes appear to be overcoming pathway bottlenecks, avoiding toxicity, controlling competing pathways, and facilitating expression and assembly of protein complexes. Finally, our fusion dataset provides powerful clues to decipher the biological activities of domains of unknown function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2782-3) contains supplementary material, which is available to authorized users.

Functional Diversity of Haloacid Dehalogenase Superfamily Phosphatases from Saccharomyces cerevisiae: BIOCHEMICAL, STRUCTURAL, AND EVOLUTIONARY INSIGHTS

The haloacid dehalogenase (HAD)-like enzymes comprise a large superfamily of phosphohydrolases present in all organisms. The Saccharomyces cerevisiae genome encodes at least 19 soluble HADs, including 10 uncharacterized proteins. Here, we biochemically characterized 13 yeast phosphatases from the HAD superfamily, which includes both specific and promiscuous enzymes active against various phosphorylated metabolites and peptides with several HADs implicated in detoxification of phosphorylated compounds and pseudouridine. The crystal structures of four yeast HADs provided insight into their active sites, whereas the structure of the YKR070W dimer in complex with substrate revealed a composite substrate-binding site. Although the S. cerevisiae and Escherichia coli HADs share low sequence similarities, the comparison of their substrate profiles revealed seven phosphatases with common preferred substrates. The cluster of secondary substrates supporting significant activity of both S. cerevisiae and E. coli HADs includes 28 common metabolites that appear to represent the pool of potential activities for the evolution of novel HAD phosphatases. Evolution of novel substrate specificities of HAD phosphatases shows no strict correlation with sequence divergence. Thus, evolution of the HAD superfamily combines the conservation of the overall substrate pool and the substrate profiles of some enzymes with remarkable biochemical and structural flexibility of other superfamily members.

Salvage of the thiamin pyrimidine moiety by plant TenA proteins lacking an active-site cysteine

The TenA protein family occurs in prokaryotes, plants and fungi; it has two subfamilies, one (TenA_C) having an active-site cysteine, the other (TenA_E) not. TenA_C proteins participate in thiamin salvage by hydrolysing the thiamin breakdown product amino-HMP (4-amino-5-aminomethyl-2-methylpyrimidine) to HMP (4-amino-5-hydroxymethyl-2-methylpyrimidine); the function of TenA_E proteins is unknown. Comparative analysis of prokaryote and plant genomes predicted that (i) TenA_E has a salvage role similar to, but not identical with, that of TenA_C and (ii) that TenA_E and TenA_C also have non-salvage roles since they occur in organisms that cannot make thiamin. Recombinant Arabidopsis and maize TenA_E proteins (At3g16990, GRMZM2G080501) hydrolysed amino-HMP to HMP and, far more actively, hydrolysed the N-formyl derivative of amino-HMP to amino-HMP. Ablating the At3g16990 gene in a line with a null mutation in the HMP biosynthesis gene ThiC prevented its rescue by amino-HMP. Ablating At3g16990 in the wild-type increased sensitivity to paraquat-induced oxidative stress; HMP overcame this increased sensitivity. Furthermore, the expression of TenA_E and ThiC genes in Arabidopsis and maize was inversely correlated. These results indicate that TenA_E proteins mediate amidohydrolase and aminohydrolase steps in the salvage of thiamin breakdown products. As such products can be toxic, TenA_E proteins may also pre-empt toxicity.

Plant, animal, and fungal micronutrient queuosine is salvaged by members of the DUF2419 protein family

Queuosine (Q) is a modification found at the wobble position of tRNAs with GUN anticodons. Although Q is present in most eukaryotes and bacteria, only bacteria can synthesize Q de novo. Eukaryotes acquire queuine (q), the free base of Q, from diet and/or microflora, making q an important but under-recognized micronutrient for plants, animals, and fungi. Eukaryotic type tRNA-guanine transglycosylases (eTGTs) are composed of a catalytic subunit (QTRT1) and a homologous accessory subunit (QTRTD1) forming a complex that catalyzes q insertion into target tRNAs. Phylogenetic analysis of eTGT subunits revealed a patchy distribution pattern in which gene losses occurred independently in different clades. Searches for genes co-distributing with eTGT family members identified DUF2419 as a potential Q salvage protein family. This prediction was experimentally validated in Schizosaccharomyces pombe by confirming that Q was present by analyzing tRNA^Asp with anticodon GUC purified from wild-type cells and by showing that Q was absent from strains carrying deletions in the QTRT1 or DUF2419 encoding genes. DUF2419 proteins occur in most Eukarya with a few possible cases of horizontal gene transfer to bacteria. The universality of the DUF2419 function was confirmed by complementing the S. pombe mutant with the Zea mays (maize), human, and Sphaerobacter thermophilus homologues. The enzymatic function of this family is yet to be determined, but structural similarity with DNA glycosidases suggests a ribonucleoside hydrolase activity.

Identification of the thiamin salvage enzyme thiazole kinase in Arabidopsis and maize

The breakdown of thiamin (vitamin B1) and its phosphates releases a thiazole moiety, 4-methyl-5-(2-hydroxyethyl)thiazole (THZ), that microorganisms and plants are able to salvage for re-use in thiamin synthesis. The salvage process starts with the ATP-dependent phosphorylation of THZ, which in bacteria is mediated by ThiM. The Arabidopsis and maize genomes encode homologs of ThiM (At3g24030 and GRMZM2G094558, respectively). Plasmid-driven expression of either plant homolog restored the ability of THZ to rescue Escherichia coli thiM deletant strains, showing that the plant proteins have ThiM activity in vivo. Enzymatic assays with purified recombinant proteins confirmed the presence of THZ kinase activity. Furthermore, ablating the Arabidopsis At3g24030 gene in a thiazole synthesis mutant severely impaired rescue by THZ. Collectively, these results show that ThiM homologs are the main source of THZ kinase activity in plants and are consequently crucial for thiamin salvage.

Identification of mitochondrial coenzyme a transporters from maize and Arabidopsis

Plants make coenzyme A (CoA) in the cytoplasm but use it for reactions in mitochondria, chloroplasts, and peroxisomes, implying that these organelles have CoA transporters. A plant peroxisomal CoA transporter is already known, but plant mitochondrial or chloroplastic CoA transporters are not. Mitochondrial CoA transporters belonging to the mitochondrial carrier family, however, have been identified in yeast (Saccharomyces cerevisiae; Leu-5p) and mammals (SLC25A42). Comparative genomic analysis indicated that angiosperms have two distinct homologs of these mitochondrial CoA transporters, whereas nonflowering plants have only one. The homologs from maize (Zea mays; GRMZM2G161299 and GRMZM2G420119) and Arabidopsis (Arabidopsis thaliana; At1g14560 and At4g26180) all complemented the growth defect of the yeast leu5Δ mitochondrial CoA carrier mutant and substantially restored its mitochondrial CoA level, confirming that these proteins have CoA transport activity. Dual-import assays with purified pea (Pisum sativum) mitochondria and chloroplasts, and subcellular localization of green fluorescent protein fusions in transiently transformed tobacco (Nicotiana tabacum) Bright Yellow-2 cells, showed that the maize and Arabidopsis proteins are targeted to mitochondria. Consistent with the ubiquitous importance of CoA, the maize and Arabidopsis mitochondrial CoA transporter genes are expressed at similar levels throughout the plant. These data show that representatives of both monocotyledons and eudicotyledons have twin, mitochondrially located mitochondrial carrier family carriers for CoA. The highly conserved nature of these carriers makes possible their reliable annotation in other angiosperm genomes.

Acyl-lipid metabolism

Acyl lipids in Arabidopsis and all other plants have a myriad of diverse functions. These include providing the core diffusion barrier of the membranes that separates cells and subcellular organelles. This function alone involves more than 10 membrane lipid classes, including the phospholipids, galactolipids, and sphingolipids, and within each class the variations in acyl chain composition expand the number of structures to several hundred possible molecular species. Acyl lipids in the form of triacylglycerol account for 35% of the weight of Arabidopsis seeds and represent their major form of carbon and energy storage. A layer of cutin and cuticular waxes that restricts the loss of water and provides protection from invasions by pathogens and other stresses covers the entire aerial surface of Arabidopsis. Similar functions are provided by suberin and its associated waxes that are localized in roots, seed coats, and abscission zones and are produced in response to wounding. This chapter focuses on the metabolic pathways that are associated with the biosynthesis and degradation of the acyl lipids mentioned above. These pathways, enzymes, and genes are also presented in detail in an associated website (ARALIP: http://aralip.plantbiology.msu.edu/). Protocols and methods used for analysis of Arabidopsis lipids are provided. Finally, a detailed summary of the composition of Arabidopsis lipids is provided in three figures and 15 tables.

Acyl-lipid metabolism

Acyl lipids in Arabidopsis and all other plants have a myriad of diverse functions. These include providing the core diffusion barrier of the membranes that separates cells and subcellular organelles. This function alone involves more than 10 membrane lipid classes, including the phospholipids, galactolipids, and sphingolipids, and within each class the variations in acyl chain composition expand the number of structures to several hundred possible molecular species. Acyl lipids in the form of triacylglycerol account for 35% of the weight of Arabidopsis seeds and represent their major form of carbon and energy storage. A layer of cutin and cuticular waxes that restricts the loss of water and provides protection from invasions by pathogens and other stresses covers the entire aerial surface of Arabidopsis. Similar functions are provided by suberin and its associated waxes that are localized in roots, seed coats, and abscission zones and are produced in response to wounding. This chapter focuses on the metabolic pathways that are associated with the biosynthesis and degradation of the acyl lipids mentioned above. These pathways, enzymes, and genes are also presented in detail in an associated website (ARALIP: http://aralip.plantbiology.msu.edu/). Protocols and methods used for analysis of Arabidopsis lipids are provided. Finally, a detailed summary of the composition of Arabidopsis lipids is provided in three figures and 15 tables.

Acyl-lipid metabolism

Acyl lipids in Arabidopsis and all other plants have a myriad of diverse functions. These include providing the core diffusion barrier of the membranes that separates cells and subcellular organelles. This function alone involves more than 10 membrane lipid classes, including the phospholipids, galactolipids, and sphingolipids, and within each class the variations in acyl chain composition expand the number of structures to several hundred possible molecular species. Acyl lipids in the form of triacylglycerol account for 35% of the weight of Arabidopsis seeds and represent their major form of carbon and energy storage. A layer of cutin and cuticular waxes that restricts the loss of water and provides protection from invasions by pathogens and other stresses covers the entire aerial surface of Arabidopsis. Similar functions are provided by suberin and its associated waxes that are localized in roots, seed coats, and abscission zones and are produced in response to wounding. This chapter focuses on the metabolic pathways that are associated with the biosynthesis and degradation of the acyl lipids mentioned above. These pathways, enzymes, and genes are also presented in detail in an associated website (ARALIP: http://aralip.plantbiology.msu.edu/). Protocols and methods used for analysis of Arabidopsis lipids are provided. Finally, a detailed summary of the composition of Arabidopsis lipids is provided in three figures and 15 tables.

Polyphosphoinositides are enriched in plant membrane rafts and form microdomains in the plasma membrane

In this article, we analyzed the lipid composition of detergent-insoluble membranes (DIMs) purified from tobacco (Nicotiana tabacum) plasma membrane (PM), focusing on polyphosphoinositides, lipids known to be involved in various signal transduction events. Polyphosphoinositides were enriched in DIMs compared with whole PM, whereas all structural phospholipids were largely depleted from this fraction. Fatty acid composition analyses suggest that enrichment of polyphosphoinositides in DIMs is accompanied by their association with more saturated fatty acids. Using an immunogold-electron microscopy strategy, we were able to visualize domains of phosphatidylinositol 4,5-bisphosphate in the plane of the PM, with 60% of the epitope found in clusters of approximately 25 nm in diameter and 40% randomly distributed at the surface of the PM. Interestingly, the phosphatidylinositol 4,5-bisphosphate cluster formation was not significantly sensitive to sterol depletion induced by methyl-beta-cyclodextrin. Finally, we measured the activities of various enzymes of polyphosphoinositide metabolism in DIMs and PM and showed that these activities are present in the DIM fraction but not enriched. The putative role of plant membrane rafts as signaling membrane domains or membrane-docking platforms is discussed.