Modification of erythrocyte membrane phospholipid composition in preterm newborns with retinopathy of prematurity: The omegaROP study

N-3 polyunsaturated fatty acids (PUFAs) may prevent retinal vascular abnormalities observed in oxygen-induced retinopathy, a model of retinopathy of prematurity (ROP). In the OmegaROP prospective cohort study, we showed that preterm infants who will develop ROP accumulate the n-6 PUFA arachidonic acid (ARA) at the expense of the n-3 PUFA docosahexaenoic acid (DHA) in erythrocytes with advancing gestational age (GA). As mice lacking plasmalogens ―That are specific phospholipids considered as reservoirs of n-6 and n-3 PUFAs― Display a ROP-like phenotype, the aim of this study was to determine whether plasmalogens are responsible for the changes observed in subjects from the OmegaROP study. Accordingly, preterm infants aged less than 29 weeks GA were recruited at birth in the Neonatal Intensive Care Unit of University Hospital Dijon, France. Blood was sampled very early after birth to avoid any nutritional influence on its lipid composition. The lipid composition of erythrocytes and the structure of phospholipids including plasmalogens were determined by global lipidomics using liquid chromatography coupled to high-resolution mass spectrometry (LC-HRMS). LC-HRMS data confirmed our previous observations by showing a negative association between the erythrocyte content in phospholipid esterified to n-6 PUFAs and GA in infants without ROP (rho = −0.485, p = 0.013 and rho = −0.477, p = 0.015 for ethanolamine and choline total phospholipids, respectively). Phosphatidylcholine (PtdCho) and phosphatidylethanolamine (PtdEtn) species with ARA, namely PtdCho16:0/20:4 (rho = −0.511, p < 0.01) and PtdEtn18:1/20:4 (rho = −0.479, p = 0.015), were the major contributors to the relationship observed. On the contrary, preterm infants developing ROP displayed negative association between PtdEtn species with n-3 PUFAs and GA (rho = −0.380, p = 0.034). They were also characterized by a positive association between GA and the ratio of ethanolamine plasmalogens (PlsEtn) with n-6 PUFA to PlsEtn with n-3 PUFAs (rho = 0.420, p = 0.029), as well as the ratio of PlsEtn with ARA to PlsEtn with DHA (rho = 0.843, p = 0.011). Altogether, these data confirm the potential accumulation of n-6 PUFAs with advancing GA in erythrocytes of infants developing ROP. These changes may be partly due to plasmalogens.

Introduction

Retinopathy of prematurity (ROP) is the leading cause of childhood blindness with an estimated incidence ranging from 6 to 34% in developed countries (Good & Hardy, 2001; Jordan, 2014). ROP is characterized by a first phase of vaso-obliteration in the central retina (phase 1), followed by the overexpression of pro-angiogenic growth factors such angiopoietins and VEGF (Sonmez, Drenser, Capone, & Trese, 2008; Sato, Shima, & Kusaka, 2011; Rivera et al., 2017) associated to neovascular events in the retina (phase 2) (Hellstrom, Smith, & Dammann, 2013). As a result, mature retinal vessels exhibit several major abnormalities, such as increased dilatation and tortuosity, as well as vascular leakage (Hellstrom et al., 2013).

Previous studies have reported that the polyunsaturated fatty acids (PUFAs) such as arachidonic acid (ARA, C20:4 n-6) and docosahexaenoic acid (DHA, C22:6 n-3) influence retinal vascularization processes in mouse models of oxygen-induced retinopathy, a mouse model of ROP (Connor et al., 2007; Sapieha et al., 2011). Moreover, human studies revealed alterations in blood levels of PUFAs in preterm newborns developing ROP (Martin et al., 2011; Lofqvist et al., 2018; Pallot et al., 2019). Particularly in the OmegaROP study, we have shown an accumulation of ARA at the expense of DHA in erythrocytes of preterm infants that will develop ROP. In cell membranes, PUFAs such as ARA and DHA are esterified on membrane phospholipids, from which they can be released by phospholipases for further intracellular metabolization and/or signaling. Within these phospholipids, plasmalogens represent a particular sub-class characterized by the presence of a vinyl-ether bond at sn-1 position of glycerol instead of an ester bond. Plasmalogens are considered as “reservoirs” for PUFAs such as ARA and DHA (Nagan & Zoeller, 2001) and are abundant in the human retina (Bretillon et al., 2008; Berdeaux et al., 2010; Acar et al., 2012). Interestingly, we have shown that plasmalogen-deficient mice exhibit retinal vascular abnormalities resembling to those observed in ROP (Saab-Aoude, Bron, Creuzot-Garcher, Bretillon, & Acar, 2013; Saab et al., 2014). Indeed, retinal vascular development in these mice was characterized by a delayed outgrowth followed by increased angiogenesis associated to the overexpression of pro-angiogenic factors such as angiopoietins (Saab et al., 2014).

In this work, we aimed to assess whether the differential accumulation of ARA and DHA in preterm infants of the OmegaROP study is associated to plasmalogens. By using high-resolution mass spectrometry (HRMS), we evaluated the concentrations of individual phospholipids species in erythrocytes of preterm infants that will or not develop ROP. In this study, blood was collected immediately after birth, in order to limit the interference with lipid nutritional intakes.

Materials and methods

Ethics statement

This study was conducted in accordance with the guidelines of the Declaration of Helsinki. The experimental procedures were approved by local ethics committee (CPP Est III, School of Medicine, Dijon, France) that waived the obtainment of a written consent as our study did not generate additional procedure to those of standard care. Instead of, an information note was given to parents and/or legal guardians. In accordance with “ethical considerations for clinical trials on medicinal products conducted with the pediatric population”, the volume of blood collected in preterm infants was limited to 0.5 mL (European Union. 2008).

