Home Potential biomarkers and molecular mechanisms in preeclampsia progression
Article Open Access

Potential biomarkers and molecular mechanisms in preeclampsia progression

  • Guohua Li , Shijia Huang , Xiaosong Liu and Qiaoling Du EMAIL logo
Published/Copyright: May 18, 2022

Abstract

This study aimed to explore potential biomarkers and molecular mechanisms in preeclampsia (PE) progression. Gene expression profiles of GSE147776 and GSE96984 were downloaded, followed by the identification of common differentially expressed genes (co-DEGs) and common differentially expressed lncRNAs (co-DElncRNAs) in PE patients between the two datasets. Key genes were identified using gene set enrichment analysis (GSEA), followed by functional enrichment analyses. Subsequently, the miRNAs of key genes and miRNA-related lncRNAs were predicted, followed by the construction of the lncRNA–miRNA–gene ceRNA network. Furthermore, the key genes associated with different gestational stages were identified. As a result, 192 co-DEGs and 16 co-DElncRNAs were revealed from the two datasets. Based on two outstanding PE-associated pathways, including glaucoma and PE, identified by GSEA, ten key genes, including IGFBP1, CORIN, and C3, were revealed. Key genes, including IL1A and IL1B, were enriched in the developmental process involved in reproduction. Furthermore, ceRNAs, such as LINC00473-miR-4476-IL1A, LINC00473-miR-1291-IL1B, and NAV2-AS4-miR-6131-REN, were identified. Moreover, REN expression was significantly upregulated in the first- and second-trimester placentae compared to C-section-term placentae. In conclusion, these key genes may serve as novel biomarkers for PE. The detection of REN expression may help in the early prediction of PE and the initiation of prophylactic medical treatment.

1 Introduction

Preeclampsia (PE) is a relatively common pregnancy disorder originating in the placenta [1]. The incidence rate of PE, which is directly related to the mortality of women and fetuses, is approximately 5–8% [2]. However, in clinical diagnosis, patients with PE are easily misdiagnosed because the clinical manifestations of PE overlap with other diseases, including hypertension, chronic kidney disease, and primary epilepsy [3]. Therefore, accurate diagnosis and timely intervention are crucial for the clinical treatment of PE.

Factually, gene expression in different molecular patterns is closely associated with PE development [4]. Bioinformatics analysis-based gene profile studies have aided the interpretation of PE progression and can be conducted to identify novel biomarkers of PE [5]. A previous study has showed that certain differentially expressed genes (DEGs), such as cytochrome P450 family 11 subfamily A member 1 (CYP11A1), participate in the pathogenesis of PE by inducing excessive autophagy [6]. Moreover, key genes such as glutathione S-transferase omega 1 (GSTO1) are considered novel biomarkers for the therapy of PE [7]. Additionally, non-coding RNAs constitute the largest group of RNA transcripts, accounting for 90% of the human genome. Aberrant regulation of specific non-coding RNAs can affect critical mechanisms in PE development, such as immunity, trophoblast proliferation, and trophoblast invasion [8]. MicroRNAs (miRNAs) and long non-coding RNAs (lncRNAs) have been reported to play crucial roles in PE development [9,10,11]. Particularly, lncRNAs may function as competing for endogenous RNAs (ceRNAs) to sponge their target miRNAs, thus mediating the post-transcriptional regulation of target genes by affecting the regulatory function of miRNAs [12]. Identification of lncRNA–miRNA–mRNA ceRNA networks may help elucidate the potential regulatory mechanisms in PE and may aid the discovery of promising biomarkers for PE [13]. Although the relationship between the molecular mechanism and pathological process of PE has been widely studied, its exact pathogenesis and the identity of associated potential diagnostic biomarkers remain uncertain.

A previous placental microarray profiling (GSE147776) study reported by Medina-Bastidas et al. [14] revealed the existence of several DEGs in pregnant women with PE and normal pregnancy (NP). However, the mechanism by which these genes interact with each other and the pathways involved in PE progression remain unknown. In this study, the common DEGs (co-DEGs) in the PE vs NP group between the GSE147776 dataset reported by Medina-Bastidas et al. and another dataset, GSE96984, were investigated, and key gene exploration based on the results of gene set enrichment analysis (GSEA) was performed. The significant functions and pathways enriched by these key genes were identified. Moreover, the common differentially expressed lncRNAs (co-DElncRNAs) in the PE vs NP group between the GSE147776 and GSE96984 datasets were identified. The miRNAs of key genes and miRNA-related lncRNAs were predicted, followed by the construction of the lncRNA–miRNA–gene ceRNA network via integration of the lncRNA–miRNA regulatory relationships and miRNA–target gene pairs. Finally, the key genes associated with different gestational stages were identified to determine the effect of the trimester on PE. The findings of this study may help provide insights into novel biomarkers and molecular mechanisms associated with PE. The flow chart of the present study is shown in Figure 1.

Figure 1 
               Experimental flow chart of the present study.
Figure 1

Experimental flow chart of the present study.

2 Materials and methods

2.1 Microarray data and preprocessing

The gene expression profile GSE147776 [14] was downloaded from the Gene Expression Omnibus (GEO) database. This dataset included seven tissue samples from patients with intrauterine growth restriction, seven from patients with PE, six from patients with PE and intrauterine growth restriction, and eight from patients with NP. The samples from patients with intrauterine growth restriction were excluded. Thus, seven tissue samples obtained from patients with PE (PE group) and eight from women with NP (NP group) were included in the training set in this analysis. The gene expression profile GSE96984, also downloaded from the GEO database, was used as a validation set in this study. This set contained data on three tissue samples obtained from PE patients (PE group) and four tissue samples obtained from women with normal pregnancy (NP group). Furthermore, using trimester and PE as keywords, only one gene expression dataset, GSE9984, was screened from the GEO database. This dataset included gene expression profiles of three gestational stages (first trimester [6 weeks + 3 days] − [8 weeks + 3 days]), second trimester ([15 weeks + 4 days] − [16 weeks + 3 days]), and C-section-term placentae). There were four placental tissue samples at each gestational stage.

The processed gene expression matrix files in the GSE147776, GSE96984, and GSE9984 datasets were downloaded. The probe IDs were converted to gene symbols based on annotation files downloaded from different platforms. Probes that did not correspond to the gene symbols were discarded. For different probes matched to one gene, the average values of the different probes were used as the final expression values.

2.2 DEG investigation and VENN plot analysis

The classic Bayesian method provided by the limma package (version 3.10.3) [15] was used to reveal the DEGs between the PE and NP groups from the GSE147776 and GSE96984 datasets, respectively. Moreover, the DElncRNAs between the PE and NP groups derived from the two datasets were identified using the same method. The threshold values for selecting DEGs and DElncRNAs were set as P-value < 0.05 and |log fold change (FC)| > 1.5. VENN plot analyses were subsequently performed to identify the co-DEGs and co-DElncRNAs exhibiting the same up- or downregulation of expression in the two datasets. Co-DEGs and co-DElncRNAs were included in the subsequent analysis.

2.3 Differential co-expression correlation analysis

