BRIEF RESEARCH REPORT article

Front. Immunol., 10 August 2026

Sec. Comparative Immunology

Volume 17 - 2026 | https://doi.org/10.3389/fimmu.2026.1838657

An evolutionary analysis of lysophospholipid acyltransferases generating membrane phospholipid diversity

  • 1. Department of Lipid Life Science, National Institute of Global Health and Medicine, Japan Institute for Health Security (JIHS), Shinjuku-ku, Tokyo, Japan

  • 2. Department of Medical Lipid Science, Graduate School of Medicine, The University of Tokyo, Bunkyo-ku, Tokyo, Japan

  • 3. Human Population Genomics laboratory, Department of Informatics, National Institute of Genetics. Yata, Mishima, Shizuoka, Japan

  • 4. Genome Medical Science Project, National Institute of Global Health and Medicine, Japan Institute for Health Security (JIHS), Shinjuku-ku, Tokyo, Japan

Abstract

All living organisms are enclosed by biological membranes composed primarily of phospholipids. The diversification of membrane lipids has shaped key membrane functions, including signaling, membrane protein regulation, and immune responses, and may also have driven organismal evolution. In humans, this diversity is mainly generated by lysophospholipid acyltransferases (LPLATs), yet the evolutionary distribution and diversity of LPLATs across organisms remain poorly understood, particularly in pathogens. Here, we systematically identified candidate LPLAT proteins belonging to two major families in representative prokaryotic and eukaryotic species and reconstructed their phylogenetic relationships. The analysis reveals a marked expansion and diversification of LPLATs in vertebrates and suggests the possibility of pathogen-specific LPLAT expansions. These findings suggest that the diversification of membrane phospholipid remodeling enzymes may accompany increasing biological complexity and pathogenicity. This study provides an evolutionary framework for understanding phospholipid metabolism and highlights pathogen LPLATs as potential targets for exploring host defense mechanisms.

1 Introduction

Biological membranes define the boundaries of cells and organelles and play central roles in numerous cellular processes, including signaling, membrane protein regulation, and immune responses (–). The principal structural components of these membranes are glycerophospholipids (hereafter referred to as phospholipids). Beyond their structural role in membrane bilayers, phospholipids can also be released as main components of alveolar surfactant and lipoproteins (, ). Mammalian phospholipids form a highly diverse molecular population, with more than a thousand molecular species arising from combinations of different polar head groups and fatty acyl chains attached to a glycerol backbone ().

The biosynthesis and remodeling of phospholipids have been extensively studied in mammals. Classical studies established that phospholipids are produced and matured through two major pathways: the Kennedy pathway () and the Lands cycle (). Central enzymes producing phospholipid diversity in both pathways are lysophospholipid acyltransferases (LPLATs), which catalyze the esterification of lysophospholipids using acyl-CoAs () (Figure 1A). In the Kennedy pathway, phosphatidic acid (PA) is produced by LPLATs from lyso-PA (LPA) and acyl-CoAs. PA is subsequently converted into several other phospholipids, phosphatidylcholine (PC), phosphatidylethanolamine (PE), phosphatidylglycerol (PG), phosphatidylinositol (PI), phosphatidylserine (PS), and cardiolipin (CL). Fatty acids of these phospholipids are reacylated by LPLATs in the Lands cycle (Figure 1B). Fourteen LPLATs have been identified in human and mouse, and recent studies combining enzyme characterization with knockout mouse models have begun to reveal their roles in the regulation of membrane phospholipid compositions. These LPLATs belong to two families: 1-acylglycerol-3-phosphate O-acyltransferase (AGPAT) and membrane-bound O-acyltransferase (MBOAT) families. AGPAT motifs 1–4 and MBOAT motifs A-D, essential for LPLAT activity, have been identified (Figure 1C) (–). LPLAT structures have been predicted by homology modeling () and AlphaFold2 analysis (), and several 3D structures have also been experimentally determined (). Conserved motifs in both AGPAT- and MBOAT-family LPLATs may form pockets for substrate binding.

Figure 1

).

