Comprehensive analysis reveals potential hub genes and therapeutic drugs in an acquired lymphedema model.

Background: Acquired lymphedema is a common and often severe complication of breast cancer surgery and radiology that seriously affects patients' quality of life. Nevertheless, the pathogenesis for acquired lymphedema is complex and remains unclear. The aim of this study is to find out possible genetic markers and potential drugs for acquired lymphedema.
Methods: First, the GSE4333 datasets, which include expression data for six female humanized hairless immunocompetent SKH-1 mice (the condition of whom mimics acquired lymphedema), were reanalyzed. According to the criteria of a fold change (FC) ≥1.4 and an adjusted P value <0.05, we identified the differentially expressed genes (DEGs) between a normal group and the lymphedema group. Next, we analyzed the Gene Ontology (GO) terms and enriched signaling pathways associated with these DEGs with an online tool DAVID. We also constructed protein-protein interaction (PPI) networks and selected meaningful gene modules for additional gene-drug interaction research. Finally, the extant drugs targeting these module genes were identified for further study of their therapeutic effects against acquired lymphedema.
Results: A total of 481 DEGs were identified that were closely associated with the immune system, inflammatory response, and extracellular matrix (ECM) structural constituent terms, among others. Moreover, we identified the top 10 significant genes in the PPI networks and identified one extant drug, fiboflapon, that targets the ALOX5AP gene. Conclusions: We ultimately identified 10 hub genes, molecular mechanisms, and one extant drug related to acquired lymphedema. The findings identified targets and a potential drug for further research on acquired lymphedema. 2022 Gland Surgery. All rights reserved.

Introduction

In 2020, female breast neoplasms exceeded lung cancer as the primary cancer type worldwide in terms of incidence (1). Axillary lymph node dissection (ALND) is performed when lymph node metastasis of breast cancer is confirmed. Due to devastation of the axillary lymphatic structure, upper limb lymphedema affects approximately 30–50% of ALND patients (2). Acquired lymphedema is a common and often severe complication of breast cancer surgery and radiology that seriously affects patients’ quality of life. In addition, many patients experience a series of physical and psychological symptoms (3-5).

There are many risk factors for acquired lymphedema, including overweight or obesity, older age, axillary radiotherapy, infection or trauma of the ipsilateral upper extremity, a broad scope of axillary surgery, and adjuvant chemotherapy. However, the conventional view of lymphatic obstruction does not adequately explain the generation of lymphedema. Presently, there are 3 main hypotheses regarding lymphedema pathogenesis (6): the lymphatic failure hypothesis, the hemodynamic hypothesis, and the interstitial hypothesis. Most researchers believe that acquired lymphedema is caused by the obstruction or interruption of lymphatic reflux of the upper limb due to axillary surgery and radiotherapy (7). A large amount of protein-rich lymph liquid is retained in the tissue space to form high protein edema. Nevertheless, the pathogenesis for lymphedema is complex and remains unclear.

Currently, acquired lymphedema diagnoses commonly rely on history, physical findings, and even imaging studies (8). Lymphedema therapy mainly involves complete decongestive physiotherapy (multilayer low-stretch bandage application, manual lymph drainage, skin care, exercises), surgical treatment, and other approaches. However, early diagnosis of lymphedema can be difficult, and lymphedema treatment efficacy is not ideal. The main reason for these challenges is a lack of clarity regarding the potential molecular mechanism of acquired lymphedema. In addition, research focusing on the genetics of acquired lymphedema pathogenesis has been relatively scarce.

