Abstract
Background
: Lymphedema is an intractable disease with no curative treatment available. Conservative treatment is the mainstay, and new drug treatment options are strongly needed. The purpose of this study was to investigate the effect of roxadustat, a prolyl-4-hydroxylase inhibitor, on lymphangiogenesis and its therapeutic effect on lymphedema in a radiation-free mouse hindlimb lymphedema model.
Methods and Results
: Male C57BL/6N mice (8–10 weeks old) were used for the lymphedema model. Mice were randomized to an experimental group receiving roxadustat or a control group. The circumferential ratio of the hindlimbs was evaluated, and lymphatic flow of the hindlimbs was compared by fluorescent lymphography up to 28 days postoperatively. The roxadustat group showed an early improvement in hindlimb circumference and stasis of lymphatic flow. The number and area of lymphatic vessels on postoperative day 7 were significantly larger and smaller, respectively, in the roxadustat group compared with the control group. Skin thickness and macrophage infiltration on postoperative day 7 were significantly reduced in the roxadustat group compared with the control group. The relative mRNA expression of hypoxia-inducible factor-1α (Hif-1α), vascular endothelial growth factor receptor-3 (VEGFR-3), vascular endothelial growth factor-C (VEGF-C), and Prospero homeobox 1 (Prox1) on postoperative day 4 was significantly higher in the roxadustat group compared with the control group.
Conclusions
: Roxadustat demonstrated a therapeutic effect in a murine model of hindlimb lymphedema through promotion of lymphangiogenesis through the activation of HIF-1α, VEGF-C, VEGFR-3, and Prox1, suggesting the potential of roxadustat as a therapeutic option in lymphedema.
Introduction
Lymphedema is a refractory condition caused by decreased lymphatic transport capacity, and characterized by local interstitial fluid accumulation, inflammation, and fatty degeneration of connective tissue. 1 Symptoms include pain, fatigue, and recurrent infections due to compromised immunity, which significantly reduces patients' activities of daily living. It affects 200 million people worldwide 2 and is a major complication after treatment of solid tumors.
Current treatments for lymphedema primarily include compression therapy, manual lymphatic drainage, and other conservative therapies, rather than rebuilding the destroyed lymphatic system. Although lymphovenous anastomosis and vascularized lymph node transplantation have been reported to be effective for lymphedema, only specialized centers are able to perform these procedures, and their effectiveness remains controversial. 3 The pathophysiology of lymphedema remains to be clarified, and there is an urgent need for new drug therapies.
Roxadustat, a novel hypoxia-inducible factor (HIF) prolyl 4-hydroxylase inhibitor, is of interest as a treatment for renal anemia in patients with chronic kidney disease. In addition, it has been reported that roxadustat exerts therapeutic effects on cutaneous wound healing,4,5 cardiac disease, 6 retinal disease, 7 spinal cord injury, 8 Parkinson's disease, 9 and malignant tumor growth10,11 by stabilizing hypoxia-inducible factor-1α (HIF-1α) expression. These findings show the great therapeutic potential of roxadustat in diseases related to HIF-1α signaling.
In murine lymphedema, lymphatic stasis and inflammation stabilize HIF-1α and result in high levels of HIF-1α expression, indicating that this is required for reparative lymphangiogenesis. 12 HIF-1α has also been shown to promote lymphangiogenesis by inducing the expression of vascular endothelial growth factor receptor-3 (VEGFR-3) and vascular endothelial growth factor-C (VEGF-C) in lymphatic endothelial cells.12,13 Although it has been reported that the deletion or inhibition of HIF-1α exacerbated edema in a murine model of tail lymphedema postoperatively,12,14 HIF-1α stabilization has not been studied in an animal model of lymphedema. Therefore, in this study, the therapeutic effects of roxadustat were investigated in a murine model of hindlimb lymphedema that we previously developed. 15
Materials and Methods
Cell culture and treatment
Endothelial Cell Basal Media MV2 (PromoCell, Heidelberg, Germany) supplemented with penicillin-streptomycin and Growth Medium MV2 Supplement Mix (PromoCell) was used to culture human dermal lymphatic endothelial cells (HDLECs) (PromoCell), and the HDLECs from passages 3–5 were used. Roxadustat (Selleck Chemicals, Houston, TX) was dissolved in dimethyl sulfoxide (DMSO) and added to the medium at final concentrations of 0, 5, 10, and 25 μM.
Proliferation assay
To assess the effect of roxadustat on HDLEC proliferation, a Cell Counting Kit-8 (CCK-8) assay (Dojindo Molecular Technologies, Japan) was performed, following a previously reported protocol. 4 In brief, cells were seeded on a 96-well plate at a density of 2 × 103 cells/well with 0, 5, 10, or 25 μM of roxadustat added to the medium and incubated for 5 days. From days 0–5, 10 μL of CCK-8 solution was added to each well, after which the samples were incubated for 60 minutes at 37°C. The absorbance was detected at 450 nm on the Infinite 200 PRO microplate reader (Tecan Japan, Kawasaki, Japan).
Scratch wound assay
HDLECs were seeded into 6-well plates at a density of 5.0 × 105 cells/well and cultured at 37°C for 24 hours. After formation of a confluent monolayer, a sterile pipette tip was used to make a linear scratch. Phosphate-buffered saline was used to wash away the cellular residue, then 2 mL of medium containing 0, 5, 10, or 25 μM roxadustat was added to each well, and the plates were incubated at 37°C. Photographs were taken at 0 and 24 hours, and the amount by which each scratch had closed and the migration rate were determined using ImageJ software (National Institutes of Health, Bethesda, MD).
Tube formation assay
Pretreatment of HDLECs with 0, 5, 10, or 25 μM roxadustat was carried out for 24 hours. Then, the HDLECs were mixed with medium and seeded into a 24-well plate precoated with 100 μL of Geltrex (Life Technologies, Grand Island, NY) at 8 × 104 cells/well, and incubated at 37°C for 6 hours. After microscopic observation, the total length of each tube was calculated using ImageJ software (National Institutes of Health). 4
Animals
All experiments involving animals in this study were approved by the Hokkaido University Institutional Animal Care and Use Committee. Male C57BL/6N mice, 8–10 weeks of age (SLC, Tokyo, Japan), were kept in a room maintained at 24°C under a 12-hours light/dark cycle with free access to food and water.
Murine model and drug administration
The murine model of hindlimb lymphedema was created using our previously established method. 15 In brief, a circumferential incision was made in the left inguinal skin, followed by resection of the inguinal lymph node and surrounding fat pad. After tying the prenodal and postnodal lymphatic vessels, the popliteal lymph node and surrounding fat pad were resected. A 1-mm–thick silicone sheet was used to fashion a 3-mm–wide rectangular splint, which was then placed in the inguinal wound and affixed to the skin and underlying muscle.
The mice were randomly assigned to either an experimental group receiving 25 mg/kg roxadustat (Selleck Chemicals) in DMSO4,7 or a control group receiving only DMSO, and were given intraperitoneal injections every 2 days, including the day of surgery, for up to 2 weeks. A total of 36 mice were used and sacrificed on day 28 (12 mice for edema assessment and fluorescence lymphatic imaging), on day 7 (12 mice for histology), and on day 4 (12 mice for RNA isolation).
Edema assessment
Lymphedema formation in the hindlimb was evaluated quantitatively on days 0, 2, 4, 7, 10, 14, 17, 21, 24, and 28 after surgery (n = 6 per group), at which time the circumference of the musculotendinous junction of the gastrocnemius muscle was measured bilaterally, and the circumference ratio was calculated according to the following formula:
Circumference ratio = (Treated hindlimb circumference/Untreated hindlimb circumference) × 100%.
Fluorescence lymphatic imaging
Lymphatic structures in the hindlimbs of the control and roxadustat groups were compared by fluorescent lymphography every week for 4 weeks postoperatively (n = 6 per group). Before imaging, isoflurane was used to anesthetize the mice, and their residual fur was removed. A 5-μL volume of indocyanine green solution (2.5 mg/mL Diagnogreen in distilled water; Daiichi Sankyo Company, Ltd., Tokyo, Japan) was subcutaneously injected into both paws with a 26-gauge needle. Imaging was performed using a near-infrared fluorescence camera system (Photodynamic Eye; Hamamatsu Photonics, Hamamatsu, Japan) 15 minutes after indocyanine green injection. The coverage of the fluorescent area (i.e., the dermal backflow at the thigh) was measured using ImageJ software (National Institutes of Health).
Histology
Six mice per group were sacrificed at 7 days after surgery, and skin samples were obtained. To minimize the impact of inflammation due to wound healing, the samples were taken from an area 6 mm distal to the inguinal wound. Then, skin sections were fixed in 4% paraformaldehyde, embedded in paraffin, and subjected to immunohistochemical or Elastica–Masson staining. Immunohistochemical staining was carried out to evaluate lymphatic endothelial cells and macrophage infiltration. Sections were incubated overnight with antibodies against LYVE-1 (Abcam, Inc., Cambridge, MA) and F4/80 (Cedarlane, Burlington, Canada).
Histofine (Nichirei, Tokyo, Japan), which was developed using 3,3′-diaminobenzidine, was used as the secondary antibody. Elastica–Masson staining was performed to evaluate skin thickness. A whole-slide scanner (Nano Zoomer Digital Pathology; Hamamatsu Photonics) was used to capture digital images of the slides, which were then visualized using NDP.view2 software (Hamamatsu Photonics). In accordance with a previous report, 16 lymphatic vessels were defined as vessels with an immunopositive endothelium and a vascular lumen.
Lymphatic vessels were identified by two examiners other than the first author after the specimens were randomized, and those identified by both examiners were designated as lymphatic vessels. LYVE-1–stained sections were scanned at low magnification (40 × ), and “hot spots” with the greatest number of lymphatic vessels were noted. Then, five of these areas in each section were observed at high magnification (200 × ), and the average number of lymphatic vessels within the hot spots in each field of view was calculated. 17
In addition, the average area of the lymphatic lumen in each field of view was calculated using ImageJ software (National Institutes of Health). Quantification of the F4/80 positive area was performed using ImageJ software (National Institutes of Health) by observing 10 randomly selected high-magnification fields of view (400 × ). The thickness of skin in the roxadustat and control groups was measured as the distance from the epidermis to the dermal-fat junction by two examiners other than the first author, using ImageJ software (National Institutes of Health) after the specimens were randomized.
RNA isolation and real-time polymerase chain reaction
To isolate the RNA, six mice per group were sacrificed at 4 days after surgery. Skin and subcutaneous tissue samples were collected from the hindlimb at an area 6 mm distal to the inguinal wound and frozen. RNeasy Fibrous Tissue Mini Kit (Qiagen, Hilden, Germany) was used to extract total RNA from the skin samples, after which a High Capacity RNA-to-cDNA Kit (Applied Biosystems, Foster City, CA) was used to reverse transcribe the RNA to cDNA. Next, Power SYBR Green PCR Master Mix and a StepOnePlus Real-Time PCR System (Applied Biosystems) were used to perform real-time quantitative reverse-transcription polymerase chain reaction (PCR), and the relative levels of PCR products were calculated using the ΔΔCt method. 15 Examinations of each sample were performed three times. Primers for reverse-transcription PCR are listed in Table 1.
Reverse-Transcription Polymerase Chain Reaction Primer Sequences and Product Size
GAPDH, glyceraldehyde-3-phosphate dehydrogenase; Hif-1α, hypoxia-inducible factor-1α; Prox1, Prospero homeobox 1; VEGF-C, vascular endothelial growth factor-C; VEGFR-3, vascular endothelial growth factor receptor-3.
Statistical analysis
Differences between the two groups were analyzed by Student's t-test, while differences among three or more groups were analyzed by one-way analysis of variance followed by a Tukey–Kramer multiple comparisons test. Results from multiple experiments are expressed as mean values ± standard error. All statistical analyses were performed using JMP software ver. 16.0.0 (SAS Institute, Inc., Cary, NC), with p < 0.05 considered to indicate statistical significance.
Results
In vitro assays
Roxadustat was found to stimulate the proliferation of HDLECs in a dose-dependent manner (Fig. 1A), and significantly enhanced their motility (p < 0.05) (Fig. 1B, C). Roxadustat also significantly enhanced HDLEC tube formation (p < 0.05) (Fig. 1D, E).

