Associations of Plasma Fatty Acid Patterns During Pregnancy With Gestational Diabetes Mellitus
Peiyun Li1,2,3 
Shan Hu1,2 
Yalun Zhu1,2 Taoping Sun1,2 Yue Huang1,2 Zihui Xu1,2 Hongjie Liu1,2 Cheng Luo1,2 Shiqiong Zhou3 Aijun Tan4* 
Liegang Liu1,2*
- 1Department of Nutrition and Food Hygiene, Hubei Key Laboratory of Food Nutrition and Safety, School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
- 2The Ministry of Education Key Laboratory of Environment and Health, School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
- 3Department of Clinical Nutrition, Wuhan Children’s Hospital (Wuhan Maternal and Child Healthcare Hospital), Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
- 4Zhuhai Center for Disease Control and Prevention, Zhuhai, China
Background: Limited studies have explored the difference of fatty acid profile between women with and without gestational diabetes mellitus (GDM), and the results were inconsistent. Individual fatty acids tend to be interrelated because of the shared food sources and metabolic pathways. Thus, whether fatty acid patters during pregnancy were related to GDM odds needs further exploration.
Objective: To identify plasma fatty acid patters during pregnancy and their associations with odds of GDM.
Methods: A hospital-based case-control study including 217 GDM cases and 217 matched controls was carried out in urban Wuhan, China from August 2012 to April 2015. All the participants were enrolled at the time of GDM screening and provided fasting blood samples with informed consent. We measured plasma concentrations of fatty acids by gas chromatography–mass spectrometry, and derived potential fatty acid patterns (FAPs) through principal components analysis. Conditional logistic regression and restricted cubic spline model were used to evaluate the associations between individual fatty acids or FAPs and odds of GDM.
Results: Twenty individual fatty acids with relative concentrations ≥0.05% were included in the analyses. Compared with control group, GDM group had significantly higher concentrations of total fatty acids, 24:1n-9, and relatively lower levels of 14:0, 15:0, 17:0, 18:0, 24:0, 16:1n-7, 20:1n-9,18:3n-6, 20:2n-6, 18:3n-3, 20:3n-3, 22:5n-3. Two novel patterns of fatty acids were identified to be associated with lower odds of GDM: (1) relatively higher odd-chain fatty acids, 14:0, 18:0, 18:3n-3, 20:2n-6, 20:3n-6 and lower 24:1n-9 and 18:2n-6 [adjusted odds ratio (OR) (95% confidence interval) (CI) for quartiles 4 vs. 1: 0.42 (0.23–0.76), P-trend = 0.002], (2) relatively higher n-3 polyunsaturated fatty acids, 24:0, 18:3n-6 and lower 16:0 and 20:4n-6 [adjusted OR (95% CI) for quartiles 4 vs. 1: 0.48 (0.26–0.90), P-trend = 0.018].
Conclusion: Our findings suggested that two novel FAPs were inversely associated with GDM odds. The combination of circulating fatty acids could be a more significant marker of GDM development than individual fatty acids or their subgroups.
Introduction
Gestational diabetes mellitus (GDM) refers to hyperglycemia diagnosed for the first time during pregnancy and is one of the most common complications of pregnancy in the world. According to the 10th edition of the Diabetes Atlas published by the International Diabetes Federation, an estimated 16.7% (21.1 million) of live births were affected by hyperglycemia in pregnancy in 2021. Of which, 80.3% were due to GDM (1). Although GDM usually resolves once the delivery ends, its impact on maternal and child health cannot be ignored. Women with GDM during pregnancy are at higher risk of adverse pregnancy outcomes and have an increased risk of developing GDM in subsequent pregnancies or type 2 diabetes later in life (2, 3). Babies exposed to hyperglycemia in utero also have a life-long higher risk of obesity and glucose intolerance (4, 5). The prevention of GDM could have far-reaching effects on the short-term and long-term health of mothers and offsprings. Hence, identifying potentially modifiable risk factors and evaluating their impact on GDM is of high priority.
Recent years, dietary fat intake has been shown to be associated with the development of insulin resistance and diabetes (6, 7). As the important composition of fat, fatty acids in tissues can reflect both the quantity and quality of dietary fat intake without being affected by recall bias and have been considered as reliable biomarkers in epidemiologic studies (8). However, up to date, limited studies have explored the difference of fatty acid profile between women with and without GDM (9). Some of them have suggested that certain fatty acids might be associated with the development of GDM, but the results were inconsistent and partly depends on the fatty acids involved (10–17). Furthermore, all the related studies only focused on traditional groups of fatty acids or a few selected individual fatty acids, whereas many circulating fatty acids of which the concentrations tend to be interrelated because of the shared food sources and metabolic pathways (18). Otherwise, fatty acid composition in previous studies were often expressed as relative amounts like mole percent or weight percent (19), which lead to the fact that a change in the level of one fatty acid might influence the amounts of others. Thus, the potential synergistic and additive effects among different fatty acids may be ignored when separately investigating the potential influence of targeted individual fatty acids or fatty acid groups on GDM risk.
To take the above complexity into account, using statistical pattern-recognition approach might be worth trying. Principal components analysis (PCA), a widely used statistical method of dimension reduction, has been applied to explore novel informative patterns such as dietary patterns or genetic patterns in complex data in the field of epidemiology (20, 21). As to biomarker panels like circulating fatty acids, PCA is also applicable and has a main effect on maximizing the information of both major and minor fatty acids. Several epidemiologic studies have identified associations of circulating or tissue fatty acid patterns (FAPs) with obesity, hypertension, cardiovascular diseases, type 2 diabetes, and cancers (22–26), which indicated that potential clinical and biological implications were existed and over that of individual fatty acids or fatty acid subgroups. However, whether FAPs during pregnancy were related to GDM odds remains unknown. We therefore conducted a matched case-control study among Chinese women to investigate the association of plasma fatty acid subgroups and individual fatty acids during pregnancy with odds of GDM initially, and further applied PCA to identify the novel FAPs associated with GDM.
Materials and Methods
Study Population
The hospital-based case-control study was carried out in urban Wuhan, China from August 2012 to April 2015. Pregnant women who screened for GDM at the outpatient clinics of the Department of Endocrinology, Tongji Hospital were invited to participant in the study. The inclusion criteria were as follows: (1) age ≥20 years; (2) gestational age at GDM screening ≥24 weeks; (3) singleton pregnancy. We excluded women who met any of the following items: history of diabetes (including but not limited to GDM), cardiovascular disease, cancer or other systemic diseases; pharmacologic treatment or dietary supplements use (e.g., fish oil, cod liver oil, pure docosahexaenoic acid supplements, albumen power, etc.) that might influence glucose or lipids metabolism; accompanied by other pregnancy complications; blood sample hemolysis or insufficiency; and incomplete basic information. Fasting blood samples (5 ml) were collected at the time of GDM screening using anticoagulant tubes and centrifuged at 3,000 rpm for 5 min. Plasma were separated from blood cells and stored at −80°C for further assay. The study was approved by the Ethics Committee of Tongji Medical College and was registered at www.clinicaltrials.gov as NCT05146401. All participants gave written informed consent before enrolling in the study.
Selection of Gestational Diabetes Mellitus Cases and Controls
According to the diagnostic criteria advocated by International Association of the Diabetes and Pregnancy Study Groups (IADPSG)/WHO (27), the diagnosis of GDM can be made if one or more glucose values are above the cut points of 5.1, 10.0, and 8.5 mmol/L at fasting, 1 and 2 h during a 75-gram oral glucose tolerance test (OGTT). Controls were randomly selected and individually matched to cases by age (± 2 years), gestational age (± 2 weeks) and parity.
Measurement of Plasma Fatty Acids
The derivatization of fatty acids in plasma was achieved by the modified direct transesterification method proposed by Lepage and Roy (28). According to the method, total plasma lipids were hydrolyzed and derived to fatty acid methyl esters (FAMEs) for further instrumental analysis. Fifty microliter of plasma was mixed with 3.6 ml of methanol-acetyl chloride (8:1; v/v) in a disposable glass tube, then 10 μl internal standard composed of non-adecanoic acid methyl ester (800 μg/ml) and 150 μl hexane were added. The sample mixture was incubated at 100°C for 1 h and cooled on ice. Afterward, 2.5 ml of 12% K2CO3 solution was slowly added to stop the reaction and neutralize the mixture. The tubes were vortexed, followed by centrifugation at 3,000 rpm for 5 min. Then, the organic layer was collected for the subsequent analysis. FAMEs were separated and analyzed by an Agilent 7890B gas chromatography (GC) coupled with an Agilent 5977A Series mass spectrometry (MS). The inlet temperature was set as 250°C. The temperature of oven was initially held at 50°C for 1 min and increased to 100°C at a rate of 25°C/min, then to 130°C at a rate of 10°C/min, and finally raised to 250°C at a rate of 5°C/min and held for 5 min. Total run time was 33 min. As to scan modes, the total ion count was used for compound identification while the selected ion monitoring was used for quantification. The dilution series of the standard mixture of 37 FAMEs and the individual standard of C22:5n-3 FAME (both from Sigma Aldrich, Poole, United Kingdom) were used for the establishment of calibration curves. The molecular weights of FAMEs, selected ions, and retention times for fatty acids detection are shown in Supplementary Table 1. The validation of analytical method was performed by reference to the European Medicines Agency Guideline on bioanalytical method validation (29). Blank samples spiked with three known concentration levels (low, medium, and high) of quality control samples were used for accuracy and precision analyses. The intra- and inter-day accuracy (expressed as recovery values) of all analytes ranged from 86.04 to 114.78% and 86.11 to 114.11%, while the intra- and inter-day precision (expressed as the coefficient of variation) of all analytes ranged from 0.63 to 7.58% and 1.15 to 12.28%, respectively. Validation results were provided in detail in the Supplementary Tables 2–4.
Identification of Fatty Acid Patterns
In the current study, 18 trace fatty acids were excluded in the analyses because they were undetectable in most samples or the relative concentrations were <0.05% on average. After exclusion, 20 individual fatty acids, including 6 saturated fatty acids (SFAs), 4 monounsaturated fatty acids (MUFAs), and 10 polyunsaturated fatty acids (PUFAs), were used to derive FAPs through PCA. Principal components were inferred as FAPs, and we determined the number of patterns to analyze based on scree plot and eigenvalues (>1). The factor loadings reflect the contributions of individual fatty acids to principal components, which means the closer the absolute value is to 1, the greater the influence of single fatty acid on the pattern. A score of each pattern was calculated by summing fatty acid concentrations weighted by the scoring coefficients of each fatty acid.
Covariates
Sociodemographic factors, lifestyle, and health information were collected by trained interviewers with standardized questionnaires. Height was measured using calibrated instrument. Prepregnancy body mass index (BMI) was calculated as self-reported prepregnancy weight divided by height squared (kg/m2). Fasting concentrations of plasma glucose (FPG), total cholesterol, triglycerides, high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C) were measured using the commercial assay kits (Biosino Bio-Technology and Science, Inc.). Fasting plasma insulin (FPI) was measured by enzyme-linked immunoassay (ELISA) kits (Mercodia Company), and insulin resistance was assessed by homeostasis model assessment (HOMA) index.
Statistical Analyses
The concentration of total fatty acids (μmol/L) was calculated as a sum of individual fatty acids (μmol/L). Individual fatty acids were expressed as mole percentage of total fatty acids. Differences of continuous variables between case and control groups were tested by the Student’s t-test or Mann-Whitney U-test, when appropriate. For categorical variables, Chi-square (χ2) test were adopted. Spearman partial correlation coefficients were used to estimate the relations between individual fatty acids and plasma glucose in OGTT, lipids and insulin resistance. To evaluate the associations between fatty acids and odds of GDM, we used conditional logistic regression. Odds ratios (ORs) with 95% confidence intervals (CIs) of GDM were showed in quartiles which were based on the distribution of fatty acids among controls. Tests for linear trend were conducted by treating median value for each quartile of individual fatty acids as continuous variables. False discovery rate (FDR) correction was applied in the multiple comparisons of fatty acids to reduce the false positive rate (30). Potential confounding factors including age, gestational age at blood collection, parity, prepregnancy BMI, family history of diabetes, smoking and alcohol use were adjusted in multivariable models.
To explore the association between combinations of plasma fatty acids and GDM, we did further analyses. We treated the scores of each pattern as continuous variables and evaluated the association of FAP score with GDM odds by ORs and 95% CIs through conditional logistic regression, with potential confounders adjusted. In addition, a restricted cubic spline model with four knots was used to assess the possible non-linear associations. The relationships of FAP score with indexes of glucose and lipid metabolism were also evaluated by Spearman partial correlation coefficients. Data were analyzed with SPSS 20.0 software package (SPSS, Inc.) or Stata version 13.0 (StataCorp). All P-values were two-sided and the threshold for statistical significance was 0.05.
Results
Characteristics of the Participants
Among the 683 confirmed participants who met the inclusion criteria, 15 women were excluded for previous diagnosis of diabetes, 17 women were excluded because of supplements use that might influence glucose or lipids metabolism, 3 women were excluded for incomplete basic information, 105 participants were excluded due to insufficient plasma samples (mainly consumed by previous studies). Hence, 253 women with GDM and 290 healthy pregnant women were eligible for further case-control matching by SPSS software. Finally, 217 GDM cases and 217 matched controls were selected in this study.
The demographic, anthropometric, reproductive, and metabolic characteristics of the participants were exhibited in Table 1 by case-control status. The two groups were comparable for age, gestational age, and parity. Compared to controls, GDM cases had significantly higher levels of prepregnancy BMI, FPG, 1- and 2-h post-glucose load, FPI, triglycerides and insulin resistance values, and are prone to have a family history of diabetes. The concentrations of total fatty acids were significantly higher in cases than in controls [median: 2624.96 (interquartile range 2073.66–3472.03) vs. 2458.38 (interquartile range 1857.25–2985.84) μmol/L, P = 0.001].
Table 1. Characteristics among women with GDM and their matched controls.
Association of Individual Fatty Acids and Gestational Diabetes Mellitus
Plasma fatty acid composition was shown in Table 2. According to the methods, only fatty acids with relative concentrations (mol%) higher than 0.05% are displayed in the results. In both case and control groups, the most abundant fatty acids were 16:0 and 18:2n-6, whereas 24:0 and 15:0 were present in low relative concentrations. When compared GDM with control subjects, significant difference was found with several fatty acids, including five SFAs (14:0, 15:0, 17:0, 18:0, 24:0), three MUFAs (16:1n-7, 20:1n-9, 24:1n-9), and five PUFAs (18:3n-6, 20:2n-6, 18:3n-3, 20:3n-3, 22:5n-3).
Table 2. Composition (mol% of total fatty acids) of plasma fatty acids among GDM cases and non-GDM controls.
Spearman partial correlations between individual fatty acids and the indexes of glucose and lipid metabolism were shown in Table 3. HOMA-IR, the index that reflects insulin resistance levels, was correlated positively with 14:0, 16:0, 16:1n-7, 20:3n-6 (β = 0.14, 0.30, 0.14, and 0.15, respectively) and negatively with 15:0, 17:0, 24:0, 20:1n-9, 18:2n-6, 18:3n-6, 20:2n-6, 18:3n-3 (β =−0.13, −0.12, −0.21, −0.24, −0.21, −0.17, −0.17, and −0.12, respectively).
Table 3. Partial correlation coefficients between plasma individual fatty acids and variables of interesta.
In the conditional logistic regression analyses, the adjusted ORs (95% CIs) of GDM across increasing quartiles of plasma total fatty acid levels were 1.00 (referent), 1.43 (0.73–2.81), 1.11 (0.56–2.22), and 2.35 (1.24–4.47), respectively. Results for associations between fatty acid subgroups, individual fatty acids and GDM odds were displayed in Table 4 (more details could be found in Supplementary Tables 5,6). After adjustment for age, gestational age at blood collection, parity, prepregnancy BMI, family history of diabetes, smoking, and alcohol use, total SFAs, MUFAs, or PUFAs were not significantly associated with GDM odds. However, odd-chain SFAs were inversely associated with GDM odds, whether individually or combined. Among the 20 individual fatty acids included in the analyses, 13 fatty acids were considered to be negatively associated with odds of GDM. It is worth noting that the abundance of the 13 fatty acids were all relatively low, with a proportion range of 0.08% (24:0) to 6.72% (18:0). In addition, the OR (95% CI) for GDM compared the highest with lowest quartiles of 24:1n-9 was 2.05 (1.12–3.76).
Table 4. Association between fatty acid subgroups, individual fatty acids and GDMa.
Association of Novel Fatty Acid Patterns and Gestational Diabetes Mellitus
The correlation matrix of individual fatty acids was shown in Figure 1. Four major FAPs were identified and explained 31.41, 18.47, 13.56, and 9.68% of the overall variability, respectively. The scree plot starts to flatten from the fifth principal component (Figure 2), which was consistent with the above results.
Figure 1. Heat map for spearman correlation matrix of 20 fatty acids measured among 217 GDM cases and matched controls.
Figure 2. The scree plot derived from principal component analysis. The variance that explained by each component (%) was exhibited.
Factor loadings of the 20 fatty acids for the four patterns were shown in Figure 3. The first pattern (FAP1) was characterized by higher relative concentrations of odd-chain fatty acids, 14:0, 18:0, 18:3n-3, 20:2n-6, 20:3n-6, and lower relative concentrations of 24:1n-9 and 18:2n-6. The second pattern (FAP2) was characterized by higher abundance of n-3 PUFAs, 24:0, 18:3n-6 and lower proportions of 16:0 and 20:4n-6. The third pattern (FAP3) was featured by higher relative concentrations of 22:6n-3, 22:4n-6, 22:5n-3 and lower relative concentrations of 18:1n-9. The fourth pattern (FAP4) had moderated high factor loadings for 16:0 and 16:1n-7.
Figure 3. Factor loadings of the 20 fatty acids for the 1st to 4th principal components. The closer the factor loading (absolute value) is to 1, the greater the influence of single fatty acid on the pattern.
According to the scoring coefficients matrix (Supplementary Table 7), the score of each pattern can be calculated for each participant. Table 5 displays the partial correlation coefficients between FAP scores and metabolic factors. Higher FAP1 and FAP2 scores were modestly correlated with lower FPG, 1-h post glucose load, total cholesterol, triglycerides and LDL-C, and FAP2 score was further negatively correlated with 2-h post-glucose load, FPI, and HOMA-IR. Moreover, FAP3 score was positively correlated with triglycerides and HDL-C, while FAP4 was positively correlated with triglycerides, FPI, and HOMA-IR.
Table 5. Partial correlation coefficients between fatty acid pattern scores and variables of interesta.
After multivariable adjustment, the FAP1 and FAP2 were inversely associated with GDM odds (Table 6). Compared with women in the lowest quartile of FAP1 score, women in the highest quartile experienced a 58% lower odds of GDM (OR, 0.42; 95% CI, 0.23–0.76). The OR (95% CI) of GDM comparing extreme quartiles of FAP2 score was 0.48 (0.26–0.90). Potential non-linear associations of FAP1 and FAP2 with odds of GDM were also found in the restricted cubic spline model (Figure 4). Nevertheless, little evidence of an association was found between quartiles of FAP3 or FAP4 scores and odds of GDM.
Table 6. ORs (95% CIs) for GDM according to quartiles of plasma fatty acid pattern scores.
Figure 4. Potential non-linear associations of FAP1 (A) and FAP2 (B) with odds of GDM assessed by restricted cubic spline model. Solid line indicated adjusted ORs for GDM, and dashed line indicated 95% CIs. Adjustments were made for age, gestational age at blood collection, parity, prepregnancy BMI, family history of diabetes, smoking and alcohol use. FAP, fatty acid pattern; GDM, gestational diabetes mellitus; ORs, odds ratios; CIs, confidence intervals.
Discussion
In this matched case-control study, we initially found that 13 individual fatty acids including 5 SFAs (14:0, 15:0, 17:0, 18:0, 24:0), 2 MUFAs (16:1n-7, 20:1n-9), and 6 PUFAs (18:3n-6, 20:2n-6, 18:3n-3, 20:3n-3, 20:5n-3, 22:5n-3) were negatively associated with odds of GDM, while 24:1n-9 was positively associated with GDM odds. Afterward, we identified two certain patterns of circulating fatty acids which were associated with lower odds of GDM. The first pattern was characterized by higher relative concentrations of odd-chain fatty acids, 14:0, 18:0, 18:3n-3, 20:2n-6, 20:3n-6 and lower relative concentrations of 24:1n-9 and 18:2n-6. The second pattern was characterized by higher abundance of n-3 PUFAs, 24:0, 18:3n-6, and lower proportions of 16:0 and 20:4n-6. To our knowledge, this is the first study focused on circulating FAPs related to odds of GDM. The above results suggested that the combinations of circulating fatty acids could be a more significant marker of GDM development than individual fatty acids or their subgroups.
For fatty acids with relatively high concentrations, small changes in the levels of other fatty acids may have little impact on their relative amounts or the relationships with diseases. However, for fatty acids with relatively low concentrations, the influence of others might be more considerable. In our study, individual fatty acids which were significantly associated with GDM odds were mostly of the second type. Thus, investigation of the FAPs with PCA may confer benefits over investigating individual fatty acids or fatty acid groups, and provide new insight into the joint effect of circulating fatty acids on GDM development.
The FAP1 that we identified contained positive loadings for odd-chain fatty acids, 14:0, 18:0, 18:3n-3, 20:2n-6, 20:3n-6 and negative loadings for 24:1n-9 and 18:2n-6. Although the main fatty acids involved in this pattern were distributed in various subgroups, internal correlation existed. Linoleic acid (18:2n-6), the most abundant n-6 PUFA in human circulating, is the precursor for the endogenous synthesis of other long-chain n-6 PUFAs, including 20:2n-6, 20:3n-6, and 20:4n-6 (8). A recent study conducted by Kim et al. found that low 18:2n-6 diet and estrogen significantly increased the expression of peroxisome proliferator-activated receptor α (PPAR-α), fatty acid desaturase 2 (Δ6 desaturase), and elongases of vary long chain fatty acids in rats liver (31). The fatty acid desaturases and elongases play a major role in the levels of circulating long chain PUFAs (32). Linoleic acid can compete with α-linolenic acid (18:3n-3) for rate-limiting enzyme like Δ6 desaturase and affect circulating level of 18:3n-3 and its conversion into very long chain n-3 PUFAs eventually (33). As a ligand-dependent transcription factor, PPAR-α can be activated by fatty acids and their derivatives (34), and have been found to be able to downregulate the expression of Δ9-16 and Δ9-18 desaturases which are the rate-limiting enzymes in synthesis of 16:1n-7 and 18:1n-9. The above evidence may partly explain the inverse correlation of 18:2n-6 with 18:0, 18:3n-3, 20:2n-6 and 20:3n-6 in our study. Several potential mechanisms might contribute to the protective role of FAP1 in GDM development. Firstly, 18:2n-6 is the direct precursor of hydroxyl conjugated linoleic acid. Research evidence has suggested that lowering dietary 18:2n-6 from 6.7 to 2.4% of calories for 12 weeks could markedly reduce the abundance of human plasma hydroxyl conjugated linoleic acid which is the important component of oxidized low density lipoprotein and has the effect of promoting inflammatory (35). Inflammation and dyslipidemia are considered to be the major mechanisms of insulin resistance. Secondly, 18:2n-6 competes with 18:3n-3 for conversion process. Study conducted by Ghafoorunissa et al. found that the substitution of one-third of dietary 18:2n-6 with 18:3n-3 resulted in lowered blood lipid levels and increased insulin sensitivity in sucrose fed rats, possibly due to the resulting high long chain n-3 PUFA levels in target tissues of insulin action (36). Thirdly, odd-chain fatty acids, a subclass of SFA with very low abundance, are proved to be biomarkers of dairy fat and dietary fiber intake. A large meta-analysis that pooled findings from 16 prospective cohort studies indicated that higher levels of 15:0 and 17:0 were associated with a lower risk of type 2 diabetes (37). The inverse association between odd-chain fatty acids and diabetes was considered to be partly related to dietary fiber and dairy fiber intake (38).
The FAP2 was characterized by higher abundance of n-3 PUFAs, 24:0, 18:3n-6 and lower proportions of 16:0 and 20:4n-6. As the most abundant fatty acid in circulation, 16:0 can be directly obtained from foods or synthesized from other fatty acids, amino acids or carbohydrates via de novo lipogenesis. The level of 16:0 in tissue was relatively stable, and changes of dietary intake seemed have little influence on its concentration because the exogenous source was counterbalanced by the endogenous biosynthesis (39). However, pathological state, chronic malnutrition, or unhealthy lifestyle factors like physical inactivity may have great impact on de novo lipogenesis and further affect 16:0 level (40, 41). Mounting evidence suggested that disruption in 16:0 homeostasis was associated with the development of diseases including cardiovascular diseases, metabolic diseases, neurodegenerative diseases and cancers (39). Higher tissue 16:0 could induce inflammatory responses and metabolic dysregulation which further result in dyslipidemia, hyperglycemia, and insulin resistance (42). These findings have also been validated in animal and cell researches (43, 44). In addition, an optional intake of 16:0 in a certain ratio with PUFA (including n-3 and n-6 families) was considered to have significant effect on human health, which was in line with our findings on FAP2 (40). The in vitro experiments showed that 22:6n-3 could reduce 16:0 induced endoplasmic reticulum stress in pancreatic β cells (45), and reverse atherosclerotic changes in human endothelial cells induced by 16:0 (46). Moreover, n-3 and n-6 PUFAs are precursors of eicosanoids which have effects in mediating inflammation and regulating immune response. Interestingly, the eicosanoids derived from n-6 PUFAs are prone to be pro-inflammatory, whereas those derived from n-3 PUFAs are anti-inflammatory (47). This suggested that n-3 PUFAs may, to some extent, improve insulin sensitivity through modulating inflammation.
Up to date, no studies have explored the potential fatty acid patters that associated with GDM odds through statistical method of dimension reduction like PCA. The results of the few studies that assessed the association of individual fatty acids with incident GDM showed some inconsistency with the present study’s findings (9). Reasons for such variability include differences in gestational age at blood sample collection, biological specimen sampled, and the method of fatty acid measurement or expression. Pooled results from two nested case-control studies conducted in China indicated that higher levels of 14:0 and 16:0 were associated with higher odds of GDM, whereas higher 18:2n-6 was associated with lower odds of GDM (14). However, the fatty acids were measured for different plasma compartments in the two studies (one study used total plasma samples and another used plasma phospholipid contents), and significant associations found in one study may not be replicated in another one. Zhang et al. found that 18:3n-3 and 22:6n-3 in early pregnancy were associated with a higher risk of GDM subtype (17). One case-control study nested in a cohort of US women suggested a beneficial role of odd-chain fatty acids, 22:5n-3 and 22:4n-6 in prevention of GDM, while 16:0, 18:3n-6, and 20:3n-6 at 10–14 weeks of pregnancy were associated with increased risk of GDM (10, 11). In normal pregnancy, plasma phospholipids in maternal circulation would increase by nearly 50% as compared to non-pregnant circulation (9), and the composition of plasma fatty acids is expected to fluctuate throughout gestation, which may lead to the discrepancy in the association between individual fatty acids and GDM risk in different trimesters on pregnancy (15, 17). Ortega-Senovilla et al. found that in the serum of women with GDM, the concentrations of most fatty acids were lower than in control women, except for 20:4n-6 and 22:6n-3, which remained the same. However, when values were expressed as a percentage of total fatty acids, different results emerged, serum from women with GDM showed significantly higher proportions of 18:2n-6, 20:4n-6, 22:6n-3, and lower proportions of 16:0, 16:1n-7, 18:1n-9, 20:5n-3 (48). The discrepancy between the above results mainly due to the difference in fatty acid expression. Only 10 individual fatty acids were measured in this study, which may lead to the possibility of lacking of comparability with other studies when fatty acids were expressed as relative concentrations. It should also be noted that the placenta plays an essential role in determining how fatty acids are transferred from maternal to embryonic circulation (49). Multiple studies have evaluated placental preference of fatty acid transport and found that the placenta places a higher preference on docosahexaenoic acid transport (9). However, GDM appears to influence the transfer of PUFAs from mother to fetus. The percentages of docosahexaenoic acid, arachidonic acid, and n-6 and n-3 PUFAs were found to be lower in the cord blood of mothers with GDM than in controls (50). Hence, the impaired transfer of some certain fatty acids through the placenta to the cord blood and fetus may be a possible factor in changing the patterns of fatty acids in GDM women.
Previous studies have explored the FAPs that related to a series of diseases such as cardiovascular disease (24), type 2 diabetes (25), metabolic syndrome (51), and cancer (52). A large nested case-cohort study indicated that a combination of fatty acids that characterized by high concentrations of 18:2n-6, odd-chain fatty acids, and very long-chain fatty acids, and low concentrations of 18:3n-6 and 16:0 was associated with lower risk of type 2 diabetes (25). Another previous study using PCA on 11 fatty acids in serum found that low 18:2n-6 factor and n-3 PUFA factor predicted metabolic syndrome development over 20 years, independent of lifestyle factors (51). Significant differences of the results on identified fatty acid patters were existed between our study and studies above, which might be attributable to several reasons. First, although metabolic syndrome, type 2 diabetes, and GDM are all metabolic diseases with insulin resistance, diverse pathogenesis, and potential risk factors still exist. Besides, circulating fatty acids could be selectively transported from pregnant women to the fetus through placenta (53), which results in the difference of fatty acid composition in the plasma of pregnant women compared with non-pregnant adults (8, 54). Furthermore, differences in study design, such as the preference of biological samples, measurement of fatty acids, sample size, and the analytical method, were all the possible reasons for the discrepancy of the results. Thus, more prospective studies are needed to validate our findings.
The major strength of our study is the application of PCA on FAPs’ exploring, making it possible to take into account the interaction of individual fatty acids. Additionally, plasma fatty acids, the objective biomarker of fatty acid status, were measured by GC–MS, which was independent of diet records and recall bias. Further, all the GDM cases included in this study were new cases and were not managed with lifestyle counseling or treated with medicine. However, several limitations should also be considered. Firstly, the sample size in the current study was relatively small. Secondly, the non-prospective design of this study disenabled us to infer the causal relations between FAPs and GDM development. The trajectory of fatty acid levels across the duration of gestation may be more informative than single-point measurement. Thus, large prospective study with longitudinal data collection is warrant in future. Thirdly, whether the FAPs identified in this study are related to certain dietary patterns or food preferences remains unclear. This could also be an important direction for future research. Fourthly, despite the adjustment for potential confounding factors in study design and statistical process, we cannot exclude the possibility that residual confounders existed.
In conclusion, in this matched case-control study, we identified two novel FAPs that were inversely associated with GDM odds. The combination of circulating fatty acids could be a more significant marker of GDM development than individual fatty acids or their subgroups. Prospective studies in other populations are needed to validate our findings, and explore how to optimize FAPs during pregnancy to achieve better health outcomes.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Tongji Medical College. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
PL, AT, and LL designed the research. PL, SH, YZ, TS, YH, ZX, HL, and SZ contributed to the data collection and analysis. PL wrote the manuscript. CL, AT, and LL edited the manuscript. LL and AT are the guarantors of this work and have full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. All authors read and approved the final manuscript.
Funding
This work was supported by the Young Scientists Fund of the National Natural Science Foundation of China (82103822 and 81803226), the Natural Science Foundation of Hubei Science and Technology Commission (2020CFB360), and the Major International (Regional) Joint Research Project (NSFC 81820108027).
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.
The reviewer XL declared a shared affiliation with the authors PL, SH, YZ, TS, YH, ZX, HL, CL, SZ, and LL to the handling editor at the time of review.
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 gratefully acknowledge the cooperation of all the pregnant women who took part in this study. We are also sincerely grateful to the staffs in Tongji Hospital for their considerable assistance with many aspects.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2022.836115/full#supplementary-material
Abbreviations
BMI, body mass index; FAME, fatty acid methyl ester; FAP, fatty acid pattern; FPG, fasting plasma glucose; FPI, fasting plasma insulin; GC–MS, gas chromatography–mass spectrometry; GDM, gestational diabetes mellitus; HDL-C, high density lipoprotein cholesterol; HOMA-IR, homeostasis model assessment of insulin resistance; LDL-C, low density lipoprotein cholesterol; MUFA, monounsaturated fatty acid; OGTT, oral glucose tolerance test; PCA, principal components analysis; PUFA, polyunsaturated fatty acid; SFA, saturated fatty acid.
References
1. International Diabetes Federation.IDF Diabetes Atlas. 10th ed. (2022). Available online at: http://www.diabetesatlas.org/ (accessed January 5, 2022).
2. Farrar D, Simmonds M, Bryant M, Sheldon TA, Tuffnell D, Golder S, et al. Hyperglycaemia and risk of adverse perinatal outcomes: systematic review and meta-analysis. BMJ. (2016) 354:i4694. doi: 10.1136/bmj.i4694
3. Song C, Lyu Y, Li C, Liu P, Li J, Ma RC, et al. Long-term risk of diabetes in women at varying durations after gestational diabetes: a systematic review and meta-analysis with more than 2 million women. Obes Rev. (2018) 19:421–9. doi: 10.1111/obr.12645
4. Damm P, Houshmand-Oeregaard A, Kelstrup L, Lauenborg J, Mathiesen ER, Clausen TD. Gestational diabetes mellitus and long-term consequences for mother and offspring: a view from Denmark. Diabetologia. (2016) 59:1396–9. doi: 10.1007/s00125-016-3985-5
5. Lowe WL Jr. Hyperglycemia and adverse pregnancy outcome follow-up study (HAPO FUS): maternal gestational diabetes mellitus and childhood glucose metabolism. Diabetes Care. (2019) 42:372–80. doi: 10.2337/dc18-1646
6. de Souza RJ, Mente A, Maroleanu A, Cozma AI, Ha V, Kishibe T, et al. Intake of saturated and trans unsaturated fatty acids and risk of all cause mortality, cardiovascular disease, and type 2 diabetes: systematic review and meta-analysis of observational studies. BMJ. (2015) 351:h3978. doi: 10.1136/bmj.h3978
7. Brown TJ, Brainard J, Song F, Wang X, Abdelhamid A, Hooper L. Omega-3, omega-6, and total dietary polyunsaturated fat for prevention and treatment of type 2 diabetes mellitus: systematic review and meta-analysis of randomised controlled trials. BMJ. (2019) 366:l4697. doi: 10.1136/bmj.l4697
8. Hodson L, Skeaff CM, Fielding BA. Fatty acid composition of adipose tissue and blood in humans and its use as a biomarker of dietary intake. Prog Lipid Res. (2008) 47:348–80. doi: 10.1016/j.plipres.2008.03.003
9. Mauro AK, Rengarajan A, Albright C, Boeldt DS. Fatty acids in normal and pathological pregnancies. Mol Cell Endocrinol. (2022) 539:111466. doi: 10.1016/j.mce.2021.111466
10. Zhu Y, Li M. Plasma phospholipid n-3 and n-6 polyunsaturated fatty acids in relation to cardiometabolic markers and gestational diabetes: a longitudinal study within the prospective NICHD fetal growth Studies. PLoS Med. (2019) 16:e1002910. doi: 10.1371/journal.pmed.1002910
11. Zhu Y, Tsai MY, Sun Q, Hinkle SN, Rawal S, Mendola P, et al. A prospective and longitudinal study of plasma phospholipid saturated fatty acid profile in relation to cardiometabolic biomarkers and the risk of gestational diabetes. Am J Clin Nutr. (2018) 107:1017–26. doi: 10.1093/ajcn/nqy051
12. Li X, Huang Y, Xing Y, Hu C, Zhang W, Tang Y, et al. Association of urinary cadmium, circulating fatty acids, and risk of gestational diabetes mellitus: a nested case-control study in China. Environ Int. (2020) 137:105527. doi: 10.1016/j.envint.2020.105527
13. Hosseinkhani S, Dehghanbanadaki H, Aazami H, Pasalar P, Asadi M, Razi F. Association of circulating omega 3, 6 and 9 fatty acids with gestational diabetes mellitus: a systematic review. BMC Endocr Disord. (2021) 21:120. doi: 10.1186/s12902-021-00783-w
14. Pan XF, Huang Y, Li X, Wang Y, Ye Y, Chen H, et al. Circulating fatty acids and risk of gestational diabetes mellitus: prospective analyses in China. Eur J Endocrinol. (2021) 185:87–97. doi: 10.1530/eje-21-0118
15. Li L, Zhu Y, Wu J, Hinkle SN, Tobias DK, Ma RCW, et al. Changes of plasma phospholipid fatty acids profiles in pregnancy in relation to the diagnosis and treatment of gestational diabetes mellitus. Clin Chem. (2021) 67:1660–75. doi: 10.1093/clinchem/hvab169
16. Tryggvadottir EA, Gunnarsdottir I, Birgisdottir BE, Hrolfsdottir L, Landberg R, Hreidarsdottir IT, et al. Early pregnancy plasma fatty acid profiles of women later diagnosed with gestational diabetes. BMJ Open Diabetes Res Care. (2021) 9:e002326. doi: 10.1136/bmjdrc-2021-002326
17. Zhang T, Jiang WR, Xia YY, Mansell T, Saffery R, Cannon RD, et al. Complex patterns of circulating fatty acid levels in gestational diabetes mellitus subclasses across pregnancy. Clin Nutr. (2021) 40:4140–8. doi: 10.1016/j.clnu.2021.01.046
18. Voortman T, Tielemans MJ, Stroobant W, Schoufour JD, Kiefte-de Jong JC, Steenweg-de Graaff J, et al. Plasma fatty acid patterns during pregnancy and child’s growth, body composition, and cardiometabolic health: the generation R study. Clin Nutr. (2018) 37:984–92. doi: 10.1016/j.clnu.2017.04.006
19. Stark KD, Van Elswyk ME, Higgins MR, Weatherford CA, Salem N Jr. Global survey of the omega-3 fatty acids, docosahexaenoic acid and eicosapentaenoic acid in the blood stream of healthy adults. Prog Lipid Res. (2016) 63:132–52. doi: 10.1016/j.plipres.2016.05.001
20. Ringner M. What is principal component analysis? Nat Biotechnol. (2008) 26:303–4. doi: 10.1038/nbt0308-303
21. Lara KM, Levitan EB, Gutierrez OM, Shikany JM, Safford MM, Judd SE, et al. Dietary patterns and incident heart failure in U.S. Adults without known coronary disease. J Am Coll Cardiol. (2019) 73:2036–45. doi: 10.1016/j.jacc.2019.01.067
22. Dahm CC, Gorst-Rasmussen A, Jakobsen MU, Schmidt EB, Tjonneland A, Sorensen TI, et al. Adipose tissue fatty acid patterns and changes in anthropometry: a cohort study. PLoS One. (2011) 6:e22587. doi: 10.1371/journal.pone.0022587
23. Yang B, Ding F, Yan J, Ye XW, Xu XL, Wang FL, et al. Exploratory serum fatty acid patterns associated with blood pressure in community-dwelling middle-aged and elderly Chinese. Lipids Health Dis. (2016) 15:58. doi: 10.1186/s12944-016-0226-3
24. Imamura F, Lemaitre RN, King IB, Song X, Lichtenstein AH, Matthan NR, et al. Novel circulating fatty acid patterns and risk of cardiovascular disease: the cardiovascular health study. Am J Clin Nutr. (2012) 96:1252–61. doi: 10.3945/ajcn.112.039990
25. Imamura F, Sharp SJ, Koulman A, Schulze MB, Kroger J, Griffin JL, et al. A combination of plasma phospholipid fatty acids and its association with incidence of type 2 diabetes: the EPIC-interact case-cohort study. PLoS Med. (2017) 14:e1002409. doi: 10.1371/journal.pmed.1002409
26. Dahm CC, Gorst-Rasmussen A, Crowe FL, Roswall N, Tjonneland A, Drogan D, et al. Fatty acid patterns and risk of prostate cancer in a case-control study nested within the European prospective investigation into cancer and nutrition. Am J Clin Nutr. (2012) 96:1354–61. doi: 10.3945/ajcn.112.034157
27. American Diabetes Association. 2. Classification and diagnosis of diabetes: standards of medical care in diabetes-2021. Diabetes Care. (2021) 44(Suppl. 1):S13–27. doi: 10.2337/dc21-S002
28. Lepage G, Roy CC. Direct transesterification of all classes of lipids in a one-step reaction. J Lipid Res. (1986) 27:114–20.
29. European Medicines Agency.Guideline on Validation of Bioanalytical Methods. EMEA/CHMP/EWP/192217/2009. (2009). Available online at: http://www.ema.europa.eu/ (accessed January 5, 2022).
30. Glickman ME, Rao SR, Schultz MR. False discovery rate control is a recommended alternative to Bonferroni-type adjustments in health studies. J Clin Epidemiol. (2014) 67:850–7. doi: 10.1016/j.jclinepi.2014.03.012
31. Kim D, Choi JE, Park Y. Low-linoleic acid diet and oestrogen enhance the conversion of alpha-linolenic acid into DHA through modification of conversion enzymes and transcription factors. Br J Nutr. (2019) 121:137–45. doi: 10.1017/s0007114518003252
32. Lankinen M, Uusitupa M, Schwab U. Genes and dietary fatty acids in regulation of fatty acid composition of plasma and erythrocyte membranes. Nutrients. (2018) 10:1785. doi: 10.3390/nu10111785
33. Goyens PL, Spilker ME, Zock PL, Katan MB, Mensink RP. Conversion of alpha-linolenic acid in humans is influenced by the absolute amounts of alpha-linolenic acid and linoleic acid in the diet and not by their ratio. Am J Clin Nutr. (2006) 84:44–53. doi: 10.1093/ajcn/84.1.44
34. Wahli W, Michalik L. PPARs at the crossroads of lipid signaling and inflammation. Trends Endocrinol Metab. (2012) 23:351–63. doi: 10.1016/j.tem.2012.05.001
35. Ramsden CE, Ringel A, Feldstein AE, Taha AY, MacIntosh BA, Hibbeln JR, et al. Lowering dietary linoleic acid reduces bioactive oxidized linoleic acid metabolites in humans. Prostaglandins Leukot Essent Fatty Acids. (2012) 87:135–41. doi: 10.1016/j.plefa.2012.08.004
36. Ghafoorunissa, Ibrahim A, Natarajan S. Substituting dietary linoleic acid with alpha-linolenic acid improves insulin sensitivity in sucrose fed rats. Biochim Biophys Acta. (2005) 1733:67–75. doi: 10.1016/j.bbalip.2004.12.003
37. Imamura F, Fretts A, Marklund M, Ardisson Korat AV, Yang WS, Lankinen M, et al. Fatty acid biomarkers of dairy fat consumption and incidence of type 2 diabetes: a pooled analysis of prospective cohort studies. PLoS Med. (2018) 15:e1002670. doi: 10.1371/journal.pmed.1002670
38. Weitkunat K, Schumann S, Nickel D, Hornemann S, Petzke KJ, Schulze MB, et al. Odd-chain fatty acids as a biomarker for dietary fiber intake: a novel pathway for endogenous production from propionate. Am J Clin Nutr. (2017) 105:1544–51. doi: 10.3945/ajcn.117.152702
39. Carta G, Murru E, Banni S, Manca C. Palmitic acid: physiological role, metabolism and nutritional implications. Front Physiol. (2017) 8:902. doi: 10.3389/fphys.2017.00902
40. Sanders FWB, Acharjee A, Walker C, Marney L, Roberts LD, Imamura F, et al. Hepatic steatosis risk is partly driven by increased de novo lipogenesis following carbohydrate consumption. Genome Biol. (2018) 19:79. doi: 10.1186/s13059-018-1439-8
41. Brouwers B, Hesselink MK, Schrauwen P, Schrauwen-Hinderling VB. Effects of exercise training on intrahepatic lipid content in humans. Diabetologia. (2016) 59:2068–79. doi: 10.1007/s00125-016-4037-x
42. Donnelly KL, Smith CI, Schwarzenberg SJ, Jessurun J, Boldt MD, Parks EJ. Sources of fatty acids stored in liver and secreted via lipoproteins in patients with nonalcoholic fatty liver disease. J Clin Invest. (2005) 115:1343–51. doi: 10.1172/jci23621
43. Hirata T, Kawai T, Hirose H, Tanaka K, Kurosawa H, Fujii C, et al. Palmitic acid-rich diet suppresses glucose-stimulated insulin secretion (GSIS) and induces endoplasmic reticulum (ER) stress in pancreatic islets in mice. Endocr Res. (2016) 41:8–15. doi: 10.3109/07435800.2015.1038352
44. Calvo-Ochoa E, Sánchez-Alegría K, Gómez-Inclán C, Ferrera P, Arias C. Palmitic acid stimulates energy metabolism and inhibits insulin/PI3K/AKT signaling in differentiated human neuroblastoma cells: the role of mTOR activation and mitochondrial ROS production. Neurochem Int. (2017) 110:75–83. doi: 10.1016/j.neuint.2017.09.008
45. Suzuki E, Matsuda T, Kawamoto T, Takahashi H, Mieda Y, Matsuura Y, et al. Docosahexaenoic acid reduces palmitic acid-induced endoplasmic reticulum stress in pancreatic B cells. Kobe J Med Sci. (2018) 64:E43–55.
46. Karbasforush S, Nourazarian AA-O, Darabi M, Rahbarghazi RA-O, Khaki-Khatibi F, Biray Avci Ç, et al. Docosahexaenoic acid reversed atherosclerotic changes in human endothelial cells induced by palmitic acid in vitro. Cell Biochem Funct. (2018) 36:203–11. doi: 10.1002/cbf.3332
47. Lands B. Historical perspectives on the impact of n-3 and n-6 nutrients on health. Prog Lipid Res. (2014) 55:17–29. doi: 10.1016/j.plipres.2014.04.002
48. Ortega-Senovilla H, Schaefer-Graf U, Herrera E. Pregnant women with gestational diabetes and with well controlled glucose levels have decreased concentrations of individual fatty acids in maternal and cord serum. Diabetologia. (2020) 63:864–74. doi: 10.1007/s00125-019-05054-x
49. Devarshi PP, Grant RW, Ikonte CJ, Hazels Mitmesser S. Maternal omega-3 nutrition, placental transfer and fetal brain development in gestational diabetes and preeclampsia. Nutrients. (2019) 11:1107. doi: 10.3390/nu11051107
50. Hai-Tao Y, Zhi-Heng G, Yi-Ru C, Yue-Ting L, Hai-Ying Z, Ya-Juan L, et al. Gestational diabetes mellitus decreased umbilical cord blood polyunsaturated fatty acids: a meta-analysis of observational studies. Prostaglandins Leukot Essent Fatty Acids. (2021) 171:102318. doi: 10.1016/j.plefa.2021.102318
51. Warensjo E, Sundstrom J, Lind L, Vessby B. Factor analysis of fatty acids in serum lipids as a measure of dietary fat quality in relation to the metabolic syndrome in men. Am J Clin Nutr. (2006) 84:442–8. doi: 10.1093/ajcn/84.1.442
52. Yang M, Ayuningtyas A, Kenfield SA, Sesso HD, Campos H, Ma J, et al. Blood fatty acid patterns are associated with prostate cancer risk in a prospective nested case-control study. Cancer Causes Control. (2016) 27:1153–61. doi: 10.1007/s10552-016-0794-6
53. De Vriese SR, Dhont M, Christophe AB. FA composition of cholesteryl esters and phospholipids in maternal plasma during pregnancy and at delivery and in cord plasma at birth. Lipids. (2003) 38:1–7. doi: 10.1007/s11745-003-1023-y
Keywords: fatty acid, gestational diabetes mellitus, pattern, case-control study, pregnancy
Citation: Li P, Hu S, Zhu Y, Sun T, Huang Y, Xu Z, Liu H, Luo C, Zhou S, Tan A and Liu L (2022) Associations of Plasma Fatty Acid Patterns During Pregnancy With Gestational Diabetes Mellitus. Front. Nutr. 9:836115. doi: 10.3389/fnut.2022.836115
Received: 15 December 2021; Accepted: 15 February 2022;
Published: 06 May 2022.
Edited by:
Rosa Casas Rodriguez, Institut de Recerca Biomèdica August Pi i Sunyer (IDIBAPS), SpainReviewed by:
Xinping Li, Huazhong University of Science and Technology, ChinaEmilio Herrera, CEU San Pablo University, Spain
Copyright © 2022 Li, Hu, Zhu, Sun, Huang, Xu, Liu, Luo, Zhou, Tan and Liu. 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: Aijun Tan, zhcdc@sina.com; Liegang Liu, lgliu@mails.tjmu.edu.cn