Next Article in Journal
Cytokinin Response Factor 9 Represses Cytokinin Responses in Flower Development
Next Article in Special Issue
Sphingosine-1-phosphate Decreases Erythrocyte Dysfunction Induced by β-Amyloid
Previous Article in Journal
Transcriptional Regulatory Network of Plant Cadmium Stress Response
Previous Article in Special Issue
Sublingual Microcirculation Specificity of Sickle Cell Patients: Morphology of the Microvascular Bed, Blood Rheology, and Local Hemodynamics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Untargeted Lipidomics of Erythrocytes under Simulated Microgravity Conditions

1
Department of Life and Environmental Sciences, Cittadella Universitaria di Monserrato, Blocco A, Room 13, 09042 Monserrato, Italy
2
Department of Mechanical, Chemical and Materials Engineering, University of Cagliari, Via Marengo 2, 09123 Cagliari, Italy
3
Department of Biomedical Science, University of Sassari, Viale San Pietro, 07100 Sassari, Italy
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(5), 4379; https://doi.org/10.3390/ijms24054379
Submission received: 18 January 2023 / Revised: 13 February 2023 / Accepted: 20 February 2023 / Published: 22 February 2023
(This article belongs to the Special Issue Roles of Erythrocytes in Human Health and Disease)

Abstract

:
Lipidomics and metabolomics are nowadays widely used to provide promising insights into the pathophysiology of cellular stress disorders. Our study expands, with the use of a hyphenated ion mobility mass spectrometric platform, the understanding of the cellular processes and stress due to microgravity. By lipid profiling of human erythrocytes, we annotated complex lipids such as oxidized phosphocholines, phosphocholines bearing arachidonic in their moiety, as well as sphingomyelins and hexosyl ceramides associated with microgravity conditions. Overall, our findings give an insight into the molecular alterations and identify erythrocyte lipidomics signatures associated with microgravity conditions. If the present results are confirmed in future studies, they may help to develop suitable treatments for astronauts after return to Earth.

1. Introduction

The term microgravity in general refers to existing residual accelerations. When gravitation is the only force acting on an object, then it results is in free fall and hence it will experience microgravity [1]. Weightlessness is the state in which a body having a certain weight is balanced by another force or remains in free fall without feeling the effects of the atmosphere, equivalent to the situation faced by an astronaut aboard a spaceship. The effects of microgravity on human physiology have been studied extensively since the time of Yuri Gagarin (in 1961) who experienced the first man-on-board orbital flight, revealing profound implications for human health [2]. Despite the great interest and commitment of the scientific community, the mechanisms by which microgravity exerts its effects on the human body are not entirely clear.
Acute changes in normal physiology are typically seen in astronauts as a response and adaptation to abnormal environments. Such peculiar alterations require the attention of doctors and scientists [3]. In addition to alterations at the genetic level [4], microgravity experienced by space travellers also induces profound alterations at the cellular level. These alterations occurring at the cellular level are reflected in a series of pathological conditions such as a reduction of bone density, muscle atrophy, endocrine disorders, cognitive disorders and cardiovascular disfunctions, body fluid and electrolyte reduction, motion sickness, immune inhibition, and anaemia [5]. For these reasons, ground-based experiments simulating factors of spaceflight conditions are needed.
Microgravity is studied in several scientific and technological fields with the aim to highlight processes that on Earth are masked by the effects of the high gravitational field. Furthermore, the study of physiological processes in microgravity conditions allows the identification of the molecular mechanisms involved in different pathologies [6].
Roughly 350 people have experienced spaceflight in the past four decades, making it difficult to develop higher levels of clinical evidence to evaluate the effectiveness of space medicine interventions [7]. This great limitation, and the importance of studying the alterations at the cellular level that affect astronauts has led various research groups to study and build instruments capable of simulating space gravitational conditions on Earth.
The most employed methods to simulate microgravity are random positioning machines (RPM) and clinostats [8]. By controlled simultaneous rotating of the two axes, the clinostat cancels the cumulative gravity vector at the centre of the device, producing an environment with an average of 10−3 g. This is accomplished by the rotation of a chamber at the centre of the device to disperse the gravity vector uniformly within a spherical volume at a constant angular speed [9].
Bioactive lipid molecules known as signalling molecules, such as fatty acid, eicosanoids, diacylglycerol, phosphatidic acid, lysophosphatidic acid, ceramide, sphingosine, sphingosine-1-phosphate, phosphatidylinositol-3 phosphate, and cholesterol, are involved in the activation or regulation of different signalling pathways leading to apoptosis. Furthermore, alterations in the lipid composition determine membrane rigidity and fluidity, and play a crucial role in membrane organization, dynamics, and function [10]. Because of their biological role, lipids have been the subject of an intense area of research since the 1960s, which unfortunately was held back due to limited instrument platforms. Nowadays, lipidomics is considered an emerging science of fundamental importance for clarifying the biochemical pathways involved in several pathologies or cellular stress adaptations [11]. Advances in mass spectrometry (MS) and data processing, as well as the incorporation of soft ionization techniques, as ESI-MS2 method, has revolutionized the use of mass spectrometry, ushering this analytical tool in the field of lipidomics [12]. The lipidomics study can be applied as untargeted and targeted approaches, each with its own advantages and limitations [13]. Untargeted lipidomics focuses on the analysis of all detectable metabolites in a sample, including unknown chemicals, while targeted lipidomics is the measurement of defined groups of metabolites. While the strength of the targeted approach is validating one or more hypotheses, untargeted lipidomics allows for the discovery of new compounds that have led to a number of breakthroughs in understanding human disease risks [14]. Untargeted analyses can be performed with or without the addition of internal rules. When internal standards are added to samples, the method can provide pseudoconcentration results for particular metabolites or for metabolites with similar physicochemical properties (e.g., lipids). While these results are not truly quantitative, they may be accurate enough for case/control comparisons [15].
Despite the great relevance of the topic, currently few studies have been carried out to investigate the behaviour of lipids in erythrocyte samples cultured under simulated gravity conditions. In 2009, Ivanova et al. investigated blood samples from Russian cosmonauts by observing significant changes in the phospholipids class [16]. An increase in the percentage of phosphatidylcholine may be clearly associated with the increase in membrane rigidity. On the other hand, changes in the physicochemical properties of the plasma membrane of erythrocytes (microviscosity and permeability) can influence the efficiency of oxygen transfer, the state of the haemoglobin, and changes in the conformation of hematoporphyrin. Furthermore, changes in the in erythrocyte structure through an ultrastructural morphological analysis can be assessed by atomic force microscopy [17]. However, the study conducted by Ivanova’s team reported data deriving mainly from studies carried out after the end of a space flight, while only few data are related to changes that occur during a space flight. Moreover, lipid and phospholipid compositions of erythrocyte membranes were assayed by thin layer chromatography followed by densitometric measurement of stained dots. This technique provides information on the entire lipid class, but hardly allows the recognition of the specific lipid compounds. Since no data are reported on this subject, we decided to exploit the potential offered by chromatographic and mass spectrometry innovations to better understand the lipid modifications suffered by erythrocytes during simulated microgravity conditions. For this reason, with the aim to better understand which metabolic and/or structural changes occur in the erythrocytes subjected to low gravity, an experimental analysis of the erythrocytes’ lipid profile and their morphology under normal- and micro-g conditions was carried out following a recent investigation on the subject [18]. In detail, human erythrocytes were cultured in simulated gravity conditions, and they were collected at different times of clinorotation. For each sample, the organic phase was collected and analysed through ion mobility Q-TOF mass spectrometer (UHPLC-IM-QTOF-MS).