To further study the genetic pathogenesis of acquired lymphoedema, we performed comprehensive bioinformatics analyses to determine the hub genes and signaling pathways in samples from a model imitating human acquired lymphedema. The treatment of lymphedema is mainly non-surgical method. Physical therapy is the primary and effective treatment. It is difficult and unethical to obtain human tissue samples from acquired lymphedema patients. So mouse lymphedema models were applied for analysis. Initially, the GSE4333 gene expression datasets (9) were downloaded from the National Center for Biotechnology Information (NCBI) databases (10). These datasets included data for female hairless immunocompetent (the capacity for a normal immune response) SKH-1 mice, which represent a model of human acquired lymphedema. Subsequently, we identified dysregulated genes between different groups. In addition, we used several methods to perform Gene Ontology (GO), signaling pathway enrichment annotation, and protein-protein interaction (PPI) analyses (11,12). Ultimately, the potential hub genes, related pathways, and existing candidate drugs were identified. The findings might provide insights into the molecular mechanisms underlying acquired lymphedema and provide theoretical support for clinical treatment. We present the following article in accordance with the STREGA reporting checklist (available at https://gs.amegroups.com/article/view/10.21037/gs-22-453/rc).

Methods

Microarray data information

We downloaded the GSE4333 expression profile (.txt format files) and correlated information built in the GPL3506 platform (.soft format file) from the NCBI-Gene Expression Omnibus (GEO) website (10,13). The GSE4333 datasets contain 3 arrays of tail skin from normal mice (no intervention) and 3 arrays of tail skin from lymphedema mice (with surgical lymphatic vessel blockage). We analyzed 6 samples from the GSE4333 datasets to identify differentially expressed genes (DEGs) between the normal and lymphedema groups.

Data preprocessing

First, the messenger RNA (mRNA) probes were retained, and the other non-mRNA probes were discarded. According to the information of the GPL3506 platform, we replaced the probe identification numbers with the official gene symbols. Finally, we used Python to explore the gene expression matrix and obtain DEGs between the normal group and the lymphedema group (14). The fold change (FC) cut-off was ≥1.4, and the adjusted P value cut-off was <0.05.

GO and pathway enrichment analysis

The GO analysis method is a common and helpful method for interpreting gene products and functional characteristics (11). The GO analysis terms were classified into biological process (BP), cellular component (CC), and molecular function (MF) categories. The Kyoto Encyclopedia of Genes and Genomes (KEGG) database is used to interpret biological functions and features of organic systems (12). We used the Database for Annotation, Visualization and Integrated Discovery (DAVID) to determine the GO and KEGG enrichment of the DEGs (15). The DEGs were used as search input for GO and KEGG pathway analyses using the official gene symbol and species information (Mus musculus).

PPI analysis and module determination

The online database Search Tool for the Retrieval of Interacting Genes (STRING; version 11.5) is a well-known tool for prediction of correlations between distinct proteins encoded by related genes (16). First, the DEGs were input into the STRING website. Then, a.tsv format file containing the PPI information was obtained. We used Cytoscape software to build PPI networks (17). Next, we used the Molecular Complex Detection (MCODE) tool to categorize the important gene modules. These modules had markedly interconnected groups in the PPI networks (18,19). We used the default parameters for all analyses. Finally, we used the above gene modules to carry out drug-gene interaction analysis.

Determination of the hub genes

We used the cytoHubba plugin of Cytoscape software to identify hub genes. This tool uses an algorithm to predict the critical nodes and subnetworks in an internal network. We used the maximal clique centrality (MCC) method to identify key genes.

Drug-gene interaction and functional analyses of potential hub genes

The Drug-Gene Interaction Database (DGIdb; https://www.dgidb.org) is an online resource for filtering drug-gene interactions. We used the module genes as search terms and determined the extant drugs targeting these genes.

Statistical analysis

We used the moderate t-test to identify DEGs. The GO and KEGG annotation enrichment results were analyzed by Fisher’s exact test. The analysis process diagram () shows the workflow of the study.

Figure 1

Framework of the data analyses. DEGs, differentially expressed genes; KEGG, Kyoto Encyclopedia of Genes and Genomes; DAVID, Database for Annotation, Visualization and Integrated Discovery; STRING, Search Tool for the Retrieval of Interacting Genes; MCODE, Molecular Complex Detection; DGIdb, The Drug-Gene Interaction Database.

Results

Identification of DEGs

According to the above criteria, a total of 481 DEGs between the normal and lymphedema groups were identified, including 307 upregulated genes and 174 downregulated genes (). A volcano map and heatmap depicting these DEGs are shown in .

Table 1

A total of 481 DEGs were identified from GSE4333, including 307 upregulated genes and 174 downregulated genes

Regulation direction	Gene name
Upregulated	Fcer1g, Dynlt1c/Dynlt1f, Tecpr1, B2m, Tlr2, Ifitm1, Cald1, Prkcz, Rptn, Nt5e, Rarres2, S100a8, Wdr18, Socs3, Csrp2, Syngr1, Ctsz, S100a9, Mfap4, Alox5ap, Cyba, Irf7, Tmem176b, Lcn2, Svep1, Srgn, Lyz2, Cdc25a, Trim30d, Erdr1, Calm1, Rac2, Ctsc, Trdn, C1qc, Txndc5, 2610528A11Rik, Ppp1r14b, Ctla2a, Ptx3, Marcks, Tmsb10, Stfa1, Gm10252/Ftl1, Sprr2a3/Sprr2a1, Ms4a6d, Stfa3, C1qa, Traf4, Capg, Sssca1, C3, Col5a2, Serpinb3a, Syngr3, Renbp, Cxcl12, Ccl6, Atxn7l1, Akap2, Cyth4, Ccl7, St13, Ifit1, Zpr1, Lox, Gnl2, 1700029K24Rik/Eno3, Slc28a2, Sec24d, Usp33, Chil4, Sepp1, Coro1a, Hdc, Lst1, Fam49b, Tpi1, Aif1, Lyn, Bgn, Sebox, Tnc, Pf4, C1qtnf3, Ctsb, Gngt2, Cbr2, Pgk2, Rnase1, Lgals1, Mfap2, H2-D1, Esd, Ndn, Tpm2, Cybb, Unc93b1, Dnaja1, Nelfe/Cfb, Fhl2, St6galnac4, Sdhd, Hexa, Cd53, Col5a1, Fth1, Col8a1, Ifitm2, Arhgdib, Ccr4, Lcp1, Plac8, Cd52, Col4a2, Plek, Dnajc13, Lgmn, Aqp1, 4930406D14Rik, Ehhadh, Ms4a6b, Selp, Col4a1, Cd44, Tmem176a, Tfpi, Ifitm3, LOC547349/5430410E06Rik/H2-L/H2-K1/H2-D1, Gpihbp1, Entpd1, Bicc1/4930533K18Rik, Ncf4, Ccl3, Tpm1, Cthrc1, Mt1, Rgs2, Ifngr1, Flot1, Arpc3, Tnfaip6, Fabp5, F2r, Gpm6b, Itm2a, Eln, Txlna, Slc35b2, Reg3g, Myeov2, Csta1, Anxa6, Col12a1, Slc38a1, Parp9, Mpp1, Oxct1, Impa1, Wdr83os, Ctla2b, Inpp5d, Ubl3, Psmb10, Zfp57, Cyc1, Vim, Gm10840/Limd2, Tmsb4x, Efemp2, Smpdl3a, Ftl1, Mob3a, Dnajb5, Ube2s, Cotl1, Szrd1, Cfh, Hspbap1, Ptma, Pnp, Fbln1, Prcp, Serping1, Anxa3, Snx32, Sat1, Mmp9, Sepw1, Fetub, Dynlt1f/Dynlt1b, Hsd17b4, Cystm1, Stim1, Igtp, Ifnar2, Ly86, Hsp90b1, Galnt10, Col27a1, Spink4, Gusb/Aebp1, Cmtm3, Clec4n, Ccnb1, Loxl1, Wnt7b, Cfp, Ly6c1, Bcl2a1a, Laptm5, Flrt2, H2-Q5, Samhd1, ar2, Psmb8, Thbs1, Krtcap2, Myh6, Slpi, 5730458M16Rik/Cd302, Gimap4, Lsp1, Itga11, Keg1, Snx5, Tmem55a, 1700071M16Rik/Pla2g7, Runx1t1, Tspan32, H2-DMa, Il10rb, Alx3, Ncf2, Msn, Ufsp2, H2-Aa, Ndufa3, Rpl7l1, Fn1, Cxcl16, Pigq, Irf1, LOC102641613/H2-K1, H2-Eb1, Nrp1, Prelp, Stab1, Ttr, Stc2, Ccr2, Gata1, Tmed9, Slc25a45, Usp36, Rps29, Lama4, Irf8, Ccl9, Sod3, AV051173, Marcksl1, Gm2518/Taf6l, Gcnt3, Psmb9, Spp1, Trdmt1, Cpe, Rsph9, Prmt1, Pgk1, Serpina1c/Serpina1b/Serpina1a, Klf13, Pkig, Mgp, Pdia4, Ube2b, C4a/C4b, Pglyrp1, Cox7b, Serf2, Cd14, Fabp2, Ddit3, Rnaset2a, Ttc17, Hnrnpa1, Cd74, Clec4b1, Mrps6/Slc5a3, 2900052L18Rik, Abca1, Sprr2j-ps, Ccdc53, Inpp5k, Rrp7a, Clps, Fcgbp, Nme1, Egfl6, Palld, Actb, Urah, Ubl5/Fbxl12, C1s1, Atp5b, H19
Downregulated	Dennd1b, Calb1, 1700025G04Rik, Mxd4, Selenbp1, Krtap13, Ocln, Mrps31, Trim35, Micall1, Ctnnbip1, Gzmb, Pcp2, Tssk3, Nle1, C78278, A330105O20Rik/Golga4, Slc16a11, Skap2, Stard8, Acoxl, Prr22, Relt, C2cd2, Cdkn2b, Nanog, Slx1b, Aff4, Aldh3a1, Tmem236, 1700007H22Rik, 1700022H16Rik, Acox3, Ppm1d, Cnn2, Pou3f2, Atp2a2, Mia3, Oog3, Apoe, Lyzl1, 2310034C09Rik, Tmem119, Rpl7, Aatk, Cilp2, 4930542C16Rik, Khdc1b, Ciart, Alox12e, Pck2, C77068, Igbp1b, Incenp, 1110064A23Rik, Kpna1, 4930550L05Rik/Fbxw2, Epcam, Gchfr, Slc22a4, Slc7a3, Arrdc2, Cox6b1, Zfp821/Ccdc7b, Wnt7a, Zmynd19, Gk5, Txndc2, Picalm, Omt2b, Mmaa, Gm8783, Lta, Fmo3, Ankrd33b, D10Wsu159e, Scn4b, Mitf, Chsy1, Sec14l4, Bbip1, Aig1, Dgat2, Nkx2-4, Phyhipl, Ufd, Ccnf, Kap, Efhb, Sdhaf3, Nop58, Il12a, Col15a1, Slurp2, Hnrnpu, 1700022F17Rik, Bglap, Il7r, Casp9, Golph3, Kifap3, Wdr1, 2310065H11Rik/2310057J18Rik, Gm10589/Pik3c2a, Tor4a, C77872, Pin1rt1, Kcnk7, Efnb1, Nf1, Smpd1, Smg5, 1700026H06Rik/Map7, Caml, Angel1, Slc12a4, Elovl6, Eif3f, Paxip1, Ostf1, Dag1, Igfbp5, Mycl, Ggps1, Tm9sf3, Ttll10, Wfdc21, Unc13c, Pygb, D1Ertd799e, Il9r, Nlrp4b, Slc2a13, Ppp5c, Slc3a2, Pds5b, 1700016D08Rik/Cluap1, Klhl32, Selplg, Wdr47, Rfpl4, Kcnd2, Cntn3, Fbxo36, Tmem225, Fam3c, Gm21685/Entpd4, Wnt1, Hbb-bs/Hbb-b2, Ppp1r12b, Tbx3, Ptger4, Rab11a, 1700125H03Rik, Tas1r3, Zxdc, Rgs14, D9Ertd167e, Cited4, Svil, Fxyd4, Nes, Pim3, Cypt12, Apoa1, Elmo2, Tmem159/4930560O18Rik, D6Ertd87e, Ahi1, 1810030O07Rik, LOC105242488/Tjap1, 2310061I04Rik, Ccdc80, Gm19935, Scn1, Col1a1, Trpv4, P4ha1, Col5a1top2A, Ccnb2, Nusap1, Iqgap3, Uhrf1, Aspm, Tubb3, Prc1, Tyms, Kiaa0101

DEGs, differentially expressed genes.

Figure 2

DEGs and heatmap from RNA-seq analysis. (A) Volcano plot showing the DEGs of GSE4333. Red dots represent significantly upregulated genes, blue dots represent significantly downregulated genes, and grey dots represent genes with no significant difference. (B) Heatmap showing the hierarchical clustering of DEGs in the normal and lymphedema groups. DEGs, differentially expressed genes.

GO and pathway enrichment analyses

We used DAVID to depict the enriched GO terms and signaling pathways of the DEGs. shows the top enrichment terms for the DEGs in the BP, CC, and MF categories. The BP terms were mainly related to immune system processes, the inflammatory response, and the innate immune response. In the CC category, the DEGs were primarily related to the extracellular exosome, extracellular region, and extracellular space components. In the MF category, the DEGs were significantly related to protein binding, extracellular matrix (ECM) structural constituent, and identical protein binding terms. According to the results of pathway analysis, the genes were highly involved in ECM-receptor interactions, Staphylococcus aureus infection, leishmaniasis, pertussis, and tuberculosis ().

Figure 3

GO analysis of DEGs between the normal and lymphedema groups. The bar plots show the top 10 significantly enriched terms for the DEGs in the BP, CC, and MF categories in terms of enrichment scores. BP, biological process; CC, cellular component; MF, molecular function; GO, Gene Ontology; DEGs, differentially expressed genes.

Figure 4

KEGG analysis of DEGs between the normal and lymphedema groups. The bar plot shows the top 10 significantly enriched pathways in terms of enrichment scores. ECM, extracellular matrix; KEGG, Kyoto Encyclopedia of Genes and Genomes; DEGs, differentially expressed genes.

PPI and module analysis

A total of 434 genes/nodes with 1,715 edges were included in the PPI networks. In contrast, 47 genes did not participate in the PPI networks (). Based on the MCODE analysis, 13 modules were identified in the network. The top 5 significant modules are shown in .

Figure 5

PPI network analysis. (A) PPI networks of the DEGs. (B) Significant gene module 1 in the PPI networks. PPI, protein-protein interaction; DEGs, differentially expressed genes.

Table 2

Top 5 significant modules

Module	Nodes	Edges	Score	Genes
Module 1	26	334	13.360	Aif1, Alox5ap, Arhgdib, Bgn, C1qc, Cd52, Cd53, Cd74, Coro1a, Cyba, Eln, Fcer1g, Inpp5d, Irf8, Laptm5, Lcp1, Lox, Ly86, Ms4a6b, Ms4a6d, Ncf2, Ncf4, Plek, Rac2, Selp, Thbs1
Module 2	7	36	6.000	Nle1, Gnl2, Nop58, Rpl7, Prmt1, Wdr18, Rpl7l1
Module 3	16	86	5.733	Actb, C1qa, Ccl3, Ccr2, Cd14, Col12a1, Col15a1, Col27a1, Col8a1, Il7r, Lcn2, Ly6c1, Lyn, Mmp9, S100a9, Vim
Module 4	16	80	5.333	Selplg, Nanog, Nes, Clec4n, Epcam, Ccl6, Lyz2, Apoe, Ccl7, Tlr2, S100a8, Cybb, Cfp, Fn1, Bglap, Nt5e
Module 5	14	52	4.000	Irf7, Sdhd, Psmb8, Atp5b, Ifitm2, Cyc1, Hsp90b1, Irf1, Cox6b1, Psmb9, Cox7b, Il10rb, Ifnar2, Ifitm1

Hub gene identification

Based on the MCC method, the connectivity degree values in the PPI network were used to determine the top 10 genes, which included Fc epsilon receptor Ig (Fcer1g), CD53 molecule (Cd53), neutrophil cytosolic factor 4 (Ncf4), Rac family small GTPase 2 (Rac2), lymphocyte antigen 86 (Ly86), lysosomal protein transmembrane 5 (Laptm5), interferon regulatory factor 8 (Irf8), neutrophil cytosolic factor 2 (Ncf2), allograft inflammatory factor 1 (Aif1), and coronin 1A (Coro1a).

Drug-gene interaction analysis

We selected the top module, module 1 (26 genes), for the analysis of drug-gene interactions (). The ALOX5AP gene, which is targeted by one existing candidate drug (fiboflapon), was found in module 1.

Discussion

Acquired lymphedema is a prevalent complication among patients who undergo ALND and axillary radiotherapy (20). Acquired lymphedema can be divided into acute acquired lymphedema and delayed acquired lymphedema. Acute acquired lymphedema, which affects numerous patients, often occurs within a few months after surgery, while delayed acquired lymphedema does not occur until a few years later. Patients’ quality of life can be negatively affected by acquired lymphedema because of chronic infections, swelling, secondary malignancies, and unsightly upper limbs. In most cases, acquired lymphedema can be treated with conservative methods, such as combined decongestive therapy, compression therapy, and patient education programs. Recently, the key conservative method of treating acquired lymphedema is complete decongestive therapy (CDT), which involves physical therapy, manual lymphatic drainage, and skin care. Lymphatic fluid accumulation in the tissues was reduced by this method (21). But it can only relieve acquired lymphedema. However, there is still no cure for lymphedema. Prevention including prophylactic use of compression sleeves and upper limb strength excise remains the key measure. Screening high-risk patients is critical for early intervention. Therefore, the current study attempted to reveal the hub genes, important pathways, and potentially effective drugs for acquired lymphedema.

The efficacy of therapeutic strategies for acquired lymphedema is still not ideal, primarily because the pathogenesis is not clear. Recent studies have found that the critical mechanisms involved in acquired lymphedema development might include inflammation, fibrosis, and adipocytokine signaling (22-24). However, few studies have explored the genetic mechanism of and potential targeted drugs for acquired lymphedema.

In this paper, the GSE4333 datasets containing microarray data for female hairless immunocompetent SKH-1 mice were utilized for gene expression analysis. We found a total of 481 DEGs, including 307 upregulated genes and 174 downregulated genes. Select bioinformatics methods, namely, GO term enrichment analysis, signaling pathway enrichment analysis, PPI network construction, hub gene selection, and drug-gene interaction analysis, were used to further characterize these DEGs.

We found interactions among proteins related to the DEGs and potential hub genes involved in the molecular mechanism through PPI network analysis. Then, we identified 10 genes with comparatively high connectivity levels (≥50): Aif1 (degree =90), Fcer1g (degree =88), Rac2 (degree =86), Irf8 (degree =80), Ly86 (degree =72), Ncf4 (degree =70), Ncf2 (degree =64), Laptm5 (degree =56), Cd53 (degree =52), and Coro1a (degree =50). Among them, Cd53, Ncf4, Rac2, Ncf2, and Coro1a were enriched in immune system processes. In addition, LY86, Laptm5, and Irf8 were related to inflammatory responses. The gene Cd53 belongs to the tetraspanin family, which is involved in the regulation of immune cell function (25,26). In our study, immune system processes were the most significantly enriched functional term. This may suggest that Cd53 acts on acquired lymphedema through an immune system process. Meanwhile, Ncf4, Rac2, Ncf2, and Coro1a participate in the same mechanism in the immune system (27-29). The gene Ly86 is also known as MD-1, and some studies have shown that MD-1 is related to inflammatory responses (30); Laptm5 has a positive feedback role in the inflammatory signaling pathway (31); and Irf8 adjusts the transcription of NLR family apoptosis inhibitory proteins related to inflammasome activation (32). In GO enrichment analysis, the annotated BP terms were associated with the inflammatory response. Therefore, the inflammatory response plays an active role in acquired lymphedema. The roles of Fcer1g and Aif1 in the mechanism of acquired lymphedema are still unclear, and further studies are needed.

In the BP annotation category, the DEGs were mostly related to the immune system and the inflammatory response. García Nores et al. suggested that in the lymphedematous region, the constituents of the skin innate immune system change during different phases. One study has shown that T-regulatory cells mediate local immunosuppression in lymphedema (33). Patients with chronic lymphedema can easily develop persistent infections and multiple neoplasms in the lymphedematous limb.

In the CC annotation category, the DEGs were markedly associated with extracellular exosomes and focal adhesion. Regarding extracellular exosomes, research has suggested that cancer exosomes trigger fibroblast-to-myofibroblast differentiation, which has major implications for the mechanisms underlying fibrotic diseases (34). Regarding focal adhesion, mutations in some transcription factors are predominantly associated with lymphedema and increased focal adhesions in human lymphatic endothelial cells (35).

In the MF category, the DEGs were significantly associated with fibronectin binding, protein binding, and ECM structural constituents. The ECM is a noncellular three-dimensional macromolecular network composed of collagens, fibronectin, laminins, several other glycoproteins, and other components. Herrera et al. (36) found that the ECM is a trigger of gradual fibrosis. Recently, additional studies have shown that chronic inflammation-induced fibrosis is important in the pathophysiology of acquired lymphedema. One basic study has shown that fibrosis gradually worsens and is irreversible in mouse tail models of acquired lymphedema (37). Rockson (38) found that inhibition of fibrosis can combat lymphoedema. With regard to protein binding, Leung (39) found that the accumulation and binding of protein molecules might damage the flow of lymph fluid and induce lymphedema and tissue oedemaedema. Regarding fibronectin binding, fibronectin is the basic component of the ECM. Yim et al. (40) found that in Korean breast cancer survivors, serum fibronectin 1 levels are associated with a high risk of lymphoedema.

According to KEGG analysis, the DEGs were highly related to the processes of ECM-receptor interaction and Staphylococcus aureus infection. Evidence has revealed that ECM receptors control the transcription of collagen by translocating to the nucleus (41,42). Some specific areas of ECM receptors contribute to tissue fibrosis by adjusting the production of ECM. Furthermore, lymphoedema associated with fluid accumulation in tissue is an important pathological feature of bacterial infections. However, how bacterial infections can be both the result and cause of lymphatic vessel dysfunction is still unclear. In recent studies, significant inhibition of both lymphatic vessel contraction and lymph velocity was observed in mice with active Staphylococcus aureus infection. When Staphylococcus aureus was cleared and inflammation was addressed, chronic inhibition of lymphatic function remained (43,44). These findings prompted us to focus on the potential mechanisms of fibrosis and bacterial inflammation.

The drug-gene interaction analysis of module 1 revealed that ALOX5AP is targeted by one existing candidate drug, fiboflapon. Fiboflapon is a 5-lipoxygenase-activating protein inhibitor and an effective oral anti-inflammatory drug (45). When the lymphatic system is injured, changes in normal tissue function occur due to inflammation. Subsequently, this can result in the lymphedema phenotype. Jiang et al. (46) investigated the interplay between inflammation and pathological tissue remodelling through inflammatory factors. In addition, Rockson (47) found that the inflammatory response is related to some gene expression pathways associated with acquired lymphedema. Hence, suppression of inflammatory reactions may delay the development of lymphedema. The findings of these studies suggest that anti-inflammatory drugs might be beneficial treatment agents for acquired lymphedema.

To date, the genes related to acquired lymphedema and drugs that can be applied to treat this condition have not been revealed. In a mouse model of acquired lymphedema, we identified 10 hub genes and one existing drug through comprehensive analyses. Among them, Cd53, Ncf4, Rac2, Ncf2, and Coro1a were enriched in immune system processes. Kataru et al. (48) found lymphatic injury could regulate immune response, which could promote tissue fibrosis. Furthermore, LY86, Laptm5, and Irf8 were related to inflammatory responses. A study showed that fluid cumulation by lymphatic injury might act as an initial signaling key to inflammation that eventually brings about the pathologic changes of lymphedema (49). So, immune system processes and immune response play important roles in acquired lymphedema. However, the functions of Fcer1g and Aif1 in the process of acquired lymphedema are still unclear. Besides, the fiboflapon-related drug maybe effective in treating acquired lymphedema. Future studies are needed to verify.

The study still had some limitations. We did not use human samples to verify the results due to the difficulty associated with acquiring humanized lymphedema samples. This kind of dataset did not exist in the GEO database. After searching in PubMed, we concluded that this is the first study to reveal the potential mechanism and therapeutic target of lymphedema using bioinformatic methods. Consequently, our findings provide insights to aid in the discovery of the molecular mechanism of acquired lymphedema.

Conclusions

Using bioinformatics methods, we identified the top 10 hub genes (Fcer1g, Cd53, Ncf4, Rac2, Ly86, Laptm5, Irf8, Ncf2, Aif1, and Coro1a), molecular mechanisms, and potential therapeutic drugs (fiboflapon) for acquired lymphedema.

The article’s supplementary files as