The differential co-expression correlation of co-DEGs was revealed using the Diffcorr package (version: 0.4.1, https://CRAN.R-project.org/package = DiffCorr) [16] in R software. Briefly, the Pearson correlation coefficient of each co-DEG was calculated, followed by a correlation test. The co-DEGs with |r (r 1r 2)| > 1 (with r 1 being the correlation coefficient of two genes in the NP group and r 2 being the correlation coefficient of two genes in the PE group) were considered as differentially co-expressed genes. The co-expression relations in both NP and PE groups were further explored with |r| > 0.8, and P < 0.05. Finally, the co-expression network in the two groups was visualized using the Cytoscape software (version 3.4.0) [17], followed by the conduction of network topology (closeness centrality and degree centrality) analysis.

2.4 GSEA based on co-DEGs and information in the disease database

GSEA is a commonly used method for analyzing gene sets generated via the conduction of genome-wide experiments [18]. It can be used to identify significant differences in gene expression between the two different biological states [19]. Based on the gene–disease associations obtained from the DisGeNET database [20], GSEA was performed in an orderly manner (according to the descending order of the log FC) for all genes in GSE147776 using the clusterProfiler package (version: 3.16.0, http://bioconductor.org/packages/release/bioc/html/clusterProfiler.html) [21]. The Benjamini & Hochberg-adjusted P < 0.05 was considered as the threshold value for conducting analysis in the present study.

2.5 Investigation of key genes of PE

The key genes of PE were investigated based on the results of GSEA and co-expression correlation analysis. Briefly, the genes that were reported in the database of PE-associated disease (revealed by GSEA), along with the exhibition of co-expression in co-expression correlation analysis, were considered key PE genes. The differential expression of these key genes in different clinical phenotypes was also analyzed. Finally, the results were visualized and illustrated using a box diagram.

2.6 Enrichment analysis of key genes

Gene Ontology-Biological Process (GO-BP) [22] and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway [23] enrichment analyses of key genes were performed using the Metascape software (parameter: minimum overlap = 3; minimum enrichment = 1.5) [24]. P < 0.05 was considered the cutoff value to obtain significant enrichment results. Moreover, clustering analysis was conducted according to the similarity of genes enriched in each term (similarity >0.3). The most statistically significant term (minimum P-value) in each cluster was selected to define the cluster. Finally, the network was constructed using the Cytoscape software (version 3.4.0) [17].

2.7 Prediction of miRNAs of key genes and construction of the lncRNA–miRNA–gene ceRNA network

The miRNAs of key genes were predicted using miRWalk 3.0 (http://mirwalk.umm.uni-heidelberg.de/) [25]. The obtained miRNA–target gene pairs were further subjected to screening using miRDB (http://mirdb.org/) [26]. The miRNA–target gene pairs with scores > 0.95 were selected for network construction. Additionally, the miRNA-related lncRNAs were predicted using the Prediction Module of DIANA-LncBase v.2 (http://carolina.imis.athena-innovation.gr/diana_tools/web/index.php?r=lncbasev2%2Findex) [27]. The lncRNA–miRNA regulatory relationships with scores > 0.7 were then obtained. The lncRNA–miRNA regulatory relationships containing co-DElncRNAs were further screened. By integrating the lncRNA–miRNA regulatory relationships and miRNA–target gene pairs, the lncRNA–miRNA–gene ceRNA network was constructed.

2.8 Identification of key genes associated with different gestational stages

To determine the effect of the trimester (different gestational ages) on PE, based on the gene expression data of three gestational stages ([6 weeks + 3 days] − [8 weeks + 3 days], second trimester [15 weeks + 4 days] − [16 weeks + 3 days], and C-section-term placentae) in GSE9984 datasets, the expression of key genes between any two gestational stages was compared using the t-test in R.

3 Results

3.1 DEG and co-DEG investigation

A total of 23,825 genes and 4,720 lncRNAs, and 18,826 genes and 19,741 lncRNAs, were annotated in GSE147776 and GSE96984, respectively. Based on principal component analysis (PCA) results, three samples in GSE147776, GSM4445727, GSM4445714, and GSM4445724, were excluded due to abnormal expression observed between the PE and NP groups (Figure A1). All samples in GSE96984 were included for DEG analysis as no abnormal samples were observed using PCA. A total of 554 DEGs (154 downregulated and 390 upregulated) and 968 DEGs (363 downregulated and 605 upregulated) between the PE and NP groups were obtained in GSE147776 and GSE96984, respectively. The heatmap of the clustering analysis conducted for DEGs indicated that the samples in the PE and NP groups could be clearly distinguished in GSE147776 (Figure 2a) and GSE96984 (Figure 2b). The volcano plot depicts the DEGs in GSE147776 (Figure 2c) and GSE96984 (Figure 2d). Additionally, 53 DElncRNAs (27 downregulated and 26 upregulated) and 571 DElncRNAs (378 downregulated and 193 upregulated) were identified between the PE and NP groups in GSE147776 and GSE96984, respectively. Heatmap and volcano plots for DElncRNAs in the two datasets are illustrated in Figure 2e–h.

Figure 2 
                  Results of DEGs and differentially expressed lncRNAs (DElncRNAs) obtained during investigation. (a) The heat map illustrating DEGs between the PE and NP groups in GSE147776; (b) the heat map illustrating DEGs between PE group and NP group in GSE96984. (c) The volcano plot illustrating the expression of genes in GSE147776; (d) the volcano plot illustrating the expression of genes in GSE96984; the green node represents downregulated genes, while the red node represents the upregulated genes; (e) the heat map illustrating DElncRNAs between the PE and NP groups in GSE147776; (f) the heat map illustrating DElncRNAs between the PE and NP groups in GSE96984. (g) The volcano plot illustrating the expression of DElncRNAs in GSE147776; (h) the volcano plot illustrating the expression of DElncRNAs in GSE96984; the green node represents downregulated lncRNAs, while the red node represents upregulated lncRNAs.
Figure 2

Results of DEGs and differentially expressed lncRNAs (DElncRNAs) obtained during investigation. (a) The heat map illustrating DEGs between the PE and NP groups in GSE147776; (b) the heat map illustrating DEGs between PE group and NP group in GSE96984. (c) The volcano plot illustrating the expression of genes in GSE147776; (d) the volcano plot illustrating the expression of genes in GSE96984; the green node represents downregulated genes, while the red node represents the upregulated genes; (e) the heat map illustrating DElncRNAs between the PE and NP groups in GSE147776; (f) the heat map illustrating DElncRNAs between the PE and NP groups in GSE96984. (g) The volcano plot illustrating the expression of DElncRNAs in GSE147776; (h) the volcano plot illustrating the expression of DElncRNAs in GSE96984; the green node represents downregulated lncRNAs, while the red node represents upregulated lncRNAs.

Furthermore, the co-DEGs and co-DElncRNAs between the GSE147776 and GSE96984 datasets were explored via VENN plot analysis. The results showed that there were 192 co-DEGs (27 downregulated and 165 upregulated) (Figure 3a) and 16 co-DElncRNAs (7 downregulated and 9 upregulated) (Figure 3b) between the two datasets.

Figure 3 
                  VENN plot illustrating common DEGs (co-DEGs) and common DElncRNAs (co-DElncRNAs) between GSE96984 and GSE147776 datasets. (a) The common downregulated and upregulated genes observed between the two datasets. (b) The common downregulated and upregulated lncRNAs observed between the two datasets.
Figure 3

VENN plot illustrating common DEGs (co-DEGs) and common DElncRNAs (co-DElncRNAs) between GSE96984 and GSE147776 datasets. (a) The common downregulated and upregulated genes observed between the two datasets. (b) The common downregulated and upregulated lncRNAs observed between the two datasets.

3.2 Differential co-expression correlation and co-expression network analysis

A total of 1,365 differential co-expression correlations were further investigated based on 192 co-DEGs using the Diffcorr package in R software. The co-expression networks in both the NP and PE groups were constructed using 1,365 differentially co-expressed elements. The results showed that there were 166 co-DEGs and 253 correlations in the co-expression network of the NP group (Figure 4a). A total of 158 co-DEGs and 891 correlations were observed in the co-expression network of the PE group (Figure 4b). The topological property analysis of the two networks showed that, compared with the NP group, the degree centrality and closeness centrality of the PE group were significantly increased (Figure 4c, P < 0.05), which further indicated that co-expression of these genes might contribute to the occurrence of diseases.

Figure 4 
                  Co-expression networks constructed by co-DEGs. (a) The co-expression network constructed by co-DEGs in the NP group. (b) The co-expression network constructed by co-DEGs in the PE group. The red node represents upregulated genes; the blue node represents downregulated genes; the blue line represents negative co-expression; the red line represents positive co-expression. (c) The degree centrality and closeness centrality between NP and PE groups. The red block represents the NP group, while the blue block represents the PE group.
Figure 4

Co-expression networks constructed by co-DEGs. (a) The co-expression network constructed by co-DEGs in the NP group. (b) The co-expression network constructed by co-DEGs in the PE group. The red node represents upregulated genes; the blue node represents downregulated genes; the blue line represents negative co-expression; the red line represents positive co-expression. (c) The degree centrality and closeness centrality between NP and PE groups. The red block represents the NP group, while the blue block represents the PE group.

3.3 Investigation of key genes based on GSEA

Based on the DisGeNET database and logFC-sequenced genes in the GSE147776 dataset, a total of 14 PE-related diseases, such as maternal hypertension, glaucoma, and PE, were revealed (Table A1). GSEA analysis showed that co-DEGs were mainly enriched in glaucoma (C0017601) and PE (C0032914) (Table A1 and Figure A2). During the analysis of enrichment of co-DEGs in the two pathways, TNF receptor superfamily member 8 (TNFRSF8), norrin cystine knot growth factor NDP (NDP), complement C3 (C3), WT1, insulin-like growth factor binding protein 1 (IGFBP1), interleukin-1A (IL1A), corin/serine peptidase (CORIN), interleukin-1B (IL1B), renin (REN), and carbonic anhydrase 12 (CA12) were identified as key genes of PE, as the expression levels of the indicated genes were significantly upregulated in the PE group (Figure 5).

Figure 5 
                  The expression of ten key genes in the NP and PE groups. The red block represents the NP group, while the blue block represents the PE group.
Figure 5

The expression of ten key genes in the NP and PE groups. The red block represents the NP group, while the blue block represents the PE group.

3.4 GO and KEGG enrichment analyses of key genes

GO-BP and KEGG analyses were performed for the indicated ten key genes. The results showed that the genes were mainly assembled in 43 GO-BP functions and 4 pathways. Clustering analysis according to the observation of similarity of key genes enriched in each term was subsequently performed. The results showed that a total of five clusters, including positive regulation of vascular endothelial growth factor production (GO: 0010575; genes: C3, IL1A, and IL1B), developmental process involved in reproduction (GO: 0003006; genes: IL1A, IL1B, NDP, REN, and WT1), endocrine process (GO: 0050886; genes: IL1B, REN, and CORIN), response to organophosphorus (GO: 0046683; genes: IL1B, REN, and WT1), and visual system development (GO: 0150063; genes: C3, NDP, and WT1) were explored.

3.5 Construction of the lncRNA–miRNA–gene ceRNA network

A total of 99 miRNA–target gene pairs were obtained, including 97 miRNAs and 8 key genes. Data on 6 co-DElncRNAs and 26 miRNAs were obtained by further predicting lncRNA–miRNA regulatory relationships. The lncRNA–miRNA–gene ceRNA network was constructed by integrating the lncRNA–miRNA regulatory relationships and miRNA–target gene pairs, including 58 interactions, 26 miRNAs, 6 lncRNAs, and 7 key genes (Figure 6). Notably, among the six lncRNAs, LINC00473 could target most miRNAs, followed by NAV2-AS4. LINC00473-miR-4476-IL1A, LINC00473-miR-1291-IL1B, and NAV2-AS4-miR-6131-REN were identified as key ceRNAs. The complementary sequences of LINC00473-miR-4476, LINC00473-miR-1291, and NAV2-AS4-miR-6131 are shown in Figure A3.

Figure 6 
                  The constructed lncRNA–miRNA–gene ceRNA network. The red diamond node represents the upregulated gene, the orange square node represents miRNA, the purple inverted triangle node represents the upregulated lncRNA, and the green inverted triangle node represents the downregulated lncRNA.
Figure 6

The constructed lncRNA–miRNA–gene ceRNA network. The red diamond node represents the upregulated gene, the orange square node represents miRNA, the purple inverted triangle node represents the upregulated lncRNA, and the green inverted triangle node represents the downregulated lncRNA.

3.6 Identification of key genes associated with different gestational stages

The expression of key genes at different gestational stages was analyzed based on the gene expression data of three gestational stages in the GSE9984 datasets. The results showed that only REN was differentially expressed among the different gestational stages. Notably, the expression of REN gradually decreased with an increase in gestation time. REN expression in the first- and second-trimester placentae was significantly higher than that in the C-section-term placentae (Figure 7). These data indicate that the detection of REN expression may be useful for the early prediction of PE.

Figure 7 
                  The expression of REN at different gestational stages. The red block represents the first-trimester group, the green block represents the second-trimester group, and the blue block represents the C-section term group.
Figure 7

The expression of REN at different gestational stages. The red block represents the first-trimester group, the green block represents the second-trimester group, and the blue block represents the C-section term group.

4 Discussion

Although PE is a pregnancy-associated complication exhibiting considerable consequences in women, the biomarkers and associated molecular mechanisms underlying PE progression remain unclear. The present bioinformatics study revealed two remarkable pathways (glaucoma and PE) using GSEA based on 14 PE-related diseases and co-DEGs. A total of ten co-DEGs, including IGFBP1, CORIN C3, IL1A, IL1B, and REN, were identified as key genes. Moreover, functional cluster analysis results showed that these key genes, including IL1A and IL1B, were mainly involved in functions involved in reproduction, such as developmental processes. Furthermore, LINC00473-miR-4476-IL1A, LINC00473-miR-1291-IL1B, and NAV2-AS4-miR-6131-REN were identified as key ceRNAs in the lncRNA–miRNA–gene ceRNA network. Additionally, we found that REN expression was significantly upregulated in first- and second-trimester placentae compared to that in the C-section-term placentae, indicating that the detection of REN expression might be useful for early prediction of PE.

As PE is a common vascular disease associated with pregnancy with a known genetic predisposition, the abnormal expression of certain genes during this process contributes to the progression of PE [28]. IGFBP1 plays essential role in embryogenesis and carcinogenesis [29]. Under normal circumstances of clinical observation, the expression of IGFBP1 in pregnant women is reportedly higher than that in non-pregnant women [30]. However, low expression of IGFBP1 has been observed in the maternal blood samples derived from PE patients compared to healthy women [31]. In contrast, an investigation based on the alteration of serum IGFBP-1 concentration in women with different stages of eclampsia showed that IGFBP-1 levels in the PE group were significantly higher than those in the control group, which was related to the severity of PE [32]. It has been proven that the expression of the IGFBP1 gene is abnormal before PE development and can be clinically detected due to abnormalities in the trophoblastic invasion, thereby highlighting potential applicability as a marker for the identification of PE [33]. Moreover, CORIN encodes an atrial natriuretic peptide-converting enzyme [34]. The CORIN-mediated atrial natriuretic peptide in the uterus under pregnancy conditions promotes spiral artery remodeling and trophoblast invasion [35]. A previous study showed that the differential expression of the CORIN gene induced by two mutations (K317E and S472g) could be observed in PE patients, which might contribute to the development of PE [36]. The potential consideration of CORIN in PE clinical treatment as a new biomarker has been investigated in a previous study [37]. Furthermore, C3 is a complement component with the highest level in serum and plays an important role in both the classical and bypass activation pathways [38]. A previous genetic case–control study based on 960 PE patients has revealed that C3 plays a potential role in the complement system in the pathogenesis of PE, and may be considered a prognostic and therapeutic biomarker of preeclamptic women [39]. It has been proven that the gene mutation of C3 (such as C3F) is associated with susceptibility to PE [40]. In the present study, IGFBP-1, CORIN, and C3 were identified as three key genes with upregulated expression between the PE and NP groups. Meanwhile, GSEA showed that all three genes were significantly enriched in pathways such as those in PE. Thus, we speculated that expression levels of IGFBP-1, CORIN, and C3 were upregulated during the development of PE, and the three genes might therefore be considered novel biomarkers for PE.

The occurrence and development of PE is closely associated with the functions of the human reproductive system [41]. The occurrence of PE represents a reproductive disadvantage unique to humans compared to other mammals [42]. Reproduction and associated factors, such as neutrophil extracellular traps, play a vital role in the development of pregnancy-related disorders in PE [43]. A previous study has indicated that the differential expression of certain genes in the reproduction pathway influences the risk of PE [44]. The interleukin-1 (IL1) family is a group of 11 cytokines (including IL1A and IL1B) that form a complex network of proinflammatory cytokines that initiate and control inflammatory responses through integrin expression on leukocytes and endothelial cells [45]. IL-1 system components, including IL1A and IL1B, have been shown to exhibit several sites of synthesis in the reproductive system, including the ovary [46]. In fact, the genetic variations of the IL1 family member IL1A are related to PE risk in the Chinese Han population [47]. A previous study showed that IL1A could be considered a therapeutic tool for monitoring PE patients because it leads to a proinflammatory cascade and results in high levels of circulating cytokines involved in PE progression [48]. Moreover, univariate analysis results suggest that the inflammatory process-related gene IL1B is involved in the host response of PE [49]. In this study, the IL-1 family members IL1A and IL1B were not only revealed as upregulated key genes in the PE and NP groups but also found to be significantly assembled during the developmental process involved in reproduction. Thus, we speculated that IL1A and IL1B might participate in the progression of PE by affecting developmental processes involved in reproduction. Additionally, LINC00473-miR-4476-IL1A and LINC00473-miR-1291-IL1B were identified as key ceRNAs in the lncRNA–miRNA–gene ceRNA network. LINC00473 reportedly contributes to PE development and may serve as a potential biomarker for PE [50]. Although miR-4476 and miR-1291 have been revealed to play a role in the development of various cancers [51,52,53], their possible role in PE development has not been reported and warrants further study. Our results suggest that these ceRNAs may be key mechanisms mediating PE development and may provide novel biomarkers for PE.

Notably, the first-trimester prediction of PE has important clinical implications, as it allows clinicians to focus on high-risk populations and aids the development of prophylactic medical treatment strategies. Herein, we found that the expression of REN in first and second-trimester placentae was significantly higher than that in the C-section-term placentae. Previous studies have confirmed that REN is associated with PE, and REN gene polymorphisms play a crucial role in PE development [54,55]. Moreover, NAV2-AS4-miR-6131-REN was identified as a key ceRNA in our study. However, the roles of NAV2-AS4 and miR-6131 in PE have not been elucidated. Considering the potential role of REN in PE, we speculate that detection of REN expression may demonstrate utility for early prediction of PE, and NAV2-AS4 and miR-6131 may also be involved in PE via regulation of REN expression.

However, there were some limitations in the present study. First, the sample size was small. Second, there was a lack of clinical verification. Finally, whether the selected lncRNA–miRNA–mRNA was dysregulated under hypoxic conditions was not explored in trophoblast cells in vitro. Thus, further experimental verification based on a larger sample size is imperative.

In conclusion, expression levels of key genes, such as IGFBP1, CORIN, C3, IL1A, IL1B, and REN, were upregulated during the development of PE, and these genes could therefore be considered novel biomarkers for PE. Moreover, IL1A and IL1B may participate in the progression of PE via developmental processes involved in reproduction. Furthermore, key ceRNAs, such as LINC00473-miR-4476-IL1A, LINC00473-miR-1291-IL1B, and NAV2-AS4-miR-6131-REN, may also contribute to PE development. Detection of REN expression may be useful for the early prediction of PE and may aid the initiation of prophylactic medical treatment.


tel: +86-21-50730190

  1. Funding information: This study was supported by the Technology Project of the Shanghai Pudong New District Health and Family Planning Commission (No. PW2019D-9).

  2. Conflict of interest: Authors state no conflict of interest.

  3. Data availability statement: The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Appendix

Figure A1 
                  The principal component analysis (PCA) results for GSE147776. (a) Results of PCA before abnormal sample exclusion; (b) results of PCA after abnormal sample exclusion.
Figure A1

The principal component analysis (PCA) results for GSE147776. (a) Results of PCA before abnormal sample exclusion; (b) results of PCA after abnormal sample exclusion.

Figure A2 
                  The results of gene set enrichment analysis performed for two significantly enriched pathways, including glaucoma and preeclampsia.
Figure A2

The results of gene set enrichment analysis performed for two significantly enriched pathways, including glaucoma and preeclampsia.

Figure A3 
                  The complementary sequences of LINC00473-miR-4476, LINC00473-miR-1291, and NAV2-AS4-miR-6131.
Figure A3

The complementary sequences of LINC00473-miR-4476, LINC00473-miR-1291, and NAV2-AS4-miR-6131.

Table A1

The enrichment result using GSEA

ID Description SetSize EnrichmentScore P. adjust CoreEnrichment
C0017601 Glaucoma 357 0.384723753 0.00401984 NDP/WT1/TTR/IL1B/CD69/CA12/IL1A/GSTT1/ROM1/NOS2/HLA-DQB1/IFT172/PTPN22/FN1/PTGER3/DAPK2/MMP3/TUB/GJB3/NOD2/TNF/SAA2/CA2/NTRK1/ANTXR1/PITX2/HLA-DRB1/TTC8/AGT/SFRP1/DDX58/ACE2/MGP/OPTN/PLK4/ARL6/COCH/NCAM1/VEGFA/EDNRB/VCAN/TMEM98/MMP1/GSTM1/PUS1/SALL1/BEST1/ABCA4/FAT4/TMEM182/COL4A1/FAS/CA4/GSTM3/IFIH1/NPPA/ADAMTS17/EDNRA/TLR4/BCOR/RPGR/ADAMTS10/WDR36/LTBP2/CDKN2A/CLU/IDUA/SPATA7/CTGF/ADRB2/PIK3R1/SPARC/IL6/TMCO1/NFKB1/PRPH2/TDRD7/CYP1B1/C1QA/BBS2/SLC4A4/TGFB2/GSTM2/RPGRIP1
C0032914 Pre-Eclampsia 146 0.465393454 0.00401984 IGFBP1/TAC3/C3/TNFRSF8/CXCL10/CD40LG/IL15/IFNG/REN/SERPIND1/NOS2/CORIN/CCL5/MMP9/SERPINE1/SELL/TAP1/HDC/NOD2/TNF/ADA/BDKRB2/VCAM1/S100B/AGT/PYGM/VEGFA/NOV/IL18/LGALS1/TAP2/PRCP/CYR61/HP/HLA-G

Notes: GSEA, gene set enrichment analysis; EnrichmentScore: positive values indicate up-regulated genes and negative values indicate down regulated genes; Core Enrichment: the core genes with the most contribution to the enrichment score.

References

[1] Phipps E, Prasanna D, Brima W, Jim B. Preeclampsia: updates in pathogenesis, definitions, and guidelines. Clin J Am Soc Nephrology. 2016;11(6):1102–13.10.2215/CJN.12081115Search in Google Scholar PubMed PubMed Central

[2] Jim B, Karumanchi SA. Preeclampsia: pathogenesis, prevention, and long-term complications. Seminars in nephrology. 2017;37(4):386–97.10.1016/j.semnephrol.2017.05.011Search in Google Scholar PubMed

[3] Rana S, Lemoine E, Granger JP, Karumanchi SA. Preeclampsia: pathophysiology, challenges, and perspectives. Circulation Res. 2019;124(7):1094–112.10.1161/CIRCRESAHA.118.313276Search in Google Scholar PubMed

[4] Liang M, Niu J, Zhang L, Deng H, Ma J, Zhou W, et al. Gene expression profiling reveals different molecular patterns in G-protein coupled receptor signaling pathways between early-and late-onset preeclampsia. Placenta. 2016;40:52–9.10.1016/j.placenta.2016.02.015Search in Google Scholar PubMed

[5] Luo S, Cao N, Tang Y, Gu W. Identification of key microRNAs and genes in preeclampsia by bioinformatics analysis. PLoS One. 2017;12(6):e0178549.10.1371/journal.pone.0178549Search in Google Scholar PubMed PubMed Central

[6] Pan T, He G, Chen M, Bao C, Chen Y, Liu G, et al. Abnormal CYP11A1 gene expression induces excessive autophagy, contributing to the pathogenesis of preeclampsia. Oncotarget. 2017;8(52):89824–36.10.18632/oncotarget.21158Search in Google Scholar PubMed PubMed Central

[7] Li X, Fang Y. Bioinformatics identification of potential genes and pathways in preeclampsia based on functional gene set enrichment analyses. Exp therapeutic Med. 2019;18(3):1837–44.10.3892/etm.2019.7749Search in Google Scholar PubMed PubMed Central

[8] Munjas J, Sopić M, Stefanović A, Košir R, Ninić A, Joksić I, et al. Non-coding RNAs in preeclampsia—molecular mechanisms and diagnostic potential. Int J Mol Sci. 2021;22(19):10652.10.3390/ijms221910652Search in Google Scholar PubMed PubMed Central

[9] Selvaraj S, Lakshmanan G, Kalimuthu K, Sekar D. Role of microRNAs and their involvement in preeclampsia. Epigenomics. 2020 Oct;12(20):1765–7. 10.2217/epi-2020-0281. [Epub 2020 Oct 20].Search in Google Scholar PubMed

[10] Hornakova A, Kolkova Z, Holubekova V, Loderer D, Lasabova Z, Biringer K, et al. Diagnostic Potential of MicroRNAs as Biomarkers in the Detection of Preeclampsia. Genet Test Mol Biomarkers. 2020;24(6):321–7.10.1089/gtmb.2019.0264Search in Google Scholar PubMed

[11] Moradi MT, Rahimi Z, Vaisi-Raygani A. New insight into the role of long non-coding RNAs in the pathogenesis of preeclampsia. Hypertens Pregnancy. 2019;38(1):41–51.10.1080/10641955.2019.1573252Search in Google Scholar PubMed

[12] Smillie CL, Sirey T, Ponting CP. Complexities of post-transcriptional regulation and the modeling of ceRNA crosstalk. Crit Rev Biochem Mol Biol. 2018;53(3):231–45.10.1080/10409238.2018.1447542Search in Google Scholar PubMed PubMed Central

[13] Jiang R, Wang T, Zhou F, Yao Y, He J, Xu D. Bioinformatics-based identification of miRNA-, lncRNA-, and mRNA-associated ceRNA networks and potential biomarkers for preeclampsia. Medicine. 2020;99(45):0000000000022985.10.1097/MD.0000000000022985Search in Google Scholar PubMed PubMed Central

[14] Medina-Bastidas D, Guzmán-Huerta M, Borboa-Olivares H, Ruiz-Cruz C, Parra-Hernández S, Flores-Pliego A, et al. Placental Microarray Profiling Reveals Common mRNA and lncRNA Expression Patterns in Preeclampsia and Intrauterine Growth Restriction. Int J Mol Sci. 2020;21(10):3597.10.3390/ijms21103597Search in Google Scholar PubMed PubMed Central

[15] Smyth GK. limma: Linear Models for Microarray Data. In: Gentleman R, Carey VJ, Huber W, Irizarry RA, Dudoit S. editors. Bioinformatics and computational biology solutions using r and bioconductor; 2005. New York, NY: Springer New York. p. 397–420.10.1007/0-387-29362-0_23Search in Google Scholar

[16] Hsu C, Juan H, Huang HC. Functional Analysis and Characterization of Differential Coexpression Networks. Sci Rep. 2015;5(1):13295.10.1038/srep13295Search in Google Scholar PubMed PubMed Central

[17] Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.10.1101/gr.1239303Search in Google Scholar PubMed PubMed Central

[18] Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic acids Res. 2016;44(W1):W90–W7.10.1093/nar/gkw377Search in Google Scholar PubMed PubMed Central

[19] Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci. 2005;102(43):15545–50.10.1073/pnas.0506580102Search in Google Scholar PubMed PubMed Central

[20] Piñero J, Bravo À, Queralt-Rosinach N, Gutiérrez-Sacristán A, Deu-Pons J, Centeno E, et al. DisGeNET: a comprehensive platform integrating information on human disease-associated genes and variants. Nucleic Acids Res. 2017;45(D1):D833–D9.10.1093/nar/gkw943Search in Google Scholar PubMed PubMed Central

[21] Yu G, Wang L, Han Y, He Q. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.10.1089/omi.2011.0118Search in Google Scholar PubMed PubMed Central

[22] Carbon S, Douglass E, Dunn N, Good BM, Harris NL, Lewis SE, et al. The Gene Ontology Resource: 20 years and still GOing strong. Nucleic Acids Res. 2019;47(D1):D330–8.10.1093/nar/gky1055Search in Google Scholar PubMed PubMed Central

[23] Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.10.1093/nar/28.1.27Search in Google Scholar PubMed PubMed Central

[24] Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523.10.1038/s41467-019-09234-6Search in Google Scholar PubMed PubMed Central

[25] Dweep H, Gretz N. miRWalk2.0: a comprehensive atlas of microRNA–target interactions. Nat Methods. 2015;12(8):697.10.1038/nmeth.3485Search in Google Scholar PubMed

[26] Chen Y, Wang X. miRDB: an online database for prediction of functional microRNA targets. Nucleic Acids Res. 2019;48(D1):D127–D31.10.1093/nar/gkz757Search in Google Scholar PubMed PubMed Central

[27] Paraskevopoulou MD, Vlachos IS, Karagkouni D, Georgakilas G, Kanellos I, Vergoulis T, et al. DIANA-LncBase v2: indexing microRNA targets on non-coding transcripts. Nucleic Acids Res. 2016;44(D1):D231–8.10.1093/nar/gkv1270Search in Google Scholar PubMed PubMed Central

[28] Gray KJ, Saxena R, Karumanchi SA. Genetic predisposition to preeclampsia is conferred by fetal DNA variants near FLT1, a gene involved in the regulation of angiogenesis. Am J Obstet Gynecol. 2018;218(2):211–8.10.1016/j.ajog.2017.11.562Search in Google Scholar PubMed PubMed Central

[29] Huang X, Zhang H, Guo X, Zhu Z, Cai H, Kong X. Insulin-like growth factor 2 mRNA-binding protein 1 (IGF2BP1) in cancer. J Hematol Oncol. 2018;11(1):88.10.1186/s13045-018-0628-ySearch in Google Scholar PubMed PubMed Central

[30] Larsson A, Palm M, Basu S, Axelsson O. Insulin-like growth factor binding protein-1 (IGFBP-1) during normal pregnancy. Gynecol Endocrinol. 2013;29(2):129–32.10.3109/09513590.2012.730574Search in Google Scholar PubMed

[31] Panjeta P, Ghalaut V, Bala J, Nanda S, Kharb S. Inverse Correlation between Insulin‑like Growth Factor‑1 and Leptin Levels in Preeclampsia. J Basic Clin Reprod Sci. 2016;5(2):94–9.10.4103/2278-960X.194480Search in Google Scholar

[32] Ingec M, Gursoy H, Yildiz L, Kumtepe Y, Kadanali S. Serum levels of insulin, IGF-1, and IGFBP-1 in preeclampsia and eclampsia. Int J Gynecol Obstet. 2004;84(3):214–9.10.1016/S0020-7292(03)00342-4Search in Google Scholar

[33] Grobman WA, Kazer RR. Serum insulin, insulin-like growth factor-I, and insulin-like growth factor binding protein-1 in women who develop preeclampsia. Obstet Gynecol. 2001;97(4):521–6.10.1097/00006250-200104000-00008Search in Google Scholar

[34] Li H, Zhang Y, Wu Q. Role of corin in the regulation of blood pressure. Curr Opnephrol Hypertens. 2017;26(2):67–73.10.1097/MNH.0000000000000297Search in Google Scholar

[35] Cui Y, Wang W, Dong N, Lou J, Srinivasan DK, Cheng W, et al. Role of corin in trophoblast invasion and uterine spiral artery remodelling in pregnancy. Nature. 2012;484(7393):246–50.10.1038/nature10897Search in Google Scholar

[36] Dong N, Zhou T, Zhang Y, Liu M, Li H, Huang X, et al. Corin mutations K317E and S472G from preeclamptic patients alert zymogen activation and cell surface targeting. J Biol Chem. 2014;289(25):17909–16.10.1074/jbc.M114.551424Search in Google Scholar

[37] Duvekot JJ, Roeters van Lennep JE Searching for new biomarkers for preeclampsia: is there a role for corin?: Mary Ann Liebert, Inc.: 140 Huguenot Street, 3rd Floor New Rochelle, NY 10801 USA; 2015.Search in Google Scholar

[38] Kimura A, Sakaguchi E, Nonaka M. Multi-component complement system of Cnidaria: C3, Bf, and MASP genes expressed in the endodermal tissues of a sea anemone, Nematostella vectensis. Immunobiology. 2009;214(3):165–78.10.1016/j.imbio.2009.01.003Search in Google Scholar

[39] Lokki AI, Kaartokallio T, Holmberg V, Onkamo P, Koskinen L, Saavalainen P, et al. Analysis of complement C3 gene reveals susceptibility to severe preeclampsia. Front Immunol. 2017;8:589.10.3389/fimmu.2017.00589Search in Google Scholar

[40] Rhim MS, Meddeb S, Kaabia O, Jalloul M, Sakouhi M, Jrzad BB, et al. C3F gene mutation is involved in the susceptibility to preeclampsia. Arch Gynecol Obstet. 2015;291(5):1023–7.10.1007/s00404-014-3515-ySearch in Google Scholar

[41] Robillard P-Y. Interest in preeclampsia for researchers in reproduction. J Reprod Immunol. 2002;53(1-2):279–87.10.1016/S0165-0378(01)00081-XSearch in Google Scholar

[42] Robillard P-Y, Hulsey TC, Dekker GA, Chaouat G. Preeclampsia and human reproduction.: An essay of a long term reflection. J Reprod Immunol. 2003;59(2):93–100.10.1016/S0165-0378(03)00040-8Search in Google Scholar

[43] Hahn S, Giaglis S, Hoesli I, Hasler P. Neutrophil NETs in reproduction: from infertility to preeclampsia and the possibility of fetal loss. Front Immunol. 2012;3:362.10.3389/fimmu.2012.00362Search in Google Scholar PubMed PubMed Central

[44] Hiby SE, Walker JJ, O'shaughnessy KM, Redman CW, Carrington M, Trowsdale J, et al. Combinations of maternal KIR and fetal HLA-C genes influence the risk of preeclampsia and reproductive success. The J Exp Med. 2004;200(8):957–65.10.1084/jem.20041214Search in Google Scholar PubMed PubMed Central

[45] Dinarello CA. Immunological and inflammatory functions of the interleukin-1 family. Annu Rev Immunol. 2009;27:519–50.10.1146/annurev.immunol.021908.132612Search in Google Scholar PubMed

[46] Gérard N, Caillaud M, Martoriati A, Goudet G, Lalmanach A-C. The interleukin-1 system and female reproduction. J Endocrinol. 2004;180(2):203–12.10.1677/joe.0.1800203Search in Google Scholar PubMed

[47] Hou Q, Liu L, Dong Y, Wu J, Du L, Dong H, et al. Genetic variations in IL1A and IL1RN are associated with the risk of preeclampsia in Chinese Han population. Sci Rep. 2014;4(1):1–6.10.1038/s41598-018-33214-3Search in Google Scholar PubMed PubMed Central

[48] Silva VRS, Soardi FC, Tanaka SCSV, da Silva-Grecco RL, Paschoini MC, Balarin MAS. Investigation of polymorphisms in preeclampsia related genes VEGF and IL1A. Arch Gynecol Obstet. 2015;291(5):1029–35.10.1007/s00404-014-3503-2Search in Google Scholar PubMed

[49] Leme Galvão LP, Menezes FE, Mendonca C, Barreto I, Alvim-Pereira C, Alvim-Pereira F, et al. Analysis of association of clinical aspects and IL1B tagSNPs with severe preeclampsia. Hypertension Pregnancy. 2016;35(1):112–22.10.3109/10641955.2015.1116554Search in Google Scholar PubMed

[50] Wu D, Xu Y, Zou Y, Zuo Q, Huang S, Wang S, et al. Long Noncoding RNA 00473 Is Involved in Preeclampsia by LSD1 Binding-Regulated TFPI2 Transcription in Trophoblast Cells. Mol Ther Nucleic Acids. 2018;12:381–92.10.1016/j.omtn.2018.05.020Search in Google Scholar PubMed PubMed Central

[51] Lin J, Ding S, Xie C, Yi R, Wu Z, Luo J, et al. MicroRNA-4476 promotes glioma progression through a miR-4476/APC/β-catenin/c-Jun positive feedback loop. Cell Death Dis. 2020;11(4):020-2474.10.1038/s41419-020-2474-4Search in Google Scholar PubMed PubMed Central

[52] Escuin D, López-Vilaró L, Bell O, Mora J, Moral A, Pérez JI, et al. MicroRNA-1291 Is Associated With Locoregional Metastases in Patients With Early-Stage Breast Cancer. Front Genet. 2020;11:562114.10.3389/fgene.2020.562114Search in Google Scholar PubMed PubMed Central

[53] Chen Y, Zhou Y, Han F, Zhao Y, Tu M, Wang Y, et al. A novel miR-1291-ERRα-CPT1C axis modulates tumor cell proliferation, metabolism and tumorigenesis. Theranostics. 2020;10(16):7193–210.10.7150/thno.44877Search in Google Scholar PubMed PubMed Central

[54] Yu S, Peng W, Zhang H, Li C, Chen X, Wei M, et al. The association between maternal and foetal REN gene polymorphisms and preeclampsia/eclampsia: A hybrid design study. Pregnancy Hypertens. 2019;18:150–5.10.1016/j.preghy.2019.09.011Search in Google Scholar PubMed

[55] Procopciuc LM, Nemeti G, Buzdugan E, Iancu M, Stamatian F, Caracostea G. Renin-angiotensin system gene variants and risk of early- and late-onset preeclampsia: A single center case-control study. Pregnancy Hypertens. 2019;18:1–8.10.1016/j.preghy.2019.08.006Search in Google Scholar PubMed

Received: 2021-06-02
Revised: 2022-01-11
Accepted: 2022-01-21
Published Online: 2022-05-18

© 2022 Guohua Li et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Biomedical Sciences
  2. Effects of direct oral anticoagulants dabigatran and rivaroxaban on the blood coagulation function in rabbits
  3. The mother of all battles: Viruses vs humans. Can humans avoid extinction in 50–100 years?
  4. Knockdown of G1P3 inhibits cell proliferation and enhances the cytotoxicity of dexamethasone in acute lymphoblastic leukemia
  5. LINC00665 regulates hepatocellular carcinoma by modulating mRNA via the m6A enzyme
  6. Association study of CLDN14 variations in patients with kidney stones
  7. Concanavalin A-induced autoimmune hepatitis model in mice: Mechanisms and future outlook
  8. Regulation of miR-30b in cancer development, apoptosis, and drug resistance
  9. Informatic analysis of the pulmonary microecology in non-cystic fibrosis bronchiectasis at three different stages
  10. Swimming attenuates tumor growth in CT-26 tumor-bearing mice and suppresses angiogenesis by mediating the HIF-1α/VEGFA pathway
  11. Characterization of intestinal microbiota and serum metabolites in patients with mild hepatic encephalopathy
  12. Functional conservation and divergence in plant-specific GRF gene family revealed by sequences and expression analysis
  13. Application of the FLP/LoxP-FRT recombination system to switch the eGFP expression in a model prokaryote
  14. Biomedical evaluation of antioxidant properties of lamb meat enriched with iodine and selenium
  15. Intravenous infusion of the exosomes derived from human umbilical cord mesenchymal stem cells enhance neurological recovery after traumatic brain injury via suppressing the NF-κB pathway
  16. Effect of dietary pattern on pregnant women with gestational diabetes mellitus and its clinical significance
  17. Potential regulatory mechanism of TNF-α/TNFR1/ANXA1 in glioma cells and its role in glioma cell proliferation
  18. Effect of the genetic mutant G71R in uridine diphosphate-glucuronosyltransferase 1A1 on the conjugation of bilirubin
  19. Quercetin inhibits cytotoxicity of PC12 cells induced by amyloid-beta 25–35 via stimulating estrogen receptor α, activating ERK1/2, and inhibiting apoptosis
  20. Nutrition intervention in the management of novel coronavirus pneumonia patients
  21. circ-CFH promotes the development of HCC by regulating cell proliferation, apoptosis, migration, invasion, and glycolysis through the miR-377-3p/RNF38 axis
  22. Bmi-1 directly upregulates glucose transporter 1 in human gastric adenocarcinoma
  23. Lacunar infarction aggravates the cognitive deficit in the elderly with white matter lesion
  24. Hydroxysafflor yellow A improved retinopathy via Nrf2/HO-1 pathway in rats
  25. Comparison of axon extension: PTFE versus PLA formed by a 3D printer
  26. Elevated IL-35 level and iTr35 subset increase the bacterial burden and lung lesions in Mycobacterium tuberculosis-infected mice
  27. A case report of CAT gene and HNF1β gene variations in a patient with early-onset diabetes
  28. Study on the mechanism of inhibiting patulin production by fengycin
  29. SOX4 promotes high-glucose-induced inflammation and angiogenesis of retinal endothelial cells by activating NF-κB signaling pathway
  30. Relationship between blood clots and COVID-19 vaccines: A literature review
  31. Analysis of genetic characteristics of 436 children with dysplasia and detailed analysis of rare karyotype
  32. Bioinformatics network analyses of growth differentiation factor 11
  33. NR4A1 inhibits the epithelial–mesenchymal transition of hepatic stellate cells: Involvement of TGF-β–Smad2/3/4–ZEB signaling
  34. Expression of Zeb1 in the differentiation of mouse embryonic stem cell
  35. Study on the genetic damage caused by cadmium sulfide quantum dots in human lymphocytes
  36. Association between single-nucleotide polymorphisms of NKX2.5 and congenital heart disease in Chinese population: A meta-analysis
  37. Assessment of the anesthetic effect of modified pentothal sodium solution on Sprague-Dawley rats
  38. Genetic susceptibility to high myopia in Han Chinese population
  39. Potential biomarkers and molecular mechanisms in preeclampsia progression
  40. Silencing circular RNA-friend leukemia virus integration 1 restrained malignancy of CC cells and oxaliplatin resistance by disturbing dyskeratosis congenita 1
  41. Endostar plus pembrolizumab combined with a platinum-based dual chemotherapy regime for advanced pulmonary large-cell neuroendocrine carcinoma as a first-line treatment: A case report
  42. The significance of PAK4 in signaling and clinicopathology: A review
  43. Sorafenib inhibits ovarian cancer cell proliferation and mobility and induces radiosensitivity by targeting the tumor cell epithelial–mesenchymal transition
  44. Characterization of rabbit polyclonal antibody against camel recombinant nanobodies
  45. Active legumain promotes invasion and migration of neuroblastoma by regulating epithelial-mesenchymal transition
  46. Effect of cell receptors in the pathogenesis of osteoarthritis: Current insights
  47. MT-12 inhibits the proliferation of bladder cells in vitro and in vivo by enhancing autophagy through mitochondrial dysfunction
  48. Study of hsa_circRNA_000121 and hsa_circRNA_004183 in papillary thyroid microcarcinoma
  49. BuyangHuanwu Decoction attenuates cerebral vasospasm caused by subarachnoid hemorrhage in rats via PI3K/AKT/eNOS axis
  50. Effects of the interaction of Notch and TLR4 pathways on inflammation and heart function in septic heart
  51. Monosodium iodoacetate-induced subchondral bone microstructure and inflammatory changes in an animal model of osteoarthritis
  52. A rare presentation of type II Abernethy malformation and nephrotic syndrome: Case report and review
  53. Rapid death due to pulmonary epithelioid haemangioendothelioma in several weeks: A case report
  54. Hepatoprotective role of peroxisome proliferator-activated receptor-α in non-cancerous hepatic tissues following transcatheter arterial embolization
  55. Correlation between peripheral blood lymphocyte subpopulations and primary systemic lupus erythematosus
  56. A novel SLC8A1-ALK fusion in lung adenocarcinoma confers sensitivity to alectinib: A case report
  57. β-Hydroxybutyrate upregulates FGF21 expression through inhibition of histone deacetylases in hepatocytes
  58. Identification of metabolic genes for the prediction of prognosis and tumor microenvironment infiltration in early-stage non-small cell lung cancer
  59. BTBD10 inhibits glioma tumorigenesis by downregulating cyclin D1 and p-Akt
  60. Mucormycosis co-infection in COVID-19 patients: An update
  61. Metagenomic next-generation sequencing in diagnosing Pneumocystis jirovecii pneumonia: A case report
  62. Long non-coding RNA HOXB-AS1 is a prognostic marker and promotes hepatocellular carcinoma cells’ proliferation and invasion
  63. Preparation and evaluation of LA-PEG-SPION, a targeted MRI contrast agent for liver cancer
  64. Proteomic analysis of the liver regulating lipid metabolism in Chaohu ducks using two-dimensional electrophoresis
  65. Nasopharyngeal tuberculosis: A case report
  66. Characterization and evaluation of anti-Salmonella enteritidis activity of indigenous probiotic lactobacilli in mice
  67. Aberrant pulmonary immune response of obese mice to periodontal infection
  68. Bacteriospermia – A formidable player in male subfertility
  69. In silico and in vivo analysis of TIPE1 expression in diffuse large B cell lymphoma
  70. Effects of KCa channels on biological behavior of trophoblasts
  71. Interleukin-17A influences the vulnerability rather than the size of established atherosclerotic plaques in apolipoprotein E-deficient mice
  72. Multiple organ failure and death caused by Staphylococcus aureus hip infection: A case report
  73. Prognostic signature related to the immune environment of oral squamous cell carcinoma
  74. Primary and metastatic squamous cell carcinoma of the thyroid gland: Two case reports
  75. Neuroprotective effects of crocin and crocin-loaded niosomes against the paraquat-induced oxidative brain damage in rats
  76. Role of MMP-2 and CD147 in kidney fibrosis
  77. Geometric basis of action potential of skeletal muscle cells and neurons
  78. Babesia microti-induced fulminant sepsis in an immunocompromised host: A case report and the case-specific literature review
  79. Role of cerebellar cortex in associative learning and memory in guinea pigs
  80. Application of metagenomic next-generation sequencing technique for diagnosing a specific case of necrotizing meningoencephalitis caused by human herpesvirus 2
  81. Case report: Quadruple primary malignant neoplasms including esophageal, ureteral, and lung in an elderly male
  82. Long non-coding RNA NEAT1 promotes angiogenesis in hepatoma carcinoma via the miR-125a-5p/VEGF pathway
  83. Osteogenic differentiation of periodontal membrane stem cells in inflammatory environments
  84. Knockdown of SHMT2 enhances the sensitivity of gastric cancer cells to radiotherapy through the Wnt/β-catenin pathway
  85. Continuous renal replacement therapy combined with double filtration plasmapheresis in the treatment of severe lupus complicated by serious bacterial infections in children: A case report
  86. Simultaneous triple primary malignancies, including bladder cancer, lymphoma, and lung cancer, in an elderly male: A case report
  87. Preclinical immunogenicity assessment of a cell-based inactivated whole-virion H5N1 influenza vaccine
  88. One case of iodine-125 therapy – A new minimally invasive treatment of intrahepatic cholangiocarcinoma
  89. S1P promotes corneal trigeminal neuron differentiation and corneal nerve repair via upregulating nerve growth factor expression in a mouse model
  90. Early cancer detection by a targeted methylation assay of circulating tumor DNA in plasma
  91. Calcifying nanoparticles initiate the calcification process of mesenchymal stem cells in vitro through the activation of the TGF-β1/Smad signaling pathway and promote the decay of echinococcosis
  92. Evaluation of prognostic markers in patients infected with SARS-CoV-2
  93. N6-Methyladenosine-related alternative splicing events play a role in bladder cancer
  94. Characterization of the structural, oxidative, and immunological features of testis tissue from Zucker diabetic fatty rats
  95. Effects of glucose and osmotic pressure on the proliferation and cell cycle of human chorionic trophoblast cells
  96. Investigation of genotype diversity of 7,804 norovirus sequences in humans and animals of China
  97. Characteristics and karyotype analysis of a patient with turner syndrome complicated with multiple-site tumors: A case report
  98. Aggravated renal fibrosis is positively associated with the activation of HMGB1-TLR2/4 signaling in STZ-induced diabetic mice
  99. Distribution characteristics of SARS-CoV-2 IgM/IgG in false-positive results detected by chemiluminescent immunoassay
  100. SRPX2 attenuated oxygen–glucose deprivation and reperfusion-induced injury in cardiomyocytes via alleviating endoplasmic reticulum stress-induced apoptosis through targeting PI3K/Akt/mTOR axis
  101. Aquaporin-8 overexpression is involved in vascular structure and function changes in placentas of gestational diabetes mellitus patients
  102. Relationship between CRP gene polymorphisms and ischemic stroke risk: A systematic review and meta-analysis
  103. Effects of growth hormone on lipid metabolism and sexual development in pubertal obese male rats
  104. Cloning and identification of the CTLA-4IgV gene and functional application of vaccine in Xinjiang sheep
  105. Antitumor activity of RUNX3: Upregulation of E-cadherin and downregulation of the epithelial–mesenchymal transition in clear-cell renal cell carcinoma
  106. PHF8 promotes osteogenic differentiation of BMSCs in old rat with osteoporosis by regulating Wnt/β-catenin pathway
  107. A review of the current state of the computer-aided diagnosis (CAD) systems for breast cancer diagnosis
  108. Bilateral dacryoadenitis in adult-onset Still’s disease: A case report
  109. A novel association between Bmi-1 protein expression and the SUVmax obtained by 18F-FDG PET/CT in patients with gastric adenocarcinoma
  110. The role of erythrocytes and erythroid progenitor cells in tumors
  111. Relationship between platelet activation markers and spontaneous abortion: A meta-analysis
  112. Abnormal methylation caused by folic acid deficiency in neural tube defects
  113. Silencing TLR4 using an ultrasound-targeted microbubble destruction-based shRNA system reduces ischemia-induced seizures in hyperglycemic rats
  114. Plant Sciences
  115. Seasonal succession of bacterial communities in cultured Caulerpa lentillifera detected by high-throughput sequencing
  116. Cloning and prokaryotic expression of WRKY48 from Caragana intermedia
  117. Novel Brassica hybrids with different resistance to Leptosphaeria maculans reveal unbalanced rDNA signal patterns
  118. Application of exogenous auxin and gibberellin regulates the bolting of lettuce (Lactuca sativa L.)
  119. Phytoremediation of pollutants from wastewater: A concise review
  120. Genome-wide identification and characterization of NBS-encoding genes in the sweet potato wild ancestor Ipomoea trifida (H.B.K.)
  121. Alleviative effects of magnetic Fe3O4 nanoparticles on the physiological toxicity of 3-nitrophenol to rice (Oryza sativa L.) seedlings
  122. Selection and functional identification of Dof genes expressed in response to nitrogen in Populus simonii × Populus nigra
  123. Study on pecan seed germination influenced by seed endocarp
  124. Identification of active compounds in Ophiopogonis Radix from different geographical origins by UPLC-Q/TOF-MS combined with GC-MS approaches
  125. The entire chloroplast genome sequence of Asparagus cochinchinensis and genetic comparison to Asparagus species
  126. Genome-wide identification of MAPK family genes and their response to abiotic stresses in tea plant (Camellia sinensis)
  127. Selection and validation of reference genes for RT-qPCR analysis of different organs at various development stages in Caragana intermedia
  128. Cloning and expression analysis of SERK1 gene in Diospyros lotus
  129. Integrated metabolomic and transcriptomic profiling revealed coping mechanisms of the edible and medicinal homologous plant Plantago asiatica L. cadmium resistance
  130. A missense variant in NCF1 is associated with susceptibility to unexplained recurrent spontaneous abortion
  131. Assessment of drought tolerance indices in faba bean genotypes under different irrigation regimes
  132. The entire chloroplast genome sequence of Asparagus setaceus (Kunth) Jessop: Genome structure, gene composition, and phylogenetic analysis in Asparagaceae
  133. Food Science
  134. Dietary food additive monosodium glutamate with or without high-lipid diet induces spleen anomaly: A mechanistic approach on rat model
  135. Binge eating disorder during COVID-19
  136. Potential of honey against the onset of autoimmune diabetes and its associated nephropathy, pancreatitis, and retinopathy in type 1 diabetic animal model
  137. FTO gene expression in diet-induced obesity is downregulated by Solanum fruit supplementation
  138. Physical activity enhances fecal lactobacilli in rats chronically drinking sweetened cola beverage
  139. Supercritical CO2 extraction, chemical composition, and antioxidant effects of Coreopsis tinctoria Nutt. oleoresin
  140. Functional constituents of plant-based foods boost immunity against acute and chronic disorders
  141. Effect of selenium and methods of protein extraction on the proteomic profile of Saccharomyces yeast
  142. Microbial diversity of milk ghee in southern Gansu and its effect on the formation of ghee flavor compounds
  143. Ecology and Environmental Sciences
  144. Effects of heavy metals on bacterial community surrounding Bijiashan mining area located in northwest China
  145. Microorganism community composition analysis coupling with 15N tracer experiments reveals the nitrification rate and N2O emissions in low pH soils in Southern China
  146. Genetic diversity and population structure of Cinnamomum balansae Lecomte inferred by microsatellites
  147. Preliminary screening of microplastic contamination in different marine fish species of Taif market, Saudi Arabia
  148. Plant volatile organic compounds attractive to Lygus pratensis
  149. Effects of organic materials on soil bacterial community structure in long-term continuous cropping of tomato in greenhouse
  150. Effects of soil treated fungicide fluopimomide on tomato (Solanum lycopersicum L.) disease control and plant growth
  151. Prevalence of Yersinia pestis among rodents captured in a semi-arid tropical ecosystem of south-western Zimbabwe
  152. Effects of irrigation and nitrogen fertilization on mitigating salt-induced Na+ toxicity and sustaining sea rice growth
  153. Bioengineering and Biotechnology
  154. Poly-l-lysine-caused cell adhesion induces pyroptosis in THP-1 monocytes
  155. Development of alkaline phosphatase-scFv and its use for one-step enzyme-linked immunosorbent assay for His-tagged protein detection
  156. Development and validation of a predictive model for immune-related genes in patients with tongue squamous cell carcinoma
  157. Agriculture
  158. Effects of chemical-based fertilizer replacement with biochar-based fertilizer on albic soil nutrient content and maize yield
  159. Genome-wide identification and expression analysis of CPP-like gene family in Triticum aestivum L. under different hormone and stress conditions
  160. Agronomic and economic performance of mung bean (Vigna radiata L.) varieties in response to rates of blended NPS fertilizer in Kindo Koysha district, Southern Ethiopia
  161. Influence of furrow irrigation regime on the yield and water consumption indicators of winter wheat based on a multi-level fuzzy comprehensive evaluation
  162. Discovery of exercise-related genes and pathway analysis based on comparative genomes of Mongolian originated Abaga and Wushen horse
  163. Lessons from integrated seasonal forecast-crop modelling in Africa: A systematic review
  164. Evolution trend of soil fertility in tobacco-planting area of Chenzhou, Hunan Province, China
  165. Animal Sciences
  166. Morphological and molecular characterization of Tatera indica Hardwicke 1807 (Rodentia: Muridae) from Pothwar, Pakistan
  167. Research on meat quality of Qianhua Mutton Merino sheep and Small-tail Han sheep
  168. SI: A Scientific Memoir
  169. Suggestions on leading an academic research laboratory group
  170. My scientific genealogy and the Toronto ACDC Laboratory, 1988–2022
  171. Erratum
  172. Erratum to “Changes of immune cells in patients with hepatocellular carcinoma treated by radiofrequency ablation and hepatectomy, a pilot study”
  173. Erratum to “A two-microRNA signature predicts the progression of male thyroid cancer”
  174. Retraction
  175. Retraction of “Lidocaine has antitumor effect on hepatocellular carcinoma via the circ_DYNC1H1/miR-520a-3p/USP14 axis”
Downloaded on 5.8.2025 from https://www.degruyterbrill.com/document/doi/10.1515/biol-2022-0053/html
Scroll to top button