Open Access

Preeclampsia: a bioinformatics approach through protein-protein interaction networks analysis

BMC Systems Biology20126:97

DOI: 10.1186/1752-0509-6-97

Received: 5 January 2012

Accepted: 23 July 2012

Published: 8 August 2012

Abstract

Background

In this study we explored preeclampsia through a bioinformatics approach. We create a comprehensive genes/proteins dataset by the analysis of both public proteomic data and text mining of public scientific literature. From this dataset the associated protein-protein interaction network has been obtained. Several indexes of centrality have been explored for hubs detection as well as the enrichment statistical analysis of metabolic pathway and disease.

Results

We confirmed the well known relationship between preeclampsia and cardiovascular diseases but also identified statistically significant relationships with respect to cancer and aging. Moreover, significant metabolic pathways such as apoptosis, cancer and cytokine-cytokine receptor interaction have also been identified by enrichment analysis. We obtained FLT1, VEGFA, FN1, F2 and PGF genes with the highest scores by hubs analysis; however, we also found other genes as PDIA3, LYN, SH2B2 and NDRG1 with high scores.

Conclusions

The applied methodology not only led to the identification of well known genes related to preeclampsia but also to propose new candidates poorly explored or completely unknown in the pathogenesis of preeclampsia, which eventually need to be validated experimentally. Moreover, new possible connections were detected between preeclampsia and other diseases that could open new areas of research. More must be done in this area to resolve the identification of unknown interactions of proteins/genes and also for a better integration of metabolic pathways and diseases.

Background

Preeclampsia is a pregnancy related disease associated with hypertension, proteinuria and increased maternal and perinatal morbidity and mortality, without known underlying mechanism and preventive treatment[1, 2]. On the other hand, the future health or possible risks of women with previous history of preeclampsia are important areas of investigation. In this direction, it is well known the increased risk of future cardiovascular disease and renal dysfunction, however, other risks are also being discussed[1, 35]. Owing to the morbidity and mortality of this pregnancy related disease and the possible multifactorial causes involved[15], several experimental procedures have been applied by researchers in the last two decades, evidently, generating an elevated number of unprocessed information.

Although some bioinformatic analysis has been performed in particular microarray assays[6, 7], an extensive data evaluation and processing has not yet been performed. Furthermore, the capabilities of bioinformatics tools for gene prioritization, network analysis, gene ontology and gene-disease relationships[8, 9], together with all available data on protein/gene expression during preeclampsia bring an interesting and valuable opportunity for an in-silico study of the disease. Therefore, the present study is focused on two main areas: I) collection and basic analysis of the genes/proteins-diseases dataset, including, protein-protein interaction network and pathway enrichment analysis and II) exploration of the related gene-diseases in order to evaluate other genetic diseases possibly related with preeclampsia.

Results

Protein-protein interaction network analysis

Preeclampsia PPI network topology reveals (Figure1) a similar behavior with respect to general topology of PPI following a power law behavior[10] and therefore scale-free properties. These types of networks have the particular feature that some nodes are highly connected compared with others on the same network. These highly connected nodes (hubs) in general, represent important proteins/genes in biological terms and therefore are treated with special attention.
https://static-content.springer.com/image/art%3A10.1186%2F1752-0509-6-97/MediaObjects/12918_2012_Article_950_Fig1_HTML.jpg
Figure 1

PPI network and topology. Left) PPI network and Right) Degree distribution. The degree distribution follows a power law distribution.

The top 50 genes with high scores and also present in the initial set (347) are shown in Table1, however, other genes were found with high scores value but there are not part of the initial gene group. As expected some of the selected genes like FN1, FLT1, F2, VEGFA, PGF, TNF, NOS and INHBA, are well known preeclampsia relates genes (see discussion) and several of them are related with signaling pathways.
Table 1

Top 50 genes obtained by analysis of the PPI network

Gene symbol

Entrez1

Score I

Gene symbol

Entrez1

Score I

4067

LYN

983.5136

2324

FLT4

426.74

2335

FN1

823.9037

6767

ST13

425.5747

29110

TBK1

672.1386

873

CBR1

424.257

8826

IQGAP1

635.9629

667

DST

414.0964

1759

DNM1

628.7167

9495

AKAP5

412.7567

5573

PRKAR1A

622.5864

8601

RGS20

408.6292

10397

NDRG1

590.672

2896

GRN

401.0426

2923

PDIA3

559.3358

4323

MMP14

400.1593

5339

PLEC

541.5909

55697

VAC14

394.4482

6303

SAT1

515.6939

5228

PGF

393.3167

2321

FLT1

508.3947

8553

BHLHE40

392.8039

2317

FLNB

505.1671

604

BCL6

385.7928

4221

MEN1

499.019

5054

SERPINE1

382.5547

7067

THRA

487.6884

4597

MVD

380.0219

2147

F2

482.7173

5569

PKIA

377.6185

3959

LGALS3BP

475.7565

140885

SIRPA

375.9565

10603

SH2B2

473.9241

3624

INHBA

367.5795

7422

VEGFA

465.5212

10272

FSTL3

367.2976

8815

BANF1

461.918

4846

NOS3

361.5554

23650

TRIM29

448.6764

3553

IL1B

357.2569

7124

TNF

441.7007

