Next Article in Journal
Urticaria: A Narrative Overview of Differential Diagnosis
Next Article in Special Issue
A Pig Model to Assess Skin Lesions after Apomorphine Application
Previous Article in Journal
In Vitro Screening Studies on Eight Commercial Essential Oils-Derived Compounds to Identify Promising Natural Agents for the Prevention of Osteoporosis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Oleic Acid-Containing Phosphatidylinositol Is a Blood Biomarker Candidate for SPG28

1
Division of Genomics, Medical Institute of Bioregulation, Kyushu University, 3-1-1, Maidashi, Higashi-ku, Fukuoka 812-8582, Japan
2
Division of Metabolomics, Medical Institute of Bioregulation, Kyushu University, 3-1-1, Maidashi, Higashi-ku, Fukuoka 812-8582, Japan
3
Department of Nutritional Sciences, Nakamura Gakuen University, 5-7-1, Befu, Jonan-ku, Fukuoka 814-0198, Japan
4
Department of Neurosurgery, Kurume University School of Medicine, 67 Asahi-Machi, Kurume, Fukuoka 830-0011, Japan
5
Department of Food and Nutrition, Beppu University Junior College, 82, Kitaishigaki, Oita 874-8501, Japan
6
Department of Neurology and Geriatric Medicine, Ehime University Graduate School of Medicine, 454, Shitsukawa, Toon 791-0295, Japan
*
Author to whom correspondence should be addressed.
Biomedicines 2023, 11(4), 1092; https://doi.org/10.3390/biomedicines11041092
Submission received: 7 March 2023 / Revised: 28 March 2023 / Accepted: 1 April 2023 / Published: 4 April 2023

Abstract

:
Hereditary spastic paraplegia is a genetic neurological disorder characterized by spasticity of the lower limbs, and spastic paraplegia type 28 is one of its subtypes. Spastic paraplegia type 28 is a hereditary neurogenerative disorder with an autosomal recessive inheritance caused by loss of function of DDHD1. DDHD1 encodes phospholipase A1, which catalyzes phospholipids to lysophospholipids such as phosphatidic acids and phosphatidylinositols to lysophosphatidic acids and lysophoshatidylinositols. Quantitative changes in these phospholipids can be key to the pathogenesis of SPG28, even at subclinical levels. By lipidome analysis using plasma from mice, we globally examined phospholipids to identify molecules showing significant quantitative changes in Ddhd1 knockout mice. We then examined reproducibility of the quantitative changes in human sera including SPG28 patients. We identified nine kinds of phosphatidylinositols that show significant increases in Ddhd1 knockout mice. Of these, four kinds of phosphatidylinositols replicated the highest level in the SPG28 patient serum. All four kinds of phosphatidylinositols contained oleic acid. This observation suggests that the amount of oleic acid-containing PI was affected by loss of function of DDHD1. Our results also propose the possibility of using oleic acid-containing PI as a blood biomarker for SPG28.

1. Introduction