2. Results and Discussion

To investigate the erythrocytes’ lipid profile after clinorotation and to describe possible variations among the different lipid categories, samples were analysed by IM-QTOF-LC/MS and representative total ion chromatograms are shown in Figure 1.
Data processing yielded 215 and 160 features for the positive (PIA) and negative ionization analysis (NIA), respectively, which were subjected to multivariate statistical analysis (MVA). Chemical composition analysis indicated that the lipid fraction was composed of lipids from the following classes: free fatty acid (FA), lysophosphatidylcholines (LysoPC), phosphatidylcholines (PC), phosphatidylethanolamines (PE), sphingomyelins (SM), ceramides (Cer), and ether-linked oxidized phosphatidylcholine (EtherOxPC). Initially, to study sample distribution, to detect outliers, and to highlight differences or common features, a PCA was performed. The unsupervised analysis of both PIA and NIA features did not indicate any sample clustering correlated to clinorotation as shown in Figure 2.
However, the arrangement of the samples in the multivariate space appeared to be influenced by the time factor. Thus, to further limit the time factor influence, for each time point, we performed a PLS-DA. The validation parameters of the PIA and NIA models built for the samples collected at 6, 9, and 24 h are reported in the caption of the resulting plots (Figure 3).
To identify metabolites that can discriminate for the two classes of samples (clinorotated vs. control samples), an OPLS-DA model of the IM-QTOF-LC/MS data was performed for each time point and for both polarities of acquisition. The OPLS-DA score plots are reported in Figure 4. In Table 1, we reported the discriminant metabolites between two classes and selected based on VIP value.
Using MS/MS fragmentation data and consulting the Metlin and Lipidomics libraries, we were able to tentatively identify the most discriminant metabolites as reported in Table 1.
Astronauts, after their return from space missions, manifest significant haematological alterations. Since the earliest space missions, symptoms such as structural alterations of red blood cells [18], anaemia [19], thrombocytopenia, [20,21], 10–17% reduction in plasma volume, and haemolysis [22] were reported. For these reasons, concern about the effects of space flight on haematological processes has been increasing. Several scientific studies allowed different theories to be proposed that may explain the alterations in the size and number of erythrocytes [23,24]. Recently, Trudel et al. showed in astronauts a degradation and 54% reduction in red blood cells [25]. Different factors can lead a human cell to programmed death, such as changes to lipid signal activity [26].
Human cells determine the characteristics of the plasma membrane bilayer by tightly controlling lipid composition and recruiting cytosolic proteins involved in structural functions or signal transduction [27]. The cell membrane is a lipid bilayer essentially formed by phospholipids, cholesterol, and glycolipids [28]. Small variations in percentage composition and molar ratio of the different classes of phospholipids and glycolipids might induce changes in the cell membrane’s fluidity and permeability. In particular, phospholipids are the main components of cell membranes and perform important biological functions.
From our results, it appears that after 6 h of clinorotation, levels of phosphocholines were increased in human erythrocytes. In particular, PC 18:1_20:4, PC 18:0_20:4, PC18:1_18:1, and PC 18:2_18:1 were found to be upregulated. Notably, PC with the arachidonic acid in their moiety were found discriminants. In particular, the proportion of sn-2-arachidonoyl-phosphatidylcholine (20:4-PC) has been shown to be inversely correlated with the activity of protein kinase B (Akt), an important kinase which promotes cell proliferation and survival. 20:4-PC reduces cell proliferation by interfering with the S-phase cell transition and by suppressing Akt downstream signalling and the expression of cyclin, such as LY294002, which is a specific inhibitor of the phosphatidylinositol-3-kinase/Akt [29]. At 9 and 24 h, erythrocytes showed other 20:4-PC upregulated: PC 18:2_20:4, PC 18:3_20:4 and PC 16:0_20:4, and PC 15:0_20:4, 15:1_20:4 and 16:0_20:4, respectively.
With the classical techniques of liquid chromatography coupled to mass spectrometry, the annotation of lipids and thus phosphocholine fatty acid composition with a good confidence interval is difficult due to the large variety of lipid species with different regiochemistry. In our study, the use of an analytical platform such as ion mobility coupled to mass spectrometry providing the collision cross section (CCS) value allows a better and more confident annotation of each metabolite. Each CCS was compared with an internal database and against the unified collision cross section compendium available on LipidMaps [30].
Additionally, a different fatty acid composition of membrane components can result in a greater sensitivity to peroxidative stress, with a consequent increase in membrane fragility. Phosphatidylcholine species containing polyunsaturated fatty acids in their moiety, particularly arachidonate, at the sn-2 position are susceptible to free radical oxidation [31]. An example is represented by 1-palmitoyl-2-arachidonoyl-sn-glycero-3-phosphatidylcholine (PC16.0_20:4), which is a common cell membrane constituent, and circulates within cholesterol particles. At 6 h of clinorotation, erythrocytes showed an upregulation of EtherOxPC 16:0_20:4, while the respective phosphocholine, PC 16:0_20:4, was found to be not discriminant. Simulated microgravity conditions increase reactive oxygen species (ROS) production in various cell types [32]. Generally, in microgravity conditions a different management of cellular resources was observed. In fact, in G0 conditions, there is a more rapid consumption of intracellular ATP, and an increase in ATP expulsion compared to cells cultured under terrestrial gravity conditions, coupled with a reducing power [8] resulting in a more oxidant environment.
Furthermore, inflammation and oxidative stress are associated with lipid peroxidation and the formation of bioactive lipids such as oxidized phosphocholines [33]. C-reactive protein (CRP), an acute-phase protein of hepatic origin that binds to specific structures expressed on the surface of dead or dying cells, promotes phagocytosis as macrophages may bind to these PC-oxidized species. Furthermore, recent studies demonstrate an enrichment of oxidized phosphatidylcholine in apoptotic cells [34]. Indeed, CRP can selectively bind on oxidized phosphatidylcholine but not on native phosphatidylcholine. In addition, oxidized phospholipids are recognized by macrophage scavenger, implying that these innate immune responses participate in cell clearance due to their proinflammatory properties [35]. Moreover, oxidized phosphatidylcholine, specifically oxidized-1-palmitoyl-2-arachidonoyl-sn-glycero-3-phosphatidylcholine (EtherOxPC 16:0_20:4), seems to be involved in ROS production. According to the study of Rouhanizadeh et al., EtherOxPC 16:0_20:4 was able to induce vascular endothelial superoxide production [36].
On the contrary, after 9 h of clinorotation, EtherOxPC(16:0_20:4) was downregulated, while PC 16:0_20:4 was upregulated, resulting non-discriminant after 24 h. These findings can lead us to hypothesize the complex adaptive response of cells.
Interestingly, several sphingomyelins were found to be downregulated for each experimental time point. This should not be surprising considering the mechanism of sphingomyelin synthesis. Indeed, sphingomyelinases (SMases) catalyse the hydrolysis of sphingomyelin to form ceramide and phosphocholine [37].
Taken together, these findings indicate that there are probably several mechanisms underlying spatial anaemia: inhibition of 20:4 PC-mediated cell proliferation and a simultaneous increase in pro-apoptotic signals.