3880

KRT19

355.4412

904

CCNT1

440.8304

9444

QKI

354.3211

3556

IL1RAP

430.4483

8773

SNAP23

348.0265

9131

AIFM1

429.2521

27250

PDCD4

345.9363

6449

SGTA

429.1609

5291

PIK3CB

345.5114

The table shown the first 50 genes selected according to the score values and arranged in descendent order. We also provide the gene symbol the associated Entrez Gene identifier.

Notes: 1) Entrez gene ID.

A total of 27 communities (k = 3) covering 161 genes were identified by communality analysis. In Figure2 (Left) we represent those communities that are superimposed to form a large connected graph. The genes involved in communities overlapping are also highly represented in the Table1 (and also the genes members of the large community). The model based clustering analysis reveal an optimal number of 8 clusters (BIC = 152192.4) with an ellipsoidal distribution with equal volume, shape and variable orientation. The genes are grouped in the clusters (C1…8) as follow: C1 (67), C2 (56), C3 (1806), C4 (59), C5 (133), C6 (23), C7 (95) and C8 (161). The C8 and C4 correspond with the highest mean scoring value: 393.3 and 348.9 respectively, and contain all the 100 genes with highest score values (part or not of the initial gene set). Furthermore, 161 genes of C8 are also the same genes detected in the communality analysis.
https://static-content.springer.com/image/art%3A10.1186%2F1752-0509-6-97/MediaObjects/12918_2012_Article_950_Fig2_HTML.jpg
Figure 2

Communality and clusters analysis. Left) Representation of the largest connected community. Red nodes represent the genes involved in communities overlapping. White nodes represent the bigger community. Right) Representation of C8 and C4 clusters and the related statistically significant biological process obtained by gene ontology enrichment analysis.

Gene ontology (GO) enrichment analyses were performed in all obtained clusters. However, for simplicity only C4 and C8 are presented (Figure2 Right). The GO analysis reveals that C8 comprise several processes related with angiogenesis, apoptosis and cell proliferation and also shared with C4 several processess involved in cell activation and biological adhesion. The relation between these processes as well as the fact that both groups are representative of the highest scored genes could indicate a particular relevance of the clusters in terms of genes-disease relationship. On the other hand, also these processes are well known involved in preeclampsia and are also consistent with the pathway enrichment analysis.

Diseases and metabolic pathway enrichment analysis

Several types of diseases were found statistically significant in the enrichment analysis; partial results are presented in Table2. Obviously, preeclampsia and even hypertension have to be present in the analysis. In the GAD database there are several disease classes and beside the presented in Table2, others like hematologic (p-value = 6.1E-06) and renal diseases (p-value = 2.5E-04) were also significant. Even when we present only the results obtained with GAD database, analysis was also performed with OMIN database confirming the ovarian cancer (p-value = 0.019) and also indicating colorectal cancer as statistically significant (p-value = 0.011), however, better results were obtained with GAD and also more consistent with literature information and pathway enrichment analysis.
Table 2

The diseases enrichment analysis

Disease classes

Disease

p-value

Cancer (p-value = 1.24E-11)

breast cancer

9.76E-07

ovarian cancer

4.30E-04

lupus erythematosus

4.4E-02

Cardiovascular (p-value = 4.28E-9)

myocardial infarction

8.7E-09

coronary artery disease

7.38E-06

Stroke

6.16E-05

heart disease, ischemic

9.7E-03

blood pressure, arterial

8.2E-03

inflammatory bowel disease

2.9E-02

Aging (p-value = 5.23E-05)

Alzheimer's Disease

2.36E-07

Longevity

2.17E-04

atherosclerosis

3.7E-02

Arthritis

3.60E-02

Reproduction (p-value = 7.9E-03)

preeclampsia

1.29E-05

pregnancy loss, recurrent

6.08E-04

The table provides the p-value obtained from disease enrichment analysis according to the Disease Class and also the sub-classes diseases following the GAD database structure. The results showed are partials as discussed in the respective section.

It is important to consider that several genes in the PPI network do not present a known relation with specific diseases, at least reported in the GAD or OMIN databases. Only around 30% of the 2400 genes were found in the databases. This difficulty means that we have to be cautious with the preeclampsia genes-diseases relationships and with the reliability of the statistical p-value, even when some important and significant inferences, can be made.

A similar situation occurred with the pathway enrichment analysis (Table3). Even when the KEGG database is the most representative of our gene space and high coherence was noticed with the physiopathology of PE, the results only cover around 50% of the initial 2400 genes. A similar procedure was also performed with the Reactome and BioCarta databases with a less covering (37% and 27% respectively) and showing a high coherence with Table3 results. These databases reveal other significant pathways like NGF, PDGF, BMP, EPO and EGFR signaling as well as apoptosis and hemostasis pathways (data not shown). Some cancerous pathways (i.e. prostatic, pancreatic and lung) were also found statistically significant in KEGG but were excluded from Table3, in order to simplify, because many of them have similar reactions with the general cancer pathway, already presented in the Table3.
Table 3

The KEGG pathway enrichment analysis

Pathway description

p-value

Pathway description

p-value

Pathways in cancer

1.10E-42

Fc gamma R-mediated phagocytosis

2.90E-07

Focal adhesion

3.50E-29