Phospholipid compositions vary widely across organisms, reflecting differences in metabolism, environmental adaptation, and cellular complexity. Although multiple mechanisms can influence phospholipid diversity, LPLAT-mediated remodeling is likely to play a central role in shaping membrane lipid compositions in most organisms. LPLATs have been studied in several model organisms, including Saccharomyces cerevisiae (), Caenorhabditis elegans (), Drosophila melanogaster (), Escherichia coli (, ), and Arabidopsis thaliana (), as well as in Homo sapiens and Mus musculus (, ). Previous studies showed that deficiency or mutation of LPLATs changed the phospholipid compositions and induced biological dysfunctions, including cell growth reduction in S. cerevisiae (), impaired epithelial cell division in C. elegans (), and abnormal germ cell development in D. melanogaster (). In mammals, several phenotypes have been reported in association with LPLAT deficiency or mutation, including visual dysfunction (), male infertility (), delayed neuronal migration in the cortex (), alleviation of neuropathic pain (), and suppression of ferroptosis in the liver () in LPLAT gene-deficient mice; intellectual disability () and lipodystrophy () in humans with LPLAT gene mutations; and suppression of cell growth () in a human cell line. Several reviews have summarized the effects of human and mouse LPLAT mutations and deficiencies on biological functions (, ). The number of LPLATs has increased during evolution. For example, humans possess 14 LPLATs, whereas yeast has only two (). Consistent with this expansion, the fatty acid composition of membrane phospholipids is more diverse in humans than in yeast, and phospholipid composition differs among species. However, systematic analyses of LPLAT distribution and evolution across diverse taxa remain limited. In addition, pathogens may expand their LPLATs to remodel their membranes in order to survive and thrive in hostile host niches, evade the host immune system, or optimize virulence. Yet, little is known about the diversity of LPLATs in pathogenic organisms. In particular, little is known about the diversity of LPLATs in pathogenic organisms.

In this study, we asked whether the diversification of membrane phospholipid composition across organisms is evolutionarily linked to LPLATs, and whether it is associated with functional changes such as pathogenicity. Here, we conducted a comprehensive survey of LPLAT proteins across representative prokaryotic and eukaryotic lineages, including bacteria and major eukaryotic supergroups such as Opisthokonta, Amoebozoa, Stramenopiles–Alveolata–Rhizaria (SAR), Archaeplastida, and Excavata, based on the classification framework proposed by Adl et al. (, ). The phylogenetic analyses reveal the diversification of LPLATs during evolution, with notable expansion in vertebrates. Furthermore, comparative analyses suggest that certain pathogenic species harbor an expanded repertoire of LPLAT homologs compared with their non-pathogenic relatives. These findings indicate that diversification of phospholipid remodeling enzymes may have accompanied organismal evolution and the emergence of pathogenicity. This study provides an evolutionary perspective on membrane phospholipid metabolism and establishes a framework for exploring the functional roles of LPLATs in biological diversity and host–pathogen interactions.

2 Methods

2.1 Identification of candidate LPLAT proteins

Orthologs of the fourteen known human and mouse LPLAT enzymes were searched in representative prokaryotic and eukaryotic organisms (Supplementary Table 1) (, ) using the following strategy.

(i) The amino acid sequences of SLC1 (NP_010231, AGPAT family) and ALE1 (NP_014818, MBOAT family) from S. cerevisiae were used individually as queries in tBLASTn searches against the NCBI RefSeq genomes database and selected organism genomes.

(ii) Additional searches were performed in the NCBI database using combinations of organism names and the keywords “1-acylglycerol-3-phosphate acyltransferase”, “AGPAT”, “membrane-bound O-acyltransferase”, or “MBOAT.” In addition, previously reported LPLATs in several organisms were included.

(iii) Candidate LPLAT amino acid sequences identified in each organism were subsequently used as queries for additional tBLASTn searches against the NCBI non-redundant (nr) database and organism-specific genome databases to identify related genomes.

Candidate proteins were selected if they contained the conserved motifs of AGPAT or MBOAT families, including the AGPAT motifs 1(HxxxxD) and 3(FPEG) (, ) or the MBOAT motif B(WHGxxxGYxxxF) () (Figure 1C). When multiple proteins with identical internal regions were identified, the protein with the longest N- or C-terminal region was selected as the representative sequence. Accession numbers (AccNo) of the proteins and LPLAT names used in this study are listed in Supplementary Tables 1–3. The number of clusters was shown in Supplementary Table 1.

2.2 Sequence alignment and phylogenetic analysis

Multiple sequence alignments were generated using Clustal Omega () with default parameters. Phylogenetic trees were inferred from the multiple sequence alignment of LPLAT homolog amino acid sequences using the maximum likelihood method implemented in PhyML (). The LG model of amino acid substitution and bootstrap analysis with 1,000 replicates were used for the phylogenetic analysis. For phylogenetic tree inference, only the alignments of the AGPAT motif 1–4 or MBOAT motif A–D regions were used. All analyses were performed using SeaView (), and the final outputs were generated within this software. All alignments were shown in Supplementary Figures 1–3.

3 Results

3.1 Identification of candidate LPLAT proteins across representative organisms

In humans and mice, LPLATs have been classified into two families: LPLAT1–10 belonging to the AGPAT family and LPLAT11–14 belonging to the MBOAT family (). Because the AGPAT and MBOAT families are evolutionarily distinct and share low sequence homology, phylogenetic trees were constructed separately for each family. To investigate their evolutionary distribution, orthologous LPLAT candidates were identified from representative prokaryotic and eukaryotic clades using the strategy described in the Methods section. All identified candidates contained conserved motifs of either the AGPAT family (Motif 1: HxxxxD and Motif 3: FPEG) (, ) or the MBOAT family (Motif B: WHGxxxGYxxxF) ().