3. Materials and Methods

3.1. Chemicals

Analytical LC-grade methanol, chloroform, acetonitrile, 2-propanol, and ammonium acetate and formiate were purchased from Sigma Aldrich (Milan, Italy). Bi-distilled water was obtained with a MilliQ purification system (Millipore, Milan, Italy). A SPLASH® LIPIDOMIX® standard component mixture was purchased from Sigma Aldrich (Milan, Italy): PC (15:0–18:1) (d7), PE (15:0–18:1) (d7), PS (15:0–18:1) (d7), PG (15:0–18:1) (d7), PI (15:0–18:1) (d7), PA (15:0–18:1) (d7), LPC (18:1) (d7), LPC 25, LPE (18:1) (d7), Chol Ester (18:1) (d7), MG (18:1) (d7), DG (15:0–18:1) (d7), TG ((15:0–18:1) (d7)-15:0)), SM (18:1) (d9), cholesterol (d7).

3.2. Cell Culture

Freshly drawn blood (Rh+) from 9 healthy adults of both sexes (men and women) was used, heparin was added and preserved in citrate-phosphate-dextrose with adenine (CPDA-1). Data are the average ± SD of three independent experiments. RBCs were separated from plasma and leukocytes by washing three times with phosphate-buffered saline (127 mM NaCl, 2.7 mM KCl, 8.1 mM Na2HPO4, 1.5 mM KH2PO4, 20 mM HEPES, 1 mM MgCl2, and pH 7.4) supplemented with 5 mM glucose (PBS glucose) to obtain packed cells. This study was conducted in accordance with Good Clinical Practice guidelines and the Declaration of Helsinki. No ethical approval has been requested as human blood samples were used only to sustain in vitro cultures and patients provided written, informed consent in ASL. 1-Sassari (Azienda Sanitaria Locale. 1-Sassari) centre before entering the study.

3.3. Microgravity Simulation

In order to study the effects caused by microgravity on human erythrocytes, the gravity simulator 3D Random Positioning Machine (RPM, Fokker Space, Netherlands) was used at the laboratory of the Department of Biomedical Sciences, University of Sassari, Sardinia, Italy. The 3D Random Positioning Machine (RPM) is a micro-weight (‘microgravity’) simulator based on the principle of ‘gravity-vector-averaging’, built by Dutch Space. The 3D RPM is constructed from two perpendicular frames that rotate independently. This setup was used to constantly change the mean value of the gravity vector to zero. In this way, the 3D RPM provides a simulated microgravity less than 10−3 g. The dimensions of the 3D RPM are limited to 1000 × 800 × 1000 mm (length × width × height). The 3D RPM is connected to a computer, and through a specific software the mode and speed of rotation were selected. Random Walk mode with an 80 degree/s (rpm) was chosen.
The red blood cell samples were carefully deposited in 2 mL tubes together with PBS-glucose (30% haematocrit, approximately 3.4 × 109 cells) in a dedicated room at 37 °C. The control group samples were placed in the static bar at 1 g to undergo the same vibrations as the samples placed in µg conditions. Both control (1 g) and case (0 g) samples were collected after different time points (0, 6, 9, 24 h). Subsequently, the red blood cells were centrifuged and resuspended in 1 mL of lysis buffer [5 mM Na2HPO4, 1 mM EDTA (pH 8.0)] and stored at −20 °C until use for lipidomic analysis or fixed for confocal microscopic analysis.

3.4. Sample Preparation for UHPLC-IM-QTOF-MS Analysis

In order to investigate changes in the lipidome, analysis by UHPLC- IM-QTOF-MS requires the extraction of lipid content from cells [38]. An amount of 50 µL of human erythrocyte solution was extracted following the Folch procedure using 0.700 mL of a methanol and chloroform mixture (2/1, v/v). Samples were vortexed every 15 min up to 1 h, when 0.350 mL of chloroform and 0.150 mL of water were subsequently added. The solution thus obtained was centrifuged at 17,700 rcf for 10 min, and 0.600 mL of the organic layer was transferred into a glass vial and dried under a nitrogen stream. The dried chloroform phase was reconstituted with 50 μL of a methanol and chloroform mixture (1/1, v/v) and 75 μL isopropanol:acetonitrile:water mixture (2:1:1 v/v/v). Quality control (QC) samples were prepared taking an aliquot of 10 μL of each sample. All samples thus prepared were injected in UHPLC-IM-QTOF-MS/MS and acquired in negative ionization mode, while for positive ionization mode they were diluted in ratio 1:10.

3.5. UHPLC-IM-QTOF-MS/MS Analysis

The chloroform phase was analysed with a 6560-drift tube ion mobility LC-QTOF-MS coupled with an Agilent 1290 Infinity II LC system. An aliquot of 4.0 μL from each sample was injected in a Luna Omega C18, 1.6 μm, 100 mm × 2.1 mm chromatographic column (Phenomenex, Castel Maggiore (BO), Italy). The column was maintained at 50 °C at a flow rate of 0.4 mL/min. The mobile phase for positive ionization mode consisted of (A) 10 mM ammonium formate solution in 60% of milliQ water and 40% of acetonitrile and (B) 10 mM ammonium formate solution containing 90% of isopropanol and 10% of acetonitrile. In positive ionization mode, the chromatographic separation was obtained with the following gradient: initially, 80% of A, then a linear decrease from 80% to 50% of A in 2.1 min, then at 30% in 10 min. Subsequently, the mobile phase A was again decreased from 30% to 1% and stayed at this percentage for 1.9 min, and then was brought back to the initial conditions in 1 min. The mobile phase for the chromatographic separation in the negative ionization mode differed only for the use of 10 mM ammonium acetate instead of ammonium formate.
An Agilent jet stream technology source was operated in both positive and negative ion modes with the following parameters: gas temperature, 200 °C; gas flow (nitrogen) 10 L/min; nebulizer gas (nitrogen), 50 psig; sheath gas temperature, 300 °C; sheath gas flow, 12 L/min; capillary voltage 3500 V for positive and 3000 V for negative; nozzle voltage 0 V; fragmentor 150 V; skimmer 65 V, octapole RF 7550 V; mass range, 50−1700 m/z; capillary voltage, 3.5 kV; collision energy 20 eV in positive and 25 eV in negative mode, mass precursor per cycle = 3. High-purity nitrogen (99.999%) was used as a drift gas with a trap fill time and a trap release time of 2000 and 500 µs, respectively. Before the analysis, the instrument was calibrated using an Agilent tuning solution at the mass range of m/z 50–1700. Samples were evaporated with nitrogen at the pressure of 48 mTorr and at the temperature of 375 °C, while an Agilent reference mass mix for mass re-calibration was continuously injected during the run schedule.
The Agilent MassHunter LC/MS Acquisition console (revision B.09.00) from The MassHunter suite was used for data acquisition.

3.6. Data Analysis

Data acquired with the Agilent 6560 DTIM Q-TOF LC-MS were pre-processed with the software MassHunter Workstation suite (Agilent Technologies, Santa Clara, CA, USA). This software (Mass Profiler 10.0) allowed us to perform mass re-calibration, DTCCSN2 re-calibration, time alignment, and deconvolution of signals, yielding a matrix containing all features present across all samples. The removal of background noise and unrelated ions was performed by a recursive feature extraction tool, yielding a matrix containing all the features present across all samples. Furthermore, to eliminate non-specific information, data matrix quality assurance was performed. This filtered matrix was then subjected to multivariate statistical analysis using SIMCA software 15.0 (Umetrics, Umeå, Sweden).
First, a principal component analysis (PCA) was carried out. This unsupervised analysis allows an observation of samples and variables distribution in the multivariate space on the basis of their similarity and dissimilarity. This was followed by partial least square-discriminant analysis (PLS-DA) with its orthogonal extension (OPLS-DA), which was used as a classificatory model to visualize and evaluate the differences between sample classes.

4. Conclusions

Spatial anaemia in astronauts has been noted since the earliest space missions, while the contributing mechanisms during space flight remained unclear. To investigate the molecular mechanisms that induce a reduction in the number of erythrocytes during spaceflight, we decided to analyse the lipid profile of human erythrocytes under microgravity conditions. Thanks to the advancement of hyphenated techniques and mass analysers, we were able to identify biologically active complex lipids susceptible to microgravity, allowing new possible hypotheses that explain the anaemia experienced by astronauts.
In more detail, lipidomic analysis of erythrocytes revealed a double mechanism that generates the reduction in the number of red blood cells. On one hand, there is an increase in the levels of 20:4 PC, reducing cellular proliferation. On the other hand, the increase in the levels of EtherOxPC 16:0_20:4 stimulates the immune response by attracting the C-reactive protein and macrophages and induces an increase in ROS production [33,34]. ROS increase caused by microgravity inevitably induces mitochondrial damage and dysfunction as indicated by the accumulation of HexCer lipid species in clinorotated erythrocytes. This accumulation acts as a pro-apoptotic signal condemning the erythrocytes to death.
In this study, we reported a set of lipid discriminants or potential biomarkers linked to microgravity exposure with the aim to explore in the future specific lipid pathways and form the foundation in the development of novel therapeutics in hopes of reducing the effects of space flights. However, further studies are needed to accurately measure lipids in erythrocytes samples to better understand the clinical effects of microgravity.

Author Contributions

Conceptualization, P.C. and G.C.; methodology, C.M., A.M. (Antonio Murgia) and A.M. (Alessia Manca); writing—original draft preparation, C.M., P.C. and G.C.; writing—review and editing, C.M., A.M. (Antonio Murgia), A.M. (Alessia Manca), P.C., A.P. and G.C.; supervision, P.C. and G.C.; funding acquisition, P.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki. No ethical approval has been requested as human blood samples were used only to perform in vitro experiments.

Informed Consent Statement

Not needed.

Data Availability Statement

On request.

Acknowledgments

We acknowledge the CeSAR (Centro Servizi d’Ateneo per la Ricerca) of the University of Cagliari, Italy for the IMQTOFMSMS experiments performed with an Agilent 6560. C.M. and A.M. (Alessia Manca) also acknowledge the PhD program in Innovation Science and Technologies available at the University of Cagliari and Life Science and Biotechnologies available at the University of Sassari, respectively.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

fatty acid (FA); phosphatidylethanolamine (PE); ether-linked phosphatidylcholine (Ether-PC); phosphatidylcholine (PC); lysophosphocholine (LPC); sphingomyelin (SM).

References

  1. Riwaldt, S.; Corydon, T.J.; Pantalone, D.; Sahana, J.; Wise, P.; Wehland, M.; Krüger, M.; Melnik, D.; Kopp, S.; Infanger, M.; et al. Role of Apoptosis in Wound Healing and Apoptosis Alterations in Microgravity. Front. Bioeng. Biotechnol. 2021, 9, 1–22. [Google Scholar] [CrossRef]
  2. Bizzarri, M.; Fedeli, V.; Piombarolo, A.; Angeloni, A. Space Biomedicine: A Unique Opportunity to Rethink the Relationships between Physics and Biology. Biomedicines 2022, 10, 2633. [Google Scholar] [CrossRef]
  3. Thirsk, R.; Kuipers, A.; Mukai, C.; Williams, D. The space-flight environment: The International Space Station and beyond. CMAJ 2009, 180, 1216–1220. [Google Scholar] [CrossRef] [Green Version]
  4. Sheyn, D.; Pelled, G.; Netanely, D.; Domany, E.; Gazit, D. The Effect of Simulated Microgravity on Human Mesenchymal Stem Cells Cultured in an Osteogenic Differentiation System: A Bioinformatics Study. Tissue Eng.-Part A 2010, 16, 3403–3412. [Google Scholar] [CrossRef] [Green Version]
  5. Hughes-Fulford, M.; Tjandrawinata, R.; Fitzgerald, J.; Gasuad, K.; Gilbertson, V. Effects of microgravity on osteoblast growth. Gravit. Space Biol. Bull. 1998, 11, 51–60. [Google Scholar]
  6. Han, X.; Qiu, L.; Zhang, Y.; Kong, Q.; Wang, H.; Wang, H.; Li, H.; Duan, C.; Wang, Y.; Song, Y.; et al. Transplantation of Sertoli-Islet Cell Aggregates Formed by Microgravity: Prolonged Survival in Diabetic Rats. Exp. Biol. Med. 2009, 234, 595–603. [Google Scholar] [CrossRef]
  7. NSBRI. National Space Biomedical Research Institute Strategic Plan 2010. Biomed Res. 2010. Available online: http://nsbri.org/wp-content/uploads/2015/08/NSBRI_strategic_plan1 (accessed on 18 January 2022).
  8. Dinarelli, S.; Longo, G.; Dietler, G.; Francioso, A.; Mosca, L.; Pannitteri, G.; Boumis, G.; Bellelli, A.; Girasole, M. Erythrocyte’s aging in microgravity highlights how environmental stimuli shape metabolism and morphology. Sci. Rep. 2018, 8, 5277. [Google Scholar] [CrossRef]
  9. Poon, C. Factors implicating the validity and interpretation of mechanobiology studies in simulated microgravity environments. Eng. Rep. 2020, 2, 1–18. [Google Scholar] [CrossRef]
  10. Maxfield, F.R.; Tabas, I. Role of cholesterol and lipid organization in disease. Nature 2005, 438, 612–621. [Google Scholar] [CrossRef]
  11. Wenk, M.R. The emerging field of lipidomics. Nat. Rev. Drug Discov. 2005, 4, 594–610. [Google Scholar] [CrossRef]
  12. Carrasco-Pancorbo, A.; Navas-Iglesias, N.; Cuadros-Rodríguez, L. From lipid analysis towards lipidomics, a new challenge for the analytical chemistry of the 21st century. Part I: Modern lipid analysis. TrAC-Trends Anal. Chem. 2009, 28, 263–278. [Google Scholar] [CrossRef]
  13. Roberts, L.D.; Souza, A.L.; Gerszten, R.E.; Clish, C.B. Targeted Metabolomics. Curr. Protoc. Mol. Biol. 2012, 98, 30.2.1–30.2.24. [Google Scholar] [CrossRef]
  14. Wang, Z.; Klipfell, E.; Bennett, B.J.; Koeth, R.; Levison, B.S.; DuGar, B.; Feldstein, A.E.; Britt, E.B.; Fu, X.; Chung, Y.-M.; et al. Gut Flora Metabolism of Phosphatidylcholine Promotes Cardiovascular Disease. Nature 2011, 472, 57–63. [Google Scholar] [CrossRef] [Green Version]
  15. Cajka, T.; Fiehn, O. Toward Merging Untargeted and Targeted Methods in Mass Spectrometry-Based Metabolomics and Lipidomics. Anal. Chem. 2016, 88, 524–545. [Google Scholar] [CrossRef]
  16. Ivanova, S.M.; Morukov, B.V.; Labetskaia, O.I.; IuV, I.; Levina, A.A.; Kozinets, G.I. Morphobiochemical assay of the red blood system in members of the prime crews of the International Space Station. Hum. Physiol. 2009, 43, 43–47. [Google Scholar]
  17. Girasole, M.; Pompeo, G.; Cricenti, A.; Congiu-Castellano, A.; Andreola, F.; Serafino, A.; Frazer, B.; Boumis, G.; Amiconi, G. Roughness of the plasma membrane as an independent morphological parameter to study RBCs: A quantitative atomic force microscopy investigation. Biochim. Biophys. Acta (BBA)-Biomembr. 2007, 1768, 1268–1276. [Google Scholar] [CrossRef] [Green Version]
  18. Manis, C.; Manca, A.; Murgia, A.; Uras, G.; Caboni, P.; Congiu, T.; Faa, G.; Pantaleo, A.; Cao, G. Understanding the Behaviour of Human Cell Types under Simulated Microgravity Conditions: The Case of Erythrocytes. Int. J. Mol. Sci. 2022, 23, 6876. [Google Scholar] [CrossRef]
  19. Alfrey, C.P.; Udden, M.M.; Leach-Huntoon, C.; Driscoll, T.; Pickett, M.H. Control of red blood cell mass in spaceflight. J. Appl. Physiol. 1996, 81, 98–104. [Google Scholar] [CrossRef]
  20. Taylor, G.R. Advances in experimental medicine. In Advances in Experimental Medicine and Biology; Plenum Press: New York, NY, USA; London, UK, 1997; Volume 225, pp. 269–271. Available online: http://linkinghub.elsevier.com/retrieve/pii/0020711X84902064 (accessed on 18 January 2022).
  21. Meehan, R.T.; Neale, L.S.; Kraus, E.T.; Stuart, C.A.; Smith, M.L.; Cintron, N.M.; Sams, C.F. Alteration in human mononuclear leucocytes following space flight. Immunology 1992, 76, 491–497. [Google Scholar]
  22. Diedrich, A.; Paranjape, S.Y.; Robertson, D. Plasma and Blood Volume in Space. Am. J. Med. Sci. 2007, 334, 80–86. [Google Scholar] [CrossRef]
  23. Erslev, A.J. Erythropoietin. N. Engl. J. Med. 1991, 324, 1339–1344. [Google Scholar] [CrossRef]
  24. Koury, M.J.; Bondurant, M.C. Erythropoietin Retards DNA Breakdown and Prevents Programmed Death in Erythroid Progenitor Cells. Science 1990, 248, 378–381. [Google Scholar] [CrossRef]
  25. Trudel, G.; Shahin, N.; Ramsay, T.; Laneuville, O.; Louati, H. Hemolysis contributes to anemia during long-duration space flight. Nat. Med. 2022, 28, 59–62. [Google Scholar] [CrossRef]
  26. Crimi, M.; Degli Esposti, M. Apoptosis-induced changes in mitochondrial lipids. Biochim. Biophys. Acta (BBA)-Mol. Cell Res. 2011, 1813, 551–557. [Google Scholar] [CrossRef]
  27. Sprong, H.; Van Der Sluijs, P.; van Meer, G. How proteins move lipids and lipids move proteins. Nat. Rev. Mol. Cell Biol. 2001, 2, 504–513. [Google Scholar] [CrossRef]
  28. Robertson, J.L. The lipid bilayer membrane and its protein constituents. J. Gen. Physiol. 2018, 150, 1472–1483. [Google Scholar] [CrossRef] [Green Version]
  29. Koeberle, A.; Shindou, H.; Koeberle, S.C.; Laufer, S.A.; Shimizu, T.; Werz, O. Arachidonoyl-phosphatidylcholine oscillates during the cell cycle and counteracts proliferation by suppressing Akt membrane binding. Proc. Natl. Acad. Sci. USA 2013, 110, 2546–2551. [Google Scholar] [CrossRef] [Green Version]
  30. Picache, J.A.; Rose, B.S.; Balinski, A.; Leaptrot, K.L.; Sherrod, S.D.; May, J.C.; McLean, J.A. Collision cross section compendium to annotate and predict multi-omic compound identities. Chem. Sci. 2019, 10, 983–993. [Google Scholar] [CrossRef] [Green Version]
  31. Subbanagounder, G.; Leitinger, N.; Schwenke, D.C.; Wong, J.W.; Lee, H.; Rizza, C.; Watson, A.D.; Faull, K.F.; Fogelman, A.M.; Berliner, J.A. Determinants of bioactivity of oxidized phospholipids: Specific oxidized fatty acyl groups at the sn-2 position. Arterioscler. Thromb. Vasc. Biol. 2000, 20, 2248–2254. [Google Scholar] [CrossRef] [Green Version]
  32. Ran, F.; An, L.; Fan, Y.; Hang, H.; Wang, S. Simulated microgravity potentiates generation of reactive oxygen species in cells. Biophys. Rep. 2016, 2, 100–105. [Google Scholar] [CrossRef] [Green Version]
  33. Ashraf, M.Z.; Kar, N.S.; Podrez, E.A. Oxidized phospholipids: Biomarker for cardiovascular diseases. Int. J. Biochem. Cell Biol. 2009, 41, 1241–1244. [Google Scholar] [CrossRef] [Green Version]
  34. Huber, J.; Vales, A.; Mitulović, G.; Blumer, M.; Schmid, R.; Witztum, J.L.; Binder, B.R.; Leitinger, N. Oxidized Membrane Vesicles and Blebs from Apoptotic Cells Contain Biologically Active Oxidized Phospholipids That Induce Monocyte-Endothelial Interactions. Arterioscler. Thromb. Vasc. Biol. 2002, 22, 101–107. [Google Scholar] [CrossRef] [Green Version]
  35. Chang, M.-K.; Binder, C.J.; Torzewski, M.; Witztum, J.L. C-reactive protein binds to both oxidized LDL and apoptotic cells through recognition of a common ligand: Phosphorylcholine of oxidized phospholipids. Proc. Natl. Acad. Sci. USA 2002, 99, 13043–13048. [Google Scholar] [CrossRef] [Green Version]
  36. Rouhanizadeh, M.; Hwang, J.; Clempus, R.E.; Marcu, L.; Lassègue, B.; Sevanian, A.; Hsiai, T.K. Oxidized-1-palmitoyl-2-arachidonoyl-sn-glycero-3-phosphorylcholine induces vascular endothelial superoxide production: Implication of NADPH oxidase. Free Radic. Biol. Med. 2005, 39, 1512–1522. [Google Scholar] [CrossRef] [Green Version]
  37. Patwardhan, G.A.; Beverly, L.J.; Siskind, L.J. Sphingolipids and mitochondrial apoptosis. J. Bioenerg. Biomembr. 2016, 48, 153–168. [Google Scholar] [CrossRef] [Green Version]
  38. Ulmer, C.Z.; Jones, C.M.; Yost, R.A.; Garrett, T.J.; Bowden, J.A. Optimization of Folch, Bligh-Dyer, and Matyash sample-to-extraction solvent ratios for human plasma-based lipidomics studies. Anal. Chim. Acta 2018, 1037, 351–357. [Google Scholar] [CrossRef]