RIG-I-like receptor signaling pathway

5.30E-07

Apoptosis

2.10E-17

Chemokine signaling pathway

5.90E-07

Neurotrophin signaling pathway

3.40E-17

Leukocyte transendothelial migration

6.20E-07

TGF-beta signaling pathway

5.10E-16

VEGF signaling pathway

8.60E-07

T cell receptor signaling pathway

8.90E-15

NOD-like receptor signaling pathway

1.20E-06

ErbB signaling pathway

1.80E-14

mTOR signaling pathway

7.00E-06

B cell receptor signaling pathway

8.70E-14

Complement and coagulation cascades

7.90E-06

Adherens junction

2.00E-12

Tight junction

1.10E-05

Fc epsilon RI signaling pathway

3.60E-12

RNA polymerase

1.60E-05

ECM-receptor interaction

6.80E-12

Ubiquitin mediated proteolysis

2.20E-05

Insulin signaling pathway

3.00E-10

GnRH signaling pathway

1.10E-04

Regulation of actin cytoskeleton

3.10E-09

Hematopoietic cell lineage

1.20E-04

Toll-like receptor signaling pathway

1.01E-08

Cytokine-cytokine receptor interaction

2.00E-04

Wnt signaling pathway

1.10E-08

Hedgehog signaling pathway

3.70E-04

Cell cycle

1.50E-08

Natural killer cell mediated cytotoxicity

3.80E-04

MAPK signaling pathway

1.50E-08

Cytosolic DNA-sensing pathway

6.70E-04

Adipocytokine signaling pathway

7.50E-08

Notch signaling pathway

2.60E-03

Progesterone-mediated oocyte maturation

2.50E-07

Renin-angiotensin system

1.80E-02

The table provides the p-value obtained from pathway enrichment analysis in the KEGG database. The pathways names correspond with the KEGG nomenclature. The results showed are partials as discussed in the respective section.

In order to simplify and enhance the understanding of the involved pathways and their relationship with the selected hubs, a fusion between both was made (Figure3). However, it is important to exalt that from the 50 hubs previously selected; only 22 present some significant pathway association with Table2. The genes: NDRG1, LGALS3BP, BANF1, SGTA, TRIM29, RGS20, PLEC, GRN, ST13, AKAP5, FSTL3, DST, PKIA, QKI, MLF2 and KRT19, for example, were not found in the KEGG database.
https://static-content.springer.com/image/art%3A10.1186%2F1752-0509-6-97/MediaObjects/12918_2012_Article_950_Fig3_HTML.jpg
Figure 3

Genes-metabolic pathways interaction. The genes and pathways were selected after hubs detection (see Table1) and enrichment analysis respectively (see Table3).

Discussion

PPI analysis

The top 50 genes selected were manually analyzed with the scientific literature in order to validate its connection with PE or even changes during pregnancy and we corroborate that several of them like FLT1, FLT4, VEGFA, PGF, TNF, FN1, F2 and NOS3 can be related to the main lines of research in preeclampsia pathogenesis hypothesis and specially angiogenesis[1116]. Moreover, PGF, INHBA and related inhibin as well as activin proteins have been considered in several predictive model of PE[17, 18]. The presence of those genes in our selection could validate the method applied and increase our confidence with respect to those genes that are poorly explored or unknown in their association with preeclampsia. In the latter group, we have for example the genes LYN, PDIA3, NDRG1 and TBK1. The PDIA3 genes encode an endoplasmic protein highly related with folding processes and in its relation with pregnancy only one article was found[19] revealing that PDIA3 gene intervene in trophoblast invasion via interleukin (IL) 11. Similarly, TBK1 could also intervene during the first moments of pregnancy to secure the implantation in relation with the nuclear factor kappa B[20]. However, there are not articles of PDIA3, TBK1 or LYN related with preeclampsia. Moreover, only one article was found describing an increased expression of NDRG1 during preeclampsia[21]. A similar problem is found with other genes such as IQGAP1, DNM1, SAT1, MEN1 and SH2B2 that also have been little explored during pregnancy.

Cluster analysis indicates that mainly C8 but also C4 probably embrace the most significant genes, at least related with centrality. All genes highest scored are part of these clusters as well as the large community graph and as we can notice the genes that lead to a community superposition are also highly scored in Table1 (Figure2). On the other hand there are other genes like ENG, VEGFB and INHA that are well known related to preeclampsia and are also part of the C8 cluster[1118]. It is important to notice that there are other genes not shown in Table1 because were not present in the initial gene set. In this group, EGFR and GRB2, are both with highest scores and there are both well related with preeclampsia[22, 23]. Thus, even when the work presented focused on the analysis of the initial group, it is possible that relevant genes were identified by PPI topology but not included in the initial subset. Moreover, the GO enrichment analysis clearly reveals that C8 and C4 clusters are related with angiogenesis, apoptosis and biological adhesion, which are also biological processes obtained with the pathways enrichment analysis. Angiogenesis and related processes (i.e. vascular growing, cell differentiation and hypoxia) are considered as central processes related to preeclampsia[24] and therefore it could support the reliability of the procedure and also the necessity to increase the study of the C8 and C4 genes.

