- Open Access
Prioritization of genes involved in endothelial cell apoptosis by their implication in lymphedema using an analysis of associative gene networks with ANDSystem
BMC Medical Genomicsvolume 12, Article number: 47 (2019)
Currently, more than 150 million people worldwide suffer from lymphedema. It is a chronic progressive disease characterized by high-protein edema of various parts of the body due to defects in lymphatic drainage. Molecular-genetic mechanisms of the disease are still poorly understood. Beginning of a clinical manifestation of primary lymphedema in middle age and the development of secondary lymphedema after treatment of breast cancer can be genetically determined. Disruption of endothelial cell apoptosis can be considered as one of the factors contributing to the development of lymphedema. However, a study of the relationship between genes associated with lymphedema and genes involved in endothelial apoptosis, in the associative gene network was not previously conducted.
In the current work, we used well-known methods (ToppGene and Endeavour), as well as methods previously developed by us, to prioritize genes involved in endothelial apoptosis and to find potential participants of molecular-genetic mechanisms of lymphedema among them. Original methods of prioritization took into account the overrepresented Gene Ontology biological processes, the centrality of vertices in the associative gene network, describing the interactions of endothelial apoptosis genes with genes associated with lymphedema, and the association of the analyzed genes with diseases that are comorbid to lymphedema.
An assessment of the quality of prioritization was performed using criteria, which involved an analysis of the enrichment of the top-most priority genes by genes, which are known to have simultaneous interactions with lymphedema and endothelial cell apoptosis, as well as by genes differentially expressed in murine model of lymphedema. In particular, among genes involved in endothelial apoptosis, KDR, TNF, TEK, BMPR2, SERPINE1, IL10, CD40LG, CCL2, FASLG and ABL1 had the highest priority. The identified priority genes can be considered as candidates for genotyping in the studies involving the search for associations with lymphedema.
Analysis of interactions of these genes in the associative gene network of lymphedema can improve understanding of mechanisms of interaction between endothelial apoptosis and lymphangiogenesis, and shed light on the role of disturbance of these processes in the development of edema, chronic inflammation and connective tissue transformation during the progression of the disease.
Lymphedema is a chronic progressive disease, resulting in a significant loss of productivity, which affects more than 150 million people worldwide . At the heart of the pathogenesis of lymphedema lies the appearance of chronic high-protein edema, characterized by the pathological accumulation of intercellular fluid with large-molecule proteins in the interstitial space due to a defect in lymphatic drainage caused by congenital malformation (primary lymphedema), lymphatic obstruction or by the destruction of lymphatic vessels (secondary lymphedema). Long existence of the high-protein edema causes chronic inflammation leading to the replacement of adipose tissue with connective tissue, increasing the volume of connective tissue matrix, which subsequently leads not only to an increase in size of body parts, but also to a secondary disturbance of lymphatic transport and drainage . Most often it affects the lower extremities, but can also affect the upper limbs, face, trunk, external genitalia, etc. . Primary lymphedema is a disease caused by dysfunction of lymphatic vessels, their aplasia, dysplasia and hypoplasia. Disturbance of drainage function of lymphatic vessels leads to accumulation of fluid rich in proteins under the epidermis. Most often, edema is localized on the lower limbs. Family forms of primary lymphedema are most often an autosomal dominant disease, but autosomal recessive types of inheritance of this pathology are also found . It is believed that developmental defects and the resulting dysfunction of the lymphatic system are the cause of primary lymphedema, as well as associated syndromes, and can be genetically determined . Secondary lymphedema is most often the result of a complex treatment of breast cancer, including removal of the axillary lymph nodes and/or radiation therapy . Recently a number of associations between gene polymorphisms, including cytokine genes , and the development of breast cancer–related lymphedema were identified [8, 9]. One of the possible causes of clinical manifestation of primary lymphedema in middle age and the development of secondary lymphedema after treatment of breast cancer may be apoptosis of the endothelium, whose role is discussed in publications on the pathogenesis of lymphedema [10, 11].
Apoptosis is a form of cell death, characterized by a number of morphological and molecular features, including exposure of phosphatidylserine on the cell membrane, blebbing of the plasma membrane, cell shrinkage, cytoskeleton rearrangement, nuclear collapse, chromatin condensation, DNA fragmentation and formation of “apoptotic bodies” . A number of works discuss the role of endothelial cell apoptosis in the pathogenesis of lymphedema [10, 11, 13, 14]. However, the association between genes involved in endothelial apoptosis and genes associated with lymphedema has not been yet analyzed.
One of the most effective and frequently used approaches for the identification of the potential associations between genes and diseases is a prioritization . Nowadays, there are many known methods for the prioritization of genes. A first, and the largest group of such methods is based on the analysis of genomic and transcriptomic data, as well as data on homologues . The second class of prioritization methods involves approaches based on the analysis of gene networks and protein-protein interactions networks. An example of such methods is GUILD . Integrated methods that combine approaches from the first and the second groups represent the third group. The well-known examples of such systems from the third group are ToppGene  and Endeavour .
Previously, we developed criteria for the prioritization of genes, by using the well-known systems and also considering the structure of associative gene networks from ANDSystem [19, 20]. ANDSystem is a computer tool, designed for the automated extraction of knowledge from the texts of scientific publications and automatic reconstruction of the associative gene networks by using the retrieved information, describing the mechanisms of diseases, as well as other complex phenotypic traits. The knowledge base of ANDSystem contains over 30 million facts describing genetic regulation, gene associations with diseases, protein-protein interactions, catalytic reactions, transport pathways, etc., extracted from more than 25 million PubMed abstracts [21, 22]. In particular, ANDSystem was used for the identification of candidate genes associated with comorbidity of preeclampsia, diabetes and obesity , asthma and tuberculosis , as well as asthma and hypertension .
In this study, our prioritization criteria were applied to identify the endothelial apoptosis-related genes, potentially involved in lymphedema. Three approaches were used to assess the quality of prioritization, based on the evaluation of the enrichment of the list of top priority genes by genes for which there was indirect evidence confirming their association with lymphedema: (1) enrichment with genes shared between lymphedema and endothelial apoptosis; (2) enrichment with genes, the names of which are significantly co-occurring in full-text articles with the key word «lymphedema»; (3) enrichment with genes differentially expressed in the murine model of lymphedema. All these quality criteria showed significant enrichment.
Among the genes with the highest priority TNF, TEK, BMPR2, SERPINE1, IL10, CD40LG, CCL2, FASLG and ABL1 can be distinguished. These genes can be used to plan experiments confirming their association with lymphedema, and to be further considered as new candidates for genotyping. The found genes can become the basis for a better understanding of the mechanisms of interaction between endothelial apoptosis and lymphedema.
The list of genes associated with lymphedema (Additional file 1: Table S1) was extracted from the CTD , Malacards , KEGG , HPO  and DisGeNET  databases, available for January 2018. This list was additionally expanded with genes associated with lymphedema, extracted from ANDSystem .
The over-represented Gene Ontology (GO) biological processes were identified using the DAVID 6.8 tool  with the following parameters: the organism – «Homo sapiens», Gene_Ontology – «GOTERM_BP_DIRECT».
The betweenness centrality of a node in a gene network was estimated using the networkx package of the Python programming language . This indicator characterizes the number of shortest pathways between all pairs of vertices of analyzed graph passing through a given vertex and reflects the functional significance of gene in gene network.
The Mann-Whitney criterion was calculated using the mannwhitneyu function of the scipy.stats package of the Python programming language .
The list of human genes involved in the Gene Ontology biological process «apoptotic process» was obtained using the AmiGO system  by the «GO:0006915» query and Organism filter set to «Homo sapiens».
The list of human genes involved in the Gene Ontology biological process «endothelial cell apoptotic process» was obtained with AmiGO using the «GO:0072577» query and with Organism filter set to «Homo sapiens».
Criteria 1 and 2 were the well-known prioritization methods: ToppGene  and Endeavour version 3.71 , respectively. These resources allow to perform ranking of a test set of genes by a training set of genes according to specific criteria characterizing the proximity of genes from the test set to the genes from the training sample. The methods of these resources use the genetic information (co-localization in the genome), the functional properties of genes (participation in the same GO categories), etc., as well as properties of the vertices of the graph of protein-protein networks. A list of genes associated with lymphedema, described above, was used as a training set for each of these methods. As a test set, the list of genes involved in endothelial cell apoptosis, described above, was used. For the ToppGene the “all Feature” parameter was selected in the “Training parameters” section, and the ranking of genes was based on the value of the “Rank” output parameter. In case of the Endeavour system all settings were set to default, the gene ranking was based on “P-value”. Thus, the lowest ranks had genes with the lowest “P-value”, while genes with the highest “P-value” obtained the highest ranks.
Criterion 3 was calculated as the proportion of the biological processes, over-represented for a set of genes associated with lymphedema (Additional file 2: Table S2), to all Gene Ontology biological processes where the analyzed gene was involved. Information on the involvement of gene in the Gene Ontology biological process was obtained from the AmiGO system . Ranks for genes were determined with the sorting of the list of genes by descending the proportion of over-represented Gene Ontology biological processes. Thus, the lowest ranks were assigned to genes with the largest proportion of over-represented Gene Ontology biological processes.
Criterion 4 was based on the use of the cross-talk centrality (CTC), calculated using the «Intelligent Filtration» function of ANDSystem. Within this criterion, CTCs were calculated separately for genes from the gene network (criterion 4A) and for their products (criterion 4B). Thus, two indices (CTCgene and CTCprotein) were determined for each gene. The centrality of CTCgene was calculated using the following formula:
where Ni – is a number of interactions of the i-th gene with members of the associative gene network of lymphedema; M – is a total number of nodes from the associative gene network of lymphedema. CTCprotein was determined for the i-th protein or miRNA in a same way. CTCgene and CTCprotein ranks were determined by sorting the values in descending order the same way as it was done for the criterion 3.
Criterion 5 represents the number of diseases comorbid to lymphedema, which are associated with the analyzed gene. Comorbid diseases are considered to be simultaneously present in one patient more often than can be expected for accidental reasons [35, 36]. The list of diseases comorbid to lymphedema was manually created by analyzing the publications corresponding to the following query to the PubMed database: “lymphedema and (comorbid or comorbidity)”. A total of 80 publications were manually analyzed and six comorbid diseases were found (Additional file 3: Table S3). All interactions between the analyzed genes and these six comorbid diseases were established using ANDSystem. In addition, the rank by criterion 5 was determined by sorting the list of genes in descending order of the number of comorbid diseases associated with these genes.
In the cases when several genes had the same value of the criterion by which the ranking was performed, the rank of these genes was calculated similarly to the Spearman rank correlation coefficient . I.e., such genes were assigned with the same rank, equal to the average arithmetic ranks of these genes, according to their position, in the sorted list of genes.
Criterion 6 was calculated as the average value of ranks obtained from the criteria 1–5.
The statistical significance of the co-occurrence of gene names from the test set and lymphedema in full-text articles was estimated as follows:
(1) The number of full-text publications presented in the PubMed Central system (https://www.ncbi.nlm.nih.gov/pmc/) was calculated according to the following parameters:
(a) the official name of the analyzed gene is co-occurring with the “lymphedema” word,
(b) the official name of the analyzed gene is mentioned in the article,
(c) the word “lymphedema” is mentioned in the article.
(2) Using the hypergeometric distribution, implemented in the hypergeom.sf function of the scipy.stats package of the Python programming language , the statistical significance of the co-occurrence of names of the analyzed genes and lymphedema was assessed.
The correlation between the ranks of criteria 1–6 and the significance of co-occurrence of the analyzed genes names and lymphedema in the full-text articles was carried out using the Spearman’s rank correlation coefficients. The calculations were performed using the Social Science Statistics system (http://www.socscistatistics.com/tests/spearman/default2.aspx).
Data on gene expression in the mouse surgical model of acute lymphedema induction were collected from the GEO resource (https://www.ncbi.nlm.nih.gov/geo/) using the «GSE4333» id . Differentially expressed genes (p-value < 0.05) for the following groups: (1) «Lymphedema Tail Skin» (test) versus «Normal Tail Skin» (control 1); and (2) «Lymphedema Tail Skin» versus «Surgical Sham Control Tail Skin» (control 2), were identified using the GEO2R tool (http://www.ncbi.nlm.nih.gov/geo/geo2r/). Mapping of mouse genes on human orthologous genes was carried out using the HomoloGene database (https://www.ncbi.nlm.nih.gov/homologene).
Results and discussion
Associative gene network of lymphedema
Analysis of information from the CTD , Malacards , KEGG , HPO  and DisGeNET  databases and from ANDSystem, revealed 69 genes associated with lymphedema (Additional file 1: Table S1). The obtained list of genes was applied to ANDSystem for the reconstruction of the associative gene network of lymphedema. The reconstructed network was automatically expanded in ANDSystem with the products of these genes. The obtained gene network contained 69 genes and 78 proteins, as well as 709 interactions between them (Fig. 2), including the following interaction types: association (460 interactions), protein-protein interaction (128 interactions), gene expression (78 interactions), regulation of gene expression (30 interactions), co-expression (6 interactions), regulation of transport (2 interactions), regulation of degradation (2 interactions), catalysis (2 interactions) and regulation of activity (1 interaction). Twenty one genes from the network were interacting only with their own products, while the remaining 48 genes and 57 proteins were highly connected with each other (Fig. 2).
Analysis of the overrepresentation of Gene Ontology (GO) biological processes for a set of genes associated with lymphedema revealed 40 processes significantly enriched with these genes (Additional file 2: Table S2). Among the most over-represented processes were lymphangiogenesis, endothelial cell proliferation, ERK1/ERK2 cascade and VEGF signaling pathway. The role of these biological processes in the pathogenesis of lymphedema is actively discussed in the literature [8, 40,41,42]. For example, Saito et al., 2013 studied the possibility of therapeutic lymphangiogenesis, using hepatocyte growth factor for the treatment of lymphedema . Miaskowski et al., 2013, described the important role of VEGF signaling pathway in lymphangiogenesis associated with inflammation in lymphedema .
It should be noted that among the significantly over-represented GO biological processes were also processes of positive and negative regulation of the apoptosis (Additional file 2: Table S2); this is consistent with recent studies suggesting the important role of apoptosis in the pathogenesis of lymphedema . The associative gene network of lymphedema involved 24 genes/proteins from the GO biological process entitled “apoptotic process” (Fig. 2). Analysis of the centrality of nodes from the gene network of lymphedema (Additional file 4: Table S4) showed that the average value of the betweenness centrality for apoptosis genes (equal to 303) statistically significantly exceeds the value calculated for all the nodes of the gene network (equal to 128), with p-value< 10− 6, according to the Mann-Whitney criterion. The high centrality of the apoptotic genes in the lymphedema gene network can be indirect evidence of the key role of the apoptosis process in the pathogenesis of lymphedema.
Prioritization of endothelial apoptosis genes by their potential association with lymphedema
Considering the importance of apoptosis of endothelial cells in the pathogenesis of lymphedema [10, 11, 13, 14], we performed a prioritization of genes involved in the «endothelial cell apoptotic process» (GO: 0072577) using a training set that included genes associated with lymphedema. The prioritization process was based on the use of well-known methods (ToppGene and Endeavour), as well as on our original approaches, allowing to consider the structure of the gene network of lymphedema, over-represented biological processes and associations of genes with diseases comorbid to lymphedema. It appeared that two genes GATA2 and KDR were simultaneously presented in the list of genes from the «endothelial cell apoptotic process», as well as in the list of genes associated with lymphedema [43, 44]. During the prioritization, these two genes were excluded from the training set and were used as control genes. Thus, the prioritization was performed for 64 genes from the testing set and 67 genes from the training set, associated with lymphedema.
According to the criterion 6 (Additional file 5: Table S5), calculated as an average value of criteria 1–5, it turned out that the KDR control gene was on the first place, while the GATA2 gene was on twelfth. Hypergeometric distribution showed that the enrichment of the top 12 genes by these two control genes is statistically significant (p-value = 0.03). It should be noted that a similar analysis, carried out by individual criteria 1–5, showed the absence of any statistically significant enrichment. Thus, the criterion 6, which takes into account five previous criteria, can be considered as the best approach to the genes prioritization.
In the Table 1 are shown top 10 of the highest priority genes according to the criterion 6. Among these top 10 genes SNPs associated with lymphedema are known only for KDR gene , which had the first place in the table. For the other genes, to our knowledge, there was no information about any SNPs associated with this disease and these genes can be considered as new candidates for lymphedema susceptibility.
Also, it was interesting to consider each remaining criterion separately, even though none of them showed any statistical significance during the evaluation of the enrichment of the top of genes, ranked by these criteria, by the control genes. Thus, both methods of prioritization (ToppGene and Endeavour) determined the KRD control gene in the first place in the resulting list. The top 5 of the most priority genes according to both ToppGene (Additional file 6: Table S6) and Endeavour (Additional file 7: Table S7) were KDR, ABL1 and TEK, presented in Table 1. Using the criterion 3 the top 10 included the control gene KDR and the TEK gene (Additional file 8: Table S8). At the same time, according to criterion 4A (Additional file 9: Table S9), the top 10 included eight genes (KDR, TNF, IL10, CCL2, CD40LG, ABL1, BMPR2 and SERPINE1) from the Table 1. While for the criterion 4B (Additional file 10: Table S10) there were seven such genes in a list (TNF, IL-10, KDR, CD40LG, CCL2, ABL1 and BMPR2). According to the criterion 5, which considers the comorbidity of diseases, the TNF gene, associated with six comorbid diseases, was in the first place (Additional file 11: Table S11), the second priority was given to CCL2, IL10 and SERPINE1 genes, associated with five comorbid diseases. These genes also appeared to be in the top 10 list (Table 1).
Verification of the prioritization results using the full-text articles
Jenssen et al., 2001 proposed an approach for automated identification of potential interactions between biological objects, based on an assessment of the statistical significance of the co-occurrence of terms in scientific publications . Because this approach is not implemented in ANDSystem the estimation of the co-occurrence between the analyzed genes and lymphedema can serve as an indirect evidence of the correctness of our prioritization results. It can be expected that the most priority genes would have the highest frequency of mentioning in articles together with lymphedema. As the ANDSystem knowledge base was created by automated analysis of PubMed, to obtain results that could be even more free of the ANDSystem data, we performed a manual analysis of the full-text articles from PubMed Central (Additional file 12: Table S12). It was found that of the 64 of analyzed genes, the 21 genes were statistically significantly more often mentioned in publications together with the “lymphedema” word (p-value< 0,05 with FDR correction). 13 out of these 21 genes appeared to be presented among the top 20 genes, ranked according to criterion 6 (Table 2), while the most significant association with lymphedema was observed for the GATA2 and KDR control genes (p-value = 10− 195 and p-value = 10− 136, respectively).
According to the hypergeometric distribution, the statistical significance of such enrichment of the top 20 genes had a p-value< 0.0004. It should be noted that only criterion 5, which considers the comorbidity, also had such a low p-value, while the p-values of other criteria were higher. At the same time, the highest Spearman’s rank correlation coefficient between the genes ranks and the significance of co-occurrence of the analyzed genes names and lymphedema in the full-text articles (r = 0,529), was observed for the criterion 6 (Table 2).
Agreement between the prioritization results and gene expression data from the mouse lymphedema model
Lymphedema is characterized by a chronic stasis of lymph in the tissues. In Tabibiazar et al., 2006, the authors performed gene expression profiling by an array in the experimental model of acute postoperative lymphedema associated with lymphatic stagnation in the tails of SKH-1 mice . The following conditions were examined: (1) lymphedema tail skin caused by surgical lymphatic vessel blockage (test); (2) normal tail skin with no intervention (control 1); (3) surgical sham control tail skin - surgical incision with no lymphatic vessel blockage (control 2). Using the data from Tabibiazar et al., 2006  (GEO series GSE4333) 10201 mouse microarray probes differentially expressed between test and control 1, and 2219 probes for test versus control 2 with p-value< 0.05 (Additional file 13: Table S13) were detected by GEO2R (http://www.ncbi.nlm.nih.gov/geo/geo2r/). After the mapping of these differentially expressed probes on orthologues human genes from the 10,201 probes 3289 differentially expressed human genes were remained, while of the 2219 probes remained 735 genes. The remaining probes were removed from the analysis due to the lack of any intersections with human genes. After combining these two sets 3494 differentially expressed genes were obtained. It appeared that 13 of 64 analyzed genes were in the combined set.
According to criterion 6, among the top 20 of the highest priority genes there were 8 of 13 differentially expressed genes (KDR, PLCG1, SERPINE1, CD40LG, IL10, CCL2, PDPK1 and THBS1). It was shown by hypergeometric distribution that such enrichment is statistically significant (p-value = 0.012, Table 3). However, the greatest enrichment value was observed for criterion 4B. This criterion is based on an assessment of the CTC centrality of the nodes of the gene network, corresponding to proteins.
Interactions of the top 10 candidates with participants of the lymphedema gene network
The associative gene networks reconstructed by ANDSystem were used for the analysis of interactions of the top 10 candidate genes, obtained by criterion 6, with the genes/proteins associated with lymphedema (Additional file 14: Table S14, Fig. 3).
The highest priority, according to Table 1, had the control gene KDR, which encodes Vascular Endothelial Growth Factor Receptor 2. KDR is one of the main regulators of endothelial cell growth , it provides endothelial survival , proliferation, migration , sprouting  and tubular morphogenesis [50, 51]. It is known that a number of polymorphisms in this gene are associated with lymphedema with extreme risk estimates OR > 2, p-value< 0.05 . In the mouse tail lymphedema model, using the quantitative real-time polymerase chain reaction, it was shown that the messenger RNA (mRNA) expression of KDR was increased in the group treated by the low-level laser therapy compared to the control . The associative gene network describing the interactions of the KDR gene consisted of 14 nodes (Fig. 3a). In particular, the KDR protein had protein-protein interactions with VEGF-C, VEGF-D, VEGFR-3 and CBL [53,54,55,56]. It is known that mutations in VEGFC, VEGFR-3 and CBL genes are associated with lymphedema [4, 57,58,59,60,61,62,63]. Serum level of VEGF-D was shown to be significantly higher in the group of patients with primary lymphedema compared with controls . It is also known from the literature that the KDR and SHP2 genes have a positive correlation of their expression levels . Mutations in the SHP2 gene are associated with lymphedema . HGF in combination with VEGF-A can activate VEGFR-2 . Mutations in the HGF gene are also associated with lymphedema [68, 69], and increased expression of HGF improves lymphedema [70, 71].
The second place in Table 1 belonged to the Tumor Necrosis Factor alpha (TNF) gene, which encodes cytokine involved in systemic inflammation . Figure 3B shows the 20 nodes associated with this gene, in particular, TNF was associated with IL-6. It is known from the literature that this protein is capable to activate the expression of the IL-6 gene , increase stability and secretion of IL-6 [74, 75]. In turn, IL-6 has been shown to be increased in models of lymphedema . Another example is the HMGB1 protein, which was also connected to TNF in the associative gene network. It is known that HMGB1 can induce expression of TNF . Also studies show the 2.4-fold increased level of HMGB1 in secondary lymphedema . Besides, TNF was found to be associated with adiponectin in the network, a level of which can be reduced in response to the TNF . Shimizu et al., 2013 showed that adiponectin can promote a lymphatic vessel formation resulting in amelioration of lymphedema .
The potential role of TNF in the molecular mechanisms of lymphedema is discussed in the literature. It is known that TNF can induce endothelial cells apoptosis [81,82,83]. Földi et al., 2000 showed a decrease in the level of expression of TNF after the complex decongestive physiotherapy in patients with peripheral leg lymphedema . Ji, 2007, reports the importance of the role of TNF in the pathophysiological changes of the lymphatic endothelial cells (LECs) and inflammatory lymphangiogenesis . Anuradha et al., 2012, shows the association of inflammatory cytokines (IL-1β, IL-12, and TNF) with pathogenesis of lymphatic filarial infection, which can lead to lymphedema . Jeong et al., 2013 by using the western blot analysis, showed that hyaluronidase treatment of lymphedema in a mouse tail model declined expression of TNF .
The third place in Table 1 had the TEK gene – a receptor tyrosine kinase, involved in a signaling pathway related to the embryonic vascular development . It is known that TEK can reduce proliferation of endothelial cell and apoptosis , while mutations in this gene can cause autosomal dominantly inherited forms of venous malformations [90, 91]. The reconstructed gene network contained four interactions of the TEK gene with genes/proteins (FOXC2, NEMO, Shp2 and CBL) associated with lymphedema (Fig. 3c). From the literature it is known that the loss of FOXC2 can result in TEK downregulation  meanwhile, a complete loss or a significant gain of FOXC2 function can lead to perturbation of lymphatic vessel formation and lymphedema . According to information from the Innatedb , HPRD  and BioGrid  databases, the TEK protein can interact with the NEMO, Shp2 and CBL proteins, mutations in which are associated with lymphedema [63, 66, 97].
The bone morphogenetic protein receptor type II (BMPR2) gene – a serine/threonine receptor kinase involved in cell growth and differentiation, osteogenesis and adipogenesis  was on the fourth place. A reduced BMPR2 expression can induce mitochondrial dysfunction of endothelial cells, promoting a pro-inflammatory and pro-apoptotic state [99, 100]. Kim, Kim, 2014 showed that knockdown of Bmpr2a and Bmpr2b result in lymphatic defects in developing zebrafish . In the constructed associative gene network, BMPR2 has four interactions with the VEGFC, IL6, HGF and FOXL1 genes/proteins (Fig. 3d). It is known that silencing of BMPR2 results in down-regulation of VEGFC expression . Also, VEGFC plays an important role in the functioning of lymphatic vessels while mutations in this gene are associated with lymphedema [57,58,59,60]. Soon et al., 2015 showed that mutations in BMPR2 lead to higher levels of IL6 , which level is increased in lymphedema models . The expression of BMPR2 was up-regulated by HGF . The increase of expression of HGF improves lymphedema [71, 105], and mutations in this gene are associated with the disease [68, 69]. According to the HPRD database, the FOXL1 protein interacts with BMPR2, and mutations in the FOXL1 gene are associated with lymphedema .
Serpin Family E Member 1 (SERPINE1), which is the inhibitor of tissue plasminogen activator and urokinase  that regulates fibrinolysis , appeared to be in the fifth place of the table. The SERPINE1 gene mediates the anti-apoptotic effect in the endothelial cell . A number of studies have shown the association of SERPINE1 with the metastasis of tumors in the lymph nodes [110,111,112,113]. In the associative gene network, SERPINE1 interacts with four genes/proteins: IL6, HGF, HMGB1 and APN (Fig. 3e). It is known that SERPINE1 expression decreases with the addition of IL-6  that is increased in lymphedema models . HGF increases SERPINE1 expression , and increase of HGF expression improves lymphedema [70, 71]. It was shown that levels of HMGB1 and SERPINE1 had positive correlation , while correlation of APN and SERPINE1 was negative . The administration of APN improved the edema of injured tails in the mouse model of lymphedema .
The sixth in the Table 1 was an anti-inflammatory cytokine IL10, playing an important role in the immunoregulation and inflammation . It was shown that IL-10 significantly blocked endothelial apoptosis . However, other study showed that IL-10 had the capacity to induce macrophage apoptosis . An increased gene expression of IL10 was found in keratinocytes derived from limb affected by lymphedema  and in wounded lymphedematous skin . In the associative network IL-10 appeared to be connected with the APN and HGF genes (Fig. 3f). It is known that APN can promote amelioration of lymphedema  and is able to significantly increase IL-10 gene expression and protein secretion . HGF can increase plasma IL-10 concentration [124, 125].
The seventh line in the Table 1 was taken by the CD40LG gene, which is expressed in T cells and its product is presented on their surface . CD40LG can bind CD40 on the B cell surface and is involved in T cell proliferation and cytokine production . It is known from the literature that the product of the CD40LG gene can increased apoptosis of endothelial cells . This gene was found to be interacting only with IL-6 in the associative gene network (Fig. 3g). Sommer et al., 2009 showed that CD40LG can up-regulate expression of IL-6 .
On the eighth line was the CCL2 (C-C Motif Chemokine Ligand 2) gene, which encodes the cytokine possessing chemotactic activity for monocytes  and basophils . Down-regulation of CCL2 by miR-495 resulted in inhibited apoptosis of human umbilical vein endothelial cells . In the associative gene network this gene was interacting with 5 genes/proteins (Fig. 2h). Among the genes involved in interactions with CCL2 were APN, CDC42, HGF, IL6 and TSC2, which are associated with lymphedema [68,69,70,71, 76, 80, 133,134,135]. It is known that APN elevates mRNA and protein level of the CCL2 and stimulates release of CCL2 in primary human monocytes . Ablation of CDC42 induced an overexpression of CCL2 . Müller et al., 2012 showed that CCL2 can induce HGF . It was shown that IL6 induced CCL2  and the IL6 level positively correlated with the CCL2 . Loss of TSC2 function can result in overexpression of CCL2 .
The FASLG gene, located on the ninth row of the Table 1, is a member of tumor necrosis factor superfamily [142, 143]. FASLG is encoding transmembrane protein Fas ligand that is responsible for the induction of apoptosis triggered by binding to FAS [144,145,146]. Abnormal FASLG expression is associated with a metastasis of tumors in the lymph nodes [147,148,149]. In the reconstructed associative gene network FASLG appeared to be associated only with IL6 (Fig. 2i). It is known that FASLG stimulation can enhanced IL-6 release .
On the tenth line of the Table 1 was the ABL1 protooncogene, which is encoding a tyrosine kinase protein, participating in cell division and differentiation, adhesion, and response to oxidative stress . ABL1 is involved in mediating endothelial apoptosis . In reconstructed gene network ABL1 was associated with 7 genes/proteins (Fig. 3j): CBL, CDK5, IL6, SHIP2, Shp2, SOS1 and SOS2. It is known that mutations in these genes are associated with lymphedema [63, 68, 153,154,155,156,157]. The Innatedb, HPRD and BioGrid databases contain data describing protein-protein interactions of ABL1 with CBL, CDK5, SHIP2, Shp2, SOS1 and SOS2. Also, ABL1 positively controls IL-6 expression .
In this work, we performed a search for new potential participants of lymphedema molecular-genetic mechanisms based on the prioritization of genes involved in endothelial cell apoptosis. Six criteria, including the well-known ToppGene and Endeavour methods, as well as our original approaches [19, 20] were used. The use of original methods allowed taking into account the overrepresented Gene Ontology biological processes, structural features of the associative gene network of lymphedema and endothelial apoptosis, as well as diseases comorbid to lymphedema. Verification of the prioritization quality using three different criteria showed significant enrichment of the most priority genes with known genes, which have simultaneous interactions with lymphedema and endothelial apoptosis, as well as with genes differentially expressed in the murine model of lymphedema. Among genes, involved in endothelial cell apoptosis, TNF, TEK, BMPR2, SERPINE1, IL10, CD40LG, CCL2, FASLG and ABL1 were identified as the most promising candidates that can be used for planning the experiments concerning their possible associations with lymphedema. Besides, analysis of the function of these genes can help in understanding the molecular-genetic role of endothelial cell apoptosis in lymphedema.
False discovery rate
- GO PB:
Gene Ontology biological processes
Lymphatic endothelial cells
Farinola N, Piller NB. CYP2A6 polymorphisms: is there a role for pharmacogenomics in preventing coumarin-induced hepatotoxicity in lymphedema patients? Pharmacogenomics. 2007;8(12):151–8.
Poveshchenko AF, Nimaev VV, Lubarsky MS, Konenkov VI. Medical and genetical aspects of lymphedema. Med Genet. 2010;9(9):3–9.
Szuba A, Rockson SG. Lymphedema: classification, diagnosis and therapy. Vasc Med. 1998;3(2):145–56.
Mendola A, Schlögel MJ, Ghalamkarpour A, Irrthum A, Nguyen HL, Fastré E, et al. Mutations in the VEGFR3 signaling pathway explain 36% of familial lymphedema. Mol Syndromol. 2013;4(6):257–66.
Brouillard P, Boon L, Vikkula M. Genetics of lymphatic anomalies. J Clin Invest. 2014;124(3):898–904.
Gupta A, Moore JA. Lymphedema. JAMA Oncol. 2018;4(5):755.
Leung G, Baggott C, West C, Elboim C, Paul SM, Cooper BA, et al. Cytokine candidate genes predict the development of secondary lymphedema following breast cancer surgery. Lymphat Res Biol. 2014;12(1):10–22.
Miaskowski C, Dodd M, Paul SM, West C, Hamolsky D, Abrams G, et al. Lymphatic and angiogenic candidate genes predict the development of secondary lymphedema following breast cancer surgery. PLoS One. 2013;8(4):e60164.
Shaitelman SF, Cromwell KD, Rasmussen JC, Stout NL, Armer JM, Lasinski BB, et al. Recent progress in cancer-related lymphedema treatment and prevention. CA Cancer J Clin. 2015;65(1):55.
Olszewski WL. Pathophysiological aspects of lymphedema of human limbs: I. Lymph protein composition. Lymphat Res Biol. 2003;1(3):235–43.
Tian W, Rockson SG, Jiang X, Kim J, Begaye A, Shuffle EM, Tu AB, Cribb M, Nepiyushchikh Z, Feroze AH, Zamanian RT. Leukotriene B4 antagonism ameliorates experimental lymphedema. Sci Transl Med. 2017;9(389):eaal3920.
Kerr J, Wyllie A, Currie A. Apoptosis: a basic biological phenomenon with wide-ranging implications in tissue kinetics. Br J Cancer. 1972;26:239–57.
Wang JF, Zhang X, Groopman JE. Activation of vascular endothelial growth factor receptor-3 and its downstream signaling promote cell survival under oxidative stress. J Biol Chem. 2004;279(26):27088–97.
Ogunbiyi S, Chinien G, Field D, Humphries J, Burand K, Sawyer B, et al. Smith, for the London lymphedema consortium a. Molecular characterization of dermal lymphatic endothelial cells from primary lymphedema skin. Lymphat Res Biol. 2011;9(1):19–30.
Chen J, Bardes EE, Aronow BJ, Jegga AG. ToppGene Suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res. 2009;37(suppl_2):W305–11.
Moreau Y, Tranchevent LC. Computational tools for prioritizing candidate genes: boosting disease gene discovery. Nat Rev Genet. 2012;13(8):523.
Guney E, Oliva B. Exploiting protein-protein interaction networks for genome-wide disease-gene prioritization. PLoS One. 2012;7(9):e43557.
Tranchevent LC, Ardeshirdavani A, ElShal S, Alcaide D, Aerts J, Auboeuf D, et al. Candidate gene prioritization with Endeavour. Nucleic Acids Res. 2016;44(W1):W117–21.
Saik OV, Demenkov PS, Ivanisenko TV, Bragina EY, Freidin MB, Goncharova IA, et al. Novel candidate genes important for asthma and hypertension comorbidity revealed from associative gene networks. BMC Med Genet. 2018;11(1):15.
Yankina MA, Saik OV, Demenkov PS, Khusnutdinova EK, Rogaev EI, Lavrik IN, et al. Analysis of the interactions of neuronal apoptosis genes in the associative gene network of Parkinson's disease. Vavilovskii Zhurnal Genetiki I SelektsII. 2018;22(1):153–60.
Demenkov PS, Ivanisenko TV, Kolchanov NA, Ivanisenko VA. ANDVisio: a new tool for graphic visualization and analysis of literature mined associative gene networks in the ANDSystem. In Silico Biol. 2012;11(3, 4):149–61.
Ivanisenko VA, Saik OV, Ivanisenko NV, Tiys ES, Ivanisenko TV, Demenkov PS, et al. ANDSystem: an associative network discovery system for automated literature mining in the field of biology. BMC Syst Biol. 2015;9(2):S2.
Glotov AS, Tiys ES, Vashukova ES, Pakin VS, Demenkov PS, Saik OV, et al. Molecular association of pathogenetic contributors to pre-eclampsia (pre-eclampsia associome). BMC Syst Biol. 2015;9(2):S4.
Bragina EY, Tiys ES, Freidin MB, Koneva LA, Demenkov PS, Ivanisenko VA, et al. Insights into pathophysiology of dystropy through the analysis of gene networks: an example of bronchial asthma and tuberculosis. Immunogenetics. 2014;66(7–8):457–65.
Mattingly CJ, Rosenstein MC, Colby GT, Forrest JN Jr, Boyer JL. The comparative Toxicogenomics database (CTD): a resource for comparative toxicological studies. J Exp Zool A Comp Exp Biol. 2006;305(9):689–92.
Rappaport N, Twik M, Plaschkes I, Nudel R, Iny Stein T, Levitt J, et al. MalaCards: an amalgamated human disease compendium with diverse clinical and genetic annotation and structured search. Nucleic Acids Res. 2016;45(D1):D877–87.
Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2016;45(D1):D353–61.
Köhler S, Vasilevsky NA, Engelstad M, Foster E, McMurry J, Aymé S, et al. The human phenotype ontology in 2017. Nucleic Acids Res. 2016;45(D1):D865–76.
Piñero J, Bravo À, Queralt-Rosinach N, Gutiérrez-Sacristán A, Deu-Pons J, Centeno E, García-García J, Sanz F, Furlong LI. DisGeNET: a comprehensive platform integrating information on human disease-associated genes and variants. Nucleic Acids Res. 2016;45(D1):D833-D839.
Usmonov DB, Saik OV, Nimaev VV. The possibilities of bioinformatic analysis in the study of the pathogenesis of lymphatic dysplasia: Engineering, Computer and Information Sciences (SIBIRCON). 2017 International multi-conference on 2017 Sep 18. Hoboken: Wiley-IEEE Press; 2017. p. 512–4.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2008;4(1):44.
Hagberg A, Swart P, S Chult D. Exploring network structure, dynamics, and function using NetworkX. Los Alamos: Los Alamos National Lab.(LANL); 2008. p. 1–6.
Oliphant TE. Python for scientific computing. Comput Sci Eng. 2007;9(3):9–12.
Carbon S, Ireland A, Mungall CJ, Shu S, Marshall B, Lewis S. AmiGO hub, web presence working group. AmiGO: online access to ontology and annotation data. Bioinformatics. 2008;25(2):288–9.
Feinstein AR. The pre-therapeutic classification of co-morbidity in chronic disease. J Chronic Dis. 1970;23(7):455–68.
Puzyrev VP. Genetic bases of human comorbidity. Russ J Genet. 2015;51(4):408–17.
Spearman C. ‘Footrule’for measuring correlation. Br J Psychol 1904–1920. 1906;2(1):89–108.
Radua J, Mataix-Cols D. Meta-analytic methods for neuroimaging data explained. Biology of mood & anxiety disorders. 2012;2(1):6.
Tabibiazar R, Cheung L, Han J, Swanson J, Beilhack A, An A, et al. Inflammatory manifestations of experimental lymphatic insufficiency. PLoS Med. 2006;3(7):e254.
Saito Y, Nakagami H, Kaneda Y, Morishita R. Lymphedema and therapeutic lymphangiogenesis. Biomed Res Int. 2013;2013(804675):1–6.
Yoon YS, Murayama T, Gravereaux E, Tkebuchava T, Silver M, Curry C, et al. VEGF-C gene therapy augments postnatal lymphangiogenesis and ameliorates secondary lymphedema. J Clin Invest. 2003;111(5):717–25.
Coso S, Zeng Y, Sooraj D, Williams ED. Conserved signaling through vascular endothelial growth (VEGF) receptor family members in murine lymphatic endothelial cells. Exp Cell Res. 2011;317(17):2397–407.
Kazenwadel J, Secker GA, Liu YJ, Rosenfeld JA, Wildin RS, Cuellar-Rodriguez J, et al. Loss-of-function germline GATA2 mutations in patients with MDS/AML or MonoMAC syndrome and primary lymphedema reveal a key role for GATA2 in the lymphatic vasculature. Blood. 2011;119(5):1283–91.
Newman B, Lose F, Kedda MA, Francois M, Ferguson K, Janda M, et al. Possible genetic predisposition to lymphedema after breast cancer. Lymphat Res Biol. 2012;10(1):2–13.
Jenssen TK, Lægreid A, Komorowski J, Hovig E. A literature network of human genes for high-throughput analysis of gene expression. Nat Genet. 2001;28(1):21.
Liu Z, Qi L, Li Y, Zhao X, Sun B. VEGFR2 regulates endothelial differentiation of colon cancer cells. BMC Cancer. 2017;17(1):593.
Ou JM, Yu ZY, Qiu MK, Dai YX, Dong Q, Shen J, et al. Knockdown of VEGFR2 inhibits proliferation and induces apoptosis in hemangioma-derived endothelial cells. Eur J Histochem. 2014;58(1):2263.
Liu Y, Qiao Y, Hu C, Liu L, Zhou L, Liu B, et al. VEGFR2 inhibition by RNA interference affects cell proliferation, migration, invasion, and response to radiation in Calu-1 cells. Clin Transl Oncol. 2016;18(2):212–9.
Gaengel K, Niaudet C, Hagikura K, Laviña B, Muhl L, Hofmann JJ, et al. The sphingosine-1-phosphate receptor S1PR1 restricts sprouting angiogenesis by regulating the interplay between VE-cadherin and VEGFR2. Dev Cell. 2012;23(3):587–99.
van Tuyl M, Groenman F, Wang J, Kuliszewski M, Liu J, Tibboel D, et al. Angiogenic factors stimulate tubular branching morphogenesis of sonic hedgehog-deficient lungs. Dev Biol. 2007;303(2):514–26.
Mellberg S, Dimberg A, Bahram F, Hayashi M, Rennel E, Ameur A, et al. Transcriptional profiling reveals a critical role for tyrosine phosphatase VE-PTP in regulation of VEGFR2 activity and endothelial cell morphogenesis. FASEB J. 2009;23(5):1490–502.
Jang DH, Song DH, Chang EJ, Jeon JY. Anti-inflammatory and lymphangiogenetic effects of low-level laser therapy on lymphedema in an experimental mouse tail model. Lasers Med Sci. 2016;31(2):289–96.
Enholm B, Karpanen T, Jeltsch M, Kubo H, Stenbacz F, Prevo R, et al. Adenoviral expression of vascular endothelial growth factor-C induces lymphangiogenesis in the skin. Circ Res. 2001;88(6):623–9.
Partanen TA, Arola J, Saaristo A, Jussila L, Ora A, Miettinen M, et al. VEGF-C and VEGF-D expression in neuroendocrine cells and their receptor, VEGFR-3, in fenestrated blood vessels in human tissues. FASEB J. 2000;14(13):2087–96.
Achen MG, Roufail S, Domagala T, Catimel B, Nice EC, Geleick DM, et al. Monoclonal antibodies to vascular endothelial growth factor-D block its interactions with both VEGF receptor-2 and VEGF receptor-3. Eur J Biochem. 2000;267(9):2505–15.
Meyer RD, Sacks DB, Rahimi N. IQGAP1-dependent signaling pathway regulates endothelial cell proliferation and angiogenesis. PLoS One. 2008;3(12):e3848.
Shin M, Male I, Beane TJ, Villefranc JA, Kok FO, Zhu LJ, Lawson ND. Vegfc acts through ERK to induce sprouting and differentiation of trunk lymphatic progenitors. Development. 2016;143(20):3785-3795.
Le Guen L, Karpanen T, Schulte D, Harris NC, Koltowska K, Roukens G, Bower NI, Van Impel A, Stacker SA, Achen MG, Schulte-Merker S. Ccbe1 regulates Vegfc-mediated induction of Vegfr3 signaling during embryonic lymphangiogenesis. Development. 2014;141(6):1239-49.
Gousopoulos E, Proulx ST, Bachmann SB, Dieterich LC, Scholl J, Karaman S, et al. An important role of VEGF-C in promoting lymphedema development. J Investig Dermatol. 2017;137(9):1995–2004.
Fastré E, Lanteigne LE, Helaers R, Giacalone G, Revencu N, Dionyssiou D, et al. Splice-site mutations in VEGFC cause loss of function and nonne-Milroy-like primary lymphedema. Clin Genet. 2018;94(1):179–81.
Saaristo A, Veikkola T, Tammela T, Enholm B, Karkkainen MJ, Pajusola K, et al. Lymphangiogenic gene therapy with minimal blood vascular side effects. J Exp Med. 2002;196(6):719–30.
Mizuno S, Yamada Y, Yamada K, Nomura N, Wakamatsu N. Clinical variability in a Japanese hereditary lymphedema type I family with an FLT4 mutation. Congenit Anom. 2005;45(2):59–61.
Hanson HL, Wilson MJ, Short JP, Chioza BA, Crosby AH, Nash RM, et al. Germline CBL mutation associated with a Noonan-like syndrome with primary lymphedema and teratoma associated with acquired uniparental isodisomy of chromosome 11q23. Am J Med Genet A. 2014;164(4):1003–9.
Fink AM, Kaltenegger I, Schneider B, Fruhauf J, Jurecka W, Steiner A. Serum level of VEGF-D in patients with primary lymphedema. Lymphology. 2004;37(4):185–9.
Tang C, Luo D, Yang H, Wang Q, Zhang R, Liu G, et al. Expression of SHP2 and related markers in non–small cell lung cancer: a tissue microarray study of 80 cases. Appl Immunohistochem Mol Morphol. 2013;21(5):386–94.
Yeang CH, Ma GC, Shih JC, Yang YS, Chen CP, Chang SP, et al. Genome-wide gene expression analysis implicates the immune response and lymphangiogenesis in the pathogenesis of fetal chylothorax. PLoS One. 2012;7(4):e34901.
Sulpice E, Ding S, Muscatelli-Groux B, Bergé M, Han ZC, Plouet J, et al. Cross-talk between the VEGF-A and HGF signalling pathways in endothelial cells. Biol Cell. 2009;101(9):525–39.
Finegold DN, Schacht V, Kimak MA, Lawrence EC, Foeldi E, Karlsson JM, et al. HGF and MET mutations in primary and secondary lymphedema. Lymphat Res Biol. 2008;6(2):65–8.
Michelini S, Vettori A, Maltese PE, Cardone M, Bruson A, Fiorentino A, et al. Genetic screening in a large cohort of italian patients affected by primary lymphedema using a next generation sequencing (NGS) approach. Lymphology. 2016;49(2):57–72.
Saito Y, Nakagami H, Morishita R, Takami Y, Kikuchi Y, Hayashi H, et al. Transfection of human hepatocyte growth factor gene ameliorates secondary lymphedema via promotion of lymphangiogenesis. Circulation. 2006;114(11):1177–84.
Lee CY, Kang JY, Lim S, Ham O, Chang W, Jang DH. Hypoxic conditioned medium from mesenchymal stem cells promotes lymphangiogenesis by regulation of mitochondrial-related proteins. Stem Cell Res Ther. 2016;7(1):38.
Popa C, Netea MG, Van Riel PL, Van Der Meer JW, Stalenhoef AF. The role of TNF-α in chronic inflammatory conditions, intermediary metabolism, and cardiovascular risk. J Lipid Res. 2007;48(4):751–62.
Tomita N, Morishita R, Tomita S, Kaneda Y, Higaki J, Ogihara T, et al. Inhibition of TNF-α, induced cytokine and adhesion molecule. Nephron Exp Nephrol. 2001;9(3):181–90.
Shimada M, Andoh A, Hata K, Tasaki K, Araki Y, Fujiyama Y, et al. IL-6 secretion by human pancreatic periacinar myofibroblasts in response to inflammatory mediators. J Immunol. 2002;168(2):861–8.
Lee C, Oh JI, Park J, Choi JH, Bae EK, Lee HJ, et al. TNFα mediated IL-6 secretion is regulated by JAK/STAT pathway but not by MEK phosphorylation and AKT phosphorylation in U266 multiple myeloma cells. Biomed Res Int. 2013;2013:580135.
Cuzzone DA, Weitman ES, Albano NJ, Ghanta S, Savetsky IL, Gardenier JC, et al. IL-6 regulates adipose deposition and homeostasis in lymphedema. Am J Phys Heart Circ Phys. 2014;306(10):H1426–34.
Barkauskaite V, Ek M, Popovic K, Harris HE, Wahren-Herlenius M, Nyberg F. Translocation of the novel cytokine HMGB1 to the cytoplasm and extracellular space coincides with the peak of clinical activity in experimentally UV-induced lesions of cutaneous lupus erythematosus. Lupus. 2007;16(10):794–802.
Zampell JC, Yan A, Avraham T, Andrade V, Malliaris S, Aschen S, et al. Temporal and spatial patterns of endogenous danger signal expression after wound healing and in response to lymphedema. Am J Phys Cell Phys. 2011;300(5):C1107–21.
Wang B, Trayhurn P. Acute and prolonged effects of TNF-α on the expression and secretion of inflammation-related adipokines by human adipocytes differentiated in culture. Pflugers Arch. 2006;452(4):418–27.
Shimizu Y, Shibata R, Ishii M, Ohashi K, Kambara T, Uemura Y, et al. Adiponectin-mediated modulation of lymphatic vessel formation and lymphedema. J Am Heart Assoc. 2013;2(5):e000438.
Mariño E, Cardier JE. Differential effect of IL-18 on endothelial cell apoptosis mediated by TNF-α and Fas (CD95). Cytokine. 2003;22(5):142–8.
Lejeune FJ, Rüegg C. Recombinant human tumor necrosis factor: an efficient agent for cancer treatment. Bull Cancer. 2006;93(8):10090–100.
Acquavella N, Quiroga MF, Wittig O, Cardier JE. Effect of simvastatin on endothelial cell apoptosis mediated by Fas and TNF-α. Cytokine. 2010;49(1):45–50.
Foldi E, Sauerwald A, Hennig B. Effect of complex decongestive physiotherapy on gene expression for the inflammatory response in peripheral lymphedema. Lymphology. 2000;33(1):19–23.
Ji RC. Lymphatic endothelial cells, inflammatory lymphangiogenesis, and prospective players. Curr Med Chem. 2007;14(22):2359–68.
Anuradha R, George PJ, Kumar NP, Fay MP, Kumaraswami V, Nutman TB, et al. Circulating microbial products and acute phase proteins as markers of pathogenesis in lymphatic filarial disease. PLoS Pathog. 2012;8(6):e1002749.
Jeong HJ, Roh K, Kim G, Kim Y, Lee J, Lee M, et al. Hyaluronidase treatment of acute lymphedema in a mouse tail model. Lymphology. 2013;46(4):160–72.
Dumont DJ, Gradwohl G, Fong GH, Puri MC, Gertsenstein M, Auerbach A, et al. Dominant-negative and targeted null mutations in the endothelial receptor tyrosine kinase, tek, reveal a critical role in vasculogenesis of the embryo. Genes Dev. 1994;8(16):1897–909.
Hu HT, Huang YH, Chang YA, Lee CK, Jiang MJ, Wu LW. Tie2-R849W mutant in venous malformations chronically activates a functional STAT1 to modulate gene expression. J Investig Dermatol. 2008;128(9):2325–33.
Frigerio A, Stevenson DA, Grimmer JF. The genetics of vascular anomalies. Curr Opin Otolaryngol Head Neck Surg. 2012;20(6):527–32.
Yadav P, De Castro DK, Waner M, Meyer L, Fay A. Vascular anomalies of the head and neck: a review of genetics. Semin Ophthalmol. 2013;28(5–6):257–66.
Thomson BR, Heinen S, Jeansson M, Ghosh AK, Fatima A, Sung HK, et al. A lymphatic defect causes ocular hypertension and glaucoma in mice. J Clin Invest. 2014;124(10):4320–4.
Tavian D, Missaglia S, Maltese PE, Michelini S, Fiorentino A, Ricci M, et al. FOXC2 disease-mutations identified in lymphedema-distichiasis patients cause both loss and gain of protein function. Oncotarget. 2016;7(34):54228.
Breuer K, Foroushani AK, Laird MR, Chen C, Sribnaia A, Lo R, et al. InnateDB: systems biology of innate immunity and beyond—recent updates and continuing curation. Nucleic Acids Res. 2012;41(D1):D1228–33.
Keshava Prasad TS, Goel R, Kandasamy K, Keerthikumar S, Kumar S, Mathivanan S, Telikicherla D, Raju R, Shafreen B, Venugopal A, Balakrishnan L. Human protein reference database—2009 update. Nucleic Acids Res. 2008;37(suppl_1):D767–72.
Stark C, Breitkreutz BJ, Reguly T, Boucher L, Breitkreutz A, Tyers M. BioGRID: a general repository for interaction datasets. Nucleic Acids Res. 2006;34(suppl_1):D535–9.
Roberts CM, Angus JE, Leach IH, McDermott EM, Walker DA, Ravenscroft JC. A novel NEMO gene mutation causing osteopetrosis, lymphoedema, hypohidrotic ectodermal dysplasia and immunodeficiency (OL-HED-ID). Eur J Pediatr. 2010;169(11):1403–7.
Xiao YT, Xiang LX, Shao JZ. Bone morphogenetic protein. Biochem Biophys Res Commun. 2007;362(3):550–3.
Spiekerkoetter E, Tian X, Cai J, Hopper RK, Sudheendra D, Li CG, et al. FK506 activates BMPR2, rescues endothelial dysfunction, and reverses pulmonary hypertension. J Clin Invest. 2013;123(8):3600–13.
Diebold I, Hennigs JK, Miyagawa K, Li CG, Nickel NP, Kaschwich M, et al. BMPR2 preserves mitochondrial function and DNA during reoxygenation to promote endothelial cell survival and reverse pulmonary hypertension. Cell Metab. 2015;21(4):596–608.
Kim JD, Kim J. Alk3/Alk3b and Smad5 mediate BMP signaling during lymphatic development in zebrafish. Mol Cells. 2014;37(3):270.
Zeng P, Cai S, Zhang JN, Yi FM, Jiang WM, Wu JB. Effects of siRNA targeting BMPR-II on the biological activities of human liver cancer cells and its mechanism. Cancer Cell Int. 2014;14(1):55.
Soon E, Crosby A, Southwood M, Yang P, Tajsic T, Toshner M, et al. Bone morphogenetic protein receptor type II deficiency and increased inflammatory cytokine production. A gateway to pulmonary arterial hypertension. Am J Respir Crit Care Med. 2015;192(7):859–72.
Ye L, Lewis-Russell JM, Davies G, Sanders AJ, Kynaston H, Jiang WG. Hepatocyte growth factor up-regulates the expression of the bone morphogenetic protein (BMP) receptors, BMPR-IB and BMPR-II, in human prostate cancer cells. Int J Oncol. 2007;30(2):521–9.
Saito Y, Nakagami H, Morishita R. Transfection of human hepatocyte growth factor gene ameliorates secondary lymphedema via promotion of lymphangiogenesis. J Vasc Surg. 2007;45(3):637.
Butler MG, Dagenais SL, Garcia-Perez JL, Brouillard P, Vikkula M, Strouse P, et al. Microcephaly, intellectual impairment, bilateral vesicoureteral reflux, distichiasis, and glomuvenous malformations associated with a 16q24. 3 contiguous gene deletion and a Glomulin mutation. Am J Med Genet A. 2012;158(4):839–49.
Ye Y, Vattai A, Zhang X, Zhu J, Thaler CJ, Mahner S, et al. Role of plasminogen activator inhibitor type 1 in pathologies of female reproductive diseases. Int J Mol Sci. 2017;18(8):1651.
Chen R, Yan J, Liu P, Wang Z, Wang C. Plasminogen activator inhibitor links obesity and thrombotic cerebrovascular diseases: the roles of PAI-1 and obesity on stroke. Metab Brain Dis. 2017;32(3):667–73.
Yao H, He G, Chen C, Yan S, Lu L, Song L, et al. PAI1: a novel PP1-interacting protein that mediates human plasma's anti-apoptotic effect in endothelial cells. J Cell Mol Med. 2017;21(9):2068–76.
Eljuga D, Razumovic JJ, Bulic K, Petrovecki M, Draca N, Bulic SO. Prognostic importance of PAI-1 in node negative breast cancer patients—results after 10 years of follow up. Pathol Res Pract. 2011;207(5):290–4.
Dhanda J, Triantafyllou A, Liloglou T, Kalirai H, Lloyd B, Hanlon R, et al. SERPINE1 and SMA expression at the invasive front predict extracapsular spread and survival in oral squamous cell carcinoma. Br J Cancer. 2014;111(11):2114.
Zhang Y, Zhang YL, Chen HM, Pu HW, Ma WJ, Li XM, et al. Expression of Bmi-1 and PAI-1 in esophageal squamous cell carcinoma. World J Gastroenterol: WJG. 2014;20(18):5533.
Sang Y, Chen MY, Luo D, Zhang RH, Wang L, Li M, et al. TEL2 suppresses metastasis by down-regulating SERPINE1 in nasopharyngeal carcinoma. Oncotarget. 2015;6(30):29240.
Umemura K, Ishioka SI, Endo T, Ezaka Y, Takahashi M, Saito T. Roles of microRNA-34a in the pathogenesis of placenta accreta. J Obstet Gynaecol Res. 2013;39(1):67–74.
Imagawa S, Fujii S, Dong J, Furumoto T, Kaneko T, Zaman T, et al. Hepatocyte growth factor regulates E box–dependent plasminogen activator inhibitor type 1 gene expression in HepG2 liver cells. Arterioscler Thromb Vasc Biol. 2006;26(10):2407–13.
Nomura S, Maeda Y, Ishii K, Katayama Y, Yagi H, Fujishima N, et al. Relationship between HMGB1 and PAI-1 after allogeneic hematopoietic stem cell transplantation. J Blood Med. 2016;7:1–4.
Nomura S, Taniura T, Shouzu A, Omoto S, Inami N, Fujita S, et al. Effects of pitavastatin on plasminogen activator inhibitor-1 in hyperlipidemic patients. Int J General Med. 2012;5:535.
Hedrich CM, Bream JH. Cell type-specific regulation of IL-10 expression in inflammation and disease. Immunol Res. 2010;47(1–3):185–206.
Yin Y, Liu W, Ji G, Dai Y. The essential role of p38 MAPK in mediating the interplay of oxLDL and IL-10 in regulating endothelial cell apoptosis. Eur J Cell Biol. 2013;92(4–5):150–9.
Wang ZQ, Bapat AS, Rayanade RJ, Dagtas AS, Hoffmann MK. Interleukin-10 induces macrophage apoptosis and expression of CD16 (FcγRIII) whose engagement blocks the cell death programme and facilitates differentiation. Immunology. 2001;102(3):331–7.
Baroni A, Buommino E, Piccolo V, Chessa MA, Russo T, Cozza V, et al. Alterations of skin innate immunity in lymphedematous limbs: correlations with opportunistic diseases. Clin Dermatol. 2014;32(5):592–8.
Kimura T, Sugaya M, Blauvelt A, Okochi H, Sato S. Delayed wound healing due to increased interleukin-10 expression in mice with lymphatic dysfunction. J Leukoc Biol. 2013;94(1):137–45.
Kumada M, Kihara S, Ouchi N, Kobayashi H, Okamoto Y, Ohashi K, et al. Adiponectin specifically increased tissue inhibitor of metalloproteinase-1 through interleukin-10 expression in human macrophages. Circulation. 2004;109(17):2046–9.
Warzecha Z, Dembiński A, Ceranowicz P, Konturek SJ, Tomaszewska R, Stachura J, et al. Inhibition of cyclooxygenase-2 reduces the protective effect of hepatocyte growth factor in experimental pancreatitis. Eur J Pharmacol. 2004;486(1):107–19.
Chen PM, Liu KJ, Hsu PJ, Wei CF, Bai CH, Ho LJ, et al. Induction of immunomodulatory monocytes by human mesenchymal stem cell-derived hepatocyte growth factor through ERK1/2. J Leukoc Biol. 2014;96(2):295–303.
Mikolajczak SA, Ma BY, Yoshida T, Yoshida R, Kelvin DJ, Ochi A. The modulation of CD40 ligand signaling by transmembrane CD28 splice variant in human T cells. J Exp Med. 2004;199(7):1025–31.
Blotta MH, Marshall JD, DeKruyff RH, Umetsu DT. Cross-linking of the CD40 ligand on human CD4+ T lymphocytes generates a costimulatory signal that up-regulates IL-4 synthesis. J Immunol. 1996;156(9):3133–40.
Wu CF, Huang FD, Sui RF, Sun JX. Preeclampsia serum upregulates CD40/CD40L expression and induces apoptosis in human umbilical cord endothelial cells. Reprod Biol Endocrinol. 2012;10(1):28.
Sommer S, Pudrith CB, Colvin CJ, Coussens PM. Mycobacterium avium subspecies paratuberculosis suppresses expression of IL-12p40 and iNOS genes induced by signalling through CD40 in bovine monocyte-derived macrophages. Vet Immunol Immunopathol. 2009;128(1–3):44–52.
Ogilvie P, Paoletti S, Clark-Lewis I, Uguccioni M. Eotaxin-3 is a natural antagonist for CCR2 and exerts a repulsive effect on human monocytes. Blood. 2003;102(3):789–94.
Kuna P, Reddigari SR, Rucinski D, Oppenheim JJ, Kaplan AP. Monocyte chemotactic and activating factor is a potent histamine-releasing factor for human basophils. J Exp Med. 1992;175(2):489–93.
Liu D, Zhang XL, Yan CH, Li Y, Tian XX, Zhu N, et al. MicroRNA-495 regulates the proliferation and apoptosis of human umbilical vein endothelial cells by targeting chemokine CCL2. Thromb Res. 2015;135(1):146–54.
Takenouchi T, Okamoto N, Ida S, Uehara T, Kosaki K. Further evidence of a mutation in CDC42 as a cause of a recognizable syndromic form of thrombocytopenia. Am J Med Genet A. 2016;170(4):852–5.
Kacerovska D, Kerl K, Michal M, Filipova H, Vrtel R, Vanecek T, et al. Giant angiofibromas in tuberous sclerosis complex: a possible role for localized lymphedema in their pathogenesis. J Am Acad Dermatol. 2012;67(6):1319–26.
Navarre P, Poitras B. Lymphoedema in tuberous sclerosis: case report and review of the literature. J Pediatr Orthop. 2014;34(6):e27–32.
Neumeier M, Bauer S, Brühl H, Eisinger K, Kopp A, Abke S, et al. Adiponectin stimulates release of CCL2, −3, −4 and −5 while the surface abundance of CCR2 and −5 is simultaneously reduced in primary human monocytes. Cytokine. 2011;56(3):573–80.
Deroanne CF, Hamelryckx D, Ho TG, Lambert CA, Catroux P, Lapière CM, et al. Cdc42 downregulates MMP-1 expression by inhibiting the ERK1/2 pathway. J Cell Sci. 2005;118(6):1173–83.
Müller AM, Jun E, Conlon H, Sadiq SA. Cerebrospinal hepatocyte growth factor levels correlate negatively with disease activity in multiple sclerosis. J Neuroimmunol. 2012;251(1–2):80–6.
Thomas M, Bayha C, Klein K, Müller S, Weiss TS, Schwab M, et al. The truncated splice variant of peroxisome proliferator-activated receptor alpha, PPARα-tr, autonomously regulates proliferative and pro-inflammatory genes. BMC Cancer. 2015;15(1):488.
Cheung CM, Vania M, Ang M, Chee SP, Li J. Comparison of aqueous humor cytokine and chemokine levels in diabetic patients with and without retinopathy. Mol Vis. 2012;18:830.
Li S, Takeuchi F, Wang JA, Fuller C, Pacheco-Rodriguez G, Moss J, et al. MCP-1 overexpressed in tuberous sclerosis lesions acts as a paracrine factor for tumor development. J Exp Med. 2005;202(5):617–24.
Ruiz-García R, Mora S, Lozano-Sánchez G, Martínez-Lostao L, Paz-Artal E, Ruiz-Contreras J, Anel A, González-Granado LI, Moreno D, Allende LM. Decreased activation-induced cell death by EBV-transformed B cells from a patient with autoimmune lymphoproliferative syndrome caused by a novel FASLG mutation. Pediatr Res. 2015;78(6):603-8.
Schneider P, Bodmer JL, Holler N, Mattmann C, Scuderi P, Terskikh A, et al. Characterization of Fas (Apo-1, CD95)-fas ligand interaction. J Biol Chem. 1997;272(30):18827–33.
Kokkonen TS, Augustin MT, Mäkinen JM, Kokkonen J, Karttunen TJ. High endothelial venules of the lymph nodes express Fas ligand. J Histochem Cytochem. 2004;52(5):693–9.
Verma RK, Gunda V, Pawar SC, Sudhakar YA. Extra cellular matrix derived metabolite regulates angiogenesis by FasL mediated apoptosis. PLoS One. 2013;8(12):e80555.
Gong Q, Qiu S, Li S, Ma Y, Chen M, Yao Y, et al. Proapoptotic PEDF functional peptides inhibit prostate tumor growth—a mechanistic study. Biochem Pharmacol. 2014;92(3):425–37.
Sjöström-Mattson J, Von Boguslawski K, Bengtsson NO, Mjaaland I, Salmenkivi K, Blomqvist C. The expression of p53, bcl-2, bax, fas and fasL in the primary tumour and lymph node metastases of breast cancer. Acta Oncol. 2009;48(8):1137–43.
Li Q, Peng J, Li XH, Liu T, Liang QC, Zhang GY. Clinical significance of Fas and FasL protein expression in gastric carcinoma and local lymph node tissues. World J Gastroenterol: WJG. 2010;16(10):1274.
Knechtel G, Hofmann G, Gerger A, Renner W, Langsenlehner T, Szkandera J, et al. Analysis of common germline polymorphisms as prognostic factors in patients with lymph node-positive breast cancer. J Cancer Res Clin Oncol. 2010;136(12):1813–9.
Seitz DH, Palmer A, Niesler U, Braumüller ST, Bauknecht S, Gebhard F, et al. Altered expression of Fas receptor on alveolar macrophages and inflammatory effects of soluble Fas ligand following blunt chest trauma. Shock. 2011;35(6):610–7.
Cobbaut M, Derua R, Döppler H, Lou HJ, Vandoninck S, Storz P, et al. Differential regulation of PKD isoforms in oxidative stress conditions through phosphorylation of a conserved Tyr in the P+ 1 loop. Sci Rep. 2017;7(1):887.
Xu J, Millard M, Ren X, Cox OT, Erdreich-Epstein A. c-Abl mediates endothelial apoptosis induced by inhibition of integrins αvβ3 and αvβ5 and by disruption of actin. Blood. 2010;115(13):2709-18.
Liebl J. Cdk5 and Foxc2–a new relationship in the lymphatic vasculature. Oncotarget. 2015;6(26):21799.
Fu MR, Conley YP, Axelrod D, Guth AA, Yu G, Fletcher J, et al. Precision assessment of heterogeneity of lymphedema phenotype, genotypes and risk prediction. Breast. 2016;29:231–40.
Agollah GD, Gonzalez-Garay ML, Rasmussen JC, Tan IC, Aldrich MB, Darne C, et al. Evidence for SH2 domain-containing 5′-inositol phosphatase-2 (SHIP2) contributing to a lymphatic dysfunction. PLoS One. 2014;9(11):e112548.
Smpokou P, Tworog-Dube E, Kucherlapati RS, Roberts AE. Medical complications, clinical findings, and educational outcomes in adults with Noonan syndrome. Am J Med Genet A. 2012;158(12):3106–11.
Cordeddu V, Yin JC, Gunnarsson C, Virtanen C, Drunat S, Lepri F, et al. Activating mutations affecting the Dbl homology domain of SOS2 cause Noonan syndrome. Hum Mutat. 2015;36(11):1080–7.
Reynaud D, Pietras E, Barry-Holson K, Mir A, Binnewies M, Jeanne M, et al. IL-6 controls leukemic multipotent progenitor cell fate and contributes to chronic myelogenous leukemia development. Cancer Cell. 2011;20(5):661–73.
Reconstruction and analysis of the endothelial cell apoptotic gene network was performed with the support of the Russian Science Foundation grant «Programmed cell death induced via death receptors: Delineating molecular mechanisms of apoptosis initiation via molecular modeling» 14-44-00011.
Reconstruction and analysis of the lymphedema gene network was supported by State Budgeted Project № 0324-2019-0045 «Investigation of cellular and molecular biological mechanisms, the identification of biomarkers, the development of new approaches to the diagnosis and correction of vascular remodeling in the pathology of the lymphatic system, diabetes mellitus and malignant neoplasms».
Reconstruction of gene networks describing the interactions of the top 10 of the most promising candidate genes with lymphedema was performed with the support of the State Budgeted Project № 0324-2019-0040 «Genetic basis of biotechnology and bioinformatics».
The publication cost was covered by the Russian Science Foundation grant «Programmed cell death induced via death receptors: Delineating molecular mechanisms of apoptosis initiation via molecular modeling» 14–44-00011.
Availability of data and materials
Results are shared in the additional files.
About this supplement
This article has been published as part of BMC Medical Genomics Volume 12 Supplement 2, 2019: Selected articles from BGRS\SB-2018: medical genomics. The full contents of the supplement are available online at https://bmcmedgenomics.biomedcentral.com/articles/supplements/volume-12-supplement-2.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Lists of genes associated with lymphedema. (XLSX 12 kb)
Table S2. Results of gene ontology enrichment analysis. (XLSX 14 kb)
Table S3. Lists of diseases comorbid to lymphedema. (XLSX 11 kb)
Table S4. Values of betweenness centrality of nodes of lymphedema associative gene network. (XLSX 16 kb)
Table S5. Ranks of analyzed genes. (XLSX 23 kb)
Table S6. Results of prioritization according to criterion 1. (XLSX 30 kb)
Table S7. Results of prioritization according to criterion 2. (XLSX 33 kb)
Table S8. Results of prioritization according to criterion 3. (XLSX 17 kb)
Table S9. Results of prioritization according to criterion 4A. (XLSX 15 kb)
Table S10. Results of prioritization according to criterion 4B. (XLSX 16 kb)
Table S11. Results of prioritization according to criterion 5. (XLSX 17 kb)
Table S12. Frequency of mentioning analyzed gene names together with lymphedema in in full-text articles. (XLSX 13 kb)
Table S13. Differentially expressed genes in mouse model of lymphedema. (XLSX 471 kb)
Table S14. Interactions of top 10 most promising candidate genes with genes/proteins associated with lymphedema. (XLSX 173 kb)