Figure 1. Representative total ion chromatograms of erythrocytes cultured in terrestrial gravity conditions and used as control (a) and compared with erythrocytes cultured in simulated gravity conditions (b) both collected after six hours of experiment.
Figure 1. Representative total ion chromatograms of erythrocytes cultured in terrestrial gravity conditions and used as control (a) and compared with erythrocytes cultured in simulated gravity conditions (b) both collected after six hours of experiment.
Ijms 24 04379 g001
Figure 2. PCA score plot of (A) positive ion mode data and (B) negative ion mode data. The black circles represent the clinorotated samples, while grey circles represent control samples. The PCA analysis for PIA shows the following validation parameters: R2X = 0.808 and Q2 = 0.629; for the NIA model: R2X = 0.721, Q2 = 0.572.
Figure 2. PCA score plot of (A) positive ion mode data and (B) negative ion mode data. The black circles represent the clinorotated samples, while grey circles represent control samples. The PCA analysis for PIA shows the following validation parameters: R2X = 0.808 and Q2 = 0.629; for the NIA model: R2X = 0.721, Q2 = 0.572.
Ijms 24 04379 g002
Figure 3. PLS-DA score plot. The grey circles represent the control samples, while black circles represent clinorotated erythrocytes samples. (a) 6 h clinorotated samples vs. 6 h control samples (PIA): R2X = 0.686, R2Y = 0.747, and Q2 = 0.434; (b) 9 h clinorotated samples vs. 9 h control samples (PIA): R2X = 0.849, R2Y = 0.945, and Q2 = 0.66; (c) 24 h clinorotated samples vs. 24 h control samples (PIA): R2X = 0.741, R2Y = 0.984, and Q2 = 0.421; (d) 6 h clinorotated samples vs. 6 h control samples (NIA): R2X = 0.776, R2Y = 0.994, and Q2 = 0.599; (e) 9 h clinorotated samples vs. 9 h control samples (NIA): R2X = 0.85, R2Y = 0.971, and Q2 = 0.147; (f) 24 h clinorotated samples vs. 24 h control samples (NIA): R2X = 0.692, R2Y = 0.995, and Q2 = 0.837.
Figure 3. PLS-DA score plot. The grey circles represent the control samples, while black circles represent clinorotated erythrocytes samples. (a) 6 h clinorotated samples vs. 6 h control samples (PIA): R2X = 0.686, R2Y = 0.747, and Q2 = 0.434; (b) 9 h clinorotated samples vs. 9 h control samples (PIA): R2X = 0.849, R2Y = 0.945, and Q2 = 0.66; (c) 24 h clinorotated samples vs. 24 h control samples (PIA): R2X = 0.741, R2Y = 0.984, and Q2 = 0.421; (d) 6 h clinorotated samples vs. 6 h control samples (NIA): R2X = 0.776, R2Y = 0.994, and Q2 = 0.599; (e) 9 h clinorotated samples vs. 9 h control samples (NIA): R2X = 0.85, R2Y = 0.971, and Q2 = 0.147; (f) 24 h clinorotated samples vs. 24 h control samples (NIA): R2X = 0.692, R2Y = 0.995, and Q2 = 0.837.
Ijms 24 04379 g003
Figure 4. OPLS-DA score plot. The grey circles represent the control samples, while black circles represent clinorotated erythrocytes samples. (a) 6 h clinorotated samples vs. 6 h control samples (PIA): R2X = 0.868, R2Y= 0.862, and Q2 = 0.503; (b) 9 h clinorotated samples vs. 9 h control samples (PIA): R2X = 0.646, R2Y= 0.708, and Q2 = 0.199; (c) 24 h clinorotated samples vs. 24 h control samples (PIA): R2X = 0.741, R2Y= 0.984, and Q2 = 0.284; (d) 6 h clinorotated samples vs. 6 h control samples (NIA): R2X = 0.776, R2Y= 0.994, and Q2 = 0.482; (e) 9 h clinorotated samples vs. 9 h control samples (NIA): R2X = 0.883, R2Y= 0.998, and Q2 = 0.033; (f) 24 h clinorotated samples vs. 24 h control samples (NIA): R2X = 0.788, R2Y= 0.998, and Q2 = 0.721.
Figure 4. OPLS-DA score plot. The grey circles represent the control samples, while black circles represent clinorotated erythrocytes samples. (a) 6 h clinorotated samples vs. 6 h control samples (PIA): R2X = 0.868, R2Y= 0.862, and Q2 = 0.503; (b) 9 h clinorotated samples vs. 9 h control samples (PIA): R2X = 0.646, R2Y= 0.708, and Q2 = 0.199; (c) 24 h clinorotated samples vs. 24 h control samples (PIA): R2X = 0.741, R2Y= 0.984, and Q2 = 0.284; (d) 6 h clinorotated samples vs. 6 h control samples (NIA): R2X = 0.776, R2Y= 0.994, and Q2 = 0.482; (e) 9 h clinorotated samples vs. 9 h control samples (NIA): R2X = 0.883, R2Y= 0.998, and Q2 = 0.033; (f) 24 h clinorotated samples vs. 24 h control samples (NIA): R2X = 0.788, R2Y= 0.998, and Q2 = 0.721.
Ijms 24 04379 g004
Table 1. Discriminant metabolites of human erythrocyte samples annotated by LC/DTIM-QTOF-MS.
Table 1. Discriminant metabolites of human erythrocyte samples annotated by LC/DTIM-QTOF-MS.
LipidAdductm/z Experimentalm/z TheoreticalΔ (ppm)RT
(min)
Fatty Acid CompositionDTCCSN2
2)
VIPSignificance
Level
Regulation in G0 Cells
6 h
LysoPC16:1+H+494.3216494.32415.16.5516:1231.581.73nsup
PC 32:2+H+730.5394730.53811.37.7732:2286.751.71**up
PC 38:5+H+808.5871808.58512.08.0518:1, 20:4292.871.50**up
PC 33:5+H+738.5070738.50680.37.74 287.331.26nsdown
SM d34:0+H+705.5930705.59053.58.0116:0, 18:0285.871.10*down
PC 38:4+H+810.6021810.60071.79.2318:0, 20:4295.111.05**up
PC 36:2+H+786.6031786.60073.09.4318:1, 18:1291.561.00*up
SM d 38:1+H+759.6328759.63756.011.31 297.551.00**down
PC 36:1+H+788.6192788.61643.210.43 293.880.99*up
FA 16:2-H251.2021251.20171.63.5916:2114.951.38**down
HexCer_AP t37:1+(C2H3O2)832.6152832.61641.510.4322:1, 15:0296.271.46**up
Etn-1-P-Cer 32:1-H687.5469687.54463.47.5114:1, 18:0265.471.13***up
EtherOxPC 36:4e + 1O+(CHO2)828.5761828.57812.48.316:0, 20:4293.501.07*up
SM d40:1+(CHO2)831.6632831.6597411.66 301.491.06nsup
PE 36:3-H740.5249740.523627.5018:2, 18:1270.581.02nsup
9 h
PC 38:6+H+806.5700806.56940.748.4818:2, 20:4293.141.79*up
SM d 38:1+H+759.6328759.63756.011.31 297.551.75**down
PC 35:5+H+766.5399766.53812.38.3315:1, 20:4280.171.73***up
PE 36:1+H+746.5712746.56942.110.7118:0, 18:1281.871.71**up
PC 36:5+H+780.5541780.55380.38.1116:1, 20:4228.331.70**up
PC 32:0+H+734.572734.56943.58.9216:0, 16:0284.831.64**up
PC 36:2+H+786.6031786.60073.09.2518:1, 18:1291.561.63**up
SM d34:0+H+705.5930705.59053.58.0116:0, 18:0285.871.61**up
PC 35:5+H+766.5407766.53813.48.51 288.551.60**down
PE 38:6+H+764.524764.52252.07.8916:0, 22:6279.421.59**up
PC 36:3+H+784.5872784.58512.78.2518:1, 18:2289.421.57**up
PC 35:4+H+768.5593768.55387.08.7315:0, 20:4289.521.56**up
SM d44:5+H+835.6691835.66881.011.6114:3, 30:2305.481.52**up
SM 42:3+H+811.6713811.66883.110.49 302.651.48**up
PC 38:7+H+804.554804.5 5380.37.9318:3, 20:4292.21.45nsup
PC 34:0+H+762.6032762.60073.310.27 290.861.38nsup
PC 36:5+H+780.5541780.55380.38.11 289.421.37*up
PC 36:4+H+782.5701782.56941.09.0816:0, 20:4229.671.36**up
PC 34:3+H+756.554759.55380.38.9218:1, 16:2286.351.33*up
PE 34:2+H+716.5251716.52253.58.3716:0, 18:2273.061.23nsup
PC 36:1+H+788.6192788.61643.210.43 293.881.21*up
SM d42:2+H+813.6870813.68443.211.61 305.351.21**up
PC 38:5+H+808.5870808.58512.38.0416:0, 22:5292.871.20nsup
PC 33:5+H+738.5070738.50680.37.74 287.331.07*down
PE 36:3+H+742.5399742.53812.58.5116:0, 20:3276.611.04nsup
PE 36:4+H+740.5239740.52252.08.1916:0, 20:4276.931.03**up
SM d42:2+H+813.6870813.68443.211.61 305.031.00**up
EtherOxPC 36:4e + 1O+(CHO2)828.5761828.57812.48.316:0, 20:4293.501.56*down
PE 36:2-H742.5403742.53921.59.5518:1, 18:1271.521.16**up
PC 38:6+(CHO2)850.5613850.56041.17.6218:2, 20:4296.341.01nsup
24 h
PC 35:4+H+768.5593768.55387.08.7315:0, 20:4289.521.84**up
PE 36:4+H+740.5239740.52252.08.1916:0, 20:4276.931.84**up
PC 36:2+H+786.6031786.60073.09.2518:1, 18:1291.561.64**up
PC 32:0+H+734.572734.56943.58.9216:0, 16:0284.831.61nsup
SM d42:2+H+813.6870813.68443.211.98 305.031.46*up
SM 42:3+H+811.6713811.66883.110.49 302.651.40**up
PC 36:3+H+784.5872784.58512.78.2518:1, 18:2289.421.35*up
PE 36:1+H+746.5712746.56942.110.7118:0, 18:1281.871.34nsup
PE 38:6+H+764.524764.52252.07.8918:2, 20:4279.421.21**up
PC 36:4+H+782.5701782.56941.09.0816:0, 20:4229.671.20***up
PC 36:5+H+780.5541780.55380.38.11 289.421.20*up
PC 38:4+H+810.6021810.60071.79.2318:0, 20:4295.111.19**up
PC 34:0+H+762.6032762.60073.310.27 290.861.13*up
PC 35:5+H+766.5399766.53812.38.3315:1, 20:4280.171.11*down
PE 36:3+H+742.5399742.53812.58.5216:0, 20:3276.611.11nsup
SM d44:5+H+835.6691835.66881.011.6114:3, 30:2305.481.09**up
SM d42:2+H+813.6870813.68443.211.62 305.351.05nsup
PC 35:5+H+766.5407766.53813.48.51 288.551.03*down
LysoPC16:1+H+494.3216494.32415.16.5516:1231.581.01*up
FA 16:2-H251.2021251.20171.63.5916:2114.951.99nsup
Cer 42:1+(CHO2)694.6369694.63552.115.2718:1, 24:0278.311.02*down
HexCer_AP t37:2+(C2H3O2)830.5979830.59992.59.4222:1, 15:1295.071.01nsdown
HexCer_AP t37:1+(C2H3O2)832.6152832.61641.510.4322:1, 15:0296.271.01nsdown
ns: not significant; *** p value < 0.001; ** p value < 0.01; * p value < 0.05.
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