As shown in Supplementary Tables 1–3, no MBOAT family LPLATs were detected in several species, including the prokaryote E. coli K-12, the choanoflagellate Monosiga brevicollis, the alveolate parasite Plasmodium falciparum, and the excavate Trypanosoma cruzi. In addition, no candidate LPLAT proteins were identified in the archaeal genomes examined.

3.2 Classification of LPLATs in representative organisms

To examine overall patterns of LPLAT diversity across major evolutionary groups, representative species were selected: E. coli K-12 (EckLPLAT), D. melanogaster (DmLPLAT), H. sapiens (HsLPLAT), S. cerevisiae S288C (ScLPLAT), M. brevicollis MX1 (MbLPLAT), Entamoeba dispar SAW760 (EdLPLAT), P. falciparum 3D7 (PfLPLAT), A. thaliana (AtLPLAT), and T. cruzi (TcLPLAT) (Supplementary Tables 1–3). Phylogenetic trees were constructed using SeaView () (Figure 2).

Figure 2

Figure 2 shows phylogenetic trees depicting the evolutionary relationships among AGPAT (Figure 2A) and MBOAT (Figure 2B) LPLAT families across various species. In Figure 2A, the phylogenetic analysis identified 7 and 16 major clades with bootstrap values ≧90 and ≧50, respectively, whereas in Figure 2B, 4 and 7 major clades were identified with bootstrap values ≧90 and ≧50, respectively. Although several deep ancestral nodes exhibited low bootstrap support (<50), likely due to low sequence homology across highly divergent evolutionary clades, the distinct separation of major clusters (such as HsLPLAT1/2 vs. HsLPLAT8/9/10) remains highly supported. Further structural or synteny analyses may be required to firmly resolve these deeper evolutionary relationships. In the AGPAT family, LPLATs were grouped into three major clusters, comprising HsLPLAT1/2, HsLPLAT3/4/5/6/7, and HsLPLAT8/9/10. In contrast, the MBOAT family was divided into two major clusters, HsLPLAT11/12 and HsLPLAT13/14. Although HsLPLAT3 and HsLPLAT4 belong to the same cluster, they are located on chromosomes 21 (https://www.ncbi.nlm.nih.gov/protein/NP_064517) and 6 (https://www.ncbi.nlm.nih.gov/protein/NP_064518), respectively. Similarly, HsLPLAT8, HsLPLAT9, and HsLPLAT10, which form another cluster, are located on chromosomes 5 (https://www.ncbi.nlm.nih.gov/protein/NP_079106), 16 (https://www.ncbi.nlm.nih.gov/protein/NP_060309), and 15 (https://www.ncbi.nlm.nih.gov/protein/NP_705841), respectively. HsLPLAT13 and HsLPLAT14 also belong to the same cluster and exhibit similar substrate preferences, yet they are located on chromosomes 2 (https://www.ncbi.nlm.nih.gov/protein/NP_620154) and 6 (https://www.ncbi.nlm.nih.gov/protein/NP_001073949), respectively. The observations suggest that these gene pairs or groups are unlikely to have arisen through tandem duplication. HsLPLAT7 diverged earlier than the other LPLATs, suggesting that it may represent a distinct evolutionary lineage (Figure 2A). AtLPLAT8 and AtLPLAT9 show a close relationship with a high bootstrap value. Although both genes are located on chromosome 1 (https://www.ncbi.nlm.nih.gov/protein/NP_176493; https://www.ncbi.nlm.nih.gov/protein/NP_172724), they are not adjacent to each other, suggesting that they may not be products of tandem duplication (Figure 2B).

The reported biochemical characteristics of human and mouse LPLATs (mammalian LPLATs) were shown in Supplementary Table 4. HsLPLATs are not a single monophyletic group but instead split into distinct evolutionary lineages, suggesting that different HsLPLAT proteins may have evolved to handle different lipid substrates and serve distinct physiological roles within the cells. Although functional verification is required to determine whether non-human LPLATs possess identical enzymatic activities, several observations suggest potential similarities. For example, in S. cerevisiae, ScLPLAT1 (SLC1) incorporates C18:1 into LPA and lysophosphatidylcholine (LPC) (, ), whereas ScLPLAT2 (LPT1, ALE1, SLC4) can utilize a broad range of lysophospholipids and can incorporate highly unsaturated fatty acids in vitro (, –). In D. melanogaster, DmLPLAT4 (farjavit), DmLPLAT5 (nessy), and DmLPLAT6 (oysgedart) exhibit substrate preferences for LPI/C20:4-CoA, LPC/C20:4-CoA, LPC/C18:1-CoA as substrates, respectively (). HsLPLAT11, HsLPLAT12, and HsLPLAT13/14 belong to the same clusters as DmLPLAT4, DmLPLAT5, and DmLPLAT6, respectively, and show substrate preferences () similar to those of the corresponding DmLPLATs (). They may be orthologs with conserved functions across vertebrates and invertebrates. Each LPLAT in S. cerevisiae and D. melanogaster exhibits characteristics similar to those of the corresponding mammalian LPLATs (Figure 2). Although PfLPLAT1 was not classified into any human LPLAT group, LPA acyltransferase (LPAAT) and LPC acyltransferase (LPCAT) activities with several acyl-CoAs have been reported (). Notably, mammalian LPLAT5 and mammalian LPLAT7, which incorporate C18:0 or C18:1 into phospholipids (), were observed only in metazoans among the species examined. These enzymes may represent a vertebrate-specific expansion. Mammalian LPLAT6 and mammalian LPLAT7 esterify fatty acids at the sn-1 position of phospholipids (Figure 1A) (–). A close relationship between HsLPLAT6 and plant AtLPLAT1, AtLPLAT2, and AtLPLAT3 was observed. Further studies are needed to characterize the biochemical properties of HsLPLAT6 and AtLPLAT1–3.

3.3 Diversification of LPLATs in Metazoa

To further investigate LPLAT evolution within metazoa, phylogenetic trees were constructed using representative invertebrates—Amphimedon queenslandica (AqLPLAT), Nematostella vectensis (NvLPLAT), Magallana gigas (MgLPLAT), D. melanogaster (DmLPLAT), C. elegans (CeLPLAT), and Capitella teleta (CtLPLAT)—and vertebrates—Danio rerio (DrLPLAT), Xenopus tropicalis (XtLPLAT), Anolis carolinensis (AcLPLAT), Gallus gallus (GgLPLAT), and H. sapiens (Supplementary Tables 1–3). Figure 3 shows phylogenetic trees depicting the evolutionary relationships among AGPAT (Figure 3A) and MBOAT (Figure 3B) LPLAT families across various species. In Figure 3A, the phylogenetic analysis identified 11 and 30 major clades with bootstrap values ≧90 and ≧50, respectively, whereas in Figure 3B, 5 and 17 major clades were identified with bootstrap values ≧90 and ≧50, respectively.

Figure 3

Members of the LPLAT2–9 group in the AGPAT family and the LPLAT12–14 group in the MBOAT family were identified in D. rerio, X. tropicalis, A. carolinensis, G. gallus, and H. sapiens. However, the analysis revealed that D. rerio lacked an ortholog of HsLPLAT1 but contained two orthologs of HsLPLAT13 (DrLPLAT13 and DrLPLAT15), suggesting a lineage-specific duplication, because DrLPLAT13 and DrLPLAT15 are located on chromosomes 17 (https://www.ncbi.nlm.nih.gov/protein/XP_073783893.1/) and 20 (https://www.ncbi.nlm.nih.gov/protein/NP_001265531.1/), respectively. In contrast, G. gallus lacked the orthologs of HsLPLAT10 and HsLPLAT11. The full complement of 10 AGPAT-family LPLATs and 4 MBOAT-family LPLATs appears to have been established before the emergence of vertebrates. These findings suggest that the diversification of phospholipid remodeling enzymes had largely occurred prior to vertebrate evolution.

LPLAT1 and LPLAT2 cluster sequences may share a common ancestor and exhibit close evolutionary relationships across vertebrates and some invertebrates. LPLAT5 and LPLAT7 showed high bootstrap values of 96 and 100, respectively, and formed distinct major lineages in the phylogenetic tree. LPLAT8, 9, and 10 formed a distinct group, which may reflect relatively recent duplication events in vertebrates. Based on the substrate specificities of mouse LPLATs, these enzymes are predicted to be involved in sn-1/sn-2 remodeling; however, no clear commonality in fatty acid recognition has been identified (). Further analyses of the substrate specificities of LPLATs from additional species will be necessary.

Vertebrate LPLATs formed tight clusters, suggesting strong conservation of these specific LPLAT isoforms across vertebrates. In particular, vertebrate LPLAT11 and LPLAT12 showed high bootstrap support values of 99 and 82, respectively, suggesting that they are orthologs. LPLAT14 and LPLAT13 groups showed strong conservation among vertebrates. LPLAT14 sequences from different vertebrate species formed a tight cluster, whereas LPLAT13 sequences formed a sister group to LPLAT14, suggesting a relatively recent duplication event in vertebrates. Although DrLPLAT13 and DrLPLAT15 are located on different chromosomes, they may have arisen through a gene duplication event rather than tandem duplication. In the LPLAT12 group, the relatively long branch lengths of the invertebrate LPLATs may reflect a greater accumulation of amino acid substitutions. The distinct clustering of vertebrate and invertebrate sequences may indicate that these lineages diverged relatively early during evolution.

Overall, vertebrates displayed a greater diversity of LPLATs than invertebrates. The diversification of LPLATs may have accompanied evolutionary processes. Nevertheless, some invertebrates showed relatively high numbers of LPLAT proteins. For example, C. elegans, M. gigas, and D. melanogaster possessed comparatively large LPLAT repertoires. In particular, C. elegans contained 12 LPLAT proteins. Previous studies reported that triple knockout of CeLPLAT3 (acl-8), CeLPLAT4 (acl-9), and CeLPLAT5 (acl-10) in C. elegans reduces PI species containing C18:0 (), suggesting that these LPLATs contribute to the regulation of PI and phosphoinositide compositions.

3.4 Differences in LPLAT repertoires between pathogenic and non-pathogenic species

To examine whether LPLAT compositions are associated with pathogenicity, LPLAT repertoires were compared between pathogenic and non-pathogenic species. Specifically, this study analyzed E. coli K-12 and the pathogenic strain E. coli O157:H7 (strain G5101) (EcoLPLAT), as well as the non-pathogenic amoeba E. dispar SAW760 and the pathogenic species Entamoeba histolytica HM-1:IMSS (EhLPLAT) (Supplementary Tables 1–3). Phylogenetic analyses were performed together with HsLPLATs (Figure 4). Figure 4 shows phylogenetic trees depicting the evolutionary relationships among AGPAT (Figure 4A) and MBOAT (Figure 4B) LPLAT families across various species. In Figure 4A, the phylogenetic analysis identified 5 and 11 major clades with bootstrap values ≧90 and ≧50, respectively, whereas in Figure 4B, 5 and 5 major clades were identified with bootstrap values ≧90 and ≧50, respectively.

Figure 4

E. coli K-12 possessed two AGPAT-family LPLATs, whereas the pathogenic species E. coli O157:H7 contained three, including an additional protein (EcoLPLAT3) (Figure 4A). In Entamoeba, both Entamoeba species possessed two AGPAT-family paralogs. Unlike E. coli, however, Entamoeba species also possess MBOAT-family LPLATs, resulting in a total of four LPLAT paralogs. The very short branch lengths between EcoLPLAT and EckLPLAT suggest high conservation across these strains. EcoLPLAT1/EckLPLAT1, EcoLPLAT2/EckLPLAT2, EdLPLAT1/EhLPLAT1, EdLPLAT3/EhLPLAT3, EdLPLAT4/EhLPLAT4, EdLPLAT5/EhLPLAT5, and EdLPLAT6/EhLPLAT6 showed high bootstrap values of 88, 100, 91, 99, 100, 100, and 100, respectively. The LPLATs in E. coli and amoebae may suggest lineage-specific gene duplication events that occurred before divergence within each lineage. HsLPLAT3–7 formed a separate evolutionary branch from the bacterial and protozoan groups. Notably, the pathogenic species E. histolytica possesses an additional protein, EhLPLAT7. Within the MBOAT family tree, EhLPLAT7 forms an independent branch that clusters more closely with the Entamoeba LPLAT3/4/5/6 group than with HsLPLATs (Figure 4B). These observations suggest that certain pathogenic species may possess expanded LPLAT repertoires relative to their non-pathogenic relatives.

4 Discussion

Membrane phospholipid biosynthesis has diversified significantly during evolution. To map these relationships, we constructed phylogenetic trees based on the conserved motifs of the AGPAT and MBOAT families. Our findings suggest that the structural diversification of these LPLATs may have driven the expansion of ecological niches and the development of specialized cellular functions. This evolutionary diversity is reflected in modern mammals; for instance, the 14 human and mouse LPLAT isoforms exhibit highly distinct tissue expression patterns and physiological roles. Specifically, LPLAT3 is highly expressed in the retina and testis, and its deficiency leads to visual dysfunction () and male infertility (), respectively. Ultimately, variations in LPLAT expression, substrate preferences, and enzymatic regulation act as primary determinants of cellular lipid profiles, offering new insights into both baseline physiology and host defense mechanisms.

In this study, MBOAT-family LPLATs were not detected in prokaryotes. Although AGPAT-family LPLATs were identified in the prokaryotic bacterium E. coli, no LPLAT candidates were detected in archaeal genomes. Archaeal membrane lipids are structurally distinct from those of bacteria and eukaryotes, as they typically contain ether-linked isoprenoid chains at both the sn-1 and sn-2 positions of the glycerol backbone (). Therefore, the acyl-CoA–dependent esterification catalyzed by LPLATs may not be required in archaeal lipid biosynthesis. Similar to E. coli, MBOAT-family LPLATs were not found in M. brevicollis, P. falciparum, or T. cruzi. MBOAT-family LPLATs are characterized by their ability to incorporate arachidonic acid into PC and PI. These lipids are involved in lipid mediator production and phosphoinositides (PIPs) signaling. Such features may have evolved to enable greater diversity in membrane lipid compositions. During evolution, MBOAT-family LPLATs and metazoan-specific AGPAT-family LPLATs such as mammalian LPLAT5, mammalian LPLAT6, and mammalian LPLAT7 appear to have arisen, potentially reflecting the increasing complexity of membrane lipid remodeling. These enzymes may be particularly important for processes such as polyunsaturated fatty acids (PUFA)s incorporation and remodeling of sn-1 fatty acids.

PA serves as the central precursor for the synthesis of many other phospholipids in the de novo pathway (Kennedy pathway). PA synthesis using a broad range of fatty acids may therefore represent a conserved and essential feature of membrane phospholipid metabolism across diverse organisms. In addition, the phylogenetic analysis identified several LPLAT candidates that could not be clearly classified into mammalian LPLAT groups (Figure 2). As with most LPLATs, enzymatic in vitro analyses combined with lipidomic studies of LPLAT knockout organisms will be essential to clarify the phospholipid species regulated by these LPLATs.

This study suggests possible functional links between LPLAT diversification and specific phospholipid metabolism. HsLPLAT11, which incorporates arachidonic acid into LPI, contributes to PI metabolism and phosphoinositide signaling pathways. A similar enzyme, DmLPLAT4 (farjavit), was identified in the arthropod D. melanogaster. DmLPLAT4 (farjavit) has been reported to synthesize arachidonic-acid–containing PI (), suggesting functional similarity to HsLPLAT11. In addition, the mouse ortholog of HsLPLAT6 catalyzes the esterification of stearic acid at the sn-1 position of PI (). In this study, orthologs of HsLPLAT6 were restricted to A. thaliana (AtLPLATs) and not widely detected among other examined species. This activity may therefore represent a relatively recent evolutionary innovation that contributes to phosphoinositide regulation.

Other LPLAT groups may also have contributed to the evolution of lipid mediator pathways. HsLPLAT12 and DmLPLAT5 (nessy) incorporate arachidonic acid into LPC (, , ). Arachidonic acid is not only a structural component of membrane phospholipids but also the precursor of lipid mediators, eicosanoids (), indicating that these enzymes may have contributed to the evolution of lipid mediator production. However, because the physiological roles of lipid mediators differ widely among organisms, and their biosynthetic mechanisms cannot be directly compared across species, further studies are needed to evaluate the relationship between LPLAT evolution and lipid mediators.

Among phospholipids, platelet-activating factor (PAF) is a potent phospholipid mediator (, ). PAF is rapidly synthesized in response to extracellular stimuli (, ), and the PAF biosynthetic enzyme lyso-PAF acetyltransferase () is regulated through three distinct mechanisms (). Mammalian LPLAT9 catalyzes the lyso-PAF acetyltransferase reaction responsible for PAF biosynthesis (, ). Notably, mammalian LPLAT9 is the only member among the fourteen LPLATs whose phosphorylation and expression are induced by innate immune stimuli, such as toll-like receptor (TLR)4 or TLR9 ligands (, , ). These observations suggest that diversification of the LPLAT family may have contributed not only to membrane phospholipid remodeling but also to the evolution of lipid-mediated signaling systems involved in host defense. This group in the phylogenetic tree also includes mammalian LPLAT8 and mammalian LPLAT10, which also have distinctive characters. Mammalian LPLAT8 contributes to the synthesis of pulmonary surfactant phospholipids essential for respiration and increases the saturated fatty acids in PC (–). Mammalian LPLAT10 incorporates very long chain fatty acid into the sn-1 position of PC (, ). Although similar physiological functions cannot be assumed in non-mammal organisms, these enzymes highlight the potential importance of LPLAT-mediated phospholipid remodeling in evolutionary adaptations.

Docosahexaenoic acid (DHA, C22:6) is a very long-chain fatty acid containing six double bonds and plays multiple functional roles in mammals (, , , ). Several studies have reported LPLAT activities that generate DHA-containing phospholipids, mediated by mammalian LPLAT3 (), mammalian LPLAT4 (), and mammalian LPLAT10 (, , –). In addition to serving as membrane components, DHA-containing phospholipids function as a reservoir for lipid mediators (), similar to arachidonic acid (C20:4)-containing phospholipids (). Considering the number of double bonds in fatty acids may provide important insights into the evolution of phospholipid biosynthesis.

Other enzymes such as mammalian LPLAT13 and mammalian LPLAT14, which preferentially utilize oleoyl-CoA (, ), and DmLPLAT6 (oysgedart) () and ScLPLAT2 (LPT1, ALE1, SLC4) (, –) in this cluster also possess similar preferences. Although ScLPLAT2 has been reported to recognize PUFA substrates in vitro, yeast membranes contain relatively low levels of PUFAs (). Similar discrepancies have been reported in mouse LPLAT studies (), in which in vitro enzymatic activity does not always correlate with the phospholipid species altered in vivo following LPLAT deficiency or overexpression.

Membrane phospholipids play important roles in immune function. Membrane nanodomains (also called microdomains), which are enriched in specific phospholipids, are involved in various cellular functions, including immunoreceptor signaling. These nanodomains serve as platforms for the accumulation of signaling molecules involved in immune receptor pathways, such as Src-family kinases, G proteins, and PIPs, thereby influencing T cell receptor (TCR) and B cell receptor (BCR) signaling (). PI is a precursor of PIPs. In addition, fatty acids contained in phospholipids, such as arachidonic acid, DHA, and linoleic acid, serve as precursors for lipid mediators. These fatty acid-derived lipid mediators are produced in response to extracellular stimuli, originating from membrane phospholipids in immune cells. Thus, membrane phospholipids and LPLATs that regulate their remodeling play critical roles in the modulation of immune signaling. Furthermore, understanding pathogen-derived phospholipids and LPLATs will provide important insights into pathogen control (). Elucidating LPLATs in both pathogens and hosts will also contribute to advances in immunological research.

Interestingly, differences in LPLAT repertoires were also observed between pathogenic and non-pathogenic species within the same genus. The pathogenic bacterium E. coli O157 contains one additional AGPAT-family LPLAT, EcoLPLAT3, compared with the non-pathogenic strain E. coli K-12. Similarly, the pathogenic amoeba E. histolytica contains one additional MBOAT-family LPLAT, EhLPLAT7, compared with the non-pathogenic species E. dispar. These additional enzymes may contribute to the synthesis of phospholipids involved in pathogenicity or may compensate for other LPLATs. The phospholipid biosynthetic pathway in E. histolytica has also been proposed, highlighting the importance of future studies (). Functional studies, including gene disruption analysis, will be required to further clarify their roles in virulence. Comparing only two pairs of organisms (E. coli and Entamoeba), which are isolated data points, is insufficient to suggest a broader evolutionary or biological trend linking lipid remodeling to virulence. Expanding such analyses to additional pathogens may provide important insights into host–pathogen interactions.

To inhabit diverse environments rather than a specific one, organisms likely require a broad range of membrane properties to ensure adaptability. Furthermore, it can be inferred that diverse membrane phospholipid compositions were advantageous for the formation of complex tissues from multiple cell types and for the acquisition of specialized functions. Evaluating the effects of phospholipids on physicochemical properties is also an important issue. For example, LPLAT3-deficient mice, which show reduced levels of DHA-containing phospholipids, exhibit disruption of photoreceptor disc structures and abnormalities in sperm morphology. However, the molecular mechanisms by which phospholipid properties regulate these processes remain unclear. The observation that LPLATs are most abundant in animals, particularly in mammals, suggests a link between organismal evolution and the diversification of phospholipids and LPLATs. It can be inferred that diverse phospholipid compositions are advantageous for generating the wide range of cellular and organelle morphologies. In addition, the presence of highly unsaturated fatty acids may contribute not only to the diversification of lipid mediators but also to a potential role as a sink for lipid oxidation. Elucidating these mechanisms, including their possible involvement in the acquisition of pathogenicity, remains an important direction for future research.

5 Limitations

This study is based primarily on publicly available genome annotations, and therefore, some species may lack the information. BLAST-based searches can be influenced by the selected query sequence and the search parameters.

Phospholipid compositions in cells are regulated not only by LPLATs but also by multiple additional factors, such as acyl-CoA synthetases (), phospholipase As (, ), transacylases (), substrate availability including lipid nutrients, and intracellular localization of LPLATs. These regulatory mechanisms may vary among organisms. For example, the bacterial enzyme Aas in E. coli can synthesize PA using acyl-acyl carrier protein (acyl-ACP) rather than acyl-CoA (). Although Aas was not included in the current analysis because of its distinct enzymatic mechanism, such enzymes indicate that additional pathways may contribute to membrane phospholipid remodeling. Therefore, LPLATs alone cannot fully explain the diversity of membrane phospholipid compositions.

In this study, phylogenetic analyses were performed using newly identified candidate LPLAT proteins. Further biochemical and biological characterization of LPLATs from diverse organisms beyond humans and mice, both in vitro and in vivo, will contribute to a deeper understanding of the evolution of LPLATs and membrane phospholipids. Additionally, the phospholipid compositions of individual organisms should be analyzed to better understand the evolution of phospholipid biosynthesis. While sequence conservation represents potential enzymatic capabilities, the known variations between in vitro activity and in vivo phospholipid profiles highlight that sequence alone does not dictate cellular phenotype. Thus, these phylogenetic clusters should be viewed as predictive bases needing in vivo lipidomic studies.

Our results contained the weak nodes. While the deeper backbone nodes of the AGPAT tree lack strong bootstrap support, the functional divergence of these lineages is independently supported by distinct chromosomal distributions (e.g., chromosomes 5, 16, and 15 for HsLPLAT8/9/10) and diverging substrate preferences observed in experimental models.

6 Conclusion

Comparative analyses of LPLATs across diverse organisms suggest that membrane phospholipid remodeling diversified during organismal evolution in parallel with LPLAT diversification. Expansion of the LPLAT family, along with regulation of their tissue distribution and enzymatic activities, may have contributed to the generation of diverse membrane environments in living organisms. These findings support the idea that LPLAT diversification is closely associated with organismal evolution. The identification of pathogen-specific LPLAT expansions further indicates that these enzymes may contribute to the emergence of pathogenic traits, although further biochemical characterization of these LPLATs will be necessary. Even if certain LPLATs are not directly involved in pathogenicity, inhibitors targeting pathogen-specific LPLATs may represent potential targets for antimicrobial drug development. Thus, the evolutionary study of phospholipid metabolism provides insights not only into membrane biology and evolutionary processes but also into host defense and the control of infectious diseases.

Statements

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: Amino acid sequences from NCBI database. Only accession numbers are listed in the Supplementary Tables 2 and 3.

Author contributions

HS: Visualization, Project administration, Funding acquisition, Writing – review & editing, Formal analysis, Data curation, Supervision, Writing – original draft, Conceptualization, Methodology, Investigation. YK: Data curation, Formal analysis, Writing – review & editing.

Funding

The author(s) declared financial support was received for this work and/or its publication. This work was supported by the Japan Institute for Health Security (JIHS) Intramural Research Fund (22T001, 24A1013, and 24A2011) and ONO Medical Research Foundation.

Acknowledgments

We would like to thank T. Shimizu, T. Hashidate-Yoshida, and T. Suzuki (Japan Institute for Health Security), and K. Yanagida and K. Nagata (The Jikei University School of Medicine) for their constructive comments.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The author HS declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Generative AI statement

The author(s) declared that generative AI was used in the creation of this manuscript. Some parts of this manuscript were improved using a generative AI tool (ChatGPT), which was used solely for language editing and did not influence the scientific content.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2026.1838657/full#supplementary-material

Supplementary Figure 1

(A) Sequence alignment utilized in Figure 2A. AGPAT motifs in HsLPLAT8 were shown. (B) Sequence alignment utilized in Figure 2B. MBOAT motifs in HsLPLAT12 were shown.

Supplementary Figure 2

(A) Sequence alignment utilized in Figure 4A. AGPAT motifs in HsLPLAT8 were shown. (B) Sequence alignment utilized in Figure 4B. MBOAT motifs in HsLPLAT12 were shown.

Supplementary Figure 3

(A) Sequence alignment utilized in Figure 3A. AGPAT motifs in HsLPLAT8 were shown. (B) Sequence alignment utilized in Figure 3B. MBOAT motifs in HsLPLAT12 were shown.

Abbreviations

AccNo, accession numbers; acyl-ACP, acyl-acyl carrier protein; AGPAT, 1-acyl-glycerol-3-phosphate O-acyltransferase; CL, cardiolipin; DHA, docosahexaenoic acid; LPA, lyso-PA; LPAAT, LPA acyltransferase; LPC, lysophosphatidylcholine; LPCAT, LPC acyltransferase; LPI, lysophosphatidylinositol; LPLAT, lysophospholipid acyltransferases; MBOAT, membrane-bound O-acyltransferase; PA, phosphatidic acid; PAF, platelet-activating factor; PC, phosphatidylcholine; PE, phosphatidylethanolamine; PG, phosphatidylglycerol; PI, phosphatidylinositol; PIP, phosphoinositide; PS, phosphatidylserine; PUFA, polyunsaturated fatty acid; SAR, Stramenopiles–Alveolata–Rhizaria; and TLR, toll-like receptor. LPLAT names were described in Supplementary Table 1.

References

Summary

Keywords

cellular membrane, LPLAT, membrane diversity, membrane remodeling, pathogen, phospholipid

Citation

Shindou H and Kawai Y (2026) An evolutionary analysis of lysophospholipid acyltransferases generating membrane phospholipid diversity. Front. Immunol. 17:1838657. doi: 10.3389/fimmu.2026.1838657

Received

25 March 2026

Revised

01 July 2026

Accepted

22 July 2026

Published

10 August 2026

Volume

17 - 2026

Edited by

Pierre Pontarotti, Centre National de la Recherche Scientifique (CNRS), France

Reviewed by

Dennis B. Bela-ong, Gyeongsang National University, Republic of Korea

Cristina Colomer Winter, University of Geneva, Switzerland

Updates

Copyright

*Correspondence: Hideo Shindou,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics