1. Introduction
Hypoxic-ischemic encephalopathy (HIE) is a pathology that affects the brain and typically manifests within the first hours of life. It is characterized by a complex spectrum of neurological disorders of varying severity due to insufficient oxygen supply and reduced blood flow to the brain during childbirth. The incidence of severe HIE in the neonatal population is estimated at 0.37 to 3 cases per 1000 live births. [
1] HIE is particularly common in preterm infants with a gestational age of less than 28 weeks, with an incidence of 4 to 48 cases per 1000 pregnancies. It is worth noting that as gestational age decreases, the likelihood of HIE increases, indicating a relationship between prematurity and the development of this condition. [
2,
3]
Several factors influence the prognosis of HIE, including the degree, location, and extent of structural brain damage resulting from HIE. In general, the mortality rate associated with HIE is estimated to be between 9.1% and 9.9%, while in cases where the disease progresses to stage three, the mortality rate can increase significantly to 70-80%.[
4] It is important to note that approximately 25% of children affected by HIE suffer severe long-term neurological impairments, including conditions such as cerebral palsy, seizures, mental retardation, cognitive impairment, and epilepsy. In addition, intellectual impairment and behavioral problems, characterized by an intelligence quotient (IQ) of 70 or less, are observed in a significant proportion of HIE cases, ranging from 50 to 80%.[
5,
6,
7,
8] These impairments can have an essential impact on the affected child’s cognitive development and overall quality of life. Therefore, accurate diagnosis, timely intervention, and appropriate treatment of HIE are crucial to minimize the risk of potential long-term adverse consequences of this condition.
In the pathogenesis of HIE, phases of primary and secondary damage to nervous tissue are distinguished. Primary damage occurs at the period of asphyxia, and it is characterized by irreversible death of brain cells, the extent of which depends on the severity and duration of hypoxia. Secondary damage is activated in the reoxygenation-reperfusion phase 2-12 hours after the primary damage. Secondary damage is caused by the activation of several pathogenetic mechanisms: glutamate and calcium stress, free radical damage, aseptic inflammatory processes, and activation of apoptosis, which leads to an increase in the volume of neural damage and worsens the prognosis for life and health.[
9]
One of the promising approaches to reducing the adverse effects of central nervous system (CNS) injury is therapeutic hypothermia (TH), which is considered the most effective and safest neuroprotective method for infants who have suffered severe birth asphyxia. [
10,
11] It involves controlled chilling of the body to lower the core temperature to around 33-34 C for a specific time frame. Numerous experimental studies have shown that HT helps to reduce metabolic expenditures of the organism, secondary energy deficits in the cells, and the release of glutamate, also inhibiting free radical production and preventing inflammation and apoptosis.[
12,
13]
Early and accurate assessment of the severity of HIE remains one of the greatest challenges in neonatal care.[
14,
15] However, modern methods for assessing the risk of brain injury in neonates have limitations related to the first hours of life, and uncertainty about the severity of progressive brain injury and the potential neurological consequences remains at this early phase.[
16] For an early differential diagnosis of HIE in newborns, approaches to analyze metabolic changes in the blood of newborns using new analytical methods, in particular chromatography-mass spectrometry, might be suitable. [
17]This method allows the identification of hundreds of metabolites within a few hours, which fits into the time frame necessary to obtain analytical results and make decisions on treatment strategies. Clinical studies during the neonatal period are very scarce. Non-invasive samples (such as urine and stool) and minimally invasive samples (such as blood spots) are beneficial to study the dynamics of the molecular composition of biological fluids. Infants with congenital neonatal asphyxia typically do not urinate for a long time due to kidney damage. Therefore, the main focus of research is limited to circulating blood biomarkers.[
18] A dried blood spot is the most promising and least invasive sample for the diagnosis of neonatal pathologies.[
19] The analysis of amino acids and acylcarnitines in dried blood spots is a globally recognized approach for the screening of metabolic defects in newborns.[
20]
Recent advances in metabolomics, proteomics and transcriptomics have focused the attention of researchers on the identification of comprehensive molecular signatures. These signatures may allow precisely separate the effects of HIE in the first six hours of an infant's life and other injuries commonly caused by perinatal asphyxia, such as renal and myocardial damage. Neonatal HIE is difficult to distinguish from sepsis or inborn errors of metabolism [
18]. Sepsis may accompany HIE and alter the clinical signs. It is important to point out that incorrect diagnosis and, accordingly, the prescription of hypothermia can be harmful to patients with non-ischemic encephalopathies. [
21] This fact points to the need to seek for markers of brain damage in neonates that are specifically related to the ischemic/hypoxic genesis of the damage. [
22]
To simplify the multifactorial system in patients with HIE and to identify specific patterns of brain injury, proteomic and metabolomic analyzes are used in
in vivo and
in vitro models. Using an ischemia-hypoxia model (carotid artery ligation and exposure to a mixed gas of 92% N
2 / 8% O
2 for two hours) our recent study demonstrated the central role of glycerophospholipids, steroid biosynthesis and fatty acid metabolism in triggering systemic responses in 7-day-old rats [
23]. Other authors have performed metabolomic studies of hypoxic-ischemic injury in neonatal pigs after hypoxia and 120 minutes of reoxygenation phase. In blood plasma, during hypoxia the content of metabolites indicating a transition to anaerobic metabolism, in particular cytidine and uridine derivatives, free fatty acids and choline, increased significantly.[
24]
The persistent high morbidity and mortality in neonates with HIE underscores the urgent need for the introduction and implementation of improved diagnostic modalities. These advanced diagnostics are important not only to reliably confirm hypoxia/ischemia induced complications, but also to enable a rapid and prudent decision to initiate therapeutic hypothermia. In addition, such approaches have the potential to separate HIE and other clinical pathologies with phenotypically congruent symptoms, including but not limited to various inflammatory diseases. The integration of nuanced diagnostic measures is expected to improve the prognostic accuracy of HIE while streamlining the distinctive diagnostic processes to effectively distinguish HIE from a spectrum of alternative illnesses.
A routine neonatal screening platform using the non-derivatized MSMS kit NeoBase 2 (Perkin Elmer, Turku, Finland), which can detect 53 metabolites and more than 30 different metabolic diseases by tandem mass spectrometry (MSMS), was used.
3. Discussion
The diagnosis of HIE is a major challenge due to limitations in current clinical methodology and tools [
16,
25,
26]. Although amplitude-integrated EEG and imaging devices such as MRI are common techniques for predicting long-term outcome and assessing the severity of brain injury, they often fall short in the critical first hours of life when rapid and accurate diagnosis is most needed to make a decision whether to initiate TH [
14,
15,
27]. In addition, there are factors such as the need for continuous EEG monitoring to detect subclinical seizures that may go clinically unnoticed, and the possibility that neuroimaging may underestimate brain injury if performed too soon after the hypoxic event [
28]. In addition, comorbidities, including perinatal sepsis and metabolic disorders, may interfere with the neuroprotective effect of TH and complicate the diagnostic process [
29,
30,
31]. Although metabolic defects in inborn are rare, they may present symptoms like those under HIE [
32,
33].
In this complex diagnostic landscape, the analysis of metabolic markers proves to be a promising approach [
34]. Metabolomics offers a nuanced view of the biochemical changes that occur because of HIE and can provide earlier, and more precise indications of damage compared to conventional imaging and neurophysiological studies [
35,
36]. This method is particularly valuable in distinguishing between HIE and other disorders with similar manifestations, such as inborn metabolic defects or genetic encephalopathies, which can mimic HIE symptoms [
37]. Consideration of a detailed family history and targeted metabolic testing are critical to identifying these conditions, which require treatment approaches other than HIE, including specific treatments instead of therapeutic hypothermia. With advances in neonatal care research, the integration of metabolomic data into existing diagnostic tools holds the potential to increase the accuracy of HIE diagnoses and personalize interventions to improve outcomes for affected infants.
Although metabolomics is a promising approach for the diagnosis of HIE, it is important to note that many of the approaches being investigated in this area are still staying in experimental field and require further standardization before they can be used in clinical practice. However, there are approaches that are already being used in clinical practice, such as the NeoBase 2 MSMS kit, which is currently used in newborn screening [
38]. It allows the simultaneous analysis of a comprehensive panel of metabolites involved in different metabolic pathways, such as amino acids, carnitines, nucleosides and lysophosphatidylcholines, a total of 57 metabolites. The main limitation of this set in the context of searching for HIE markers is of course the small number of metabolites compared to standard mass spectrometry. However, we assumed that this set of metabolites would be sufficient for the differential diagnosis of HIE. In particular, the applicability of this kit for the diagnosis of various pathological conditions in addition to conventional newborn screening has already been demonstrated [
39,
40,
41]. Also, the NeoBase 2 MSMS kit takes advantage of the convenience and practicality of analyzing dried blood spots on paper, which are widely available and commonly used as a source of neonatal biological samples [
42]. This eliminates the need for invasive procedures and simplifies sample collection, transportation and storage.
Research on HIE markers is constrained by strict ethical standards associated with neonatal studies and the low incidence rate of HIE, limiting the clinical research capacity to fully validate diagnostic tests. To overcome these hurdles, experimental animal models of HIE are invaluable as they provide the opportunity to study the pathogenetic mechanisms and markers of brain injury in a controlled environment. In our study, we used the Rice-Vannucci model, a widely recognized experimental model that closely mimics the pathophysiological changes observed in neonates with HIE [
43]. We have also modeled one of the most common conditions of neonatal sepsis [
44], which complicates the diagnosis of HIE and may also be a contraindication to hypothermic therapy when sepsis and HIE are combined. Additionally, we assessed how TH impacts the profile of markers, particularly focusing on prognostic implications.
In the all three experiments, we observed a statistically significant increase in the levels of the acylcarnitines mainly short-chain acylcarnitine. Acylcarnitines are formed by the esterification of carnitine with acyl groups, which is needed to transport fatty acids through the mitochondrial membrane for subsequent beta-oxidation yielding production of energy. During hypoxia-ischemia, energy production in the cells is disturbed, which leads to a change in metabolism and redirection of bioenergetics toward anaerobic glycolysis. This in turn affects the levels of various metabolites, including acylcarnitines. The accumulation of certain acylcarnitines may reflect a disturbance in mitochondrial function and energy metabolism. Clinical studies have demonstrated that changes in the acylcarnitine profile can serve as a biomarker for HIE. A study examining a cohort of 67 children diagnosed with HIE revealed that butyrylcarnitine (C4) could serve as a potential prognostic biomarker. The observed correlations between circulating C4 levels, neuron-specific enolase concentrations and MRI abnormalities suggest that this biomarker is useful for the prognostic assessment of HIE [
45]. Analysis of umbilical cord blood showed elevated concentrations of long-chain acylcarnitines in two different groups of newborns: those who had clinical signs of hypoxia-ischemia but had a normal neurological examination, and those who were diagnosed with HIE [
46]. Another study in which newborns were diagnosed with HIE revealed elevated levels of acylcarnitines, specifically palmitoylcarnitine (C16:0), stearoylcarnitine (C18:0) and oleoylcarnitine (C18:1), which were corelated with more acidic pH values in the umbilical artery and lower Apgar scores five minutes after delivery [
47]. Dave M. and co-authors explored the tissue-specific origin of acylcarnitines using a mouse model of HIE and quantified these metabolites in skeletal muscle, brain and blood. They found that acylcarnitines C16:0 and C18:0 increase in plasma as well as in brain tissue 30 minutes after HIE and dissipate after 24 hours. On the other hand, all other acylcarnitines were elevated only in muscle and plasma. Authors made conclusions that plasma acylcarnitine profiles predominantly reflect muscle-derived changes rather than those occurring in the brain [
48]. In our study, after applying the FDR correction, acylcarnitines were excluded as reliable biomarkers for the diagnosis of HIE, although they were initially statistically significant. Due to predominant skeletal muscle origin of acylcarnitines following systemic hypoxic/ischemic events, the early post-HIE elevation in their levels may be more indicative for the presence of systemic hypoxia rather than the index of extent the cerebral injury. Consequently, acylcarnitines could serve as markers of systemic hypoxia in neonates and indirectly point to possible injury to a brain, the organ most susceptible to oxygen deficiency. However, comprehensive, long-term, multi-center studies of changes in acylcarnitine levels in neonates with HIE are needed for definitive conclusions.
In modeling a systemic inflammatory process associated with HIE, our study found change levels of metabolites responsible for several dysfunctions, including decreased oxygenation, excitotoxicity, impaired glucose utilization, impaired fatty acid metabolism, impaired protein synthesis, deterioration of cell membrane integrity, mitochondrial dysfunction, and impaired energy production. We found a decrease in short-chain acylcarnitine levels during hypoxic-ischemic injury, followed by an increase associated with increase in inflammatory response, which could indicate carnitine deficiency. Carnitine deficiency is a contributing factor to the progression of inflammatory processes [
49]. Although the analysis involving a panel of 57 metabolites provided insights into the pathophysiological processes, principal component analysis failed to distinguish between groups subjected to 24-hour LPS treatment and those modeled with HIE concurrent with systemic inflammation, as significant clustering was only observed within the control group. Further studies with a broader range of metabolites or alternative analytical methods would be required to refine the distinction between these conditions.
In the context of hypothermia after HIE, we observed a significant decrease in markers associated with HIE, such as acylcarnitines and the amino acids identified in the OPLS-DA diagnostic model. In addition, we observed significant changes in metabolic pathways responsible for ammonia detoxification, reduced lipid peroxidation activity and stabilization of cellular energy metabolism after TH. These results suggest a positive effect of TH on biochemical processes and recovery mechanisms after HIE, which is consistent with previous studies in this field. [
50,
51]. Furthermore, our results emphasize the sensitivity of the investigated markers to the extent of brain damage after HIE. Considering these results, we emphasize the potential utility of the NeoBase 2 MSMS kit for exploring neuroprotective strategies in animal models during the acute phase of the disease.
In this study, OPLS-DA and logistic regression models with high diagnostic accuracy were developed for timely diagnosis of HIE. Moreover, three- and six-hours’ models were similar in respect to included variables and performance. The models based on the DBS profile at six hours were validated using a new set of samples from experiment 3. Validation demonstrated sensitivity and specificity more than 86%. Interestingly, glycine was the only DBS compound included in logistic three- and six-hours’ regression models. Its VIP was the highest in both OPLS-DA models, which also included various amino acids. These compounds may be suggested as HIE markers in DBS. The main marker for HIE that we have identified is glycine, which is consistent with clinical data for newborns with HIE [
52,
53]. Also, an increase in glycine has been described for models of global brain ischemia [
53,
54,
55] and for patients in the acute phase of stroke in both plasma and cerebrospinal fluid [
56,
57]. We suggest that the main source of glycine could be a product of NADH synthesis reaction from serine under hypoxic conditions [
58]. NADH plays a vital role in the aerobic synthesis of ATP. Independently of the tricarboxylic acid cycle, serine catabolism within the mitochondria generates NADH through a folate-dependent mechanism. The enzyme serine hydroxymethyltransferase 2 catalyzes the conversion of serine to glycine, concurrently generating 5,10-methylene tetrahydrofolate (methylene-THF). Subsequently, methylene tetrahydrofolate dehydrogenase 2 oxidizes methylene-THF, with reduction of NAD+ and formation of 10-formyl-THF. It is also important to consider that hypoxic/ischemic conditions can disrupt neurotransmitter homeostasis, potentially leading to extracellular accumulation of neurotransmitters, including glycine, which can subsequently be released into the bloodstream [
56]. This phenomenon has been observed particularly in ischemic stroke model, with the increase in glycine concentration correlating with the extent of injury [
59].
In addition to glycine, the OPLS-DA diagnostic model also includes glutamate, leucine, isoleucine, hydroxyproline, citrulline, phenylalanine and the combined markers glutamine/lysine and valine as main markers of HIE. Glutamate in particular is known to be a key excitatory neurotransmitter in the brain and plays an important role in acute neurological conditions including HIE [
60]. In neonatal HIE, there is a critical disturbance in the energy metabolism of brain cells. The interruption of oxygen and glucose supply impedes respiration and oxidative phosphorylation and leads to a failure of ATP-dependent ions pumping. This leads to a pathological reversal of the functions of the glutamate transporter functions, resulting in a significant efflux of glutamate into the extracellular space. The result is a chain reaction of excitotoxicity that exacerbates neuronal damage and is reflected in increased plasma glutamate concentrations [
61]. In addition, hypoxia/ischemia leads to a metabolic shift toward anaerobic metabolic pathways, disabling the Krebs cycle and leading to an accumulation of cycle intermediates and byproducts of pyruvate metabolism, potentially increasing the concentrations of precursor amino acids identified in our study. The role of amino acids as biomarkers for early HIE has previously been demonstrated in both model experiments and clinical studies [
17,
46,
62,
63,
64]. Furthermore, these amino acids could serve as alternative energy substrates for the neonatal brain and muscles under hypoxic conditions [
65]. It is important to note that the enrichment of metabolic pathways that we identified 3 and 6 hours after HIE indicates the response of the brain and entire organism to hypoxic/ischemic damage, including a transition to anaerobic metabolism, compromised cellular energy production, perturbed glucose metabolism and protein biosynthesis, elevated oxidative stress, and disrupted reparation processes. The diagnostic biomarkers we have identified are intrinsically linked to these early pathophysiological events and are pathogenetically relevant. They denote the initiation of a patho-biochemical cascade that underlies the cerebral damage associated with HIE. Therefore, these markers can serve as critical indicators for the early detection of brain injury and may assist in guiding therapeutic interventions during the initial phase of HIE.
4. Materials and Methods
4.1 Animals
The animal protocols performed were reviewed and approved by the institutional animal ethics committee according to the guidelines of FELASA (Federation of European Laboratory Animal Science Associations). The experiments were performed with outbred white rats obtained from the animal facility of the A.N. Belozersky Institute of Physico-Chemical Biology. To ensure animal welfare, dams and their pups were housed in cages in a temperature-controlled environment with a temperature of 21±2°C. In addition, a light-dark cycle was established in which the light was switched on from 9:00 am to 9:00 pm. The dams had ad libitum access to food and water, and pups were examined daily for health status.
4.2. HIE Modeling
The study used the well-established Rice-Vannucci rat model to investigate hypoxic-ischemic brain injury (HIE). The experiments were conducted with seven-day-old postnatal rats of both sexes. The left carotid artery was then surgically isolated and obliterated with an electrocautery to deliberately restrict the blood supply to the brain. The rats were then exposed to a gas mixture of 8% oxygen and 92% nitrogen for a defined period of 2 hours in a multigas CO2 incubator (Binder, Tuttlingen, Germany) at a temperature of 37◦C, starting 1.5 hours after surgery, to induce hypoxia. The mortality rate after hypoxia-ischemia induction was in a range of 3–5%. The experiments were conducted to explore different impacts on rat pups with HIE, which were divided into experimental groups as follows.
4.2.1.HIE with Sampling at Different Times
Pups were randomly divided into the following experimental groups: intact rats (n = 10); animals with HIE whose blood was taken 3 hours after hypoxia (n = 13); animals with HIE whose blood was taken 6 hours after hypoxia (n = 12) (
Figure 1A). Three and six hours after the hypoxic exposure, dry blood spots were obtained.
4.2.2. HIE Model with LPS-Induced Inflammation
In this experiment, 6-day-old rat pups were injected intraperitoneally with 4 mg/kg body weight of lipopolysaccharide (LPS) 24 hours before the induction of hypoxia and ischemia. The animals were divided into three experimental groups: intact rats (n = 11), rats injected with LPS without HIE (n = 11) and rats injected with LPS with HIE (n = 27) (
Figure 1B). After 6 hours of hypoxic exposure, dry blood spots and blood plasma were collected for analysis.
4.2.3. Modeling of therapeutic hypothermia
Pups were randomly assigned to either normothermic or hypothermic recovery immediately after termination of hypoxia for a duration of 6 hours. Normothermic recovery was maintained at a temperature of 37°C, while hypothermic recovery was maintained at 30°C according to the previously described protocol (PMID: 25996893). Pups were placed in an open container in a water bath to achieve the desired temperatures. In this experimental set there were three groups of animals: intact rats (n = 13), rats with HIE and hypothermia (n = 16), and rats with HIE and normothermia (n = 14) (
Figure 1C). As in the previous experiments, dry blood spots (DBS) were collected after 6 hours of hypoxic exposure for further analysis.
4.3. Dried Blood Spot metabolome analysis (FIA-MRM-MS)
The levels of 57 low-molecular-weight compounds, including 14 amino acids, 36 carnitines, 2 nucleosides, succinylacetone and 4 lysophosphatidylcholines (
Table S7 Supplementary S2), were measured in DBS using the flow injection analysis tandem mass spectrometry (FIA-MS/MS) as a part of a routine newborn screening at the Clinical Diagnostic Laboratory of the National Medical Research Center for Obstetrics Gynecology and Perinatology named after Academician V.I. Kulakov of the Ministry of Healthcare of the Russian Federation. Experiments were performed on UHPLC QSight MD autosampler coupled to a QSight 225MD tandem mass spectrometer (MSMS) in a positive electrospray (ESI) mode using multiple reaction monitoring (MRM) for signal acquisition. The NeoBase 2 non-derivatized MSMS kit #3044-0010 (Perkin Elmer, Turku, Finland) was applied to extract and quantitatively analyze target compounds in DBS. The kit includes NeoBase 2 Internal Standards, Extraction Working Solution (EWS), NeoBase 2 Controls (Low and High), Neo MSMS Flow Solvent, NeoBase 2 Extraction Solution, NeoBase 2 Succinylacetone Assay Solution. NeoBase 2 assay was assessed based on the recommendations of the CLSI guidelines [
66]. Several studies confirmed good analytical performance of Neobase 2 kit in terms of precision, accuracy, limit of detection (LOD), lower limit of quantification (LLOQ), linearity, recovery, and carryover [
67,
68,
69,
70].
All analytical steps were performed in accordance with the manufacturer’s guidelines. In brief, a 3.2 mm disc was automatically punched from the DBS of each sample into a 96-well plate using Panthera-Pancher 9 (Perkin Elmer, Turku, Finland). DBS spots were incubated in the shaker with 125 µl of the extraction working solution (EWS) including internal standards (IS) for 30 minutes +45 °C. Then 100 µl supernatant was transferred to another plate’s well for 60-minute incubation at room temperature. The 10 µL injection into the ion source was performed under the following conditions: capillary voltage 5.5 kV, source temperature 120°C, hot source induced desolvation (HSID) temperature 320°C, drying gas flow 80 L/hr, and nebulizer gas flow 120 L/hr. FIA-MRM-MS was performed at solvent flow rate 1 µL/sec with duration run 1.5 minutes and scanning rate 1.76 sec/scan. MRM parameters are listed in the
Table S7 Supplementary S2. The low and high level controls (LC, HC) were run at the start and the end of in each plate. The mean measured QCs were within ± 3 SD of the values provided by the quality control certificate. Processing of the data and analytes concentration calculation was carried out by Simplicity 3Q MD 2.1 (Perkin Elmer, Turku, Finland) software.
Figure S9 Supplementary S2 presents an example of a typical flow injection chromatogram obtained for a LC sample provided by NeoBase™ 2 Non-derivatized MSMS kit on the QSight 225MD tandem mass spectrometer.
4.4. Statistical Methods
The groups were compared pairwise by the Mann-Whitney test, and p-values were corrected using the Benjamin-Hochberg method (false discovery rate, FDR). A p-value less than 0.05 was considered statistically significant. Compounds with statistically significant differences in levels were used for principal component analysis (PCA).
To investigate the correlation between the time elapsed since hypoxic-ischemic brain injury and the molecular composition of blood, a Spearman test was applied with criteria of significant P<0.05.
For each task, a pairwise comparative analysis of interactions between compounds within groups was performed: a correlation matrix and an adjacency change matrix were calculated for two groups. Based on the adjacency change matrix, a difference matrix was computed for hierarchical cluster analysis using the Ward’s method. The optimal number and composition of clusters were determined based on the maximum silhouette coefficient.
Pathway enrichment analysis was performed using the over-representation analysis in the MetabolanalystR package. A result with p<0.1 was considered statistically significant. If a marker set contained compounds with multiple identification variants, new variable sets were created for each identification variant, each containing one of the options. In this case, pathways were considered statistically significant if p was less than 0.1 for each marker set. Enrichment values and probability of observing were averaged across all marker sets, and the median value was calculated for the number of markers in each pathway.
For hypoxia diagnosis, classification models were built based on OPLS-DA. Sensitivity and specificity were determined based on leave-one-out cross-validation results.
For diagnostic models based on logistic regression, a new dataset was created, including metabolite concentrations, the product of concentrations for each pair of metabolites, and the squares of the concentration values for each metabolite. In the preliminary marker set, variables with a variable projection value greater than 1 in orthogonal projections analysis were selected. A stepwise variable selection was performed based on the Akaike information criterion to obtain the minimum criterion value at each stage, while observing a decrease in the criterion value. Subsequently, a stepwise exclusion of variables was carried out, with coefficients having a maximum probability of zero equality, until coefficients for variables had a probability of zero less than 0.05. The model's quality was evaluated based on leave-one-out cross-validation.
Additional validation of 6 h’ models (OPLS and logistic regression) was performed by 100-times cross-validation using train/test split (70/30) of first experiment samples (
Figure 1A, control (n=10) and 6h after hypoxia-ischemia (n=12) groups) with addition of samples from the third experiment (
Figure 1C, control (n=13) and normothermic hypoxia-ischemia (n=14) groups).
Author Contributions
Conceptualization, D.S., O.I., and N.S.; methodology, C.E.; software, C.E. and E.P.; validation, Y.S. and E.P.; formal analysis, Y.S. and K.G.; investigation, V.C., Y.S., K.G., and C.E.; resources, D.S. and E.P.; data curation, K.G. and Y.S.; writing—original draft preparation, Y.S., C.E., K.G., A.S. and N.S.; writing—review and editing, O.I., E.P., G.S., D.Z. and D.S.; visualization, K.G., N.S., and C.E.; supervision, G.S., N.S. and D.S.; project administration, D.S. and O.I.; funding acquisition, D.S. All authors have read and agreed to the published version of the manuscript.