- 1Department of Obstetrics and Gynecology, National Clinical Research Centre for Obstetric and Gynecologic Diseases, State Key Laboratory of Complex Severe and Rare Diseases, Peking Union Medical College Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
- 2State Key Laboratory of Complex Severe and Rare Diseases, Medical Research Center, Peking Union Medical College Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
Preeclampsia (PE) is one of the leading causes of maternal and perinatal morbidity and mortality. However, it is still uncertain how PE affects neonate metabolism. We conducted an untargeted metabolomics analysis of cord blood to explore the metabolic changes in PE neonates. Umbilical cord serum samples from neonates with preeclampsia (n = 29) and non-preeclampsia (non-PE) (n = 32) pregnancies were analyzed using the UHPLC-QE-MS metabolomic platform. Different metabolites were screened, and pathway analysis was conducted. A subgroup analysis was performed among PE neonates to compare the metabolome between appropriate-for-gestational-age infants (n = 21) and small-for-gestational-age (SGA) infants (n = 8). A total of 159 different metabolites were detected in PE and non-PE neonates. Creatinine, N4-acetylcytidine, sphingomyelin (D18:1/16:0), pseudouridine, uric acid, and indolelactic acid were the most significant differential metabolites in the cord serum of PE neonates. Differential metabolite levels were elevated in PE neonates and were involved in the following metabolic pathways: glycine, serine, and threonine metabolism; sphingolipid, glyoxylate, and dicarboxylate metabolism; and arginine biosynthesis. In PE neonates, SGA neonates showed increased levels of hexacosanoyl carnitine and decreased abundance of 3-hydroxybutyric acid and 3-sulfinoalanine. Taurine-related metabolism and ketone body-related pathways were mainly affected. Based on the UHPLC-QE-MS metabolomics analysis, we identified the metabolic profiles of PE and SGA neonates. The abundance of metabolites related to certain amino acid, sphingolipid, and energy metabolism increased in the umbilical cord serum of PE neonates.
Introduction
Preeclampsia (PE) is a major multisystem complication in pregnancy, characterized clinically by new hypertension after 20 weeks of gestation with proteinuria or other clinical signs of impaired end-organ function, including abnormalities of the kidneys, liver, brain, and platelets (1). PE is one of the main causes of iatrogenic premature births (2) and perinatal deaths worldwide, accounting for approximately 10–20% of perinatal mortality (3). It is commonly accepted that the pathophysiology of PE is attributed to the abnormal formation of spiral arteries in the maternal placenta, giving rise to placental oxidative stress and leading to inappropriate and exaggerated maternal responses involving endothelial dysfunction and systemic inflammation (4). However, several studies have shown that preeclampsia has adverse effects on preterm birth and fetal growth restriction in the perinatal period and has long-term effects on offspring in adulthood, manifested as conditions such as an increased risk of hypertension, cardiovascular disease, neurological diseases, and Alzheimer’s disease (5).
Metabolomics, a newly developed subject following genomics, transcriptomics, and proteomics, is an important component of systems biology. Metabolomics refers to the quantitative measurement of metabolic changes caused by pathophysiological changes and the analysis of changes in all endogenous molecule substances at a certain moment to determine the overall state of the organism. Metabolites hold huge potential sources of biomarkers and can also reveal the metabolic process of diseases, assisting in exploring their pathogenesis.
Several studies have demonstrated the metabolome perturbation of PE and discovered characteristic metabolites in the blood, urine, and placenta of pregnant women with PE, including lipids (6), fatty acids (7, 8), metabolites associated with lipid transport (9, 10), amino acids and related metabolites (11, 12), and purine-related metabolites (13). These studies have mainly focused on pregnant women with preeclampsia; however, there are relatively few studies on newborns. A previous study reported an increased level of certain metabolites in the cord blood of PE neonates, including the urea cycle and carnitine synthesis (14). A tryptophan-targeted metabolomics study showed changes in tryptophan metabolites in PE fetal plasma, and 3-hydroxyanthranilic acid was only present in fetal circulation and not in the maternal circulation (15). Youssef et al. (16) found that the lipidomics of PE fetuses was altered, manifested as increased triglycerides, cholesterol, and lipoprotein; however, the mechanism by which maternal PE affects neonatal metabolism is not currently well understood, and the metabolic profile of neonates with preeclampsia remains to be refined.
In this study, we conducted a non-targeted metabolomic analysis of umbilical cord serum of PE and non-PE pregnancies to elucidate the metabolic changes in PE neonates. Moreover, the metabolome of small-for-gestational-age (SGA) infants was compared with that of appropriate-for-gestational-age (AGA) infants in PE neonates.
Materials and Methods
Study Population
This case-control study was approved by the Institutional Review Board of the Peking Union Medical College Hospital. All women provided written informed consent. Cord serum samples were collected from women whose pregnancies were complicated by PE (n = 29) and those without PE (n = 32). All participants were singleton pregnancies that delivered their babies at Peking Union Medical College Hospital, China, between January 2020 and August 2021.
PE was diagnosed according to the criteria published in Practice Bulletin No. 202 of the American College of Obstetricians and Gynecologists in 2020. Women in the PE group had the following characteristics after 20 weeks of gestation: (1) Systolic blood pressure (SBP) ≥ 140 mmHg and/or diastolic blood pressure (DBP) ≥ 90 mmHg on two occasions at least 4 h apart; (2) proteinuria ≥ 300 mg/24 h or severe features (thrombocytopenia, impaired liver function, persistent epigastric pain, renal insufficiency, pulmonary edema, new-onset headache, visual disturbances). Given that neonatal metabolic status is related to gestational age, 10 preterm non-PE neonates were selected for the non-PE group to match the gestational age of the two groups. The cause of preterm birth in these 10 cases was a breech presentation complicated by premature rupture of membranes.
Women with fetal abnormalities, chronic kidney disease, pre-pregnancy diabetes, immune system disease, or other chronic medical diseases were excluded.
Sample Collection
Cord blood samples were obtained right immediately after delivery and were transferred to the laboratory. Samples were centrifuged at 3,000 rpm for 10 min at 4°C. The supernatant serum was extracted, divided into 100 μL aliquots, and stored at –80°C.
Metabolites Extraction
The 50 μL samples were transferred to EP tubes. After adding 200 μL of extract solution (acetonitrile:methanol = 1:1, containing isotopically-labeled internal standard mixture), the samples were vortexed for 30 s, sonicated for 10 min in an ice-water bath, and incubated for 1 h at –40°C to precipitate proteins. The samples were centrifuged at 12,000 rpm for 15 min at 4°C. The resulting supernatant was transferred into a fresh glass vial for further analysis. A quality control (QC) sample was prepared by mixing an equal aliquot of the supernatant from all samples.
UHPLC-QE-MS Analysis
LC-MS/MS analyses were performed using a UHPLC system (Vanquish, Thermo Fisher Scientific) with a UPLC BEH amide column (2.1 mm × 100 mm, 1.7 μm) coupled to a Q Exactive HFX mass spectrometer (Orbitrap MS, Thermo). The mobile phase consisted of 25 mmol/L ammonium acetate and 25 mmol/L ammonia hydroxide in water (pH = 9.75) (A) and acetonitrile (B). The auto-sampler temperature was 4°C, and the injection volume was 3 μL.
The QE HFX mass spectrometer was used for continuously evaluation of the full-scan MS spectrum. The ESI conditions were set as follows: sheath gas flow rate, 30 Arb; Aux gas flow rate, 25 Arb; capillary temperature, 350°C; full MS resolution as 60,000, MS/MS resolution as 7,500, collision energy, 10/30/60 in NCE mode; spray voltage as 3.6 kV (UHPLC-QE-MS(+), positive) or –3.2 kV (UHPLC-QE-MS(–), negative), respectively.
Raw data were converted to the mzXML format and processed with an in-house program developed using R and based on XCMS for peak detection, extraction, alignment, and integration. The following steps were taken: Filtering, removing, replacing, normalizing the features.
Statistical Analysis
The metabolomic data were analyzed by univariate analysis (UVA) and multivariate analysis (MVA) on the UHPLC-QE-MS (+) and UHPLC-QE-MS (–) platforms. For MVA, the dataset was imported into SIMCA (Version 14.1, Umetrics, Umea, Sweden), and the principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (OPLS-DA) models were created. The PCA model was used to analyze the overall distribution of each sample, and the OPLS-DA model was used to analyze the differences in metabolomics between the two groups. For the OPLS-DA model, cross-validation was performed to evaluate the model to avoid bias, and a permutation test was performed to evaluate whether the model was overfitting. Based on the OPLS-DA model, variable importance to projection (VIP) values are reported for each variable (17).
Student’s t-test and Mann-Whitney U-test were applied for UVA to all metabolites, and the false discovery rate (FDR) was performed to adjust the p-value (18). Differential metabolites were selected using MVA and UVA: VIP > 1 and p < 0.05 (adjusted p-value). The metabolites were annotated using the Human Metabolome Database (HMDB)1 (19) and online Kyoto Encyclopedia of Genes and Genomes (KEGG)2 (20) database. MetaboAnalyst 5.03 (21, 22) was used for metabolomic data analysis, including hierarchical cluster analysis and heat map visualization.
Based on the summarized differential metabolites screened by the positive and negative modes, MetaboAnalyst (version 5.0) was used for comprehensive pathway analysis (including enrichment and topology analyses). The enrichment of metabolites in the pathways and the impact factors were analyzed. To select metabolic pathways with a high correlation with preeclampsia, the screening criteria were p < 0.05 or impact > 0.2.
Figures were created using GraphPad (version 8.0.2) and SIMCA (Version 14.1).
Results
Baseline Maternal and Neonatal Characteristics
The characteristics of mothers and newborns are shown in Table 1. Between the PE and non-PE groups, there were no significant differences in maternal age, gravidity, parity, the proportion of assisted reproductive pregnancy, and gestational diabetes mellitus (GDM). There was no significant difference in gestational age at delivery, premature birth rate, cesarean section rate, neonatal Apgar score, or incidence of neonatal complications between the two groups. The SBP, DBP, and body mass index (BMI) of the PE group were significantly higher than those of the non-PE group. The birth weight of the PE neonates was significantly lower than that of the non-PE neonates, and the incidence of SGA was also significantly higher in the PE group.
Umbilical Cord Serum Metabolomics Analysis in Non-preeclampsia and Preeclampsia Group
After data processing, 4,698 and 3,621 metabolic characteristics were detected on the positive [UHPLC-QE-MS (+)] and negative [UHPLC-QE-MS(–)] platforms, respectively. PCA of data from the two platforms showed that PE and non-PE neonatal metabolomic profiles were separated to some extent (Figure 1). OPLS-DA analysis on both platforms showed that the metabolomic profile could significantly distinguish PE and non-PE groups (Figure 2). In the positive and negative platforms, the P-values of CV-ANOVA were 1.8 × 10–5 and 6.9 × 10–9, respectively, indicating good models. Permutation tests on the two platforms showed no overfitting in the OPLS-DA models. Combined MVA and UVA analysis (VIP > 1 and p < 0.05), 1,298 and 838 metabolic characteristics were found to be different in neonatal cord blood on UHPLC-QE-MS (+) and UHPLC-QE-MS (–) platforms, respectively. Among these, 113 and 46 differential metabolites were annotated based on several databases on the UHPLC-QE-MS (+) and UHPLC-QE-MS (–) platforms, respectively.
Figure 1. Principle component analysis (PCA) plots show some separation of metabolic profiles between PE (in blue) and non-PE (in green) for (a)UHPLC-QE-MS(+) and (b)UHPLC-QE-MS(–). (A) POS [UHPLC-QE-MS(+)], 3 components model:R2X = 0.655, Q2 = 0.537. (B) NEG [UHPLC-QE-MS(–)], 9 components model:R2X = 0.503, Q2 = 0.141.
Figure 2. Orthogonal partial least squares discriminant analysis (OPLS-DA) plots show significant separation of metabolic profiles between PE (in blue) and non-PE (in green) for (a1,b1). (A) UHPLC-QE-MS (+) OPLS-DA model, R2X = 0.209, R2Y = 0.78, Q2 = 0.441. (B) UHPLC-QE-MS(–) OPLS-DA model, R2X = 0.175, R2Y = 0.805, Q2 = 0.544. (C) A permutation test of UHPLC-QE-MS (+) model. The Y-axis intercepts were: R2Y (0, 0.6), Q2(0, –0.89). (D) A permutation test of UHPLC-QE-MS(–) model. The Y-axis intercepts were: R2Y (0, 0.7), Q2 (0, –1.01).
Hierarchical cluster analysis was performed on the screened differential metabolites. The heatmap (Figure 3) showed that differential metabolites could be clustered into two types, enriched in umbilical cord serum of PE and non-PE groups. In the positive platform, 106 metabolites increased in PE (including 59 lipid-related metabolites, 16 amino acid-related metabolites, three nucleotide metabolites, and others), and seven decreased in PE (including three amino acid-related metabolites and others). Among these, creatinine, N4-acetylcytidine, and sphingomyelin (D18:1/16:0) showed significant differences and were enriched in the PE group. In the negative platform, 42 metabolites increased (including 13 amino acid-related metabolites, seven amino acid-related metabolites, four nucleotide metabolites, and others), and four metabolites decreased (including one lipid metabolite, one carbohydrate metabolite, and others). Among these, pseudouridine, uric acid, and indolelactic acid showed significant differences and were enriched in the PE group. Important differential metabolites are listed in Table 2.
Figure 3. The heatmap showing the hierarchical clustering characteristics of metabolites of PE and non-PE neonates. The color blue indicates decreasing expression, and red indicates increasing expression. The color intensity reflects the corresponding abundance difference. (A) Differential metabolite characteristics of the UHPLC-QE-MS (+) model. (B) Differential metabolite characteristics of the UHPLC-QE-MS (–) model.
Pathway Analysis of Preeclampsia Neonates
To explore the biological functions of differential metabolites, 159 differential metabolites from both UHPLC-QE-MS (+) and UHPLC-QE-MS (–) platforms were integrated into pathway analysis. The 10 pathways with the highest correlation are listed in Table 3. The topological impact factors and enrichment analysis p-values of the corresponding pathways are shown in the bubble graph (Figure 4). Screening by the criteria: p < 0.05 or impact > 0.2, four metabolic pathways with high influence on PE were selected: glycine, serine, and threonine metabolism, sphingolipid metabolism, glyoxylate and dicarboxylate metabolism, and arginine biosynthesis. The abundance of 10 metabolites in the four pathways above increased in the PE group (Figure 5).
Figure 4. Metabolic pathway analysis results related to differential metabolites between PE and non-PE neonates, presented by bubble plots. Each bubble represents a metabolic pathway. The x-axis of the bubble and its size represent the impact factors of the path in topology analysis. The y-axis of the bubble and its color represent the p-value of enrichment analysis [represented by the negative natural logarithm, -ln(P)]. The color intensity reflects the corresponding p-value and the enrichment degree.
Figure 5. Abundance of differential metabolites in four metabolic pathways related to PE neonates. (A) Sphingolipid metabolism. (B) Glyoxylate and dicarboxylate metabolism. (C) Glycine, serine, and threonine metabolism. (D) Arginine biosynthesis. Mean standard deviation and p-values are shown in the graph. *represents p < 0.05, **represents p < 0.01, ***represents p < 0.001.
Metabolomics Characteristics of Small-for-Gestational-Age in Preeclampsia Neonates
None of the newborns in the non-PE group developed SGA, and the incidence of SGA in the PE group was 27.6%. The PE group was further divided into two subgroups: appropriate for gestational age (AGA) subgroup of PE neonates (PE-AGA group, n = 21) and SGA subgroup of PE neonates (PE-SGA group, n = 8).
MVA showed that metabolomics could clearly distinguish PE-SGA newborns from AGA newborns on positive and negative ion platforms, as shown in Figure 6. Therefore, UVA and MVA were used to screen for metabolites, and the criteria were as follows: VIP > 1 and p < 0.05. Eighteen and seven differential metabolites were found on the UHPLC-QE-MS (+) and UHPLC-QE-MS(–) platforms, respectively (Figure 7).
Figure 6. OPLS-DA plots show a significant separation of metabolic profiles between PE-SGA (in blue) and PE-AGA (in green) for (a,b). (A) UHPLC-QE-MS(+) OPLS-DA model, R2X = 0.169, R2Y = 0.894. (B) UHPLC-QE-MS(–) OPLS-DA model, R2X = 0.199, R2Y = 0.869.
Figure 7. The volcano map shows the p-values, fold changes, and VIP-values of all metabolic characteristics. The color green indicates decreasing expression, and red indicates increasing expression. The spot size reflects the corresponding VIP-value. The major metabolites are highlighted in the maps. (A) Metabolite characteristics of the UHPLC-QE-MS (+) model. (B) Metabolite characteristics of the UHPLC-QE-MS (–) model.
In the UHPLC-QE-MS (+) platform, 14 metabolites increased in PE-SGA (including eight lipid-related metabolites, two amino acid-related metabolites, and others), and four metabolites decreased in PE-SGA (including two lipid-related metabolites, one amino acid-related metabolite, and others). Hexacosanoyl carnitine was the major significant difference. In the UHPLC-QE-MS (–) platform, three metabolites increased in PE-SGA (including two lipid-related metabolites and one amino acid-related metabolite), and four metabolites decreased in PE-SGA (including one lipid-related metabolite, one amino acid-related metabolite, and one carbohydrate metabolite). 3-hydroxybutyric acid and 3-sulfinoalanine were the most significantly different metabolites.
Metabolic Pathway Analysis of Small-for-Gestational-Age in Preeclampsia Neonates
A total of 25 different metabolites from the positive and negative platforms were integrated. Six metabolites related to SGA were distributed in eight metabolic pathways, analyzed by MetaboAnalyst 5.0 (specific results are shown in Figure 8). Topological analysis showed that taurine and hypotaurine metabolism had the greatest influence. The metabolic pathways with the highest enrichment degree were those of synthesis and degradation of ketone bodies. The abundance of 3-hydroxybutyric acid and 3-sulfinoalanine, involved in the two pathways as mentioned above, was reduced in SGA neonates (Figure 9).
Figure 8. Metabolic pathway analysis results of SGA newborns, presented by bubble graph. Each bubble represents a metabolic pathway. The x-axis of the bubble and its size represent the impact factors of the pathway in the topology analysis. The y-axis of the bubble and its color represent the p-value of enrichment analysis [represented by the negative natural logarithm, -ln(P)]. The color intensity reflects the corresponding p-value and the enrichment degree.
Figure 9. The decreased abundance of differential metabolites in two important metabolic pathways related to PE-SGA neonates. (A) Synthesis and degradation of ketone bodies. (B) Taurine and hypotaurine metabolism. Mean, standard deviation, and p-values are shown in the graph. *represents p < 0.05.
Discussion
Preeclampsia is a complex pregnancy disorder that causes damage to fetal growth and development. PE neonates have an average of 5% lower birth weight as compared to children born after an uncomplicated pregnancy, and PE is a significant contributor to preterm birth (23). Preterm birth is the world’s leading cause of neonatal morbidity and mortality. As for short term outcomes, PE is associated with higher rates of infant respiratory distress syndrome, intraventricular hemorrhage, sepsis, bronchopulmonary dysplasia, and neurodevelopmental disability in childhood (24). For long term outcomes, due to abnormal placentation, there is elevated cardiovascular risk in offspring after intrauterine exposure of preeclampsia. Most studies found in offspring exposed to intrauterine preeclampsia, at the age of 9–17 years old, higher systolic blood pressure, other found higher diastolic blood pressure and some both (25). Therefore, it is of great significance to explore the metabolic profiles of neonates with PE and SGA neonates. Based on the UHPLC-QE-MS platform, our study observed metabolic alterations in cord serum involved in the metabolism of amino acids, lipids, and nucleosides in neonates with PE.
Firstly we investigated the metabolites associated with PE newborns, including 6 metabolites that have highest VIP value. We found out that creatinine and uric acid were significantly higher in PE newborns which is consistent with previously reported (14). Moreover, these metabolites also have long been considered to be significantly higher in pregnant women with PE (11). As mentioned before, it is bothering that cord blood might reflect on maternal values and not only the newborns. Creatinine and uric acid were generally considered to be associated with kidney function. It was reported that neonatal creatinine could reflects the level of maternal renal function (26), and predict the onset of neonatal acute kidney injury to some extent (27). Uric acid have effects on fetal growth and kidney development, and lead to a tendency of low-birth-weight infants (28). High levels of uric acid can lead to hypertension-like pathological changes such as arteriole thickening, arteriosclerosis and vasoconstriction in rats, reflecting the correlation between uric acid and hypertension (29). The possible molecular mechanism might due to uric acid can hinder vascular endothelial cell proliferation and function (30). Therefore, the offspring of PE tend to have an increased risk of hypertension in the long term. It may be partly related to the vascular changes resulting from exposure to high levels creatinine and uric acid in fetal and neonatal periods. However, the mechanism of long-term impact in adulthood is unclear.
N4-acetylcytidine and pseudouridine also increased in our PE newborns. Previous studies have found that these metabolites were elevated in the blood of patients with chronic kidney disease, uremia, and pediatric acute kidney injury, respectively (31, 32). These metabolites, cleared by the kidney (33), are related to the decrease of glomerular filtration rate (34) and may reflect the renal function of patients to some extent.
Indoleolactic acid is a tryptophan metabolite of an essential amino acid during pregnancy that meets the needs of fetal growth and development (35). Consistent with our study, Jääskeläinen et al. also found that indoleolactic acid increased in umbilical cord blood of PE newborns (14). In addition, Morita et al. found that PE can lead to an elevated concentration of tryptophan metabolites levels, including indoleolactic acid, in neonatal serum (15). Previous studies have shown that phenyllactate dehydrogenase and acyl-CoA dehydrogenase in the placenta can convert indoleolactic acid to indolepropionic acid and that PE may inhibit the activity of both enzymes, leading to the accumulation of indoleolactic acid (36). In addition, indoleolactic acid has been found to help maintain cell stability in human umbilical vein endothelial cells, antagonizing the abnormal endothelial function of PE. Besides, indoleolactic acid show an anti-inflammatory effect in the intestine of premature infants, and its low level may be implicated in irritable bowel syndrome and NEC (37). Therefore, increased indoleolactic acid seems to be a self-protective factor for newborns with PE.
Further, we use metabolites to conduct pathway analysis and PE newborns were associated with the following pathways.
(1) Glycine, serine, and threonine metabolism. Glycine is located at the central position in the glycine, serine, and threonine pathways. It can also be interconverted with l-serine and l-threonine. Sarcosine is metabolized to glycine by sarcosine dehydrogenase. A latest study of neonates also found elevated glycine in SGA neonates (38). From the perspective of nutrition, it is speculated that this amino acid alteration, we found, may reflect the enhanced decomposition of protein to make up for the insufficient energy supply of PE placenta. Glycine may alleviate inflammatory injury induced by ischemia and reperfusion and inhibit interleukin production (39, 40). It was found that glycine in FGR newborns was relatively increased, and glycine concentration in umbilical cord blood was inversely proportional to oxygen content (41). The change in glycine levels in cord blood may reflect the adaptation of newborns with PE in the ischemia and hypoxia state. It is generally believed that abnormal vasculogenesis in the PE placenta may lead to relative ischemia and hypoxia in the fetal supply. This is the first study to determine the characteristics of the glycine metabolic pathway in PE newborns.
(2) Arginine biosynthesis. Our data showed that arginine biosynthesis was altered in PE neonates due to L-citrulline and L-glutamine’s increased abundance. Glutamine is a precursor of citrulline, synthesized in the gut. Citrulline is mainly converted to arginine in the kidneys in adults (42). In preterm neonates, citrulline maybe directly converted in the guts due to imperfect gut-kidney organ circulation (43). A previous study found that citrulline was slightly higher in PE cord blood than in non-PE cord blood (44). Some studies have shown that citrulline antagonizes hypertension and reduces arterial tension by activating the citrulline-nitric oxide cycle (45). However, a recent systematic review showed no strong evidence of citrulline in fighting hypertension, and more clinical studies are needed (46). In addition, glutamine is a potential marker of brain injury. In previous studies, among newborns with hypoxic-ischemic brains, those with poor outcomes had higher glutamine concentrations (47). Glutamine is also associated with severe brain edema (48). Therefore, glutamine may help in the early recognition and treatment of neonatal PE brain injury.
(3) Sphingolipid metabolism. N-Acylsphingosine, a ceramide, is a sphingolipid with antiangiogenic properties. Ceramide, sphingomyelin, and sphinganine play important roles in the pathogenesis of cardiovascular diseases and regulation of endothelial cell function (49). A study in newborns found an increase in ceramide in feces before the onset of NEC, and abnormal sphingolipid metabolism may change the permeability of intestinal cells, leading to inflammatory diseases (50). Ceramide also promotes apoptosis (51). The ceramide concentration is often higher in patients with Alzheimer’s disease, leading to a wide range of neuronal apoptosis and cognitive impairments (52). However, the role of ceramide in neonatal cognition and brain development remains unclear.
PE tend to have a high incidence of low birth weight and SGA (53). Our study also reflected this feature. As mentioned above, uric acid is an important risk factor for low birth weight infants (28). The free activity of uric acid in maternal-fetal circulation can hinder the proliferation and function of endothelial cells, causing interference in the fetal growth and organ development (30, 54). Another neonatal metabolomics study has found a correlation between aberrant sphingolipids and SGA (55). Growth restriction is associated with placental cell apoptosis, which may be regulated by sphingolipids (56, 57). The abnormal sphingolipid metabolism, observed in our study, may increase the signaling of apoptosis, limiting the fetal growth. On the other hand, it is commonly accepted that the main pathophysiological mechanism of SGA is the reduced uterine-placenta blood perfusion. Poor formation of spiral arteries, associated with PE, can also lead to the similar damage, increasing the incidence of fetal growth restriction (58).
In order to explore the potential molecular mechanism underling PE associated adverse neonatal outcomes, we investigated the metabolites associated with SGA neonates. We observed abnormalities in amino acid, lipid, and ketone body metabolisms in SGA neonates. Carnitine is a carrier of lipid transport, and lipids move from the cytoplasm to the mitochondria, where they are metabolized. Mitochondrial dysfunction makes it difficult to consume fat, resulting in carnitine accumulation (59). Beken et al. also found abnormal carnitine in SGA neonates, including reduced propionyl carnitine and methylglutaryl carnitine (38). Another study reported an increased carnitine concentration in neonates with a high body fat rate (60). Combining our study with other studies, we speculate that SGA neonates have abnormal lipid metabolism and fat accumulation. Individuals born SGA are at high risk of long-term obesity (61). Although long-term growth is affected by many factors, such as the environment and diet, metabolic abnormalities originating from newborns may be one of the reasons.
Previous studies have indicated that fetuses exhibit growing energy requirements during fetal growth in normal pregnancy, including increased glucose and fatty acids consumption. Fetal ketogenesis is a response to the potential risk of hypoglycemia (62) and is also an important backup supplement for fetal energy supply, resulting in the accumulation of acetoacetic and hydroxybutyric acids (63). A positive correlation between birth weight and hydroxybutyric acid concentration in umbilical cord blood was also observed in studies by Mansel et al. (64) and Lowe et al. (65). In this study, the low level of 3-hydroxybutyric acid indicates that the metabolic pathway of ketone bodies in SGA neonates was impaired. Inadequate energy supply and placental nutritional disorders can result in adaptive changes in newborns that manifest as weight loss.
In the taurine metabolic pathway, 3-sulfinoalanine produces taurine. Taurine metabolism plays an important role in facilitating substance exchange and nutrient supply in the placenta. One study found a positive association between maternal taurine intake and neonatal height in pregnancy (66). Previous studies have shown that low taurine levels can limit neonatal neurodevelopment, retinal development, and intestinal absorption function (67). Animal experiments have shown that taurine supplementation may have a protective effect on neural development in neonatal rats (68). However, little research has been conducted on human neonates due to limitations of ethical requirements. The latest meta-analysis of taurine supplementation did not confirm its role and showed no significant effects on the growth and development of newborns (69). On the other hand, taurine is also an important biomarker of individual neurodegeneration (70). Due to its antioxidant property, high levels of taurine can reduce cognitive impairment caused by neurodegeneration lesions (71, 72). It is observed, children with SGA have a relatively elevated risk of cognitive impairment, learning difficulties, and even cerebral palsy in the long term (73). Low abundance of taurine in our study may partially explain this problem. Therefore, detection of taurine level may be beneficial for timely evaluation and prevention of neurological diseases, not only for newborns, but also for offspring.
The results of our study reflect the value of metabolomics in evaluating the pathological mechanisms of neonates with PE. These metabolites have diverse functions. Not only can a single metabolite reflect the inherent characteristics of the disease, but a variety of related metabolites can also better reflect the changes in a certain biochemical mechanism. These findings are insufficient to fully reflect all the changes involved in the pathogenesis of PE neonates; however, they still provide a basis and perspective for further research on PE neonates. This is the first metabolomic study of SGA neonates with PE, although it was not the primary objective of our study. This study recorded and statistically analyzed metabolic characteristics and recent complications of PE neonates. However, the long-term effects of metabolic characteristics on neonates remain unclear. Further studies on the risk of related metabolic or cardiovascular diseases based on metabolomic profiles may be needed in the future.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by the Institutional Review Board of Peking Union Medical College Hospital. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author Contributions
XW, JL, and YS conceived and designed the study. XW and JL collected the serum materials and wrote the manuscript. JL and XH performed the experiments. XW, JL, and XH analyzed and modified the data. YS provided advice, assistance, review and editing. All authors read and approved the manuscript.
Funding
Funding was provided by the National Natural Science Foundation of China (Nos. 81871191 and 81900723).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank all participants and their attending obstetricians, nurses for their support and selfless dedication.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fped.2022.869381/full#supplementary-material
Footnotes
References
1. Lyall F, Belfort M. Pre-Eclampsia: Etiology and Clinical Practice. Cambridge: Cambridge University Press (2011).
2. Basso O, Rasmussen S, Weinberg CR, Wilcox AJ, Irgens LM, Skjaerven R. Trends in fetal and infant survival following preeclampsia. JAMA. (2006) 296:1357–62. doi: 10.1001/jama.296.11.1357
3. Lawn JE, Blencowe H, Pattinson R, Cousens S, Kumar R, Ibiebele I, et al. Stillbirths: Where? When? Why? How to make the data count? Lancet. (2011) 377:1448–63. doi: 10.1016/S0140-6736(10)62187-3
5. Lu HQ, Hu R. Lasting effects of intrauterine exposure to preeclampsia on offspring and the underlying mechanism. AJP Rep. (2019) 9:e275–91. doi: 10.1055/s-0039-1695004
6. Austdal M, Skråstad RB, Gundersen AS, Austgulen R, Iversen AC, Bathen TF. Metabolomic biomarkers in serum and urine in women with preeclampsia. PLoS One. (2014) 9:e91923. doi: 10.1371/journal.pone.0091923
7. Liu Y, Zu L, Cai W, Cheng Q, Hua T, Peng L, et al. Metabolomics revealed decreased level of omega-3 PUFA-derived protective eicosanoids in pregnant women with pre-eclampsia. Clin Exp Pharmacol Physiol. (2019) 46:705–10. doi: 10.1111/1440-1681.13095
8. Jain S, Jayasimhulu K, Clark JF. Metabolomic analysis of molecular species of phospholipids from normotensive and preeclamptic human placenta electrospray ionization mass spectrometry. Front Biosci J Virtual Library. (2004) 9:3167–75. doi: 10.2741/1470
9. Tan B, Ma Y, Zhang L, Li N, Zhang J. The application of metabolomics analysis in the research of gestational diabetes mellitus and preeclampsia. J Obstet Gynaecol Res. (2020) 46:1310–8. doi: 10.1111/jog.14304
10. Chen T, He P, Tan Y, Xu D. Biomarker identification and pathway analysis of preeclampsia based on serum metabolomics. Biochem Biophys Res Commun. (2017) 485:119–25. doi: 10.1016/j.bbrc.2017.02.032
11. Kenny LC, Broadhurst D, Brown M, Dunn WB, Redman CW, Kell DB, et al. Detection and identification of novel metabolomic biomarkers in preeclampsia. Reprod Sci. (2008) 15:591–7. doi: 10.1177/1933719108316908
12. Austdal M, Thomsen LCV, Tangerås LH, Skei B, Mathew S, Bjørge L, et al. Metabolic profiles of placenta in preeclampsia using HR-MAS MRS metabolomics. Placenta. (2015) 36:1455–62. doi: 10.1016/j.placenta.2015.10.01
13. Senyavina NV, Khaustova SA, Grebennik TK, Pavlovich SV. Analysis of purine metabolites maternal serum for evaluating the risk of gestosis. Bull Exp Biol Med. (2013) 155:682–4. doi: 10.1007/s10517-013-2225-y
14. Jääskeläinen T, Kärkkäinen O, Jokkala J, Litonius K, Heinonen S, Auriola S, et al. A non-targeted LC-MS profiling reveals elevated levels of carnitine precursors and trimethylated compounds in the cord plasma of pre-eclamptic infants. Sci Rep. (2018) 8:1–2. doi: 10.1038/s41598-018-32804-5
15. Morita I, Kawamoto M, Yoshida H. Difference in the concentration of tryptophan metabolites between maternal and umbilical foetal blood. J Chromatogr. (1992) 576:334–9. doi: 10.1016/0378-434792)80208-8
16. Youssef L, Simoes RV, Miranda J, Garcia-Martín ML, Paules C, Crovetto F, et al. Paired maternal and fetal metabolomics reveal a differential fingerprint in preeclampsia versus fetal growth restriction. Sci Rep. (2021) 111:14422. doi: 10.1038/s41598-021-93936-9
17. Jansson J, Willing B, Lucio M, Fekete A, Dicksved J, Halfvarson J, et al. Metabolomics reveals metabolic biomarkers of Crohn’s disease. PLoS One. (2009) 4:e6386. doi: 10.1371/journal.pone.0006386
18. Zhu D, Hero AO, Qin ZS, Swaroop A. High throughput screening of co-expressed gene pairs with controlled false discovery rate (FDR) and minimum acceptable strength (MAS). J Comput Biol. (2005) 12:1029–45. doi: 10.1089/cmb.2005.12.1029
19. Wishart DS, Knox C, Guo AC, Eisner R, Young N, Gautam B, et al. HMDB: a knowledgebase for the human metabolome. Nucleic Acids Res. (2009) 37:D603–10. doi: 10.1093/nar/gkn810
20. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. (2000) 28:27–30. doi: 10.1093/nar/28.1.27
21. Pang Z, Chong J, Zhou G, de Lima Morais DA, Chang L, Barrette M, et al. MetaboAnalyst 5.0: narrowing the gap between raw spectra and functional insights. Nucleic Acids Res. (2021) 49:W388–96. doi: 10.1093/nar/gkab382
22. Chong J, Wishart DS, Xia J. Using MetaboAnalyst 4.0 for comprehensive and integrative metabolomics data analysis. Curr Protoc Bioinformatics. (2019) 68:e86. doi: 10.1002/cpbi.86
23. Bokslag A, van Weissenbruch M, Mol BW, de Groot CJ. Preeclampsia; short and long-term consequences for mother and neonate. Early Hum Dev. (2016) 102:47–50. doi: 10.1016/j.earlhumdev.2016.09.007
24. Goldenberg RL, Culhane JF, Iams JD, Romero R. Epidemiology and causes of preterm birth. Lancet. (2008) 371:75–84. doi: 10.1016/S0140-6736(08)60074-4
25. Saigal S, Doyle LW. An overview of mortality and sequelae of preterm birth from infancy to adulthood. Lancet. (2008) 371:261–9. doi: 10.1016/S0140-673608)60136-1
26. Askenazi DJ, Ambalavanan N, Goldstein SL. Acute kidney injury in critically ill newborns: what do we know? What do we need to learn? Pediatr Nephrol. (2009) 24:265–74. doi: 10.1007/s00467-008-1060-2
27. Allegaert K, Smits A, Mekahli D, van den Anker JN. Creatinine at birth correlates with gestational age and birth weight: another factor of the imbroglio in early neonatal life. Neonatology. (2020) 117:637–40. doi: 10.1159/000510658
28. Sagen N, Haram K, Nilsen ST. Serum urate as a predictor of fetal outcome in severe preeclampsia. Acta Obstet Gynecol Scand. (1984) 63:71–5. doi: 10.3109/00016348409156277
29. Sanchez-Lozada LG, Tapia E, Avila-Casado C, Soto V, Franco M, Santamaría J, et al. Mild hyperuricemia induces glomerular hypertension in normal rats. Am J Physiol Renal Physiol. (2002) 283:F1105–10. doi: 10.1152/ajprenal.00170.2002
30. Feig DI, Nakagawa T, Karumanchi SA, Oliver WJ, Kang DH, Finch J, et al. Hypothesis: uric acid, nephron number, and the pathogenesis of essential hypertension. Kidney Int. (2004) 66:281–7. doi: 10.1111/j.1523-1755.2004.00729.x
31. Wang S, Xiao C, Liu C, Li J, Fang F, Lu X, et al. Identification of biomarkers of sepsis-associated acute kidney injury in pediatric patients based on UPLC-QTOF/MS. Inflammation. (2020) 43:629–40. doi: 10.1007/s10753-019-01144-5
32. Sekula P, Goek O-N, Quaye L, Barrios C, Levey AS, Romisch-Margl W, et al. A metabolome-wide association study of kidney function and disease in the general population. J Am Soc Nephrol. (2016) 274:1175–88. doi: 10.1681/ASN.2014111099
33. Dzúrik R, Lajdová I, Spustová V, Opatrny K. Pseudouridine excretion in healthy subjects and its accumulation in renal failure. Nephron. (1992) 61:64–7. doi: 10.1159/000186836
34. Niewczas MA, Mathew AV, Croall S, Byun J, Major M, Sabisetti VS, et al. Circulating modified metabolites and a risk of ESRD in patients with type 1 diabetes and chronic kidney disease. Diabetes Care. (2017) 40:383–90. doi: 10.2337/dc16-0173
35. Badawy AA. Tryptophan metabolism disposition and utilization in pregnancy. Biosci Rep. (2015) 35:e00261. doi: 10.1042/BSR20150197
36. Agus A, Planchais J, Sokol H. Gut microbiota regulation of tryptophan metabolism in health and disease. Cell Host Microbe. (2018) 23:716–24. doi: 10.1016/j.chom.2018.05.003
37. Meng D, Sommella E, Salviati E, Campiglia P, Ganguli K, Djebali K, et al. Indole-3-lactic acid, a metabolite of tryptophan, secreted by Bifidobacterium longum subspecies infantis is anti-inflammatory in the immature intestine. Pediatr Res. (2020) 88:209–17. doi: 10.1038/s41390-019-0740-x
38. Beken S, Abali S, Yildirim Saral N, Guner B, Dinc T, Albayrak E, et al. Early postnatal metabolic profile in neonates with different birth weight status: a pilot study. Front Pediatr. (2021) 9:646860. doi: 10.3389/fped.2021.646860
39. Liu A, Fang H, Dirsch O, Jin H, Dahmen U. Early release of macrophage migration inhibitory factor after liver ischemia and reperfusion injury in rats. Cytokine. (2012) 57:150–7. doi: 10.1016/j.cyto.2011.11.009
40. Ito K, Ozasa H, Noda Y, Koike Y, Arii S, Horikawa S. Effect of non-essential amino acid glycine administration on the liver regeneration of partially hepatectomized rats with hepatic ischemia/reperfusion injury. Clin Nutr. (2008) 27:773–80. doi: 10.1016/j.clnu.2008.06.012
41. Economides DL, Nicolaides KH, Gahl WA, Bernardini I, Bottoms S, Evans M. Cordocentesis in the diagnosis of intrauterine starvation. Am J Obstet Gynecol. (1989) 161:1004–8. doi: 10.1016/0002-9378(89)90772-2
42. Wijnands KA, Castermans TM, Hommen MP, Meesters DM, Poeze M. Arginine and citrulline and the immune response in sepsis. Nutrients. (2015) 7:1426–63. doi: 10.3390/nu7031426
43. Bourdon A, Rougé C, Legrand A, Robert CD, Piloquet H, Vodovar M, et al. Urinary citrulline in very low birth weight preterm infants receiving intravenous nutrition. Br. J. Nutr. (2011) 108:1150–4. doi: 10.1017/S0007114511006660
44. Alacam H, Dikmen ZG, Yaman H, Cakir E, Deren O, Akgul EO, et al. The role of asymmetric dimethyl arginine and oxidant/antioxidant system in preeclampsia. Fetal Pediatr Pathol. (2011) 30:387–93. doi: 10.3109/15513815.2011.618869
45. Mirenayat MS, Moradi S, Mohammadi H, Rouhani MH. Effect of L-Citrulline supplementation on blood pressure: a systematic review and meta-analysis of clinical trials. Curr Hypertens Rep. (2018) 20:98. doi: 10.1007/s11906-018-0898-3
46. Barkhidarian B, Khorshidi M, Shab-Bidar S, Hashemi B. Effects of L-citrulline supplementation on blood pressure: a systematic review and meta-analysis. Avicenna J Phytomed. (2019) 9:10–20.
47. Shibasaki J, Aida N, Morisaki N, Tomiyasu M, Nishi Y, Toyoshima K. Changes in brain metabolite concentrations after neonatal hypoxic-ischemic encephalopathy. Radiology. (2018) 288:840–8. doi: 10.1148/radiol.2018172083
48. Poveda MJ, Bernabeu A, Concepción L, Roa E, de Madaria E, Zapater P, et al. Brain edema dynamics in patients with overt hepatic encephalopathy A magnetic resonance imaging study. Neuroimage. (2010) 522:481–7. doi: 10.1016/j.neuroimage.2010.04.260
49. Venkataraman K, Lee YM, Michaud J, Thangada S, Ai Y, Bonkovsky HL, et al. Vascular endothelium as a contributor of plasma sphingosine 1-phosphate. Circ Res. (2008) 102:669–76.
50. Rusconi B, Jiang X, Sidhu R, Ory DS, Warner BB, Tarr PI. Gut sphingolipid composition as a prelude to necrotizing enterocolitis. Sci Rep. (2018) 81:10984. doi: 10.1038/s41598-018-28862-4
51. Jana A, Hogan EL, Pahan K. Ceramide and neurodegeneration: susceptibility of neurons and oligodendrocytes to cell damage and death. J Neurol Sci. (2009) 278:5–15. doi: 10.1016/j.jns.2008.12.010
52. Loft LMI, Moseholm KF, Pedersen KKW, Jensen MK, Koch M, Cronjé HT. Sphingomyelins and ceramides: possible biomarkers for dementia? Curr Opin Lipidol. (2022) 33:57–67. doi: 10.1097/MOL.0000000000000804
53. Bakker R, Steegers EA, Hofman A, Jaddoe VW. Blood pressure in different gestational trimesters, fetal growth, and the risk of adverse birth outcomes: the generation R study. Am J Epidemiol. (2011) 174:797–806. doi: 10.1093/aje/kwr151
54. Chang FM, Chow SN, Huang HC, Hsieh FJ, Chen HY, Lee TY, et al. The placental transfer and concentration difference in maternal and neonatal serum uric acid at parturition: comparison of normal pregnancies and gestosis. Biol Res Pregnancy Perinatol. (1987) 8:35–9.
55. Horgan RP, Broadhurst DI, Walsh SK, Dunn WB, Brown M, Roberts CT, et al. Metabolic profiling uncovers a phenotypic signature of small for gestational age in early pregnancy. J Proteome Res. (2011) 10:3660–73. doi: 10.1021/pr2002897
56. Smith SC, Baker PN, Symonds EM. Increased placental apoptosis in intrauterine growth restriction. Am J Obstet Gynecol. (1997) 177:1395–401. doi: 10.1016/s0002-9378(97)70081-4
57. Maceyka M, Payne SG, Milstien S, Spiegel S. Sphingosine kinase, sphingosine-1-phosphate, and apoptosis. Biochim Biophys Acta. (2002) 1585:193–201. doi: 10.1016/s1388-1981(02)00341-4
58. Everett TR, Lees CC. Beyond the placental bed: placental and systemic determinants of the uterine artery Doppler waveform. Placenta. (2012) 33:893–901. doi: 10.1016/j.placenta.2012.07.011
59. Zhou X, Han TL, Chen H, Baker PN, Qi H, Zhang H. Impaired mitochondrial fusion, autophagy, biogenesis and dysregulated lipid metabolism is associated with preeclampsia. Exp Cell Res. (2017) 359:195–204. doi: 10.1016/j.yexcr.2017.07.029
60. Kadakia R, Scholtens DM, Rouleau GW, Talbot O, Ilkayeva OR, George T, et al. Cord blood metabolites associated with newborn adiposity and hyperinsulinemia. J Pediatr. (2018) 203:144–9.e1. doi: 10.1016/j.jpeds.2018.07.056
61. Vicens-Calvet E, Espadero RM, Carrascosa A, Spanish Sga Collaborative Group, Small for Gestational Age. Longitudinal study of the pubertal growth spurt in children born small for gestational age without postnatal catch-up growth. J Pediatr Endocrinol Metab. (2002) 15:381–8. doi: 10.1515/jpem.2002.15.4.381
62. Hawdon JM, Ward Platt MP. Metabolic adaptation in small for gestational age infants. Arch Dis Child. (1993) 68:262–8. doi: 10.1136/adc.68.3_Spec_No.262
63. Bon C, Raudrant D, Golfier F, Poloce F, Champion F, Pichot J, et al. Feto-maternal metabolism in human normal pregnancies: study of 73 cases. Ann Biol Clin (Paris). (2007) 65:609–19. doi: 10.1684/abc.2007.0168
64. Mansell T, Vlahos A, Collier F, Ponsonby AL, Vuillermin P, Ellul S, et al. The newborn metabolome: associations with gestational diabetes, sex, gestation, birth mode, and birth weight. Pediatr Res. (2021) 15:1.
65. Lowe WL, Bain JR, Nodzenski M, Reisetter AC, Muehlbauer MJ, Stevens RD, et al. Maternal BMI and glycemia impact the fetal metabolome. Diabetes Care. (2017) 40:902–10. doi: 10.2337/dc16-2452
66. Jung YM, Choi MJ. Relation of taurine intake during pregnancy and newborns’. Growth. Adv Exp Med Biol. (2019) 1155:283–92. doi: 10.1007/978-981-13-8023-5_27
67. Wharton BA, Morley R, Isaacs EB, Cole TJ, Lucas A. Low plasma taurine and later neurodevelopment. Arch Dis Child Fetal Neonatal Ed. (2004) 89:F497–8. doi: 10.1136/adc.2003.048389
68. Zhu XY, Ma PS, Wu W, Zhou R, Hao YJ, Niu Y, et al. Neuroprotective actions of taurine on hypoxic-ischemic brain damage in neonatal rats. Brain Res Bull. (2016) 124:295–305. doi: 10.1016/j.brainresbull.2016.06.010
69. Cao SL, Jiang H, Niu SP, Wang XH, Du S. Effects of taurine supplementation on growth in low birth weight infants: a systematic review and meta-analysis. Indian J Pediatr. (2018) 85:855–60. doi: 10.1007/s12098-018-2609-0
70. Menzie J, Pan C, Prentice H, Wu JY. Taurine and central nervous system disorders. Amino Acids. (2014) 46:31–46. doi: 10.1007/s00726-012-1382-z
71. Javed H, Khan A, Vaibhav K, Moshahid Khan M, Ahmad A, Ejaz Ahmad M, et al. Taurine ameliorates neurobehavioral, neurochemical and immunohistochemical changes in sporadic dementia of Alzheimer’s type (SDAT) caused by intracerebroventricular streptozotocin in rats. Neurol Sci. (2013) 34:2181–92. doi: 10.1007/s10072-013-1444-3
72. Taranukhin AG, Taranukhina EY, Saransaari P, Podkletnova IM, Pelto-Huikko M, Oja SS. Neuroprotection by taurine in ethanol-induced apoptosis in the developing cerebellum. J Biomed Sci. (2010) 17(Suppl. 1):S12. doi: 10.1186/1423-0127-17-S1-S12
Keywords: preeclampsia, neonate, SGA, metabolomics, pathway
Citation: Wang X, Liu J, Hui X and Song Y (2022) Metabolomics Applied to Cord Serum in Preeclampsia Newborns: Implications for Neonatal Outcomes. Front. Pediatr. 10:869381. doi: 10.3389/fped.2022.869381
Received: 04 February 2022; Accepted: 11 March 2022;
Published: 25 April 2022.
Edited by:
Offer Erez, Soroka Medical Center, IsraelReviewed by:
Cristina Loddo, Azienda Ospedaliero-Universitaria Cagliari, ItalyOrna Staretz Chacham, Ben-Gurion University of the Negev, Israel
Copyright © 2022 Wang, Liu, Hui and Song. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yingna Song, ZGVhcnlpbmduYUAxNjMuY29t