Abstract
Although artesunate (ART) is generally accepted as a safe and well-tolerated first-line treatment of severe malaria, cases of severe side effects and toxicity of this compound are also documented. This study applied larval zebrafishes to determine the acute toxicity and efficacy of ART and performed RNA-sequencing analyses to unravel the underlying signaling pathways contributing to ART’s activities. Results from acute toxicity assay showed that a single-dose intravenous injection of ART from 3.6 ng/fish (1/9 maximum nonlethal concentration) to 41.8 ng/fish (lethal dose 10%) obviously induced pericardial edema, circulation defects, yolk sac absorption delay, renal edema, and swim bladder loss, indicating acute cardiotoxicity, nephrotoxicity, and developmental toxicity of ART. Efficacy assay showed that ART at 1/2 lowest observed adverse effect level (LOAEL) exerted cardioprotective effects on zebrafishes with verapamil-induced heart failure. Artesunate significantly restored cardiac malformation, venous stasis, cardiac output decrease, and blood flow dynamics reduction. No adverse events were observed with this treatment, indicating that ART at doses below LOAEL was effective and safe. These results indicate that ART at low doses was cardioprotective, but revealed cardiotoxicity at high doses. RNA-sequencing analysis showed that gene expression of
Introduction
As complementary and alternative medicine, traditional Chinese medicine (TCM) gains increasing popularity and attracts growing attention around the world.
Traditional Chinese medicine herbs exert only few side effects and adverse events were less frequently reported than by conventional Western medicines due to its particular pharmacological properties and appropriate use by TCM clinical physicians. 17- 19 The quintessence of TCM is the combination and interaction of different herbs in one formula, in which possible side effect exerted by one herb are neutralized by another herb in the herbal mixture, while the efficacy of each herb is synergistically enhanced. In fact, many herbal components have certain toxicities, causing damages on nervous, liver, renal, respiratory, and reproductive systems. 20 Without the interactions in formula, the use of single herbal components, such as ART, may bear some risks in clinical application and side effects may occur in cases of overdose. Recently, concerns about the safety of ART have been raised due to the potential toxicity observed during its treatment. 21- 23 Therefore, in-depth knowledge on the pharmacological efficacy-safety profile of ART is urgently needed for its clinical application. Here, we hypothesized that ART’s dosage may be a key factor that accounts for its efficacy-safety relationship, that is, ART may induce toxicity at high doses, but exerts safe efficacy at low doses. To verify this hypothesis, we applied a larval Zebrafish model to study the efficacy-toxicity relationship of ART.
Zebrafish (
Materials and Methods
Chemicals and Reagents
Artesunate (Batch number: B20992) with a purity of ≥ 98% was purchased as a powder from Shanghai Yuanye Biotechnology Co, Ltd (Shanghai, China). Verapamil hydrochloride (Batch number: K1629079) with a purity of ≥ 98% as a powder was purchased from Shanghai Aladdin Reagent Co, Ltd (Shanghai, China). Dimethyl sulfoxide (DMSO) was obtained from Sigma-Aldrich (Taufkirchen, Germany). Trizol reagent was purchased from Takara Biotechnology Co, Ltd (Dalian, China). Artesunate and verapamil hydrochloride were dissolved with DMSO and ultrapure water for experimental use, respectively.
Zebrafish Handling
Wide-type AB strain of Zebrafish was purchased from China Zebrafish Resource Center (CZRC), Institute of Hydrobiology, Chinese Academy of Sciences (Wuhan, China) and bred by Hunter Biotechnology, Inc. (Hangzhou, China). All zebrafishes were accredited by the Association for Assessment and Accreditation of Laboratory Animal Care (AAA LAC) International (SYXK2012-0171). After natural pair-mating and reproduction by adult zebrafishes, larval zebrafishes at the 2 days postfertilization stage (dpf) were generated and housed in a light-controlled aquaculture facility with a standard 14:10 hours day/night photoperiod and fed with live brine shrimp twice a day and dry flake once a day. The water temperature was maintained at 28°C (0.2% instant ocean salt, pH 6.9-7.2, conductivity 480-510 μS/cm, and hardness 53.7∼71.6 mg/L CaCO3).
Determination of Maximum Nonlethal Concentration and Lethal Concentration 10% of ART
Totally 270 larval zebrafishes (2 dpf) were grouped (30 fishes per group) and distributed into 6-well plates in 3 mL fresh fish water for each well. Eight concentrations of ART (0, 5, 10, 20, 30, 40, 60, 80 ng/fish) were used to intravenous (IV) inject zebrafishes in each group. The injection volume was 2 nL/fish. Dead zebrafishes were daily recorded and removed. The dose-mortality response curve was generated by using Origin 8.0 (OriginLab, Northampton, Massachusetts). Maximum nonlethal concentration (MNLC) and Lethal concentration 10% (LC10) were determined by logistic regression calculation.
Acute Toxicity Assay
To evaluate the acute toxicity of ART, 180 larval zebrafishes (2 dpf) were grouped (n = 30) and IV injected with 1/9 MNLC, 1/3 MNLC, MNLC, and LC10 of ART as well as ultrapure water (normal control [NC] group) and DMSO (solvent control group) for 3 days, respectively. The injection volume was 2 nL/fish. After treatment, the heart, brain, mandible, eyes, liver, intestine, spine, muscle, color, blood circulation, and all abnormal symptoms of each Zebrafish were observed under the microscope. The occurrence of edema, hemorrhage, and thrombosis were also observed. All abnormal phenotypes were statistically recorded by double-blind evaluation.
Efficacy Assay of ART
To evaluate the cardioprotective effect of ART, acute heart failure model of zebrafish was established by in administration of verapamil hydrochloride. Totally 210 larval zebrafishes (2 dpf) were grouped (n = 30) and IV injected with 0, 1.25, 2.5, and 5 ng/fish of ART, respectively, for 3 days. The injection volume was 2 nL/fish. The heart area was observed under the dissecting microscopy (SZX7; Olympus, Tokyo, Japan), and the cardiac output, venous stasis, and blood flow rate were calculated and obtained with Heart Blood Flow Analysis System (Zebralab3.3 PB2084C; ViewPoint, Ltd, Lyon, France). The bodies of zebrafishes were collected and 2 replicates were applied for RNA-sequencing analysis.
RNA-Sequencing Analysis
Total RNA was isolated from zebrafishes by using Trizol reagent according to the instructions of the manufacturer and treated with RNase-free DNase I to remove the residual DNA. The quality and integrity of the RNA sample were confirmed using an Agilent 2100 Bioanalyzer (Agilent, Waldbronn, Germany). The concentration of total RNA was quantified using a NanoDrop 2000 spectrophotometer (Thermo Scientific, Massachusetts). For complementary DNA (cDNA) library preparation, 3 μg total RNA were captured by Dynabeads Oligo (dT) (Life Technologies, New York) and sheared to fragments of ∼200 bp. Reverse transcription was performed by Superscript III cDNA Synthesis Kit. The cDNA was end-repaired, A-tailed and ligated to Illumina sequencing adapters and amplified by polymerase chain reaction (PCR). Library preparation was performed by TruSeq RNA LT V2 Sample Prep Kit (Illumina, San Diego, California). The sequencing library was qualified by Qubit 2.0 (Life technologies) and Agilent 2100 Bioanalyzer. The DNA sequencing was performed on the Illumina XTen Sequencing System (Illumina). Real-time image analysis and base calling were performed by HiSeq Control Software (HCS 1.5).
The raw data (raw reads) obtained from the Illumina HiSeq platform were processed by removing the adaptors, sequences with uncertain bases, low-quality sequences, and sequences of less than 50 bp to generate clean data (clean reads). The clean reads from the Fastq files were mapped to the Zebrafish reference genome using the Spliced Transcripts Alignment to Reference (STAR) software. Differential expression analysis of the different groups was performed with biological replicates using the DESeq software. Value of
Real Time PCR (RT-PCR) Assay
Real Time PCR assay was performed to verify the expression of candidate genes by using an ABI QuantStudio 7 Flex Real-Time PCR System (Applied Biosystems, California). The total RNA of zebrafishes was extracted using Trizol reagent and synthesized to cDNA via reverse transcription. Real Time PCR reaction system had a 20.0 μL volume: 10.0 μL SYBR Premix Ex Taq II (Tli RnaseH Plus), 0.8 μL PCR Forward Primer, 0.8 μL PCR Reverse Primer, 2.0 μL template cDNA, 0.4 μL ROX Reference Dye, and 6.0 μL ddH2O. The reaction condition was set to 95°C for 30 seconds initial denaturation, 40 cycles of 95°C for 5 seconds denaturation, 60°C for 35 seconds annealing, and 72°C for 40 seconds extension. At the end of each reaction, a melting curve analysis was performed. Data were normalized to the expression of 18 S and presented by using 2−ΔΔCT method (Table 1).
Primer Sequences Used for Real Time PCR Analysis.
Results
Acute Toxicity of ART
ART-induced dose-mortality curve was shown in Figure 1 (upper). No death of zebrafishes was observed within ART doses from 5 to 30 ng/fish, and no survival was seen with ART at dose of 80 ng/fish. Accordingly, MNLC and lethal dose 10% (LD10) were estimated as 32.4 and 41.8 ng/fish. During the double-blind observation, the lowest observed adverse effect level (LOAEL) of ART was estimated as 2.5 ng/fish. For acute toxicity assay, 1/9 of MNLC (3.6 ng/fish), 1/3 of MNLC (10.8 ng/fish), MNLC, and LD10 were adopted.

Acute toxicity observation of ART in zebrafishes. Upper: dose-mortality response curve with MNLC and LD10; Middle: microscopic observation (J: mandible; E: eye; H: heart; L: liver; Y: yolk sac; In: intestine; Sb: swim bladder) on A: (normal control group); B: (solvent control group); C: (1/9 MNLC ART treated group); D: (1/3 MNLC ART treated group); E: (MNLC ART treated group); and F: (LD10 ART treated group); Lower: observation of renal edema, muscle degeneration (muscle texture abnormality), and cardiac abnormality. LD10 indicates lethal dose 10%, MNLC, maximum nonlethal concentration.
As shown in Figure 1 (middle and lower), morphological abnormalities of zebrafishes were observed with ART treatments. In both NC and solvent control groups, abnormal phenotypes were not found in the zebrafishes. If treated with increasing ART doses, the incidences of abnormal manifestations, such as pericardium edema, slowed circulation defect (venous stasis), and yolk sac absorption delay, were found mainly in the cardiovascular system of zebrafishes, indicating dose-dependent acute toxicity of ART, especially cardiotoxicity. The hearts in these zebrafishes suffered from dilation and edema and were filled with turbid fluid, which seemed to be bigger than those of normal zebrafishes. Particularly, almost all fishes (incidence of >96%) suffered from circulation defects upon ART treatment at a dose of LD10. Artesunate at doses from MNLC to LD10 caused severe renal edema and swim bladder loss, indicating acute nephrotoxicity and developmental toxicity. Besides, ART at LD10 exerted toxicity in almost all organs in zebrafishes, including brain, mandible, eyes, muscle, and body length (Table 2).
Artesunate-Induced Abnormal Manifestations Their Occurrence Rate in Zebrafishes.a
a “−” indicates 0/30.
Cardioprotective Effect of ART
To determine the effective and safe dose range and verify the cardiotoxicity of ART, another acute toxicity assay was conducted with ART treatment at doses of 1.25, 2.5, and 5.0 ng/fish. As shown in Figure 2A, typical manifestation of cardiotoxicity, such as venous congestion and stasis, was observed in zebrafishes treated with 2.5 and 5.0 ng/fish ART. Thus, 1.25 ng/fish (1/2 LOAEL) was safe for the cardiovascular system and selected for the efficacy assay. As shown in Figure 2B and lower, heart dilatation, venous stasis, cardiac output decrease, and blood flow dynamics reduction were observed in zebrafishes treated with verapamil hydrochloride, of which treatment induced heart failure in 100% zebrafishes. No observable toxicity or death was found in zebrafishes after treatment of verapamil hydrochloride. After 1.25 ng/fish ART treatment, those verapamil-induced heart failure associated symptoms in zebrafishes were significantly restored toward normal levels (Figure 2B and lower panel), indicating a profound cardioprotective effect of ART at 1.25 ng/fish.

Cardiotoxicity and cardioprotective effect of ART in zebrafishes. (A) observation of circulation defect (venous congestion and stasis); (B) observation of heart failure phenotype. Lower: calculation of heart area, cardiac output, venous stasis, and blood flow rate. NC: normal control group; Model: acute heart failure group. Data (mean ± SD) with different lowercase letters are significantly different with each other at LSD multiple comparisons. LSD indicates Fisher's least significant difference; SD, standard deviation.
RNA-Sequencing-Based Transcriptional Profile of ART
Nearly 56 to 76 million RNA-sequencing reads for zebrafishes in the NC, model, and ART groups were obtained. For expression measurements, more than 92% of total reads mapped to reference database genes. Statistical analysis revealed 630 (NC vs model) and 368 (model vs ART) differentially expressed genes (

Heat map overview of the top 50 significantly expressed genes in zebrafishes.

Volcano plot and KEGG pathways for the comparison between the different treatment groups as assayed by RNA-sequencing analysis. KEGG Kyoto Encyclopedia of Genes and Genomes.
As shown in Table 3, besides
Significantly Expressed Genes Which Are Related to Cardiotoxicity or Cardioprotection.
Abbreviation: SE for log2 (fold change).
Relative Messenger RNA Expression of Candidate Genes in Zebrafishes
Real time PCR was performed to assess the relative messenger RNA expression of selected candidate genes in zebrafishes. As shown in Figure 5, when compared with NC,

Relative messenger RNA expression of candidate genes in zebrafishes. Data (mean ± SD) with different lowercase letters are significantly different with each other at multiple comparisons. SD indicates standard deviation.
Discussion
During the Vietnam War, numerous Vietnamese soldiers were dying from malaria and the Vietnamese government asked China for help. Then, the Chinese government launched a research project to search for antimalaria drugs from TCM. As a result, artemisinin was screened out as a strong antimalaria compound from
To verify this hypothesis, we applied the larval Zebrafish model to conduct acute toxicity and efficacy assays. Larval Zebrafish represents a well-established animal model, which suitably serves for toxicity assessment. It offers many advantages compared to traditional
Our results from the acute toxicity assay showed that a single-dose IV injection of ART from 1/9 MNLC to above doses exerted obvious toxicity in the cardiovascular system, resulting in pericardial edema and circulation defects. Besides this obvious cardiotoxicity, ART at nonlethal doses (≤ MNLC) also exerted acute nephrotoxicity and developmental toxicity. Since the average body weight of the larval zebrafishes was 25 mg, the toxic dose range (1/9 MNLC to LD10) can be estimated as 0.144 to 0.932 mg/kg in Zebrafish, which could be extrapolated to 0.0144 to 0.0932 mg/kg in human beings. 64 The recommended daily doses of ART in clinical practice are 1.2 to 2.4 mg/kg human body weight, which is overwhelmingly higher than the above toxic doses extrapolated by the larval Zebrafish assay. The assay discovered such toxicity within 3 days after single injection of ART, indicating that single IV injection of ART potentially induces acute toxicity in its clinical doses for at least 3 days. However, it is generally known that ART is well tolerated and appears to be a safe option even during pregnancy. This controversy to our result can be explained by the following reason: Zebrafish is a very sensitive model suitable for the discovery of potential toxicity or side effects of therapeutics, and ART’s toxicity might have been underestimated due to its invisibility in the clinic. On the other hand, ART at 1/2 LOAEL exerted cardioprotective effects on zebrafishes with heart failure, which significantly restored cardiac malformation, venous stasis, cardiac output decrease, and blood flow dynamics reduction. No adverse events were observed with this treatment, indicating that ART at doses below LOAEL was effective and safe. For the first time, the above results revealed the efficacy-toxicity relationship of ART and verified our hypothesis that ART's dosage decides, whether it is toxic or effective. Such a biphasic dose–response relationship has been described as “hormesis,” which refers to a phenomenon that is characterized by a “low-dose activation and high-dose inhibition”. 65,66 In the context of ART, there was also a hormetic duality in response by zebrafishes in reply to an impetus that spurred favorable effects at a low dose and harmful effects at higher doses.
The RNA-sequencing result showed that 7 genes (
Furthermore, the gene expression of a transmenbrane modulator of Wnt signaling pathway, Frizzled class receptor 7a, was significantly upregulated in heart failure zebrafishes and significantly downregulated by ART at 1/2 LOAEL, indicating that Wnt signaling pathway might mediate the cardioprotective effect of ART. As a classic signaling pathway, Wnt signaling plays an essential role in cardiac disease and serves as an important target for cardiac medication. 67 Artesunate has been reported to exert antimalaria and anticancer effects by inhibiting Wnt signaling pathway, 68,69 indicating Wnt signaling as a target for the treatments of ART. Interestingly, the Wnt signaling pathway also mediates the cardiotoxic effects of other therapeutic agents, such as doxorubicin (DOX). DOX is an effective chemotherapeutic agent in treating human neoplasms such as leukemia, lymphomas, and solid tumors. 70 Acute cardiotoxicity has been recognized as a severe complication with DOX chemotherapy, often resulting in congestive heart failure and irreversible degenerative cardiomyopathy. 71 The cardiotoxicity of DOX was mediated by Wnt signaling in cardiomyocytes. 72 Secreted frizzled-related proteins (sFRP) are ligands that bind to Fz receptors such as Frizzled class receptor 7a. 73 It can relieve DOX-induced cardiotoxicity by modulation of Wnt signaling. 74,75 Therefore, the fact that ART exerts a sFRP-like effect on the cardiovascular disease by targeting Frizzled class receptor 7a in Wnt signaling indicates a therapeutic potential of ART for relieving DOX-induced cardiotoxicity in cancer treatment. It has been reported that some other herbal medicines have successfully alleviated cardiotoxicity in DOX-treated mice with cancer, providing a potential alternative treatment for cancer patients under DOX treatment. 76 It remains to be investigated in the future, whether the potential cardioprotective effects of ART can be exploited to reduce DOX-mediated cardiotoxicity in combination treatment protocols.
Supplemental Material
Supplemental Material, Table_S1 - Cardiotoxicity and Cardioprotection by Artesunate in Larval Zebrafish
Supplemental Material, Table_S1 for Cardiotoxicity and Cardioprotection by Artesunate in Larval Zebrafish by Chuanrui Zheng, Letian Shan, Peijian Tong and Thomas Efferth in Dose-Response
Supplemental Material
Supplemental Material, Table_S2 - Cardiotoxicity and Cardioprotection by Artesunate in Larval Zebrafish
Supplemental Material, Table_S2 for Cardiotoxicity and Cardioprotection by Artesunate in Larval Zebrafish by Chuanrui Zheng, Letian Shan, Peijian Tong and Thomas Efferth in Dose-Response
Footnotes
Authors’ Note
Chuanrui Zheng and Letian Shan contributed equally to this manuscript; they conceived of and designed this study, and wrote the article. Peijian Tong improved the design. Thomas Efferth improved the design, draft, and language editing of this manuscript. All listed authors approved the manuscript for publication and agreed to be accountable for all aspects of this work.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was funded by the National Natural Science Foundation of China (Grant No: 81774331).
Supplemental Material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
