- 1Department of Fetal Medicine & Prenatal Diagnosis Center, Shanghai Key Laboratory of Maternal Fetal Medicine, Shanghai Institute of Maternal-Fetal Medicine and Gynecologic Oncology, Shanghai First Maternity and Infant Hospital, School of Medicine, Tongji University, Shanghai, China
- 2Department of Gynecology and Obstetrics, The First Affiliated Hospital of Dalian Medical University, Dalian, Liaoning, China
Background: Fetal growth restriction (FGR) occurs in 10% of pregnancies worldwide. Placenta dysfunction, as one of the most common causes of FGR, is associated with various poor perinatal outcomes. The main objectives of this study were to screen potential diagnostic biomarkers for FGR and to evaluate the function of immune cell infiltration in the process of FGR.
Methods: Firstly, differential expression genes (DEGs) were identified in two Gene Expression Omnibus (GEO) datasets, and gene set enrichment analysis was performed. Diagnosis-related key genes were identified by using three machine learning algorithms (least absolute shrinkage and selection operator, random forest, and support vector machine model), and the nomogram was then developed. The receiver operating characteristic curve, calibration curve, and decision curve analysis curve were used to verify the validity of the diagnostic model. Using cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT), the characteristics of immune cell infiltration in placental tissue of FGR were evaluated and the candidate key immune cells of FGR were screened. In addition, this study also validated the diagnostic efficacy of TREM1 in the real world and explored associations between TREM1 and various clinical features.
Results: By overlapping the genes selected by three machine learning algorithms, four key genes were identified from 290 DEGs, and the diagnostic model based on the key genes showed good predictive performance (AUC = 0.971). The analysis of immune cell infiltration indicated that a variety of immune cells may be involved in the development of FGR, and nine candidate key immune cells of FGR were screened. Results from real-world data further validated TREM1 as an effective diagnostic biomarker (AUC = 0.894) and TREM1 expression was associated with increased uterine artery PI (UtA-PI) (p-value = 0.029).
Conclusion: Four candidate hub genes (SCD, SPINK1, TREM1, and HIST1H2BB) were identified, and the nomogram was constructed for FGR diagnosis. TREM1 was not only associated with a variety of key immune cells but also correlated with increased UtA-PI. The results of this study could provide some new clues for future research on the prediction and treatment of FGR.
1 Introduction
Fetal growth restriction (FGR) affects 10% of all pregnancies worldwide and is a major cause of poor perinatal outcomes (1, 2). It is a condition in which a fetus does not attain its genetically conferred growth potential because of underlying pathologies, including maternal, fetal, infection, and placental abnormalities. FGR is commonly defined as the estimated fetal weight (EFW) or abdominal circumference (AC) less than the tenth percentile for the gestational age (3, 4).
The placenta is an important organ during pregnancy as the interface between fetal development and maternal circulation (5). Placenta dysfunction is a common cause of FGR, and accounts for 25%~30% of all FGR cases, which can increase the risk of iatrogenic preterm birth, very low birth weight, and poor long-term neurological prognosis (6, 7). Many major clinical problems in human pregnancy, such as FGR and preeclampsia (PE), although classically presenting in the third trimester, have their origins in the first trimester when immune regulation at the mother-fetus interface is abnormal, namely placenta dysfunction (8). Several studies have shown that immune regulation at the mother-fetus interface in embryo implantation, decidualization, and placentation (9, 10) plays an important role in the process, and is an important factor affecting the outcome of pregnancy. Furthermore, with the development of high-throughput sequencing technology, it is possible to explore the relationship between the expression of pathogenic genes and FGR caused by placental dysfunction at the transcriptome and epigenetics levels. In vitro and in vivo studies, differential expression of multiple genes have been detected through cell models and placental tissue, which may contribute to abnormal placental function through various ways, including immune regulation, and are strongly associated with pregnancy-related diseases (11–13). Additionally, given its high incidence and associated mortality rate, timely diagnosis and prediction of FGR are linked to improved outcomes (4). Various invasive diagnostic methods, such as chorionic villus sampling, have been proven effective in diagnosing and analyzing aberrant gene regulatory networks in fetuses suspected of having FGR (2). These results add new perspectives to the study of the mechanisms of FGR on one hand and provide possible clues for molecular or drug therapies targeting candidate genes and diagnostic genes to prevent and predict placenta-related diseases on the other.
Despite the exciting results of genes research in the human placenta, there is still a lack of effective diagnostic biomarkers for FGR, and our current knowledge of how numerous genes contribute to human placental development and function in pregnancy processes is still very limited. The critical role of immune factors in embryo implantation and placenta formation has also received extensive attention (14).
In this study, we used three machine learning algorithms to explore the potential biomarkers and underlying pathways involved in the development of FGR and conducted immune-related analysis of FGR placental tissue. In addition, real-world data were collected to validate the diagnostic efficacy of the biomarkers of interest identified in this study and to explore its relationship with multiple clinical parameters.
2 Methods
2.1 Data processing and identification of differentially expressed genes
GSE147776 and GSE203507 are data sets from GEO database for gene expression microarray analysis of placental tissues. GSE147776, based on the GPL20844 platform, contains 13 samples of FGR and 8 samples of normal pregnancy. GSE203507, based on the GPL16791 platform, contains 21 samples of FGR and 10 samples of normal pregnancy (Supplementary Table 1). By applying R packages “limma” and “sva”, we standardized the data in datasets, performed batch effect correction, and screened differentially expressed genes (DEGs). The cut-off criteria were adjusted p < 0.05 and | log fold change (FC)| > 1. The heatmap and volcano diagram were obtained using the “ggplots” package.
2.2 Enrichment analysis
Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed by using the “clusterProfiler” and “pathview” R packages to predict the potential function of DEGs between two groups. The potential underlying molecular mechanisms of DEGs were further investigated by applying gene set enrichment analysis (GSEA). In addition, we performed disease ontology (DO) enrichment analysis of DEGs using the “Dose” R package (15).
2.3 Hub gene identification and construction and evaluation of diagnostic model
In this experiment, hub genes were defined as overlapping genes from three machine learning algorithms, including: least absolute shrinkage and selection operator (LASSO) regression curve with “glmnet” package (16), random forest (RF) with the “randomForest” package (17) and support vector machine model (SVM-RFE) with the “e1071” package (18). Based on the selected hub genes, we constructed a nomogram model through the “rms” package to predict the incidence of FGR. The calibration curve and receiver operating characteristic (ROC) curve were drawn to evaluate the suitability of our nomogram for clinical use (19, 20).
2.4 Gene set enrichment analysis and immunological correlation analysis
The “clusterProfiler” package was also used to conduct GSEA enrichment analysis to explore potential molecular mechanisms for TRMDGs. The “CIBERSORT” package was used to analyze differences in 22 immune cell infiltration levels between FGR and control groups. Box plots can show the differences in immune cells between the different groups. LASSO regression analysis was performed to identify the FGR group of candidate key immune cells from 22 types of immune cells. In addition, the correlation relationship between TREM1 and the immune cell infiltration level was analyzed and a lollipop plot was drawn to summarize the correlation between each immune cell and TREM1.
2.5 Validation based on the real-world data
To validate the expression changes and the prognostic value of TREM1, placenta tissues were collected from 62 singleton pregnancies with FGR and 24 singleton controls undergoing routine prenatal monitoring at the Department of Fetal Medicine & Prenatal Diagnosis Center of Shanghai First Maternity and Infant Hospital in China from January 2021 to December 2023. The clinical characteristics of patients are presented in Table 1 and Supplementary Table 2. FGR was defined as birthweight below the 10th percentile for gestational age. Cases complicated with other fetal major structural and chromosomal anomalies were excluded. All placental tissues were collected during cesarean sections or spontaneous vaginal delivery. Tissue samples of approximately 1 cm3 from various regions near the cord attachment on the maternal side were dissected. These tissues were rinsed in saline solution to eliminate maternal contamination. A portion of the tissues was fixed in a 4% paraformaldehyde solution, embedded in paraffin, and subjected to immunohistochemistry staining (21). The remaining samples were removed the maternal side and fetal side of the placenta, snap-frozen in liquid nitrogen and stored at −80°C for subsequent RT-PCR analysis (22) and Western blot (23). Then, we also underwent immunohistochemical (IHC) staining. Detailed methods are found in the Supplementary Materials.
Partial clinical information from patients in the FGR group and the control group also were collected and the following variables for this study were extracted: age, gender of offspring, fetal weight, features of umbilical artery (UA) and uterine artery (UtA) Doppler examination in pregnancy, and the history of PE. The pulsatility index (PI) was estimated by the ultrasonography software.
TREM1 differentially expressed analysis between groups was performed using the “limma” R package. Spearman’s correlation analysis was used to describe the correlation between TREM1 expression and clinical traits. Furthermore, a ROC curve was used to assess the predictive value of TREM1 for FGR in the real world. This study was approved by the Ethics Committee of Shanghai First Maternity and Infant Hospital (Ethical number KS2133), and written informed consent was obtained from all participants before the collection of clinical samples and data.
2.6 Statistics analysis
All statistical analyses were conducted using R packages [R software (version 4.2.0)]. A p-value < 0.05 was considered statistically significant (p-value < 0.001 = ∗∗∗, p-value < 0.01 = ∗∗, and p-value < 0.05 = ∗). The process and study design are presented in a flowchart (Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1381795/fimmu-15-1381795-HTML/image_m/fimmu-15-1381795-g001.jpg)
Figure 1. The flow chart for the whole design. DEGs, differentially expressed genes. CIBERSORT, cell-type identification by estimating relative subsets of RNA transcripts. GO, Gene Ontology; DO, disease ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; GSEA, gene set enrichment analysis; SVM-RFE, support vector machine-recursive feature elimination; LASSO, least absolute shrinkage and selection operator; IHC, immunohistochemical.
3 Results
3.1 Identification of differentially expressed genes
After standardization and batch effect removal of two GEO datasets (Figures 2A, B), a total of 290 DEGs were screened (Figure 2C). Top 20 up and down-regulated genes were shown in Figure 2D. Among them, 133 genes were up-regulated, and 157 genes were down-regulated in the FGR group (Supplementary Table 3). The protein-protein interaction network of DEGs was analysised through the public database STRING and we visualized the result using Cytoscape, a network visualization software (Supplementary Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1381795/fimmu-15-1381795-HTML/image_m/fimmu-15-1381795-g002.jpg)
Figure 2. Combining different datasets and identification of DEGs. (A, B) Boxplots of mRNA expression distribution before and after removing batch effects. (C, D) DEG heatmap and volcano plot between FGR and control group. DEGs, differentially expressed genes; FGR, fetal growth restriction.
3.2 Enrichment analysis
Enrichment analysis of 290 DEGs was performed. Specifically, in the DO enrichment analysis, we found that DEGs may be related to the occurrence and development of female reproductive system diseases, asthma, pre−eclampsia and other diseases (Supplementary Figure 2A). The specific p-values of the DO enrichment analysis were shown in Supplementary Table 4. On the other hand, biological process (BP) analyses showed that leukocyte migration, regulation of cell-cell adhesion, and lymphocyte mediated immunity are the biological activities in which the DEGs are primarily involved. The result of cellular component (CC) analysis showed that the DEGs are involved in the composition of lumenal side of endoplasmic reticulum membrane, basal plasma membrane, and basolateral plasma membrane. Immune receptor activity, glycosaminoglycan binding, and cytokine receptor binding were mainly enriched according to the molecular function (MF) analysis (Supplementary Figure 2B). KEGG pathways analysis also indicated that DEGs are enriched in several cytokine interactions and immune-related pathways (Supplementary Figure 2C). GSEA analysis was performed, and the enriched pathways in the FGR group are presented in Supplementary Figure 1D.
3.3 Screening for key genes based on machine learning models
The 290 pivotal gene expression profiles were used to construct prediction functions using three machine learning models: Firstly, a total of 290 genes were screened from the DEGs via 10-fold cross-validation by SVM-RFE algorithm as diagnostic markers (Figure 3A). LASSO regression selected 18 predicted genes from among the statistically significant univariate variables out of candidate variables (Figure 3B). And 50 genes were screened from DEGs using RF and top 30 genes relative relevance was ranked from high to low (Figure 3C). Finally, the three algorithms identified four key genes (SCD, SPINK1, TREM1 and HIST1H2BB) with overlap (Figure 3D). In addition, we also found a correlation between the expression levels of four hub genes in the FGR group, and the specific correlation coefficients and p-values are displayed in Figure 3E. We observed significant intercorrelations among the expression levels of these four hub genes in the FGR placenta.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1381795/fimmu-15-1381795-HTML/image_m/fimmu-15-1381795-g003.jpg)
Figure 3. Hub gene identification and correlation analysis. (A) based on SVM-RFE to screen biomarkers. (B) LASSO logistic regression algorithm to screen diagnostic markers. (C) top-10 genes according to their discriminant ability in the RF algorithm. (D) The Venn diagram showed the intersection of diagnostic markers obtained by the results of three algorithms. (E) Correlation between 4 diagnostic genes. *P< 0.05, **P< 0.01, ***P< 0.001. SVM-RFE, support vector machine-recursive feature elimination; LASSO, least absolute shrinkage and selection operator; RF, random forest.
3.4 Construction of diagnostic model and developing a nomogram
The expression levels of four key genes in the FGR group and the control group were extracted, and a multivariate Logistic regression model was constructed. The results of Logistic regression model suggested that four key genes have predictive efficacy for the occurrence of FGR (SCD, p-value < 0.001; SPINK1, p-value = 0.001; TREM1, p-value = 0.001; HIST1H2BB, p-value = 0.008) (Supplementary Table 5). To graphically evaluate each individual (SCD, SPINK1, TREM1 and HIST1H2BB), a nomogram was developed based on 4 key genes to predict the likelihood of FGR in each fetus (Figure 4A). According to the results of the decision curve analysis (DCA), the nomogram model offered a better clinical benefit. Furthermore, the calibration plot indicated that the nomogram operated in line with the ideal model (Figures 4B, C). The result of ROC curve showed that the prediction model based on 4 key genes had high accuracy in diagnosing FGR and the area under the curve (AUC) is 0.971 (Figure 4D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1381795/fimmu-15-1381795-HTML/image_m/fimmu-15-1381795-g004.jpg)
Figure 4. Hub genes for FGR diagnosis. (A) Nomogram is constructed to predict the occurrence of FGR. (B) DCA curves. (C) the calibration curves. (D) the diagnostic efficacy verification based on the ROC curve of nomogram and 4 hub genes respectively. OR, odds ratio; CI, confidence interval; FGR, fetal growth restriction; DCA decision curve analysis; ROC, receiver operating characteristic curve; AUC, area under the curve.
3.5 GSEA analysis and immune characteristics analysis
GSEA results revealed the potential mechanisms of 4 hub genes and the details were presented in Supplementary Figures 2A–D. The results of immune cell infiltration analysis showed that the infiltration level of 22 immune cells was significantly different between the FGR group and the control group (Supplementary Figures 5A, D. Specifically, the infiltration level of T cells follicular helper (p-value = 0.016), NK cells activated (p-value < 0.001) and Macrophages M1 (p-value = 0.001) increased significantly in the FGR group. The level of infiltration of Macrophages M2 (p-value = 0.006) and Dendritic cells activated (p-value = 0.001) in placental tissues of FGR was significantly lower than that of control group (Supplementary Figures 4A–E). In order to further screen out key immune cells that may play an important role, we performed LASSO regression analysis (Supplementary Figure 5B) on 22 immune cells. B cells naive, B cells memory, T cells gamma delta, NK cells resting, NK cells activated, Macrophages M0, Macrophages M1, Macrophages M2 and Dendritic cells activated were selected as candidate key immune cells in the FGR group (Supplementary Figure 5C). The correlation between the level of infiltration of 22 types of immune cells in FGR placental tissue were also investigated and the results were shown in Supplementary Figure 5E.
3.6 Immunologic correlation analysis of TREM1
In this study, we further explored the correlation between the expression level of TREM1 and 9 key immune cells and found that TREM1 was positively correlated with T cells follicular helper (Correlation coefficients: 0.36, p-value = 0.008) and Macrophages M1 (Correlation coefficients: 0.35, p-value = 0.012); and negatively correlated with Dendritic cells activated (Correlation coefficients: -0.31, p-value = 0.026) and Macrophages M2 (Correlation coefficients: -0.32, p-value = 0.023) (Figure 5A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1381795/fimmu-15-1381795-HTML/image_m/fimmu-15-1381795-g005.jpg)
Figure 5. Immunologic and clinical correlation analysis of TREM1 and validation based on real-world data. (A) the correlation between TREM1 and immune cell infiltration in the FGR group. (B) the ROC curve of the diagnostic efficacy verification of TREM1. (C) The placental mRNA expression of TREM1 in placental tissue of normal (n=24) and FGR (n=62) fetus by real-time quantitative PCR. (D) IHC in the control group. (E) IHC in the FGR group. (F) The localization and expression of TREM1 in controls (n=6) and FGR (n=6) by IHC. The blue arrow points to extravillous trophoblast and the red arrow points to syntrophoblast. The blue spots represent the nuclei and brown-particle represents positive expression of TREM1. (G) Representative Western blot of TREM1. (H) The protein level of TREM1 in controls (n=24) and FGR (n=24) by western blot. (I) The relationship between TREM1 expression and uterine artery Doppler in the subgroup analysis. AUC, area under the curve; FGR, fetal growth restriction; IHC, immunohistochemistry; AOD, The average optical density.
3.7 Validation of TREM1 in the real-world
Samples from two groups were matched based on maternal age, gestational age at delivery, gravidity and parity and the FGR singleton placentas exhibited a statistically significant increase in mRNA expression of the TREM1 gene than the normal placentas, as quantified by RT-qPCR (p-value = 1.729e−08), which was consistent with our results above (Figure 5C). Western blot analysis further confirmed this result (p-value = 1.678e−11) (Figures 5G, H). IHC revealed localization of the TREM1 protein on syntrophoblast and extravillous trophoblast (EVT). Analysis of average optical density (AOD) demonstrated a significantly higher expression intensity of TREM-1 in the FGR group compared to the normal group (p-value = 0.044) (Figures 5D–F).
We further explored the differences in TREM1 expression in different subgroups, including age, gender of offspring, history of PE, and features of UA and UtA Doppler examination pregnancy (Supplementary Figures 6A–D). The results indicated that elevated TREM1 expression in the placenta of FGR patients with UtA spectrum abnormalities (p-value = 0.043) (Figure 5H). Therefore, patients in the FGR group were divided into three groups: normal UtA, UtA with bilateral notch, and elevated UtA-PI. Then the analysis of differences in TREM1 expression among the three groups were performed (Figure 5I). We found that the TREM1 expression in patients with elevated UtA-PI is significantly higher than that in patients with normal uterine artery spectrum (p-value = 0.029). In addition, ROC curve results indicated TREM1 likelihood as valuable biomarkers in the real world (AUC = 0.894) (Figure 5B).
4 Discussion
FGR caused by placenta dysfunction is a risk factor for multiple adverse pregnancy outcomes including fetal death (24). However, to our knowledge, there is still a lack of biomarkers that can effectively predict or treat placenta-derived FGR. Therefore, in this study, we obtained DEGs by comparing gene expression levels in placental tissues of FGR group and control group. Then, three machine learning algorithms (RF model, SVM-RFE and LASSO) were used to screen out four key genes and construct a prognostic model. Whether using supervised or unsupervised techniques, machine learning algorithms can handle complex non-linear relationships, efficiently process large-scale data, automatically identify patterns and regularities, rapidly iterate and optimize models, and provide evaluations of predictive performance and accuracy (25, 26). Nomogram is also more quantitative and intuitive, which is convenient for clinicians to use (27). TREM1 has been shown to play a key role in both innate and adapted immune responses and is considered as a potential pathogenic gene for a variety of diseases (28, 29). Therefore, we took TREM1 as the focus of follow-up studies to further explore the relationship between TREM1 expression and immune cell infiltration in the FGR group. TREM1, which was identified as the potential biomarker, is highly correlated with the degree of infiltration of four immune cells (T cells follicular helper cells, Macrophages M1, DC and Macrophages M2). In addition, we also verified its expression levels and diagnostic efficacy using real-world data and explored its relationship with various clinical features.
The RF model, SVM-RFE and LASSO were subsequently screened for four key diagnostic biomarkers (HIST1H2BB, SCD, SPINK1 and TREM1). HIST1H2BB (Histone Cluster 1H2B Family Member B) encodes histone H2B type 1-b, which is a replication-dependent histone. This gene participates in the packaging of Telomere ends and RNA polymerase I promoter opening, DNA repair, transcription regulation. DNA replication plays an important role in chromosomal stability and has been shown to be a potential biomarker for high-grade serous ovarian cancer (30). However, its implication in placental diseases remains undisclosed. SPINK1 encodes a trypsin inhibitor secreted from pancreatic acinar cells into pancreatic juice, which plays an important role in various digestive systems (31). While mRNA expressing of SPINK1 was noted to be altered in mouse placenta development, its specific function in the placenta remained unknown (32). No associations with other pregnancy complications have been reported for SPINK1. SCD (Stearoyl-CoA desaturase) encodes a critical enzyme in fatty acid metabolism, catalyzing the rate-limiting step in the formation of monounsaturated fatty acids (33, 34). A previous study revealed that the inhibition of SCD attenuated the impact of oleic acid, resulting in the downregulation of migration and proliferation in human extravillous trophoblast (EVT) cells (35). Elevated placental mRNA and protein expressions of SCD were observed in gestational diabetes mellitus (GDM) pregnancies due to promoting the synthesis of palmitic acid (PA) into palmitoleic acid (POA) with anti-inflammatory effect, suggesting a potential association between SCD and the promotion of fetal growth (36, 37). Placental inflammation is often one of the main causes of FGR (38), so we hypothesized that the up-regulation of SCD expression in the placenta of FGR may be a compensatory effect against the inflammatory response. TREM1 encodes an immunoglobulin (Ig) superfamily transmembrane protein plays a key role in regulating innated and adaptive immunity. Specifically, TREM1 can multimerizes and forms a complex with transmembrane adapter TYROBP/DAP12, a SYK-mediated cascade of tyrosine phosphorylation, activating multiple downstream mediators, such as BTK, MAPK1 and MAPK3 to promote the release of pro-inflammatory cytokines and chemokines (28, 39, 40). In addition, upregulation of TREM1 expression can activate signaling pathways including Toll-like receptor (TLR) and NOD-like receptor engagement, which play a role in immune regulation (41, 42). Recently, as more researchers have begun to explore the role TREM1 plays in obstetric disorders, TREM1 have been reported to be over-expressed in the third-trimester placenta and maternal serum of PE in several studies (43–45). FGR and PE share similar pathogenesis of inadequate placentation, inflammation, and maternal vascular dysfunction (46). Furthermore, previous study also mentioned that birth weight was negatively correlated with sTREM-1 in preterm infants’ peripheral blood (47). However, no studies have reported the association between FGR and TREM1 to authors’ knowledge. Thus, in our study, we verified the mRNA and protein expression of TREM1 in second and third trimester of FGR placentas compared to gestational age-matched controls for the first time by RT-PCR, Western blot and IHC. The mRNA and protein levels of TREM1 were relatively low and stable in controls from 20 to 38 gestational weeks, while were significantly higher in FGR placentas. We also explored that the elevated mRNA levels of TREM1 was unrelated to the pregnant complicated with PE, or the maternal age, fetal gender, and umbilical artery doppler. The sole association between TREM1 overexpression and the elevated UtA-PI, along with its localization within the syntrophoblast and EVT in IHC, suggests a potential significant involvement of TREM1 in the remodeling processes of uterine spiral arteries. Studies have explored TREM1 can be both expressed in TEV-1 cells, an EVT cell line and BeWo cells, mimicking the function and phenotype of villous trophoblast (VT) cells (48). EVTs plays a dominant role in the remodeling of spiral arteries by invading spiral artery wall and replacing the smooth muscle, and finally reducing blood vessel resistance and plasticity (49, 50). However, the study manifested that the overexpression of TREM1 gene promotes migration and invasion of TEV-1 cells through activation of the NF-κB pathway (43). BeWo is the most extensively used as a cell culture model to mimic in vivo essentialization of placental villous trophoblast (51). The study suggested that although TREM1 did not appear to play a role in BeWo cell fusion, it is required for the induction of human chorionic gonadotropin hormone (hCG), a placental-specific protein associated with syncytialization (44). Although the upregulation of TREM1 in FGR placenta and its positive regulatory effects on EVTs seem to be contradictory, the event of placentation involves an intricate coordination of multiple cells, products, formed structures and immune system (52). TREM1 serves as a unique insight into understanding of this complex coordination under physiological and pathological conditions.
The interactions between the trophoblast cells and the maternal various immune cells have an impact on the outcome of the pregnancy (53). In the third trimester placental pathology of FGR, the expression of CD68+ Macrophages was higher, while the ratio of Macrophages M2 with anti-inflammatory effect was decreased (54), which is consistent with our findings through immune characteristics analysis. According to the results of the previous studies mentioned above and combined with the results of this study, we assume that the decrease in the number of Macrophages M2 is related to the increased production of various pro-inflammatory cytokines (55, 56), which may be related to the pathogenesis of FGR. Our study also showed increased T cells follicular helper, natural killers (NK) cells activated and decreased dendritic cells (DCs) activated in FGR placenta. in pregnant women of FGR. T cells helper cells are a special type of CD4+T cells, which can play a key role in adapted immune response by promoting B cell activation and antibody production and regulating immune memory (57). And a Single-cell RNA sequencing demonstrated that Aire+ cell depletion in pregnancy, which was thought to be a sign of causing FGR during early mouse pregnancy, results in expansion of T follicular helper cells (58). Although the specific mechanism is not clear, this result suggested that T cells follicular helper cells may be one of the causes of FGR. Studies manifested higher proportion of NK cells in umbilical cord blood (59) and reduced activation of peripheral blood DCs (60) in FGR pregnancies. DC, as an antigen presenting cell (APC), plays an important role in the development of FGR. On the one hand, DC can recognize and present fetal tissue-specific antigens in the process of placenta formation, thereby inducing maternal immune tolerance to the fetus. On the other hand, in the process of fetal growth and development, DC plays an important role in maintaining the balance of maternal immune system by regulating the activation degree of various T cells to avoid adverse effects of excessive immune response on fetal growth and development (59).
Interestingly, TREM1 was significantly positively correlated with T cells follicular helper cells and Macrophages M1, and negatively correlated with DC and Macrophages M2 infiltration levels. In addition, based on LASSO regression results, four immune cells (T cells helper cells, DCs, Macrophages M1 and Macrophages M2) overlapped with characteristic immune cells in the FGR group, suggesting that TREM1 may play a role in regulating the immune response to FGR by influencing characteristic immune cells. Overall, multiple infiltrating immune cells collectively contribute to the development of FGR, with TREM1 potentially playing a crucial role in this process.
A limitation of this study is that first trimester placenta was not available to examine the localization and over-expression of TREM-1, which cannot provide more information about how TREM-1 influence the development of human placenta in early pregnancy. Furthermore, although we found that TREM-1 was over expressed in FGR with increased UtA-PI, uterine artery measurements in these clinical cases were all at an advanced stage of the disease. There is a lack of evidence regarding the relationship between UtA-PI changes and TREM1 expression in the first trimester. This will be our next research plan to explore the predictive value of TREM-1 for early screening of FGR that mediated by placental malperfusion.
This study concluded that SCD, SPINK1, TREM1 and HIST1H2BB were diagnostic indicators of FGR by using three machine learning algorithms, and a nomogram was constructed to assess the probability of each patient developing FGR. Furthermore, the results of this study suggested that a variety of immune cells may have a role in the onset and progression of FGR, and TREM1 was not only associated with four candidate key immune cells of FGR, but also correlated with increased UtA-PI, which suggested that TREM1 may be involved in remodeling processes of uterine spiral arteries through immune regulation, thus affecting the development of FGR, and provide some new clues for our future research on the prediction and treatment of FGR.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving humans were approved by Ethics Committee of Shanghai First Maternity and Infant Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
XW: Conceptualization, Data curation, Formal analysis, Investigation, Supervision, Writing – original draft, Writing – review & editing. ZL: Formal analysis, Investigation, Methodology, Software, Writing – original draft, Writing – review & editing. LC: Data curation, Investigation, Methodology, Software, Supervision, Validation, Writing – original draft, Writing – review & editing. DS: Investigation, Validation, Writing – review & editing. QS: Validation, Writing – review & editing. LZ: Validation, Writing – review & editing. FZ: Validation, Writing – review & editing. LS: Funding acquisition, Resources, Supervision, Visualization, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Key Research and Development Program of China (grant number 2022YFC2704700); the National Natural Science Foundation of China (grant number 82071656) and Shanghai Municipal Science and Technology Commission (grant number 21Y11907500; 23DZ2303400)
Acknowledgments
The authors thank the staff of the Department of Fetal Medicine & Prenatal Diagnosis Center at the Shanghai First Maternity and Infant Hospital for their diligent work and all the patients for their selfless participation.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2024.1381795/full#supplementary-material
Supplementary Figure 1 | Interaction network analysis differentially expressed genes.
Supplementary Figure 2 | Functional analysis of DEGs. (A) DO analysis. (B) GO enrichment analysis. (C) KEGG pathway enrichment analysis. (D) GSEA analysis. DEGs, differentially expressed genes. DO, disease ontology; GO, gene ontology; KEGG, Kyoto encyclopedia of genes and genomes; GSEA, gene set enrichment analysis.
Supplementary Figure 3 | The GSEA analysis. (A) SCD. (B) SPINK1. (C) TREM1. (D) HIST1H2BB.
Supplementary Figure 4 | The Violin diagram of the difference in 5 immune cells infiltration levels.
Supplementary Figure 5 | Evaluation and visualization of immune cell infiltration. (A) PCA. (B, C) The LASSO technique identifies candidate key immune cells in the FGR group. (D) Boxplot diagram of the proportion of 22 types of immune cells. (E) Correlation analysis of 22 types of immune cells. PCA, Principal component analysis; LASSO, least absolute shrinkage and selection operator; FGR, fetal growth restriction.
Supplementary Figure 6 | Differences of TREM1 expression and various clinical factors.
Supplementary Table 1 | Information on microarray datasets obtained from GEO. GEO, Gene Expression Omnibus.
Supplementary Table 2 | Detailed clinical characteristics of the FGR and control groups.
Supplementary Table 3 | Details of DEGs. DEGs, differentially expressed genes.
Supplementary Table 4 | Details of the disease enrichment analysis (DO).
Supplementary Table 5 | The OR value and 95% CI of 4 hub genes for FGR diagnosis. OR, odds ratio; CI, confidence interval.
Supplementary Table 6 | Primers of PCR.
References
1. Egaña-Ugrinovic G, Sanz-Cortes M, Figueras F, Couve-Perez C, Gratacós E. MRI insular cortical morphometry and its association with neurobehavior in late-onset small-for-gestational-age fetuses. Ultrasound obstetrics gynecology: Off J Int Soc Ultrasound Obstetrics Gynecology. (2014) 44:322–9. doi: 10.1002/uog.13360
2. Gaccioli F, Aye I, Sovio U, Charnock-Jones DS, Smith GCS. Screening for fetal growth restriction using fetal biometry combined with maternal biomarkers. Am J obstetrics gynecology. (2018) 218:S725–s37. doi: 10.1016/j.ajog.2017.12.002
3. Mylrea-Foley B, Napolitano R, Gordijn S, Wolf H, Lees CC, Stampalija T. Do differences in diagnostic criteria for late fetal growth restriction matter? Am J obstetrics gynecology MFM. (2023) 5:101117. doi: 10.1016/j.ajogmf.2023.101117
4. Gordijn SJ, Beune IM, Thilaganathan B, Papageorghiou A, Baschat AA, Baker PN, et al. Consensus definition of fetal growth restriction: a Delphi procedure. Ultrasound obstetrics gynecology: Off J Int Soc Ultrasound Obstetrics Gynecology. (2016) 48:333–9. doi: 10.1002/uog.15884
5. Hoo R, Nakimuli A, Vento-Tormo R. Innate immune mechanisms to protect against infection at the human decidual-placental interface. Front Immunol. (2020) 11:2070. doi: 10.3389/fimmu.2020.02070
6. Cohen E, Baerts W, van Bel F. Brain-sparing in intrauterine growth restriction: Considerations for the neonatologist. Neonatology. (2015) 108:269–76. doi: 10.1159/000438451
7. Figueras F, Cruz-Martinez R, Sanz-Cortes M, Arranz A, Illa M, Botet F, et al. Neurobehavioral outcomes in preterm, growth-restricted infants with and without prenatal advanced signs of brain-sparing. Ultrasound obstetrics gynecology: Off J Int Soc Ultrasound Obstetrics Gynecology. (2011) 38:288–94. doi: 10.1002/uog.9041
8. Smith GC. First-trimester determination of complications of late pregnancy. Jama. (2010) 303:561–2. doi: 10.1001/jama.2010.102
9. Knöfler M, Haider S, Saleh L, Pollheimer J, Gamage T, James J. Human placenta and trophoblast development: key molecular mechanisms and model systems. Cell Mol Life sciences: CMLS. (2019) 76:3479–96. doi: 10.1007/s00018-019-03104-6
10. Yockey LJ, Iwasaki A. Interferons and proinflammatory cytokines in pregnancy and fetal development. Immunity. (2018) 49:397–412. doi: 10.1016/j.immuni.2018.07.017
11. Jin M, Xu Q, Li J, Xu S, Tang C. Micro-RNAs in human placenta: tiny molecules, immense power. Molecules (Basel Switzerland). (2022) 27:5943–61. doi: 10.3390/molecules27185943
12. Song YP, Lv JW, Zhang ZC, Qian QH, Fan YJ, Chen DZ, et al. Effects of gestational arsenic exposures on placental and fetal development in mice: the role of cyr61 m6A. Environ Health Perspect. (2023) 131:97004. doi: 10.1289/ehp12207
13. An X, Qin J, Hu X, Zhou Y, Fu B, Wei H. Overexpression of lipocalin 2 in PBX1-deficient decidual NK cells promotes inflammation at the maternal-fetal interface. Am J Reprod Immunol (New York NY: 1989). (2023) 89:e13676. doi: 10.1111/aji.13676
14. Moffett A, Shreeve N. Local immune recognition of trophoblast in early human pregnancy: controversies and questions. Nat Rev Immunol. (2023) 23:222–35. doi: 10.1038/s41577-022-00777-2
15. Yu G, Wang LG, Yan GR, He QY. DOSE: an R/Bioconductor package for disease ontology semantic and enrichment analysis. Bioinf (Oxford England). (2015) 31:608–9. doi: 10.1093/bioinformatics/btu684
16. Friedman J, Hastie T, Tibshirani R. Regularization paths for generalized linear models via coordinate descent. J Stat software. (2010) 33:1–22. doi: 10.18637/jss.v033.i01
17. Wu Y, Seufert I, Al-Shaheri FN, Kurilov R, Bauer AS, Manoochehri M, et al. DNA-methylation signature accurately differentiates pancreatic cancer from chronic pancreatitis in tissue and plasma. Gut. (2023) 72:2344–53. doi: 10.1136/gutjnl-2023-330155
18. Noble WS. What is a support vector machine? Nat Biotechnol. (2006) 24:1565–7. doi: 10.1038/nbt1206-1565
19. Chen S, Liu Y, Yang J, Liu Q, You H, Dong Y, et al. Development and validation of a nomogram for predicting survival in male patients with breast cancer. Front Oncol. (2019) 9:361. doi: 10.3389/fonc.2019.00361
20. Wolbers M, Koller MT, Witteman JC, Steyerberg EW. Prognostic models with competing risks: methods and application to coronary risk prediction. Epidemiol (Cambridge Mass). (2009) 20:555–61. doi: 10.1097/EDE.0b013e3181a39056
21. Xiao G, Li Y, Wang M, Li X, Qin S, Sun X, et al. FBXW7 suppresses epithelial-mesenchymal transition and chemo-resistance of non-small-cell lung cancer cells by targeting snai1 for ubiquitin-dependent degradation. Cell proliferation. (2018) 51:e12473. doi: 10.1111/cpr.12473
22. Osipyants AI, Poloznikov AA, Smirnova NA, Hushpulian DM, Khristichenko AY, Chubar TA, et al. L-ascorbic acid: A true substrate for HIF prolyl hydroxylase? Biochimie. (2018) 147:46–54. doi: 10.1016/j.biochi.2017.12.011
23. Dong Z, Liang X, Zhang Q, Luo S, Liu H, Wang X, et al. Folic acid deficiency enhances the tyr705 and ser727 phosphorylation of mitochondrial STAT3 in in vivo and in vitro models of ischemic stroke. Trans stroke Res. (2021) 12:829–43. doi: 10.1007/s12975-020-00860-7
24. Graafmans WC, Richardus JH, Macfarlane A, Rebagliato M, Blondel B, Verloove-Vanhorick SP, et al. Comparability of published perinatal mortality rates in Western Europe: the quantitative impact of differences in gestational age and birthweight criteria. BJOG: an Int J obstetrics gynaecology. (2001) 108:1237–45. doi: 10.1111/j.1471-0528.2001.00291.x
25. Uddin S, Khan A, Hossain ME, Moni MA. Comparing different supervised machine learning algorithms for disease prediction. BMC Med Inf decision making. (2019) 19:281. doi: 10.1186/s12911-019-1004-8
26. Greener JG, Kandathil SM, Moffat L, Jones DT. A guide to machine learning for biologists. Nat Rev Mol Cell Biol. (2022) 23:40–55. doi: 10.1038/s41580-021-00407-0
27. Zhou X, Ning Q, Jin K, Zhang T, Ma X. Development and validation of a preoperative nomogram for predicting survival of patients with locally advanced prostate cancer after radical prostatectomy. BMC Cancer. (2020) 20:97. doi: 10.1186/s12885-020-6565-5
28. Bouchon A, Dietrich J, Colonna M. Cutting edge: inflammatory responses can be triggered by TREM-1, a novel receptor expressed on neutrophils and monocytes. J Immunol (Baltimore Md: 1950). (2000) 164:4991–5. doi: 10.4049/jimmunol.164.10.4991
29. Liu X, Chen W, Zhu G, Yang H, Li W, Luo M, et al. Single-cell RNA sequencing identifies an Il1rn(+)/Trem1(+) macrophage subpopulation as a cellular target for mitigating the progression of thoracic aortic aneurysm and dissection. Cell Discovery. (2022) 8:11. doi: 10.1038/s41421-021-00362-2
30. Valle BL, Rodriguez-Torres S, Kuhn E, Díaz-Montes T, Parrilla-Castellar E, Lawson FP, et al. HIST1H2BB and MAGI2 methylation and somatic mutations as precision medicine biomarkers for diagnosis and prognosis of high-grade serous ovarian cancer. Cancer Prev Res (Philadelphia Pa). (2020) 13:783–94. doi: 10.1158/1940-6207.Capr-19-0412
31. Qian YY, Zou WB, Chen JM, Li ZS, Liao Z. Identification of a novel SPINK1 deletion in a teenager with idiopathic chronic pancreatitis. Digestive liver disease: Off J Ital Soc Gastroenterol Ital Assoc Study Liver. (2017) 49:941–3. doi: 10.1016/j.dld.2017.04.024
32. Ishida M, Ohashi S, Kizaki Y, Naito J, Horiguchi K, Harigaya T. Expression profiling of mouse placental lactogen II and its correlative genes using a cDNA microarray analysis in the developmental mouse placenta. J Reprod Dev. (2007) 53:69–76. doi: 10.1262/jrd.18002
33. Wang J, Yu L, Schmidt RE, Su C, Huang X, Gould K, et al. Characterization of HSCD5, a novel human stearoyl-CoA desaturase unique to primates. Biochem Biophys Res Commun. (2005) 332:735–42. doi: 10.1016/j.bbrc.2005.05.013
34. Paton CM, Ntambi JM. Biochemical and physiological function of stearoyl-CoA desaturase. Am J Physiol Endocrinol Metab. (2009) 297:E28–37. doi: 10.1152/ajpendo.90897.2008
35. Yang C, Lim W, Bazer FW, Song G. Oleic acid stimulation of motility of human extravillous trophoblast cells is mediated by stearoyl-CoA desaturase-1 activity. Mol Hum Reprod. (2017) 23:755–70. doi: 10.1093/molehr/gax051
36. Balachandiran M, Bobby Z, Dorairajan G, Jacob SE, Gladwin V, Vinayagam V, et al. Placental accumulation of triacylglycerols in gestational diabetes mellitus and its association with altered fetal growth are related to the differential expressions of proteins of lipid metabolism. Exp Clin Endocrinol diabetes: Off journal German Soc Endocrinol [and] German Diabetes Assoc. (2021) 129:803–12. doi: 10.1055/a-1017-3182
37. Ferchaud-Roucher V, Barner K, Jansson T, Powell TL. Maternal obesity results in decreased syncytiotrophoblast synthesis of palmitoleic acid, a fatty acid with anti-inflammatory and insulin-sensitizing properties. FASEB journal: Off Publ Fed Am Societies Exp Biol. (2019) 33:6643–54. doi: 10.1096/fj.201802444R
38. Kim CJ, Romero R, Chaemsaithong P, Kim JS. Chronic inflammation of the placenta: definition, classification, pathogenesis, and clinical significance. Am J obstetrics gynecology. (2015) 213:S53–69. doi: 10.1016/j.ajog.2015.08.041
39. Carrasco K, Boufenzer A, Jolly L, Le Cordier H, Wang G, Heck AJ, et al. TREM-1 multimerization is essential for its activation on monocytes and neutrophils. Cell Mol Immunol. (2019) 16:460–72. doi: 10.1038/s41423-018-0003-5
40. Ormsby T, Schlecker E, Ferdin J, Tessarz AS, Angelisová P, Köprülü AD, et al. Btk is a positive regulator in the TREM-1/DAP12 signaling pathway. Blood. (2011) 118:936–45. doi: 10.1182/blood-2010-11-317016
41. Arts RJ, Joosten LA, Dinarello CA, Kullberg BJ, van der Meer JW, Netea MG. TREM-1 interaction with the LPS/TLR4 receptor complex. Eur Cytokine network. (2011) 22:11–4. doi: 10.1684/ecn.2011.0274
42. Bouchon A, Facchetti F, Weigand MA, Colonna M. TREM-1 amplifies inflammation and is a crucial mediator of septic shock. Nature. (2001) 410:1103–7. doi: 10.1038/35074114
43. Xie Y, Li X, Lv D, He M, Sun Y, Lin X, et al. TREM-1 amplifies trophoblastic inflammation via activating NF-κB pathway during preeclampsia. Placenta. (2021) 115:97–105. doi: 10.1016/j.placenta.2021.09.016
44. Lim R, Barker G, Lappas M. TREM-1 expression is increased in human placentas from severe early-onset preeclamptic pregnancies where it may be involved in syncytialization. Reprod Sci (Thousand Oaks Calif). (2014) 21:562–72. doi: 10.1177/1933719113503406
45. Turhan U, Ertaş S. A promising novel biomarker for early-onset preeclampsia: soluble trigger receptor expressed on myeloid cells-1 (sTREM-1). J maternal-fetal neonatal medicine: Off J Eur Assoc Perinatal Medicine Fed Asia Oceania Perinatal Societies Int Soc Perinatal Obstet. (2022) 35:1623–8. doi: 10.1080/14767058.2020.1846706
46. Burton GJ, Jauniaux E. Pathophysiology of placental-derived fetal growth restriction. Am J obstetrics gynecology. (2018) 218:S745–s61. doi: 10.1016/j.ajog.2017.11.577
47. Garofoli F, Borghesi A, Mazzucchelli I, Tzialla C, Di Comite A, Tinelli C, et al. Preterm newborns are provided with triggering receptor expressed on myeloid cells-1. Int J immunopathology Pharmacol. (2010) 23:1297–301. doi: 10.1177/039463201002300439
49. Lyall F. Mechanisms regulating cytotrophoblast invasion in normal pregnancy and pre-eclampsia. Aust New Z J obstetrics gynaecology. (2006) 46:266–73. doi: 10.1111/j.1479-828X.2006.00589.x
50. Burton GJ. Oxygen, the Janus gas; its effects on human placental development and function. J Anat. (2009) 215:27–35. doi: 10.1111/j.1469-7580.2008.00978.x
51. Orendi K, Gauster M, Moser G, Meiri H, Huppertz B. The choriocarcinoma cell line BeWo: syncytial fusion and expression of syncytium-specific proteins. Reprod (Cambridge England). (2010) 140:759–66. doi: 10.1530/rep-10-0221
52. Ortega MA, Fraile-Martínez O, García-Montero C, Sáez MA, Álvarez-Mon MA, Torres-Carranza D, et al. The pivotal role of the placenta in normal and pathological pregnancies: A focus on preeclampsia, fetal growth restriction, and maternal chronic venous disease. Cells. (2022) 11:568–600. doi: 10.3390/cells11030568
53. Sanguansermsri D, Pongcharoen S. Pregnancy immunology: decidual immune cells. Asian Pacific J Allergy Immunol. (2008) 26:171–81.
54. Bezemer RE, Schoots MH, Timmer A, Scherjon SA, Erwich J, van Goor H, et al. Altered levels of decidual immune cell subsets in fetal growth restriction, stillbirth, and placental pathology. Front Immunol. (2020) 11:1898. doi: 10.3389/fimmu.2020.01898
55. Deer E, Herrock O, Campbell N, Cornelius D, Fitzgerald S, Amaral LM, et al. The role of immune cells and mediators in preeclampsia. Nat Rev Nephrol. (2023) 19:257–70. doi: 10.1038/s41581-022-00670-0
56. Schonkeren D, van der Hoorn ML, Khedoe P, Swings G, van Beelen E, Claas F, et al. Differential distribution and phenotype of decidual macrophages in preeclamptic versus control pregnancies. Am J Pathol. (2011) 178:709–17. doi: 10.1016/j.ajpath.2010.10.011
57. Walker LSK. The link between circulating follicular helper T cells and autoimmunity. Nat Rev Immunol. (2022) 22:567–75. doi: 10.1038/s41577-022-00693-5
58. Gillis-Buck E, Miller H, Sirota M, Sanders SJ, Ntranos V, Anderson MS, et al. Extrathymic Aire-expressing cells support maternal-fetal tolerance. Sci Immunol. (2021) 6:eabf1968–1992. doi: 10.1126/sciimmunol.abf1968
59. Huang L, Li P, Feng T, Xiong F. Changes of dendritic cell and natural killer cell on the cord blood with idiopathic fetal growth restriction. J maternal-fetal neonatal medicine: Off J Eur Assoc Perinatal Medicine Fed Asia Oceania Perinatal Societies Int Soc Perinatal Obstet. (2022) 35:7526–31. doi: 10.1080/14767058.2021.1951214
Keywords: fetal growth restriction, machine-learning, immune cell infiltration, placenta, TREM1 (triggering receptor expressed on myeloid cells)
Citation: Wei X, Liu Z, Cai L, Shi D, Sun Q, Zhang L, Zhou F and Sun L (2024) Integrated transcriptomic analysis and machine learning for characterizing diagnostic biomarkers and immune cell infiltration in fetal growth restriction. Front. Immunol. 15:1381795. doi: 10.3389/fimmu.2024.1381795
Received: 16 February 2024; Accepted: 20 August 2024;
Published: 04 September 2024.
Edited by:
Jelmer R. Prins, University Medical Center Groningen, NetherlandsReviewed by:
Maulikkumar P. Patel, Washington University in St. Louis, United StatesSemanti Ghosh, Swami Vivekananda University, India
Gunjan Saini, Purdue University, United States
Copyright © 2024 Wei, Liu, Cai, Shi, Sun, Zhang, Zhou and Sun. 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: Luming Sun, luming_sun@163.com
†These authors have contributed equally to this work