Future combination of centrality indexes and specific clusters selection together with machine learning procedures or genetic algorithm optimization based on groups differentiability or external prediction subset, could reduce even more the final gene space and we are currently working in this direction.

Diseases and pathway analysis

We manually evaluate through scientific literature the connection of PE with the identified metabolic pathways and diseases. The relationship between preeclampsia and inflammatory, immunologic, angiogenic and hemodynamic responses are very well documented[2228] and therefore are expected not only in the gene space but also in the metabolic pathways analysis. The metabolic pathway analysis (Table3) indicates a strong significance of the cancer pathway that is consistent with the disease analysis (Table2) and also with the VEGF signaling, that is present in several related processes like cytokine-cytokine receptor, angiogenesis and also cancer pathway. With the procedure carried out in the present study a simplification of the metabolic pathways was possible, however, more needs to be done in this area in order to better integrate not only the hubs genes but also the comprehensive data created by the PPI interaction.

Considering the associated metabolic pathways we can notice that several signaling pathways are statistically significant and some of them are connected with PIK3CB, LYN and linked in several cases with TNF (Figure3) that is a central gene affecting several processes and also widely studied in its relation with preeclampsia[28]. Similarly, SH2B2 that encodes an adapter protein of several tyrosine kinase receptors is also connected with metabolic pathways indirectly related with apoptosis. Contrasting PDIA3 and NDRG1 are not present in the Figure3 connected with any metabolic pathway, but the recent relation found between PDIA3 and IL11 could open a relationship with cytokine-cytokine receptors[19], specifically through the hematopoietins receptor family. We can also notice (Figure3 and Table3) that in the pathways associated with central proteins (the hubs), those highly connected are the cytokine-cytokine receptors, focal adhesion and apoptosis pathways and they contain almost the complete space of genes mainly studied in preeclampsia.

The elevated number of signaling pathways that we found statistically significant, together with the hubs distribution detected in Table1, seems to point out the idea of a signaling related disease similar to cancer, however, the apoptosis and angiogenic mechanism had been related previously with PE and are also highly represented in our study.

On the other hand, the relationship between PE and cardiovascular diseases is well known. Women with previous history of PE or even pregnancy hypertension present an increased risk of future development of cardiovascular disease or hypertension[2931]. This is clearly expressed in our results (Table2). Moreover, our results also suggest a significant relationship between cancer and PE (at this point it is not possible to say if this correlation is positive or negative) indicated by diseases enrichment analysis and also by the metabolic pathways. However, the experimental and epidemiological evidence that could support the influence of PE in future cancer development is for now inconclusive, even when could be reasonably expected by the wide kind of genes that are actually shared between both diseases (i.e. angiogenesis related genes). Several articles report a reduced risk of future cancer development in women with previous history of PE, but others could not find any statistical relationship[3235].

There are neither clinical signs in common between PE and Alzheimer and nor epidemiological studies. However, the connection between both diseases has been questioned before (at least in late-onset Alzheimer's disease)[36] and was significant in the present study. Therefore, although it is not possible to validate our findings experimentally with the current information, it actually opens new possibilities for future works. A similar problem concerns to other ageing related factors like atherosclerosis, arthritis or longevity (that for obvious reasons will be difficult to explore in PE during long-term studies).

The present study points out several advantages of the bioinformatics analysis but some limitations were also found. The detection of genes that are very well known to be involved in PE thought the applied methodology as well as the consistence of the results across metabolic information could support the novel candidates found, however, it is necessary to reduce even more the genes space applying other methodologies as well as to design new experimental experiences. On the other hand, the limitation of the human protein interaction information suggests that also orthologous genes should be needed in order to increase the PPI covering of the initial dataset and to increase the capabilities of the metabolic pathways and disease enrichment analysis.

Conclusions

In the present study we applied several bioinformatics tools in order to create a comprehensive initial database of genes statistically related to preeclampsia and a further expansion through the construction of related protein-protein interaction network. Using several centrality indexes for hubs detection, some well know preeclampsia related genes like FLT1, TNF, VEGFA and PGF were detected as well as other genes with high scores, like PDIA3, NDRG1, TBK1, LYN, IQGAP1, DNM1, SAT1, MEN1 and SH2B2 that have been poorly explored or unknown in the current state of the art of preeclampsia physiopathology.

Through disease enrichment analysis we corroborated the well know connection between PE and cardiovascular disease, but we also found a possible link of PE with aging and cancer. Moreover we also found the cancer pathway, focal adhesion, ECM-receptor interaction, apoptosis and cytokine-cytokine receptors interactions metabolic pathways highly represented by the PPI network that are in agreement with some of the preeclampsia-diseases relationship found, as well as with the central topics of the current preeclampsia investigations.

Methods

Dataset construction

Experimental proteomic/genomic data comparing normal (N) and preeclamptic (PE) pregnancies was obtained analysing the Gene Expression Omnibus (GEO) database[37]. The datasets considered are represented in Table4.
Table 4

GEO data used for dataset construction

GEO

Sample

Method

Collection (Weeks)

Tissue

Ref.

GSE6573

10(PE), 10(N)

Affymetrix

At delivery

Placenta

[25]

GSE47071

10(PE), 4(N)

Agilent

See notes

Placenta

[6]

GSE25906

23(PE), 37(N)

Illumina

27-38

Placenta

[26]

GSE147222

12(PE), 11(N)

Affymetrix

See notes

Placenta

[27]

GSE10588

17(PE), 26(N)

ABI Human

At delivery

Placenta

[38]

The Table shown the GEO data sources, references and complementary information used in the study.

Notes: 1) The 10 samples were collected for placenta biopsy during delivery. Five of the preeclamptic samples were collected before 31 weeks of gestation and the rest after this time. 2) The 23 samples were collected for placenta biopsy during delivery, however, because in GEO appear two platform types, both were analysed.