Effects of roxadustat on HDLECs in vitro.
Edema assessment
Edema disappeared grossly on days 21 and 28 in the roxadustat and control groups, respectively (Fig. 2A). The hindlimb circumference ratio reached a maximum on day 4 in both groups, and then gradually approached 100%. As in our previous study, 15 the hindlimb circumference ratio reached 100% on day 28 in the control group and on day 21 in the roxadustat group. The circumference ratio was significantly lower in the roxadustat group compared with the control group from postoperative days 2–24 (Fig. 2B).

Effects of roxadustat on edema.
Fluorescence lymphatic imaging
The control group showed diffuse dermal backflow throughout the treated hindlimb from 1 to 4 weeks postoperatively. In contrast, the roxadustat group showed dermal backflow, but with lower fluorescence compared with the control group (Fig. 3A). Compared with the control group, the coverage of the fluorescent area at the thigh at 2 weeks postoperatively was significantly reduced in the roxadustat group (*p < 0.05) (Fig. 3B).

Effects of roxadustat on lymphatic flow.
Histology
The roxadustat group had more lymphatic vessels (Fig. 4A, B) and a smaller lymphatic lumen area (Fig. 4A, C) compared with control group. There were also significantly more LYVE-1–positive lymphatic vessels in the roxadustat group (8.5 ± 0.5) compared with the control group (5.2 ± 0.4) (Fig. 4B, p < 0.05).