Manis, C.; Murgia, A.; Manca, A.; Pantaleo, A.; Cao, G.; Caboni, P. Untargeted Lipidomics of Erythrocytes under Simulated Microgravity Conditions. Int. J. Mol. Sci. 2023, 24, 4379. https://doi.org/10.3390/ijms24054379

AMA Style

Manis C, Murgia A, Manca A, Pantaleo A, Cao G, Caboni P. Untargeted Lipidomics of Erythrocytes under Simulated Microgravity Conditions. International Journal of Molecular Sciences. 2023; 24(5):4379. https://doi.org/10.3390/ijms24054379

Chicago/Turabian Style

Manis, Cristina, Antonio Murgia, Alessia Manca, Antonella Pantaleo, Giacomo Cao, and Pierluigi Caboni. 2023. "Untargeted Lipidomics of Erythrocytes under Simulated Microgravity Conditions" International Journal of Molecular Sciences 24, no. 5: 4379. https://doi.org/10.3390/ijms24054379

APA Style

Manis, C., Murgia, A., Manca, A., Pantaleo, A., Cao, G., & Caboni, P. (2023). Untargeted Lipidomics of Erythrocytes under Simulated Microgravity Conditions. International Journal of Molecular Sciences, 24(5), 4379. https://doi.org/10.3390/ijms24054379

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