Hereditary spastic paraplegia (HSP) is a genetic neurological disorder with progressive spasticity of the lower limbs and muscle weakness caused by axonal damage in the pyramidal tract [1,2,3,4,5]. The prevalence of HSP is estimated to be 1.2 to 9.6 per 100,000 [6]. HSP is genetically heterogeneous and is classified into 83 types, at least by 2022, with different responsible genes [7]. HSP is thought to be caused by the obstacles of mitochondrial metabolism, lipid metabolism, membrane transport, axonal transport and myelin formation, and its pathogenesis is extremely diverse [4,6]. Although diagnoses of specific types of HSPs are generally difficult due to the subtle or no differences in their symptoms among types, the molecular mechanisms of pathogenesis can be at least partially different depending on the causative genes. Due to the late onset of HSP [8], it would be useful if we could identify their types with biomarkers before HSP symptoms develop. In this study, we quantitated lipids to identify molecule species showing quantitative changes associated with DDHD1 dysfunction as a potential biomarker for the diagnosis of HSPs.
Spastic paraplegia type 28 (SPG28;OMIM#609340) is an autosomal recessive type of HSP caused by DDHD1 dysfunction [9,10]. DDHD1 encodes phospholipase A1 (PLA1), which catalyzes phosphatidic acid (PA) and phosphatidylinositol (PI) to lysophosphatidic acid (LPA) and lysophosphatidylinositol (LPI), respectively. PLA1 specifically hydrolyzes the sn-1 ester bonds of PAs and PIs, producing lysophospholipid with a hydrocarbon group at a sn-2 position, while PLA2, another type of phospholipase A, specifically hydrolyzes the sn-2 ester bonds of PAs and PIs [11,12]. PLA1 can be classified into two types: those that function extracellularly (extracellular PLA1) and those that function intracellularly (intracellular PLA1, iPLA1). DDHD1 is included in the latter type. In mammals, there are three families of iPLA1: DDHD1, p125 and DDHD2 [13,14,15,16]. In humans, DDHD1 is known to be highly expressed in the central nervous system and testes [10,17]. p125 and DDHD2 are expressed in a wide range of tissues outside the nervous system [13,15,16,18]. These three iPLA1s are known to have distinct subcellular localizations [13,16,19]. Since iPLA1 is highly conserved among vertebrates, animal models are useful for observing changes in substrates due to PLA1 deficiency [13].
SPG28 is a very rare genetic disorder that has been reported only from six pedigrees thus far [10,20,21,22]. We have previously reported a novel homozygous 4 bp deletion of the DDHD1 gene as the responsible variant for SPG28 [22]. We have established Ddhd1 knockout mice carrying a 5 bp deletion at a very similar position as in the patients and have replicated neurological and histochemical symptoms as seen in the SPG28 patients [23]. Apart from our report, there are two studies that have generated and analyzed Ddhd1 knockout mice (Figure 1 and Table 1) [17,24]. Both strains lack function of the Ddhd1 gene by inserting a lacZ trapping cassette and a neo cassette into the intron of Ddhd1, inducing frameshift [25]. One study reported abnormal sperm motility and mitochondrial morphology [17], while another study reported a significant increase in PI 18:1/20:4 and a significant decrease in LPI 20:4 in the cerebrum of Ddhd1 knockout mice [24]. Our Ddhd1 knockout mice also consistently showed significant changes in certain phospholipids in the brain tissues [23,24]. However, quantitative changes of phospholipids have not been examined in other tissues, such as peripheral blood.
HSP is clinically as well as genetically heterogeneous. Neurological examination of patients often fails to distinguish their types without DNA testing [1]. Whole exome sequencing (WES) is commonly used as DNA testing to conduct molecular diagnoses of specific types of HSP [3]. Therefore, if biomarkers specific to particular types of HSP and applicable to patients’ peripheral blood are available, they lead to cheaper and faster molecular diagnosis without DNA testing. Furthermore, biomarkers are potentially useful for examining pathogenic mechanisms, disease progression and drug efficacy.
Biomarkers have been established for some monogenetic diseases, such as phenylketonuria (PKU) and alkaptonuria (AKU). These diseases are caused by defects in genes encoding enzymes PHA and HGD as in the case of DDHD1 in SPG28 [26,27]. L-phenylalanine and homogentistic acid are substrates of PHA and HGA, respectively, and have been established as biomarkers for these diseases [28,29]. In monogenic diseases caused by defects in enzyme genes, it is most straightforward to examine the substrates of the enzymes as candidate biomarkers. As quantitative changes of specific phospholipids in Ddhd1 knockout mice have been reported [23,24], lipid molecules are good candidates for the biomarker of SPG28. Since there are many species of phospholipids with different hydrocarbon groups, it is essential to distinguish molecular species with different side chains by high-resolution lipidome analyses. To address the need for SPG28 biomarkers, here we searched for lipid molecules whose abundances are altered by DDHD1 dysfunction.

2. Materials and Methods

2.1. Animals

We previously established Ddhd1 knockout mouse strain (Ddhd1 (−/−)) by introducing a 5-base-pair deletion in the second exon of the Ddhd1 gene [23]. The deletion results in a premature termination within the second exon of the DDHD1 gene, resulting in the expression of a protein completely lacking the DDHD domain. They were fed a standard pellet diet (CLEA Japan, Inc., Tokyo, Japan) and filtered water. The animals were kept under condition of a 12:12 h light:dark cycle and housed in groups of two to five animals per cage. The cages were changed once per week. Mouse husbandry and all mouse experiments were carried out at The Animal Center at Nakamura Gakuen University. All mice were bred in SPF (specific pathogen free) area. Although we did not examine sperm motility abnormalities in our Ddhd1 (−/−), we maintained the strain by mating Ddhd1 (+/−) females and Ddhd1 (+/−) males because they were considered to be infertile based on a previous report [17]. The six individuals in each genotype were chosen completely at random and used for blood collection. Animals used in the experiments were observed by visual inspection to be in normal health status. The background of all mice used in the experiments was C57BL/6J. The animal experiments were conducted with the approval of the Animal Ethics Committee of Nakamura Gakuen University (protocol code: 2016-1; date of approval: 4 April 2016). The study is reported in accordance with the recommendation of the Animal Research: Reporting of In Vivo Experiments (ARRIVE) guidelines.

2.2. Genotyping

Tails were cut with lengths of 2 mm at 2 weeks of age for genotyping. DNA was extracted by heating the tail tissues in 50 mM sodium hydroxide solution at 95 °C for 10 min. After heating, 1M Tris-HCl (pH 8.0) was added to neutralize, and the supernatant obtained by centrifugation at 16,000× g for 10 min was used as a template for PCR. The second exon of Ddhd1 was amplified using primers 5′-AAGGTACATCTGGCTTGAAG-3′ and 5′-GAGCTGTGGGTATAGTTGTG-3′. Genotypes were determined by Sanger sequencing of PCR products.

2.3. Sample Preparation of Mouse Plasma

Blood samples were collected from the heart of 6-month-old Ddhd1 (+/+), Ddhd1 (+/−) and Ddhd1 (−/−) (n = 6 each). Blood collection was performed at 6 months of age in all mice. Both male and female mice were used in a 1:1 ratio to experimental groups. The mice were anesthetized by injecting 10 mL per kg of body weight of mixed anesthesia. Mixed anesthesia was composed by adjusting medetomidine hydrochloride (0.3 mg/kg) (Kyoritsu Seiyaku, Tokyo, Japan), midazolam (4 mg/kg) (Teva Takeda Pharma Ltd., Nagoya, Japan) and butorphanol tartrate (5 mg/kg) (Meiji Seika Pharma Co, Ltd., Tokyo, Japan) with saline (Otsuka Pharmaceutical Co., Ltd., Tokyo, Japan). After complete disappearance of the reflection, the chest was opened, and approximately 1 mL of blood was collected from the left ventricle. The syringe used for blood collection was moistened beforehand with sodium heparin (Mochida Pharmaceutical Co., Ltd., Tokyo, Japan). The blood was collected in a tube containing 1 µL of sodium heparin and centrifuged at 3000× g for 20 min at 4 °C. The supernatant was transferred to a new tube and stored as plasma at −80 °C until use.

2.4. Human Sample Preparation

Peripheral blood was collected from the SPG28 patient carrying the homozygous 4 bp deletion in the DDHD1 gene, originally described in Miura et al., 2016 [22], and from unrelated five control individuals free of HSP. The SPG28 patient was a male aged 61 at the time of blood collection. All control individuals were males aged 30, 47, 60, 75 and 76 at the time of the blood collection. Serum was isolated by centrifugation at 3000× g for 10 min after blood collection. Peripheral blood was taken following informed consent. The study was approved by the ethics committees of Kyushu University (protocol code: 2020-444; data of approval: 9 October 2020), Faculty of Medicine and Kurume University School of Medicine (protocol code: 21057; the data of approval: 16 June 2021). All methods were carried out in accordance with relevant guidelines and regulations, including the Declaration of Helsinki.

2.5. Lipid Extraction

We performed lipid extraction from mouse plasma and human sera using an acidic methanol extraction, as described previously, with minor modifications [23]. Briefly, each 50 μL sample of mouse plasma and human sera was mixed with 1 mL of 20 mM acetic acid in methanol containing the internal standards: PC 15:0/18:1 (d7), 2126 pmol; PA 15:0/18:1 (d7), 102 pmol; PI 15:0/18:1 (d7), 118 pmol; LPC 18:1 (d7), 473 pmol; LPI 17:1, 1714 pmol; and LPA 17:0, 571 pmol. The samples were vortexed for 1 min, and then, sonication was performed for 5 min at room temperature. The samples were centrifuged at 16,000× g at 4 °C for 5 min, and 700 μL of the supernatant was transferred in clean tubes. The supernatants were evaporated using a nitrogen gas, and the dried extracts were stored at −80 °C until lipidome analysis.

2.6. Lipidome Analysis for Mouse Plasma and Human Sera

The levels of PC, PI, LPC, LPI and LPA in mouse plasma and human sera were quantified using a supercritical fluid supercritical fluid chromatography (SFC) system (Shimadzu Co., Kyoto, Japan) coupled to a triple–quadrupole mass spectrometer (QqQMS) equipped with an electrospray ionization ion source (Shimadzu Co., Kyoto, Japan) as described previously [30]. The SFC system was equipped with a binary pump (i.e., CO2 pump and a pump for the modifier and make-up solvent), autosampler, temperature-controlled column oven, and back pressure regulator (BPR). The SFC/QqQMS systems and data acquisition were controlled using Shimadzu LabSolutions ver. 5.99 SP2 software.
All samples were resuspended in 200 μL methanol/chloroform (1:1, v/v), and 2 μL was injected onto an ACQUITY UPC2™ Torus diethylamine (DEA) column (3.0 mm i.d. × 100 mm, 1.7 μm particle size, Waters Co., Milford, MA, USA). The temperature of the column oven and autosampler was set at 50 and 4 °C, respectively. BPR was set at 10 MPa. The mobile phase was composed of supercritical carbon dioxide (A) and methanol/water (95/5, v/v) with 0.1% (w/v) ammonium acetate (B). The mobile phase B was also used for modifier and make-up solvent. The flow rate of the mobile phase and make-up pump was set at 1.0 and 0.1 mL/min, respectively. The chromatographic separation was performed using a gradient of increasing the mobile phase B concentration as follows: 0–1 min: gradient was held at 1% B; 1–24 min; linear gradient of 1% to 75% B; 24–26 min; gradient was held at 75% B; 26–26.1 min; gradient was returned to 1% B; 26.1–30 min; the initial conditions were restored and the column was allowed to equilibrate for 3.9 min.
SFC/QqQMS analyses under multiple-reaction monitoring (MRM) were carried out in the positive-ion mode and in the negative-ion mode (polarity switching mode). The QqQMS parameters were set as follows: electrospray voltage of 4.0 kV in the positive-ion mode and −3.5 kV in the negative-ion mode; heating gas flow rate, 10 L/min; drying gas flow rate, 10 L/min; nebulizing gas flow rate, 3 L/min; heat block temperature, 400 °C; desolvation temperature, 250 °C; detector voltage, 2.3 kV; dwell time, 2 ms; pause time, 2 ms, and polarity switching time of 15 ms. The information of optimized MRM parameters for the lipid molecules is shown in Table 2, Table 3, Table 4 and Table 5 and Tables S1–S10. Identification and quantification of the lipid molecules were performed using Multi-ChromatoAnalysT (Beforce Co., Fukuoka, Japan). The quantitative content of identified lipid molecules was calculated using the ratio of the peak area of each analyte to that of the internal standard of its representative lipid class.
When the objective to examine which molecules among several molecular species changed significantly between Ddhd1 (+/+) and Ddhd1 (−/−), the analysis was performed by Student’s t test followed by correction using Benjamini–Hochberg method.

3. Results

3.1. PC, PA and PI in Mouse Plasma

Phosphatidylcholine (PC), PA and PI in mouse plasma from Ddhd1 (+/+) and Ddhd1 (−/−) were analyzed by lipidome analysis. PI and PA are substrates of DDHD1. PC is a precursor of PA, of which the conversion is catalyzed by phospholipase D [31,32]. PCs, PAs and PIs are expected to accumulate due to the loss of DDHD1 enzymatic activity in Ddhd1 (−/−) mice. A total of 88 different species of PCs and 64 different species of PIs were identified in mouse plasma (Tables S1 and S2). PA was not detected in any of the molecular species. Although there was no significant change in the total amount of PC in Ddhd1 (−/−) (Figure 2A), a significant increase was observed in the total amount of PIs in Ddhd1 (−/−), as we expected (p = 7.2 × 10−3) (Figure 2B). We also found significant changes in nine species of PIs, 18:1/20:4, 20:1/20:4, 18:1/20:3, 15:0/24:0, 16:0/18:1, 18:1/18:1, 19:0/20:4, 16:0/20:1 and 18:1/22:5 in Ddhd1 (−/−) (Table 2 and Table S3). These quantitative changes were consistently in the same direction of increase in Ddhd1 (−/−).

3.2. LPC, LPA and LPI in Mouse Plasma

We also quantitated PLA1 metabolites, lysophosphatidylcholine (LPC), LPA and LPI in mouse plasma. We detected 72, 8 and 18 species of LPCs, LPAs and LPIs, respectively (Tables S4–S6). Although we observed no significant changes in LPC and LPI in Ddhd1 (−/−), total LPA was significantly elevated in Ddhd1 (−/−) (p = 1.1 × 10−3) (Figure 3). Among the eight species of LPAs detected in mouse plasma, only LPA 24:0 was significantly elevated in Ddhd1 (−/−) (p = 0.0004) (Table 3). Although we have previously reported a significant decrease in LPI 20:4 (sn-2) in Ddhd1 (−/−) cerebra [16], this significant change was not replicated in Ddhd1 (−/−) plasma from mice (Table 3 and Table S7).

3.3. Focused lipidome Analysis of Human Sera

We performed lipidome analysis on human sera by focusing on nine molecule species significantly changed in the mouse lipidome analyses. Since simple statistical analyses are not applicable due to the unavailability of biological replication of SPG28 patients, we ranked the levels of the seven lipid molecules in the six human samples, including the SPG28 patient. The total amount of PIs observed in the SPG28 patient was within the distribution of the ones in the controls (Figure 4A). Of the nine PIs that showed significant change in Ddhd1 (−/−) mice, seven PIs, 16:0/18:1, 16:0/20:1, 18:1/18:1, 18:1/20:3, 18:1/20:4, 18:1/22:5 and 20:1/20:4, were also detected in human sera. Among these, four PIs, 18:1/18:1, 18:1/20:3 18:1/24:0 and 18:1/22:5, showed the highest amount in the SPG28 patient (Figure 4B). The amounts of all four PIs in SPG28 consistently exceeded the 95% CI calculated from the control samples (Table 4 and Table S8). This indicates that these four PIs were elevated in the SPG28 sera. Notably, all four of these PIs contain oleic acid (18:1) in their side chains. LPA 24:0 was not detected in human sera. The SPG28 serum showed the lowest value for total LPAs, in contrast to the elevation of LPAs observed in mice (Figure 4C and Table S9).

3.4. PIs and LPIs Containing Oleic Acid

Since the four PIs showing the highest values in SPG28 commonly contain oleic acid, we examined the amounts of PIs containing oleic acid in mouse plasma and human sera. Oleic acid-containing PIs were significantly elevated in Ddhd1 (−/−) mice (p = 7.0 × 10−4) and showed the highest values in the SPG28 patient (Figure 5). We examined the ranks of all oleic acid-containing PIs from the results of the lipidome analysis for human sera. Twelve oleic acid-containing PIs were detected in human sera. Nine of them showed the highest values in SPG28. Furthermore, the levels of these nine oleic acid-containing PIs exceeded the 95% CI of the control samples (Table 5 and Table S10). This indicates the tendency of an increase in oleic acid-containing PIs in SPG28 sera.

4. Discussion

We observed a significant increase in the total amount of PI in Ddhd1 (−/−) mice that is consistent with the function of PLA1 converting PI into LPI. The increase in total PI is considered to be due to the accumulation of the substrate PI by the failure of the metabolism of PI to LPI due to the loss of DDHD1 function, which is consistent with the prediction based on the metabolic pathway of DDHD1. PI is one of the essential phospholipid components of the plasma membrane. PI is biased toward the inner side of the plasma membrane [33]. Since DDHD1 is a type of intracellular PLA1 (iPLA), the intracellular amount of PI can be greatly affected by the function of DDHD1 [13]. We observed consistent trends of increase in four PIs, 16:0/18:1, 18:1/18:1, 18:1/20:3 and PI 18:1/22:5, both in Ddhd1 (−/−) plasma from mice and in the human SPG28 serum [23]. Although differences in the content of several phospholipids have been reported between plasma and sera, we assumed the correspondence of quantitative changes in phospholipids between plasma and sera [34]. Notably, all of the PIs that have the highest levels in SPG28 contain oleic acid (18:1) in the side chain, an unsaturated fatty acid with 18 carbons and one double bond. Significant increases in PI 18:1/20:4 have also been observed in the brains of a different strain of Ddhd1 knockout mice previously published [24]. There was no significant increase in oleic acid-containing PI in Ddhd1 (+/−) mice, suggesting that sufficient enzymatic activity of DDHD1 remained (Figure 4A). This is consistent with the mode of inheritance of SPG28 being autosomal recessive. Since the increased levels of oleic acid-containing PIs in the peripheral blood are reasonably attributable to the deficiency of DDHD1, our results indicate the possibility of oleic acid-containing PIs as a blood biomarker for SPG28.
LPA is a phospholipid that is abundant in plasma. LPA is also known to be a potential biomarker for ovarian cancer since LPA in plasma is shown to be specifically increased in ovarian cancer patients [35,36,37]. Although LPA is the product of a reaction catalyzed by DDHD1, the total amount of LPA was increased in Ddhd1 (−/−) mice. Among the eight species of LPAs identified, only LPA 24:0 showed a significant increase, suggesting that the increase in total LPAs is mainly attributable to LPA 24:0. LPA is yielded not only by the PLA1-mediated pathway but also by other pathways involving other enzymes such as lysophospholipase D (LysoPLD) catalyzing lysophoshatidylcholine, lysophosphatidyllethanolamine (LPE) or lysophosphatidylserine (LPS) into LPA [38]. LysoPLD is a secreted enzyme that is a major enzyme-producing PLA in plasma [39]. The excess of LPAs observed in our current analyses is, therefore, attributable to the compensatory pathway such as the one mediated by LysoPLD. Increased activity of LysoPLD and accumulation of LPA are known to occur in female-specific cancers, hepatitis C virus and inflammatory diseases [40,41]. The activity change in LysoPLD is thought to be involved in these diseases in multiple aspects. Increased activity of LysoPLD has been shown to be triggered primarily by inflammatory signals [42]. We therefore speculate that the loss of motoneurons in the pyramidal tract in Ddhd1 knockout mice can act as inflammatory signals, triggering an increase in LysoPLD activity [23]. Since no significant decrease in LPC 24:0 was observed in the Ddhd1 (−/−) mice (Figure S1), the observed increase in LPA 24:0 may be supplied by the metabolism of LPE or LPS, which are out of the scope of our focused lipidome analyses. LPA was at the lowest level in the SPG28 patient in contrast to the observation in mice. This is partly attributable to the differences in biological age between the human and mice at examination. The age of 6 months in the Ddhd1 (−/−) mice is prior to the onset of HSP-like phenotypes and is likely in the middle of the inflammatory process, which may trigger LysoPLD activation without neurological phenotypes. In contrast, since the human sera were collected from the SPG28 patient at the age of 61 after the progression of HSP symptoms, it is possible that the motoneuron degradation prior to the HSP symptoms had reached a plateau and the inflammatory signals might have been weakened.
The significant reduction in LPI 20:4 (sn-2) found in Ddhd1 (−/−) cerebra was not observed in Ddhd1 (−/−) plasma from mice (Table S6) [23]. This discrepancy can be partly attributed to the difference in tissues used for lipid extraction. Plasma is cell-free. DDHD1 is one of the iPLAs and is mainly involved in the intracellular metabolism of phospholipids [43]. Therefore, it is likely that the lipidome analysis using plasma, which is cell-free, does not directly reflect the effects of intracellular DDHD1 metabolism. Since DDHD1 seems to be functionally defective in blood cells, a decrease in LPI 20:4 (sn-2) as in the cerebra is likely to occur.
SPG39, SPG54 and SPG56 are also caused by defects in enzymes involved in PI metabolism, PNPLA6, DDHD2 and CYP2U1, respectively [10]. Notably, DDHD2 is responsible for SPG54 functions such as in iPLA1 as well as in DDHD1 [44,45]. Since abnormal PI metabolism is likely to be the common pathological mechanism for these types of HSPs, peripheral PI levels are good candidates for biomarkers of other HSPs sharing pathogenesis with SPG28.
The major problem with the current study is that statistical analyses were not applicable to human samples due to the unavailability of additional SPG28 samples, since the disease is extremely rare. Although further studies of specificity and sensitivity of oleic acid-containing PIs using additional SPG28 patients are required for establishment as a biomarker, oleic acid-containing PIs may still be useful to group SPGs sharing similar pathogenic pathways and may still help in clarifying their pathogenic mechanisms, potentially contributing to the treatment. Identification of the blood biomarker is useful not only for diagnosis of HSP types in clinical practice but also as a platform for therapeutic drug development. For example, we should be able to evaluate the effect of a drug by measuring oleic acid-containing PI after administering an SPG28 therapeutic agent to patients. To achieve these goals, it is essential to establish standard values for the amount of oleic acid-containing PI. Measurements of oleic acid-containing PI in larger numbers of SPG28 patients and establishment of the standard values are important tasks for this goal.

5. Conclusions

Lipidome analysis of SPG28 model mice and a SPG28 patient was performed to search for a biomarker in blood that is useful for the diagnosis of SPG28. Lipidome analysis of mouse plasma revealed a significant increase in total PI and LPA in Ddhd1 (−/−) mice. We examined the changes in each PI and identified nine PIs that were significantly elevated in Ddhd1 (−/−) mice. Next, we examined the level of these nine PIs and total PA in human serum. Seven of the nine PIs found to be elevated in Ddhd1 (−/−) mice were also detected in human sera. Four of these has the highest level in the SPG28 patient, and all of them were oleic acid-containing PIs (PI 18:1/18:1, PI 18:1/20:3, PI 18:1/20:4 and PI 18:1/22:5). LPA showed lower levels in the SPG28 patient, different from the results observed in Ddhd1 (−/−) mice. From our observation, we propose that oleic acid-containing PI is a promising blood biomarker for SPG28.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biomedicines11041092/s1, Figure S1: LPC 24:0 in plasma from mice Table S1: PCs detected in mouse plasma; Table S2: PIs detected in mouse plasma; Table S3: PIs identified from mice plasma; Table S4: LPCs detected in mouse plasma; Table S5: LPIs detected in mouse plasma; Table S6; LPAs detected in mouse plasma; Table S7: LPAs identified from mice plasma; Table S8: Results in human sera for PIs significantly changed in Ddhd1 (−/−); Table S9: LPAs in human sera; Table S10: Oleic acid-containing PIs in human sera.

Author Contributions

T.M., S.M. and H.S. designed the experiments. K.M. provided the mouse facilities. T.M., K.M., G.H. and H.S. prepared the samples. M.T., Y.I. and T.B. performed lipidome analysis. R.F., S.M. and H.S. supervised the project. T.M. and H.S. wrote the manuscript with contributions from all authors. All authors have read and agreed to the published version of the manuscript.

Funding

S.M. received funding from Grants-in-Aid of MEXT, Japan (#23K06853) to support this research. This work was partly performed in the Cooperative Research Project Program of the Medical Institute of Bioregulation, Kyushu University.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and was approved by the Ethics Committee of Kyushu University (protocol code: 2020-444; date of approval: 9 October 2020) and Kurume University (protocol code: 21057; date of approval: 16 June 2021). The animal study protocol was approved by the Institutional Review Board (or Ethics Committee) of Nakamura Gakuen University (protocol code 2016-1; date of approval: 4 April 2016).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request, while all the analyzed data are included within this article and its supplementary information files.

Conflicts of Interest

The authors declare that they have no competing interests.

References

  1. Tesson, C.; Koht, J.; Stevanin, G. Delving into the complexity of hereditary spastic paraplegias: How unexpected phenotypes and inheritance modes are revolutionizing their nosology. Hum. Genet. 2015, 134, 511. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Klebe, S.; Stevanin, G.; Depienne, C. Clinical and genetic heterogeneity in hereditary spastic paraplegias: From SPG1 to SPG72 and still counting. Rev. Neurol. 2015, 171, 505–530. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Fink, J. Hereditary Spastic Paraplegia: Clinical Principles and Genetic Advances. Semin. Neurol. 2014, 34, 293–305. [Google Scholar] [CrossRef] [PubMed]
  4. Blackstone, C. Hereditary spastic paraplegia. Handb. Clin. Neurol. 2018, 148, 633–652. [Google Scholar] [CrossRef]
  5. Blackstone, C. Cellular Pathways of Hereditary Spastic Paraplegia. Annu. Rev. Neurosci. 2012, 35, 25–47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Fink, J.K. Hereditary spastic paraplegia: Clinico-pathologic features and emerging molecular mechanisms. Acta Neuropathol. 2013, 126, 307–328. [Google Scholar] [CrossRef] [Green Version]
  7. Saputra, L.; Kumar, K.R. Challenges and Controversies in the Genetic Diagnosis of Hereditary Spastic Paraplegia. Curr. Neurol. Neurosci. Rep. 2021, 21, 15. [Google Scholar] [CrossRef]
  8. Elsayed, L.E.O.; Eltazi, I.Z.; Ahmed, A.E.; Stevanin, G. Insights into Clinical, Genetic, and Pathological Aspects of Hereditary Spastic Paraplegias: A Comprehensive Overview. Front. Mol. Biosci. 2021, 8, 690899. [Google Scholar] [CrossRef]
  9. Bouslam, N.; Benomar, A.; Azzedine, H.; Bouhouche, A.; Namekawa, M.; Klebe, S.; Durr, A.; Ruberg, M.; Brice, A.; Yahyaoui, M.; et al. Mapping of a new form of pure autosomal recessive spastic paraplegia (SPG28). Ann. Neurol. 2005, 57, 567–571. [Google Scholar] [CrossRef]
  10. Tesson, C.; Nawara, M.; Salih, M.A.; Rossignol, R.; Zaki, M.S.; Al Balwi, M.; Schule, R.; Mignot, C.; Obre, E.; Bouhouche, A.; et al. Alteration of Fatty-Acid-Metabolizing Enzymes Affects Mitochondrial Form and Function in Hereditary Spastic Paraplegia. Am. J. Hum. Genet. 2012, 91, 1051–1064. [Google Scholar] [CrossRef] [Green Version]
  11. Aloulou, A.; Rahier, R.; Arhab, Y.; Noiriel, A.; Abousalham, A. Phospholipases: An Overview. Methods Mol. Biol. 2018, 1835, 69–105. [Google Scholar] [CrossRef]
  12. Richmond, G.S.; Smith, T.K. Phospholipases A1. Int. J. Mol. Sci. 2011, 12, 588–612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Tani, K.; Kogure, T.; Inoue, H. The intracellular phospholipase A1 protein family. Biomol. Concepts 2012, 3, 471–478. [Google Scholar] [CrossRef] [PubMed]
  14. Higgs, H.N.; Glomset, J.A. Purification and Properties of a Phosphatidic Acid-preferring Phospholipase A1 from Bovine Testis: Examination of the Molecular Basis of Its Activation. J. Biol. Chem. 1996, 271, 10874–10883. [Google Scholar] [CrossRef] [Green Version]
  15. Tani, K.; Mizoguchi, T.; Iwamatsu, A.; Hatsuzawa, K.; Tagaya, M. p125 Is a Novel Mammalian Sec23p-interacting Protein with Structural Similarity to Phospholipid-modifying Proteins. J. Biol. Chem. 1999, 274, 20505–20512. [Google Scholar] [CrossRef] [Green Version]
  16. Nakajima, K.-I.; Mizoguchi, T.; Nagahama, M.; Tagaya, M.; Tani, K.; Sonoda, H.; Aoki, J.; Arai, H. A Novel Phospholipase A1 with Sequence Homology to a Mammalian Sec23p-interacting Protein, p125. J. Biol. Chem. 2002, 277, 11329–11335. [Google Scholar] [CrossRef] [Green Version]
  17. Baba, T.; Kashiwagi, Y.; Arimitsu, N.; Kogure, T.; Edo, A.; Maruyama, T.; Nakao, K.; Nakanishi, H.; Kinoshita, M.; Frohman, M.A.; et al. Phosphatidic Acid (PA)-preferring Phospholipase A1 Regulates Mitochondrial Dynamics. J. Biol. Chem. 2014, 289, 11497. [Google Scholar] [CrossRef] [Green Version]
  18. Mizoguchi, T.; Nakajima, K.-I.; Hatsuzawa, K.; Nagahama, M.; Hauri, H.-P.; Tagaya, M.; Tani, K. Determination of Functional Regions of p125, a Novel Mammalian Sec23p-Interacting Protein. Biochem. Biophys. Res. Commun. 2000, 279, 144–149. [Google Scholar] [CrossRef] [PubMed]
  19. Shimoi, W.; Ezawa, I.; Nakamoto, K.; Uesaki, S.; Gabreski, G.; Aridor, M.; Yamamoto, A.; Nagahama, M.; Tagaya, M.; Tani, K. p125 Is Localized in Endoplasmic Reticulum Exit Sites and Involved in Their Organization. J. Biol. Chem. 2005, 280, 10141–10148. [Google Scholar] [CrossRef] [Green Version]
  20. Liguori, R.; Giannoccaro, M.P.; Arnoldi, A.; Citterio, A.; Tonon, C.; Lodi, R.; Bresolin, N.; Bassi, M.T. Impairment of brain and muscle energy metabolism detected by magnetic resonance spectroscopy in hereditary spastic paraparesis type 28 patients with DDHD1 mutations. J. Neurol. 2014, 261, 1789–1793. [Google Scholar] [CrossRef] [Green Version]
  21. Wu, C.; Fan, D. A Novel Missense Mutation of the DDHD1 Gene Associated with Juvenile Amyotrophic Lateral Sclerosis. Front. Aging Neurosci. 2016, 8, 291. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Miura, S.; Morikawa, T.; Fujioka, R.; Kosaka, K.; Yamada, K.; Hattori, G.; Motomura, M.; Taniwaki, T.; Shibata, H. A novel frameshift mutation of DDHD1 in a Japanese patient with autosomal recessive spastic paraplegia. Eur. J. Med. Genet. 2016, 59, 413–416. [Google Scholar] [CrossRef] [PubMed]
  23. Morikawa, T.; Ohishi, H.; Kosaka, K.; Shimojo, T.; Nagano, A.; Taniguchi, I.; Fujioka, R.; Moriyama, K.; Unoki, M.; Takahashi, M.; et al. Ddhd1 knockout mouse as a model of locomotive and physiological abnormality in familial spastic paraplegia. Biosci. Rep. 2021, 41, BSR20204171. [Google Scholar] [CrossRef] [PubMed]
  24. Inloes, J.M.; Jing, H.; Cravatt, B.F. The Spastic Paraplegia-Associated Phospholipase DDHD1 Is a Primary Brain Phosphatidylinositol Lipase. Biochemistry 2018, 57, 5759–5767. [Google Scholar] [CrossRef]
  25. Skarnes, W.C.; Rosen, B.; West, A.P.; Koutsourakis, M.; Bushell, W.; Iyer, V.; Mujica, A.O.; Thomas, M.; Harrow, J.; Cox, T.; et al. A conditional knockout resource for the genome-wide study of mouse gene function. Nature 2011, 474, 337–344. [Google Scholar] [CrossRef] [Green Version]
  26. Scriver, C.R. ThePAH gene, phenylketonuria, and a paradigm shift. Hum. Mutat. 2007, 28, 831–845. [Google Scholar] [CrossRef]
  27. Fernández-Cañón, J.M.; Granadino, B.; De Bernabé, D.B.; Renedo, N.; Fernández-Ruiz, E. The molecular basis of alkaptonuria. Nat. Genet. 1996, 14, 19–24. [Google Scholar] [CrossRef]
  28. Mistry, J.B.; Bukhari, M.; Taylor, A.M. Alkaptonuria. Rare Dis. 2013, 1, e27475. [Google Scholar] [CrossRef] [Green Version]
  29. Blau, N. Genetics of Phenylketonuria: Then and Now. Hum. Mutat. 2016, 37, 508–515. [Google Scholar] [CrossRef] [Green Version]
  30. Takeda, H.; Izumi, Y.; Takahashi, M.; Paxton, T.; Tamura, S.; Koike, T.; Yu, Y.; Kato, N.; Nagase, K.; Shiomi, M.; et al. Widely-targeted quantitative lipidomics method by supercritical fluid chromatography triple quadrupole mass spectrometry. J. Lipid Res. 2018, 59, 1283–1293. [Google Scholar] [CrossRef] [Green Version]
  31. McDermott, M.; Wang, Y.; Wakelam, M.; Bankaitis, V. Mammalian phospholipase D: Function, and therapeutics. Prog. Lipid Res. 2020, 78, 101018. [Google Scholar] [CrossRef] [PubMed]
  32. Yamashita, A.; Kumazawa, T.; Koga, H.; Suzuki, N.; Oka, S.; Sugiura, T. Generation of lysophosphatidylinositol by DDHD domain containing 1 (DDHD1): Possible involvement of phospholipase D/phosphatidic acid in the activation of DDHD1. Biochim. Biophys. Acta 2010, 1801, 711–720. [Google Scholar] [CrossRef]
  33. Zachowski, A. Phospholipids in animal eukaryotic membranes: Transverse asymmetry and movement. Biochem. J. 1993, 294, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Ishikawa, M.; Maekawa, K.; Saito, K.; Senoo, Y.; Urata, M.; Murayama, M.; Tajima, Y.; Kumagai, Y.; Saito, Y. Plasma and Serum Lipidomics of Healthy White Adults Shows Characteristic Profiles by Subjects’ Gender and Age. PLoS ONE 2014, 9, e91806. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Xu, Y.; Shen, Z.; Wiper, D.W.; Wu, M.; Morton, R.E.; Elson, P.; Kennedy, A.W.; Belinson, J.; Markman, M.; Casey, G. Lysophosphatidic Acid as a Potential Biomarker for Ovarian and Other Gynecologic Cancers. JAMA 1998, 280, 719–723. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Sedláková, I.; Vávrová, J.; Tošner, J.; Hanousek, L. Lysophosphatidic acid (LPA)—A perspective marker in ovarian cancer. Tumor Biol. 2010, 32, 311–316. [Google Scholar] [CrossRef]
  37. Li, Y.-Y.; Zhang, W.-C.; Zhang, J.-L.; Zheng, C.-J.; Zhu, H.; Yu, H.-M.; Fan, L.-M. Plasma levels of lysophosphatidic acid in ovarian cancer versus controls: A meta-analysis. Lipids Health Dis. 2015, 14, 72. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Aoki, J.; Taira, A.; Takanezawa, Y.; Kishi, Y.; Hama, K.; Kishimoto, T.; Mizuno, K.; Saku, K.; Taguchi, R.; Arai, H. Serum Lysophosphatidic Acid Is Produced through Diverse Phospholipase Pathways. J. Biol. Chem. 2002, 277, 48737–48744. [Google Scholar] [CrossRef] [Green Version]
  39. Ninou, I.; Magkrioti, C.; Aidinis, V. Autotaxin in Pathophysiology and Pulmonary Fibrosis. Front. Med. 2018, 5, 180. [Google Scholar] [CrossRef] [Green Version]
  40. Watanabe, N.; Ikeda, H.; Nakamura, K.; Ohkawa, R.; Kume, Y.; Aoki, J.; Hama, K.; Okudaira, S.; Tanaka, M.; Tomiya, T.; et al. Both Plasma Lysophosphatidic Acid and Serum Autotaxin Levels are Increased in Chronic Hepatitis C. J. Clin. Gastroenterol. 2007, 41, 616–623. [Google Scholar] [CrossRef]
  41. Magkrioti, C.; Galaris, A.; Kanellopoulou, P.; Stylianaki, E.-A.; Kaffe, E.; Aidinis, V. Autotaxin and chronic inflammatory diseases. J. Autoimmun. 2019, 104, 102327. [Google Scholar] [CrossRef]
  42. Liu, S.; Murph, M.; Panupinthu, N.; Mills, G.B. ATX-LPA receptor axis in inflammation and cancer. Cell Cycle 2009, 8, 3695–3701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Inoue, H.; Baba, T.; Sato, S.; Ohtsuki, R.; Takemori, A.; Watanabe, T.; Tagaya, M.; Tani, K. Roles of SAM and DDHD domains in mammalian intracellular phospholipase A1 KIAA0725p. Biochim. Biophys. Acta 2012, 1823, 930–939. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Schuurs-Hoeijmakers, J.H.; Geraghty, M.T.; Kamsteeg, E.-J.; Ben-Salem, S.; de Bot, S.T.; Nijhof, B.; van de Vondervoort, I.I.; van der Graaf, M.; Nobau, A.C.; Otte-Höller, I.; et al. Mutations in DDHD2, Encoding an Intracellular Phospholipase A1, Cause a Recessive Form of Complex Hereditary Spastic Paraplegia. Am. J. Hum. Genet. 2012, 91, 1073–1081. [Google Scholar] [CrossRef] [Green Version]
  45. Doi, H.; Ushiyama, M.; Baba, T.; Tani, K.; Shiina, M.; Ogata, K.; Miyatake, S.; Fukuda-Yuzawa, Y.; Tsuji, S.; Nakashima, M.; et al. Late-onset spastic ataxia phenotype in a patient with a homozygous DDHD2 mutation. Sci. Rep. 2014, 4, 7132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. (A) The structure of Ddhd1 gene (NM_176845.5). Exons including the region of the DDHD domain that is critical for enzymatic activity are shown in green. (B) Ddhd1 knockout mice established previously. Strain 1 and 2 deleted the Ddhd1 gene by inserting a cassette of lacZ and neomycin resistance gene inside an intron [17,24]. In strain 3, Ddhd1 is disrupted by inducing 5 bp deletion to induce premature termination [23]. The map is not drawn to scale.
Figure 1. (A) The structure of Ddhd1 gene (NM_176845.5). Exons including the region of the DDHD domain that is critical for enzymatic activity are shown in green. (B) Ddhd1 knockout mice established previously. Strain 1 and 2 deleted the Ddhd1 gene by inserting a cassette of lacZ and neomycin resistance gene inside an intron [17,24]. In strain 3, Ddhd1 is disrupted by inducing 5 bp deletion to induce premature termination [23]. The map is not drawn to scale.
Biomedicines 11 01092 g001
Figure 2. Comparison of total amount of PC and PI. Total amount of PC (A) and PI (B) in plasma from mice (n = 6 each). Ddhd1 (+/+) mice and Ddhd1 (−/−) mice are shown in opened and filled columns, respectively. The unit of the vertical axis is pmol/mg. Error bars indicate mean ± SD. ns: not significant, *: p < 0.05. All data were analyzed by Student’s t test.
Figure 2. Comparison of total amount of PC and PI. Total amount of PC (A) and PI (B) in plasma from mice (n = 6 each). Ddhd1 (+/+) mice and Ddhd1 (−/−) mice are shown in opened and filled columns, respectively. The unit of the vertical axis is pmol/mg. Error bars indicate mean ± SD. ns: not significant, *: p < 0.05. All data were analyzed by Student’s t test.
Biomedicines 11 01092 g002
Figure 3. Comparison of total amount of LPC, LPI and LPA. Amount of LPC (A), LPI (B) and LPA (C) in plasma from mice. LPA is shown in only total amount due to the inseparability of sn-1 and sn-2 in our analyses. Ddhd1 (+/+) mice and Ddhd1 (−/−) mice are shown in opened and filled columns, respectively. The unit of the vertical axis is pmol/mg. Error bars indicate mean ± SD. ns: not significant, *: p < 0.05. All data were analyzed by Student’s t test.
Figure 3. Comparison of total amount of LPC, LPI and LPA. Amount of LPC (A), LPI (B) and LPA (C) in plasma from mice. LPA is shown in only total amount due to the inseparability of sn-1 and sn-2 in our analyses. Ddhd1 (+/+) mice and Ddhd1 (−/−) mice are shown in opened and filled columns, respectively. The unit of the vertical axis is pmol/mg. Error bars indicate mean ± SD. ns: not significant, *: p < 0.05. All data were analyzed by Student’s t test.
Biomedicines 11 01092 g003
Figure 4. Lipidome analysis using human sera of the SPG28 patient and five controls. (A) Total PIs in human sera. (B) Focused analyses of human sera on PI species significantly changed in Ddhd1 (−/−) plasma from mice. PIs arranged in the order of decreasing q values obtained from lipidome analysis using plasma from mice. Individually enlarged charts are shown below. Red arrows indicate the PI species with the highest levels in the SPG28 patient. (C) Total LPAs in human sera. The data points of SPG28 are represented by red squares. The unit of the vertical axis is pmol/mg.
Figure 4. Lipidome analysis using human sera of the SPG28 patient and five controls. (A) Total PIs in human sera. (B) Focused analyses of human sera on PI species significantly changed in Ddhd1 (−/−) plasma from mice. PIs arranged in the order of decreasing q values obtained from lipidome analysis using plasma from mice. Individually enlarged charts are shown below. Red arrows indicate the PI species with the highest levels in the SPG28 patient. (C) Total LPAs in human sera. The data points of SPG28 are represented by red squares. The unit of the vertical axis is pmol/mg.
Biomedicines 11 01092 g004
Figure 5. PIs and LPAs containing oleic acid. (A) Total oleic acid-containing PIs in plasma from mice. ns: not significant. The data were analyzed by Dunnett’s multiple comparison test. (B) Total oleic acid-containing PIs in human sera. The data point of SPG28 is represented by a red square. The unit of the vertical axis is pmol/mg. Error bars indicate mean ± SD. *: p < 0.05.
Figure 5. PIs and LPAs containing oleic acid. (A) Total oleic acid-containing PIs in plasma from mice. ns: not significant. The data were analyzed by Dunnett’s multiple comparison test. (B) Total oleic acid-containing PIs in human sera. The data point of SPG28 is represented by a red square. The unit of the vertical axis is pmol/mg. Error bars indicate mean ± SD. *: p < 0.05.
Biomedicines 11 01092 g005
Table 1. Observed phenotypes of three different Ddhd1 knockout mice [17,23,24].
Table 1. Observed phenotypes of three different Ddhd1 knockout mice [17,23,24].
StrainPhenotypes
Strain 1
(Baba, 2014 [17])
・Shortened mitochondrial sheath in sperm
・Impaired sperm motility
Strain 2
(Inloes, 2018 [24])
・Accumulation of PI 18:1/20:4
・Decrease of LPI 20:4
Strain 3
(Morikawa, 2021 [23])
・Gait disturbance-like symptpms
・Decrease of LPI 20:4 (sn-2) in cerebra
・Axonal decrease in the pyramidal tract
Table 2. PIs identified from mice plasma.
Table 2. PIs identified from mice plasma.
PIsRetention Time (min)Ddhd1 (+/+)Ddhd1 (−/−)p-Valueq-Value
Average±SDAverage±SD
PI 18:1/20:417.19 ± 00.6450.1311.4950.2710.00010.006*
PI 20:1/20:417.08 ± 0.010.01410.0050.04100.0110.00050.017*
PI 18:1/20:317.13 ± 00.1530.0430.4350.1290.00090.020*
PI 18:1/22:517.2 ± 0.010.02580.0130.07820.0270.00290.046*
PI 19:1/20:417.13 ± 0.010.007770.0030.02390.0090.00340.046*
PI 16:0/20:116.94 ± 0.010.01410.0060.03670.0120.00410.046*
PI 15:0/24:017.28 ± 00.05240.0300.16680.0620.00410.046*
PI 18:1/18:117.02 ± 00.1850.0780.3700.0850.00520.046*
PI 16:0/18:117.03 ± 00.2140.0650.4370.1260.00570.046*
PI 18:1/22:617.28 ± 00.08460.0490.2380.0920.00790.051
PI 16:0/18:217.15 ± 00.8290.2191.610.4930.00880.051
PI 18:2/20:117.03 ± 0.010.01310.0060.02840.0090.0100.053
PI 16:0/18:017.01 ± 0.020.01960.0060.03750.0110.0110.055
PI 17:0/20:317.09 ± 0.020.02050.0060.03430.0080.0120.055
PI 18:1/20:517.3 ± 00.2010.0940.3970.1110.0130.055
PI 18:1/18:217.14 ± 00.3760.1500.7510.2390.0140.056
PI 16:0/16:017.08 ± 0.010.01940.0070.02960.0040.0160.061
PI 16:0/20:417.2 ± 01.200.2332.250.7880.0170.061
PI 20:1/22:616.23 ± 3.930.001580.0010.005710.0030.0210.069
PI 16:0/22:617.28 ± 00.3080.1640.7490.3250.0220.070
PI 15:0/22:517.24 ± 00.02430.0070.04710.0170.0220.070
PI 18:1/18:317.23 ± 0.010.008160.0060.01870.0070.0230.070
PI 16:0/20:317.13 ± 00.3180.0980.6550.2690.0250.070
PI 17:0/18:217.1 ± 0.010.03270.0100.05800.0200.0290.079
PI 16:1/18:217.26 ± 0.010.02280.0180.05570.0230.0320.081
PI 16:0/22:517.2 ± 00.08880.0480.2270.1150.0320.081
PI 16:1/20:417.31 ± 0.010.01050.0060.03380.0200.0350.082
PI 18:0/18:016.92 ± 0.010.02720.0070.04190.0120.0370.084
PI 18:0/20:417.07 ± 05.600.7567.111.1740.0370.084
PI 18:0/20:216.99 ± 0.010.09750.0280.15710.0480.0370.084
PI 18:0/18:116.93 ± 00.2260.0610.3370.0840.0370.084
PI 18:0/20:317.01 ± 0.011.320.3501.910.4340.0410.084
PI 15:0/22:216.98 ± 00.02650.0100.04690.0170.0420.084
PI 18:2/23:017.16 ± 0.020.001200.0010.002680.0010.0420.084
PI 18:0/20:116.84 ± 0.020.01560.0060.02790.0110.0460.084
PI 18:2/20:417.29 ± 0.020.02740.0140.05430.0230.0480.086
PI 18:0/18:217.02 ± 01.090.2651.550.3830.0500.086
PI 16:0/18:317.24 ± 0.010.01680.0070.03410.0170.0580.097
PI 18:0/18:317.12 ± 0.010.01670.0070.02650.0070.0590.097
PI 18:2/25:017.06 ± 0.050.002250.0010.003280.0010.0590.097
PI 19:0/20:417 ± 00.08090.0290.1310.0460.0640.10
PI 16:1/18:117.15 ± 00.01820.0120.03930.0190.0660.10
PI 16:0/20:517.31 ± 00.3890.1700.7640.3700.0670.10
PI 16:0/16:117.16 ± 0.010.01980.0110.03890.0180.0750.11
PI 22:1/22:616.32 ± 3.950.0009250.0010.002800.0020.0780.11
PI 19:0/22:617.12 ± 0.020.004920.0030.00950.0050.0800.11
PI 18:0/22:617.16 ± 00.7120.2861.1220.3770.0810.11
PI 18:2/19:016.97 ± 00.01800.0090.03230.0140.0900.12
PI 18:2/18:217.25 ± 00.07990.0510.16490.0880.0910.12
PI 15:0/22:017.41 ± 0.020.01160.0100.02320.0110.110.14
PI 17:0/20:417.13 ± 00.09690.0140.13050.0410.120.14
PI 17:0/18:116.98 ± 0.010.00920.0030.01490.0070.120.15
PI 20:0/20:416.96 ± 0.010.03020.0230.05680.0270.120.15
PI 15:0/22:117 ± 0.220.007410.0040.01260.0060.130.15
PI 18:0/22:416.97 ± 00.03250.0110.04780.0170.130.15
PI 19:0/20:316.96 ± 0.010.01540.0090.02680.0130.130.15
PI 18:0/22:517.09 ± 0.010.1810.0830.2780.1030.130.15
PI 15:0/24:117.37 ± 0.020.00670.0060.01400.0090.160.18
PI 16:1/18:017.02 ± 0.010.01750.0070.02510.0080.160.18
PI 18:0/20:517.18 ± 01.100.2291.330.3240.220.23
PI 17:0/20:517.25 ± 00.02590.0110.03650.0150.230.24
PI 20:4/22:117.21 ± 0.010.004810.0020.00650.0030.280.29
PI 20:0/22:614.22 ± 6.350.001530.0020.002330.0010.460.47
PI 20:5/24:117.19 ± 0.030.003250.0010.003050.0020.840.84
The unit of measurement is pmol/mg. The q value was calculated by Student’s t-= test followed by correction by the Benjamini–Hochberg method. * q < 0.05.
Table 3. LPAs identified from mice plasma.
Table 3. LPAs identified from mice plasma.
LPAsRetention TimeDdhd1 (+/+) Ddhd1 (−/−) p-Valueq-Value
(min)AverageSDAverageSD
LPA 24:022.89 ± 0.041.170.473.120.450.00010.0004*
LPA 18:022.89 ± 0.030.3290.040.4370.070.0150.061
LPA 22:622.89 ± 0.131.480.181.380.100.290.78
LPA 22:5 22.84 ± 0.060.9050.100.8790.040.61>0.99
LPA 20:222.89 ± 0.020.2750.070.2670.040.83>0.99
LPA 20:422.32 ± 0.060.08250.020.07860.030.84>0.99
LPA 20:322.89 ± 0.030.7780.180.7960.110.86>0.99
LPA 18:222.89 ± 0.030.5460.240.5310.100.90>0.99
The unit of measurement is pmol/mg. The q value was calculated by Student’s t test followed by correction by the Benjamini–Hochberg method. * q < 0.05.
Table 4. Results in human sera for PIs significantly changed in Ddhd1 (−/−).
Table 4. Results in human sera for PIs significantly changed in Ddhd1 (−/−).
PIsRetention Time (min)SPG28Control 1Control 2Control 3Control 4Control 5Rank of SPG28 among
All Samples
Average of ConrolsSD of Conrols95% CISPG28
out of 95% CI?
PI 18:1/20:4 *17.20 ± 0.013.332.601.581.622.330.79711.790.6340.556(1.23–2.34)Yes
PI 20:1/20:417.10 ± 0.020.09630.05340.05490.06840.13230.021220.06600.0370.0321(0.0340–0.0981)No
PI 18:1/20:3 *17.17 ± 0.010.5890.4890.2530.2570.3530.12310.2950.1210.106(0.189–0.401)Yes
PI 18:1/22:5 *17.21 ± 0.010.1430.1340.09280.07300.08920.032510.08430.03270.0287(0.0556–0.113)Yes
PI 16:0/20:116.97 ± 0.010.09710.06610.05860.10230.12610.017830.07420.03740.0328(0.0414–0.107)No
PI 18:1/18:1 *17.04 ± 0.012.852.271.832.331.820.66511.780.5980.524(1.26–2.31)Yes
PI 16:0/18:117.06 ± 0.016.086.405.287.254.591.8535.071.851.62(3.45–6.70)No
PIs are arranged in describing order of q value obtained from lipidome analysis in mice. The unit of measurement is pmol/mg. * PIs that have the highest level in the SPG28 patient.
Table 5. Oleic acid-containing PIs in human sera.
Table 5. Oleic acid-containing PIs in human sera.
PIsRetention Time (min)SPG28Control 1Control 2Control 3Control 4Control 5Rank of SPG28
among All Samples
Average of Conrols±SD of Conrols95% CISPG28 out of 95% CI?
PI 16:1/18:117.17 ± 0.010.3740.3620.3060.2520.4150.081920.2830.1140.100(0.183–0.384)No
PI 16:0/18:117.06 ± 0.016.086.405.287.254.591.8535.071.851.62(3.45–6.70)No
PI 17:0/18:117.00 ± 0.010.2880.3160.1790.2120.2110.095120.2030.07100.0622(0.141–0.265)No
PI 18:3/18:117.28 ± 0.020.1250.06200.08260.02810.03160.017210.04430.02430.0213(0.0230–0.0656)Yes
PI 18:1/18:217.17 ± 0.013.251.932.331.592.190.80311.770.5440.477(1.29–2.24)Yes
PI 18:1/18:117.04 ± 0.012.852.271.832.331.820.66511.780.5980.524(1.26–2.31)Yes
PI 18:0/18:116.95 ± 0.0111.29.4710.09.286.503.1717.682.572.25(5.44–9.93)Yes
PI 18:1/20:517.32 ± 0.020.1070.07140.09160.04440.03620.020010.05270.02560.0224(0.0303–0.0751)Yes
PI 18:1/20:417.20 ± 0.013.332.601.581.622.330.79711.790.6340.556(1.23–2.34)Yes
PI 18:1/20:317.17 ± 0.010.5890.4890.2530.2570.3530.12310.2950.1210.106(0.189–0.401)Yes
PI 18:1/22:617.30 ± 0.010.2500.1380.1320.1610.1560.050210.1270.04000.0351(0.0922–0.162)Yes
PI 18:1/22:517.21 ± 0.010.1430.1340.09280.07300.08920.032510.08430.03270.0287(0.0556–0.113)Yes
Total Oleic acid-containing PI 28.624.322.223.118.77.70119.26.035.29(13.9–24.5)Yes
Total PI 13413814212212667.0311927.123.7(95.2–143)No
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Morikawa, T.; Takahashi, M.; Izumi, Y.; Bamba, T.; Moriyama, K.; Hattori, G.; Fujioka, R.; Miura, S.; Shibata, H. Oleic Acid-Containing Phosphatidylinositol Is a Blood Biomarker Candidate for SPG28. Biomedicines 2023, 11, 1092. https://doi.org/10.3390/biomedicines11041092

AMA Style

Morikawa T, Takahashi M, Izumi Y, Bamba T, Moriyama K, Hattori G, Fujioka R, Miura S, Shibata H. Oleic Acid-Containing Phosphatidylinositol Is a Blood Biomarker Candidate for SPG28. Biomedicines. 2023; 11(4):1092. https://doi.org/10.3390/biomedicines11041092

Chicago/Turabian Style

Morikawa, Takuya, Masatomo Takahashi, Yoshihiro Izumi, Takeshi Bamba, Kosei Moriyama, Gohsuke Hattori, Ryuta Fujioka, Shiroh Miura, and Hiroki Shibata. 2023. "Oleic Acid-Containing Phosphatidylinositol Is a Blood Biomarker Candidate for SPG28" Biomedicines 11, no. 4: 1092. https://doi.org/10.3390/biomedicines11041092

APA Style

Morikawa, T., Takahashi, M., Izumi, Y., Bamba, T., Moriyama, K., Hattori, G., Fujioka, R., Miura, S., & Shibata, H. (2023). Oleic Acid-Containing Phosphatidylinositol Is a Blood Biomarker Candidate for SPG28. Biomedicines, 11(4), 1092. https://doi.org/10.3390/biomedicines11041092

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop