BioMed Research International
Volume 2021 (2021), Article ID 6677271, 16 pages
https://doi.org/10.1155/2021/6677271

Research Article

Blood Plasma Metabolic Profile of Newborns with Hypoxic-Ischaemic Encephalopathy by GC-MS

Yanjuan Jia,1,2 Xiaoni Jia,³ Hui Xu,1,2 Lan Gao,⁴ Chaojun Wei,1,2 Yonghong Li,1,2 Xia Liu,³ Xiaoling Gao,1,2 and Li Wei⁵

¹NHC Key Laboratory of Diagnosis and Therapy of Gastrointestinal Tumor, Gansu Provincial Hospital, Lanzhou, China, 730000
²The Institute of Clinical Research and Translational Medicine, Gansu Provincial Hospital, Lanzhou, China, 730000
³The Neonatal Department, Qingyang People's Hospital, Qingyang, China, 73000
⁴The Department of Life Science, Lanzhou University, Lanzhou, China, 730000
⁵The Laboratory Center, Gansu Provincial Hospital, Lanzhou, China, 730000

Correspondence should be addressed to Xiaoling Gao; gaoxl008@hotmail.com and Li Wei; weil18@lzu.edu.cn

Yanjuan Jia and Xiaoni Jia contributed equally to this work.

Received 26 October 2020; Revised 10 May 2021; Accepted 27 May 2021; Published 24 June 2021

Academic Editor: Cem Karadeniz

Copyright © 2021 Yanjuan Jia et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Background. Early diagnosis of hypoxic-ischaemic encephalopathy (HIE) is crucial in preventing neurodevelopmental disabilities and reducing morbidity and mortality. The study was to investigate the plasma metabolic signatures in the peripheral blood of HIE newborns and explore the potential diagnostic biomarkers. Method. In the present study, 24 newborns with HIE and 24 healthy controls were recruited. The plasma metabolites were measured by gas chromatography-mass spectrometry (GC-MS), and the raw data was standardized by the EigenMS method. Significantly differential metabolites were identified by multivariate statistics. Pathway enrichment was performed by bioinformatics analysis. Meanwhile, the diagnostic value of candidate biomarkers was evaluated. Result. The multivariate statistical models showed a robust capacity to distinguish the HIE cases from the controls. 52 metabolites were completely annotated. 331 significantly changed pathways were enriched based on seven databases, including 33 overlapped pathways. Most of them were related to amino acid metabolism, energy metabolism, neurotransmitter biosynthesis, pyrimidine metabolism, the regulation of HIF by oxygen, and GPCR downstream signaling. 14 candidate metabolites showed great diagnostic potential on HIE. Among them, alpha-ketoglutaric acid has the potential to assess the severity of HIE in particular. Conclusion. The blood plasma metabolic profile could comprehensively reflect the metabolic disorders of the whole body under hypoxia-ischaemic injury. Several candidate metabolites may serve as promising biomarkers for the early diagnosis of HIE. Further validation based on large clinical samples and the establishment of guidelines for the clinical application of mass spectrometry data standardization methods are imperative in the future.

1. Introduction

Hypoxic-ischaemic encephalopathy (HIE) is mainly caused by perinatal hypoxia-ischaemia, being responsible for nearly 23% of neonatal deaths worldwide [1, 2]. The incidence of HIE is 1.5 out of 1000 live births annually, which is higher in developing countries [3]. Approximately 25% to 30% of survivors may be at risk of lifelong neurological impairment and disabilities [4]. The efficiency of many current existing methods on HIE diagnosis, such as Apgar series score, cord blood gas, cranial ultrasound, and magnetic resonance imaging (MRI), have proven to be controversial and modest due to individual subjectivity, methodological shortcomings, and the complexity of the encephalopathy [5–12]. For instance, MRI, as one of the main popular diagnostic tools, often shows false-negative results due to its low sensitivity [13]. What is worse, unstable newborns may not be able to tolerate the transport or the duration of brain MRI scanning [14]. Indeed, no fast and reliable biomarker is currently available to help clinicians accurately assess the degree of the hypoxic-ischaemic (HI) injury [15].

In recent years, metabolomics is widely used to screen biomarkers for early diagnosis of disease and pathogenesis research [16]. The metabolic profiles have been established in biological fluid or tissue of several animal models for HIE and umbilical cord blood (UCB) obtained from HIE infants [17]. Several altered metabolites and perturbed metabolic pathways have been identified that may provide some information for establishing a panel of biomarkers for the early diagnosis of HIE. Actually, HI injury could seriously affect the integrity and permeability of the blood-brain barrier [18].Thus, various small molecules, including some metabolites, may be secreted into the peripheral circulation, and detection of the metabolic profile may play a vital role in the comprehensive understanding of HIE. So far, the blood plasma metabolic profile for HIE newborns compared with that for the normal remains unclear. Herein, we aimed to systematically analyze the plasma metabolic changes in the plasma of HIE newborns using metabolomics to explore a set of noninvasive biomarkers for the early diagnosis of HIE and early neuroprotection treatment.

2. Materials and Methods

2.1. Ethics Approval and Consent to Participate