Selection of the patients

As described previously (Pallot et al., 2019), all preterm infants born before 29 weeks GA and hospitalized in the neonatal intensive care unit of the Dijon University Hospital, Dijon, France, between 31 July 2015 and 31 January 2018 were included in the study. A 0.5-mL blood sample was collected by venipuncture in a heparinized tube within the first 48 h of life. Red blood cells were immediately separated from plasma by centrifugation at 1860 × g at 4°C. The red blood cell pellet was then washed three times with an isotonic saline solution. Samples were stored at −80°C until lipidomic analyses.

ROP screening was performed with the wide-field RETCAM II^® camera (Clarity Medical Systems; Pleasanton, CA, United States) using a lid speculum after topical anesthesia by chlorhydrate oxybuprocaine, 1.6 mg/0.4 mL (Théa, Clermont-Ferrand, France). Pupillary dilation was previously performed using one drop of tropicamide, 2 mg/0.4 mL, Théa, Clermont-Ferrand, France). The procedure was completed by a trained nurse and all fundus photographs were analyzed by two trained pediatrics-specialized ophthalmologists. Screening began at 4–6 weeks of life but never before 31 weeks of postconceptional age (PCA). Fundus imaging was repeated every other week until 39 weeks PCA if no ROP was detected, and every week or up to twice a week in case of ROP. ROP staging was determined according to the International Classification of ROP (International Committee for the Classification of Retinopathy of, 2005). Subjects were classified in the group suffering from ROP (ROP group) or in the group of unaffected controls (no-ROP group). Within the ROP group, subjects were classified into type 1 ROP and type 2 ROP. The major risks of developing ROP, namely gestational age (GA), weight at birth, duration of mechanical ventilation, sepsis, use of erythropoietin, red blood cell transfusion and cerebral hemorrhage were documented.

Characterization and quantification of individual phospholipid species

Total lipids were extracted from red blood cells according to Moilanen and Nikkari by using chloroform/methanol (1:1, v:v) (Moilanen & Nikkari, 1981). The phosphorus content of the total lipid extracts was determined according to the method developed by Bartlett and Lewis (Bartlett & Lewis, 1970). The samples were then diluted at a concentration of 500 μg/mL in chloroform/methanol 1:1 (v/v). Quality control (QC) were prepared by pooling 10 µL of each resuspended lipid extract. The concentrations of individual phospholipids species of erythrocytes total lipids were determined by Hydrophilic Interaction Liquid Chromatography coupled to High Resolution Mass Spectrometry (HILIC-HR-MS).