Effects of roxadustat on lymphatic vessels.
The area of lymphatic lumen in the roxadustat group (334 ± 38 μm2) was significantly less than that in the control group (1385 ± 163 μm2) (Fig. 4C, p < 0.05). Furthermore, the F4/80-positive area was significantly smaller in the roxadustat group (1.9% ± 0.3%) compared with the control group (3.1% ± 0.4%) (*p < 0.05) (Fig. 5A, B). Skin thickness in the roxadustat group (240 ± 18 μm) was significantly less than that in the control group (296 ± 12 μm) (*p < 0.05) (Fig. 5C, D).

Effects of roxadustat on macrophage infiltration and skin thickness.
Real-time PCR
The relative mRNA expression levels of Hif-1α, VEGFR-3, VEGF-C, and Prospero homeobox 1 (Prox1) were significantly higher in the roxadustat group (1.6 ± 0.2, 4.3 ± 0.7, 2.5 ± 0.3, and 2.8 ± 0.4, respectively) compared with the control group (Fig. 6, p < 0.05).

Real-time reverse-transcription polymerase chain reaction analysis of the skin and subcutaneous tissues demonstrated that Hif-1α, VEGF-C, VEGFR-3, and Prox1 were expressed at significantly higher levels in the roxadustat group compared with the control group at postoperative day 4 (*p < 0.05). Error bars = standard error of the mean. Hif-1α, hypoxia-inducible factor-1α; Prox1, Prospero homeobox 1; VEGF-C, vascular endothelial growth factor-C; VEGFR-3, vascular endothelial growth factor receptor-3.
Discussion
In this study, we used a novel radiation-free murine model of hindlimb lymphedema, which we had previously developed. 15 The lack of appropriate animal models is one of the reasons the elucidation of the pathogenesis of lymphedema has been hampered. Many studies of lymphedema in animal models have been conducted mainly in mouse-tail12,14 and hindlimb edema models,18,19 but because the mouse-tail edema model does not involve lymph node excision and many murine models of hindlimb edema use radiation, a suitable clinical lymphedema animal model for molecular biological studies did not exist. Our lymphedema model does not use radiation, thus eliminating its effects. In addition, the lymph node excision model more closely mimics clinical lymphedema.
Although there have been several reports on the most appropriate circumferential measurement site for hindlimb lymphedema, including a point located at 5 or 6 mm from the heel,15,18 we took circumferential measurements at the musculotendinous junction of the gastrocnemius muscle to ensure anatomic consistency at each time point, 17 taking into account the growth of the mice during the experiment, and evaluated it in comparison with the untreated side.
The circumferential ratio, which indicates the degree of edema, improved faster in the roxadustat group compared with the control group. Fluorescent imaging evaluation of lymphatic flow showed significantly lower fluorescence after 2 weeks in the roxadustat group, whereas dermal backflow, which is indicative of lymphatic stasis, was diffuse throughout the hindlimb over the 4-week postoperative period in the control group, suggesting that lymphatic stasis was improved in the roxadustat group.
Yamamoto et al 20 reported that the degree of dermal backflow was positively correlated with clinical severity, and the present results are consistent with theirs. The formation of lymphatic collateral vessels and enhanced function of the draining collecting vessels were considered as possible reasons for the reduced dermal backflow in the roxadustat group, but neither could be confirmed in our study due to the low resolution of the fluorescence camera system.
Microscopically, dilation of lymphatic vessels has been reported at the early stages of lymphedema, 21 and our previous study involving a murine model of hindlimb lymphedema revealed cutaneous and subcutaneous lymphatic dilation. 15 In this study, we found that on postoperative day 7, the number of cutaneous and subcutaneous lymphatic vessels was significantly larger, and their lumen area was significantly smaller in the roxadustat group compared with the control group. This suggests that lymphangiogenesis was more pronounced in the roxadustat group, and that the dilation of lymphatic vessels, which indicates lymphatic stasis, was improved.
In addition, skin thickness in the roxadustat group was significantly less than that in the control group. Furthermore, significantly less macrophage infiltration was found in the roxadustat group than in the control group, which is consistent with previous reports. 22
HIF-1 is a transcriptional activator of various genes that play a role in the adaptive response of cells to hypoxia. Dysfunction in systems that regulate HIF-1 activity has been implicated in the pathogenesis of malignancies and other diseases. 23 HIF-1 has a heterodimeric structure comprising oxygen-sensitive HIF-1α and constitutively expressed HIF-1β subunits. Under normoxic conditions, the ubiquitin–proteasome pathway rapidly degrades the HIF-1α subunit after hydroxylation of proline residues by HIF prolyl-4-hydroxylases. 7
The HIF system targets a wide range of genes, and in particular, increased HIF expression through prolyl hydroxylation domain protein inhibitors such as roxadustat has been reported to have therapeutic effects in cutaneous wound healing,4,5 cardiac disease, 6 retinal disease, 7 spinal cord injury, 8 Parkinson's disease, 9 and malignant tumors,10,11 and is attracting attention as a new pharmacological target.
HIF-1α promotes lymphangiogenesis by inducing VEGFR-3 and VEGF-C expression in lymphatic endothelial cells12,13 and by increasing VEGFR-3 expression through the transcriptional activation of Prox1. 24 It has also been reported that inhibition or deletion of HIF-1α exacerbates edema in a murine model of tail lymphedema postoperatively.12,14 These findings suggest that HIF-1α expression may be involved in postoperative lymphatic regeneration. In this study, Hif-1α, VEGFR-3, VEGF-C, and Prox1 mRNA expression levels were significantly elevated in the roxadustat group, suggesting that roxadustat stabilizes Hif-1α expression and promotes lymphangiogenesis through the increased expression of VEGFR-3, VEGF-C, and Prox1.
VEGF-C is a critical regulator of lymphangiogenesis and controls numerous lymphatic cell processes through VEGFR-3 signaling, including differentiation, migration, and survival. 25 It has been shown that the local administration of VEGF-C or gene therapy involving adenoviral vectors significantly increases lymphangiogenesis, and decreases swelling in animal models of both primary and secondary lymphedema.26,27 However, the application of VEGF-C therapy in the clinical setting has yet to be realized because VEGF-C acts as a key regulator of the tumor microenvironment, potentially increasing the risk of recurrence and metastasis.25,28
Previous studies have reported that roxadustat is noncarcinogenic, 29 and it was found not to promote the initiation, progression, or metastasis of tumors in a VEGF-sensitive spontaneous breast cancer model. 30 Meanwhile, a phase 2 study is being conducted to evaluate roxadustat for the treatment of anemia in patients receiving chemotherapy for nonmyeloid malignancies (NCT04076943).
In addition, there are reports that roxadustat inhibits tumor growth by enhancing macrophage phagocytosis of malignant tumors, 10 and by increasing sensitivity to chemotherapy and inhibiting tumor growth in a murine model of malignancy. 11 Although roxadustat is reported to have antitumor effects as described above, an increase in VEGF-C mRNA expression was also found in this study, and thus its effects on malignant tumors warrant further study.
This study has some limitations that merit consideration. First, the murine model of hindlimb lymphedema examined in this study does not address chronic lymphedema and is based on acute lymphedema. Human lymphedema progresses slowly and does not improve spontaneously, whereas in animal models it often resolves spontaneously. In this model, gross edema resolves in ∼1 month. Therefore, although this model anatomically approximates human lymphedema, it is not a perfect mimic, including the course of the disease. Second, although data in numerous animal studies have demonstrated that roxadustat does not induce carcinogenesis 29 or promote cancer, 30 further safety verification through clinical trials is needed to ensure its safety against malignant tumors in the clinical context.
Conclusions
In summary, roxadustat exerted a therapeutic effect in a murine model of hindlimb lymphedema by promoting lymphangiogenesis through activation of HIF-1α, VEGFR-3, VEGF-C, and Prox1. This suggests that roxadustat induces reparative lymphangiogenesis against impaired lymphatic function and thus shows promise as a therapeutic option for lymphedema.
Footnotes
Acknowledgments
The authors thank Kohei Oashi and Daisuke Iwasaki for their advice regarding this study.
Authors' Contributions
Y.H. contributed to conceptualization, investigation, data curation, and writing—original draft preparation. M.O. provided methodology and funding acquisition. E.F. contributed to methodology and validation. K.I. performed supervision and data curation. T.M. assisted with validation and visualization. M.H. performed data curation and formal analysis. Y.Y. contributed to conceptualization and project administration. T.M. assisted with conceptualization, investigation, supervision, and validation.
Author Disclosure Statement
No competing financial interests exist.
Funding Information
This study was funded by the Japan Society for the Promotion of Science (grant no. JP20K09841).