Each experiment was analyzed independently in order to reduce the number of genes. In our case we considered an adjusted p-value ≤0.05 and a fold expression ≥2 as discussed elsewhere[6, 7, 2527, 38, 39]. Initially the p-value was obtained by a bootstrapping procedure with 1000 or 10000 iterations (depending on size of the sample) obtaining 645 statistically significant modulated genes, however, applying the false discovery rate (FDR) correction by the Benajmini-Hochberg method[40], this sample was reduced to 330 genes.

In addition, several text mining exploration tools were used to complement the GEO results. There are several tools to perform a text data mining analysis but several of them require extra information (i.e. chromosome region) instead phenotype or diseases notation (i.e. diseases name or related keywords). In our case we choose those methods that do not require previous genetic knowledge of the disease[8]. Moreover, the text mining procedures usually could provide several false positive associations and therefore those tools which also combine text-mining with other data sources in the analysis are preferred[8, 41]. Considering these aspects, we used the following tools: PolySearch[42], Candid[43] and Phenopred[44]. Candid and PhenoPred use several heterogeneous data sources to overcome bias while in PolySearch analyse was restricted to PubMed publications. Obviously many other algorithms could be used in alternative. In order to reduce the risk of including biased relationships, the top 10–20 genes/proteins with the highest scores were selected and individually analyzed considering the preeclampsia related scientific publications. Some of the top genes were also present in the previous dataset (GEO), therefore, the final dataset contained 347 genes.

Protein-protein interaction network (PPI)

The proteins associated with the previous 347 genes were identified and cross-referenced with the IRefIndex (v1.16)[45] and a signaling curated databases[46] that were used to create the protein-protein interaction (PPI) network. The IRefIndex database provide an index of protein interactions available in several databases like: BIND, BioGRID, DIP, HPRD that simplify the task consuming process of inter-database mapping and lead to a comprehensive covering of the available known protein interactions space. On the other hand, this PPI database is easily integrated in Cytoscape. Additionally, many diseases are related with signaling pathways modifications and therefore the inclusion of this interaction database considerable improve the PPI space. The interaction search was restricted to Homo Sapiens and includes all kind of experimental procedure as well as some predictive interactions (mainly from the OPHID database). The curation of the final database was performed both, manually and using home-made software to remove duplicate interactions and unify isoforms notation with unique genes. We obtained our final PPI network with 3279 interactions and 2400 nodes.

Some of the proteins present in our initial dataset had not any known experimental interaction (at least in humans) and therefore the 2400 nodes cover only 234 (67.45%) genes of the initial set (347). The network visualization and network topology indexes, calculated in the hubs detection process, were carried out using Cytoscape 2.8.2 and CytoHubba[47, 48].

Several methodologies are available for hubs and essential genes identification, and all of them with the respective advantage and limitations[47, 4954]. Some strategies are the use of genetic algorithm or machine learning procedures[49, 50], however, the centrality approaches are by far the most applied procedures even by simplicity and because several studies had being pointed out its applicability[47, 51, 52]. Therefore, several centrality indexes were evaluated: Betweenness, bottleneck, density of maximum neighborhood component (DMNC), node degree, edge percolation component (EPC), eccentricity, maximal clique centrality (MCC), maximum neighborhood component (MNC), radiality and stress[47]. On the other hand, to obtain a scoring index we created the measurement (Score I) as follow:
Score I = i = 1 N c 1 max I c i I c i max I c i min I c i · 100
(1)

Where Ici is the values of centrality indexes and i = 1…Nc, and is the number of calculated centrality indexes (Nc = 10). As we can note, score I is the sum of all the indexes percent value after individual normalization and therefore is restricted to a maximal value of 100×Nc which simplify even better the top genes selection. With the normalized centrality indexes we also performed a model-based clustering analysis using R-package[53] in order to study hubs distribution with respect to centrality ranks. We also performed a communality (or cliques) network analysis by clique percolation method using CFinder[54]. The communality analysis provides a better topology description of the network including the location of highly connected sub-graphs (cliques) and/or overlapping modules that usually correspond with relevant biological information.

Pathway and diseases enrichment analysis

The pathways and diseases enrichment analysis were performed through the DAVID bioinformatics resource 6.7[55], exploring the well know databases: KEGG, BioCarta and Reactome (pathways related) as well as OMIN and Genetic Association Database (GAD) (diseases related analysis). This online resource (DAVID) integrate, in a faster computational analysis, a wide range of enrichment analysis thought different databases providing also a substantial statistical description. The analysis was carried out considering the complete gene space of the PPI network. We also used DAVID in order to perform a gene ontology enrichment analysis in the obtained clusters.

Declarations

Acknowledgements

This study was supported by ``Fundação para a Ciência e a Tecnologia'' (FCT), grant: SFRH/BPD/72391/2010.