Liquid chromatography analyses were performed using a Dionex UltiMate™ 3000 LC pump from Thermo Scientific (San Jose, CA, United States) equipped with an autosampler. The injection volume was 10 μL. Separation of lipid classes was achieved under HILIC conditions using an Accucore HILIC column (150 mm × 2.1 mm i. d., 2.6 µm, Thermo). The column was maintained at 40°C. The mobile phases consisted of (A) ACN/H₂O (99:1, v/v) containing 10 mM ammonium acetate and (B) ACN/H₂O (50:50, v/v) containing 10 mM ammonium acetate. The solvent-gradient system of the analytical pump was as follows: 0 min 100% A, 10 min 92% A, 40 min 50% A, 41–60 min 100% A. The flow rate was set to 500 μL.min−1. In order to guarantee analytical performance, Quality QC were used every eight test sample. HR-MS analyses of phospholipids were carried out using the Orbitrap FusionTM (Thermo Scientific, United States) Mass Spectrometer equipped with an EASY-MAX NGTM Ion Source (H-ESI). H-ESI source parameters were optimized and set as follows: Ion transfer tube temperature of 285°C, sheath gas flow rate of 35 au, auxiliary gas flow rate of 25 au, sweep gas of one au, and vaporizer temperature of 370°C. Positive and negative ions were monitored alternatively by switching polarity approach with a spray voltage set to 3500 V in positive and negative ion modes. The Orbitrap mass analyzer was employed to obtain all mass spectra in full scan mode with a mass range of 200–2000 Da, and a target resolution of 120,000 (FWHM at m/z 200). All MS data were recorded using a max injection time of 50 ms, automated gain control (AGC) at 4.105 and one microscan. An Intensity Threshold filter of 1.103 counts was applied. For MS/MS analyses, High-energy Collisional Dissociation (HCD) was employed for the fragmentation of PL species with optimized stepped collision energy of 30% (±5%). The linear ion trap (LIT) was used to acquire spectra for fragment ions in data-dependent mode. The AGC target was set to 2.104 with a max injection time of 50 ms. All MS and MS/MS data were acquired in the profile mode. The Orbitrap Fusion was controlled by XcaliburTM 4.1 software. The identification of PL species was performed, using the data of high accuracy and the information collected from fragmentation spectra (tolerance 5 ppm for MS1 and 20 ppm for MS2), with the help of the LipidSearchTM software and the LIPID MAPS^® database (https://www.lipidmaps.org/).

Relative quantification of the abundances of choline glycerophospholipids (ChoGpl) [including phosphatidyl cholines (PtdCho) and plasmenyl cholines (choline plasmalo gens or PlsCho)] and ethanolamine glycerophospholipids (EtnGpl) [including phosphatidylethanolamines (PtdEtn) and plasmenylethanolamines (ethanolamine plasmalogens or PlsEtn)] molecular species between samples was performed in the high resolution MS1 mode (positive for ChoGpl, negative for EtnGpl) by normalization of targeted phospholipid ion peak areas to the PtdCho14:0/14:0 or PtdEtn14:0/14:0 internal standards, respectively. Due to the lack of available lipid standards representing individual molecular species of EtnGpl and ChoGpl, the abundances of ChoGpl and EtnGpl molecular species was reported as the percentage of the total ChoGpl or EtnGpl ion abundance, respectively (after normalization on the PtdCho14:0/14:0 and PtdEtn14:0/14:0 internal standards).

Statistical analyses

Statistical analysis was performed using GraphPad Prism v6.05 (GraphPad Software, San Diego, CA, United States), XLSTAT v2018.02.50494 (Addinsoft, Paris, France), and R Project v4.0.2 (Revolutions). Quantitative data were expressed as median and interquartile range [IQR]. The groups were compared using the nonparametric Mann-Whitney test for quantitative variables and the Chi-2 test or Fisher exact test for qualitative variables. The Benjamini-Hochberg false discovery adjusted p-value was applied to correct for multiple testing. Linear regression analyses were used to determine the R-squared (R²) correlation coefficient values. Spearman correlations were carried out to compare the levels of individual or total levels of phospholipids and ratios as a function of GA. A p-value lower than 0.05 was considered as statistically significant and the tests were two-tailed.

Results

Characteristics of the population

The characteristics of the population are presented in Table 1. As described before (Pallot et al., 2019), fifty-eight preterm infants born before 29 weeks of GA were included in the study. Six infants died and the mortality rate was 11.5%. The final population included five sets of heterozygote twins. Blood samples were obtained from 52 infants at a median time of 12 h and a maximum time of 48 h after birth. No difference was observed between the ROP and the no-ROP groups for sampling time, gender, use of erythropoietin (EPO) and cerebral hemorrhage. ROP was associated with significantly higher sepsis and red blood cell transfusion (p = 0.023 and p = 0.025, respectively), lower GA and birth weight (p < 0.001 and p = 0.001, respectively) and higher duration of mechanical ventilation (p < 0.001). The mean follow-up for ROP screening was 11.3 ± 4.5 weeks of life.

TABLE 1

TABLE 1. Main characteristics of the OmegaROP population.

The mean number of screening examinations was 3.6 ± 2.0 per infant. The incidence of ROP was 51.9%, including three cases of type 1 ROP (11.1%) and 24 cases of type 2 ROP (88.9%). Subjects with type 1 ROP underwent laser therapy on both eyes. No intravitreal injection of bevacizumab was performed. Twenty-six ROP cases were observed in zone 2 (96.3%) and one ROP case was observed in zone 3 (3.7%). We did not observe any ROP in zone 1. One case of ROP was stage 1 (3.7%), 19 cases were stage 2 (70.4%) and seven cases were stage 3 (25.9%). Four subjects were classified as a “plus” stage and three of them underwent laser treatment.

Individual phospholipid species distribution in erythrocytes

No significant difference was observed between groups in the proportions of individual EtnGpl and ChoGpl species as well as in total levels of PtdCho, PlsCho, PtdEtn, and PlsEtn (Figure 1 and Supplementary Table S1). The predominant species were PtdCho16:0/18:1, PtdCho16:0/16:0 and PtdCho16:0/20:4 + PtdCho16:1/20:3 for ChoGpl and PtdEtn16:0/18:1, PtdEtn16:0/20:4 + PtdEtn16:1/20:3, and PtdEtn18:0/20:4 + 16:0/22:4 for EtnGpl. The most represented plasmalogen species were PlsCho16:0/16:1, PlsCho18:0/16:0 and PlsCho18:1/18:2 for the choline subgroup and PlsEtn18:1/22:4 and PlsEtn18:0/22:4 the ethanolamine subgroup. The wide interquartile ranges confirmed the high interindividual variability previously observed for the fatty acid concentrations (Pallot et al., 2019).

FIGURE 1

FIGURE 1. Concentrations of major individual ChoGpl species (A) and major EtnGpl species (B) of red blood cells of preterm infants without or with retinopathy of prematurity (% of total ChoGpl or % of total EtnGpl). No significant difference was observed between no-ROP and ROP groups for any specie. Abbreviations of individual phospholipid species are as follows: Position on the glycerol backbone as shown as sn-1/sn-2 of the fatty alcohol radicals (abbreviated as number of carbons: Number of double bonds). PtdCho: Phosphatidylcholine; PtdEtn: Phosphatidylethanolamine; ROP: Retinopathy of prematurity.

Associations between gestational age and erythrocyte phospholipid species

As in our previous work on the OmegaROP cohort and considering the results of the principal component analysis showing an interaction between GA and lipid data (Pallot et al., 2019), we checked for Spearman correlations between GA and individual phospholipid species (Table 2). In order to identify the phospholipid origin of the differential accumulation of ARA and/or DHA in subjects with or without ROP, we focused our attention on phospholipids carrying n-6 and n-3 PUFAs. Even if the correlations are weak, total ChoGpl and EtnGpl carrying n-6 PUFAs were negatively associated with GA in the no-ROP group (rho = −0.485, p = 0.013 and rho = −0.477, p = 0.015 for PtdCho + PlsCho and PtdEtn + PlsEtn, respectively). Within the choline subgroup, this negative association was weak but significant only for total PtdCho species esterified with n-6 PUFAs (rho = −0.509, p = 0.009), but not for PlsCho. PtdCho species carrying ARA are likely to contribute to this finding since total PtdCho species carrying ARA were also negatively associated with GA in the no-ROP group (rho = −0.495, p = 0.011, Panel B Figure 2), and particularly the individual PtdCho16:0/20:4, PtdCho18:1/20:4, and PtdCho18:0/20:4 species (Table 3). PC carrying ARA represented 28.8% of total PtdCho species (Panel A on Figure 2). Within EtnGpl, only PlsEtn carrying n-6 PUFAs were negatively associated with GA (rho = −0.587, p = 0.002). Although the correlation is not strong, PlsEtn18:1/20:4 may be a significant contributor to this observation (rho = −0.395, p = 0.049). No significant association with GA was observed for PtdEtn species.

TABLE 2

TABLE 2. Spearman correlations between gestational age and erythrocyte ChoGpl and EtnGpl esterified with n-6 and/or n-3 PUFAs in preterm infants with or without retinopathy of prematurity.

FIGURE 2

FIGURE 2. PtdCho with ARA are associated with GA in preterm infants without ROP. (A) PtdCho with ARA (20:4) accounted for 28.8% of total retinal PtdCho species. (B) The sum of PtdCho species carrying ARA was negatively associated to GA in erythrocytes of preterm infants of the no-ROP group (rho = −0.495, p = 0.011). Abbreviations of individual PtdCho species are as follows: Position on the glycerol backbone as shown as sn-1/sn-2 of the fatty alcohol radicals (abbreviated as number of carbons: Number of double bonds). PtdCho: Phosphatidylcholine; ROP: Retinopathy of prematurity.

TABLE 3

TABLE 3. Spearman correlations between gestational age and erythrocyte individual phospholipids species esterified with n-6 and/or n-3 PUFAs in preterm infants with or without retinopathy of prematurity.

In the no-ROP group, no significant association between GA and phospholipids carrying n-6 PUFAs was observed (Tables 2 and 3). Only a weak significant negative association was observed between GA and total EtnGpl esterified with n-3 PUFAs (rho = −0.380, p = 0.034). Probably as a consequence, the ratios of total PtdEtn with n-6 PUFAs to total PtdEtn with n-3 PUFAs, total PlsEtn with n-6 PUFAs to total PlsEtn with n-3 PUFAs, and total PtdEtn with ARA to total PtdEtn with DHA were impacted and positively associated with GA in the ROP group (rho = 0.420, 0.886, and 0.843, p = 0.029, 0.006, and 0.011, respectively; Table 2).

Discussion

This study characterizes the phospholipid composition of erythrocytes in preterm infants born before 29 weeks GA. The clinical characteristics of our population were comparable to those of several studies. Indeed, the incidence of ROP was high (51.9%) and in agreement with other very-low-GA populations (Austeng, Kallen, Ewald, Jakobsson, & Holmstrom, 2009).

Our data show a weak but negative association between ChoGpl and EtnGpl carrying n-6 PUFAs with GA in the no-ROP group, while no significant association was observed in the ROP group. These findings are in the line with our previous reports related to the OmegaROP study (Pallot et al., 2019). Considering that ChoGpl and EtnGpl species represent more than 90% of total phospholipids in erythrocytes (Acar et al., 2007; Acar et al., 2009; Acar et al., 2012), we assumed that the changes observed in their concentrations would be a reliable indicator of the whole phospholipid pool in erythrocytes. However, further analyses on phosphatidylserines and phosphatidylinositols could be of interest to draw a more complete picture of erythrocyte phospholipidome alterations in ROP.

In preterm infants who will not develop ROP, erythrocyte relative levels of phospholipids carrying ARA decrease as GA increases, while no change was observed in preterm infants who will develop ROP. This negative correlation mostly relies on PtdCho species and more specifically on individual PtdCho16:0/20:4, PtdCho18:1/20:4, and PtdCho18:0/20:4 species. This finding suggests an in utero accumulation of ARA in erythrocytes of preterm infants that will develop ROP relatively to those that will not develop ROP. This GA-related change is in line with the finding of Bernhard and collaborators that revealed a remodeling in plasma PtdCho species esterified with ARA in early preterm infants (Bernhard et al., 2016). They showed that the plasma PtdCho ARA to DHA ratio in preterm infants of less than 33 weeks GA largely exceeds that of term infants, thus likely contributing to the prematurity-related impaired overall health. The data we have obtained in 25- to 28-weeks GA newborns is then in agreement with this observation. However, it is important to keep in mind that the biological materials used in these two studies were different, making that further analyses on plasma PtdCho of our subjects would be of interest.

Whether the relative accumulation of ARA in erythrocytes is associated to the onset of ROP remains to be further investigated. Indeed, on one hand it is very clear that together with DHA, ARA is essential during gestation and early postnatal life for an optimal development of the infant (Tai, Wang, & Chen, 2013; Koletzko et al., 2020; Bernhard, Poets, & Franz, 2021). On the other hand, ARA is known to play a pivotal role in the promotion of inflammation, especially through its eicosanoid derivatives such as prostaglandin E2 (PGE2). PGE2 is known to be involved in abnormal angiogenesis and in the pathogenesis of proliferative retinopathies such as ROP (Yanni, Barnett, Clark, & Penn, 2009; Hartnett & Penn, 2012; Schoenberger et al., 2012; Xie et al., 2021). Accordingly, other studies have shown that inflammation is a significant risk factor for developing ROP (Holm et al., 2017; Rivera et al., 2017). Finally, it is recognized that ARA competes with DHA, the latter being known to inhibit retinal neovascularization in a mouse model of ROP (Connor et al., 2007). DHA and ARA are the most prevalent PUFAs in the human retina (Bretillon et al., 2008; Acar et al., 2012). As a subtle equilibrium in their needs may exist within this tissue, we cannot exclude that minor variations in the ARA to DHA ratio may impact retinal physiology.

In accordance with our previous observation on the same cohort (Pallot et al., 2019), our data display a high interindividual variability. This is especially true for ARA whose levels in total lipids were of 13.48 ± 7.29% and 13.86 ± 7.63% of total fatty acids in the no-ROP and ROP groups, respectively (Pallot et al., 2019). The present study shows that ARA levels were low especially in PtdCho, and independently from the development of ROP. Although extensive analyses, we were unable to connect these very low ARA levels to any other variable, including GA, blood sampling time after birth or any other parameter related to the infants’ health. Since we have used analytical QCs and as we have repeated our analyses several times, we believe that these changes are of physiologic origin and not the result of analytical bias. The probable physiological origin of these changes is reinforced by similar observations showing very-low levels of ARA (Glen et al., 1994; Khan et al., 2002) or even ARA levels close to zero (Li et al., 2022) in erythrocytes of human subjects with neurologic disorders or of infants with severe malnutrition (Smit et al., 1997). Such data may suggest a marked dysregulation of erythrocyte membrane fatty acids in some infants. It can be the result of either increased phospholipase A2 activity and/or enhanced lipid peroxidation. Free radical oxidation during technical steps of sample handling and lipid analysis can also be at the origin of the degradation of PUFAs such as AA, but also that of plasmalogens. To avoid such biases in our study, we took care to take a number of precautions starting from immediate processing of the blood sample after venipuncture, to the isolation and washing of red cells at a temperature of 4°C, their immediate storage at −80°C, and the use of the Moilanen & Nikkari procedure for lipid extraction. This last methodology doesn’t use acidic conditions that are known to be deleterious for plasmalogens. Another hypothesis to explain the low levels of ARA would rely on a reduced bioavailability of ARA in maternal blood and/or its abnormal placental transfer. Further investigation is needed to better understand the origin of these modifications.

Circulating lipids are also subject to a postnatal remodeling. Indeed, several studies documented an increase in the plasma and erythrocyte levels of linoleic acid (the dietary precursor of ARA) and a decrease in the levels of ARA during the first weeks of life as a consequence of a LA-rich nutrition (Bernhard et al., 2014; Bockmann et al., 2021). Considering that the birth to blood sampling delay in our study exceeded 24 h in some cases, we cannot exclude the influence of the source of lipid (parenteral versus enteral nutrition) on the lipid data, especially in infants with lower GA for which the parenteral lipid supply is higher. Nevertheless, we did not find any association between erythrocyte lipids (and especially PtdCho16:0/20:4) and the sampling time in the subjects of our cohort (data not shown).

Interestingly, the n-6 to n-3 PUFA ratio in erythrocytes seems to display a similar pattern in the no-ROP group when considering PlsEtn species, as PlsEtn carrying n-6 PUFAs were negatively associated with GA (rho = −0.587, p = 0.002), with PlsEtn18:1/20:4 being the most significant contributor to this observation (rho = −0.395, p = 0.049). Considering that erythrocyte lipid composition could represent a reliable indicator of the lipid composition of the retina in newborns (Carlson, Carver, & House, 1986; Makrides, Neumann, Byard, Simmer, & Gibson, 1994; Sarkadi-Nagy et al., 2004), we may speculate that retinal lipids display similar modifications in ARA, DHA and plasmalogen levels of newborns developing ROP. On the contrary, infants that will develop ROP seem to accumulate n-6 PUFAs in their plasmalogens as the ratio of PlsEtn carrying n-6 PUFAs to PlsEtn carrying n-3 PUFAs is strongly positively associated to GA (rho = 0.886, p = 0.006), ARA and DHA being the major contributors to this observation (rho = 0.843, p = 0.011). These data are in line with the ROP-like phenotype observed in plasmalogen-deficient mice (Saab et al., 2014) and then reinforce the idea that plasmalogens may participate to the pathophysiology of this disease. It would be however interesting to check whether these mice display changes in the levels of PUFAs such as ARA and/or DHA in their erythrocytes in order to strengthen the conclusions of the present paper. While PUFAs status of the newborn correlates with maternal status (Bernhard et al., 2016), and as it can be influenced by dietary supplementation after birth (Schulzke, Patole, & Simmer, 2011), plasmalogen content of tissues only relies on fetal de novo synthesis as no materno-fetal transfer of ether-lipids has been demonstrated so far (Das, Holmes, Wilson, & Hajra, 1992). In our study, blood samples were collected within the first 48 h of life. Considering that the erythrocyte life-span ranges from 35 to 50 days in preterm infants (Pearson, 1967), the modifications observed in the present study might be assigned to differences in situ lipid metabolisms even if, again, a nutritional influence cannot be excluded for PUFAs.

Several limitations must be acknowledged in our study. First, the correlation coefficients found in this study are weak and deserve more investigations. Second, our population included a limited number of patients and subjects. Third, our population only included three preterm infants with severe ROP, while this population was specifically concerned by treatments targeting vascular events. Furthermore, while ChoGpl and EtnGpl species are the main contributors to erythrocytes phospholipidome, further lipidomic analyses including phosphatidylserine and phosphatidylinositol individual phospholipid species could provide useful complementary information.

Taken together, our results seem to confirm the alterations of erythrocyte PUFA profile in preterm infants developing ROP and suggest that plasmalogens may contribute to them. Considering the importance of plasmalogens in the cellular bioavailability of PUFAs and their involvement in the vascular development of the retina, our study suggests that investigating the relationships between plasmalogen metabolism and ROP could be of particular interest to decipher the pathophysiological mechanisms driving ROP.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ebi.ac.uk/metabolights/MTBLS4219/descriptors

Ethics statement

The studies involving human participants were reviewed and approved by CPP Est III, School of Medicine, Dijon, France. Written informed consent from the participants’ legal guardian/next of kin was not required to participate in this study in accordance with the national legislation and the institutional requirements.

Author contributions

CP, DS, CC-G, LB, NA designed the study. RK, CP, SC, JM, DS, CC, NT, DM, AD, SG, OB, NA acquired, analyzed and interpreted the data. RK, CP, JM, DS, AB, CC-G, LB, OB, NA drafted the manuscript. All authors reviewed the manuscript.

Funding

The authors thank National Research Institute for Agriculture, Food and Environment (INRAE); Regional Council of Burgundy (PARI Grant); European Regional Development Fund (FEDER); Agence Nationale de la Recherche (ANR-11-LABX-0021-01); VISIO Foundation; Fondation de France/Fondation de l’Oeil; Groupe Lipides et Nutrition (GLN) for their financial support.

Acknowledgments

The authors would like to thank all medical staff of the Neonatal Intensive Care Unit of Dijon University Hospital for their unique contribution to the study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

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/fcell.2022.921691/full#supplementary-material

References

Acar, N., Berdeaux, O., Gregoire, S., Cabaret, S., Martine, L., Gain, P., et al. (2012). Lipid composition of the human eye: Are red blood cells a good mirror of retinal and optic nerve fatty acids? PLoS ONE 7 (4), e35102. doi:10.1371/journal.pone.0035102

PubMed Abstract | CrossRef Full Text | Google Scholar

Acar, N., Berdeaux, O., Juaneda, P., Gregoire, S., Cabaret, S., Joffre, C., et al. (2009). Red blood cell plasmalogens and docosahexaenoic acid are independently reduced in primary open-angle glaucoma. Exp. Eye Res. 89 (6), 840–853. doi:10.1016/j.exer.2009.07.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Acar, N., Gregoire, S., Andre, A., Juaneda, P., Joffre, C., Bron, A. M., et al. (2007). Plasmalogens in the retina: In situ hybridization of dihydroxyacetone phosphate acyltransferase (DHAP-AT) the first enzyme involved in their biosynthesis and comparative study of retinal and retinal pigment epithelial lipid composition. Exp. Eye Res. 84 (1), 143–151. doi:10.1016/j.exer.2006.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Austeng, D., Kallen, K. B., Ewald, U. W., Jakobsson, P. G., and Holmstrom, G. E. (2009). Incidence of retinopathy of prematurity in infants born before 27 weeks' gestation in Sweden. Arch. Ophthalmol. 127 (10), 1315–1319. doi:10.1001/archophthalmol.2009.244

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartlett, E. M., and Lewis, D. H. (1970). Spectrophotometric determination of phosphate esters in the presence and absence of orthophosphate. Anal. Biochem. 36 (1), 159–167. doi:10.1016/0003-2697(70)90343-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Berdeaux, O., Juaneda, P., Martine, L., Cabaret, S., Bretillon, L., Acar, N., et al. (2010). Identification and quantification of phosphatidylcholines containing very-long-chain polyunsaturated fatty acid in bovine and human retina using liquid chromatography/tandem mass spectrometry. J. Chromatogr. A 1217 (49), 7738–7748. doi:10.1016/j.chroma.2010.10.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernhard, W., Poets, C. F., and Franz, A. (2021). Parenteral nutrition for preterm infants: Correcting for arachidonic and docosahexaenoic acid may not suffice. Arch. Dis. Child. Fetal Neonatal Ed. 106 (6), 683. doi:10.1136/archdischild-2021-321871

CrossRef Full Text | Google Scholar

Bernhard, W., Raith, M., Koch, V., Kunze, R., Maas, C., Abele, H., et al. (2014). Plasma phospholipids indicate impaired fatty acid homeostasis in preterm infants. Eur. J. Nutr. 53 (7), 1533–1547. doi:10.1007/s00394-014-0658-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernhard, W., Raith, M., Koch, V., Maas, C., Abele, H., Poets, C. F., et al. (2016). Developmental changes in polyunsaturated fetal plasma phospholipids and feto-maternal plasma phospholipid ratios and their association with bronchopulmonary dysplasia. Eur. J. Nutr. 55 (7), 2265–2274. doi:10.1007/s00394-015-1036-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Bockmann, K. A., von Stumpff, A., Bernhard, W., Shunova, A., Minarski, M., Frische, B., et al. (2021). Fatty acid composition of adipose tissue at term indicates deficiency of arachidonic and docosahexaenoic acid and excessive linoleic acid supply in preterm infants. Eur. J. Nutr. 60 (2), 861–872. doi:10.1007/s00394-020-02293-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Bretillon, L., Thuret, G., Gregoire, S., Acar, N., Joffre, C., Bron, A. M., et al. (2008). Lipid and fatty acid profile of the retina, retinal pigment epithelium/choroid, and the lacrimal gland, and associations with adipose tissue fatty acids in human subjects. Exp. Eye Res. 87 (6), 521–528. doi:10.1016/j.exer.2008.08.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Carlson, S. E., Carver, J. D., and House, S. G. (1986). High fat diets varying in ratios of polyunsaturated to saturated fatty acid and linoleic to linolenic acid: A comparison of rat neural and red cell membrane phospholipids. J. Nutr. 116 (5), 718–725. doi:10.1093/jn/116.5.718

PubMed Abstract | CrossRef Full Text | Google Scholar

Connor, K. M., SanGiovanni, J. P., Lofqvist, C., Aderman, C. M., Chen, J., Higuchi, A., et al. (2007). Increased dietary intake of omega-3-polyunsaturated fatty acids reduces pathological retinal angiogenesis. Nat. Med. 13 (7), 868–873. doi:10.1038/nm1591

PubMed Abstract | CrossRef Full Text | Google Scholar

Das, A. K., Holmes, R. D., Wilson, G. N., and Hajra, A. K. (1992). Dietary ether lipid incorporation into tissue plasmalogens of humans and rodents. Lipids 27 (6), 401–405. doi:10.1007/BF02536379

PubMed Abstract | CrossRef Full Text | Google Scholar

European Union (2008). Ethical considerations for clinical trials on medicinal products conducted with the paediatric population. Eur. J. Health Law 15 (2), 223–250.

PubMed Abstract | CrossRef Full Text | Google Scholar

Glen, A. I., Glen, E. M., Horrobin, D. F., Vaddadi, K. S., Spellman, M., Morse-Fisher, N., et al. (1994). A red cell membrane abnormality in a subgroup of schizophrenic patients: Evidence for two diseases. Schizophr. Res. 12 (1), 53–61. doi:10.1016/0920-9964(94)90084-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Good, W. V., and Hardy, R. J. (2001). The multicenter study of early treatment for retinopathy of prematurity (ETROP). Ophthalmology 108 (6), 1013–1014. doi:10.1016/s0161-6420(01)00540-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Hartnett, M. E., and Penn, J. S. (2012). Mechanisms and management of retinopathy of prematurity. N. Engl. J. Med. 367 (26), 2515–2526. doi:10.1056/NEJMra1208129

PubMed Abstract | CrossRef Full Text | Google Scholar

Hellstrom, A., Smith, L. E., and Dammann, O. (2013). Retinopathy of prematurity. Lancet 382 (9902), 1445–1457. doi:10.1016/S0140-6736(13)60178-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Holm, M., Morken, T. S., Fichorova, R. N., VanderVeen, D. K., Allred, E. N., Dammann, O., et al. (2017). Systemic inflammation-associated proteins and retinopathy of prematurity in infants born before the 28th week of gestation. Invest. Ophthalmol. Vis. Sci. 58 (14), 6419–6428. doi:10.1167/iovs.17-21931

PubMed Abstract | CrossRef Full Text | Google Scholar

International Committee for the Classification of Retinopathy of, P (2005). The international classification of retinopathy of prematurity revisited. Arch. Ophthalmol. 123 (7), 991–999. doi:10.1001/archopht.123.7.991

PubMed Abstract | CrossRef Full Text | Google Scholar

Jordan, C. O. (2014). Retinopathy of prematurity. Pediatr. Clin. North Am. 61 (3), 567–577. doi:10.1016/j.pcl.2014.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Khan, M. M., Evans, D. R., Gunna, V., Scheffer, R. E., Parikh, V. V., Mahadik, S. P., et al. (2002). Reduced erythrocyte membrane essential fatty acids and increased lipid peroxides in schizophrenia at the never-medicated first-episode of psychosis and after years of treatment with antipsychotics. Schizophr. Res. 58 (1), 1–10. doi:10.1016/s0920-9964(01)00334-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Koletzko, B., Bergmann, K., Brenna, J. T., Calder, P. C., Campoy, C., Clandinin, M. T., et al. (2020). Should formula for infants provide arachidonic acid along with DHA? A position paper of the European academy of paediatrics and the child health foundation. Am. J. Clin. Nutr. 111 (1), 10–16. doi:10.1093/ajcn/nqz252

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, N., Yang, P., Tang, M., Liu, Y., Guo, W., Lang, B., et al. (2022). Reduced erythrocyte membrane polyunsaturated fatty acid levels indicate diminished treatment response in patients with multi- versus first-episode schizophrenia. NPJ Schizophr. 8 (1), 7. doi:10.1038/s41537-022-00214-2

CrossRef Full Text | Google Scholar

Lofqvist, C. A., Najm, S., Hellgren, G., Engstrom, E., Savman, K., Nilsson, A. K., et al. (2018). Association of retinopathy of prematurity with low levels of arachidonic acid: A secondary analysis of a randomized clinical trial. JAMA Ophthalmol. 136 (3), 271–277. doi:10.1001/jamaophthalmol.2017.6658

PubMed Abstract | CrossRef Full Text | Google Scholar

Makrides, M., Neumann, M. A., Byard, R. W., Simmer, K., and Gibson, R. A. (1994). Fatty acid composition of brain, retina, and erythrocytes in breast- and formula-fed infants. Am. J. Clin. Nutr. 60 (2), 189–194. doi:10.1093/ajcn/60.2.189

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, C. R., Dasilva, D. A., Cluette-Brown, J. E., Dimonda, C., Hamill, A., Bhutta, A. Q., et al. (2011). Decreased postnatal docosahexaenoic and arachidonic acid blood levels in premature infants are associated with neonatal morbidities. J. Pediatr. 159 (5), 743–742. doi:10.1016/j.jpeds.2011.04.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Moilanen, T., and Nikkari, T. (1981). The effect of storage on the fatty acid composition of human serum. Clin. Chim. Acta. 114 (1), 111–116. doi:10.1016/0009-8981(81)90235-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Nagan, N., and Zoeller, R. A. (2001). Plasmalogens: Biosynthesis and functions. Prog. Lipid Res. 40 (3), 199–229. doi:10.1016/s0163-7827(01)00003-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Pallot, C., Mazzocco, J., Meillon, C., Semama, D. S., Chantegret, C., Ternoy, N., et al. (2019). Alteration of erythrocyte membrane polyunsaturated fatty acids in preterm newborns with retinopathy of prematurity. Sci. Rep. 9 (1), 7930. doi:10.1038/s41598-019-44476-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Pearson, H. A. (1967). Life-span of the fetal red blood cell. J. Pediatr. 70 (2), 166–171. doi:10.1016/s0022-3476(67)80410-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Rivera, J. C., Holm, M., Austeng, D., Morken, T. S., Zhou, T. E., Beaudry-Richard, A., et al. (2017). Retinopathy of prematurity: Inflammation, choroidal degeneration, and novel promising therapeutic strategies. J. Neuroinflammation 14 (1), 165. doi:10.1186/s12974-017-0943-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Saab, S., Buteau, B., Leclere, L., Bron, A. M., Creuzot-Garcher, C. P., Bretillon, L., et al. (2014). Involvement of plasmalogens in post-natal retinal vascular development. PLoS ONE 9 (6), e101076. doi:10.1371/journal.pone.0101076

PubMed Abstract | CrossRef Full Text | Google Scholar

Saab-Aoude, S., Bron, A. M., Creuzot-Garcher, C. P., Bretillon, L., and Acar, N. (2013). A mouse model of in vivo chemical inhibition of retinal calcium-independent phospholipase A2 (iPLA2). Biochimie 95 (4), 903–911. doi:10.1016/j.biochi.2012.12.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Sapieha, P., Stahl, A., Chen, J., Seaward, M. R., Willett, K. L., Krah, N. M., et al. (2011). 5-Lipoxygenase metabolite 4-HDHA is a mediator of the antiangiogenic effect of omega-3 polyunsaturated fatty acids. Sci. Transl. Med. 3 (69), 69ra12. doi:10.1126/scitranslmed.3001571

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarkadi-Nagy, E., Wijendran, V., Diau, G. Y., Chao, A. C., Hsieh, A. T., Turpeinen, A., et al. (2004). Formula feeding potentiates docosahexaenoic and arachidonic acid biosynthesis in term and preterm baboon neonates. J. Lipid Res. 45 (1), 71–80. doi:10.1194/jlr.M300106-JLR200

PubMed Abstract | CrossRef Full Text | Google Scholar

Sato, T., Shima, C., and Kusaka, S. (2011). Vitreous levels of angiopoietin-1 and angiopoietin-2 in eyes with retinopathy of prematurity. Am. J. Ophthalmol. 151 (2), 353–357. doi:10.1016/j.ajo.2010.08.037

PubMed Abstract | CrossRef Full Text | Google Scholar

Schoenberger, S. D., Kim, S. J., Sheng, J., Rezaei, K. A., Lalezary, M., Cherney, E., et al. (2012). Increased prostaglandin E2 (PGE2) levels in proliferative diabetic retinopathy, and correlation with VEGF and inflammatory cytokines. Invest. Ophthalmol. Vis. Sci. 53 (9), 5906–5911. doi:10.1167/iovs.12-10410

PubMed Abstract | CrossRef Full Text | Google Scholar

Schulzke, S. M., Patole, S. K., and Simmer, K. (2011). Long-chain polyunsaturated fatty acid supplementation in preterm infants. Cochrane Database Syst. Rev. 2, CD000375. doi:10.1002/14651858.CD000375.pub4

PubMed Abstract | CrossRef Full Text | Google Scholar

Smit, E. N., Dijkstra, J. M., Schnater, T. A., Seerat, E., Muskiet, F. A., Boersma, E. R., et al. (1997). Effects of malnutrition on the erythrocyte fatty acid composition and plasma vitamin E levels of Pakistani children. Acta Paediatr. 86 (7), 690–695. doi:10.1111/j.1651-2227.1997.tb08569.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sonmez, K., Drenser, K. A., Capone, A., and Trese, M. T. (2008). Vitreous levels of stromal cell-derived factor 1 and vascular endothelial growth factor in patients with retinopathy of prematurity. Ophthalmology 115 (6), 1065–1070. doi:10.1016/j.ophtha.2007.08.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Tai, E. K., Wang, X. B., and Chen, Z. Y. (2013). An update on adding docosahexaenoic acid (DHA) and arachidonic acid (AA) to baby formula. Food Funct. 4 (12), 1767–1775. doi:10.1039/c3fo60298b

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, T., Zhang, Z., Cui, Y., Shu, Y., Liu, Y., Zou, J., et al. (2021). Prostaglandin E2 promotes pathological retinal neovascularisation via EP4R-EGFR-Gab1-AKT signaling pathway. Exp. Eye Res. 205, 108507. doi:10.1016/j.exer.2021.108507

PubMed Abstract | CrossRef Full Text | Google Scholar

Yanni, S. E., Barnett, J. M., Clark, M. L., and Penn, J. S. (2009). The role of PGE2 receptor EP4 in pathologic ocular angiogenesis. Invest. Ophthalmol. Vis. Sci. 50 (11), 5479–5486. doi:10.1167/iovs.09-3652

PubMed Abstract | CrossRef Full Text | Google Scholar