Detection of characteristic sub pathway network for angiogenesis based on the comprehensive pathway network
© Huang and Li; licensee BioMed Central Ltd. 2010
Published: 18 January 2010
Pathways in biological system often cooperate with each other to function. Changes of interactions among pathways tightly associate with alterations in the properties and functions of the cell and hence alterations in the phenotype. So, the pathway interactions and especially their changes over time corresponding to specific phenotype are critical to understanding cell functions and phenotypic plasticity.
With prior-defined pathways and incorporated protein-protein interaction (PPI) data, we counted PPIs between corresponding gene sets of each pair of distinct pathways to construct a comprehensive pathway network. Then we proposed a novel concept, characteristic sub pathway network (CSPN), to realize the phenotype-specific pathway interactions. By adding gene expression data regarding a given phenotype, angiogenesis, active PPIs corresponding to stimulation of interleukin-1 (IL-1) and tumor necrosis factor α (TNF-α) on human umbilical vein endothelial cells (HUVECs) respectively were derived. Two kinds of CSPN, namely the static or the dynamic CSPN, were detected by counting active PPIs.
A comprehensive pathway network containing 37 signalling pathways as nodes and 263 pathway interactions were obtained. Two phenotype-specific CSPNs for angiogenesis, corresponding to stimulation of IL-1 and TNF-α on HUVEC respectively, were addressed. From phenotype-specific CSPNs, a static CSPN involving interactions among B cell receptor, T cell receptor, Toll-like receptor, MAPK, VEGF, and ErbB signalling pathways, and a dynamic CSPN involving interactions among TGF-β, Wnt, p53 signalling pathways and cell cycle pathway, were detected for angiogenesis on HUVEC after stimulation of IL-1 and TNF-α respectively. We inferred that, in certain case, the static CSPN maintains related basic functions of the cells, whereas the dynamic CSPN manifests the cells' plastic responses to stimulus and therefore reflects the cells' phenotypic plasticity.
The comprehensive pathway network helps us realize the cooperative behaviours among pathways. Moreover, two kinds of potential CSPNs found in this work, the static CSPN and the dynamic CSPN, are helpful to deeply understand the specific function of HUVEC and its phenotypic plasticity in regard to angiogenesis.
The advent of high-throughput technologies has encouraged the appearance of a lot of analysing tools, which were intended for interpreting gene expression data and extracting biological insight. No longer being limited to produce a ranked list of differentially expressed genes, most of those analysing tools are designed from the perspective of system biology and identify gene sets over-represented in biological processes [1–6]. A gene set could be a group of genes corresponding to a pathway or a group of genes sharing the same Gene Ontology term. Gene Set Enrichment Analysis  and other similar tools are capable to determine pathways which are enriched in the gene list of specific phenotype. Yet these methods fail to furnish the interactions among pathways, not to mention changes of the pathway interactions. Actually, the interactions among pathways are helpful for better understanding of cooperation of pathways which are critical in the functioning of each related pathway. Likewise, changes of the pathway interactions are important for comprehending alterations in the properties and functions of the cell and therefore the phenotypic plasticity, which is cell's ability to alter phenotype in response to specific environmental stimulus [1, 7, 8].
In recent years, several efforts have been made to surmount the limitation of those methods. Global pathway crosstalk network (GPCN) based method  calculates possible crosstalks among pathways and then analyzes microarray data based on these crosstalks to infer sub networks where enriched pathways tied closely to each other as potential modules. Pathway dependency structure (PDS) based method , in the first place, determines pathways involved in two transitions of cancer progression, and successively constructs a pathway network for each transition according to the interdependences of pathways. However, there are still problems to overcome. (i) Although whether interactions among pathways exist could be determined and p-values could be given, the types of interactions are lack in consideration and discussion. GPCN based method did not consider sharing components , which means that different pathways have common components and interact with each other via them, while embedded pathway network based method just considered sharing components . (ii) Generally, the components and their interrelations within distinct pathways are manually collected from literature or other resources by expert curators. Consequently, the importance of context that some genes in specific pathway may be not involved under particular condition must be considered. However, GPCN based method and PDS based method both did not take into account this importance, though gene level analysis of the latter considered it. As a result, these two methods both missed interactions between pairs of pathways at least one of which was not enriched [5, 6]. And GPCN based method even missed interactions which were significant in current case but insignificant in comprehensive pathway network . (iii) The alteration of pathway interactions remains to be studied. GPCN based method did not consider the alteration of interactions . And since PDS based method had no whole possible interdependence network as the base of analysis, nodes of networks constructed with this method were so different that a big part of the real alteration of interactions among involved pathways was still concealed .
Actually, from current literatures, we found that pathway may interact with one other by means of five ways: sharing components , components physically interacting with counterparts from the other pathway [9, 10], components being enzymatic targets of the other pathway , components relating with components of the other pathway via transcription , one pathway generating signal which modulates or competes for a key modulator or mediator of the other . At length, the first way means distinct pathways have common components like genes or gene products and these shared components produce interactions among pathways; the second way depends on the physical protein-protein interactions (PPIs) like binding and phosphorylation; the third way could be reflected by the indirect PPI by two successive enzymes in the metabolic pathway; the fourth way could be represented by indirect PPIs between transcriptional factors and transcribed gene products, which reflect the corresponding gene regulation; the fifth way involves all kinds of PPIs mentioned above.
So we incorporated distinct PPIs to derive interactions among pathways. By counting PPIs between corresponding gene sets of each pair of pathways, pathway interactions could be assessed. Furthermore, since no gene is involved in specific pathway in any case and no PPI works under any condition, we incorporated gene expression data to infer active PPIs, which have at least one differentially expressed corresponding gene, in certain case for distinct phenotype, such as angiogenesis. And then we constructed phenotype-specific pathway networks with the derived sets of active PPIs. Two sub phenotype-specific pathway networks were detected as angiogenesis-related characteristic sub pathway networks (CSPNs) to manifest the specific function of the cell and its phenotypic plasticity with regard to angiogenesis.
Angiogenesis studied in this work, as the process of new blood vessel formation, is fundamental to reproduction, development, and repair. Strictly regulated by many inducers and inhibitors, it could be aroused when the local equilibrium between them is changed. At the phenotype level, angiogenesis is critical to inflammatory responses and cancers . In our work, two series of gene expression data were used to infer CSPNs for angiogenesis. Those data are about human umbilical vein endothelial cells (HUVEC) after stimulations of Interleukin-1 (IL-1) and tumor necrosis factor α (TNF-α), both are known proangiogenic factors. Actually, after being treated with different proangiogenic factors or same factor but with various dosages, some pathway interactions of the cell keep constant, others alter. Those constant interactions are likely required immediately and continually by the cell for basic function in regard to particular stimulus like IL-1 and TNF-α for angiogenesis, whereas those changed interactions are probably reflect the phenotype plasticity concerning specific stimulus such as IL-1 and TNF-α for angiogenesis.
Preparation of PPI data
To complement substantial information loss of restricting the pathway interactions to those by means of physical PPIs like binding or phosphorylation, we considered all the five kinds of pathway interactions. We gathered the PPI data from Human Protein Reference Database  and Kyoto Encyclopedia of Genes and Genomes . HPRD PPIs are mainly physical PPIs. The generalized PPI in KEGG includes three types of interactions: PPI by two successive enzymes in a known pathway, PPI like binding or phosphorylation, PPI between transcription factor and transcribed gene product via gene expression, which actually reflects the gene expression regulation between transcription factor and its regulated gene .
Preparation of pathway data
Generally, established pathways from various databases can be approximately divided into three categories: metabolic pathways, signalling pathways, and pathways of specific diseases or drugs . At length, metabolic pathways focus on the metabolism and metabolic products of cells, having much less relations with PPI. In contrast, signalling pathways mainly involve genes and gene products. Also, the environmental information processing is of biological essence to cells, while the cellular information processing is of clinical consequence: a breakdown in it underlies many diseases . In KEGG, pathways designed for the study of specific diseases are often integrated from some other pathways in former two categories, and pathways from 'Drug development' category are drug structure maps on the structure relationships of drugs. Thus, considering the ways in which pathway interacts with each other are of proteins, we gathered 48 human pathways from 'Environmental Information Processing' and 'Cellular Processes' categories of KEGG as the reference of pathways to elucidate the cooperation of pathways in cells and cells' phenotypic plasticity.
Construction of pathway network and detection of CSPNs
Preparation of time course Microarray data
Since angiogenesis could be aroused by different stimulus, the CSPNs for angiogenesis would be more credible if more data about different stimulus were used. Finally, two series of time course gene expression data, GSE973 and GSE9055, about human umbilical vein endothelial cells (HUVEC) after stimulations of IL-1 and TNF-α were obtained from the Gene Expression Omnibus (GEO). GSE973  has four samples: HUVEC were stimulated with IL-1 for 0, 0.5, 1, 2.5 and 6 hours; and GSE9055  has twenty-five samples: Using TNF-α, samples were collected every 15 min to arrays. We selected samples of 0 h 15 m, 1 h 15 m, 3 h 15 m, 8 h 00 m of GSE9055 in order to compare these two series appropriately. At each time point of these two series, only one sample (microarray) is offered. Consequently, we could not directly make use of coexpression to infer active PPI at teach time point. Also, t test and Bayesian analysis are not appropriate for discovering differentially expressed genes. Finally, the differentially expressed genes were obtained with fold change method. Taking the quality of gene expression data and the shortage of fold change method itself into account, we should admit that genes with high expression level but relatively small fold change may also reflect significant biological meaning, even if their fold change values are smaller than 2. Hence, we inferred that all the immediate neighbours of those differentially expressed genes in the PPI network were underlying relevant genes, for their approved interactions with those differentially expressed genes support the presumption of their relevance to some extent. Then we considered those PPIs which have at least one differentially expressed corresponding gene as the active PPIs. Successively, we exploited this kind of active interaction in construction of phenotype-specific pathway network.
Comprehensive pathway network construction
The top 16 Hub pathways (degree > 18) identified in Comprehensive pathway network
B cell receptor signaling pathway
ErbB signaling pathway
Fc epsilon RI signaling pathway
VEGF signaling pathway
GnRH signaling pathway
Insulin signaling pathway
T cell receptor signaling pathway
Phosphatidylinositol signaling system
Natural killer cell mediated cytotoxicity
Regulation of actin cytoskeleton
MAPK signaling pathway
mTOR signaling pathway
Types of pathway interaction
Comparison of degrees of different pathway networks
Number of Nodes
Number of Edges
Comprehensive pathway network based on all types of pathway interactions
Comprehensive pathway network only based on sharing components
Comprehensive pathway network based on interactions without sharing components
Phenotype-specific pathway network in IL-1 case for angiogenesis, based on comprehensive pathway network
Phenotype-specific pathway network in TNF-α case for angiogenesis, based on comprehensive pathway network
Phenotype-specific pathway network in IL-1 case for angiogenesis, based on enriched pathway pipeline
Phenotype-specific pathway network construction
Quality of pathway interaction in phenotype-specific pathway network
If we take into account the interactions among pathways from a narrow perspective of enriched pathways, we will miss important interactions. There are two typical enriched pathway based approaches, both of which have this problem [5, 6]. Therefore, we argue that it is more appropriate to construct pathway network for relevant case through assessing the significance of pathway interaction, no matter whether the pathways are enriched in current case or not. Actually, this approach is much better for comprehending the cooperation between distinct pathways, since the enrichment analysis only reveals to which degree the whole but not the relevant part of gene set of pathway takes part in current case. We built a phenotype-specific pathway network in IL-1 case based on enriched pathway pipeline to illustrate its shortcoming (Table 2). There were only 15 interactions involving 10 pathways in this network. Important interactions such as interaction between Wnt and TGF-β signalling pathways disappeared.
Detecting CSPNs for angiogenesis
From two constructed phenotype-specific pathway networks, static CSPN and dynamic CSPN for angiogenesis were detected. On the one hand, firstly, the common static sub pathway network within which interactions keep constant were identified. Then a part of this sub network involving B cell receptor, T cell receptor, Toll-like receptor, MAPK, VEGF and ErbB pathways was detected as static CSPN (Figure 2). The high degrees of those pathways and the almost full connected nature of this part of network were maintained at each time point in both cases. These pathways are highly connected in comprehensive pathway network too, which means their wide and deep involvement in different phenotypes. Also, B cell receptor, T cell receptor and Toll-like receptor signalling pathways are important in immunity; MAPK pathway is involved in various cellular functions; VEGF pathway is considered to be crucial in signal transduction in both physiologic and pathologic angiogenesis; ErbB pathway couples binding of extracellular growth factor ligands to intracellular signalling pathways regulating diverse biologic responses [13, 18]. Therefore, we derived that, conforming to interactions in the static CSPN, pathways in static CSPN cooperated with each other to receive, process and relay extracellular stimulus for angiogenesis. On the other hand, there was another group of pathways including TGF-β, Wnt, p53 pathways and cell cycle pathway, with interactions among them was detected as dynamic CSPN (Figure 2). Actually, those interactions changed observably over time in both series. Their relatively low degrees in comprehensive pathway interaction network might mean that they only involved in specific phenotypes. In addition, these four pathways all take part in the control of cell fate [13, 19]. We inferred that the cooperation of them played a key role in manifesting response of HUVEC to stimulus like IL-1 and TNF-α via controlling development of HUVEC, and reflected the phenotypic plasticity of HUVEC in regard to angiogenesis.
In detail, interactions among pathways within static CSPN are the immediate and constant response of HUVEC to the stimulation of IL-1 or TNF-α signal. The B cell receptor, T cell receptor and Toll-like receptor signalling pathways represented the innate immune system to depress the excitement from stimulus . They cooperated with each other to interfere with the interactions among MAPK, ErbB and VEGF pathways, which likely promoted angiogenesis [18, 21, 22]. Also, within dynamic CSPN, p53 pathway was possibly aroused to prevent angiogenesis via negatively regulating cell cycle pathway . But, cooperating with each other to interact with cell cycle pathway, TGF-β and Wnt pathways inside dynamic CSPN also were stirred up to promote angiogenesis [24, 25].
We compared the dynamic CSPN at each time point too (Additional file 2: Comparison of dynamic CSPN at different time points). In IL-1 case, p53 only interacted with cell cycle pathway at second and fourth time points. The first interaction might respond to the excitement from static CSPN, since MAPK, ErbB and VEGF pathway immediately cooperated with each other to promote angiogenesis after treatment. The fourth interaction probably reacted to the agitation from dynamic CSPN itself, for TGF-β and Wnt pathways collaborated to propel angiogenesis after the spread of signals from static CSPN. Actually, Wnt pathway interacted with VEGF pathway at the third time point (Additional file 3: Pathway interaction at different time point(s) in IL-1 case), which likely invoked the collaboration between Wnt and TGF-β pathways as well as their regulation to cell cycle pathway. In TNF-α case, at last three time points, p53 all interacted with cell cycle pathway, as well as TGF-β and Wnt pathways cooperated to propel angiogenesis by interacting with cell cycle pathway (see Additional file 4). Virtually, all the interactions between pathways within dynamic CSPN and pathways from static CSPN were appeared after the first time point, which possibly implied the transmission of signals from static CSPN .
We inferred that the difference of evolvement of dynamic CSPN between IL-1 case and TNF-α case probably reflected the diversity of traits or dosages of stimulus, though uncertainty of the cause remained before more reliable data could be gathered for this pipeline. Also, there was cooperation between static and dynamic CSPNs. The static CSPN was likely the core one, because of its immediate and sustained response to stimulus, while the dynamic CSPN was probably the complementary one, by reason of its same effect to stimulus as the static CSPN as well as the delay of its response. It is reported that crosstalk of VEGF and Notch pathways is crucial to tumor angiogenesis , resulting in blockade of VEGF pathway sometimes doesn't work . In angiogenesis of HUVEC, it is also possible that blocking VEGF pathway make no effect for prevention of angiogenesis, on account of the work of TGF-β and Wnt pathways in dynamic CSPN.
In this work, we mainly proposed a methodology to consider the alteration of interactions among relevant pathways for specific phenotype, and then to detect two kinds CSPNs from constructed phenotype specific pathway networks of different cases for the same certain phenotype. We took angiogenesis on HUVEC as an example. In the past years, many useful methods, which could identify gene sets over-represented in biological processes, were offered [2, 3]. For example, LMMA reconstructed the gene association network specific for angiogenesis, and extracted significant pathways from the network . For moving forward to the next stage, the interactions among pathways and the alterations of these interactions require further investigation. Recently, GPCN based  and PDS based  methods were proposed to follow up this subject. But three problems still existed: the types of interactions were lack in consideration and discussion; the importance of context needed to be taken into account; the alteration of pathway interactions remained to be studied.
To solve these problems, we considered all five types of pathway interactions to construct a comprehensive pathway network, and inferred relevant genes of each pathway and active PPIs with Microarray data. Then, phenotype-specific pathway networks were constructed. From these networks, alteration of pathway interactions over time was distinct in each network, and alteration of pathway interactions with different stimulus was clear after comparison of corresponding networks. A static CSPN and a dynamic CSPN were also detected from these phenotype-specific pathway networks. The effect and cooperating mechanism of distinct pathways within both CSPNs, the influence of the diversity of traits or dosages of stimulus to the evolvement of dynamic CSPN, as well as the complementary role of dynamic CSPN to static CSPN, were elaborated.
Actually, it will enhance the reliability of "active PPI" by requiring corresponding genes firstly are coexpressed within distinct case (series in GEO). We will go further in the future to refine our work from this respect. In addition, the detection of CSPNs for specific phenotype requires ruling out the effect of other factors such as the difference of stimulus and experimental bias. As a result, the quantity and quality of corresponding time-course microarray data are asked to be appropriate to obtain reliable CSPNs. But now even the demand of quantity can not be met. We also will go further in the future to provide related active PPIs for each pathway interaction and then predict the pathway interaction type.
The comprehensive pathway network helps us realize the cooperative behaviour among pathways by providing most of possible interactions among pathways. By combining phenotype-specific gene expression data, two kinds of potential CSPNs detected in this work, the static CSPN and the dynamic CSPN, are helpful to deeply understand the specific function and phenotypic plasticity of angiogenesis.
We thank Dr. Gu J, Zheng RH and Hu Z in our lab for the kind help. This work is supported by NSFC (Nos. 30873464, 60934004 and 60721003) and FANEDD (No. 200366).
This article has been published as part of BMC Bioinformatics Volume 11 Supplement 1, 2010: Selected articles from the Eighth Asia-Pacific Bioinformatics Conference (APBC 2010). The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/11?issue=S1.
- Wang X, Dalkic E, Wu M, Chan C: Gene module level analysis: identification to networks and dynamics. Curr Opin Biotechnol 2008, 19(5):482–491. 10.1016/j.copbio.2008.07.011PubMed CentralView ArticlePubMedGoogle Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al.: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA 2005, 102: 15545–15550. 10.1073/pnas.0506580102PubMed CentralView ArticlePubMedGoogle Scholar
- Li S, Wu L, Zhang Z: Constructing biological networks through combined literature mining and microarray analysis: a LMMA approach. Bioinformatics 2006, 22(17):2143–2150. 10.1093/bioinformatics/btl363View ArticlePubMedGoogle Scholar
- Wu X, Jiang R, Zhang M, Li S: Network-based global inference of human disease genes. Mol Syst Biol 2008, 4: 189. 10.1038/msb.2008.27PubMed CentralView ArticlePubMedGoogle Scholar
- Li Y, Agarwal P, Rajagopalan D: A global pathway crosstalk network. Bioinformatics 2008, 24(12):1442–1447. 10.1093/bioinformatics/btn200View ArticlePubMedGoogle Scholar
- Edelman EJ, Guinney J, Chi JT, Febbo PG, Mukherjee S: Modeling cancer progression via pathway dependencies. PLoS Comput Biol 2008, 4(2):e28. 10.1371/journal.pcbi.0040028PubMed CentralView ArticlePubMedGoogle Scholar
- Feinberg AP: Phenotypic plasticity and the epigenetics of human disease. Nature 2007, 447(7143):433–440. 10.1038/nature05919View ArticlePubMedGoogle Scholar
- Swynghedauw B: Phenotypic plasticity of adult myocardium: molecular mechanisms. J Exp Biol 2006, 209(12):2320–2327. 10.1242/jeb.02084View ArticlePubMedGoogle Scholar
- Lu LJ, Sboner A, Huang YJ, Lu HX, Gianoulis TA, Yip KY, Kim PM, Montelione GT, Gerstein MB: Comparing classical pathways and modern networks: towards the development of an edge ontology. Trends Biochem Sci 2007, 32(7):320–331. 10.1016/j.tibs.2007.06.003View ArticlePubMedGoogle Scholar
- Guo X, Wang XF: Signaling crosstalk between TGF-/BMP and other pathways. Cell Res 2009, 19: 71–88. 10.1038/cr.2008.302PubMed CentralView ArticlePubMedGoogle Scholar
- Folkman J: Angiogenesis in cancer, vascular, rheumatoid and other disease. Nat Med 1995, 1(1):27–31. 10.1038/nm0195-27View ArticlePubMedGoogle Scholar
- Keshava Prasad TS, Goel R, Kandasamy K, Keerthikumar S, Kumar S, Mathivanan S, Telikicherla D, Raju R, Shafreen B, Venugopal A, et al.: Human Protein Reference Database - 2009 update. Nucleic Acids Res 2009, 37: D767–72. 10.1093/nar/gkn892PubMed CentralView ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res 2000, 28: 27–30. 10.1093/nar/28.1.27PubMed CentralView ArticlePubMedGoogle Scholar
- Fraser ID, Germain RN: Navigating the network: signaling crosstalk in hematopoietic cells. Nat Immunol 2009, 10(4):327–331. 10.1038/ni.1711PubMed CentralView ArticlePubMedGoogle Scholar
- Castillo-Davis CI, Hartl DL: GeneMerge--post-genomic analysis, data mining, and hypothesis testing. Bioinformatics 2003, 9(7):891–892. 10.1093/bioinformatics/btg114View ArticleGoogle Scholar
- Mayer H, Bilban M, Kurtev V, Gruber F, Wagner O, Binder BR, de Martin R: Deciphering regulatory patterns of inflammatory gene expression from interleukin-1-stimulated human endothelial cells. Arterioscler Thromb Vasc Biol 2004, 24(7):1192–1198. 10.1161/01.ATV.0000131263.06296.77View ArticlePubMedGoogle Scholar
- Series GSE9055, Time course gene expression of HUVEC after TNF-α treatment, GEO[http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE9055]
- Yoo SA, Kwok SK, Kim WU: Proinflammatory role of vascular endothelial growth factor in the pathogenesis of rheumatoid arthritis: prospects for therapeutic intervention. Mediators Inflamm 2008, 2008: 129873. 10.1155/2008/129873PubMed CentralView ArticlePubMedGoogle Scholar
- Roarty K, Baxley SE, Crowley MR, Frost AR, Serra R: Loss of TGF-β or Wnt5a results in an increase in Wnt/β-catenin activity and redirects mammary tumor phenotype. Breast Cancer Res 2009, 11(2):R19. 10.1186/bcr2244PubMed CentralView ArticlePubMedGoogle Scholar
- Pasare C, Medzhitov R: Control of B-cell responses by Toll-like receptors. Nature 2005, 438(7066):364–8. 10.1038/nature04267View ArticlePubMedGoogle Scholar
- Le XF, Mao W, Lu C, Thornton A, Heymach JV, Sood AK, Bast RC Jr: Specific blockade of VEGF and HER2 pathways results in greater growth inhibition of breast cancer xenografts that overexpress HER2. Cell Cycle 2008, 7(23):3747–58.PubMed CentralView ArticlePubMedGoogle Scholar
- Curry JM, Eubank TD, Roberts RD, Wang Y, Pore N, Maity A, Marsh CB: M-CSF signals through the MAPK/ERK pathway via Sp1 to induce VEGF production and induces angiogenesis in vivo. PLoS One 2008, 3(10):e3405. 10.1371/journal.pone.0003405PubMed CentralView ArticlePubMedGoogle Scholar
- Teodoro JG, Parker AE, Zhu X, Green MR: p53-mediated inhibition of angiogenesis through up-regulation of a collagen prolyl hydroxylase. Science 2006, 313(5789):968–71. 10.1126/science.1126391View ArticlePubMedGoogle Scholar
- Ferrari G, Cook BD, Terushkin V, Pintucci G, Mignatti P: Transforming growth factor-beta 1 (TGF-beta1) induces angiogenesis through vascular endothelial growth factor (VEGF)-mediated apoptosis. J Cell Physiol 2009, 219(2):449–58. 10.1002/jcp.21706PubMed CentralView ArticlePubMedGoogle Scholar
- Suomalainen M, Thesleff I: Patterns of Wnt pathway activity in the mouse incisor indicate absence of Wnt/beta-catenin signaling in the epithelial stem cells. Dev Dyn, in press.Google Scholar
- Li JL, Harris AL: Crosstalk of VEGF and Notch pathways in tumour angiogenesis: therapeutic implications. Front Biosci 2009, 14: 3094–110. 10.2741/3438View ArticleGoogle Scholar
- Thurston G, Kitajewski J: VEGF and Delta-Notch: interacting signalling pathways in tumour angiogenesis. Br J Cancer 2008, 99(8):1204–9. 10.1038/sj.bjc.6604484PubMed CentralView ArticlePubMedGoogle Scholar
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.