Authors’ Affiliations

(1)
Department of Biological Sciences, Biochemistry, Faculty of Pharmacy, University of Porto, Portugal/Institute for Molecular and Cell Biology (IBMC)
(2)
Department of Obstetrics and Gynecology, Faculty of Medicine, University of Porto, Department of Obstetrics and Gynecology, Sao Joao Hospital./INEB – Institute of Biomedical Engineering, University of Porto

References

  1. Leeman L, Fontaine P: Hypertensive disorders of pregnancy. Am Fam Physician. 2008, 78 (1): 93-100.
  2. National Collaborating Centre for Women's and Children's Health: Hypertension in pregnancy. The management of hypertensive disorders during pregnancy. 2010, London (UK): National Institute for Health and Clinical Excellence (NICE), 46-Clinical guideline; no. 107
  3. Sheppard SJ, Khalil RA: Risk factors and mediators of the vascular dysfunction associated with hypertension in pregnancy. Cardiovasc Hematol Disord Drug Targets. 2010, 10 (1): 33-52.View Article
  4. Borna S, Neamatipoor E, Radman N: Risk of coronary artery disease in women with history of pregnancies complicated by Preeclampsia and LBW. J Matern Fetal Neonatal Med. 2011, 19-Epub ahead of print
  5. Carty DM, Delles C, Dominiczak AF: Preeclampsia and future maternal health. J Hypertens. 2010, 28 (7): 1349-1355. 10.1097/HJH.0b013e32833a39d0.View Article
  6. Nishizawa H, Pryor-Koishi K, Kato T, Kowa H, Kurahashi H, Udagawa Y: Microarray analysis of differentially expressed fetal genes in placental tissue derived from early and late onset severe pre-eclampsia. Placenta. 2007, 28 (5–6): 487-497.View Article
  7. Founds SA, Conley YP, Lyons-Weiler JF, Jeyabalan A, Hogge WA, Conrad KP: Altered global gene expression in first trimester placentas of women destined to develop Preeclampsia. Placenta. 2009, 30 (1): 15-24. 10.1016/j.placenta.2008.09.015.View Article
  8. Tranchevent LC, Capdevila FB, Nitsch D, Moor BD, Causmaecker PD, Moreau Y: A guide to web tools to prioritize candidate genes. Brief Bioinform. 2010, 12 (1): 22-32.View Article
  9. Özgür A, Vu T, Erkan G, Radev DR: Identifying gene-disease associations using centrality on a literature mined gene-interaction network. Bioinformatics. 2008, 24: i277-i285. 10.1093/bioinformatics/btn182.View Article
  10. Stelzl U, Worm U, Lalowski M, Haenig C, Brembeck FH, Goehler H, Stroedicke M, Zenkner M, Schoenherr A, Koeppen S, Timm J, Mintzlaff S, Abraham C, Bock N, Kietzmann S, Goedde A, Toksoz E, Droege A, Krobitsch S, Korn B, Birchmeier W, Lehrach H, Wanker EE: A human protein-protein interaction network: a resource for annotating the proteome. Cell. 2005, 122 (6): 957-968. 10.1016/j.cell.2005.08.029.View Article
  11. Ward K: Searching for genetic factors underlying pre-eclampsia: recent progress and persistent challenges. Minerva Ginecol. 2008, 60 (5): 399-419.
  12. Nejatizadeh A, Stobdan T, Malhotra N, Pasha MA: The genetic aspects of pre-eclampsia: achievements and limitations. Biochem Genet. 2008, 46 (7–8): 451-479.View Article
  13. Laresgoiti-Servitje E, Gómez-López N, Olson DM: An immunological insight into the origins of pre-eclampsia. Hum Reprod Update. 2010, 16 (5): 510-524. 10.1093/humupd/dmq007.View Article
  14. Winn VD, Gormley M, Fisher SJ: The impact of Preeclampsia on gene expression at the maternal-fetal interface. Pregnancy Hypertens. 2011, 1 (1): 100-108. 1
  15. Ishida J, Matsuoka T, Saito-Fujita T, Inaba S, Kunita S, Sugiyama F, Yagami K, Fukamizu A: Pregnancy-associated homeostasis and dysregulation: lessons from genetically modified animal models. J Biochem. 2011, 150 (1): 5-14. 10.1093/jb/mvr069.View Article
  16. Rajakumar A, Chu T, Handley DE, Bunce KD, Burke B, Hubel CA, Jeyabalan A, Peters DG: Maternal gene expression profiling during pregnancy and preeclampsia in human peripheral blood mononuclear cells. Placenta. 2011, 32 (1): 70-78. 10.1016/j.placenta.2010.10.004.View Article
  17. Kuc S, Wortelboer EJ, van Rijn BB, Franx A, Visser GH, Schielen PC: Evaluation of 7 serum biomarkers and uterine artery Doppler ultrasound for first-trimester prediction of preeclampsia: a systematic review. Obstet Gynecol Surv. 2011, 66 (4): 225-239. 10.1097/OGX.0b013e3182227027.View Article
  18. Spencer K, Yu CK, Savvidou M, Papageorghiou AT, Nicolaides KH: Prediction of pre-eclampsia by uterine artery Doppler ultrasonography and maternal serum pregnancy-associated plasma protein-A, free beta-human chorionic gonadotropin, activin A and inhibin A at 22 + 0 to 24 + 6 weeks' gestation. Ultrasound Obstet Gynecol. 2006, 27 (6): 658-663. 10.1002/uog.2676.View Article
  19. Sonderegger S, Yap J, Menkhorst E, Weston G, Stanton PG, Dimitriadis E: Interleukin (IL)11 mediates protein secretion and modification in human extravillous trophoblasts. Hum Reprod. 2011, 26 (10): 2841-2849. 10.1093/humrep/der259.View Article
  20. King AE, Critchley HO, Kelly RW: The NF-kappaB pathway in human endometrium and first trimester decidua. Mol Hum Reprod. 2001, 7 (2): 175-183. 10.1093/molehr/7.2.175.View Article
  21. Choi SJ, Oh SY, Kim JH, Sadovsky Y, Roh CR: Increased expression of N-myc downstream-regulated gene 1 (NDRG1) in placentas from pregnancies complicated by intrauterine growth restriction or preeclampsia. Am J Obstet Gynecol. 2007, 196 (1): 45.e1-7. 10.1016/j.ajog.2006.08.029.View Article
  22. Anteby EY, Ayesh S, Shochina M, Hamani Y, Schneider T, Al-Shareef W, Hochberg A, Ariel I: Growth factor receptor-protein bound 2 (GRB2) upregulation in the placenta in preeclampsia implies a possible role for ras signalling. Eur J Obstet Gynecol Reprod Biol. 2005, 118 (2): 174-181. 10.1016/j.ejogrb.2004.04.029. 1View Article
  23. Faxén M, Nasiell J, Blanck A, Nisell H, Lunell NO: Altered mRNA expression pattern of placental epidermal growth factor receptor (EGFR) in pregnancies complicated by preeclampsia and/or intrauterine growth retardation. Am J Perinatol. 1998, 15 (1): 9-13. 10.1055/s-2007-993890.View Article
  24. Verlohren S, Stepan H, Dechend R: Angiogenic growth factors in the diagnosis and prediction of pre-eclampsia. Clin Sci (Lond). 2012, 122 (2): 43-52. 10.1042/CS20110097.View Article
  25. Herse F, Dechend R, Harsem NK, Wallukat G, Janke J, Qadri F, Hering L, Muller DN, Luft FC, Staff AC: Dysregulation of the circulating and tissue-based renin-angiotensin system in preeclampsia. Hypertension. 2007, 49 (3): 604-611. 10.1161/01.HYP.0000257797.49289.71.View Article
  26. Tsai S, Hardison NE, James AH, Motsinger-Reif AA, Bischoff SR, Thames BH, Piedrahita JA: Transcriptional profiling of human placentas from pregnancies complicated by preeclampsia reveals disregulation of sialic acid acetylesterase and immune signalling pathways. Placenta. 2011, 32 (2): 175-182. 10.1016/j.placenta.2010.11.014.View Article
  27. Winn VD, Gormley M, Paquet AC, Kjaer-Sorensen K, Kramer A, Rumer KK, Haimov-Kochman R, Yeh RF, Overgaard MT, Varki A, Oxvig C, Fisher SJ: Severe preeclampsia-related changes in gene expression at the maternal-fetal interface include sialic acid-binding immunoglobulin-like lectin-6 and pappalysin-2. Endocrinology. 2009, 150 (1): 452-462.View Article
  28. Xie C, Yao MZ, Liu JB, Xiong LK: A meta-analysis of tumor necrosis factor-alpha, interleukin-6, and interleukin-10 in preeclampsia. Cytokine. 2011, 6 (3): 550-559.View Article
  29. Portelinha A, Cerdeira AS, Belo L, Braga J, Tejera E, Pinto A, Pinto F, Areias MJ, Patrício B, Rebelo I: Haemostatic factors in women with history of preeclampsia. Thromb Res. 2009, 124 (1): 52-56. 10.1016/j.thromres.2008.10.005.View Article
  30. Williams D: Long-term complications of preeclampsia. Semin Nephrol. 2011, 31 (1): 111-122. 10.1016/j.semnephrol.2010.10.010.View Article
  31. Roberts JM, Hubel CA: Pregnancy: a screening test for later life cardiovascular disease. Womens Health Issues. 2010, 20–5: 304-307.View Article
  32. Ma H, Henderson KD, Sullivan-Halley J, Duan L, Marshall SF, Ursin G, Horn-Ross PL, Largent J, Deapen DM, Lacey JV, Bernstein L: Pregnancy-related factors and the risk of breast carcinoma in situ and invasive breast cancer among postmenopausal women in the California Teachers Study cohort. Breast Cancer Res. 2010, 12 (3): R35-10.1186/bcr2589.View Article
  33. Calderon-Margalit R, Friedlander Y, Yanetz R, Deutsch L, Perrin MC, Kleinhaus K, Tiram E, Harlap S, Paltiel O: Preeclampsia and subsequent risk of cancer: update from the Jerusalem Perinatal Study. Am J Obstet Gynecol. 2009, 200 (1): 63.e1-5. 10.1016/j.ajog.2008.06.057.View Article
  34. Vatten LJ, Romundstad PR, Jenum PA, Eskild A: Angiogenic balance in pregnancy and subsequent breast cancer risk and survival: a population study. Cancer Epidemiol Biomarkers Prev. 2074, 2009: 18-
  35. Nechuta S, Paneth N, Velie EM: Pregnancy characteristics and maternal breast cancer risk: a review of the epidemiologic literature. Cancer Causes Control. 2010, 21 (7): 967-989. 10.1007/s10552-010-9524-7.View Article
  36. van Dijk M, van Bezu J, Poutsma A, Veerhuis R, Rozemuller AJ, Scheper W, Blankenstein MA, Oudejans CB: The pre-eclampsia gene STOX1 controls a conserved pathway in placenta and brain upregulated in late-onset Alzheimer's disease. J Alzheimers Dis. 2010, 19 (2): 673-679.
  37. Barrett T, Troup DB, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Muertter RN, Holko M, Ayanbule O, Yefanov A, Soboleva A: NCBI GEO: archive for functional genomics data sets-10 years on. Nucleic Acids Res. 2011, 39 (Database issue): D1005-D1010.View Article
  38. Sitras V, Paulssen RH, Grønaas H, Leirvik J, Hanssen TA, Vårtun A, Acharya G: Differential placental gene expression in severe preeclampsia. Placenta. 2009, 30 (5): 424-433. 10.1016/j.placenta.2009.01.012.View Article
  39. Guda P, Chittur SV, Guda C: Comparative analysis of protein-protein interactions in cancer-associated genes. Genomics Proteomics Bioinformatics. 2009, 7 (1–2): 25-36.View Article
  40. Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Statist Soc Ser B. 1995, 57: 289-300.
  41. Mahdavi MA, Lin YH: False positive reduction in protein-protein interaction predictions using gene ontology annotations. BMC Bioinforma. 2007, 23: 8-262.
  42. Cheng D, Knox C, Young N, Stothard P, Damaraju S, Wishart DS: PolySearch: a webbased text mining system for extracting relationships between human diseases, genes, mutations, drugs and metabolites. Nucleic Acids Res. 2008, 36: W399-W405. 10.1093/nar/gkn296.View Article
  43. Hutz JE, Kraja AT, McLeod HL, Province MA: CANDID: a flexible method for prioritizing candidate genes for complex human traits. Genet Epidemiol. 2008, 32 (8): 779-790. 10.1002/gepi.20346.View Article
  44. Radivojac P, Peng K, Clark WT, Peters BJ, Mohan A, Boyle SM, Mooney SD: An integrated approach to inferring gene-disease associations in humans. Proteins. 2008, 72 (3): 1030-1037. 10.1002/prot.21989. 15View Article
  45. Turner B, Razick S, Turinsky AL, Vlasblom J, Crowdy EK, Cho E, Morrison K, Donaldson IM, Wodak SJ: iRefWeb: interactive analysis of consolidated protein interaction data and their supporting evidence. 2010, Oxford: Database, 12-
  46. Cui Q, Ma Y, Jaramillo M, Bari H, Awan A, Yang S, Zhang S, Liu L, Lu M, O'Connor-McCourt M, Purisima EO, Wang E: A map of human cancer signaling. Mol Syst Biol. 2007, 3: 152-View Article
  47. Chung-Yen L, Chia-Hao C, Hsin-Hung W, Shu-Hwa C, Chin-Wen H, Ming-Tat K: Hubba: hub objects analyzer—a framework of interactome hubs identification for network biology. Nucleic Acids Res. 2008, 36: w438-w443. 10.1093/nar/gkn257.View Article
  48. Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T: Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011, 27 (3): 431-432. 10.1093/bioinformatics/btq675.View Article
  49. Kong W, Mou X, Hu X: Exploring matrix factorization techniques for significant genes identification of Alzheimer’s disease microarray gene expression data. BMC Bioinforma. 2011, 12 (Suppl 5): S7-10.1186/1471-2105-12-S5-S7.View Article
  50. Mishra H, Singh N, Misra K, Lahiri T: An ANN-GA model based promoter prediction in Arabidopsis thaliana using tilling microarray data. Bioinformation. 2011, 6 (6): 240-243. 10.6026/97320630006240.View Article
  51. Rio G, Koschützki D, Coello G: How to identify essential genes from molecular networks?. BMC Syst Biol. 2009, 3: 102-10.1186/1752-0509-3-102.View Article
  52. Pang K, Sheng H, Mab X: Understanding gene essentiality by finely characterizing hubs in the yeast protein interaction network. Biochem Biophys Res Commun. 2010, 401: 112-116. 10.1016/j.bbrc.2010.09.021.View Article
  53. Fraley C, Raftery AE: MCLUST: software for model-based cluster analysis. J Classif. 1999, 16: 297-306. 10.1007/s003579900058.View Article
  54. Adamcsek B, Palla G, Farkas IJ, Derényi I, Vicsek T: CFinder: locating cliques and overlapping modules in biological networks. Bioinformatics. 2006, 22: 1021-1023. 10.1093/bioinformatics/btl039.View Article
  55. Huang DW, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37 (1): 1-13. 10.1093/nar/gkn923.View Article

Copyright

© Tejera et al.; licensee BioMed Central Ltd. 2012

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.