The study was conducted according to the Declaration of Helsinki and approved by the Ethics Committee of the Gansu Provincial Hospital and Qingyang People’s Hospital (Gansu, China). The informed consent forms were signed by all included newborns’ parents before admission, agreeing on the use of clinical test data and the rest of human biosamples for scientific research by the hospital. The blood samples used in the study were the rest of the sample after the blood routine tests or the measurement of ABO blood type.

2.2. Study Design

76 newborn infants were recruited between November 2017 and March 2019. 28 were excluded because of sample hemolysis, severe jaundice, or insufficient quantity. 48 infants were left for the study. The HIE cases were enrolled from the Neonatal Intensive Care Unit of Qingyang People’s Hospital. Newborns with neural tube defects, alcohol and drug embryopathies, serious cardiac abnormality hydrops, important gastrointestinal malformation, congenital malformations, and hereditary metabolic diseases were excluded. Inclusion criteria for HIE group were as follows: (1) there were severe signs of fetal distress or an abnormal obstetric history related to fetal distress; (2) , ; (3) 5 min [8]; and (4) obvious nervous system symptom. Meanwhile, magnetic resonance imaging (MRI) was used to separate the severity of HIE [11, 12]. The injury was scored by using the modified Barkovich score based on the basal ganglia and thalami in combination with cortical involvement, the watershed areas, and the posterior limb of the internal capsule [19–22]. There were 14 mild HIE and 10 moderate/severe HIE included in the study. The control group was enrolled from the Obstetrics Department of Qingyang People’s Hospital. Inclusion criteria for the control group were as follows: (1) , ; (2) 5 min . All cases had demographic and clinical details collected.

2.3. Sample Collection and Storage

2 mL blood was drawn from the femoral or radial artery into the EDTA-K₂ tube after birth. Then, the samples were transported to the laboratory center of Qingyang People’s Hospital for a blood routine test by well-trained professionals within 20 minutes. Next, the samples were processed according to the following standard operating procedure: samples were centrifuged at at room temperature for 10 minutes, and then, two aliquots of 200 μL of plasma were transferred to new microtubes and stored at -80°C. All samples were transported to the Gansu Provincial Hospital for mass spectrometry detection using dry ice.

2.4. Sample Extraction

50 μL plasma was taken into the 1.5 mL Eppendorf (EP) tubes, extracted with 200 μL methanol, added 5 μL of adonitol (0.5 mg/mL stock in dH₂O) as internal standard, vortex mixing for the 30 s, and treated with ultrasound for 10 min (incubated in ice water), then centrifuged for 15 min at 12000 rpm, 4°C. After that, 180 μL supernatant of each sample was transferred into fresh 1.5 mL EP tubes, and 40 μL from each sample was taken into fresh 1.5 mL EP tubes as a quality control (QC) sample. Furthermore, the extracts were dried completely in a vacuum concentrator without heating, then 30 μL methoxyamination hydrochloride (20 mg/mL in pyridine) incubated for 30 min at 80°C was added, and 40 μL of the BSTFA regent (1% TMCS, ) to the sample aliquots was added, incubated for 1.5 h at 70°C. 5 μL of FAMEs (standard mixture of fatty acid methyl esters, C8-C30, 1 mg/mL, in chloroform) was added to the QC sample when cooling to the room temperature. All samples were analyzed by the gas chromatograph system coupled with a Pegasus HT time-of-flight mass spectrometer (GC-TOF-MS).

2.5. GC-MS Analysis

GC-TOF-MS analysis was performed using an Agilent 7890 gas chromatograph system equipped with a Pegasus HT time-of-flight mass spectrometer. The system has a DB-5MS capillary column coated with 5% diphenyl cross-linked with 95% dimethylpolysiloxane ( inner diameter, 0.25 μm film thickness; J&W Scientific, Folsom, CA, USA). A 1 μL aliquot of the analyte was injected in the splitless mode. Helium was used as the carrier gas with a purge flow rate of 3 mL min⁻¹ at the front inlet and a gas flow rate of 1 mL min⁻¹ through the column. The initial temperature was kept at 50°C for 1 min, then raised to 310°C at a rate of 20°C min⁻¹, then kept for 6 min at 310°C. The temperature of injection, transfer line, and ion source were 280, 280, and 250°C, respectively. The energy in the electron impact mode was -70 eV. The mass spectrometry data were obtained in the full-scan mode with the m/z range of 50-500 at a rate of 12.5 spectra per second after a solvent delay of 4.83 min.

2.6. Data Processing and Normalization

The raw peak exaction, data baseline filtering and baseline calibration, peak alignment, deconvolution analysis, peak identification, and peak area integration were carried out by using Chroma TOF 4.3X software of LECO Corporation and the LECO-Fiehn Rtx5 database [23]. The annotation of mass spectral data is the first step for the data interpretation. Mass spectrum match and retention index match were considered simultaneously in metabolite identification. Peaks detected in <50% of QC samples or in QC samples were removed [24]. The missing values were filled up by half of the minimum value. Referring to the method for the normalization of blood metabolomics data [25], the EigenMS method has been selected with higher entropy across groups for data normalization (Supplementary Table 1). The log2 transformation of the peak area and the normalization were performed on the MetaboGroupS platform (https://www.omicsolution.org/wukong/MetaboGroupS/) for further analysis.

2.7. Statistical Analysis and Pathway Enrichment

Student’s test and Fisher’s exact test were applied to analyze demographic and clinical data by SPSS version 22 (SPSS Inc., Chicago, IL, USA). Principal component analysis (PCA) was carried out to reduce the dimension of the data and visualize the distribution and the grouping of the samples. 95% confidence interval in the PCA score plot was used as the threshold to identify potential outliers in the dataset. Orthogonal partial least squares discriminant analysis (OPLS-DA) was applied to improve understanding of the variables responsible for the classification and obtain the contribution of each variable to the model. Both analysis approaches were conducted by SIMCA software (V15.0.2, Sartorius Stedim Data Analytics AB, Umea, Sweden). A 7-fold cross-validation was performed to check the robustness and predictive ability of the OPLS-DA model. The performance of the model was evaluated by the goodness of fit parameters and . A permutation test was performed to validate the OPLS-DA model by resampling the model 200 times under the random permutation. Significantly differential metabolites were identified based on the combination of the first principal component of variable importance in projection (VIP) with a two-tailed Student’s -test () between the control and HIE groups.

The biological functions and pathway analysis of significantly changed metabolites were enriched by the MetaboAnalyst 4.0 and the “Wu Kong” platform (https://www.omicsolution.org/wkomics/main/) based on several databases including the Kyoto Encyclopedia of Genes and Genomes (KEGG; http://www.genome.jp/kegg/), the PubChem Database (PCDB; https://www.ncbi.nlm.nih.gov/pccompound/), the Human Metabolome Database (HMDB; http://www.hmdb.ca/), and the Small Molecule Pathway Database (SMPDB; http://smpdb.ca/). The pathway enrichment analysis and pathway topology analysis of the differential metabolites were performed by Fisher’s exact test and out-degree centrality. Visualization of hierarchical clustering and the correlation of the differential metabolites were conducted by R software (V3.6.1). The enrichment results were visualized by Cytoscape software (V 3.7.1). The receiver operating characteristic (ROC) was performed on the OmicShare platform (https://www.omicshare.com/tools/Home/Soft/roc).

3. Results

3.1. Study Population

The baseline clinical characteristics of the study population are presented in Table 1. There was no significant difference in neonatal sex, method of delivery, nationality, and maternal age between the HIE group and the control group. Gestational age, Apgar 1, Apgar 5, birth weight, leucocyte count, and lymphocyte count manifested a significant difference between the two groups (Table 1).

Table 1: Baseline characteristics of the participants ().

Control ()
HIE ()
value

Gestational age (weeks, mean)

0.004

Gender (male, , %)
13 (54.2%)
16 (66.7)
0.376

Birth weight (g)

<0.001

Method of delivery

0.080

    Vaginal
9 (37.5)
15 (62.5)

    C-section
15 (62.5)
9 (37.5)

Apgar 1 (mean (min–max))
9 (9-10)
5 (2-8)
<0.001

Apgar 5 (mean (min–max))
10 (9-10)
7 (5-8)
<0.001

Maternal ethnicity

    Han
24
24
1.00

Maternal age (years)

0.180

Inflammatory markers

    Leucocyte count (10⁹ cells/L)

0.047

    Neutrophil count (10⁹ cells/L)

0.551

    Lymphocyte count (10⁹ cells/L)

0.001

    Monocyte count (10⁹ cells/L)

0.795

Other index

    RBC (10¹² cells/L)

0.892

    Hb (g/L)

0.634

    HCT (%)

0.278

    PLT (10⁹/L)

0.869

    MCV (fL)

0.138

    MCH (pg)

0.410

    MCHC (g/L)

0.407

    MPV (fL)

0.047

    CV (%)

0.942

    SD (fL)

0.361

    PDW

0.041

    PCT

0.637

Data presented as or percentage or median. HIE: hypoxic-ischaemic encephalopathy.

3.2. Identification of Differential Metabolites

The QC samples were distributed in the range of ±2std, suggesting the good stability of the system (Supplementary Figure 1). A total of 480 peaks were detected in both groups. After a series of data management including relative standard deviation denoising, missing values filled by half of the minimum value, and data normalization by the internal standard approach, finally, there were 343 peaks retained. Then, the EigenMS method was preferred to standardize the data based on the entropy across groups among the seven methods (Supplementary Table 1) [25]; the higher the entropy across groups, the better the model. The PCA score plot was made to examine the distribution of the HIE group and the control group (, ; Figure 1(a)). The model showed a clear separation of the HIE group from the control group. OPLS-DA score plots further confirmed that this model was efficient and reliable. Goodness-of-fit parameters of the model were high, with an explanation capacity of 0.987 and a prediction capacity of 0.964 () (Figure 1(b)).

(a)

(b)

(c)

Figure 1: Multivariate statistical analysis on spectral data of neonatal plasma in groups HIE and CON: (a) PCA scattered plot; (b) OPLS-DA scattered plot; (c) permutation test. HIE: hypoxic-ischaemic encephalopathy; CON: healthy group.

Moreover, the permutation tests showed that the and values from the constructed model (top right) were always higher than the simulated values of the permutation test (bottom left), suggesting that the OPLS-DA model was robust and not overfitted (Figure 1(c)). According to the and , 71 differential metabolites were identified, of which 46 metabolites were upregulated and 25 were downregulated (Figure 2(a)). The correlation among the differential metabolites was calculated by Pearson analysis and was shown in Figure 2(b).

(a)

(b)

Figure 2: Heat map of hierarchical clustering and correlation analysis for group HIE vs. CON: (a) hierarchical clustering; (b) correlation analysis. Red indicates positive correlation and blue indicates negative correlation. The nonsignificant correlation was marked with fork number. HIE: hypoxic-ischaemic encephalopathy; CON: healthy group.

3.3. Annotation and Classification of Metabolites

The identified metabolites were annotated based on the databases of HMDB, KEGG, and PCDB. 52 differential metabolites were annotated out of a total of 71, including 33 upregulated metabolites and 19 downregulated metabolites (Supplementary Table 2). For their chemical properties, the metabolite subclasses of these potential biomarkers were grouped into the following categories: “amino acids, peptides, and analogs,” “carbohydrates and carbohydrate conjugates,” “fatty acids and conjugates,” “dicarboxylic acids and derivatives,” “benzoic acids and derivatives,” “beta hydroxy acids and derivatives,” and others (Supplementary Table 2).

3.4. Perturbed Metabolic Pathways in HIE

Among the identified potential biomarkers, the metabolite pathway overrepresentation test was performed to comprehensively understand the biological process and altered pathways in HIE in comparison with the control group. The results showed that 331 pathways were significantly perturbed in HIE based on seven databases (, Supplementary Table 3), and the top 50 pathways were shown in the Sankey network (Figure 3(a)). Among them, 32 signaling pathways were jointly enriched in more than two databases (Supplementary Table 4), including GABA synthesis, release, reuptake and degradation, neurotransmitter uptake and metabolism in glial cells, pyrimidine metabolism, regulation of hypoxia-inducible factor (HIF) by oxygen, tricarboxylic acid (TCA) cycle, TCA cycle and respiratory electron transport, and tryptophan metabolism. Meanwhile, several pathways related to downregulated metabolites mainly included fatty acid biosynthesis, phospholipases, sphingomyelin metabolism, fatty acid activation, neurotransmitter release cycle, and free fatty acid receptors. Based on the SMPDB, the top 50 of enrichment results were shown in Figure 3(b). Based on the KEGG, significantly abnormal pathways including alanine, aspartate and glutamate metabolism, nitrogen metabolism, citrate cycle, taurine, and hypotaurine metabolism, arginine and proline metabolism, and D-glutamine and D-glutamate metabolism were marked in Figure 3(c). More than 37% of the differential metabolites were involved in the amino acid metabolism and energy metabolism and were reassembled into a brief plot (Figure 4(a)). Finally, integrating the results of pathway enrichment analysis with the published results on several animal models and UCB [26, 27], 16 metabolites were selected as a candidate biomarker set, including alanine, pyruvic acid, alpha-ketoglutaric acid, glutamic acid, succinic acid, phenylalanine1, cystine, valine, isoleucine, hydroxylamine, taurine, L-malic acid, and myo-inositol, glutamine, tyrosine, and arachidic acids. The relative expression levels of candidate biomarkers in both groups were shown in the radar plot (Figure 4(b)). Given the genome-scale network model of human metabolism, the tissue locations, subcellular locations, dysfunctional enzymes, drug action pathways, and diseases associated with these candidate biomarkers were enriched (Figure 5(a)). These blood metabolites were mainly located in the placenta, muscle, and myelin sheath (Figure 5(a)).

(a)

(b)

(c)

Figure 3: Pathway analysis of altered metabolites in neonatal plasma for group HIE vs. CON: (a) Sankey network based on Reactome, SMPDB, HMDB, WikiPathways, EHMN, HumanCyc, INOH, and KEGG; (b) barplot based on the SMPDB; (c) dot plot based on the KEGG. HIE: hypoxic-ischaemic encephalopathy; CON: healthy group.

(a)

(b)

Figure 4: Assembly of the alterations of amino acid metabolic pathways in neonatal plasma for group HIE vs. CON. The colors correspond to the level of metabolites: red indicates increased levels, while blue indicates decreased levels. (a) Overview of the altered pathways of amino acid and energy metabolism; (b) radar chart comparing two groups on sixteen variables. TCA: tricarboxylic acid; HIE: hypoxic-ischaemic encephalopathy; CON: healthy group.

(a)

(b)

Figure 5: The enrichment analysis and prediction of diagnostic value of candidate differential metabolites. (a) Enrichment analysis of the differential metabolites. Five parts involved in the map including tissue locations, subcellular locations, metabolite-associated drug action pathways, metabolite-associated diseases, and metabolite-associated dysfunction enzymes. The same shape represented the same group. The colors correspond to the level of differential metabolites: red indicates increased levels, while blue indicates decreased levels. (b) Diagnostic value of 16 candidate differential metabolites. AUC: area under the curve; ROC: receiver operating characteristic.

3.5. Diagnostic Capacity of the Candidate Metabolites

Through the logistic regression analysis, the diagnostic performance of each candidate biomarker was further evaluated (Figure 5(b)). The result showed that alanine 1, glutamic acid, glutamine 1, L-malic acid, succinic acid, pyruvic acid, and taurine showed an excellent performance with the area under the curve (AUC) more than 0.95 and with high specificity and sensitivity (>0.85). Arachidic acid, phenylalanine 1, alpha-ketoglutaric acid, and tyrosine 1 exhibited a good capacity () and with specificity and sensitivity (>0.70). Valine, myo-inositol, and isoleucine showed a moderate capacity (), suggesting that these metabolites could have potential diagnostic value for HIE except hydroxylamine and cystine (). Meanwhile, alpha-ketoglutaric acid and hydroxylamine had good diagnostic value in discriminating mild HIE from severe HIE () (Supplementary Figure 2).

4. Discussions

HIE is characterized by depression of the level of consciousness, often accompanied by a series of complex clinical symptoms, such as abnormal muscle tone and power, respiratory depression, cerebral nerve dysfunction, and seizures. The metabolic profiles from animal models of HIE demonstrated the complex multifactorial nature of neonatal HIE [28–30]. Our study firstly reported the blood plasma metabolic profile of newborns with or without HIE. Significantly altered metabolites and the perturbed pathway have been found in HIE. And the plasma levels of some metabolites showed a potential role as a biomarker for the rapid diagnosis of HIE at an early stage.

4.1. Energy Metabolism and Amino Acid Metabolism

In our study, more than 1/3 of changed metabolites in HIE was related to mitochondrial energy metabolism and amino acid metabolism. Abnormalities in energy metabolism mainly involved changes in the intermediates of the TCA cycle and changes in the energy sources. In the TCA cycle, three intermediates including alpha-ketoglutaric acid, succinic acid, and L-malic acid were predominantly elevated in HIE, which was consistent with the results of urine metabolomics [31]. It has been suggested that the elevated succinic acid may act as a marker of brain injury under HI [32]. Also, succinic acid is an important vascular growth factor and might contribute to the repair of brain injury [33]. As an alternate energy source, the increased branched-chain amino acid like valine and isoleucine, and carnitine have been found in our study, which may meet the emergency needs of the system in case of energy reduction [26]. Also, the reduced creatinine levels found in the blood may be associated with impaired energy transfer systems [27].

Some metabolites like amino acid usually have a dual role acting as a source of energy metabolism and as precursors for neurotransmitter synthesis. Similar to the results of UCB metabolomics, metabolic pathways associated with the synthesis of serotonin and dopamine such as tryptophan metabolism, tyrosine metabolism, and phenylalanine metabolism were altered in our research [34], suggesting the perturbation of brain neurotransmitter synthesis [35, 36]. Increased inhibitory neurotransmitters like taurine and alanine and elevated excitatory neurotransmitters like glutamic acid have been identified in higher abundance in the HIE group. This excitotoxicity may result from the depolarization of the cell membrane and the release of a large amount of glutamate under free radical damage [3] and further cause more cell damage [9, 37]. The results of the diagnostic analysis showed that 14 candidate metabolites that were disordered in amino acid metabolism and energy metabolism had high specificity and sensitivity in the diagnosis of HIE. In particular, alpha-ketoglutaric acid had a good ability to distinguish the severity of HIE. It may be a good and broad-spectrum candidate biomarker since it was involved in 181 changed pathways.

4.2. Other Metabolism

Pyrimidine metabolism and ketone body metabolism were aberrant in blood plasma metabolomics, which were also present in the UCB metabolomics [34]. The decrease of the arachidonic acid level discovered in the plasma may increase the free radical damage and reduce neuronal repair capacity in HIE because this metabolite is an inflammatory mediator and vasodilator which can defend against oxidative stress in the brain and participates in the growth and repair of neurons [38, 39]. As a lipid metabolism-related molecule, the elevated myo-inositol in PB plasma may indicate the impairment in lipid metabolism in HIE, which was also present at urinary metabolomics [31]. Besides, metabolite set enrichment analysis also revealed several disturbed pathways closely related to brain dysfunction [40, 41], including GABA synthesis, release, reuptake and degradation, fatty acid biosynthesis, and regulation of development of the central nervous system [42]. This further suggested that each metabolic pathway does not function independently but interacts with each other to perform its biological function, so the regulatory mechanism of this complex metabolic network needs to be further investigated.

4.3. Heterogeneity Analysis on the Metabolic Profile of HIE

In terms of different species, compared with the finding in mouse [27, 28], piglet [30, 43–45], and human primate [29], we found that only several coexpressed metabolites were present in the PB plasma of HIE. For the different biosample types, there were also few coexpressed metabolites among human PB, UCB, and urine [26, 31, 34, 46]. To sum up, this may result from (1) the intrinsic heterogeneity of biosamples determined by species, sample type, sampling time, and inclusion and exclusion criteria; (2) the difference in specificity and sensitivity of the different analytical platform including GC/MS, LC/MS, H-NMR, and self-established platform; and (3) the different means of data processing. Normalization is a challenging and essential part of data analysis. Different methods of standardization directly affected the results of data analysis [25]. In our study, the standard normalization method and EigenMS method were used to normalize the data. For the standard normalization method, the results showed that the PCA could not separate both groups, and only 11 significantly differential metabolites were matched from the KEGG database. In contrast, the EigenMS method is more suitable for processing data in our research, and 52 differential metabolites were identified. 6 significantly coexpressed differential metabolites were identified between both methods (Supplementary Figure 3). This suggested that the formulation of guidelines for standardized processing of mass spectrometry data is also an urgent problem to be solved in the future.

4.4. Limitation

The study has its limitations. First, the sample size is relatively small, and the result of diagnostic analysis did not rule out the interference of confounding factors including the weight, gender, mode of delivery, and gestational age. So, multicenter large-scale validation is required to verify the results. Second, factors of preanalysis such as the fasting before sampling were not well-controlled because of the urgency and uncontrollability of time during the delivery process. Third, the long-term follow-up with scientific outcome assessment will be necessary in future clinical research. In our study, we only conducted a telephone follow-up every 6 months for the included cases due to the economic limitation and the imperfect medical system. The outcome data were summarized as follows: 3 cases died, 1 case refused to return visit, 1 case failed to connect, and the remaining 19 cases had no symptoms of neurodevelopmental retardation. The third stage during the progression of HIE may last for months and even years. Finally, the integrated analysis of multiomics will be essential to explore useful biomarkers and new therapeutic targets for HIE diagnosis and treatment.

5. Conclusions

Generally, HIE is a complex disease accompanied by multiple organ injury. The blood plasma metabolic profile could comprehensively reflect the changes of the whole body under HI injury. Most of the altered pathways were associated with amino acid metabolism, energy metabolism, biosynthesis of neurotransmitters, tryptophan, and pyrimidine metabolism. 14 candidate metabolites showed good diagnostic potential. But, the verification of large and independent cohorts is necessary. Meanwhile, the clinical application and standardization of mass spectrometry is a great challenge and dilemma at present.

Data Availability

The data is available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

XLG and LW conceived and designed the experiments and wrote the final manuscript. YJJ and XNJ performed the experiment and data analysis and drafted the manuscript. LG contributed to the design and wrote the draft. HX, CJW, YHL, and XL assisted in the data collection. All the authors of this paper have read and approved the final version submitted. Yanjuan Jia and Xiaoni Jia are listed as co-first authors.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (grant number 81830372), Gansu Youth Science and Technology Fund Project (grant number 18JR3RA057), and Science and Technology Bureau of Chengguan District of Lanzhou (No. 2019-1-7).

Supplementary Materials

Supplementary Table 1: entropy results of each group of samples with seven normalization methods. Supplementary Table 2: key differential metabolites detected by gas chromatography-mass spectrometry (GC-MS) for HIE vs. CON. Supplementary Table 3: the metabolic pathways based on seven databases. Supplementary Table 4: the overlapped pathways based on seven databases. Supplementary Figure 1: the quality control diagram. Supplementary Figure 2: the prediction of diagnostic value of differential candidate metabolites on the severity of HIE. Supplementary Figure 3: comparison of differential metabolites screened by different data standardization methods. (Supplementary Materials)

References

R. E. Black, S. Cousens, H. L. Johnson et al., “Global, regional, and national causes of child mortality in 2008: a systematic analysis,” The Lancet, vol. 375, no. 9730, pp. 1969–1987, 2010.

S. Oza, J. E. Lawn, D. R. Hogan, C. Mathers, and S. N. Cousens, “Neonatal cause-of-death estimates for the early and late neonatal periods for 194 countries: 2000-2013,” Bulletin of the World Health Organization, vol. 93, no. 1, pp. 19–28, 2015.

X. Qin, J. Cheng, Y. Zhong et al., “Mechanism and treatment related to oxidative stress in neonatal hypoxic-ischemic encephalopathy,” Frontiers in Molecular Neuroscience, vol. 12, p. 88, 2019.

L.-L. Xiong, L.-L. Xue, Y. Jiang et al., “Suppression of PDGF induces neuronal apoptosis after neonatal cerebral hypoxia and ischemia by inhibiting P-PI3K and P-AKT signaling pathways,” Brain Research, vol. 1719, pp. 77–88, 2019.

N. Merchant and D. Azzopardi, “Early predictors of outcome in infants treated with hypothermia for hypoxic-ischaemic encephalopathy,” Developmental Medicine and Child Neurology, vol. 57, pp. 8–16, 2015.

P. Yeh, K. Emary, and L. Impey, “The relationship between umbilical cord arterial pH and serious adverse neonatal outcome: analysis of 51,519 consecutive validated samples,” BJOG : An International Journal of Obstetrics and Gynaecology, vol. 119, no. 7, pp. 824–831, 2012.

J. Slattery, A. Morgan, and J. Douglas, “Early sucking and swallowing problems as predictors of neurodevelopmental outcome in children with neonatal brain injury: a systematic review,” Developmental Medicine and Child Neurology, vol. 54, no. 9, pp. 796–806, 2012.

American Academy of Pediatrics, “The Apgar score,” Pediatrics, vol. 117, no. 4, pp. 1444–1447, 2006.

M. Douglas-Escobar and M. Weiss, “Hypoxic-ischemic encephalopathy: a review for the clinician,” JAMA Pediatrics, vol. 169, no. 4, pp. 397–403, 2015.

G. De Bernardo, R. De Santis, M. Giordano, D. Sordino, G. Buonocore, and S. Perrone, “Predict respiratory distress syndrome by umbilical cord blood gas analysis in newborns with reassuring Apgar score,” Italian Journal of Pediatrics, vol. 46, no. 1, p. 20, 2020.

A. Barkovich, B. Hajnal, D. Vigneron et al., “Prediction of neuromotor outcome in perinatal asphyxia: evaluation of MR scoring systems,” AJNR. American Journal of Neuroradiology, vol. 19, no. 1, pp. 143–149, 1998.

F. Cowan, M. Rutherford, F. Groenendaal et al., “Origin and timing of brain lesions in term infants with neonatal encephalopathy,” The Lancet, vol. 361, no. 9359, pp. 736–742, 2003.

N. Rollins, T. Booth, M. Morriss, P. Sanchez, R. Heyne, and L. Chalak, “Predictive value of neonatal MRI showing no or minor degrees of brain injury after hypothermia,” Pediatric Neurology, vol. 50, no. 5, pp. 447–451, 2014.

H. van Laerhoven, T. de Haan, M. Offringa, B. Post, and J. van der Lee, “Prognostic tests in term neonates with hypoxic-ischemic encephalopathy: a systematic review,” Pediatrics, vol. 131, no. 1, pp. 88–98, 2013.

E. Graham, A. Everett, J. Delpech, and F. Northington, “Blood biomarkers for evaluation of perinatal encephalopathy: state of the art,” Current Opinion in Pediatrics, vol. 30, no. 2, pp. 199–203, 2018.

T. Wang, K. Shao, Q. Chu et al., “Automics: an integrated platform for NMR-based metabonomics spectral processing and data analysis,” BMC Bioinformatics, vol. 10, no. 1, p. 83, 2009.

E. Locci, G. Bazzano, R. Demontis, A. Chighine, V. Fanos, and E. d’Aloja, “Exploring perinatal asphyxia by metabolomics,” Metabolites, vol. 10, no. 4, p. 141, 2020.

W. Lee, A. Michael-Titus, and D. Shah, “Hypoxic-ischaemic encephalopathy and the blood-brain barrier in neonates,” Developmental Neuroscience, vol. 39, no. 1-4, pp. 49–58, 2017.

T. E. Inder, R. W. Hunt, C. J. Morley et al., “Randomized trial of systemic hypothermia selectively protects the cortex on MRI in term hypoxic-ischemic encephalopathy,” The Journal of Pediatrics, vol. 145, no. 6, pp. 835–837, 2004.

M. Rutherford, L. A. Ramenghi, A. D. Edwards et al., “Assessment of brain tissue injury after moderate hypothermia in neonates with hypoxic-ischaemic encephalopathy: a nested substudy of a randomised controlled trial,” The Lancet Neurology, vol. 9, no. 1, pp. 39–45, 2010.

T. Alderliesten, L. de Vries, M. Benders, C. Koopman, and F. Groenendaal, “MR imaging and outcome of term neonates with perinatal asphyxia: value of diffusion-weighted MR imaging and ¹H MR spectroscopy,” Radiology, vol. 261, no. 1, pp. 235–242, 2011.

R. K. Lawrence and T. E. Inder, “Anatomic changes and imaging in assessing brain injury in the term infant,” Clinics in Perinatology, vol. 35, no. 4, pp. 679–693, 2008.

T. Kind, G. Wohlgemuth, D. Y. Lee et al., “FiehnLib: mass spectral and retention index libraries for metabolomics based on quadrupole and time-of-flight gas chromatography/mass spectrometry,” Analytical Chemistry, vol. 81, no. 24, pp. 10038–10048, 2009.

W. Dunn, D. Broadhurst, P. Begley et al., “Procedures for large-scale metabolic profiling of serum and plasma using gas chromatography and liquid chromatography coupled to mass spectrometry,” Nature Protocols, vol. 6, no. 7, pp. 1060–1083, 2011.

S. Wang, X. Chen, D. du et al., “MetaboGroup S: a group entropy-based web platform for evaluating normalization methods in blood metabolomics data from maintenance hemodialysis patients,” Analytical Chemistry, vol. 90, no. 18, pp. 11124–11130, 2018.

S. Chigusa, T. Moroi, and Y. Shoji, “State-of-the-art calculation of the decay rate of electroweak vacuum in the standard model,” Physical Review Letters, vol. 119, no. 21, article 211801, 2017.

J. Liu, R. A. Sheldon, M. R. Segal et al., “¹H nuclear magnetic resonance brain metabolomics in neonatal mice after hypoxia-ischemia distinguished normothermic recovery from mild hypothermia recoveries,” Pediatric Research, vol. 74, no. 2, pp. 170–179, 2013.

B. J. Blaise, L. Schwendimann, V. Chhor et al., “Persistently altered metabolic phenotype following perinatal excitotoxic brain injury,” Developmental Neuroscience, vol. 39, no. 1-4, pp. 182–191, 2017.

P. T. Chun, R. J. McPherson, L. C. Marney et al., “Serial plasma metabolites following hypoxic-ischemic encephalopathy in a nonhuman primate model,” Developmental Neuroscience, vol. 37, no. 2, pp. 161–171, 2015.

R. Solberg, J. Kuligowski, L. Pankratov et al., “Changes of the plasma metabolome of newly born piglets subjected to postnatal hypoxia and resuscitation with air,” Pediatric Research, vol. 80, no. 2, pp. 284–292, 2016.

E. Locci, A. Noto, M. Puddu et al., “A longitudinal 1H-NMR metabolomics analysis of urine from newborns with hypoxic-ischemic encephalopathy undergoing hypothermia therapy. Clinical and medical legal insights,” PLoS One, vol. 13, no. 4, article e0194267, 2018.

E. Chouchani, V. Pell, E. Gaude et al., “Ischaemic accumulation of succinate controls reperfusion injury through mitochondrial ROS,” Nature, vol. 515, no. 7527, pp. 431–435, 2014.

D. Hamel, M. Sanchez, F. Duhamel et al., “G-protein-coupled receptor 91 and succinate are key contributors in neonatal postcerebral hypoxia-ischemia recovery,” Arteriosclerosis, Thrombosis, and Vascular Biology, vol. 34, no. 2, pp. 285–293, 2014.

N. M. Denihan, J. A. Kirwan, B. H. Walsh et al., “Untargeted metabolomic analysis and pathway discovery in perinatal asphyxia and hypoxic-ischaemic encephalopathy,” Journal of Cerebral Blood Flow and Metabolism, vol. 39, no. 1, pp. 147–162, 2019.

L.-Y. Liu, H. J. Zhang, L. Y. Luo et al., “Blood and urinary metabolomic evidence validating traditional Chinese medicine diagnostic classification of major depressive disorder,” Chinese Medicine, vol. 13, no. 1, p. 53, 2018.

P. Zheng, H. C. Gao, Q. Li et al., “Plasma metabonomics as a novel diagnostic approach for major depressive disorder,” Journal of Proteome Research, vol. 11, no. 3, pp. 1741–1748, 2012.

N. Chirinian and N. Mann, “Therapeutic hypothermia for management of neonatal asphyxia: what nurses need to know,” Critical Care Nurse, vol. 31, no. 3, pp. e1–12, 2011.

Z. Wang, C. Liang, G. Li, C. Yu, and M. Yin, “Neuroprotective effects of arachidonic acid against oxidative stress on rat hippocampal slices,” Chemico-Biological Interactions, vol. 163, no. 3, pp. 207–217, 2006.

F. Darios and B. Davletov, “Omega-3 and omega-6 fatty acids stimulate cell membrane expansion by acting on syntaxin 3,” Nature, vol. 440, no. 7085, pp. 813–817, 2006.

M. Jan, “Cerebral palsy: comprehensive review and update,” Annals of Saudi Medicine, vol. 26, no. 2, pp. 123–132, 2006.

S. Bae, H. Kim, S. Shin, J. Kim, Y. Jeong, and J. Moon, “Decipher reliable biomarkers of brain aging by integrating literature-based evidence with interactome data,” Experimental & Molecular Medicine, vol. 50, no. 4, pp. 1–15, 2018.

B. Shah and K. Catt, “GPCR-mediated transactivation of RTKs in the CNS: mechanisms and consequences,” Trends in Neurosciences, vol. 27, no. 1, pp. 48–53, 2004.

Á. Sánchez-Illana, R. Solberg, I. Lliso et al., “Assessment of phospholipid synthesis related biomarkers for perinatal asphyxia: a piglet study,” Scientific Reports, vol. 7, no. 1, article 40315, 2017.

J. Kuligowski, R. Solberg, Á. Sánchez-Illana et al., “Plasma metabolite score correlates with hypoxia time in a newly born piglet model for asphyxia,” Redox Biology, vol. 12, pp. 1–7, 2017.

D. Sachse, A. L. Solevåg, J. P. Berg, and B. Nakstad, “The role of plasma and urine metabolomics in identifying new biomarkers in severe newborn asphyxia: a study of asphyxiated newborn pigs following cardiopulmonary resuscitation,” PLoS One, vol. 11, no. 8, article e0161123, 2016.

S. N. Reinke, B. H. Walsh, G. B. Boylan et al., “1H NMR derived metabolomic profile of neonatal asphyxia in umbilical cord serum: implications for hypoxic ischemic encephalopathy,” Journal of Proteome Research, vol. 12, no. 9, pp. 4230–4239, 2013.